[DanielRCalvete]

Hello,

I am running cavitation cases with interPhaseChangeFoam, and the tool explaned by Andrea is very useful to me. I am trying to implement it in OF 2.3, but I have a doubt about it:

I see in this code that Pressure Forces is integrated multiplying by rho(p) which is 1 if the pDimensions = Pascals according to this piece of code:

 Foam::scalar Foam::forceMultiPhase::rho(const volScalarField& p) const
{
    if (p.dimensions() == dimPressure)
    {
        return 1.0;
    }
    else
    {
        if (rhoName_ != "rhoInf")
        {
            FatalErrorIn("forceMultiPhase::rho(const volScalarField& p)")
                << "Dynamic pressure is expected but kinematic is provided."
                << exit(FatalError);
        }

        return rhoRef_;
    }
}

Maybe I am missunderstanding something in the code of Andrea, but since the pressure forces is the integration of pressure through the surface, I think it has not sense to multiply by rho*alpha or rho*(1- alpha) to obtain this forces.

What I propose is the next piece of code in forceMultiPhace.C in calcForcesMoment() class (notice that is for OF2.3):

	 vectorField fN
            (
                Sfb[patchI]*(p.boundaryField()[patchI] - pRef)
*
            );
	// Viscous Forces

            vectorField fT

	    (
		((Sfb[patchI]*limitedAlpha1.boundaryField()[patchI]) & devRhoReffb1[patchI])
		+ ((Sfb[patchI]*(scalar(1)-limitedAlpha1.boundaryField()[patchI])) & devRhoReffb2[patchI]) 
	
	    );

I will appreciate a lot your comments about. Thank you in advance.

DRC