[usergk]

Hello,

I am trying to implement combustion in OpenFOAM and I get the error below. This happens mid-way during the simulation, after a few time steps.

Any idea why this could be occurring? I am relatively new to using OpenFOAM, so any information would be appreciated.

Thanks!
gk

[24] #0 Foam::error:rintStack(Foam::Ostream&) in "/home/gk/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
[24] #1 Foam::sigFpe::sigFpeHandler(int) in "/home/gk/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
[24] #2 in "/lib/libc.so.6"
[24] #3 in "/lib/libm.so.6"
[24] #4 pow in "/lib/libm.so.6"
[24] #5 Foam::ODEChemistryModel<Foam:siChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam: :janafThermo<Foam:erfectGas> > > >:mega(Foam::Reaction<Foam::sutherlandTransport< Foam::specieThermo<Foam::janafThermo<Foam:erfect Gas> > > > const&, Foam::Field<double> const&, double, double, double&, double&, int&, double&, double&, int&) const in "/home/gk/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libchemistryModel.so"
[24] #6 Foam::ODEChemistryModel<Foam:siChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam: :janafThermo<Foam:erfectGas> > > >::tc() const in "/home/gk/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libchemistryModel.so"
[24] #7
[24] in "/home/gk/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/pFoam"
[24] #8 __libc_start_main in "/lib/libc.so.6"
[24] #9
[24] in "/home/gk/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/pFoam"
[node69:00818] *** Process received signal ***
[node69:00818] Signal: Floating point exception (8)
[node69:00818] Signal code: (-6)
[node69:00818] Failing at address: 0x58f800000332
[node69:00818] [ 0] /lib/libc.so.6(+0x33af0) [0x2b0cdf0abaf0]
[node69:00818] [ 1] /lib/libc.so.6(gsignal+0x35) [0x2b0cdf0aba75]
[node69:00818] [ 2] /lib/libc.so.6(+0x33af0) [0x2b0cdf0abaf0]
[node69:00818] [ 3] /lib/libm.so.6(+0x13e81) [0x2b0cdebf0e81]
[node69:00818] [ 4] /lib/libm.so.6(pow+0x15) [0x2b0cdec02765]
[node69:00818] [ 5] /home/gk/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libchemistryModel.so(_ZNK4Foam17ODEChemistryModelI NS_17psiChemistryModelENS_19sutherlandTransportINS _12specieThermoINS_11janafThermoINS_10perfectGasEE EEEEEE5omegaERKNS_8ReactionIS8_EERKNS_5FieldIdEEdd RdSI_RiSI_SI_SJ_+0x285) [0x2b0cdd978ff5]
[node69:00818] [ 6] /home/gk/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libchemistryModel.so(_ZNK4Foam17ODEChemistryModelI NS_17psiChemistryModelENS_19sutherlandTransportINS _12specieThermoINS_11janafThermoINS_10perfectGasEE EEEEEE2tcEv+0x57e) [0x2b0cdd98424e]
[node69:00818] [ 7] pFoam() [0x426bf3]
[node69:00818] [ 8] /lib/libc.so.6(__libc_start_main+0xfd) [0x2b0cdf096c4d]
[node69:00818] [ 9] pFoam() [0x421119]
[node69:00818] *** End of error message ***

[usergk]

Can anyone explain why this floating point error may be occurring?

Thanks,
gk

[jabecker]

Hello usergk,

Did you ever figure out the issue? Your help would be greatly appreciated. It seems like the error is generally caused by having an initial condition set wrong, but I cannot seem to figure out why mine is not correct. When I have things set as zero, I want them to be set as zero. Mine comes from a slightly different thermophysicalProperties folder, but I would imagine it is caused from the same problem. My error reads:

Selecting thermodynamics package hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
#0 Foam::error:rintStack(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::hPsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::calculate() in "/opt/openfoam170/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#4 Foam::hPsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::hPsiThermo(Foam::fvMesh const&) in "/opt/openfoam170/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#5 Foam::basicPsiThermo::addfvMeshConstructorToTable< Foam::hPsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > > >::New(Foam::fvMesh const&) in "/opt/openfoam170/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#6 Foam::basicPsiThermo::New(Foam::fvMesh const&) in "/opt/openfoam170/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#7
in "/opt/openfoam170/applications/bin/linux64GccDPOpt/chtMultiRegionFoam"
#8 __libc_start_main in "/lib/libc.so.6"
#9
in "/opt/openfoam170/applications/bin/linux64GccDPOpt/chtMultiRegionFoam"
Floating point exception

I really hope you can help! Thanks in advance,
Jen

[usergk]

Hi Jen,

Some of the possible ways:
Initial conditions be more accurate
Mesh may need to be improved
Change fvSchemes/fvSolution files - There are some posts on the forum that talk about this (esp when k/epsilon blows up)

gk

[jabecker]

Hello gk,

Thank you for your quick response! I have another question: Do you work with chtMultiRegionFoam? This is the program I am currently working with and it is giving me quite the headache. I found out that the issue was indeed a zero condition in the initial folder that should not be set as zero, but this is not what I had set it as. My initial conditions for my topAir get completely changed when I use ./Allrun and this is where the error comes from. The only way I can get the program to get rid of this error is by running ./Allrun and then completely changing my 0/topAir files all over again. It seems really strange, but I cannot figure out why it is replacing all of the folders.

Your help again would be greatly appreciated. Thanks,
Jen

[usergk]

Sorry, I don't use this particular solver. You can try searching on the forum for similar issues (assuming you haven't already).

Best,
gk