reactingFoam in openfoam1.7.1

Farshad_Noravesh · November 28, 2010, 12:22

Hi,
I see that reactingFoam has always had a janaf temperature problem and it has not been solved even in openfoam1.7. although some attempt has been done for steadyreactingFoam which only and only works in openfoam1.5 , i wonder how people use nonpremixed turbulent combustion(NTC) solver such as reactingFoam. Is there any alternative method for (NTC).
here is my error in openfoam1.7 while simulating a burner inside a 2D furnace:
DILUPBiCG: Solving for hs, Initial residual = 0.998633, Final residual = 7.53168e-14, No Iterations 3
[1]
[1]
[1] --> FOAM FATAL ERROR:
[1] attempt to use janafThermo<equationOfState> out of temperature range 200 -> 5000; T = 5644.91
[1]
[1] From function janafThermo<equationOfState>::checkT(const scalar T) const
[1] in file /home/farshad/OpenFOAM/OpenFOAM-1.7.1/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 63.
[1]
FOAM parallel run aborting
[1]
[1] #0 Foam::error:

rintStack(Foam::Ostream&) in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #1 Foam::error::abort() in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #2 Foam::specieThermo<Foam::janafThermo<Foam:

erfect Gas> >::Hs(double) const in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
[1] #3 Foam::specieThermo<Foam::janafThermo<Foam:

erfect Gas> >::THs(double, double) const in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
[1] #4 Foam::hsPsiMixtureThermo<Foam::reactingMixture<Foa m::sutherlandTransport<Foam::specieThermo<Foam::ja nafThermo<Foam:

erfectGas> > > > >::calculate() in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
[1] #5 Foam::hsPsiMixtureThermo<Foam::reactingMixture<Foa m::sutherlandTransport<Foam::specieThermo<Foam::ja nafThermo<Foam:

erfectGas> > > > >::correct() in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
[1] #6
[1] in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/reactingFoam"
[1] #7 __libc_start_main in "/lib64/libc.so.6"
[1] #8
[1] in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/reactingFoam"
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 1 in communicator MPI_COMM_WORLD
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun has exited due to process rank 1 with PID 9244 on
node farshad exiting without calling "finalize". This may
have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).
--------------------------------------------------------------------------

Farshad_Noravesh · November 28, 2010, 13:26

Hi,
I see that firefoam has the same problem which reactingfoam does.please tell me if you have any solution to debug.
max(fu) = 1
min(fu) = 0
GAMG: Solving for ILambda_0_0, Initial residual = 0.345337, Final residual = 9.43977e-16, No Iterations 1
GAMG: Solving for ILambda_1_0, Initial residual = 0.278267, Final residual = 1.31212e-15, No Iterations 1
GAMG: Solving for ILambda_2_0, Initial residual = 0.242452, Final residual = 1.49772e-15, No Iterations 1
GAMG: Solving for ILambda_3_0, Initial residual = 0.22723, Final residual = 1.88724e-15, No Iterations 1
GAMG: Solving for ILambda_4_0, Initial residual = 0.966266, Final residual = 7.47684e-05, No Iterations 13
GAMG: Solving for ILambda_5_0, Initial residual = 0.985391, Final residual = 4.27461e-05, No Iterations 6
GAMG: Solving for ILambda_6_0, Initial residual = 0.989419, Final residual = 4.79087e-06, No Iterations 3
GAMG: Solving for ILambda_7_0, Initial residual = 0.990994, Final residual = 5.93108e-05, No Iterations 1
GAMG: Solving for ILambda_8_0, Initial residual = 0.991156, Final residual = 6.93516e-15, No Iterations 1
GAMG: Solving for ILambda_9_0, Initial residual = 0.990052, Final residual = 1.14939e-14, No Iterations 1
GAMG: Solving for ILambda_10_0, Initial residual = 0.986653, Final residual = 1.05384e-14, No Iterations 1
GAMG: Solving for ILambda_11_0, Initial residual = 0.934546, Final residual = 1.03974e-14, No Iterations 1
GAMG: Solving for ILambda_12_0, Initial residual = 0.390461, Final residual = 9.31151e-05, No Iterations 7
GAMG: Solving for ILambda_13_0, Initial residual = 0.381106, Final residual = 7.07213e-05, No Iterations 6
GAMG: Solving for ILambda_14_0, Initial residual = 0.380427, Final residual = 8.74956e-05, No Iterations 5
GAMG: Solving for ILambda_15_0, Initial residual = 0.384545, Final residual = 6.21265e-05, No Iterations 4
Radiation solver iter: 1
DILUPBiCG: Solving for hs, Initial residual = 0.75152, Final residual = 8.35865e-08, No Iterations 2

--> FOAM FATAL ERROR:
attempt to use janafThermo<equationOfState> out of temperature range 200 -> 6000; T = 199.989

From function janafThermo<equationOfState>::checkT(const scalar T) const
in file /home/farshad/OpenFOAM/OpenFOAM-1.7.1/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 63.

FOAM aborting

#0 Foam::error:

rintStack(Foam::Ostream&) in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam::specieThermo<Foam::janafThermo<Foam:

erfect Gas> >::Hs(double) const in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#3 Foam::specieThermo<Foam::janafThermo<Foam:

erfect Gas> >::THs(double, double) const in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#4 Foam::hsPsiMixtureThermo<Foam::veryInhomogeneousMi xture<Foam::sutherlandTransport<Foam::specieThermo <Foam::janafThermo<Foam:

erfectGas> > > > >::calculate() in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#5 Foam::hsPsiMixtureThermo<Foam::veryInhomogeneousMi xture<Foam::sutherlandTransport<Foam::specieThermo <Foam::janafThermo<Foam:

erfectGas> > > > >::correct() in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#6
in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/fireFoam"
#7 __libc_start_main in "/lib64/libc.so.6"
#8
in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/fireFoam"
Aborted (core dumped)
[farshad@farshad firefoam]$

nakul · November 29, 2010, 01:15

Hi
I am also uisng reactingFOAM and I also get similar error. I was told in one of the threads that problem is with Boundary conditions.

You may also try changing your BC and if it works for you then please share your insights.

http://www.cfd-online.com/Forums/ope...ctingfoam.html
http://www.cfd-online.com/Forums/ope...ctingfoam.html

Farshad_Noravesh · November 29, 2010, 06:32

Dear Nakul,

Thanks for your help. my error seems to be removed by changing fluxCorrectedVelocity as well as totalpressure and using inletOulet . Now my problem is that the reaction does not start even after 0.2 second.
Is it possible to post a case which has used ignition as well. although i am simulating high temperature air combustion(T_air=1300) and i dont think i need ignition. I am confused about this strange and sensitive reactingFoam.

nakul · November 29, 2010, 06:37

Hi,
It would be better if you provide a little bit more details about your test case. At present it is difficult to say what might be the problem with your case.

Farshad_Noravesh · November 29, 2010, 07:34

Dear Nakul,
Thanks for your help. I really appreciate. The attached folder is my case which is the simulation of a nonpremixed turbulent combustion of a burner inside a 2D HTAC(high temperature air combusiton). I dont know if my k and epsilon boundary conditions are appropriate for reactingFoam.
Kind Regards,
Farshad

nakul · November 29, 2010, 07:56

Hi,
Since you are not getting results with zeroGradient BC. Try "total pressure" BC at outlet for pressure.

k and epsilon can have similar BC types as that of U and T.

You might try fluxCorrected U at exhaust(outlet) with totalpressure BC for P. In that case k, epsilon, T and other scalars can remain as zeroGradient at outlet. Main problem is with outlet BC only. If they are applied correctly I hope that solver should work fine.

Do tell me what results you get. You just need to play around with BC, if you are using the case setup that you uploaded. Just as a precaution run "checkMesh" command to see if your mesh is satisfactory.

Farshad_Noravesh · November 29, 2010, 12:29

Hi Nakul,
I did exactly what you said but the temperature still does not increase more than 1 degree. please give some more help.
Kind Regards,
Farshad

nakul · November 30, 2010, 00:46

Hi

What results did you get for checkMesh?

Actually I am myself trying to sove this problem only. Sorry, I can't help you furhter at this moment. If I would find something I would definitely post it here.

Meanwhile you may also try looking for answers. If you come across something do share it here.
Thanx.

Farshad_Noravesh · November 30, 2010, 02:26

Hi,
My furnace is burning at 2200 now. I only changed mass fraction of oxygen to 1 at inlet and nitrogen to 0. Do you know how can i modify reaction equations to consider nitrogen as well? i see that CH4+02->..
but i dont see nitrogen . How can i include N2 and what are the coefficients? my checkMesh is ok.
Thanks

nakul · November 30, 2010, 03:38

Hi
To include N2 you would have to find a detailed mechanism for this reaction. The elementary reactions of this reaction would also involve N2 and its oxides. You may google for the oxidation mechanism of methane. The coefficents would be specified along with the mechanism.

In addition coefficients for large no. of elementary reactions are given in "$tutorials/combustion/dieselFOAM/chemkin" folder.

Then you may have to specify precise values of mass fractions. As far as my experience goes reactingFOAM becomes unstable as no of reactions increase. But you may give it a try.

November 28, 2010, 12:22	reactingFoam in openfoam1.7.1	#1
Farshad_Noravesh Member Farshad Join Date: Oct 2010 Posts: 76 Rep Power: 16	Hi, I see that reactingFoam has always had a janaf temperature problem and it has not been solved even in openfoam1.7. although some attempt has been done for steadyreactingFoam which only and only works in openfoam1.5 , i wonder how people use nonpremixed turbulent combustion(NTC) solver such as reactingFoam. Is there any alternative method for (NTC). here is my error in openfoam1.7 while simulating a burner inside a 2D furnace: DILUPBiCG: Solving for hs, Initial residual = 0.998633, Final residual = 7.53168e-14, No Iterations 3 [1] [1] [1] --> FOAM FATAL ERROR: [1] attempt to use janafThermo<equationOfState> out of temperature range 200 -> 5000; T = 5644.91 [1] [1] From function janafThermo<equationOfState>::checkT(const scalar T) const [1] in file /home/farshad/OpenFOAM/OpenFOAM-1.7.1/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 63. [1] FOAM parallel run aborting [1] [1] #0 Foam::error:rintStack(Foam::Ostream&) in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so" [1] #1 Foam::error::abort() in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so" [1] #2 Foam::specieThermo<Foam::janafThermo<Foam:erfect Gas> >::Hs(double) const in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so" [1] #3 Foam::specieThermo<Foam::janafThermo<Foam:erfect Gas> >::THs(double, double) const in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so" [1] #4 Foam::hsPsiMixtureThermo<Foam::reactingMixture<Foa m::sutherlandTransport<Foam::specieThermo<Foam::ja nafThermo<Foam:erfectGas> > > > >::calculate() in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so" [1] #5 Foam::hsPsiMixtureThermo<Foam::reactingMixture<Foa m::sutherlandTransport<Foam::specieThermo<Foam::ja nafThermo<Foam:erfectGas> > > > >::correct() in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so" [1] #6 [1] in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/reactingFoam" [1] #7 __libc_start_main in "/lib64/libc.so.6" [1] #8 [1] in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/reactingFoam" -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 1 in communicator MPI_COMM_WORLD with errorcode 1. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them. -------------------------------------------------------------------------- -------------------------------------------------------------------------- mpirun has exited due to process rank 1 with PID 9244 on node farshad exiting without calling "finalize". This may have caused other processes in the application to be terminated by signals sent by mpirun (as reported here). --------------------------------------------------------------------------

November 28, 2010, 13:26	firefoam has the same problem as reactingFoam does	#2
Farshad_Noravesh Member Farshad Join Date: Oct 2010 Posts: 76 Rep Power: 16	Hi, I see that firefoam has the same problem which reactingfoam does.please tell me if you have any solution to debug. max(fu) = 1 min(fu) = 0 GAMG: Solving for ILambda_0_0, Initial residual = 0.345337, Final residual = 9.43977e-16, No Iterations 1 GAMG: Solving for ILambda_1_0, Initial residual = 0.278267, Final residual = 1.31212e-15, No Iterations 1 GAMG: Solving for ILambda_2_0, Initial residual = 0.242452, Final residual = 1.49772e-15, No Iterations 1 GAMG: Solving for ILambda_3_0, Initial residual = 0.22723, Final residual = 1.88724e-15, No Iterations 1 GAMG: Solving for ILambda_4_0, Initial residual = 0.966266, Final residual = 7.47684e-05, No Iterations 13 GAMG: Solving for ILambda_5_0, Initial residual = 0.985391, Final residual = 4.27461e-05, No Iterations 6 GAMG: Solving for ILambda_6_0, Initial residual = 0.989419, Final residual = 4.79087e-06, No Iterations 3 GAMG: Solving for ILambda_7_0, Initial residual = 0.990994, Final residual = 5.93108e-05, No Iterations 1 GAMG: Solving for ILambda_8_0, Initial residual = 0.991156, Final residual = 6.93516e-15, No Iterations 1 GAMG: Solving for ILambda_9_0, Initial residual = 0.990052, Final residual = 1.14939e-14, No Iterations 1 GAMG: Solving for ILambda_10_0, Initial residual = 0.986653, Final residual = 1.05384e-14, No Iterations 1 GAMG: Solving for ILambda_11_0, Initial residual = 0.934546, Final residual = 1.03974e-14, No Iterations 1 GAMG: Solving for ILambda_12_0, Initial residual = 0.390461, Final residual = 9.31151e-05, No Iterations 7 GAMG: Solving for ILambda_13_0, Initial residual = 0.381106, Final residual = 7.07213e-05, No Iterations 6 GAMG: Solving for ILambda_14_0, Initial residual = 0.380427, Final residual = 8.74956e-05, No Iterations 5 GAMG: Solving for ILambda_15_0, Initial residual = 0.384545, Final residual = 6.21265e-05, No Iterations 4 Radiation solver iter: 1 DILUPBiCG: Solving for hs, Initial residual = 0.75152, Final residual = 8.35865e-08, No Iterations 2 --> FOAM FATAL ERROR: attempt to use janafThermo<equationOfState> out of temperature range 200 -> 6000; T = 199.989 From function janafThermo<equationOfState>::checkT(const scalar T) const in file /home/farshad/OpenFOAM/OpenFOAM-1.7.1/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 63. FOAM aborting #0 Foam::error:rintStack(Foam::Ostream&) in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::error::abort() in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so" #2 Foam::specieThermo<Foam::janafThermo<Foam:erfect Gas> >::Hs(double) const in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so" #3 Foam::specieThermo<Foam::janafThermo<Foam:erfect Gas> >::THs(double, double) const in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so" #4 Foam::hsPsiMixtureThermo<Foam::veryInhomogeneousMi xture<Foam::sutherlandTransport<Foam::specieThermo <Foam::janafThermo<Foam:erfectGas> > > > >::calculate() in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so" #5 Foam::hsPsiMixtureThermo<Foam::veryInhomogeneousMi xture<Foam::sutherlandTransport<Foam::specieThermo <Foam::janafThermo<Foam:erfectGas> > > > >::correct() in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so" #6 in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/fireFoam" #7 __libc_start_main in "/lib64/libc.so.6" #8 in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/fireFoam" Aborted (core dumped) [farshad@farshad firefoam]$

November 30, 2010, 02:26	reactions for N2	#10
Farshad_Noravesh Member Farshad Join Date: Oct 2010 Posts: 76 Rep Power: 16	Hi, My furnace is burning at 2200 now. I only changed mass fraction of oxygen to 1 at inlet and nitrogen to 0. Do you know how can i modify reaction equations to consider nitrogen as well? i see that CH4+02->.. but i dont see nitrogen . How can i include N2 and what are the coefficients? my checkMesh is ok. Thanks

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
help with reactingFoam	piccinini	OpenFOAM	8	October 4, 2013 06:08
reactingFoam wedge handling wrong U	dhondupant	OpenFOAM Bugs	1	December 9, 2010 08:34
reactingFoam - turbulent reacting flow	hamburgFoam	OpenFOAM	0	December 7, 2009 13:57
information about reactingFOAM	arezo.namazi	OpenFOAM	0	October 26, 2009 05:11
ReactingFoam solver	muthukaalai	OpenFOAM Running, Solving & CFD	1	June 16, 2008 14:36

November 29, 2010, 01:15		#3
nakul Senior Member Nakul Join Date: Apr 2010 Location: India Posts: 147 Rep Power: 16	Hi I am also uisng reactingFOAM and I also get similar error. I was told in one of the threads that problem is with Boundary conditions. You may also try changing your BC and if it works for you then please share your insights. http://www.cfd-online.com/Forums/ope...ctingfoam.html http://www.cfd-online.com/Forums/ope...ctingfoam.html

November 29, 2010, 06:32		#4
Farshad_Noravesh Member Farshad Join Date: Oct 2010 Posts: 76 Rep Power: 16	Dear Nakul, Thanks for your help. my error seems to be removed by changing fluxCorrectedVelocity as well as totalpressure and using inletOulet . Now my problem is that the reaction does not start even after 0.2 second. Is it possible to post a case which has used ignition as well. although i am simulating high temperature air combustion(T_air=1300) and i dont think i need ignition. I am confused about this strange and sensitive reactingFoam.

November 29, 2010, 06:37		#5
nakul Senior Member Nakul Join Date: Apr 2010 Location: India Posts: 147 Rep Power: 16	Hi, It would be better if you provide a little bit more details about your test case. At present it is difficult to say what might be the problem with your case.

November 29, 2010, 07:56		#7
nakul Senior Member Nakul Join Date: Apr 2010 Location: India Posts: 147 Rep Power: 16	Hi, Since you are not getting results with zeroGradient BC. Try "total pressure" BC at outlet for pressure. k and epsilon can have similar BC types as that of U and T. You might try fluxCorrected U at exhaust(outlet) with totalpressure BC for P. In that case k, epsilon, T and other scalars can remain as zeroGradient at outlet. Main problem is with outlet BC only. If they are applied correctly I hope that solver should work fine. Do tell me what results you get. You just need to play around with BC, if you are using the case setup that you uploaded. Just as a precaution run "checkMesh" command to see if your mesh is satisfactory.

November 29, 2010, 12:29	still the same problem	#8
Farshad_Noravesh Member Farshad Join Date: Oct 2010 Posts: 76 Rep Power: 16	Hi Nakul, I did exactly what you said but the temperature still does not increase more than 1 degree. please give some more help. Kind Regards, Farshad

November 30, 2010, 00:46		#9
nakul Senior Member Nakul Join Date: Apr 2010 Location: India Posts: 147 Rep Power: 16	Hi What results did you get for checkMesh? Actually I am myself trying to sove this problem only. Sorry, I can't help you furhter at this moment. If I would find something I would definitely post it here. Meanwhile you may also try looking for answers. If you come across something do share it here. Thanx.

November 30, 2010, 03:38		#11
nakul Senior Member Nakul Join Date: Apr 2010 Location: India Posts: 147 Rep Power: 16	Hi To include N2 you would have to find a detailed mechanism for this reaction. The elementary reactions of this reaction would also involve N2 and its oxides. You may google for the oxidation mechanism of methane. The coefficents would be specified along with the mechanism. In addition coefficients for large no. of elementary reactions are given in "$tutorials/combustion/dieselFOAM/chemkin" folder. Then you may have to specify precise values of mass fractions. As far as my experience goes reactingFOAM becomes unstable as no of reactions increase. But you may give it a try.