move update of rhoeps to pEqn.H

This commit is contained in:
danielque
2021-12-23 15:29:31 +01:00
parent b8ccfb9986
commit 4bdb5f06d4
2 changed files with 3 additions and 1 deletions

View File

@ -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())

View File

@ -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);