[luca]

Hi all,

I try to build in OpenFOAM a pressure based algorithm for prediciting flows at all speed. I formulate pressure equation in term of pressure-correction like in original SIMPLE procedure.
After I've solved pressure-correction equation i need to update fluxes in order to satisfy local continuity equation.
I use this empirical formulation after I've corrected volScalarField "rho", "p" and "U":

--
p.correctBoundaryConditions();

phi = fvc::interpolate(rho*U) & mesh.Sf();
phi.boundaryField() = ( rho.boundaryField()*U.boundaryField() ) & mesh.Sf().boundaryField();
--

Is it right?!

This procedure works but I'm not sure about that..


Thanks in advance at all OpenFOAM users