[Nicolas_OF]

Hi,

I'm new to OpenFoam (2.1.0). I try to simulate underexpanded axisymmetric supersonic flow of Argon. I used tutorial model : LadenburgJet60Psi.My simulation seems to work when difference between inlet and chamber pressure is small (factor of 10). But when this ratio becomes of the order of 1000, I get the famous error "floating exception point". I attempt different parameters, boundaries and I don't success to find where is my error.

Can Openfoam solve this kind of problem ?


Thank you and sorry for my English.




Here my parameters :

BlockMesh :

convertToMeters 0.001;

vertices
(
( 0 0 0)
(18 0 0)
( 0 1 -0.008726655121)
(18 1 -0.008726655121)
( 0 9 -0.017453310242)
(18 9 -0.017453310242)
( 0 1 0.008726655121)
(18 1 0.008726655121)
( 0 9 0.017453310242)
(18 9 0.017453310242)
);

blocks
(
hex (0 1 3 2 0 1 7 6) (60 10 1) simpleGrading (1 1 1)
hex (2 3 5 4 6 7 9 8) (60 10 1) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
inlet
{
type patch;
faces
(
(0 2 6 0)
);
}

outlet
{
type patch;
faces
(
(1 3 7 1)
(3 5 9 7)
);
}

freestreamInlet
{
type patch;
faces
(
(2 4 8 6)
);
}

freestream
{
type patch;
faces
(
(4 8 9 5)
);
}

wedge1
{
type wedge;
faces
(
(0 2 3 1)
(2 4 5 3)
);
}

wedge2
{
type wedge;
faces
(
(0 1 7 6)
(6 7 9 8)
);
}
);

mergePatchPairs
(
);





ThermoPhysicalProperties :

thermoType ePsiThermo<pureMixture<sutherlandTransport<specieT hermo<hConstThermo<perfectGas>>>>>;

mixture
{
specie
{
nMoles 1;
molWeight 39.948;
}
thermodynamics
{
Cp 520;
Hf 30;
}
transport
{
As 2.054e-06;
Ts 162.4;
}
}


p :


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

internalField uniform 101;

boundaryField
{
inlet
{
type fixedValue;
value uniform 271724;
}

outlet
{
type waveTransmissive;
field p;
phi phi;
rho rho;
psi psi;
fieldInf 101;
gamma 1.67;
lInf 0.025;
value uniform 101;
}

freestream
{
type totalPressure;
value uniform 101;
p0 uniform 101;
U U;
phi phi;
rho none;
psi psi;
gamma 1.67;
}

freestreamInlet
{
type zeroGradient;
}

wedge1 {type wedge;}
wedge2 {type wedge;}
}


T :

dimensions [0 0 0 1 0 0 0];

internalField uniform 98.0;

boundaryField
{
inlet
{
type fixedValue;
value uniform 299.1;
}

outlet
{
type zeroGradient;
}

freestream
{
type totalTemperature;
value uniform 97;
T0 uniform 97;
U U;
phi phi;
rho none;
psi psi;
gamma 1.67;
}

freestreamInlet
{
type fixedValue;
value uniform 97.0;
}

wedge1 {type wedge;}
wedge2 {type wedge;}
}


U :

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

internalField uniform (0 0 0);

boundaryField
{
inlet
{
type fixedValue;
value uniform (515.6 0 0);
}

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

freestream
{
type zeroGradient;
}

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

wedge1 {type wedge;}
wedge2 {type wedge;}
}