From 82ca430a02e1074df3bccf803f13236fb3daa1dc Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
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 ed59cb37579..d15f78aff58 100644
--- 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_))
);
}
--
GitLab