8 - ((gamma*psiv + (1.0 - gamma)*psil) -
psi)*
pSat
21 p.boundaryField().updateCoeffs();
25 phiv -= phiGradp/rhof;
42 pEqn.solve(
mesh.solver(
p.name() +
"Final"));
46 pEqn.solve(
mesh.solver(
p.name()));
51 phiv += (phiGradp + pEqn.flux())/rhof;
55 Info<<
"Predicted p max-min : " <<
max(
p).value()
62 + ((gamma*psiv + (1.0 - gamma)*psil) -
psi)*pSat,
72 - ((gamma*psiv + (1.0 - gamma)*psil) -
psi)*pSat
75 p.correctBoundaryConditions();
77 Info<<
"Phase-change corrected p max-min : " <<
max(
p).value()
85 if (
piso.found(
"removeSwirl"))
87 label swirlCmpt(readLabel(
piso.lookup(
"removeSwirl")));
89 Info<<
"Removing swirl component-" << swirlCmpt <<
" of U" <<
endl;
90 U.field().replace(swirlCmpt, 0.0);
93 U.correctBoundaryConditions();