piso algorithm in phaseChange phonomena
 

hi dear foamer
for incompressible flow there is no state equation for pressure
so for updating pressure we should derive equation from continuity and momentum equation. so in piso algorithm as ∇ · u = 0 use the following equation :

fvScalarMatrix pEqn ( fvm::laplacian(rUA, p) == fvc::div(phi) );

however in phase change phenomena ∇ · u = 0 is just globally true but in local ∇ · u = source term

so should i change the pressure equation or not ?