zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

BUG: pressure FO - corrected handling of enumeration in calcPressure function. See #1356

Browse Source
This commit is contained in:
Andrew Heather
2019-07-09 11:44:52 +01:00
parent 1159aaa810
commit 5bd7c4b5bb
1 changed files with 32 additions and 40 deletions
Download Patch File Download Diff File Expand all files Collapse all files

20
src/functionObjects/field/pressure/pressure.C
View File

@ -146,16 +146,15 @@ Foam::tmp<Foam::volScalarField> Foam::functionObjects::pressure::calcPressure
const tmp<volScalarField>& tp const tmp<volScalarField>& tp
) const ) const
{ {
switch (mode_) if (mode_ & TOTAL)
{
case TOTAL:
{ {
return return
tp tp
+ dimensionedScalar("pRef", dimPressure, pRef_) + dimensionedScalar("pRef", dimPressure, pRef_)
+ rhoScale(p, 0.5*magSqr(lookupObject<volVectorField>(UName_))); + rhoScale(p, 0.5*magSqr(lookupObject<volVectorField>(UName_)));
} }
case ISENTROPIC:
if (mode_ & ISENTROPIC)
{ {
const basicThermo* thermoPtr = const basicThermo* thermoPtr =
p.mesh().lookupObjectPtr<basicThermo>(basicThermo::dictName); p.mesh().lookupObjectPtr<basicThermo>(basicThermo::dictName);
@ -177,13 +176,8 @@ Foam::tmp<Foam::volScalarField> Foam::functionObjects::pressure::calcPressure
return tp()*(pow(1 + (gamma - 1)/2*sqr(Mb), gamma/(gamma - 1))); return tp()*(pow(1 + (gamma - 1)/2*sqr(Mb), gamma/(gamma - 1)));
} }
default:
{ return tp + dimensionedScalar("pRef", dimPressure, pRef_);
return
tp
+ dimensionedScalar("pRef", dimPressure, pRef_);
}
}
} }
@ -207,11 +201,9 @@ Foam::tmp<Foam::volScalarField> Foam::functionObjects::pressure::coeff
return tpCoeff; return tpCoeff;
} }
else
{
return std::move(tp); return std::move(tp);
} }
}
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * // // * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //