From c45f1897506eac8a6dac912687d6bbc0cc609392 Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Mon, 7 Feb 2022 18:48:15 +0000 Subject: [PATCH] dynamicPressureFvPatchScalarField: Updated to use the patch density when evaluating the external total pressure from the static pressure and tangential velocity. This improves numerical stability and solution smoothness for compressible cases. --- .../dynamicPressure/dynamicPressureFvPatchScalarField.C | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/finiteVolume/fields/fvPatchFields/derived/dynamicPressure/dynamicPressureFvPatchScalarField.C b/src/finiteVolume/fields/fvPatchFields/derived/dynamicPressure/dynamicPressureFvPatchScalarField.C index c60e1422bd..2976730668 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/dynamicPressure/dynamicPressureFvPatchScalarField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/dynamicPressure/dynamicPressureFvPatchScalarField.C @@ -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("rho", "rho")), psiName_(dict.lookupOrDefault("psi", "none")), - gamma_(psiName_ != "none" ? dict.lookup("gamma") : 1), + gamma_(dict.lookupOrDefault("gamma", 1)), p0_("p0", dict, p.size()) { if (dict.found("value"))