[vianne]

Dear all,

I have an urgent problem and I hope someone can help
I'm trying to simulate the reactingFoam-tutorial from OpenWiki, but I get the following error message running reactingFoam (blockMesh is okay). I use exactly the files uploaded, and I also tried my own case, but I always get the same error message. I would appreciate any suggestions...

Reading chemistry properties
Reading g
Reading thermophysicalProperties
Selecting psiChemistryModel ODEChemistryModel<gasThermoPhysics>
Selecting thermodynamics package hPsiMixtureThermo<reactingMixture<gasThermoPhysics >>
Selecting chemistryReader chemkinReader
Selecting chemistrySolver ode
Selecting ODE solver SIBS
ODEChemistryModel: Number of species = 5 and reactions = 1
Reading field U
Reading/calculating face flux field phi
Creating turbulence model.
Selecting turbulence model type laminar
Creating field DpDt
Courant Number mean: 0.0388322 max: 0.378462
Starting time loop
Courant Number mean: 0.010219 max: 0.0995953
deltaT = 2.63158e-05
Time = 2.63158e-05
Solving chemistry
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/steffi/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/steffi/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/steffi/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/steffi/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/reactingFoam"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/steffi/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/reactingFoam"
#6 main in "/home/steffi/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/reactingFoam"
#7 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#8 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/i386/elf/start.S:122
Floating point exception

[vianne]

Oh, I just forgot to tell you: I use Ubuntu Linux 9.10, OpenFOAM 1.6, Qt 4.5.2. And I really use the exact files without any changes given in the OpenWiki tutorial for reactingFoam. Has anyone else perhaps tried out this tutorial and found out, why this error is appearing?
When I use unset FOAM_SIGFPE, the simulation is running, but there's only mixing and the reaction doesn't take place...

[kalle]

Hi,

First, the reactingFoam tutorial is a bit of a strange set up. It tries to simulate heptane burning with pure oxygen. Preferrably, the setting for the oxidizer inlet should be 0.23 O2 and 0.77 N2. I recall people having trouble with this before. Now, this error is maybe not related to that.

Didn't you make any changes to the tutorial? When I try it on my OF1.5 I need to change a header of a file (well, not that important), add 'RASproperties', and fix the link to the chemkin-files.

Regards,
Kalle