CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

XiFoam solver

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: 7
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: 7
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: 37
Rep Power: 4
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: 4
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: 23
Rep Power: 5
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


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
Quarter Burner mesh with periosic condition SamCanuck FLUENT 2 August 31, 2011 11:34
Working directory via command line Luiz CFX 4 March 6, 2011 21:02
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08
compressible two phase flow in CFX4.4 youngan CFX 0 July 1, 2003 23:32
Setting a B.C using UserFortran in 4.3 tokai CFX 10 July 17, 2001 16:25

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