CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

totalPressure in rhoSonicFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   April 16, 2009, 15:41
Default totalPressure in rhoSonicFoam
Join Date: Mar 2009
Posts: 40
Rep Power: 10
mihir1310 is on a distinguished road
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:rintStack(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>:pVS<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:perator/<Foam::Vector<double>, double, 3>(Foam::VectorSpace<
Foam::Vector<double>, double, 3> const&, double) at ~/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/l
#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
#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
#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
#8 Foam::tmp<Foam::GeometricField<Foam::Vector<double >, Foam::fvPatchField, Foam::volMesh>
> Foam:perator/<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
#9 main at ~/OpenFOAM/OpenFOAM-1.4.1/applications/solvers/compressible/rhoSonicFoam/createF
#10 __libc_start_main in "/lib/"
#11 __gxx_personality_v0 in "/home/msamel/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGcc

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 !!!!!!!!!!!
mihir1310 is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem when running rhoSonicFoam solver gaottino OpenFOAM Running, Solving & CFD 3 March 1, 2011 05:59
variable Time step in rhoSonicFoam mihir1310 OpenFOAM Running, Solving & CFD 2 March 31, 2009 11:22
TotalPressure %3d data to specify floooo OpenFOAM Running, Solving & CFD 0 December 8, 2008 04:58
RhoSonicFoam srinath OpenFOAM Running, Solving & CFD 0 May 30, 2008 01:21
RhoSOnicFoam shuo OpenFOAM Running, Solving & CFD 2 November 6, 2006 01:25

All times are GMT -4. The time now is 14:06.