[nimasam]

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 ?