The relationships for calculating char particle burning rates are presented and discussed in detail by Smith [ 342]. The particle reaction rate, (kg/m -s), can be expressed as
|=||bulk diffusion coefficient (m/s)|
|=||mean reacting gas species concentration in the bulk (kg/m )|
|=||mean reacting gas species concentration at the particle surface (kg/m )|
|=||chemical reaction rate coefficient (units vary)|
|=||apparent reaction order (dimensionless)|
In Equation 14.3-1, the concentration at the particle surface, , is not known, so it should be eliminated, and the expression is recast as follows:
This equation has to be solved by an iterative procedure, with the exception of the cases when or . When , Equation 14.3-2 can be written as
In the case of , if there is a finite concentration of reactant at the particle surface, the solid depletion rate is equal to the chemical reaction rate. If there is no reactant at the surface, the solid depletion rate changes abruptly to the diffusion-controlled rate. In this case, however, FLUENT will always use the chemical reaction rate for stability reasons.
FLUENT Model Formulation
A particle undergoing an exothermic reaction in the gas phase is shown schematically in Figure 14.3.1. and are the temperatures in Equation 22.9-86.
Based on the analysis above, FLUENT uses the following equation to describe the rate of reaction of a particle surface species with the gas phase species . The reaction stoichiometry of reaction in this case is described by
and the rate of reaction is given as
|=||rate of particle surface species depletion (kg/s)|
|=||particle surface area (m )|
|=||mass fraction of surface species in the particle|
|=||effectiveness factor (dimensionless)|
|=||rate of particle surface species reaction per unit area (kg/m -s)|
|=||bulk partial pressure of the gas phase species (Pa)|
|=||diffusion rate coefficient for reaction|
|=||kinetic rate of reaction (units vary)|
|=||apparent order of reaction|
The effectiveness factor, , is related to the surface area, and can be used in each reaction in the case of multiple reactions. is given by
The kinetic rate of reaction is defined as
The rate of the particle surface species depletion for reaction order is given by
For reaction order ,
Extension for Stoichiometries with Multiple Gas Phase Reactants
When more than one gas phase reactant takes part in the reaction, the reaction stoichiometry must be extended to account for this case:
To describe the rate of reaction of a particle surface species in the presence of gas phase species , it is necessary to define the diffusion-limited species for each solid particle reaction, i.e., the species for which the concentration gradient between the bulk and the particle surface is the largest. For the rest of the species, the surface and the bulk concentrations are assumed to be equal. The concentration of the diffusion-limited species is shown as and in Figure 14.3.1, and the concentrations of all other species are denoted as . For stoichiometries with multiple gas phase reactants, the bulk partial pressure in Equations 14.3-4 and 14.3-8 is the bulk partial pressure of the diffusion-limited species, for reaction .
The kinetic rate of reaction is then defined as
|=||bulk partial pressure of gas species|
|=||reaction order in species|
When this model is enabled, the constant (Equation 14.3-6) and the effectiveness factor (Equation 14.3-4) are entered in the Reactions panel, as described in Section 14.3.2.
Reactions involving only particle surface reactants can be modeled, provided that the particle surface reactants and products exist on the same particle.
The reaction rate for this case is given by Equation 14.3-9.
Solid Decomposition Reactions
The decomposition reactions of particle surface species can be modeled.
The reaction rate for this case is given by Equations 14.3-4- 14.3-10, where the diffusion-limited species is now the gaseous product of the reaction. If there are more than one gaseous product species in the reaction, it is necessary to define the diffusion-limited species for the particle reaction as the species for which the concentration gradient between the bulk and the particle surface is the largest.
Solid Deposition Reactions
The deposition reaction of a solid species on a particle can be modeled with the following assumptions:
The theoretical analysis and Equations 14.3-4- 14.3-10 are applied for the surface reaction rate calculation, with the mass fraction of the surface species set to unity in Equations 14.3-4, 14.3-8, and 14.3-9.
In FLUENT, for the particle surface species to be deposited on a particle, a finite mass of the species must already exist in the particle. This allows for activation of the deposition reaction selectively to particular injection particles. It follows that, to initiate the solid species deposition reaction on a particle, the particle must be defined in the Set Injection Properties panel (or Set Multiple Injection Properties panel) to contain a small mass fraction of the solid species to be deposited. See Section 14.3.3 for details on defining the particle surface species mass fractions.
Gaseous Solid Catalyzed Reactions on the Particle Surface
Reactions of gaseous species catalyzed on the particle surface can also be modeled following Equations 14.3-4- 14.3-10 for the surface reaction rate calculation, with the mass fraction of the surface species set to unity in Equations 14.3-4, 14.3-8, and 14.3-9.
The catalytic particle surface reaction option is enabled in FLUENT when Particle Surface is selected as the Reaction Type in the Reactions panel and there are no solid species in the reaction stoichiometry. The solid species acting as a catalyst for the reaction is defined in the Reactions panel. The catalytic particle surface reaction will proceed only on those particles containing the catalyst species. See Section 14.3.3 for details on defining the particle surface species mass fractions.