CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   reactingFoam in openfoam1.7.1 (http://www.cfd-online.com/Forums/openfoam/82500-reactingfoam-openfoam1-7-1-a.html)

Farshad_Noravesh November 28, 2010 12:22

reactingFoam in openfoam1.7.1
 
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::printStack(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::perfect Gas> >::Hs(double) const in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
[1] #3 Foam::specieThermo<Foam::janafThermo<Foam::perfect 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::perfectGas> > > > >::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::perfectGas> > > > >::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

firefoam has the same problem as reactingFoam does
 
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::printStack(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::perfect Gas> >::Hs(double) const in "/home/farshad/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#3 Foam::specieThermo<Foam::janafThermo<Foam::perfect 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::perfectGas> > > > >::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::perfectGas> > > > >::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

1 Attachment(s)
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

still the same problem
 
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

reactions for N2
 
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.


All times are GMT -4. The time now is 03:07.