thermoSingleLayer: Provide better stabilization for the energy equation as delta -> 0

2016-05-25 18:12:07 +01:00
parent ef3c7c045d
commit 608da26c5d
1 changed files with 4 additions and 4 deletions
--- a/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C
+++ b/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C
@ -265,12 +265,12 @@ tmp<fvScalarMatrix> 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_))
    );
 }