For the sectional method the coagulation rate expression is simply that given by Equation 19–15 through Equation 19–18 . The collision kernels as given by equation Equation 19–50 , Equation 19–51 , and Equation 19–52 are directly used by Ansys Chemkin, with sizes i and j corresponding to representative particles for the sections under consideration.
It has been postulated [38] that whether a collision between two particles results in the formation of an aggregate depends on the size of the particles involved in the collision. In other words, the size of the particles involved in a collision influences whether the particles will stick together or not. Models have been developed to estimate the lower and upper limits of such "collision efficiency." Note that these models are applicable only to the free-molecular regime.
Assuming that a moving particle does not experience friction (due to collisions) from the surrounding gas, the lower limit for the free-molecular regime is established using the Lennard-Jones potential. It is written as
(19–68) |
In Equation 19–68, is the Lennard Jones potential well depth and is the dimensionless surface temperature. The latter accounts for the additional energy involved in a collision due to the potential and is given by
(19–69) |
The potential well depth is a function of particle size and is computed using the Hamaker () constant as
(19–70) |
where , , and are functions of particle size (diameter ) and the minimum separation distance ().
(19–71) |
The lower limit given by Equation 19–68 is assumed to be valid for "sufficiently small" particles. On the other hand, for "sufficiently large" particles, it is argued that the particles moving away after a collision experience friction and there is more chance that they get trapped in the potential energy well. The upper limit for the collision efficiency is thus considered to be unity.
As suggested in the literature [17], Ansys Chemkin uses a mixing function to blend the upper and lower limits. This mixing function and the effective collision efficiency are then expressed as
(19–72) |
(19–73) |
As mentioned before, the upper limit used in Equation 19–73 is taken to be unity.
In Equation 19–72, and are the parameters that control the blending effect of the mixing function. Note that the mixing function approaching 1 means the upper limit of the collision efficiency is the dominant part in the blended effective value whereas its value of 0 means the lower limit is dominant. The latter implies that the value of should be such that for small particles the Tanh term becomes -1. Thus, one may write
(19–74) |
In Equation 19–74, is a constant and is the diameter at which the mixing function approaches zero. A small integer value of (around 5) is then sufficient to yield such that the Tanh function approaches -1 as the particle diameter approaches . It should be noted that Equation 19–74 is just one of the many possible choices of setting . Noting that the value of the mixing function is 0.5 for = , represents the threshold on either side of which one of the limits of the collision efficiency starts to dominate. For a given , the larger the value of the quicker is the either limit reached by the mixing function and hence the effective collision efficiency. In Ansys Chemkin, the default value of is chosen to be 7 whereas is the particle diameter corresponding to the smallest section as specified by the user.
For simulations involving sooting flames [17], has been correlated with the maximum flame temperature and the corresponding values are of the order of 10 nm. Ansys Chemkin expects that a value of is given by the user. The default values of the Hamaker constant and the minimum separation distance used by Chemkin are also based on those used in [17]. These default values are only a guideline and may not be applicable to particulate systems other than soot.
The final general expression for the coagulation kernel in the free molecular regime involves multiplication of Equation 19–50 by the collision efficiency as calculated in Equation 19–73 and by a constant factor () accounting for the van der Waals enhancement factor. Thus,
(19–75) |