CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

XiFoam solver: changes needed to switch from air-combustion to oxy-combustion?

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

LinkBack Thread Tools Display Modes
Old   May 15, 2011, 03:06
Default reactingFoam with variable diffusivity
Yashar Afarin
Join Date: May 2010
Location: Toronto- Canada
Posts: 40
Rep Power: 9
yashar.afarin is on a distinguished road
Send a message via Skype™ to yashar.afarin
Hi nakul,

thanks for your posts.
could you please send me some information about implementing reactingFoam with variable diffusivity. I am Modeling methane turbulent flame and in part of my research I want to assess Variable diffusivity effect on my results.

sincerely yours
yashar.afarin is offline   Reply With Quote

Old   May 16, 2011, 04:40
Senior Member
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 9
nakul is on a distinguished road

Implementing the feature of variable diffusivity right from the scratch may be quite long and difficult. But I think that you may find some thing on the net in this regard. I haven't checked but if you come across something do post it here. Stuff you may find on net may not be even for reactingFoam but you would get an idea of editing existing thermodynamic models which is the fundamental way of implementing variable diffusivity.

Alternatively you may assume a constant value for Lewis no. and multiply it with alpha which can be "variable" in OF. This way you would be able to have variable diffusivity in a much simpler way as far as reactingFoam is concerned.

nakul is offline   Reply With Quote

Old   August 27, 2013, 08:34
Default XiFoam output files
Join Date: Jul 2013
Posts: 38
Rep Power: 6
ni-openfoam-user is on a distinguished road
As a continuation of the previous discussion regarding calculation of temperature when using XiFoam. I am quite new to OpenFoam and as per my understanding, please correct if I am mistaken, for each file in the 0 directory there is an associated IOobject contained within the code.

For example within the /0 folder in the XiFoam tutorial there are the files: alphat, b, epsilon, k, mut, p, Su, T, Tu, U, Xi

Therefore it is my understanding that within the code there is an IOoject for each of these parameters.

For example within createFields.H of the XiFoam solver (contained in applications/solvers/combustion/XiFoam) we have:

volvectorField U

I have been able to track down the associated IOobject for each of the parameters which are saved to file apart from b, p and T.

I notice in createFields.H, there is:

volScalarField& p = thermo.p();


volScalarField& b = composition.Y(b);

and in previous releases of OpenFoam, createFields.H also contained: const volScalarField& T = thermo.T();

Can you please comment on how/where b, p and T are calculated and saved to file and where their associated IOobject, if it exists, is located.

Many thanks,

ni-openfoam-user is offline   Reply With Quote

Old   August 28, 2013, 08:21
Default Calculate Temperature in hhuMixtureThermo
New Member
Vito Raso
Join Date: Apr 2013
Location: Bari(Italy)
Posts: 7
Rep Power: 6
Vito31388 is on a distinguished road
Originally Posted by ehsanyasari View Post
I checked your solver. I presume that "T" is calcualted according to the "hs" and "Cp", but in a little bit more complex way. Because "Cp" is also dependent to temperature.
What I found is that:
1-in "hsEqn.H": thermo.correct(); is implemented in :
in online version in line 00136.

2- in "correct()" function you can find another function called "calcualte()" which also implemented in the same file as previous one and in line 00068.

3-in this function you also can find function "THs".

The function "THs" which calculate the temperature based on "Hs" and "Cp" ,in a iterative method, is implemented in :


Please let me know if you think it is correct ?
Hi Ehsan your post was very helpful but I've a doubt about temperature calculation in library hhuMixtureThermo.C; I'd like to ask you if temperature is calculated with TH() function or based on Tb and Tu using products and reactans JANAF coefficients i.e.
T= Tu(b) + Tb(1-b).
Thanks in advance, Vito
Vito31388 is offline   Reply With Quote

Old   October 17, 2014, 12:21
New Member
Join Date: Aug 2012
Location: Germany
Posts: 25
Rep Power: 7
sharifi is on a distinguished road

I'm working with XiFoam and I have the same understanding problem about calculating T.
In XiFoam there is an equation for "h" which I'm pretty sure Its the total enthalpy, so that means for an adiabatic case it would not be used.

In my understanding, theoretically If we have the adiabatic temperature of the flame we can calculate the tempreture as function of "b" every where.
I can imagine that we can find out the adiabatic temperature using formation enthalpy of reactant and product and the "cp" and "mass fraction" of fuel.

I could not see the links that nakul sent, and its still an open question for me.
It would be very nice if you guys could help me to understand that.

sharifi is offline   Reply With Quote

Old   April 13, 2016, 10:05
Default XiFoam-Hydrogen
New Member
Reza khodadadi
Join Date: Apr 2011
Location: London, UK
Posts: 24
Rep Power: 8
reza_65 is on a distinguished road
Hi guys,

I was using rhoreactingFoam for ddt of Hydrogen, and now I want to switch to the XiFoam. Can you please tell me how can I have the Hydrogen properties and coefficients for simulating Hydrogen DDT?
Moreover, I used to set high temperature circle for ignition in setFields, can you please tell me whether I can use the same way or no I should use the ignition in combustion properties?

reza_65 is offline   Reply With Quote


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
Working directory via command line Luiz CFX 4 March 6, 2011 21:02
air bubble is disappear increasing time using vof xujjun CFX 9 June 9, 2009 07:59
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08
A role of combustion modeling in flow solver? Yoon Main CFD Forum 0 November 26, 2006 13:00
compressible two phase flow in CFX4.4 youngan CFX 0 July 1, 2003 23:32

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