totalPressure in rhoSonicFoam
 

hello im trying to simulate a supersonic free jet . I tried to impose totalPressure boudary condition at the inlet :-

type totalPressure;
U U;
value uniform 0;
phi phiv;
rho rho;
psi psi;
gamma 1..4;
p0 uniform 1e5;


However i get the following error
#0 Foam::error::printStack(Foam::Ostream&) at OSspecific/Unix/printStack.C:203
#1 Foam::sigFpe::sigFpeHandler(int) at OSspecific/Unix/signals/sigFpe.C:130
#2 Uninterpreted: [0xffffe420]
#3 void VectorSpaceOps<3, 0>::opVS<Foam::Vector<double>, Foam::VectorSpace<Foam::Vector<dou
ble>, double, 3>, double, Foam::divideOp3<double, double, double> >(Foam::Vector<double>&, F
oam::VectorSpace<Foam::Vector<double>, double, 3> const&, double const&, Foam::divideOp3<dou
ble, double, double>) at ~/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude/VectorSpaceM.H:34
#4 Foam::Vector<double> Foam::operator/<Foam::Vector<double>, double, 3>(Foam::VectorSpace<
Foam::Vector<double>, double, 3> const&, double) at ~/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/l
nInclude/VectorSpaceI.H:512
#5 void Foam::divide<Foam::Vector<double> >(Foam::Field<Foam::Vector<double> >&, Foam::ULis
t<Foam::Vector<double> > const&, Foam::UList<double> const&) at ~/OpenFOAM/OpenFOAM-1.4.1/sr
c/OpenFOAM/lnInclude/FieldFunctions.C:509
#6 void Foam::divide<Foam::fvPatchField, Foam::Vector<double> >(Foam::FieldField<Foam::fvPa
tchField, Foam::Vector<double> >&, Foam::FieldField<Foam::fvPatchField, Foam::Vector<double>
> const&, Foam::FieldField<Foam::fvPatchField, double> const&) at ~/OpenFOAM/OpenFOAM-1.4.1
/src/OpenFOAM/lnInclude/FieldFieldFunctions.C:626
#7 void Foam::divide<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>(Foam::Geometr
icField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<Foam
::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Fo
am::fvPatchField, Foam::volMesh> const&) at ~/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude
/GeometricFieldFunctions.C:563
#8 Foam::tmp<Foam::GeometricField<Foam::Vector<double >, Foam::fvPatchField, Foam::volMesh>
> Foam::operator/<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>(Foam::GeometricFi
eld<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<do
uble, Foam::fvPatchField, Foam::volMesh> const&) at ~/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/l
nInclude/GeometricFieldFunctions.C:563
#9 main at ~/OpenFOAM/OpenFOAM-1.4.1/applications/solvers/compressible/rhoSonicFoam/createF
ields.H:96
#10 __libc_start_main in "/lib/libc.so.6"
#11 __gxx_personality_v0 in "/home/msamel/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGcc
DPDebug/rhoSonicFoam"

which obviously a FPE error due to division by 0. It mostly occurs while calculating rhoE . what can be the error in initializing the totalPressure BC in the above manner ???
with a fixedValue bc for pressure, the code runs absolutely perfect !!!!!!!!!!!