CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Enthalpy equation for multicomponent mixtures (http://www.cfd-online.com/Forums/openfoam/68347-enthalpy-equation-multicomponent-mixtures.html)

 alberto_cuoci September 16, 2009 05:11

Enthalpy equation for multicomponent mixtures

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:

div(p*phi/rho)-p*div(phi/rho)

I expected just the first term:

div(p*phi/rho)

However for multicomponent mixtures we have:

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

So, I expected this term in the enthalpy equation:

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

Thank you very much.
Alberto

 dominik_christ October 8, 2009 10:59