[ANSYS, Inc. Logo] return to home search
next up previous contents index

8.3.4 Using the Non-Premixed Model with the Inert Model

To model the effect of dilution on combustion without the expense of using two mixture fractions, ANSYS FLUENT allows the introduction of an inert stream into the domain. Unlike a secondary mixture fraction, the inert does not chemically equilibrate with the primary fuel and oxidizer - instead, its composition remains constant after mixing. However the inert stream does affect the solution due to its influence on enthalpy, specific heat, and density of the mixture. The equation for conservation of inert is written as:


 \frac{\partial \rho Y_I}{\partial t} + \nabla \cdot ({\vec ... ...ft(\frac{\mu_t}{{\rm Sc}_t} \nabla \left(Y_I \right) \right) (8.3-5)


where      
  $Y_I$ = inert stream tracer
  Sc $_t$ = turbulent Schmidt number
  $\mu_t$ = turbulent viscosity
  $\rho$ = density

Equation  8.3-5 has no sources or sinks, because the problem is reduced to tracking a conserved scalar when it is assumed that the inert components have the same turbulent diffusivities.



Mixture Composition


The mixture properties are computed from the mean ( $\bar{f}$) and variance ( $\bar{f'}$) of the mixture fraction in the cell, the reaction progress variable ( $c$, when the partially premixed model is enabled), the cell enthalpy ( $H$, for non-adiabatic flows), and the inert tracer ( $Y_I$). The mixture is modeled as a blend of inert and active species, but the PDF tables need to be accessed with conditioned variables. Conditioning is necessary to take into account the volume taken up by the inert fraction, yet still be able to utilize previously built tables by straightforward lookup. The mean mixture fraction and mixture fraction variance used to access the PDF table is given by:


$\displaystyle \bar{f_I}$ $\textstyle =$ $\displaystyle \frac{\bar{f}}{1-Y_I}$ (8.3-6)
$\displaystyle \bar{f'^2_I}$ $\textstyle =$ $\displaystyle \frac{\bar{f'^2}}{(1-Y_I)^2}$ (8.3-7)

The reaction progress variable $c$ is not conditioned, however the cell enthalpy must be conditioned to account for the inert enthalpy. The inert enthalpy and active enthalpy are obtained from the following relationships:


 H(T) = (1-Y_I) H^{pdf}(T) + Y_I H^I(T) (8.3-8)

where $H$ is the enthalpy of the cell at temperature $T$, $H^{pdf}$ is the enthalpy of the active mixture-fraction stream and $H^I$ is the enthalpy of the inert stream. Here it is assumed that the inert and the active streams have the same temperature, but different enthalpies. To calculate the temperature in the cell, Equation  8.3-8 is solved for the temperature and for $H^{pdf}$, which gives the partitioning of the energy between the inert and active streams.

The inert enthalpy is defined as


 H^I(T) = \sum_{i=0}^{N_I} Y_i \int_{T_0}^{T} C_{p_i} dT (8.3-9)

where $Y_i$ refers to the mass fraction of specie $i$ defined in the inert stream, $T_0$ is the reference temperature, $C_{p_i}$ the specific heat of specie $i$, and $N_I$ is the number of inert species.

The inert and PDF enthalpies are defined further in this equation in the separate User's Guide.

Property Evaluation

The specific heat of the mixture is evaluated by mixing the inert and active streams in the following way:


 C_p(T) = (1-Y_I) C_p^{pdf}(T) + Y_I C_p^I(T) (8.3-10)

The density of the mixture is calculated by using a harmonic average of the densities of the active and inert streams, weighted by the inert tracer:


 \rho(T) = \left[ \frac{(1-Y_I)}{ \rho^{pdf}(T)} +\frac{Y_I}{\rho^I(T)} \right]^{-1} (8.3-11)

Here, the inert density ( $\rho^I$) is calculated from the ideal gas law.

For information on how to set up the inert model, see this section in the separate User's Guide.


next up previous contents index Previous: 8.3.3 Using the Non-Premixed
Up: 8.3 Restrictions and Special
Next: 8.4 The Laminar Flamelet
Release 12.0 © ANSYS, Inc. 2009-01-23