[rafigh]

Dear FOAMers
I'm new in OpenFoam and would like to model a pressure driven flow through a plug nozzle using sonicfoam.
geometry, grid and velocity plot for nozzle pressure 4 bar and ambient pressure 3.5 bar , attached. when i decrease ambient pressure to 3 bar solver diverging.
1.png

2.png

3.png

BC

for pdimensions [1 -1 -2 0 0 0 0];

internalField uniform 350000;

boundaryField
{
gasInlet
{

type totalPressure;
p0 uniform 400000;
U U;
phi phi;
rho none;
psi none;
gamma 1.3;
value uniform 400000;
}

airInlet
{

type totalPressure;
p0 uniform 350000;
U U;
phi phi;
rho none;
psi none;
gamma 1.3;
value uniform 350000;
}

outlet
{
type waveTransmissive;
field p;
phi phi;
rho rho;
psi thermosi;
gamma 1.3;
fieldInf 350000;
lInf 1;
value uniform 350000;
}

symPlane
{
type symmetry;
}

freestreem
{
type zeroGradient;
}

spikeWall
{
type zeroGradient;
}

frontAndBack
{
type empty;
}
}
///////////////////////////////
for U

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
gasInlet
{

type pressureInletOutletVelocity;
phi phi;
rho rho;
value uniform (0 0 0);
}

airInlet
{

type pressureInletOutletVelocity;
phi phi;
rho rho;
value uniform (0 0 0);
}

outlet
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0 0 0);

}



symPlane
{
type symmetry;
}

freestreem
{
type supersonicFreestream;
pInf 350000;
TInf 300;
UInf (650 0 0);
gamma 1.3;
value uniform (0 0 0);
}

spikeWall
{
type fixedValue;
value uniform (0 0 0);
}

frontAndBack
{
type empty;
}
}



Many thanks for your help.
regards

[ksfbangbangda]

Hello, I want to use waveTransmissive BC, but I don't konw how to use it? can you give me a case? 1125449167@qq.com

[xshmuel]

Hi I am quite new to OpenFOAM and I was wonderign how you got your geometry of the plug nozzle so nice? i just learned how to use the blockmesh tool butthis seems like quite a cumbersome job to get such a nice gemetry as you got. so how did you manage?