|September 16, 2009, 05:11||
Enthalpy equation for multicomponent mixtures
Join Date: Aug 2009
Posts: 3Rep Power: 5
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::Sp(fvc::div(phi), h)
- fvm::laplacian(turbulence->alphaEff(), h)
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.
|October 8, 2009, 10:59||
Join Date: Mar 2009
Posts: 12Rep Power: 6
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...
|Thread||Thread Starter||Forum||Replies||Last Post|
|Full energy equation for enthalpy||hartinger||OpenFOAM Running, Solving & CFD||18||April 23, 2012 07:31|
|Calculation of the Governing Equations||Mihail||CFX||5||July 25, 2008 17:29|
|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|