thermoSingleLayer: revert change to q function made shortly before the release of OpenFOAM-2.3.1

Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1502
This commit is contained in:
Henry
2015-01-28 18:02:56 +00:00
parent eb4ed526cc
commit ad03953433

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-2014 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -257,16 +257,9 @@ tmp<fvScalarMatrix> thermoSingleLayer::q(volScalarField& hs) const
{ {
dimensionedScalar Tstd("Tstd", dimTemperature, 298.15); dimensionedScalar Tstd("Tstd", dimTemperature, 298.15);
volScalarField htcst(htcs_->h()); volScalarField boundedAlpha(max(alpha_, ROOTVSMALL));
volScalarField htcwt(htcw_->h()); volScalarField htcst(htcs_->h()*boundedAlpha);
volScalarField htcwt(htcw_->h()*boundedAlpha);
const volScalarField mask(pos(delta_ - deltaSmall_));
forAll(mask, i)
{
htcst[i] *= max(mask[i], ROOTVSMALL);
htcwt[i] *= max(mask[i], ROOTVSMALL);
}
htcst.correctBoundaryConditions(); htcst.correctBoundaryConditions();
htcwt.correctBoundaryConditions(); htcwt.correctBoundaryConditions();