thermoSingleLayer::q: Use a local "alpha" without hydrophilic/phobic adjustment

Resolves bug-report https://bugs.openfoam.org/view.php?id=2605
This commit is contained in:
Henry Weller 2017-07-10 12:08:28 +01:00
parent b5331fdc86
commit c989277b5e

View file

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -262,15 +262,17 @@ void thermoSingleLayer::updateSubmodels()
tmp<fvScalarMatrix> thermoSingleLayer::q(volScalarField& hs) const tmp<fvScalarMatrix> thermoSingleLayer::q(volScalarField& hs) const
{ {
const volScalarField alpha(pos(delta_ - deltaSmall_));
return return
( (
// Heat-transfer to the primary region // Heat-transfer to the primary region
- fvm::Sp(htcs_->h()/Cp_, hs) - 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 // Heat-transfer to the wall
- fvm::Sp(htcw_->h()/Cp_, hs) - fvm::Sp(htcw_->h()/Cp_, hs)
+ htcw_->h()*(hs/Cp_ + alpha_*(Tw_- T_)) + htcw_->h()*(hs/Cp_ + alpha*(Tw_- T_))
); );
} }