5 tmp<fvScalarMatrix> p_rghEqnComp;
53 max(alpha1, scalar(0))*(psi1/
rho1)
54 +
max(alpha2, scalar(0))*(psi2/
rho2)
73 *(p_rghEqnComp &
p_rgh);
74 phi += p_rghEqnIncomp.flux();
79 U.correctBoundaryConditions();
83 (
p_rgh +
gh*(alpha1*rho10 + alpha2*rho20))
84 /(1.0 -
gh*(alpha1*psi1 + alpha2*psi2)),
88 rho1 = rho10 + psi1*
p;
89 rho2 = rho20 + psi2*
p;