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

wired error in sprayFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 13, 2012, 21:11
Default wired error in sprayFoam
  #1
Member
 
Join Date: Jun 2012
Posts: 65
Rep Power: 5
conceptone is on a distinguished road
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 !
Attached Files
File Type: docx sprayCloudProperties_copy.docx (26.5 KB, 27 views)
conceptone is offline   Reply With Quote

Old   June 16, 2012, 05:19
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,251
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
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":
Quote:
Code:
[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: how to add a new gas

Last edited by wyldckat; June 16, 2012 at 07:44. Reason: see "PS:"
wyldckat is offline   Reply With Quote

Old   June 16, 2012, 14:27
Default
  #3
Member
 
Join Date: Jun 2012
Posts: 65
Rep Power: 5
conceptone is on a distinguished road
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 is offline   Reply With Quote

Old   June 16, 2012, 17:27
Default
  #4
Member
 
Join Date: Jun 2012
Posts: 65
Rep Power: 5
conceptone is on a distinguished road
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.
conceptone is offline   Reply With Quote

Old   June 16, 2012, 17:31
Default
  #5
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,251
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
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: Debug version of OpenFOAM-1.6 - and yes, you'll only need to build the debug version on your personal installation!
wyldckat is offline   Reply With Quote

Old   June 16, 2012, 22:18
Default
  #6
Member
 
Join Date: Jun 2012
Posts: 65
Rep Power: 5
conceptone is on a distinguished road
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 is offline   Reply With Quote

Old   June 16, 2012, 22:29
Default
  #7
Member
 
Join Date: Jun 2012
Posts: 65
Rep Power: 5
conceptone is on a distinguished road
Luckly, I found the correct version by myself

XieXie!
conceptone is offline   Reply With Quote

Old   June 20, 2012, 22:10
Default
  #8
Member
 
Join Date: Jun 2012
Posts: 65
Rep Power: 5
conceptone is on a distinguished road
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).
------------------------------------------
conceptone is offline   Reply With Quote

Old   June 21, 2012, 16:41
Default
  #9
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,251
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Seriously, stop posting on two threads the same question ... You've already been answered here: how to add a new gas #13
wyldckat 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
Multiple Injection points in sprayFoam kripalpariyangat OpenFOAM Pre-Processing 11 June 29, 2015 13:25
aachenBomb with sprayFOAM not working jenzkeller OpenFOAM Bugs 8 June 8, 2014 18:01
sprayFoam crashes lukasfischer OpenFOAM Running, Solving & CFD 3 July 14, 2013 11:08
Standard Drag Model does not exist for sprayFoam arash1 OpenFOAM Running, Solving & CFD 0 April 8, 2012 17:52
sprayFoam results depend on the type of chemistryReader (chemkinReader/foamChemistry arash1 OpenFOAM Running, Solving & CFD 1 April 2, 2012 12:47


All times are GMT -4. The time now is 02:44.