CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM

Enthalpy equation for multicomponent mixtures

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

Like Tree1Likes
  • 1 Post By alberto_cuoci

LinkBack Thread Tools Search this Thread Display Modes
Old   September 16, 2009, 05:11
Default Enthalpy equation for multicomponent mixtures
New Member
Alberto Cuoci
Join Date: Aug 2009
Posts: 3
Rep Power: 16
alberto_cuoci is on a distinguished road
Hello everybody,

I would like to write a laminar, steady-state solver for multicomponent, reacting mixtures. I have some questions about the enthalpy equation which is used in several solvers for multicomponent mixtures (in steady state conditions, SIMPLE algorithm):

fvm::div(phi, h)
- fvm::Sp(fvc::div(phi), h)
- fvm::laplacian(turbulence->alphaEff(), h)
fvc::div(phi/fvc::interpolate(rho)*fvc::interpolate(p, "div(U,p)"))
- p*fvc::div(phi/fvc::interpolate(rho))

1. Thermal diffusivity

In laminar conditions the turbulence->alphaEff() term should become k/Cp (where k is the thermal conductivity). Is it correct?

2. Pressure term
I cannot understand why the term associated to the pressure is written as:


I expected just the first term:


3. Temperature gradient

In the enthalpy equation the div(k grad(T)) term is written as laplacian(k/Cp grad(h)) since grad(h)=Cp grad(T) .
However for multicomponent mixtures we have:

grad(h)=Cp grad(T) + sum h_i grad(Y_i)

where i is the species index and the sum is over all the species.Therefore:

grad(T) = 1/Cp ( grad(h) - sum h_i grad(Y_i) )

So, I expected this term in the enthalpy equation:

- fvm::laplacian(k/Cp, h - sum ( h_i grad(Y_i) ) )

Thank you very much.
dongchao yang likes this.
alberto_cuoci is offline   Reply With Quote

Old   October 8, 2009, 10:59
New Member
Dominik Christ
Join Date: Mar 2009
Posts: 28
Rep Power: 17
dominik_christ is on a distinguished road
Concerning your second point:

Isn't basically phi=rho*U at the faces?
Then I would understand to be phi/rho = U at the faces.
Then p*div(phi/rho) = p* div(U) which would be the effect of compression...

dominik_christ is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
Full energy equation for enthalpy hartinger OpenFOAM Running, Solving & CFD 18 April 23, 2012 07:31
dp/dt, udp/dx terms in static enthalpy equation CFDtoy Main CFD Forum 0 April 13, 2007 10:41
Energy equation and Enthalpy Martin FLUENT 1 January 25, 2001 17:11
Enthalpy equation and Cp dependency upon specis Mads Bang Main CFD Forum 1 July 30, 2000 06:43

All times are GMT -4. The time now is 09:12.