Initial condition in OpenFoam
 

Hello,

I want to initialize hydrostatic pressure in my simulation. I have tried it to do it by 'funkySetFields -time 0'. It works for zero time but at first iteration it suddenly changes.

I have pressure inlet and pressure outlet as a BC's. In 0/p I have given fixed values for pressure uniform 101320(for outlet) and uniform 118600 (for inlet).

Can anyone help me in this problem ?

Best regards,
Gitesh

What solver are you using?

Can you provide a little more information about what you have tried? Step by step?

Did you use setFields to initialize water surface?

Hello mgdenno,

OK. So my system is air is feeding by pipe into water tank. I am using twoPhaseEulerFoam solver.

I am using 'funkySetFields' for case initialization. It works ok. But at after first iteration it calculate different things. You can see in pictures (pressure) in attached file.

Best regards,
Gitesh

Hi Gitesh,

I am not at all familiar with twoPhaseEulerFoam, but which variables did you initialize? pressure? phase?

Looks like maybe you only initialized the pressure but not the phase?

Also which way is gravity acting?

MD

Hello,

I also initialized phase. It works ok but the problem in pressure.

Gravity is in +ve x direction.

BR,
Gitesh

This is not a correct setup for an incompressible code. You should specify a velocity at the inlet. If you do that, you also do not need to initialize the pressure gradient.

Best,

Hello alberto,

I want to initialize hydrostatic pressure for water. Also I have to give pressure inlet value because air is coming from there. Moreover I have pressure outlet there from where air will be out.

So, I have used funkySetFields for hydrostatic pressure initialization. So, is there any other bc from where we can directly calculate pressure gradient with out initialization ?

Thank you !!

BR,
Gitesh

Setting the pressure at the outlet is fine. Setting it also at the inlet is going to give you troubles. You should specify the velocity of air. If you do this, initializing the hydrostatic pressure becomes irrelevant, since at the first iteration the code will find the appropriate pressure field.

Best,

Hello Albarto,

Thank you !!

So, you mean I have to give velocity of air at inlet and to initialize pressure gradient by funkySetFields ?

BR,
Gitesh

My suggestion is to:

1) Specify velocity at the inlet, and set the condition on the pressure to zeroGradient there.

2) Fix the value of the pressure at the outlet, and set the velocity to zeroGradient or inletOutlet, depending on your case.

3) You can use a uniform initialization for the pressure in this case, since it does not really matter. The solver will find the correct pressure filed at the first time step.

Best,

hi alberto,
i am working on 2D multi element airfoils, i have a Question about 2th suggestion.
"Fix the value of the pressure at the outlet, and set the velocity to zeroGradient or inletOutlet, depending on your case."
what do you mean with depending on your case,in my case which one is better? inletOutlet or fixedValue ?

hi
inletOutlet is an "outlet" BC not inlet.this is the same zeroGradient only with a difference.
When there is a backflow in outlet,inletOutlet BC sets a value on the cells that have backflow by the value we have specified.
So if you think may you have backflow use inletOulet otherwise zeroGradient suffices.
Hope it helps.

hi ehsan,
i understand what you said, thank you very much.
if my domain is large enough around my airfoil, it doesn’t need to use the inletOutlet boundary condition, is it right?

yes,inletOutlet is used in internal flows in common.
I didn't received your case you said that had sent.

hi Dear ehsan,
i put my multi element airfoil,
thank you very much for your help.

hi friends! i have the problem of hydrostatic pressure initializing in a tank full of water like Gitesh P!
and what is funkysetField? how can i set hydrostatic pressure condition in my domain ? i shoul say my inlet is not the whole left side of my geometry , it is a little nozzle in there (left side of geometry)
so what can i do?