From 2883311a23e5e6ca8f341217f7e768f940a6773d Mon Sep 17 00:00:00 2001
From: Andrew Heather <a.heather@opencfd.co.uk>
Date: Thu, 15 Jun 2017 14:51:27 +0100
Subject: [PATCH] ENH: rhoPimple[DyM]Foam - updated treatment for dpdt
---
.../compressible/rhoPimpleFoam/createFields.H | 18 ++----------------
.../cfdTools/compressible/createDpdt.H | 18 ++++++++++++++++++
2 files changed, 20 insertions(+), 16 deletions(-)
create mode 100644 src/finiteVolume/cfdTools/compressible/createDpdt.H
diff --git a/applications/solvers/compressible/rhoPimpleFoam/createFields.H b/applications/solvers/compressible/rhoPimpleFoam/createFields.H
index deaba4dd1fd..9d2307d4183 100644
--- a/applications/solvers/compressible/rhoPimpleFoam/createFields.H
+++ b/applications/solvers/compressible/rhoPimpleFoam/createFields.H
@@ -35,8 +35,7 @@ volVectorField U
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
- mesh,
- true
+ mesh
);
#include "compressibleCreatePhi.H"
@@ -57,20 +56,7 @@ autoPtr<compressible::turbulenceModel> turbulence
)
);
-Info<< "Creating field dpdt\n" << endl;
-volScalarField dpdt
-(
- IOobject
- (
- "dpdt",
- runTime.timeName(),
- mesh,
- IOobject::READ_IF_PRESENT,
- IOobject::AUTO_WRITE
- ),
- mesh,
- dimensionedScalar("dpdt", p.dimensions()/dimTime, 0)
-);
+#include "createDpdt.H"
#include "createK.H"
diff --git a/src/finiteVolume/cfdTools/compressible/createDpdt.H b/src/finiteVolume/cfdTools/compressible/createDpdt.H
new file mode 100644
index 00000000000..7c730069247
--- /dev/null
+++ b/src/finiteVolume/cfdTools/compressible/createDpdt.H
@@ -0,0 +1,18 @@
+Info<< "Creating field dpdt\n" << endl;
+volScalarField dpdt
+(
+ IOobject
+ (
+ "dpdt",
+ runTime.timeName(),
+ mesh,
+ IOobject::NO_READ,
+ IOobject::NO_WRITE
+ ),
+ fvc::ddt(p)
+);
+
+if (!thermo.dpdt())
+{
+ dpdt == dimensionedScalar("0", dpdt.dimensions(), 0);
+}
--
GitLab