[clementolivier]

Hi everyone,

I am a new user on OpenFoam. I would like to run a compressible and biphasic case. It seems that, it is possible using compressibleInterFoam. However, I also want to change the two fluids' equation of state. Do I have to modify the C code of the solver, or there is a simpler way ?

Thanks

Best regards

Clément

[nimasam]

you should change source code, there is no easier task

[clementolivier]

Thank you for your answer, I will get back to work !

[olivierG]

hello,

This is already possible to change the fluid equation of state.
Take a look in the depthCharge tutorial, in TransportProperties : would can choose constant, linear, perfectFluid, ...

Regards,
olivier

[clementolivier]

Thank you for your answer olivierG.
I can see how to change the constants of the model (nu,rho,rho0,psi), but I do not quite see how to change the formula of the equation of state.

Do you mean that the argument of transportModel specifies it ? If it is the case, is there any exhaustive list of all allowed arguments ?

[olivierG]

hello,

First what the state equation you would use ?
And which version of openfoam do you use ?

To get the list of state equation, put "dummy" instead of perfectFluid : you will get a list.
To find more info, you should take a look at the source code. this part was in the solver for the 2.1.x version, but has moved with the 2.2.x

regards,
olivier

[clementolivier]

Thank you again for your help

I want to use the equation of state called the Stifeffened Gas Equation of State for the liquid, and the Perfect Gas Equation of State for the gas.

I am currently usig OpenFoam v2.1.1.

I am sorry but I have done a grep -rin "perfectFluid" * in the directory tutorials/multiphase/compressibleInterFoam, and it is nowhere written. Can you explain more precisely where whould I type "dummy".

I am right now reading the C++ code in order to fix my problem.

[olivierG]

hello,

1) in your tuto case, set "dummy" (or anythink), then launch the solver: you will get an error message, which will say dummy is not a equation of state, and will list the possible one.

2) For the 2.1.1, the definition is in the solver source, never in the tutorial. take a look at applications/solvers/multiphase/compressibleInterFoam/phaseEquationsOfState/ in your installation directories.

regards,

olivier

[clementolivier]

Thanks,

I still do not see where should I set "dummy" (or anything).

Moreover the directory applications/solvers/multiphase/compressibleInterFoam/phaseEquationsOfState/ does not exist in my version. I will get I last version and it will be easier for me to communicate.

regards,

Clément

[olivierG]

hello,

Sorry, i wasn't clear.

copy the tuto (depthCharge3D).
Open depthCharge3D/constant/transportProperties
You have 2 phase, air and water.
In the definition of water, you should have:

equationOfState
{
type perfectFluid;
rho0 1000;
R 3000;
}

Just change "perfectFluid" by "dummy", and run the solver.

regards,
olivier

[clementolivier]

Thank you very much,

Indeed in the last version, I can see everything you are talking about.

I think,this will fix my problem.

Thank you again

Regards,

Clément

[aylalisa]

Hallo Olivier,

I am not clear about equationOfState models.
What meanding has R?

Originally Posted by olivierG

hello,

Sorry, i wasn't clear.

copy the tuto (depthCharge3D).
Open depthCharge3D/constant/transportProperties
You have 2 phase, air and water.
In the definition of water, you should have:

equationOfState
{
type perfectFluid;
rho0 1000;
R 3000;
}

Just change "perfectFluid" by "dummy", and run the solver.

regards,
olivier

Aylalisa

[olivierG]

hello,

R is the specific gaz constant = R/M (M=molar mass).
R is used to calculate the compressibility psi = 1/RT.
For liquid, you still need to set compressibility and use the same law (you can guess this from the name "perfectLiquid"). But R is really bigger, to act as a very low compressible media.

regards,
olivier

[sh.d]

hi
In compressible liquid/gas flows, the pEqn reads :

(
fvm::ddt(psi,p)+fvc::div(phi)+fvm::div(phid,p,".." )-fvm::laplacian(..,p)
)
Can someone tell me what is the form of the differential equation?

[Tobi]

Hello sh.d,

the equation for that formula should be something like that:



I think its not correct but gives you a hint to read the lines.