From 608da26c5deedb4df918bc4bc4db97935db60504 Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Wed, 25 May 2016 18:12:07 +0100 Subject: [PATCH] thermoSingleLayer: Provide better stabilization for the energy equation as delta -> 0 --- .../thermoSingleLayer/thermoSingleLayer.C | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C b/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C index ed59cb375..d15f78aff 100644 --- a/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C +++ b/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C @@ -265,12 +265,12 @@ tmp thermoSingleLayer::q(volScalarField& hs) const return ( // Heat-transfer to the primary region - - fvm::Sp(alpha_*htcs_->h()/Cp_, hs) - + alpha_*htcs_->h()*(hs/Cp_ - T_ + TPrimary_) + - fvm::Sp(htcs_->h()/Cp_, hs) + + htcs_->h()*(hs/Cp_ + alpha_*(TPrimary_ - T_)) // Heat-transfer to the wall - - fvm::Sp(alpha_*htcw_->h()/Cp_, hs) - + alpha_*htcw_->h()*(hs/Cp_ - T_ + Tw_) + - fvm::Sp(htcw_->h()/Cp_, hs) + + htcw_->h()*(hs/Cp_ + alpha_*(Tw_- T_)) ); }