[BakedAlmonds]

Hi,

My main aim of using reactingFoam was to do a pollution study.
Therefore, I tried with the reactingFoam tutorial but faced an error.

What I did:
- from chemistryProperties: turn off chemistry
- from combustionProperties: change active to false
- from reactions: delete reactions equations
- from 0/CH4 file: change fuel value to uniform 0

reactingFoam ran well until the last step. When I change to an input without CH4, it starts to show this error.

Selecting chemistryReader foamChemistryReader
  #0  Foam: : error: : printStack(Foam: : Ostream&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/pla
  tforms/linux64Gcc48DPInt32Opt/lib/libOpenFOAM.so"
  #1  Foam: : sigFpe: : sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/lin
  ux64Gcc48DPInt32Opt/lib/libOpenFOAM.so"
  #2  ? in "/lib64/libc.so.6"
  #3  Foam: : heThermo<Foam: : psiReactionThermo, Foam: : SpecieMixture<Foam: : reactingMi
  xture<Foam: : sutherlandTransport<Foam: : species: : thermo<Foam: : janafThermo<Foam: : pe
  rfectGas<Foam: : specie> >, Foam: : sensibleEnthalpy> > > > >: : he(Foam: : Field<double
  > const&, Foam: : Field<double> const&, int) const in "/opt/OpenFOAM/OpenFOAM-v3.0
  +/platforms/linux64Gcc48DPInt32Opt/lib/libreactionThermophysicalModels.so"
  #4  Foam: : heThermo<Foam: : psiReactionThermo, Foam: : SpecieMixture<Foam: : reactingMi
  xture<Foam: : sutherlandTransport<Foam: : species: : thermo<Foam: : janafThermo<Foam: : pe
  rfectGas<Foam: : specie> >, Foam: : sensibleEnthalpy> > > > >: : heThermo(Foam: : fvMesh
   const&, Foam: : word const&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gc
  c48DPInt32Opt/lib/libreactionThermophysicalModels.so"
  #5  Foam: : psiReactionThermo: : addfvMeshConstructorToTable<Foam: : hePsiThermo<Foam: 
  : psiReactionThermo, Foam: : SpecieMixture<Foam: : reactingMixture<Foam: : sutherlandTr
  ansport<Foam: : species: : thermo<Foam: : janafThermo<Foam: : perfectGas<Foam: : specie> >
  , Foam: : sensibleEnthalpy> > > > > >: : New(Foam: : fvMesh const&, Foam: : word const&)
   in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libreacti
  onThermophysicalModels.so"
  #6  Foam: : autoPtr<Foam: : psiReactionThermo> Foam: : basicThermo: : New<Foam: : psiReact
  ionThermo>(Foam: : fvMesh const&, Foam: : word const&) in "/opt/OpenFOAM/OpenFOAM-v3
  .0+/platforms/linux64Gcc48DPInt32Opt/lib/libreactionThermophysicalModels.so"
  #7  Foam: : psiReactionThermo: : New(Foam: : fvMesh const&, Foam: : word const&) in "/op
  t/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libreactionThermo
  physicalModels.so"
  #8  Foam: : psiChemistryModel: : psiChemistryModel(Foam: : fvMesh const&, Foam: : word c
  onst&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/lib
  chemistryModel.so"
  #9  Foam: : chemistryModel<Foam: : psiChemistryModel, Foam: : sutherlandTransport<Foam
  : : species: : thermo<Foam: : janafThermo<Foam: : perfectGas<Foam: : specie> >, Foam: : sens
  ibleEnthalpy> > >: : chemistryModel(Foam: : fvMesh const&, Foam: : word const&) in "/o
  pt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libchemistryMode
  l.so"
  #10  Foam: : EulerImplicit<Foam: : chemistryModel<Foam: : psiChemistryModel, Foam: : sut
  herlandTransport<Foam: : species: : thermo<Foam: : janafThermo<Foam: : perfectGas<Foam: : 
  specie> >, Foam: : sensibleEnthalpy> > > >: : EulerImplicit(Foam: : fvMesh const&, Foa
  m: : word const&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Op
  t/lib/libchemistryModel.so"
  #11  Foam: : psiChemistryModel: : addfvMeshConstructorToTable<Foam: : EulerImplicit<Fo
  am: : chemistryModel<Foam: : psiChemistryModel, Foam: : sutherlandTransport<Foam: : spec
  ies: : thermo<Foam: : janafThermo<Foam: : perfectGas<Foam: : specie> >, Foam: : sensibleEn
  thalpy> > > > >: : New(Foam: : fvMesh const&, Foam: : word const&) in "/opt/OpenFOAM/O
  penFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libchemistryModel.so"
  #12  Foam: : autoPtr<Foam: : psiChemistryModel> Foam: : basicChemistryModel: : New<Foam: 
  : psiChemistryModel>(Foam: : fvMesh const&, Foam: : word const&) in "/opt/OpenFOAM/Op
  enFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libchemistryModel.so"
  #13  Foam: : psiChemistryModel: : New(Foam: : fvMesh const&, Foam: : word const&) in "/o
  pt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libchemistryMode
  l.so"
  #14  Foam: : combustionModels: : psiChemistryCombustion: : psiChemistryCombustion(Foam
  : : word const&, Foam: : fvMesh const&, Foam: : word const&) in "/opt/OpenFOAM/OpenFOA
  M-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libcombustionModels.so"
  #15  Foam: : combustionModels: : laminar<Foam: : combustionModels: : psiChemistryCombust
  ion>: : laminar(Foam: : word const&, Foam: : fvMesh const&, Foam: : word const&) in "/op
  t/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libcombustionMode
  ls.so"
  #16  Foam: : combustionModels: : psiCombustionModel: : adddictionaryConstructorToTable
  <Foam: : combustionModels: : laminar<Foam: : combustionModels: : psiChemistryCombustion>
   >: : New(Foam: : word const&, Foam: : fvMesh const&, Foam: : word const&) in "/opt/Open
  FOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libcombustionModels.so"
   
  #17  Foam: : combustionModels: : psiCombustionModel: : New(Foam: : fvMesh const&, Foam: : 
  word const&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/l
  ib/libcombustionModels.so"
  #18  ? in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/rea
  ctingFoam"
  #19  __libc_start_main in "/lib64/libc.so.6"
  #20  ? in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/rea
  ctingFoam"
  Floating point exception
  [ofuser@boot2docker aaaa]$ reactingFoam > log.reacting
  #0  Foam: : error: : printStack(Foam: : Ostream&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/pla
  tforms/linux64Gcc48DPInt32Opt/lib/libOpenFOAM.so"
  #1  Foam: : sigFpe: : sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/lin
  ux64Gcc48DPInt32Opt/lib/libOpenFOAM.so"
  #2  ? in "/lib64/libc.so.6"
  #3  Foam: : heThermo<Foam: : psiReactionThermo, Foam: : SpecieMixture<Foam: : reactingMi
  xture<Foam: : sutherlandTransport<Foam: : species: : thermo<Foam: : janafThermo<Foam: : pe
  rfectGas<Foam: : specie> >, Foam: : sensibleEnthalpy> > > > >: : he(Foam: : Field<double
  > const&, Foam: : Field<double> const&, int) const in "/opt/OpenFOAM/OpenFOAM-v3.0
  +/platforms/linux64Gcc48DPInt32Opt/lib/libreactionThermophysicalModels.so"
  #4  Foam: : heThermo<Foam: : psiReactionThermo, Foam: : SpecieMixture<Foam: :  reactingMi
  xture<Foam: : sutherlandTransport<Foam: : species: : thermo<Foam: : janafThermo<Foam: : pe
  rfectGas<Foam: : specie> >, Foam: : sensibleEnthalpy> > > > >: : heThermo(Foam: : fvMesh
   const&, Foam: : word const&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gc
  c48DPInt32Opt/lib/libreactionThermophysicalModels.so"
  #5  Foam: : psiReactionThermo: : addfvMeshConstructorToTable<Foam: : hePsiThermo<Foam: 
  : psiReactionThermo, Foam: : SpecieMixture<Foam: : reactingMixture<Foam: : sutherlandTr
  ansport<Foam: : species: : thermo<Foam: : janafThermo<Foam: : perfectGas<Foam: : specie> >
  , Foam: : sensibleEnthalpy> > > > > >: : New(Foam: : fvMesh const&, Foam: : word const&)
   in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libreacti
  onThermophysicalModels.so"
  #6  Foam: : autoPtr<Foam: : psiReactionThermo> Foam: : basicThermo: : New<Foam: : psiReact
  ionThermo>(Foam: : fvMesh const&, Foam: : word const&) in "/opt/OpenFOAM/OpenFOAM-v3
  .0+/platforms/linux64Gcc48DPInt32Opt/lib/libreactionThermophysicalModels.so"
  #7  Foam: : psiReactionThermo: : New(Foam: : fvMesh const&, Foam: : word const&) in "/op
  t/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libreactionThermo
  physicalModels.so"
  #8  Foam: : psiChemistryModel: : psiChemistryModel(Foam: : fvMesh const&, Foam: : word c
  onst&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/lib
  chemistryModel.so"
  #9  Foam: : chemistryModel<Foam: : psiChemistryModel, Foam: : sutherlandTransport<Foam
  : : species: : thermo<Foam: : janafThermo<Foam: : perfectGas<Foam: : specie> >, Foam: : sens
  ibleEnthalpy> > >: : chemistryModel(Foam: : fvMesh const&, Foam: : word const&) in "/o
  pt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libchemistryMode
  l.so"
  #10  Foam: : EulerImplicit<Foam: : chemistryModel<Foam: : psiChemistryModel, Foam: : sut
  herlandTransport<Foam: : species: : thermo<Foam: : janafThermo<Foam: : perfectGas<Foam: : 
  specie> >, Foam: : sensibleEnthalpy> > > >: : EulerImplicit(Foam: : fvMesh const&, Foa
  m: : word const&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Op
  t/lib/libchemistryModel.so"
  #11  Foam: : psiChemistryModel: : addfvMeshConstructorToTable<Foam: : EulerImplicit<Fo
  am: : chemistryModel<Foam: : psiChemistryModel, Foam: : sutherlandTransport<Foam: : spec
  ies: : thermo<Foam: : janafThermo<Foam: : perfectGas<Foam: : specie> >, Foam: : sensibleEn
  thalpy> > > > >: : New(Foam: : fvMesh const&, Foam: : word const&) in "/opt/OpenFOAM/O
  penFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libchemistryModel.so"
  #12  Foam: : autoPtr<Foam: : psiChemistryModel> Foam: : basicChemistryModel: : New<Foam: 
  : psiChemistryModel>(Foam: : fvMesh const&, Foam: : word const&) in "/opt/OpenFOAM/Op
  enFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libchemistryModel.so"
  #13  Foam: : psiChemistryModel: : New(Foam: : fvMesh const&, Foam: : word const&) in "/o
  pt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libchemistryMode
  l.so"
  #14  Foam: : combustionModels: : psiChemistryCombustion: : psiChemistryCombustion(Foam
  : : word const&, Foam: : fvMesh const&, Foam: : word const&) in "/opt/OpenFOAM/OpenFOA
  M-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libcombustionModels.so"
  #15  Foam: : combustionModels: : laminar<Foam: : combustionModels: : psiChemistryCombust
  ion>: : laminar(Foam: : word const&, Foam: : fvMesh const&, Foam: : word const&) in "/op
  t/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libcombustionMode
  ls.so"
  #16  Foam: : combustionModels: : psiCombustionModel: : adddictionaryConstructorToTable
  <Foam: : combustionModels: : laminar<Foam: : combustionModels: : psiChemistryCombustion>
   >: : New(Foam: : word const&, Foam: : fvMesh const&, Foam: : word const&) in "/opt/Open
  FOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libcombustionModels.so"
   
  #17  Foam: : combustionModels: : psiCombustionModel: : New(Foam: : fvMesh const&, Foam: : 
  word const&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/l
  ib/libcombustionModels.so"
  #18  ? in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/rea
  ctingFoam"
  #19  __libc_start_main in "/lib64/libc.so.6"
  #20  ? in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/rea
  ctingFoam"

Can anyone assist me on this? I can attached the whole folder if required.
Really appreciate if anyone could help.

[ybapat]

Hello,

I think you have changed fuel at inlet to 0 from fuel inlet boundary. In that case all species have 0 mass fraction on that boundary. This will lead to floating point exception.

-Yogesh

[BakedAlmonds]

If i get what you are trying to say correctly, I have to have an input to satisfy the law of mass conservation. So I should have at least an input which will match the mass flow to the outlet?

I will try changing fuel to 0 and air to 1 then. Will update if it works when i have the chance to try it out.

[BakedAlmonds]

Originally Posted by ybapat

Hello,

I think you have changed fuel at inlet to 0 from fuel inlet boundary. In that case all species have 0 mass fraction on that boundary. This will lead to floating point exception.

-Yogesh

Hi yogesh, my guess of changing the input to air doesnt solve the problem. Apparently my interpretation of what you wanted to say is not true.
Could you further elaborate on the point you are saying from the quote?
I would like to turn off the fuel and input a species of my choice to study the dispersion, thats why i changed the fuel to 0. I used setFields instead of having a constant input.

[ybapat]

Hello,

My guess was that you are getting floating point exceptions as no species is entering from inlet as you have mass fractions for all species =0 at the inlet. But even after setting mass fraction of air at the inlet you are getting the same error then something else might be causing the problem.

Regards,
-Yogesh