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

transonicEntrainmentPressureFvPatchScalarField: Apply a fixed total-pressure condition for outflow

Browse Source
This commit is contained in:
Henry Weller
2021-12-06 11:33:48 +00:00
parent 31fc0249cc
commit 02a40ce77b
1 changed files with 7 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
src/finiteVolume/fields/fvPatchFields/derived/transonicEntrainmentPressure/transonicEntrainmentPressureFvPatchScalarField.C
View File

@ -173,8 +173,13 @@ void Foam::transonicEntrainmentPressureFvPatchScalarField::updateCoeffs()
const scalar Mb = 0.5;
refValue() =
p0_/pow(scalar(1) - (0.5*gM1ByG)*psip*Unp*mag(Unp), 1/gM1ByG);
valueFraction() = 1 - min((Ma - Mb)/(1 - Mb), scalar(1));
p0_
/pow
(
1 - (0.5*gM1ByG)*psip*negPart(Unp)*mag(Unp),
1/gM1ByG
);
valueFraction() = 1 - min(max(Ma - Mb, scalar(0))/(1 - Mb), scalar(1));
mixedFvPatchField<scalar>::updateCoeffs();
}