[conceptone]

hi,everyone,i am a new openfoam user. I am running the sprayFoam case, where i implement the specie C12H26. I have changed the chem.inp and therm.dat. But when the case runs a few steps, it cracks. It stopped at 0.00151s.The error is like below:

Courant Number mean: 3.71158e-05 max: 0.0695979
deltaT = 1.25e-06
Time = 0.001515

Solving cloud sprayCloud

--> Cloud: sprayCloud
Added 6 new parcels

[6] #0 Foam::error:rintStack(Foam::Ostream&) in "/usr/apps1/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64Gcc44DPOpt/lib/libOpenFOAM.so"
[6] #1 Foam::sigFpe::sigHandler(int) in "/usr/apps1/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64Gcc44DPOpt/lib/libOpenFOAM.so"
[6] #2 __restore_rt at sigaction.c:0
[6] #3 __ieee754_exp at interp.c:0
[6] #4 exp in "/lib64/libm.so.6"
[6] #5 Foam::LiquidEvaporation<Foam::ReactingCloud<Foam:: ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam: :SprayParcel<Foam::ReactingParcel<Foam::ThermoParc el<Foam::KinematicParcel<Foam:article> > > > > > > > >::calculate(double, int, double, double, double, double, double, double, Foam::Field<double>&) const in "/usr/apps1/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64Gcc44DPOpt/lib/liblagrangianSpray.so"
[6] #6 void Foam::ReactingParcel<Foam::ThermoParcel<Foam::Kine maticParcel<Foam:article> > >::calc<Foam::ReactingParcel<Foam::ThermoParcel<Fo am::KinematicParcel<Foam:article> > >::TrackingData<Foam::SprayCloud<Foam::ReactingClo ud<Foam::ThermoCloud<Foam::KinematicCloud<Foam::Cl oud<Foam::SprayParcel<Foam::ReactingParcel<Foam::T hermoParcel<Foam::KinematicParcel<Foam:article> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::Ki nematicParcel<Foam:article> > >::TrackingData<Foam::SprayCloud<Foam::ReactingClo ud<Foam::ThermoCloud<Foam::KinematicCloud<Foam::Cl oud<Foam::SprayParcel<Foam::ReactingParcel<Foam::T hermoParcel<Foam::KinematicParcel<Foam:article> > > > > > > > > >&, double, int) in "/usr/apps1/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64Gcc44DPOpt/bin/sprayFoam"

......

Is there anybody who can help me? I guess it should be caused by particle or parcel parameters set in spraycloudProperties file, so I attached it. It can be open by word.
Thanks !

[wyldckat]

Greetings Conceptone,

An M$ Office Word file? Are you serious? Why didn't you attach the normal "spraycloudProperties" file? Perhaps packaged in zip format?
Lucky you that LibreOffice can open this blasphemous file format

Anyway, ranting aside: the error indicates that something went wrong when doing an exponential calculation, namely using "exp":

[6] #1  Foam::sigFpe::sigHandler(int)
...
[6] #4  exp in "/lib64/libm.so.6"

The "sigFpe" is a signal function for Floating Point Exception: http://en.wikipedia.org/wiki/SIGFPE - basically because "exp" tried to do something it shouldn't try to do.

The solution - use divide-and-conquer:

1. Since you added a new model, you should validate it first. This is also a good exercise for you, so go search in the "OpenFOAM-*/applications" folder for applications that use chemistry and see if there is any that uses the chemistry set more explicitly. I think the "utilities" and "test" sub-folders are the best place to look.
2. After you've validated that the equations and values were properly added to the "chem.inp" and "therm.dat" files, try a simpler version of this case you're running. Use a really smaller profile for initial particles.
3. If this fails, you'll have to build the debug version of OpenFOAM, so you can get closer to where "exp" gave you the problem. Here you can also try the "gdbOF" plugin: http://openfoamwiki.net/index.php/Contrib_gdbOF

Good luck!
Bruno


PS: to close the loop on this thread, it looks like this was initially being followed here: http://www.cfd-online.com/Forums/ope...d-new-gas.html

[conceptone]

Oh man, thanks your kind help! I'll check these things again..I'm just a new user coming through the period..

PSerhaps we should have a dinner together

of course I'm kidding!

[conceptone]

Oh by the way is it true that if I want to run the debug version the openfoam should be installed on my own computer? I mean my openfoam is installed on the service computer, so I can't modify the bashrc file since I don't have the managing right.

[wyldckat]

You can easily do your own personal build of OpenFOAM - just follow the instructions: http://www.openfoam.org/download/source.php

As for the debug build, here's an old thread on this subject: http://www.cfd-online.com/Forums/ope...oam-1-6-a.html - and yes, you'll only need to build the debug version on your personal installation!

[conceptone]

Thanks man!I tried the source pack from http://www.openfoam.org/download/source.php

But it seems it is only for version 2.1.1 while mine is 2.1.0 because when I type in source command "source $HOME/.cshrc" I got the error message from Linux:
elephant 166% source $HOME/.cshrc
/usr/apps1/openfoam-2.1.0/OpenFOAM-2.1.1/bin/foamEtcFile: Command not found.
/usr/apps1/openfoam-2.1.0/OpenFOAM-2.1.1/bin/foamCleanPath: Command not found.
/usr/apps1/openfoam-2.1.0/OpenFOAM-2.1.1/bin/foamCleanPath: Command not found.
/usr/apps1/openfoam-2.1.0/OpenFOAM-2.1.1/bin/foamCleanPath: Command not found.
/usr/apps1/openfoam-2.1.0/OpenFOAM-2.1.1/etc/config/settings.csh: No such file or directory

So is there any version for 2.1.0?

[conceptone]

Luckly, I found the correct version by myself

XieXie!

[conceptone]

Hi Bruno, the linux system manager has helped me to build the debug version(it seemed that it's not gdb) so I am not sure where should I find the so-called detail debug information. So I copy the message from log.sprayFoam file, can you help me about the error? Thanks.

[6] #0 Foam::error:rintStack(Foam::Ostream&) at /usr/apps1/openfoam-2.1.1/OpenFOAM-2.1.1/src/OSspecific/POSIX/printStack.C:201
[6] #1 Foam::sigFpe::sigHandler(int) at /usr/apps1/openfoam-2.1.1/OpenFOAM-2.1.1/src/OSspecific/POSIX/signals/sigFpe.C:117
[6] #2 __restore_rt at sigaction.c:0
[6] #3 __ieee754_exp at interp.c:0
[6] #4 exp in "/lib64/libm.so.6"
[6] #5 Foam::exp(double) at /usr/apps1/openfoam-2.1.1/OpenFOAM-2.1.1/src/OpenFOAM/lnInclude/Scalar.H:260
[6] #6 Foam::NSRDSfunc1::f(double, double) const at /usr/apps1/openfoam-2.1.1/OpenFOAM-2.1.1/src/thermophysicalModels/thermophysicalFunctions/lnInclude/NSRDSfunc1.H:108
[6] #7 Foam::C12H26:v(double, double) const at /usr/apps1/openfoam-2.1.1/OpenFOAM-2.1.1/src/thermophysicalModels/properties/liquidProperties/C12H26/C12H26I.H:34
[6] #8 Foam::LiquidEvaporation<Foam::ReactingCloud<Foam:: ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam: :SprayParcel<Foam::ReactingParcel<Foam::ThermoParc el<Foam::KinematicParcel<Foam:article> > > > > > > > >::calculate(double, int, double, double, double, double, double, double, double, double, Foam::Field<double> const&, Foam::Field<double>&) const at /usr/apps1/openfoam-2.1.1/OpenFOAM-2.1.1/src/lagrangian/intermediate/lnInclude/LiquidEvaporation.C:172
[6] #9 void Foam::ReactingParcel<Foam::ThermoParcel<Foam::Kine maticParcel<Foam:article> > >::calcPhaseChange<Foam::ReactingParcel<Foam::Ther moParcel<Foam::KinematicParcel<Foam:article> > >::TrackingData<Foam::SprayCloud<Foam::ReactingClo ud<Foam::ThermoCloud<Foam::KinematicCloud<Foam::Cl oud<Foam::SprayParcel<Foam::ReactingParcel<Foam::T hermoParcel<Foam::KinematicParcel<Foam:article> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::Ki nematicParcel<Foam:article> > >::TrackingData<Foam::SprayCloud<Foam::ReactingClo ud<Foam::ThermoCloud<Foam::KinematicCloud<Foam::Cl oud<Foam::SprayParcel<Foam::ReactingParcel<Foam::T hermoParcel<Foam::KinematicParcel<Foam:article> > > > > > > > > >&, double, int, double, double, double, double, double, double, double, int, double, Foam::Field<double> const&, Foam::Field<double>&, double&, double&, double&, Foam::Field<double>&) at /usr/apps1/openfoam-2.1.1/OpenFOAM-2.1.1/src/lagrangian/intermediate/lnInclude/ReactingParcel.C:517
...
[pn1:29034] [21] sprayFoam [0x472099]
[pn1:29034] *** End of error message ***
--------------------------------------------------------------------------
mpiexec noticed that process rank 6 with PID 29034 on node pn1 exited on signal 8 (Floating point exception).
------------------------------------------

[wyldckat]

Seriously, stop posting on two threads the same question ... You've already been answered here: http://www.cfd-online.com/Forums/ope...tml#post367702 #13