CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

reactingFoam in openfoam1.7.1

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 28, 2010, 12:22
Unhappy reactingFoam in openfoam1.7.1
  #1
Member
 
Farshad
Join Date: Oct 2010
Posts: 76
Rep Power: 6
Farshad_Noravesh is on a distinguished road
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 is offline   Reply With Quote

Old   November 28, 2010, 13:26
Unhappy firefoam has the same problem as reactingFoam does
  #2
Member
 
Farshad
Join Date: Oct 2010
Posts: 76
Rep Power: 6
Farshad_Noravesh is on a distinguished road
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]$
Farshad_Noravesh is offline   Reply With Quote

Old   November 29, 2010, 01:15
Default
  #3
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
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
Problem With reactingFOAM
nakul is offline   Reply With Quote

Old   November 29, 2010, 06:32
Default
  #4
Member
 
Farshad
Join Date: Oct 2010
Posts: 76
Rep Power: 6
Farshad_Noravesh is on a distinguished road
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.
Farshad_Noravesh is offline   Reply With Quote

Old   November 29, 2010, 06:37
Default
  #5
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
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.
nakul is offline   Reply With Quote

Old   November 29, 2010, 07:34
Default
  #6
Member
 
Farshad
Join Date: Oct 2010
Posts: 76
Rep Power: 6
Farshad_Noravesh is on a distinguished road
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
Attached Files
File Type: zip 2dfurnace1.zip (20.0 KB, 37 views)
Farshad_Noravesh is offline   Reply With Quote

Old   November 29, 2010, 07:56
Default
  #7
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
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.
nakul is offline   Reply With Quote

Old   November 29, 2010, 12:29
Unhappy still the same problem
  #8
Member
 
Farshad
Join Date: Oct 2010
Posts: 76
Rep Power: 6
Farshad_Noravesh is on a distinguished road
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
Farshad_Noravesh is offline   Reply With Quote

Old   November 30, 2010, 00:46
Default
  #9
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
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.
nakul is offline   Reply With Quote

Old   November 30, 2010, 02:26
Question reactions for N2
  #10
Member
 
Farshad
Join Date: Oct 2010
Posts: 76
Rep Power: 6
Farshad_Noravesh is on a distinguished road
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
Farshad_Noravesh is offline   Reply With Quote

Old   November 30, 2010, 03:38
Default
  #11
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
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.
nakul is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
help with reactingFoam piccinini OpenFOAM 8 October 4, 2013 05: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 13:36


All times are GMT -4. The time now is 21:39.