
The laws that you activate depend upon the particle type that you select. In the Set Injection Properties panel you will specify the Particle Type, and FLUENT will use a given set of heat and mass transfer laws for the chosen type. All particle types have predefined sequences of physical laws as shown in the table below:
Particle Type  Description  Laws Activated 
Inert  inert/heating or cooling  1, 6 
Droplet  heating/evaporation/boiling  1, 2, 3, 6 
Combusting  heating;
evolution of volatiles/swelling; heterogeneous surface reaction 
1, 4, 5, 6 
Multicomponent  multicomponent droplets/particles  7 
In addition to the above laws, you can define your own laws using a userdefined function. See the separate UDF Manual for information about userdefined functions.
You can also extend combusting particles to include an evaporating/boiling material by selecting Wet Combustion in the Set Injection Properties panel.
FLUENT's physical laws (Laws 1 through 6), which describe the heat and mass transfer conditions listed in this table, are explained in detail in the sections that follow.
Inert Heating or Cooling (Law 1/Law 6)
The inert heating or cooling laws (Laws 1 and 6) are applied when the particle temperature is less than the vaporization temperature that you define, , and after the volatile fraction, , of a particle has been consumed. These conditions may be written as
Law 1:
Law 6:
where is the particle temperature, is the initial mass of the particle, and is its current mass.
Law 1 is applied until the temperature of the particle/droplet reaches the vaporization temperature. At this point a noninert particle/droplet may proceed to obey one of the masstransfer laws (2, 3, 4, and/or 5), returning to Law 6 when the volatile portion of the particle/droplet has been consumed. (Note that the vaporization temperature, , is an arbitrary modeling constant used to define the onset of the vaporization/boiling/volatilization laws.)
When using Law 1 or Law 6, FLUENT uses a simple heat balance to relate the particle temperature, , to the convective heat transfer and the absorption/emission of radiation at the particle surface:
where  
=  mass of the particle (kg)  
=  heat capacity of the particle (J/kgK)  
=  surface area of the particle (m )  
=  local temperature of the continuous phase (K)  
=  convective heat transfer coefficient (W/m K)  
=  particle emissivity (dimensionless)  
=  StefanBoltzmann constant (5.67 x 10 W/m K )  
=  radiation temperature, 
Equation 22.911 assumes that there is negligible internal resistance to heat transfer, i.e., the particle is at uniform temperature throughout.
is the incident radiation in W/m :
(22.912) 
where is the radiation intensity and is the solid angle.
Radiation heat transfer to the particle is included only if you have enabled the P1 or discrete ordinates radiation model and you have activated radiation heat transfer to particles using the Particle Radiation Interaction option in the Discrete Phase Model panel.
Equation 22.911 is integrated in time using an approximate, linearized form that assumes that the particle temperature changes slowly from one time value to the next:
As the particle trajectory is computed, FLUENT integrates Equation 22.913 to obtain the particle temperature at the next time value, yielding
where is the integration time step and
(22.915) 
and
(22.916) 
FLUENT can also solve Equation 22.913 in conjunction with the equivalent mass transfer equation using a stiff coupled solver. See Section 22.11.7 for details.
The heat transfer coefficient, , is evaluated using the correlation of Ranz and Marshall [ 296, 297]:
where
=  particle diameter (m)  
=  thermal conductivity of the continuous phase (W/mK)  
Re  =  Reynolds number based on the particle diameter and  
the relative velocity (Equation 22.23)  
Pr  =  Prandtl number of the continuous phase ( ) 
Finally, the heat lost or gained by the particle as it traverses each computational cell appears as a source or sink of heat in subsequent calculations of the continuous phase energy equation. During Laws 1 and 6, particles/droplets do not exchange mass with the continuous phase and do not participate in any chemical reaction.
Droplet Vaporization (Law 2)
Law 2 is applied to predict the vaporization from a discrete phase droplet. Law 2 is initiated when the temperature of the droplet reaches the vaporization temperature, , and continues until the droplet reaches the boiling point, , or until the droplet's volatile fraction is completely consumed:
The onset of the vaporization law is determined by the setting of , a modeling parameter that has no physical significance. Note that once vaporization is initiated (by the droplet reaching this threshold temperature), it will continue to vaporize even if the droplet temperature falls below . Vaporization will be halted only if the droplet temperature falls below the dew point. In such cases, the droplet will remain in Law 2 but no evaporation will be predicted. When the boiling point is reached, the droplet vaporization is predicted by a boiling rate, Law 3, as described in a section that follows.
Mass Transfer During Law 2
During Law 2, the rate of vaporization is governed by gradient diffusion, with the flux of droplet vapor into the gas phase related to the gradient of the vapor concentration between the droplet surface and the bulk gas:
where
=  molar flux of vapor (kgmol/m s)  
=  mass transfer coefficient (m/s)  
=  vapor concentration at the droplet surface (kgmol/m )  
=  vapor concentration in the bulk gas (kgmol/m ) 
Note that FLUENT's vaporization law assumes that is positive (evaporation). If conditions exist in which is negative (i.e., the droplet temperature falls below the dew point and condensation conditions exist), FLUENT treats the droplet as inert ( ).
The concentration of vapor at the droplet surface is evaluated by assuming that the partial pressure of vapor at the interface is equal to the saturated vapor pressure, , at the particle droplet temperature, :
where is the universal gas constant.
The concentration of vapor in the bulk gas is known from solution of the transport equation for species or from the PDF lookup table for nonpremixed or partially premixed combustion calculations:
where is the local bulk mole fraction of species , is the local absolute pressure, and is the local bulk temperature in the gas. The mass transfer coefficient in Equation 22.920 is calculated from the Sherwood number correlation [ 296, 297]:
where  =  diffusion coefficient of vapor in the bulk (m /s)  
Sc  =  the Schmidt number,  
=  particle (droplet) diameter (m) 
The vapor flux given by Equation 22.920 becomes a source of species in the gas phase species transport equation, (see Section 22.14) or from the PDF lookup table for nonpremixed combustion calculations.
The mass of the droplet is reduced according to
where  =  molecular weight of species (kg/kgmol)  
=  mass of the droplet (kg)  
=  surface area of the droplet (m ) 
FLUENT can also solve Equation 22.924 in conjunction with the equivalent heat transfer equation using a stiff coupled solver. See Section 22.11.7 for details.
Defining the Vapor Pressure and Diffusion Coefficient
You must define the vapor pressure as a polynomial or piecewise linear function of temperature ( ) during the problem definition. Note that the vapor pressure definition is critical, as is used to obtain the driving force for the evaporation process (Equations 22.920 and 22.921). You should provide accurate vapor pressure values for temperatures over the entire range of possible droplet temperatures in your problem. Vapor pressure data can be obtained from a physics or engineering handbook (e.g., [ 279]).
You must also input the diffusion coefficient, , during the setup of the discrete phase material properties. Note that the diffusion coefficient inputs that you supply for the continuous phase are not used in the discrete phase model.
Heat Transfer to the Droplet
Finally, the droplet temperature is updated according to a heat balance that relates the sensible heat change in the droplet to the convective and latent heat transfer between the droplet and the continuous phase:
where  =  droplet heat capacity (J/kgK)  
=  droplet temperature (K)  
=  convective heat transfer coefficient (W/m K)  
=  temperature of continuous phase (K)  
=  rate of evaporation (kg/s)  
=  latent heat (J/kg)  
=  particle emissivity (dimensionless)  
=  StefanBoltzmann constant (5.67 x 10 W/m K )  
=  radiation temperature, , where is the radiation intensity 
Radiation heat transfer to the particle is included only if you have enabled the P1 or discrete ordinates radiation model and you have activated radiation heat transfer to particles using the Particle Radiation Interaction option in the Discrete Phase Model panel.
The heat transferred to or from the gas phase becomes a source/sink of energy during subsequent calculations of the continuous phase energy equation.
Droplet Boiling (Law 3)
Law 3 is applied to predict the convective boiling of a discrete phase droplet when the temperature of the droplet has reached the boiling temperature, , and while the mass of the droplet exceeds the nonvolatile fraction, ( ):
and
When the droplet temperature reaches the boiling point, a boiling rate equation is applied [ 188]:
where  =  heat capacity of the gas (J/kgK)  
=  droplet density (kg/m )  
=  thermal conductivity of the gas (W/mK) 
Equation 22.928 was derived assuming steady flow at constant pressure. Note that the model requires in order for boiling to occur and that the droplet remains at fixed temperature ( ) throughout the boiling law.
When radiation heat transfer is active, FLUENT uses a slight modification of Equation 22.928, derived by starting from Equation 22.925 and assuming that the droplet temperature is constant. This yields
or
Using Equation 22.917 for the Nusselt number correlation and replacing the Prandtl number term with an empirical constant, Equation 22.930 becomes
In the absence of radiation, this result matches that of Equation 22.928 in the limit that the argument of the logarithm is close to unity. FLUENT uses Equation 22.931 when radiation is active in your model and Equation 22.928 when radiation is not active. Radiation heat transfer to the particle is included only if you have enabled the P1 or discrete ordinates radiation model and you have activated radiation heat transfer to particles using the Particle Radiation Interaction option in the Discrete Phase Model panel.
The droplet is assumed to stay at constant temperature while the boiling rate is applied. Once the boiling law is entered it is applied for the duration of the particle trajectory. The energy required for vaporization appears as a (negative) source term in the energy equation for the gas phase. The evaporated liquid enters the gas phase as species , as defined by your input for the destination species (see Section 22.14).
Devolatilization (Law 4)
The devolatilization law is applied to a combusting particle when the temperature of the particle reaches the vaporization temperature, , and remains in effect while the mass of the particle, , exceeds the mass of the nonvolatiles in the particle:
and
where is the mass fraction of the evaporating/boiling material if Wet Combustion is selected (otherwise, ). As implied by Equation 22.932, the boiling point, , and the vaporization temperature, , should be set equal to each other when Law 4 is to be used. When wet combustion is active, and refer to the boiling and evaporation temperatures for the combusting material only.
FLUENT provides a choice of four devolatilization models:
Each of these models is described, in turn, below.
Choosing the Devolatilization Model
You will choose the devolatilization model when you are setting physical properties for the combustingparticle material in the Materials panel, as described in Section 22.14.2. By default, the constant rate model (Equation 22.934) will be used.
The Constant Rate Devolatilization Model
The constant rate devolatilization law dictates that volatiles are released at a constant rate [ 27]:
where  =  particle mass (kg)  
=  fraction of volatiles initially present in the particle  
=  initial particle mass (kg)  
=  rate constant (s ) 
The rate constant is defined as part of your modeling inputs, with a default value of 12 s derived from the work of Pillai [ 284] on coal combustion. Proper use of the constant devolatilization rate requires that the vaporization temperature, which controls the onset of devolatilization, be set appropriately. Values in the literature show this temperature to be about 600 K [ 27].
The volatile fraction of the particle enters the gas phase as the devolatilizing species , defined by you (see Section 22.14). Once in the gas phase, the volatiles may react according to the inputs governing the gas phase chemistry.
The Single Kinetic Rate Model
The single kinetic rate devolatilization model assumes that the rate of devolatilization is firstorder dependent on the amount of volatiles remaining in the particle [ 16]:
where  =  particle mass (kg)  
=  mass fraction of volatiles initially present in the particle  
=  mass fraction of evaporating/boiling material (if wet combustion  
is modeled)  
=  initial particle mass (kg)  
=  kinetic rate (s ) 
Note that , the fraction of volatiles in the particle, should be defined using a value slightly in excess of that determined by proximate analysis. The kinetic rate, , is defined by input of an Arrhenius type preexponential factor and an activation energy:
FLUENT uses default rate constants, and , as given in [ 16].
Equation 22.935 has the approximate analytical solution:
which is obtained by assuming that the particle temperature varies only slightly between discrete time integration steps.
FLUENT can also solve Equation 22.937 in conjunction with the equivalent heat transfer equation using a stiff coupled solver. See Section 22.11.7 for details.
The Two Competing Rates (Kobayashi) Model
FLUENT also provides the kinetic devolatilization rate expressions of the form proposed by Kobayashi [ 185]:
where and are competing rates that may control the devolatilization over different temperature ranges. The two kinetic rates are weighted to yield an expression for the devolatilization as
where  =  volatile yield up to time  
=  initial particle mass at injection  
=  yield factors  
=  ash content in the particle 
The Kobayashi model requires input of the kinetic rate parameters, , , , and , and the yields of the two competing reactions, and . FLUENT uses default values for the yield factors of 0.3 for the first (slow) reaction and 1.0 for the second (fast) reaction. It is recommended in the literature [ 185] that be set to the fraction of volatiles determined by proximate analysis, since this rate represents devolatilization at low temperature. The second yield parameter, , should be set close to unity, which is the yield of volatiles at very high temperature.
By default, Equation 22.940 is integrated in time analytically, assuming the particle temperature to be constant over the discrete time integration step. FLUENT can also solve Equation 22.940 in conjunction with the equivalent heat transfer equation using a stiff coupled solver. See Section 22.11.7 for details.
The CPD Model
In contrast to the coal devolatilization models presented above, which are based on empirical rate relationships, the chemical percolation devolatilization (CPD) model characterizes the devolatilization behavior of rapidly heated coal based on the physical and chemical transformations of the coal structure [ 110, 111, 128].
General Description
During coal pyrolysis, the labile bonds between the aromatic clusters in the coal structure lattice are cleaved, resulting in two general classes of fragments. One set of fragments has a low molecular weight (and correspondingly high vapor pressure) and escapes from the coal particle as a light gas. The other set of fragments consists of tar gas precursors that have a relatively high molecular weight (and correspondingly low vapor pressure) and tend to remain in the coal for a long period of time during typical devolatilization conditions. During this time, reattachment with the coal lattice (which is referred to as crosslinking) can occur. The high molecular weight compounds plus the residual lattice are referred to as metaplast. The softening behavior of a coal particle is determined by the quantity and nature of the metaplast generated during devolatilization. The portion of the lattice structure that remains after devolatilization is comprised of char and mineralcompoundbased ash.
The CPD model characterizes the chemical and physical processes by considering the coal structure as a simplified lattice or network of chemical bridges that link the aromatic clusters. Modeling the cleavage of the bridges and the generation of light gas, char, and tar precursors is then considered to be analogous to the chemical reaction scheme shown in Figure 22.9.2.
The variable represents the original population of labile bridges in the coal lattice. Upon heating, these bridges become the set of reactive bridges, . For the reactive bridges, two competing paths are available. In one path, the bridges react to form side chains, . The side chains may detach from the aromatic clusters to form light gas, . As bridges between neighboring aromatic clusters are cleaved, a certain fraction of the coal becomes detached from the coal lattice. These detached aromatic clusters are the heavymolecularweight tar precursors that form the metaplast. The metaplast vaporizes to form coal tar. While waiting for vaporization, the metaplast can also reattach to the coal lattice matrix (crosslinking). In the other path, the bridges react and become a char bridge, , with the release of an associated light gas product, . The total population of bridges in the coal lattice matrix can be represented by the variable , where .
Reaction Rates
Given this set of variables that characterizes the coal lattice structure during devolatilization, the following set of reaction rate expressions can be defined for each, starting with the assumption that the reactive bridges are destroyed at the same rate at which they are created ( ):
where the rate constants for bridge breaking and gas release steps, and , are expressed in Arrhenius form with a distributed activation energy:
where , , and are, respectively, the preexponential factor, the activation energy, and the distributed variation in the activation energy, is the universal gas constant, and is the temperature. The ratio of rate constants, , is set to 0.9 in this model based on experimental data.
Mass Conservation
The following mass conservation relationships are imposed:
where is the fraction of broken bridges ( ). The initial conditions for this system are given by the following:
where is the initial fraction of char bridges, is the initial fraction of bridges in the coal lattice, and is the initial fraction of labile bridges in the coal lattice.
Fractional Change in the Coal Mass
Given the set of reaction equations for the coal structure parameters, it is necessary to relate these quantities to changes in coal mass and the related release of volatile products. To accomplish this, the fractional change in the coal mass as a function of time is divided into three parts: light gas ( ), tar precursor fragments ( ), and char ( ). This is accomplished by using the following relationships, which are obtained using percolation lattice statistics:
The variables , , , and are the statistical relationships related to the cleaving of bridges based on the percolation lattice statistics, and are given by the following equations:
is the ratio of bridge mass to site mass, , where
where and are the side chain and cluster molecular weights respectively. is the lattice coordination number, which is determined from solidstate nuclear magnetic esonance (NMR) measurements related to coal structure parameters, and is the root of the following equation in (the total number of bridges in the coal lattice matrix):
In accounting for mass in the metaplast (tar precursor fragments), the part that vaporizes is treated in a manner similar to flash vaporization, where it is assumed that the finite fragments undergo vapor/liquid phase equilibration on a time scale that is rapid with respect to the bridge reactions. As an estimate of the vapor/liquid that is present at any time, a vapor pressure correlation based on a simple form of Raoult's Law is used. The vapor pressure treatment is largely responsible for predicting pressuredependent devolatilization yields. For the part of the metaplast that reattaches to the coal lattice, a crosslinking rate expression given by the following equation is used:
where is the amount of mass reattaching to the matrix, is the amount of mass in the tar precursor fragments (metaplast), and and are rate expression constants.
CPD Inputs
Given the set of equations and corresponding rate constants introduced for the CPD model, the number of constants that must be defined to use the model is a primary concern. For the relationships defined previously, it can be shown that the following parameters are coal independent [ 110]:
These constants are included in the submodel formulation and are not input or modified during problem setup.
There are an additional five parameters that are coalspecific and must be specified during the problem setup:
The first four of these are coal structure quantities that are obtained from NMR experimental data. The last quantity, representing the char bridges that either exist in the parent coal or are formed very early in the devolatilization process, is estimated based on the coal rank. These quantities are entered in the Materials panel as described in Section 22.14.2. Values for the coaldependent parameters for a variety of coals are listed in Table 22.9.1.
Coal Type  
Zap (AR)  3.9  .63  277  40  .20 
Wyodak (AR)  5.6  .55  410  42  .14 
Utah (AR)  5.1  .49  359  36  0 
Ill6 (AR)  5.0  .63  316  27  0 
Pitt8 (AR)  4.5  .62  294  24  0 
Stockton (AR)  4.8  .69  275  20  0 
Freeport (AR)  5.3  .67  302  17  0 
Pocahontas (AR)  4.4  .74  299  14  .20 
Blue (Sandia)  5.0  .42  410  47  .15 
Rose (AFR)  5.8  .57  459  48  .10 
1443 (lignite, ACERC)  4.8  .59  297  36  .20 
1488 (subbituminous, ACERC)  4.7  .54  310  37  .15 
1468 (anthracite, ACERC)  4.7  .89  656  12  .25 
Particle Swelling During Devolatilization
The particle diameter changes during devolatilization according to the swelling coefficient, , which is defined by you and applied in the following relationship:
where  =  particle diameter at the start of devolatilization  
=  current particle diameter 
The term is the ratio of the mass that has been devolatilized to the total volatile mass of the particle. This quantity approaches a value of 1.0 as the devolatilization law is applied. When the swelling coefficient is equal to 1.0, the particle diameter stays constant. When the swelling coefficient is equal to 2.0, the final particle diameter doubles when all of the volatile component has vaporized, and when the swelling coefficient is equal to 0.5 the final particle diameter is half of its initial diameter.
Heat Transfer to the Particle During Devolatilization
Heat transfer to the particle during the devolatilization process includes contributions from convection, radiation (if active), and the heat consumed during devolatilization:
Radiation heat transfer to the particle is included only if you have enabled the P1 or discrete ordinates radiation model and you have activated radiation heat transfer to particles using the Particle Radiation Interaction option in the Discrete Phase Model panel.
By default, Equation 22.966 is solved analytically, by assuming that the temperature and mass of the particle do not change significantly between time steps:
(22.967) 
where
(22.968) 
and
(22.969) 
FLUENT can also solve Equation 22.966 in conjunction with the equivalent mass transfer equation using a stiff coupled solver. See Section 22.11.7 for details.
Surface Combustion (Law 5)
After the volatile component of the particle is completely evolved, a surface reaction begins which consumes the combustible fraction, , of the particle. Law 5 is thus active (for a combusting particle) after the volatiles are evolved:
and until the combustible fraction is consumed:
When the combustible fraction, , has been consumed in Law 5, the combusting particle may contain residual "ash'' that reverts to the inert heating law, Law 6 (described previously).
With the exception of the multiple surface reactions model, the surface combustion law consumes the reactive content of the particle as governed by the stoichiometric requirement, , of the surface "burnout'' reaction:
where is defined in terms of mass of oxidant per mass of char, and the oxidant and product species are defined in the Set Injection Properties panel.
FLUENT provides a choice of four heterogeneous surface reaction rate models for combusting particles:
Each of these models is described in detail below. You will choose the surface combustion model when you are setting physical properties for the combustingparticle material in the Materials panel, as described in Section 22.14.2. By default, the diffusionlimited rate model will be used.
The DiffusionLimited Surface Reaction Rate Model
The diffusionlimited surface reaction rate model which is the default model in FLUENT, assumes that the surface reaction proceeds at a rate determined by the diffusion of the gaseous oxidant to the surface of the particle:
where  =  diffusion coefficient for oxidant in the bulk (m /s)  
=  local mass fraction of oxidant in the gas  
=  gas density (kg/m )  
=  stoichiometry of Equation 22.972 
Equation 22.973 is derived from the model of Baum and Street [ 27] with the kinetic contribution to the surface reaction rate ignored. The diffusionlimited rate model assumes that the diameter of the particles does not change. Since the mass of the particles is decreasing, the effective density decreases, and the char particles become more porous.
The Kinetic/Diffusion Surface Reaction Rate Model
The kinetic/diffusionlimited rate model assumes that the surface reaction rate is determined either by kinetics or by a diffusion rate. FLUENT uses the model of Baum and Street [ 27] and Field [ 106], in which a diffusion rate coefficient
and a kinetic rate
are weighted to yield a char combustion rate of
where is the surface area of the droplet ( ), is the partial pressure of oxidant species in the gas surrounding the combusting particle, and the kinetic rate, , incorporates the effects of chemical reaction on the internal surface of the char particle (intrinsic reaction) and pore diffusion. In FLUENT, Equation 22.976 is recast in terms of the oxidant mass fraction, , as
The particle size is assumed to remain constant in this model while the density is allowed to decrease.
When this model is enabled, the rate constants used in Equations 22.974 and 22.975 are entered in the Materials panel, as described in Section 22.14.
The Intrinsic Model
The intrinsic model in FLUENT is based on Smith's model [ 342], assuming the order of reaction is equal to unity. Like the kinetic/diffusion model, the intrinsic model assumes that the surface reaction rate includes the effects of both bulk diffusion and chemical reaction (see Equation 22.977). The intrinsic model uses Equation 22.974 to compute the diffusion rate coefficient, , but the chemical rate, , is explicitly expressed in terms of the intrinsic chemical and pore diffusion rates:
is the effectiveness factor, or the ratio of the actual combustion rate to the rate attainable if no pore diffusion resistance existed [ 198]:
(22.979) 
where is the Thiele modulus:
is the density of oxidant in the bulk gas (kg/m ) and is the effective diffusion coefficient in the particle pores. Assuming that the pore size distribution is unimodal and the bulk and Knudsen diffusion proceed in parallel, is given by
where is the bulk molecular diffusion coefficient and is the porosity of the char particle:
(22.982) 
and are, respectively, the apparent and true densities of the pyrolysis char.
(in Equation 22.981) is the tortuosity of the pores. The default value for in FLUENT is , which corresponds to an average intersecting angle between the pores and the external surface of 45 [ 198].
is the Knudsen diffusion coefficient:
where is the particle temperature and is the mean pore radius of the char particle, which can be measured by mercury porosimetry. Note that macropores ( Å) dominate in lowrank chars while micropores ( Å) dominate in highrank chars [ 198].
(in Equations 22.978 and 22.980) is the specific internal surface area of the char particle, which is assumed in this model to remain constant during char combustion. Internal surface area data for various pyrolysis chars can be found in [ 341]. The mean value of the internal surface area during char combustion is higher than that of the pyrolysis char [ 198]. For example, an estimated mean value for bituminous chars is 300 m /g [ 53].
(in Equations 22.978 and 22.980) is the intrinsic reactivity, which is of Arrhenius form:
where the preexponential factor and the activation energy can be measured for each char. In the absence of such measurements, the default values provided by FLUENT (which are taken from a least squares fit of data of a wide range of porous carbons, including chars [ 341]) can be used.
To allow a more adequate description of the char particle size (and hence density) variation during combustion, you can specify the burning mode , relating the char particle diameter to the fractional degree of burnout (where ) by [ 340]
where is the char particle mass and the subscript zero refers to initial conditions (i.e., at the start of char combustion). Note that where the limiting values 0 and correspond, respectively, to a constant size with decreasing density (zone 1) and a decreasing size with constant density (zone 3) during burnout. In zone 2, an intermediate value of , corresponding to a decrease of both size and density, has been found to work well for a variety of chars [ 340].
When this model is enabled, the rate constants used in Equations 22.974, 22.978, 22.980, 22.981, 22.983, 22.984, and 22.985 are entered in the Materials panel, as described in Section 22.14.
The Multiple Surface Reactions Model
Modeling multiple particle surface reactions follows a pattern similar to the wall surface reaction models, where the surface species is now a "particle surface species''. For the mixture material defined in the Species Model panel, the particle surface species can be depleted or produced by the stoichiometry of the particle surface reaction (defined in the Reactions panel). The particle surface species constitutes the reactive char mass of the particle, hence, if a particle surface species is depleted, the reactive "char'' content of the particle is consumed, and in turn, when a surface species is produced, it is added to the particle "char'' mass. Any number of particle surface species and any number of particle surface reactions can be defined for any given combusting particle.
Multiple injections can be accommodated, and combusting particles reacting according to the multiple surface reactions model can coexist in the calculation, with combusting particles following other char combustion laws. The model is based on oxidation studies of char particles, but it is also applicable to gassolid reactions in general, not only to char oxidation reactions.
See Section 14.3 for information about particle surface reactions.
Limitations
Note the following limitations of the multiple surface reactions model:
Heat and Mass Transfer During Char Combustion
The surface reaction consumes the oxidant species in the gas phase; i.e., it supplies a (negative) source term during the computation of the transport equation for this species. Similarly, the surface reaction is a source of species in the gas phase: the product of the heterogeneous surface reaction appears in the gas phase as a userselected chemical species. The surface reaction also consumes or produces energy, in an amount determined by the heat of reaction defined by you.
The particle heat balance during surface reaction is
where is the heat released by the surface reaction. Note that only a portion ( ) of the energy produced by the surface reaction appears as a heat source in the gasphase energy equation: the particle absorbs a fraction of this heat directly. For coal combustion, it is recommended that be set at 1.0 if the char burnout product is CO and 0.3 if the char burnout product is CO [ 39].
Radiation heat transfer to the particle is included only if you have enabled the P1 or discrete ordinates radiation model and you have activated radiation heat transfer to particles using the Particle Radiation Interaction option in the Discrete Phase Model panel.
By default, Equation 22.986 is solved analytically, by assuming that the temperature and mass of the particle do not change significantly between time steps. FLUENT can also solve Equation 22.986 in conjunction with the equivalent mass transfer equation using a stiff coupled solver. See Section 22.11.7 for details.
Multicomponent Particle Definition (Law 7)
Multicomponent particles are described in FLUENT as a mixture of species within droplets/particles. The particle mass is the sum of the masses of the components
(22.987) 
The density of the particle can be either constant, or volumeaveraged:
(22.988) 
For particles containing more than one component it is difficult to assign the whole particle to one process like boiling or heating. Therefore it can be only modeled by a law integrating all processes of relevance in one equation. The source terms for temperature and component mass are the sum of the sources from the partial processes:
(22.990) 
The equation for the particle temperature consists of terms for radiation, convective heating (Equation 22.911) and vaporization. Radiation heat transfer to the particle is included only if you have enabled P1 or DiscreteOrdinates (DO) radiation and you have activated radiation heat transfer to the particles using the Particle Radiation Interaction option in the Discrete Phase Model panel.
The mass of the particle components is only influenced by the vaporization (Equation 22.920), where is the molecular weight of species . The mass transfer coefficient of component is calculated from the Sherwood correlation (Equation 22.923). The concentration of vapor at the particle surface depends on the saturation pressure of the component.
Raoult's Law
The correlation between the vapor concentration of a species over the surface and its mole fraction in the condensed phase (under the assumption of an ideal gas law) is described by Raoult's law:
(22.991) 
You can define your own law other than Raoult's Law for vapor concentration at the particle surface using a userdefined function.
Click here to go to the FLUENT UDF Manual for details.