A detailed reaction mechanism for sulfur oxidation has been proposed by Kramlich [ 187]. The mechanism consists of 20 reversible reactions and includes 12 species (S, , SH, SO, , , H, , OH, , O and ). However, such a detailed mechanism is not practical for CFD applications unless used in a flamelet based model. An attempt was made therefore to reduce the number of reactions and species to a manageable level while predicting the major species concentration levels within acceptable tolerance. The need to retain the minor intermediate species SH and SO has been the major constraint and as such, an eight-step mechanism has been identified using PSR (a Fortran program for modeling well stirred reactors with gas reactions) and SENKIN (a Fortran program for predicting homogeneous gas phase chemical kinetics). The rationale behind the selection of an eight-step mechanism was to remove species S and from the detailed mechanism. Table 20.2.1 lists the reduced mechanism, which has been developed in-house, with the modified rate constants. For reduction calculations O and OH concentrations have been calculated through partial equilibrium assumptions based on and concentrations respectively. was used as the dilutant. Since each reaction of the eight-step reduced mechanism is reversible, for each adjacent pair of reactions given in Table 20.2.1, the second reaction is in fact the reverse reaction of the first.
The reduced mechanism given in Table 20.2.1 closely follows the concentration levels but slightly overpredicts the concentrations at temperatures below 1500 K. Above 1500 K, both mechanisms are in close agreement for and concentration predictions. However, SO and SH are not well correlated by the reduced mechanism when compared against the predictions using the original detailed mechanism.
A major concern in these mechanisms is the presence of H radical and the method in which to calculate its concentration in the flow field. At present, the concentration of H radical is assumed to be proportional to the O radical concentration, which can be evaluated from one of the existing methods in FLUENT; viz. Partial Equilibrium (Section 20.1.3) or Equilibrium (Section 20.1.3). The user is then given the option to vary the proportionality constant. Although this assumption is open to debate, the lack of simple relation to calculate the H radical concentration in a flame has prompted the present choice.
Present implementation allows the user to either include or remove SO from the calculations. Also, depending on the form of fuel sulfur release (e.g., or SO ) the species may or may not be present for the calculation. The user is also given the extended option of partitioning the intermediate fuel sulfur species to and . However, there is no literature to guide the user on how to select a correct partition fraction.
A is in m /gmol-s, E is J/gmol (assumed 1 cal = 4.18585 J), A units for the thirteenth reaction is m /gmol -s, and A units for the fifteenth reaction is m /gmol -s.
In addition, the following two reactions were included in FLUENT to complete the SOx mechanism, with the rate constants taken from Hunter's work [ 150].
M = argon, nitrogen, oxygen
k = x exp(+4185.85/RT) m /gmol /sec
where R = 8.313 J/gmol-K
k = x exp(-39765.575/RT) m /gmol/sec
The reverse rate of Equation 20.2-7 was determined through the equilibrium constant for that equation.