[wang219910611]

Hi everyone,

I was trying to do myCavitatingFoam for the last few days. But I met some problem when I was trying to write the 0pEqn.

Following is my 0pEqn:

{
volScalarField pgl=(rhov*rhol*(rhov-rhol))/(magSqr(rhov)*psil-magSqr(rhol)*psiv);
if ( rho>rhol)

p=pSat+(B/N)*(pow((rho/rhol),N)-1.0);

else if (rho<rhov)

p=pSat*pow((rho/rhov),r);

else

p=pSat+pgl*log10(rhov*psil*(rhol+alphav*(rhov-rhol))/(rhol*(rhov*psil-alphav*(rhov*psil-rhol*psiv))));

}

Apparently, this is not in OpenFoam style. And the system told me that I should use a "bool" function instead.

Anybody has any idea of this? Any help will be appreciated.

[alexeym]

Hi,

surely people can try to guess but...

0. what is rho?
1. what is rhol?
2. what is rhov?

In fact the question is are they scalars or fields?

[wang219910611]

Originally Posted by alexeym

Hi,

surely people can try to guess but...

0. what is rho?
1. what is rhol?
2. what is rhov?

In fact the question is are they scalars or fields?

Thanks for the reply. I have checked their dimension again and found the mistake.

rho is the density, volScalarField
rhov is the gas density, dimensionedScalar(which should be volScalarField)
rhol is the liquid density, dimensionedScalar(which should be volScalarField)

Now the compile goes well. I'm going to test the solver with case now.