[nathan.kim]

Hi,
While I am validating the cavitatingFOAM to solve acoustic cavitation under an ultrasonic cavitation device, I have a question on the source code of pEqn.C.

I understand it tries to solve pEqn using variables in the previous time step and update the variables accordingly. So I expect the red marked statement should be

rho = max ( ... )

instead of '=='.

If anybody can tell me if this is correct, before I really change it, it will be so thankful.

BG

---------------(quote begins)-------------------
while (pimple.correctNonOrthogonal())
{
fvScalarMatrix pEqn
(
...
);

pEqn.solve(mesh.solver(p.select(pimple.finalInnerI ter())));

if (pimple.finalNonOrthogonalIter())
{
phi += (phiGradp + pEqn.flux())/rhof;
}
}

Info<< "Predicted p max-min : " << max(p).value()
<< " " << min(p).value() << endl;

rho == max
(
psi*p
+ alphal*rhol0
+ ((alphav*psiv + alphal*psil) - psi)*pSat,
rhoMin
);

#include "alphavPsi.H"

p =
(
rho
- alphal*rhol0
- ((alphav*psiv + alphal*psil) - psi)*pSat
)/psi;

p.correctBoundaryConditions();
-------------(quote ends)-----------------------------

[nathan.kim]

I found the operator '==' overrides C/C++ standard operator in another search. The operator is redefined in fvMatrix.