From fb57ecbfdf0db3ff4362eed7365724f3bb72247c Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Wed, 14 Dec 2016 10:42:49 +0000 Subject: [PATCH] functionObjects::forceCoeffs: Corrected setting of rhoRef for compressible flow Resolves bug-report https://bugs.openfoam.org/view.php?id=2387 --- src/functionObjects/forces/forceCoeffs/forceCoeffs.C | 3 +++ src/functionObjects/forces/forces/forces.C | 2 +- 2 files changed, 4 insertions(+), 1 deletion(-) diff --git a/src/functionObjects/forces/forceCoeffs/forceCoeffs.C b/src/functionObjects/forces/forceCoeffs/forceCoeffs.C index cab68d60f..24adeedd3 100644 --- a/src/functionObjects/forces/forceCoeffs/forceCoeffs.C +++ b/src/functionObjects/forces/forceCoeffs/forceCoeffs.C @@ -156,6 +156,9 @@ bool Foam::functionObjects::forceCoeffs::read(const dictionary& dict) // Free stream velocity magnitude dict.lookup("magUInf") >> magUInf_; + // Reference (free stream) density + dict.lookup("rhoInf") >> rhoRef_; + // Reference length and area scales dict.lookup("lRef") >> lRef_; dict.lookup("Aref") >> Aref_; diff --git a/src/functionObjects/forces/forces/forces.C b/src/functionObjects/forces/forces/forces.C index 440b5645b..872efa78b 100644 --- a/src/functionObjects/forces/forces/forces.C +++ b/src/functionObjects/forces/forces/forces.C @@ -639,7 +639,7 @@ bool Foam::functionObjects::forces::read(const dictionary& dict) // Reference density needed for incompressible calculations if (rhoName_ == "rhoInf") { - rhoRef_ = readScalar(dict.lookup("rhoInf")); + dict.lookup("rhoInf") >> rhoRef_; } // Reference pressure, 0 by default