externalWallHeatFluxTemperature: Use the surface emissivity for both emission and absorption

Resolves bug-report https://bugs.openfoam.org/view.php?id=2953
This commit is contained in:
Henry Weller
2018-05-23 14:14:01 +01:00
parent 7601e8e86f
commit 6c560795cb

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-2017 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2018 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -405,14 +405,14 @@ void Foam::externalWallHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
scalarField lambdaTa4(pow4((1 - TpLambda)*Ta)); scalarField lambdaTa4(pow4((1 - TpLambda)*Ta));
hp += emissivity_*sigma.value()*(pow4(Ts) - lambdaTa4)/Tp; hp += emissivity_*sigma.value()*(pow4(Ts) - lambdaTa4)/Tp;
hpTa += sigma.value()*(emissivity_*lambdaTa4 + pow4(Ta)); hpTa += emissivity_*sigma.value()*(lambdaTa4 + pow4(Ta));
} }
else else
{ {
// ... if there is no solid wall thermal resistance use // ... if there is no solid wall thermal resistance use
// the current wall temperature // the current wall temperature
hp += emissivity_*sigma.value()*pow3(Tp); hp += emissivity_*sigma.value()*pow3(Tp);
hpTa += sigma.value()*pow4(Ta); hpTa += emissivity_*sigma.value()*pow4(Ta);
} }
} }