Species Transport Equations
When you choose to solve conservation equations for chemical species, FLUENT predicts the local mass fraction of each species, , through the solution of a convection-diffusion equation for the th species. This conservation equation takes the following general form:
where is the net rate of production of species by chemical reaction (described later in this section) and is the rate of creation by addition from the dispersed phase plus any user-defined sources. An equation of this form will be solved for species where is the total number of fluid phase chemical species present in the system. Since the mass fraction of the species must sum to unity, the th mass fraction is determined as one minus the sum of the solved mass fractions. To minimize numerical error, the th species should be selected as that species with the overall largest mass fraction, such as N when the oxidizer is air.
Mass Diffusion in Laminar Flows
In Equation 14.1-1, is the diffusion flux of species , which arises due to concentration gradients. By default, FLUENT uses the dilute approximation, under which the diffusion flux can be written as
Here is the diffusion coefficient for species in the mixture.
For certain laminar flows, the dilute approximation may not be acceptable, and full multicomponent diffusion is required. In such cases, the Maxwell-Stefan equations can be solved; see Section 8.9.2 for details.
Mass Diffusion in Turbulent Flows
In turbulent flows, FLUENT computes the mass diffusion in the following form:
where is the turbulent Schmidt number ( where is the turbulent viscosity and is the turbulent diffusivity). The default is 0.7. Note that turbulent diffusion generally overwhelms laminar diffusion, and the specification of detailed laminar diffusion properties in turbulent flows is generally not warranted.
Treatment of Species Transport in the Energy Equation
For many multicomponent mixing flows, the transport of enthalpy due to species diffusion
can have a significant effect on the enthalpy field and should not be neglected. In particular, when the Lewis number
for any species is far from unity, neglecting this term can lead to significant errors. FLUENT will include this term by default. In Equation 14.1-4, is the thermal conductivity.
Diffusion at Inlets
For the pressure-based solver in FLUENT, the net transport of species at inlets consists of both convection and diffusion components. (For the density-based solvers, only the convection component is included.) The convection component is fixed by the inlet species mass fraction specified by you. The diffusion component, however, depends on the gradient of the computed species field at the inlet. Thus the diffusion component (and therefore the net inlet transport) is not specified a priori. See Section 14.1.5 for information about specifying the net inlet transport of species.
The Generalized Finite-Rate Formulation for Reaction Modeling
The reaction rates that appear as source terms in Equation 14.1-1 are computed in FLUENT by one of three models:
The generalized finite-rate formulation is suitable for a wide range of applications including laminar or turbulent reaction systems, and combustion systems with premixed, non-premixed, or partially-premixed flames.
The Laminar Finite-Rate Model
The laminar finite-rate model computes the chemical source terms using Arrhenius expressions, and ignores the effects of turbulent fluctuations. The model is exact for laminar flames, but is generally inaccurate for turbulent flames due to highly non-linear Arrhenius chemical kinetics. The laminar model may, however, be acceptable for combustion with relatively slow chemistry and small turbulent fluctuations, such as supersonic flames.
The net source of chemical species due to reaction is computed as the sum of the Arrhenius reaction sources over the reactions that the species participate in:
where is the molecular weight of species and is the Arrhenius molar rate of creation/destruction of species in reaction . Reaction may occur in the continuous phase between continuous-phase species only, or at wall surfaces resulting in the surface deposition or evolution of a continuous-phase species.
Consider the th reaction written in general form as follows:
|=||number of chemical species in the system|
|=||stoichiometric coefficient for reactant in reaction|
|=||stoichiometric coefficient for product in reaction|
|=||symbol denoting species|
|=||forward rate constant for reaction|
|=||backward rate constant for reaction|
Equation 14.1-6 is valid for both reversible and non-reversible reactions. (Reactions in FLUENT are non-reversible by default.) For non-reversible reactions, the backward rate constant, , is simply omitted.
The summations in Equation 14.1-6 are for all chemical species in the system, but only species that appear as reactants or products will have non-zero stoichiometric coefficients. Hence, species that are not involved will drop out of the equation.
For a non-reversible reaction (that is, the Include Backward Reaction button is disabled), the molar rate of creation/destruction of species in reaction ( in Equation 14.1-5) is given by
|=||molar concentration of species in reaction (kgmol/m )|
|=||rate exponent for reactant species in reaction|
|=||rate exponent for product species in reaction|
For a reversible reaction, the molar rate of creation/destruction of species in reaction is given by
Note that the rate exponent for the reverse reaction part in Equation 14.1-8 is always the product species stoichiometric coefficient ( ).
See Section 14.1.4 for information about inputting the stoichiometric coefficients and rate exponents for both global forward (non-reversible) reactions and elementary (reversible) reactions.
represents the net effect of third bodies on the reaction rate. This term is given by
where is the third-body efficiency of the th species in the th reaction. By default, FLUENT does not include third-body effects in the reaction rate calculation. You can, however, opt to include the effect of third-body efficiencies if you have data for them.
The forward rate constant for reaction , , is computed using the Arrhenius expression
|=||pre-exponential factor (consistent units)|
|=||temperature exponent (dimensionless)|
|=||activation energy for the reaction (J/kgmol)|
|=||universal gas constant (J/kgmol-K)|
You (or the database) will provide values for , , , , , , , and, optionally, during the problem definition in FLUENT.
If the reaction is reversible, the backward rate constant for reaction , , is computed from the forward rate constant using the following relation:
where is the equilibrium constant for the th reaction, computed from
where denotes atmospheric pressure (101325 Pa). The term within the exponential function represents the change in Gibbs free energy, and its components are computed as follows:
where and are the standard-state entropy and standard-state enthalpy (heat of formation). These values are specified in FLUENT as properties of the mixture material.
FLUENT can use one of three methods to represent the rate expression in pressure-dependent (or pressure fall-off) reactions. A "fall-off'' reaction is one in which the temperature and pressure are such that the reaction occurs between Arrhenius high-pressure and low-pressure limits, and thus is no longer solely dependent on temperature.
There are three methods of representing the rate expressions in this fall-off region. The simplest one is the Lindemann [ 212] form. There are also two other related methods, the Troe method [ 122] and the SRI method [ 357], that provide a more accurate description of the fall-off region.
Arrhenius rate parameters are required for both the high- and low-pressure limits. The rate coefficients for these two limits are then blended to produce a smooth pressure-dependent rate expression. In Arrhenius form, the parameters for the high-pressure limit ( ) and the low-pressure limit ( ) are as follows:
The net rate constant at any pressure is then taken to be
where is defined as
and is the concentration of the bath gas, which can include third-body efficiencies. If the function in Equation 14.1-17 is unity, then this is the Lindemann form. FLUENT provides two other forms to describe , namely the Troe method and the SRI method.
In the Troe method, is given by
The parameters , , , and are specified as inputs.
In the SRI method, the blending function is approximated as
In addition to the three Arrhenius parameters for the low-pressure limit ( ) expression, you must also supply the parameters , , , , and in the expression.
| Chemical kinetic mechanisms are usually highly non-linear and form a set of stiff coupled equations. See Section
14.1.7 for solution procedure guidelines. Also, if you have a chemical mechanism in CHEMKIN [
178] format, you can import this mechanism into
FLUENT as described in Section
The Eddy-Dissipation Model
Most fuels are fast burning, and the overall rate of reaction is controlled by turbulent mixing. In non-premixed flames, turbulence slowly convects/mixes fuel and oxidizer into the reaction zones where they burn quickly. In premixed flames, the turbulence slowly convects/mixes cold reactants and hot products into the reaction zones, where reaction occurs rapidly. In such cases, the combustion is said to be mixing-limited, and the complex, and often unknown, chemical kinetic rates can be safely neglected.
FLUENT provides a turbulence-chemistry interaction model, based on the work of Magnussen and Hjertager [ 229], called the eddy-dissipation model. The net rate of production of species due to reaction , , is given by the smaller (i.e., limiting value) of the two expressions below:
|where||is the mass fraction of any product species,|
|is the mass fraction of a particular reactant,|
|is an empirical constant equal to 4.0|
|is an empirical constant equal to 0.5|
In Equations 14.1-26 and 14.1-27, the chemical reaction rate is governed by the large-eddy mixing time scale, , as in the eddy-breakup model of Spalding [ 350]. Combustion proceeds whenever turbulence is present ( ), and an ignition source is not required to initiate combustion. This is usually acceptable for non-premixed flames, but in premixed flames, the reactants will burn as soon as they enter the computational domain, upstream of the flame stabilizer. To remedy this, FLUENT provides the finite-rate/eddy-dissipation model, where both the Arrhenius (Equation 14.1-8), and eddy-dissipation (Equations 14.1-26 and 14.1-27) reaction rates are calculated. The net reaction rate is taken as the minimum of these two rates. In practice, the Arrhenius rate acts as a kinetic "switch'', preventing reaction before the flame holder. Once the flame is ignited, the eddy-dissipation rate is generally smaller than the Arrhenius rate, and reactions are mixing-limited.
FLUENT allows multi-step reaction mechanisms (number of reactions
) with the eddy-dissipation and finite-rate/eddy-dissipation models, these will likely produce incorrect solutions. The reason is that multi-step chemical mechanisms are based on Arrhenius rates, which differ for each reaction. In the eddy-dissipation model, every reaction has the same, turbulent rate, and therefore the model should be used only for one-step (reactant
product), or two-step (reactant
product) global reactions. The model cannot predict kinetically controlled species such as radicals. To incorporate multi-step chemical kinetic mechanisms in turbulent flows, use the EDC model (described below).
| The eddy-dissipation model requires products to initiate reaction (see Equation
14.1-27). When you initialize the solution for steady flows,
FLUENT sets all species mass fractions to a maximum of the user specified initial value and 0.01. This is usually sufficient to start the reaction. However, if you converge a mixing solution first, where all product mass fractions are zero, you may then have to patch products into the reaction zone to ignite the flame. See Section
14.1.7 for details.
The Eddy-Dissipation Model for LES
When the LES turbulence model is used, the turbulent mixing rate, in Equations 14.1-26 and 14.1-27, is replaced by the subgrid-scale mixing rate. This is calculated as
|=||subgrid-scale mixing rate (s )|
|=||= strain rate tensor (s )|
The Eddy-Dissipation-Concept (EDC) Model
The eddy-dissipation-concept (EDC) model is an extension of the eddy-dissipation model to include detailed chemical mechanisms in turbulent flows [ 228]. It assumes that reaction occurs in small turbulent structures, called the fine scales. The length fraction of the fine scales is modeled as [ 127]
where denotes fine-scale quantities and
|=||volume fraction constant = 2.1377|
The volume fraction of the fine scales is calculated as . Species are assumed to react in the fine structures over a time scale
where is a time scale constant equal to 0.4082.
In FLUENT, combustion at the fine scales is assumed to occur as a constant pressure reactor, with initial conditions taken as the current species and temperature in the cell. Reactions proceed over the time scale , governed by the Arrhenius rates of Equation 14.1-8, and are integrated numerically using the ISAT algorithm [ 290]. ISAT can accelerate the chemistry calculations by two to three orders of magnitude, offering substantial reductions in run-times. Details about the ISAT algorithm may be found in Sections 18.2.4 and 18.2.5. ISAT is very powerful, but requires some care. See Section 18.3.9 for details on using ISAT efficiently.
The source term in the conservation equation for the mean species , Equation 14.1-1, is modeled as
where is the fine-scale species mass fraction after reacting over the time .
The EDC model can incorporate detailed chemical mechanisms into turbulent reacting flows. However, typical mechanisms are invariably stiff and their numerical integration is computationally costly. Hence, the model should be used only when the assumption of fast chemistry is invalid, such as modeling the slow CO burnout in rapidly quenched flames, or the NO conversion in selective non-catalytic reduction (SNCR).
See Section 14.1.7 for guidelines on obtaining a solution using the EDC model.