Consider the th wall surface reaction written in general form as follows:
where , , and represent the gas phase species, the bulk (or solid) species, and the surface-adsorbed (or site) species, respectively. , , and are the total numbers of these species. , , and are the stoichiometric coefficients for each reactant species , and , , and are the stoichiometric coefficients for each product species . is the overall reaction rate constant.
The summations in Equation 14.2-1 are for all chemical species in the system, but only species involved as reactants or products will have non-zero stoichiometric coefficients. Hence, species that are not involved will drop out of the equation.
The rate of the th reaction is
where represents molar concentrations on the wall. It is assumed that the reaction rate does not depend on concentrations of the bulk (solid) species. From this, the net molar rate of production or consumption of each species is given by
The forward rate constant for reaction ( ) is computed using the Arrhenius expression. For example,
|where||=||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 .
Wall Surface Reaction Boundary Conditions
As shown in Equations 14.2-2- 14.2-5, the goal of surface reaction modeling is to compute concentrations of gas species and site species at the wall; i.e., and . Assuming that, on a reacting surface , the mass flux of each gas species is balanced with its rate of production/consumption, then
The mass fraction is related to concentration by
is the net rate of mass deposition or etching as a result of surface reaction; i.e.,
is the site species concentration at the wall, and is defined as
where is the site density and is the site coverage of species .
Using Equations 14.2-7 and 14.2-8, expressions can be derived for the mass fraction of species at the wall and for the net rate of creation of species per unit area. These expressions are used in FLUENT to compute gas phase species concentrations, and if applicable, site coverages, at reacting surfaces using a point-by-point coupled stiff solver.
Including Mass Transfer To Surfaces in Continuity
In the surface reaction boundary condition described above, the effects of the wall normal velocity or bulk mass transfer to the wall are not included in the computation of species transport. The momentum of the net surface mass flux from the surface is also ignored because the momentum flux through the surface is usually small in comparison with the momentum of the flow in the cells adjacent to the surface. However, you can include the effect of surface mass transfer in the continuity equation by activating the Mass Deposition Source option in the Species Model panel.
Wall Surface Mass Transfer Effects in the Energy Equation
Species diffusion effects in the energy equation due to wall surface reactions are included in the normal species diffusion term described in Section 14.1.1.
If you are using the pressure-based solver, you can neglect this term by turning off the Diffusion Energy Source option in the Species Model panel. (For the density-based solvers, this term is always included; you cannot turn it off.) Neglecting the species diffusion term implies that errors may be introduced to the prediction of temperature in problems involving mixing of species with significantly different heat capacities, especially for components with a Lewis number far from unity. While the effect of species diffusion should go to zero at Le = 1, you may see subtle effects due to differences in the numerical integration in the species and energy equations.
Modeling the Heat Release Due to Wall Surface Reactions
The heat release due to a wall surface reaction is, by default, ignored by FLUENT. You can, however, choose to include the heat of surface reaction by activating the Heat of Surface Reactions option in the Species Model panel and setting appropriate formation enthalpies in the Materials panel.
Slip Boundary Formulation for Low-Pressure Gas Systems
Most semiconductor fabrication devices operate far below atmospheric pressure, typically only a few millitorrs. At such low pressures, the fluid flow is in the slip regime and the normally used no-slip boundary conditions for velocity and temperature are no longer valid.
Knudsen number , Kn, which is defined as the ratio of mean free path to a characteristic length scale of the system, is used to describe various flow regimes. Since the mean free path increases as the pressure is lowered, the high end of Kn values represents free molecular flow and the low end the continuum regime. The range in between these two extremes is called the slip regime ( ).
In the slip regime, the gas-phase velocity at a solid surface differs from the velocity at which the wall moves, and the gas temperature at the surface differs from the wall temperature. Maxwell's models are adopted for these physical phenomena in FLUENT for their simplicity and effectiveness.
Here, and represents the velocity component that is parallel and normal to the wall, respectively. The subscripts , and indicate gas, wall and cell-center velocities. is the distance from cell center to the wall. is the momentum accommodation coefficient of the gas mixture and its value is calculated as mass-fraction weighted average of each gas species in the system.
The mean free path, , is computed as follows:
is the Lennard-Jones characteristic length of species . is the Boltzmann constant, .
Equations 14.2-12 and 14.2-13 indicate that while the gas velocity component normal to the wall is the same as the wall normal velocity, the tangential components slip. The values lie somewhere between the cell-center and the wall values. These two equations can be combined to give a generalized formulation:
is the thermal accommodation coefficient of the gas mixture and is calculated as .
| The low-pressure slip boundary formulation is available only with the pressure-based solver.