[peyman.havaej]

Dear All

I use the cavitating flow solver in the OpenFOAM. The cavitatingFoam solver employes a Homogeneous Equilibrium Model (HEM), in which the pressure equation depends on the compressibility of each phase and the density of the mixture is calculated according to an equilibrium barotropic equation.
In the OpenFOAM v-6, there are minor differences in the pressure and density equation between cavitatingFoam and cavitatingDyMFoam (mesh motion).
For example, in pEqn:
cavitatingFoam:

fvScalarMatrix pEqn
        (
            fvm::ddt(psi, p)
          - (rhol0 + (psil - psiv)*pSat)*fvc::ddt(alphav) - pSat*fvc::ddt(psi)
          + fvc::div(phi, rho)
          + fvc::div(phiGradp)
          - fvm::laplacian(rhorAUf, p)
        );

cavitatingDyMFoam:

fvScalarMatrix pEqn
        (
            fvc::ddt(rho)
          + psi*correction(fvm::ddt(p))
          + fvc::div(phi, rho)
          + fvc::div(phiGradp)
          - fvm::laplacian(rhorAUf, p)
        );

for the density of mixture:
cavitatingFoam:

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

cavitatingDyMFoam:

    volScalarField rho0(rho - psi*p);
...
...
...
    rho == max(rho0 + psi*p, rhoMin);

While I can understand the implementation of both solvers, I think pEqn and density correction in the static and Dynamic mesh solvers are not the same.

For clarification,

psi*correction(fvm::ddt(p))

is not equal to the

 fvm::ddt(psi, p)

Moreover, the density correction obviously is not the same. Please correct me if I am wrong.

Do you know what is the main reason behind this difference in the implementation of cavitatingFoam and cavitatingDyMFoam?

Thanks.
Whit Kind Regards,
Peyman

[Yang16024]

Hi, peyman

I am studying the cavitatingFoam solver. But I have a question with this pEqn.
It would be very much appreciated if you can tell me the meaning of term phiGradp in the pEqn and why this term is in the pEqn?
Thank you very much.