Problem with rhoSimpleFoam
 

I am trying to simulate a simple Plane configuration. M=0.9 & alpha-1.
No.Of cells : 459487 (tetrahedral).
Solver : rhoSimpleFoam

I attached the "system" folder.
when i run the simulation i receive the following error.

idea@IDeA:~/OpenFOAM/idea-1.6/run/Problems/RAEplane$ rhoSimpleFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.6-f802ff2d6c5a
Exec : rhoSimpleFoam
Date : Mar 16 2010
Time : 16:15:11
Host : IDeA
PID : 13488
Case : /home/idea/OpenFOAM/idea-1.6/run/Problems/RAEplane
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting RAS turbulence model laminar

Starting time loop

Time = 1

#0 Foam::error::printStack(Foam::Ostream&) in "/home/idea/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/idea/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib/libc.so.6"
#3 Foam::DILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/home/idea/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 Foam::DILUPreconditioner::DILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::dictionary const&) in "/home/idea/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#5 Foam::lduMatrix::preconditioner::addasymMatrixCons tructorToTable<Foam::DILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::dictionary const&) in "/home/idea/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#6 Foam::lduMatrix::preconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) in "/home/idea/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/idea/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#8 Foam::fvMatrix<Foam::Vector<double> >::solve(Foam::dictionary const&) in "/home/idea/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhoSimpleFoam"
#9 main in "/home/idea/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhoSimpleFoam"
#10 __libc_start_main in "/lib/libc.so.6"
#11 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
Floating point exception


Anyone please help me.

Regards
mecbe
:confused:

Hi mecbe,

A so early crash often indicates and error in the internal field values at the start.
Make sure to use absolute pressure and absolute temperature. Also the turbulence properties in the internal field should be non-zero.

Maybe this helps

stawrogin

Hi all,

where can I check if absolute or relative pressure is used?
Moreover, if I want to use a reference pressure of 101325 [Pa], where have I to set it?

You have to explore the code to see whether it uses absolute or relative pressure..

I am not sure about specifying free stream values..

Check this link. It may help
http://www.aero.polimi.it/freecase/?...:ONERA_M6_wing

Regards
mecbe