CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

dieselFoam/sprayFoam floating point error when initializing case

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 16, 2012, 16:31
Default dieselFoam/sprayFoam floating point error when initializing case
  #1
New Member
 
James
Join Date: Jan 2012
Posts: 11
Rep Power: 5
jmart7 is on a distinguished road
Hello everyone,

So I am trying to simulate a multiple spray injector case based on the aachenBomb for sprayFoam but using a different geometry, but when the case is being initialized I get a floating point exception after the kEpsilon RAS model is selected. Looking at the error, I can see that it has do do with calculating the kEpsilon coefficients when using mutkWallFunction in the boundaries.

I have tried dieselFoam and I get the same error. Here is the output, any help would be greatly appreciated! Thanks!

Create mesh for time = 0

Reading chemistry properties


Reading g
Reading thermophysical properties

Selecting psiChemistryModel ODEChemistryModel<gasThermoPhysics>
Selecting thermodynamics package hsPsiMixtureThermo<reactingMixture<gasThermoPhysic s>>
Selecting chemistryReader chemkinReader
ODEChemistryModel: Number of species = 5 and reactions = 1
Selecting ODE solver SIBS
Creating component thermo properties:
multi-component carrier - 5 species
liquids - 1 components
solids - 0 components

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
#0 Foam::error:rintStack(Foam::Ostream&) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 __restore_rt at sigaction.c:0
#3 Foam::compressible::RASModels::mutkWallFunctionFvP atchScalarField::calcMut() const in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#4 Foam::compressible::RASModels::mutkWallFunctionFvP atchScalarField::updateCoeffs() in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#5 Foam::fvPatchField<double>::evaluate(Foam::UPstrea m::commsTypes) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/sprayFoam"
#6 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::evaluate() in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/sprayFoam"
#7 Foam::compressible::RASModels::kEpsilon::kEpsilon( Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&, Foam::word const&, Foam::word const&) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#8 Foam::compressible::RASModel::adddictionaryConstru ctorToTable<Foam::compressible::RASModels::kEpsilo n>::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&, Foam::word const&) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#9 Foam::compressible::RASModel::New(Foam::GeometricF ield<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&, Foam::word const&) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#10 Foam::compressible::turbulenceModel::addturbulence ModelConstructorToTable<Foam::compressible::RASMod el>::NewturbulenceModel(Foam::GeometricField<doubl e, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&, Foam::word const&) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#11 Foam::compressible::turbulenceModel::New(Foam::Geo metricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&, Foam::word const&) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so"
#12 main in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/sprayFoam"
#13 __libc_start_main in "/lib64/tls/libc.so.6"
#14 Foam::UOPstream::write(char) in "/share/home/01803/jaime/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/sprayFoam"
jmart7 is offline   Reply With Quote

Old   August 7, 2012, 08:04
Default
  #2
Member
 
Adam
Join Date: Jun 2011
Posts: 32
Rep Power: 6
Smed is on a distinguished road
Any luck resolving this? I'm getting a similar error.
Smed is offline   Reply With Quote

Old   September 16, 2013, 07:54
Default
  #3
Member
 
Frank Ubber
Join Date: Aug 2013
Posts: 32
Rep Power: 4
kornickel is on a distinguished road
---POST DELETED---

Sorry, wrong thread
kornickel is offline   Reply With Quote

Old   September 16, 2013, 08:26
Unhappy
  #4
Senior Member
 
dkxls's Avatar
 
Armin
Join Date: Feb 2011
Location: Helsinki, Finland
Posts: 152
Rep Power: 10
dkxls will become famous soon enough
sprayFoam has basically been unusable ever since it was introduced in 2.0.x.

Even if in some cases it would have produced some results (i.e. didn't crash), these results were most likely wrong - at least for engine spray simulation!

I fixed most of the issues concerning the thermodynamics and contribute them to upstream OpenFOAM (2.2.x), but there are still several issues, especially concerning the breakup models.
For more information, see the list of bugs I reported/fixed in the last couple of weeks (and there are still more to come...):
http://www.openfoam.org/mantisbt/vie...eporter_id=294

Also have a look here:
spray droplet distortion parameter

To get started with sprayFoam, get a recent git checkout, set the atomization model to "none" (it's broken altogether right now) and use maybe the ReitzDiwakar or PilchErdman breakup model (at least don't use ReitzKHRT or TAB/ETAB).

I would very much appreciate more people testing this stuff...
dkxls is offline   Reply With Quote

Old   September 16, 2013, 08:34
Default
  #5
Member
 
Frank Ubber
Join Date: Aug 2013
Posts: 32
Rep Power: 4
kornickel is on a distinguished road
---POST DELETED---

Sorry, wrong thread
kornickel is offline   Reply With Quote

Old   September 16, 2013, 08:41
Default
  #6
Senior Member
 
dkxls's Avatar
 
Armin
Join Date: Feb 2011
Location: Helsinki, Finland
Posts: 152
Rep Power: 10
dkxls will become famous soon enough
Oho!
I just read the header and didn't exactly read your message, so yes you should probably start a new thread.

And then again, my post is not related to RANS turbulence models at all, so ... also questionable if it really belongs here.
dkxls is offline   Reply With Quote

Reply

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
MPI Error - simpleFoam - Floating Point Exception scott OpenFOAM Running, Solving & CFD 3 April 13, 2012 16:34
Floating point exception error Alan OpenFOAM Running, Solving & CFD 10 April 6, 2012 14:02
floating point exception? syd STAR-CCM+ 3 April 20, 2010 18:55
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 12:24
Floating point error Sohail Ahmed FLUENT 0 March 12, 2004 01:02


All times are GMT -4. The time now is 04:34.