[ma-tri-x]

Hi everyone!

I want to calculate a bubble with compressibleInterFoam with version 2.2.1.
In the TermoPhysicalPropertiesair or ""water you can specify the equationOfState that is used. However it seems that adiabaticPerfectFluid (provided with the source code) is not a valid option. If I type "bla" into the entry and look at the list in the error message, only incompressiblePerfectFluid, PerfectGas, PerfectFluid, icoPolynomial are shown as valid. Is this an error in the code or am I wrong? I would like to use the equation from adiabaticPerfectFluid.

Thanks!

[Tscar]

Hi,
did you find solution to this?

[ma-tri-x]

Yes, it was my Master-Thesis. But I changed to OpenFoam extend 3.1. And changed the solver to my needs. Never worked with commercial OpenFoam anymore since then.

[Tscar]

Thank you for your reply. May I ask you for a case where you used the adiabaticPerfectFluid model? It would surely save me a few hours.

[ma-tri-x]

Sorry, as I said, I used openfoam extend, which doesn't have this kind of structure. I didn't solve it with openfoam 221. What I did was taking the compressibleInterFoam from openfoam extend and then modified the equations of state p(rho) such that in pEqn.H and createFields.H they resemble an adiabatic one for the gas and the compressible Tait Equation for the liquid. See my article http://www.sciencedirect.com/science...45793015003795

Regards,
Max

[chubb87]

For me it is possible to use adiabaticPerfectFluid with standard OpenFoam version (2.4.0), which is the same as Tait equation I think.

Regards

[mg.mithun]

Hi Max,
Can you share your thesis with me
mg.mithun@gmail.com

[bigorange]

Originally Posted by ma-tri-x

Yes, it was my Master-Thesis. But I changed to OpenFoam extend 3.1. And changed the solver to my needs. Never worked with commercial OpenFoam anymore since then.

Hi, ma-tri-x! What do you mean by saying " commercial OpenFoam"? Do you mean that the OpenFoam extend is a better choice if one wants to use OpenFoam by modifing the exist solvers easily?

Regards,
big orange

[ma-tri-x]

Originally Posted by mg.mithun

Hi Max,
Can you share your thesis with me
mg.mithun@gmail.com

Hi mg.mithun!

Have I already replied to your request? If not, here's a link to download the thesis:
https://www.dropbox.com/s/rm9bctajnx...ation.pdf?dl=0

Best regards,
Max

[ma-tri-x]

Originally Posted by bigorange

Hi, ma-tri-x! What do you mean by saying " commercial OpenFoam"? Do you mean that the OpenFoam extend is a better choice if one wants to use OpenFoam by modifing the exist solvers easily?

Regards,
big orange

I am not putting judgement. The basis of the solvers is the same, as I see it. And the package is weight too complex to judge what branch of the software works better in handling own code development, especially when you're not a professional but only skilled c++ programmer.
The point is, that the solver code differs in terms of putting the fluids into categories (resulting into some different handling of the top-level classes in the object-oriented structure) from extend to openfoam. No judgement.

[bigorange]

Originally Posted by ma-tri-x

I am not putting judgement. The basis of the solvers is the same, as I see it. And the package is weight too complex to judge what branch of the software works better in handling own code development, especially when you're not a professional but only skilled c++ programmer.
The point is, that the solver code differs in terms of putting the fluids into categories (resulting into some different handling of the top-level classes in the object-oriented structure) from extend to openfoam. No judgement.

Hi, ma-tri-x!

Thanks for your reply! I understand what you mean!

[Open_Foam]

Hi ma-tri-x,
I found your thesis very interesting but I was wondering how you added the Tait equation. Is there any source code available online for that. If not, can you share any test case for the same to my email (ram22071989@gmail.com) ?
:-)

[ma-tri-x]

Just replace the linear Equation of state in pEqn.H:

rho1 = rho10*pow((p+B1)/(pRef+B1),1.0/gamma1.value());

[Open_Foam]

Thank you for your reply but is this the only thing required to implement TAIT equation of state for water as suggested in your thesis shared above?
Thank you again

[Didu]

I'm also trying to use the adiabaticPerfectFluid, in theory is the best, but how can I change the value of the constant K (Bulk Modulus)? since K=rho/psi