From 4bdb5f06d4a755aa74c3ee69b8932d77f28af92b Mon Sep 17 00:00:00 2001
From: danielque <daniel.queteschiner@jku.at>
Date: Thu, 23 Dec 2021 15:29:31 +0100
Subject: [PATCH] move update of rhoeps to pEqn.H

---
 .../solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C        | 1 -
 applications/solvers/cfdemSolverRhoPimple/pEqn.H               | 3 +++
 2 files changed, 3 insertions(+), 1 deletion(-)

diff --git a/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C b/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C
index 984b52f6..8f1c2f7c 100644
--- a/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C
+++ b/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C
@@ -133,7 +133,6 @@ int main(int argc, char *argv[])
             {
                 // besides this pEqn, OF offers a "pimple consistent"-option
                 #include "pEqn.H"
-                rhoeps=rho*voidfraction;
             }
 
             if (pimple.turbCorr())
diff --git a/applications/solvers/cfdemSolverRhoPimple/pEqn.H b/applications/solvers/cfdemSolverRhoPimple/pEqn.H
index e107c5a4..d1c99774 100644
--- a/applications/solvers/cfdemSolverRhoPimple/pEqn.H
+++ b/applications/solvers/cfdemSolverRhoPimple/pEqn.H
@@ -2,6 +2,7 @@ rho = thermo.rho();
 rho = max(rho, rhoMin);
 rho = min(rho, rhoMax);
 rho.relax();
+rhoeps = rho*voidfraction;
 
 volScalarField rAU(1.0/UEqn.A());
 surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rhoeps*rAU));
@@ -86,6 +87,8 @@ U.correctBoundaryConditions();
 fvOptions.correct(U);
 K = 0.5*magSqr(U);
 
+rhoeps = rho*voidfraction;
+
 if (thermo.dpdt())
 {
     dpdt = fvc::ddt(voidfraction,p);