[headwing]

Hi,
I'm in the middle of a masters thesis concerning a biological CFD problem.

OpenFOAM is the third solver I've tried, and so far, seems promising.

However, when I try to run simpleFoam in my case directory, the following error results...

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.7.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 1.7.0-113391ee57bd
Exec   : simpleFoam
Date   : Jan 14 2011
Time   : 12:07:07
Host   : matt-desktop
PID    : 3763
Case   : /home/matt/OpenFOAM/matt-1.7.0/run/Smokarran-adult
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

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

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kOmegaSST
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so"
#1  Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so"
#2   in "/lib/libc.so.6"
#3  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::mag<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/openfoam170/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#4  Foam::incompressible::RASModels::kOmegaSST::kOmegaSST(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/opt/openfoam170/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#5  Foam::incompressible::RASModel::adddictionaryConstructorToTable<Foam::incompressible::RASModels::kOmegaSST>::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/opt/openfoam170/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#6  Foam::incompressible::RASModel::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/opt/openfoam170/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#7  
 in "/opt/openfoam170/applications/bin/linux64GccDPOpt/simpleFoam"
#8  __libc_start_main in "/lib/libc.so.6"
#9  
 in "/opt/openfoam170/applications/bin/linux64GccDPOpt/simpleFoam"
Floating point exception

Apparently this is a division-by-zero sort of error...

Does the answer seem obvious to anyone, or have any of you had this sort of thing before?

NOTE: I should probably mention that I'm using OpenFOAM 1.7.0 on CAELinux2010 (it came pre-compiled) and that the "2Dairfoil" tutorial seems to run fine (it completed 500 iterations without a hitch) when I invoke simpleFoam in that case directory.

Thanks for your help!

Matt

[nakul]

Hi,

This is a common problem. There is something wrong with your BCs or Initial Conditions. Check them again.