zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

dynamicPressureFvPatchScalarField: Updated to use the patch density

Browse Source 
when evaluating the external total pressure from the static pressure and
tangential velocity.  This improves numerical stability and solution smoothness
for compressible cases.
This commit is contained in:
Henry Weller
2022-02-07 18:48:15 +00:00
parent a36082e62f
commit c45f189750
1 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
src/finiteVolume/fields/fvPatchFields/derived/dynamicPressure/dynamicPressureFvPatchScalarField.C
View File

@ -76,12 +76,12 @@ inline void Foam::dynamicPressureFvPatchScalarField::UpdateCoeffs
operator==
(
p0p*(scalar(1) + psip*K0)
/pow(scalar(1) + psip*gM1ByG*Kp, 1/gM1ByG));
p0p/pow(scalar(1) + psip*gM1ByG*(Kp - K0), 1/gM1ByG)
);
}
else
{
operator==(p0p*(scalar(1) + psip*K0)/(scalar(1) + psip*Kp));
operator==(p0p/(scalar(1) + psip*(Kp - K0)));
}
}
}
@ -136,7 +136,7 @@ Foam::dynamicPressureFvPatchScalarField::dynamicPressureFvPatchScalarField
fixedValueFvPatchScalarField(p, iF, dict, false),
rhoName_(dict.lookupOrDefault<word>("rho", "rho")),
psiName_(dict.lookupOrDefault<word>("psi", "none")),
gamma_(psiName_ != "none" ? dict.lookup<scalar>("gamma") : 1),
gamma_(dict.lookupOrDefault<scalar>("gamma", 1)),
p0_("p0", dict, p.size())
{
if (dict.found("value"))