From a271fd0aa0cd2c0abb0efbdd5c72da3eeaf58bc3 Mon Sep 17 00:00:00 2001 From: danielque Date: Thu, 23 Dec 2021 15:22:40 +0100 Subject: [PATCH] create rhoeps field in createFields.H this was a major issue: rhoeps must be created outside the run loop and have old time values, otherwise any time derivative will be zero, leading to incorrect results --- .../solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C | 1 - applications/solvers/cfdemSolverRhoPimple/createFields.H | 3 +++ 2 files changed, 3 insertions(+), 1 deletion(-) diff --git a/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C b/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C index f7d3d217..984b52f6 100644 --- a/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C +++ b/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C @@ -122,7 +122,6 @@ int main(int argc, char *argv[]) #include "rhoEqn.H" } - volScalarField rhoeps("rhoeps",rho*voidfraction); // --- Pressure-velocity PIMPLE corrector loop while (pimple.loop()) { diff --git a/applications/solvers/cfdemSolverRhoPimple/createFields.H b/applications/solvers/cfdemSolverRhoPimple/createFields.H index 23ddbab4..6761ccc8 100644 --- a/applications/solvers/cfdemSolverRhoPimple/createFields.H +++ b/applications/solvers/cfdemSolverRhoPimple/createFields.H @@ -65,6 +65,9 @@ Info<< "Reading thermophysical properties\n" << endl; dimensionedScalar("zero", dimensionSet(1,-1,-3,0,0,0,0), 0.0) ); + volScalarField rhoeps("rhoeps", rho*voidfraction); + rhoeps.oldTime(); // switch on saving old time + Info<< "\nCreating fluid-particle heat flux field\n" << endl; volScalarField Qsource (