diff --git a/src/regionModels/surfaceFilmModels/kinematicSingleLayer/kinematicSingleLayer.C b/src/regionModels/surfaceFilmModels/kinematicSingleLayer/kinematicSingleLayer.C index 9debb2d99..059ae34a7 100644 --- a/src/regionModels/surfaceFilmModels/kinematicSingleLayer/kinematicSingleLayer.C +++ b/src/regionModels/surfaceFilmModels/kinematicSingleLayer/kinematicSingleLayer.C @@ -204,7 +204,7 @@ tmp kinematicSingleLayer::pp() void kinematicSingleLayer::correctAlpha() { - alpha_ == pos0(delta_ - deltaSmall_); + alpha_ == pos(delta_ - deltaSmall_); } diff --git a/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C b/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C index f6e7dd901..2912f029f 100644 --- a/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C +++ b/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C @@ -268,15 +268,17 @@ void thermoSingleLayer::updateSubmodels() tmp thermoSingleLayer::q(volScalarField& hs) const { + const volScalarField alpha(pos(delta_ - deltaSmall_)); + return ( // Heat-transfer to the primary region - fvm::Sp(htcs_->h()/Cp_, hs) - + htcs_->h()*(hs/Cp_ + alpha_*(TPrimary_ - T_)) + + htcs_->h()*(hs/Cp_ + alpha*(TPrimary_ - T_)) // Heat-transfer to the wall - fvm::Sp(htcw_->h()/Cp_, hs) - + htcw_->h()*(hs/Cp_ + alpha_*(Tw_- T_)) + + htcw_->h()*(hs/Cp_ + alpha*(Tw_- T_)) ); }