[P.Satish]

Hi!

I am new to openfoam and I have a very small query.

In cfd solvers, there is an option for Operation pressures, where I can enter 1 atm so at the boundaries I can enter the guage pressure as the operating pressure is 1atm.

Where is the option of operating pressure in OpenFOAM.

In unsteady solver in fvSolution folder, there its clearly mentioned for solvers
p{
}
U{
}
PISO
{
nCorrectors 2;
nNonOrthogonal Correctors 0;
pRefCell 0;
pRefValue 0;
}

so in the pressure velocity coupling PISO, there pressure reference value is showing 0, so what ever the pressure boundary value I enter at inlet and outlet makes the pressure absolute one. Please correct me whether the above pRefValue is same as the operating pressure in Fluent or not.

When it comes to steady state case

when I had gone to the similar folder

it shows

SIMPLE
{ nNonOrthogaonalCorrectors 0;
}

there is no pref value entered there?

do I have to type the same there, like what is there in PISO in unsteady solver?

or by default Reference pressure is zero internally and what ever the value of pressure I enter at inlet/outlet boundaries becomes absolute pressure in steady state cases?

Please help me in clarifying my query.

[alberto]

Hi,

in OpenFOAM pRef is used only if your boundary conditions are set in such a way not to prescribe any reference for the pressure.

So, for example, if you open a duct open to the athmosphere at the outlet, and you have the velocity specified at the inlet, you should specify p equal to zero at the outlet (in incompressible flows, if you specify 1 it doesn't make any difference, you simply translate the values of the pressure in your whole domain of the same difference).

If you have a simulation where p has only Neumann (zeroGradient) conditions, the pRef value will be used, setting p = pRefVal in the cell of index pRefCell.

The values of pRefCell and pRefValue are always specified in the PISO/SIMPLE/PIMPLE subdictionary in fvSolution.

Best,

[P.Satish]

Alberto Thanks for ur reply.

You are saying the pref will be used only when we give the condition pgradient = 0 at the inlet/outlet. and it wont use if we directly specify the outlet/inlet pressure as 0? Please correct me if I am wrong.

In fluent, if we specify pref=1atm, then if we entering a guage perssure of
0 atm at outlet then we can directly enter the value 0.

even if I am giving the perssure boundary condition p=0 at outlet, there wont be any wrong in specifying the pref value = 0 in tne fvSolution file, right?

when I am using a steady state case, do I have to clearly enter "nCorrectors" as 1, or by default it takes as one inside?

what is the value given for "nNonOrthogonal Correctors" for a tetra mesh?

[alberto]

Originally Posted by P.Satish

Alberto Thanks for ur reply.

You are saying the pref will be used only when we give the condition pgradient = 0 at the inlet/outlet. and it wont use if we directly specify the outlet/inlet pressure as 0? Please correct me if I am wrong.

This actually depends on the solver. In some solver (i.e. pisoFoam) the reference is set using directly

Eqn.setReference(pRefCell, pRefValue);

In others the reference is set only if required.

In fluent, if we specify pref=1atm, then if we entering a guage perssure of
0 atm at outlet then we can directly enter the value 0.

even if I am giving the perssure boundary condition p=0 at outlet, there wont be any wrong in specifying the pref value = 0 in tne fvSolution file, right?

Right, since you're solving for incompressible flows. As long as you do not over-specify the boundary (setting p and U at the same side, for example), it will be fine. The only difference could be a translation of the values of p.

when I am using a steady state case, do I have to clearly enter "nCorrectors" as 1, or by default it takes as one inside?

Yes, the dictionary entry is required.

what is the value given for "nNonOrthogonal Correctors" for a tetra mesh?

It depends on the mesh quality. Use checkMesh. If the skewness is OK, you probably do not need correctors. If the skeweness is high, you need them. The generally suggested value is between 1 and 3.

Best,

[pizicai]

Hi, alberto, I have confusion about this : " setting p = pRefVal in the cell of index pRefCell."

what cell is the pRefCell? pRefCell is usually set to 0, does this mean the pressure of NO.0 cell is zero? but how come I know NO.0 cell locates at the outlet in your case?

if this question is not described clear enough to you, please tell me.

Thanks.

[Fine]

Hi Yalong

You probably figured this out by now...just in case:

An alternative to pRefCell is to specify pRefPoint i.e.
pRefPoint (0 0 0);

Or, if you want to specify pRefCell you can open your grid in parafoam. In the menu go to Edit and then go to "Find data". There you can enter the cell ID as 0 and it will be displayed on your mesh with a coloured dot.

Set your view in paraFoam to solid colour and surface with edges...then you can see the cell more clearly

good luck
f

[maysmech]

Hi all,
In a fluidized bed which the top is set zero constant pressure, it is predicted to have positive pressure in the inlet, because of pressure drop in the direction from inlet to outlet. the mean time averaged of pressure is positive in inlet but the instantaneous pressure fluctuates between large positive and large negative values from time step to time step.
I saw this trend in many solvers like pisoFoam, DPMFoam, ... which the instantaneous pressure (P) sign is not similar to what we predict. Again, This is for cases which we set outlet as constant pressure. Any idea?

[crazzy.pirate43]

Originally Posted by maysmech

Hi all,
In a fluidized bed which the top is set zero constant pressure, it is predicted to have positive pressure in the inlet, because of pressure drop in the direction from inlet to outlet. the mean time averaged of pressure is positive in inlet but the instantaneous pressure fluctuates between large positive and large negative values from time step to time step.
I saw this trend in many solvers like pisoFoam, DPMFoam, ... which the instantaneous pressure (P) sign is not similar to what we predict. Again, This is for cases which we set outlet as constant pressure. Any idea?

Hi maysmech,

It was long time for this post, but I have the same problem now as in my case the Mean pressure is always from positive at the inlet of the channel ro zero at outlet while the instantinious pressure fluctuate between positive and negative values at the inlet to zero at the outlet? So i was wondering if you found a solution fro the your problem?????

Thanks inadvance