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

Flame gets flushed away in reactingFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 25, 2011, 05:00
Default Flame gets flushed away in reactingFoam
  #1
Member
 
Tibo
Join Date: Jun 2011
Posts: 68
Rep Power: 5
megacrout is on a distinguished road
Hi,

I am simulating a premixed flame with reactingFoam. I use methane with oxygen-rich air and a semi-developped chemistry (4 reactions, one reversible) with pretty unrealistic but still plausible kinetic data (however the problem is unlikely to be found here).

I set the inlet temperature to 800K but the mixture does not ignite. So I set the initial internal temperature to 2000K. The flame does ignite but the heat does not seem to diffuse fast enough so that the 800K mixture slowly replaces the hot gas and the flame disappears. Using a smaller inlet velocity does not solve the problem. I actually first tried with 1m/s and this should work.

I reduced the mesh down to 5x5mm cells. Courant numbers are perfectly fine.

Does it have sth to do with the solver? I looked pretty intensively at the equations and donīt see where the problem could lie (premixed gases should be fine).

Thanks for your help.

Tibo
megacrout is offline   Reply With Quote

Old   October 27, 2011, 05:11
Default
  #2
New Member
 
Dhondu Pant
Join Date: Nov 2009
Posts: 6
Rep Power: 6
dhondupant is on a distinguished road
Quote:
Originally Posted by megacrout View Post
Hi,

I am simulating a premixed flame with reactingFoam. I use methane with oxygen-rich air and a semi-developped chemistry (4 reactions, one reversible) with pretty unrealistic but still plausible kinetic data (however the problem is unlikely to be found here).

I set the inlet temperature to 800K but the mixture does not ignite. So I set the initial internal temperature to 2000K. The flame does ignite but the heat does not seem to diffuse fast enough so that the 800K mixture slowly replaces the hot gas and the flame disappears. Using a smaller inlet velocity does not solve the problem. I actually first tried with 1m/s and this should work.

I reduced the mesh down to 5x5mm cells. Courant numbers are perfectly fine.

Does it have sth to do with the solver? I looked pretty intensively at the equations and donīt see where the problem could lie (premixed gases should be fine).

Thanks for your help.

Tibo
Have you tried 1D solvers with your 4 reactions chemical kinetic so far (like Cantera) or is your flame no laminar flame at all?
dhondupant is offline   Reply With Quote

Old   November 2, 2011, 04:51
Default
  #3
Member
 
Tibo
Join Date: Jun 2011
Posts: 68
Rep Power: 5
megacrout is on a distinguished road
My flame is laminar but Iīm going to study turbulent flames as soon as it works with laminar ones. Plus I need to solve it 2D (I need both axial and radial profiles). So Iīd prefer sticking to reactingFoam.

Do you know how to solve this problem? or at least where this problem could come from?

Btw, I actually got good results with other kinetic data, so it seems that I was wrong assuming the problem could not be here. however, I do not know how to modify my data so that I get results out of it...
megacrout is offline   Reply With Quote

Old   November 2, 2011, 11:27
Default
  #4
New Member
 
Jens Keller
Join Date: Nov 2011
Location: Karlsruhe
Posts: 5
Rep Power: 4
jenzkeller is on a distinguished road
Quote:
Originally Posted by megacrout View Post
My flame is laminar but Iīm going to study turbulent flames as soon as it works with laminar ones. Plus I need to solve it 2D (I need both axial and radial profiles). So Iīd prefer sticking to reactingFoam.

Do you know how to solve this problem? or at least where this problem could come from?

Btw, I actually got good results with other kinetic data, so it seems that I was wrong assuming the problem could not be here. however, I do not know how to modify my data so that I get results out of it...
Well I had to model the laminar "Bennett"-flame using a self modified OpenFOAM solver.

In my case using any reacting mechanism with less number of reactions as in the GRI 1.2 was giving me bad results compared to experimental data.

So consider using much more chemical reactions for a laminar flame as for a turbulent flame.

You should also know that reactinFoam means physically Sc=1 using a unity Lewis number assumption which is absolutely inappropriate for laminar flames.

You might want to modify reactingFOAM too, by giving each species its own Lewis number.

cheers
jenzkeller is offline   Reply With Quote

Old   November 2, 2011, 12:33
Default
  #5
Member
 
Tibo
Join Date: Jun 2011
Posts: 68
Rep Power: 5
megacrout is on a distinguished road
Hi Jens,

Thanks for your answer too.

First of all, I am comparing OpenFOAM results with data from simulations obtained with other programms, not experimental data. So getting "real" data is not quite the problem even though non-conservation of inert species or non-sustainability of the flame obviously arenīt ok.
I managed simulating one-step reactions but it just does not seem to work with an intermediate number of reactions (for which I compute the kinetic data myself).

I read that detail about Sc=1 but I donīt understand where it comes from and - more important - what it implies. As for Le=1, I am not quite sure it still is the fact (from OF1.7.1 on) because of the commented line in hsEqn.H.

So, apart from the number of reactions and Sc numbers, do you know what I should have a look at to improve my results (i.e. get my flame to keep burning)?

Thank you very much.
megacrout is offline   Reply With Quote

Old   November 3, 2011, 04:24
Default
  #6
New Member
 
Jens Keller
Join Date: Nov 2011
Location: Karlsruhe
Posts: 5
Rep Power: 4
jenzkeller is on a distinguished road
Quote:
Originally Posted by megacrout View Post
So, apart from the number of reactions and Sc numbers, do you know what I should have a look at to improve my results (i.e. get my flame to keep burning)?
Well you can decrease the Courant-number which should not be above Co=0.2 or you can try to use a coarser mesh for an initial flame simulation.

cheers
jenzkeller is offline   Reply With Quote

Old   November 7, 2011, 05:03
Default
  #7
Member
 
Tibo
Join Date: Jun 2011
Posts: 68
Rep Power: 5
megacrout is on a distinguished road
Quote:
Originally Posted by megacrout View Post
Hi Jens,
I read that detail about Sc=1 but I donīt understand where it comes from and - more important - what it implies. As for Le=1, I am not quite sure it still is the fact (from OF1.7.1 on) because of the commented line in hsEqn.H.
I get it now: the Schmidt number is Sc = mu / rho / D. Setting it to unity means D = mu/rho which explains why the specieīs mass equation in reactingFoam (YEqn.H) contains mu instead of D to decribe the mass diffusion of the species!

However, I keep thinking a unity Lewis number does not apply to reactingFoam anymore. Can anyone help?
megacrout is offline   Reply With Quote

Old   November 7, 2011, 05:33
Default
  #8
New Member
 
Jens Keller
Join Date: Nov 2011
Location: Karlsruhe
Posts: 5
Rep Power: 4
jenzkeller is on a distinguished road
Quote:
Originally Posted by megacrout View Post
I get it now: the Schmidt number is Sc = mu / rho / D. Setting it to unity means D = mu/rho which explains why the specieīs mass equation in reactingFoam (YEqn.H) contains mu instead of D to decribe the mass diffusion of the species!

However, I keep thinking a unity Lewis number does not apply to reactingFoam anymore. Can anyone help?
Exactly.

I can advise you to read the first pages of Thierry Poinsot's "Theoretical and Numerical Combustion"

Than you will understand and see the simplifications done in the energy and species equation of OpenFOAM.

If you explicitly want to calculate Le=1 within your species equation just simply replace thermo->mu() with thermo->alpha().
jenzkeller is offline   Reply With Quote

Reply

Tags
ignition, reactingfoam

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
ReactingFoam Laminar flame simulation transport properties gbansal OpenFOAM Running, Solving & CFD 3 November 23, 2012 16:03
reactingFoam wedge handling wrong U dhondupant OpenFOAM Bugs 1 December 9, 2010 07:34
reactingFoam: adiabatic flame temperature rFkhemek OpenFOAM 2 July 31, 2009 06:42
reactingFoam: New kinetic model errors rFkhemek OpenFOAM 0 July 11, 2009 15:13
Crosswind flame with reactingFoam torvic OpenFOAM Running, Solving & CFD 1 September 10, 2007 17:48


All times are GMT -4. The time now is 13:47.