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

Combustion Modeling-Using chemFoam in reactive flow

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

Like Tree4Likes
  • 1 Post By Hanzo
  • 1 Post By Hanzo
  • 1 Post By Hanzo
  • 1 Post By mhsn

Reply
 
LinkBack Thread Tools Display Modes
Old   October 24, 2012, 10:41
Default Combustion Modeling-Using chemFoam in reactive flow
  #1
New Member
 
Mhsn
Join Date: Oct 2012
Posts: 24
Rep Power: 4
mhsn is on a distinguished road
Dear Foamers,

I like to model premixed combustion in a narrow channel using detailed chemistry. I know that chemFoam is a solver that can implement detailed reaction mechanisms; however, it is only for one cell Does anyone know how I can use the capability of this solver in other combustion solvers, for instance in XiFoam?

Any help would be appreciated!
Thanks
mhsn is offline   Reply With Quote

Old   October 31, 2012, 04:18
Default
  #2
Member
 
Join Date: Nov 2010
Location: Tokyo / Japan
Posts: 40
Rep Power: 5
Hanzo is on a distinguished road
Hi mhsn,

I think chemFoam is just a validation solver for the chemistry provided in OpenFoam. One can use it to analyse reactions and compare with other solvers like chemkin etc:

http://www.openfoam.org/version2.0.0/chemistry.php

I think, the detailed chemistry can be used by the other solvers as well (e.g. you can set up a reactingFoam case using the same chemistry functionality like in the chemFoam case).

Or is there some functionality I am not considering?
Hanzo is offline   Reply With Quote

Old   October 31, 2012, 13:11
Default
  #3
New Member
 
Mhsn
Join Date: Oct 2012
Posts: 24
Rep Power: 4
mhsn is on a distinguished road
Quote:
Originally Posted by Hanzo View Post
Hi mhsn,

I think, the detailed chemistry can be used by the other solvers as well (e.g. you can set up a reactingFoam case using the same chemistry functionality like in the chemFoam case).
Hanzo, Thanks for your reply. Actually, what you're saying about using the same chemistry functionality like in the chemFoam case is my question. I'm very new to OF and not sure how I can do it! Can you be more specific or link me to something that has been done like that so that I can get some clues how to start?

Thanks
mhsn is offline   Reply With Quote

Old   October 31, 2012, 23:30
Default
  #4
Member
 
Join Date: Nov 2010
Location: Tokyo / Japan
Posts: 40
Rep Power: 5
Hanzo is on a distinguished road
Quote:
Originally Posted by mhsn View Post
Can you be more specific or link me to something that has been done like that so that I can get some clues how to start?

Thanks
I think, using the following is a nice start:

http://www.tfd.chalmers.se/~hani/kur...actingFoam.pdf

It's a tutorial written by Andreas Lundstrom and he introduces two own reactions, shows how to write that down in chemkin files and also explains how to set up the 7 high-low coeffs for this simple case.

Later, you could use the chemkinToFoam tool to convert the input into the nativeOpenFoam chemistry format, because, I think, this format is much more flexible and easier to read.

Last but not least: there are some minor typos in the tutorial when it comes to compute some values, so don't get confused by this

Also, the Tutorial in the OpenFoam Wiki:
http://openfoamwiki.net/index.php/Tu..._firstTutorial

however, here the reactions are already quite complex
mhsn likes this.
Hanzo is offline   Reply With Quote

Old   November 1, 2012, 16:21
Default
  #5
New Member
 
Mhsn
Join Date: Oct 2012
Posts: 24
Rep Power: 4
mhsn is on a distinguished road
Thanks for the information
I will take a look at these tutorials to see how it works.
mhsn is offline   Reply With Quote

Old   November 2, 2012, 02:00
Default
  #6
Member
 
Join Date: Nov 2010
Location: Tokyo / Japan
Posts: 40
Rep Power: 5
Hanzo is on a distinguished road
You are welcome!

I forgot, here is the link to the case files for the Lundstrom tutorial.

http://www.tfd.chalmers.se/~hani/kur...utorial.tar.gz
mhsn likes this.
Hanzo is offline   Reply With Quote

Old   November 9, 2012, 16:47
Default
  #7
New Member
 
Mhsn
Join Date: Oct 2012
Posts: 24
Rep Power: 4
mhsn is on a distinguished road
Hanzo,
the tutorial you sent to me helped a lot. Now I can run my test case using chemkin mechanisms. But, I have another issue now. What if I want to use transport properties from chemkin files? Do you know how I can read the transport properties like I do for them file in chemkin format?
Thanks
mhsn is offline   Reply With Quote

Old   November 13, 2012, 04:40
Default
  #8
Member
 
Join Date: Nov 2010
Location: Tokyo / Japan
Posts: 40
Rep Power: 5
Hanzo is on a distinguished road
Quote:
Originally Posted by mhsn View Post
Hanzo,
the tutorial you sent to me helped a lot. Now I can run my test case using chemkin mechanisms. But, I have another issue now. What if I want to use transport properties from chemkin files? Do you know how I can read the transport properties like I do for them file in chemkin format?
Thanks
sorry, what do you mean with transport properties from chemkin files?

you only have the reaction parameters and the thermal properties stored in the chemkin files.
Hanzo is offline   Reply With Quote

Old   November 13, 2012, 10:20
Default
  #9
New Member
 
Mhsn
Join Date: Oct 2012
Posts: 24
Rep Power: 4
mhsn is on a distinguished road
Basically, there are three different chemkin files for each mechanism: 1) for reactions 2) for thermal properties and 3) for transport properties such as diffusion coefficients, thermal conductivity, viscosity and such other properties.
OF has readers for reactions and thermal properties, and I'm not sure if there is any way that I can use chemkin transport file! Do you have any idea how these properties can be implemented for reactingFoam?
mhsn is offline   Reply With Quote

Old   November 13, 2012, 21:27
Default
  #10
Member
 
Join Date: Nov 2010
Location: Tokyo / Japan
Posts: 40
Rep Power: 5
Hanzo is on a distinguished road
Quote:
Originally Posted by mhsn View Post
Basically, there are three different chemkin files for each mechanism: 1) for reactions 2) for thermal properties and 3) for transport properties such as diffusion coefficients, thermal conductivity, viscosity and such other properties.
I never used Chemkin but now I see your point.

Quote:
Originally Posted by mhsn View Post
OF has readers for reactions and thermal properties, and I'm not sure if there is any way that I can use chemkin transport file! Do you have any idea how these properties can be implemented for reactingFoam?
When it comes to the transport properties I usually have a look on the equation the solver tries to solve. In reactingFoam if you have a look in the Yeqn.H

http://foam.sourceforge.net/docs/cpp/a02406_source.html
Code:
00024             fvScalarMatrix YiEqn 
00025             ( 
00026                 fvm::ddt(rho, Yi) 
00027               + mvConvection->fvmDiv(phi, Yi) 
00028               - fvm::laplacian(turbulence->muEff(), Yi) 
00029              == 
00030                 combustion->R(Yi) 
00031             );
The Diffusion coefficient is set to muEff. This comes from the underlying assumption of unity turbulent Lewis number which is commonly used in turbulent combustion (http://www.openfoam.org/mantisbt/view.php?id=277) .

To input your own diffusion coefficients you should do something similar like here
Non-constant Diffusion Coefficient

By the way, you are working a lot with the chemistry in openFoam. Maybe you can have a small look on my problem -> Solver Validation chemFoam using analytical Solution
I still have some conversion issues there
mhsn likes this.
Hanzo is offline   Reply With Quote

Old   November 15, 2012, 16:38
Default
  #11
New Member
 
Mhsn
Join Date: Oct 2012
Posts: 24
Rep Power: 4
mhsn is on a distinguished road
Hi Hanzo, thanks for your reply. Actually, because for the work I'm going to do there are plenty of species (around 50) with different mechanisms, that should not be the very right way to follow that. I really like a way like how thermal properties files are read in chemkin format.
About your question, I checked that post. It is really weird that you have discrepancy by factor 7! I think you may have to double check all units that you've used to see if they are consistent with other units used in the solver! That may be the reason!

Thanks
Hanzo likes this.
mhsn 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
Modeling the mixing of air and kerosene in a flow channel StefanG CFX 3 June 11, 2012 20:21
Modeling Turbulent Reactive Flow sanjibdsharma OpenFOAM 43 August 26, 2010 11:23
Modeling liquid-phase reactive flow sanjibdsharma OpenFOAM 0 October 22, 2009 05:42
mass flow in is not equal to mass flow out saii CFX 2 September 18, 2009 08:07
What is the difference between liquid reactive flow and gas reactive flow? James Main CFD Forum 6 May 15, 2009 12:14


All times are GMT -4. The time now is 18:57.