zhyang-dev/CFDEMcoupling-PFM
1
0
Fork 0
mirror of https://github.com/ParticulateFlow/CFDEMcoupling-PFM.git synced 2025-12-08 06:37:44 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

move recalculation of rho in pEqn.H

Browse Source 
just move calculation further down
add average value to output of rho
This commit is contained in:
danielque
2021-12-23 15:37:12 +01:00
parent 54db2db656
commit e16e405828
1 changed files with 8 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
applications/solvers/cfdemSolverRhoPimple/pEqn.H
View File

@ -67,14 +67,6 @@ else
// Explicitly relax pressure for momentum corrector // Explicitly relax pressure for momentum corrector
p.relax(); p.relax();
// Recalculate density from the relaxed pressure
rho = thermo.rho();
rho = max(rho, rhoMin);
rho = min(rho, rhoMax);
rho.relax();
Info<< "rho max/min : " << max(rho).value()
<< " " << min(rho).value() << endl;
if (modelType=="A") if (modelType=="A")
{ {
U = HbyA - rAU*(voidfraction*fvc::grad(p)-Ksl*Us); U = HbyA - rAU*(voidfraction*fvc::grad(p)-Ksl*Us);
@ -87,8 +79,16 @@ U.correctBoundaryConditions();
fvOptions.correct(U); fvOptions.correct(U);
K = 0.5*magSqr(U); K = 0.5*magSqr(U);
// Recalculate density from the relaxed pressure
rho = thermo.rho();
rho = max(rho, rhoMin);
rho = min(rho, rhoMax);
rho.relax();
rhoeps = rho*voidfraction; rhoeps = rho*voidfraction;
Info<< "rho max/min/ave : " << max(rho).value()
<< " " << min(rho).value() << " " << average(rho).value() << endl;
if (thermo.dpdt()) if (thermo.dpdt())
{ {
dpdt = fvc::ddt(voidfraction,p); dpdt = fvc::ddt(voidfraction,p);