CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

error with reactingFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 18, 2016, 07:15
Default error with reactingFoam
  #1
Member
 
Almond Wong
Join Date: May 2016
Posts: 68
Rep Power: 10
BakedAlmonds is on a distinguished road
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.

Code:
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.

Last edited by BakedAlmonds; June 21, 2016 at 07:37.
BakedAlmonds is offline   Reply With Quote

Old   June 20, 2016, 08:04
Default
  #2
Senior Member
 
Yogesh Bapat
Join Date: Oct 2010
Posts: 102
Rep Power: 15
ybapat is on a distinguished road
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
ybapat is offline   Reply With Quote

Old   June 21, 2016, 07:36
Default
  #3
Member
 
Almond Wong
Join Date: May 2016
Posts: 68
Rep Power: 10
BakedAlmonds is on a distinguished road
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 is offline   Reply With Quote

Old   June 21, 2016, 09:29
Default
  #4
Member
 
Almond Wong
Join Date: May 2016
Posts: 68
Rep Power: 10
BakedAlmonds is on a distinguished road
Quote:
Originally Posted by ybapat View Post
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.
BakedAlmonds is offline   Reply With Quote

Old   June 22, 2016, 02:21
Default
  #5
Senior Member
 
Yogesh Bapat
Join Date: Oct 2010
Posts: 102
Rep Power: 15
ybapat is on a distinguished road
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
ybapat is offline   Reply With Quote

Reply

Tags
reactingfoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
How to switch off combustion and reaction in reactingFoam shenzhou1987 OpenFOAM Running, Solving & CFD 16 October 26, 2017 15:31
calculation of thermal conductivity in reactingFoam sto16 OpenFOAM Running, Solving & CFD 11 April 24, 2017 23:46
calculate flame speed using reactingFoam IColin OpenFOAM Running, Solving & CFD 0 February 4, 2014 15:14
reactingFoam wedge handling wrong U dhondupant OpenFOAM Bugs 1 December 9, 2010 07:34
reactingFoam - turbulent reacting flow hamburgFoam OpenFOAM 0 December 7, 2009 12:57


All times are GMT -4. The time now is 15:19.