[Farshad_Noravesh]

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]

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]

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]

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]

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]

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]

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]

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]

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]

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]

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.