From 6ebaa63d1cd5c5e569fc028c64b8a0d472cdee5a Mon Sep 17 00:00:00 2001
From: Andrew Heather <a.heather@opencfd.co.uk>
Date: Tue, 21 Mar 2017 13:27:47 +0000
Subject: [PATCH] BUG: forceCoeffs - corrected handling of pDyn for
compressible cases. Fixes #436
---
src/functionObjects/forces/forceCoeffs/forceCoeffs.C | 10 +++++++++-
1 file changed, 9 insertions(+), 1 deletion(-)
diff --git a/src/functionObjects/forces/forceCoeffs/forceCoeffs.C b/src/functionObjects/forces/forceCoeffs/forceCoeffs.C
index f4abd527f7e..03bac52926d 100644
--- a/src/functionObjects/forces/forceCoeffs/forceCoeffs.C
+++ b/src/functionObjects/forces/forceCoeffs/forceCoeffs.C
@@ -3,7 +3,7 @@
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
- \\/ M anipulation | Copyright (C) 2015-2016 OpenCFD Ltd.
+ \\/ M anipulation | Copyright (C) 2015-2017 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@@ -250,6 +250,14 @@ bool Foam::functionObjects::forceCoeffs::read(const dictionary& dict)
// Free stream velocity magnitude
dict.lookup("magUInf") >> magUInf_;
+ // If case is compressible we must read rhoInf (store in rhoRef_) to
+ // calculate the reference dynamic pressure
+ // - note: for incompressible, rhoRef_ is already initialised
+ if (rhoName_ != "rhoInf")
+ {
+ dict.lookup("rhoInf") >> rhoRef_;
+ }
+
// Reference length and area scales
dict.lookup("lRef") >> lRef_;
dict.lookup("Aref") >> Aref_;
--
GitLab