Abstract and Applied Analysis

Volume 2013 (2013), Article ID 490123, 6 pages

http://dx.doi.org/10.1155/2013/490123

Research Article

The Fundamental Aspects of TEMOM Model for Particle Coagulation due to Brownian Motion—Part II: In the Continuum Regime

He Qing and Xie Mingliang

State Key Laboratory of Coal Combustion, Huazhong University of Science and Technology, Wuhan 430074, China

Received 24 August 2013; Revised 13 November 2013; Accepted 3 December 2013

Academic Editor: Shuyu Sun

Copyright © 2013 He Qing and Xie Mingliang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The fundamental aspects of the Taylor-series expansion method of moment (TEMOM) model proposed to model the aerosol population balance equation due to Brownian coagulation in the continuum regime is shown in this study, such as the choice of the expansion point u, the relationship between asymptotic behavior and analytical solution, and the error of the high-order moment equations. All these analyses will contribute to the buildup of the theoretical system of the TEMOM model.

1. Introduction

The population balance equations (PBE) are used to describe the evolution process of aerosol particles in a wide range of physical, chemical, and environmental subjects, such as nucleation, coagulation, diffusion, convection, and so on. When the Brownian coagulation plays a dominant role in such cases where aerosol particles at a high concentration are concerned or where suspended particles have evolved for a long time [1], the PBE for a monovariants system can be written as [2]

in which is the number density concentrations of the particles with volume from to at time , is the collision frequency function between particles with volumes and .

Because of its strong nonlinear partial integral-differential structure, the direct solution is very complicated and only a limited number of analytical solutions exist for simple coagulation kernel [3]. So several methods are proposed to solve this equation numerically, such as the sectional method (SM) [4], the Monte Carlo method (MCM) [5], and the method of moment (MM) [6]. With lower computational cost compared to the SM and MCM, the moment method has been widely used and become a powerful tool for investigating evolution processes of aerosol particles [7, 8].

By multiplying both the sides of (1) then integrating over the entire particle size distribution (PSD) [9], the growth rate of the particle moment can be obtained as follows:

where the moment is defined as

One main difficulty of the moment method is the closure of the moment equations. There exist several methods to overcome this bottleneck, including but not limited to the quadrature method of moment (QMOM) [10], the direct quadrature method of moment (DQMOM) [11], and the Taylor-series expansion method of moment (TEMOM) [3, 12].

It should be pointed out that the TEMOM has no prior assumption for the PSD using the Taylor-series expansion to achieve the closure and is considered as a promising approach to approximate the PBE for its relative simplicity of implementation and high accuracy [13]. Based on TEMOM model, the important information about the PSD, namely, the particle number density, particle mass, and geometric standard deviation, can be obtained for Brownian coagulation over the entire size regimes, and its results have a great agreement with other moment methods [3, 12–16]. But in these works, some fundamental problems are not clarified, for example, why the expansion point is set to be instead of other formulas; why the Taylor-series are truncated just at the first three terms; and what about the errors estimation of the present TEMOM model. In the present study, mainly as a methodological introduction, we would like to demonstrate the theoretical analysis to answer these questions for Brownian coagulation in the continuum regime.

2. Brief Review of TEMOM Model and Its Solutions in the Continuum Regime

At the initial time, the particle size maybe small in the free molecule regime. As time advances, the particle volume will grow due to coagulation between particles, and the particle size will transform to the near continuum regime via the transition regime and finally will tend to the continuum regime [17]. Therefore, the characteristic of particle evolution in the continuum regime is important to understand the coagulation mechanism. The collision frequency function β in the continuum regime is

where the constant , is the Boltzmann’s constant, is temperature, and is gas viscosity. Substituting (4) into (2), a set of moment equations including integral and fractional moments can be obtained. Using a Taylor-series expansion at for , the fractional moments can be approximated by the combination of integral moments as follows:

then the closure of the moment equations will be achieved automatically without any prior assumption about the particle size spectrum. The minimum number of moments for closing the equations is the first three-order moments , , and , which represents or are proportional to the total particle number concentration, and the total particle mass concentration, the total scattering light, respectively. According to the results derived by Yu et al. [3] and Xie and Wang [18], the equations can be rearranged as follows:

where the dimensionless moment is defined as

It is clear that remains constant due to the rigorous mass conservation requirement, and its initial conditions can be noted as , , and . Equations in (6) are a set of ordinary differential equations and can be solved directly. The main process is described briefly as follows. Because of the same structures of the first and the third equations in (6), the following relationship can be obtained:

This equation can be dissolved directly as follows:

where is the integration constant, and the dimensionless moment can be expressed as

Then substituting (10) into the first formula in (6), the relationship between and can be obtained:

in which is the integration constant. Then the second moment and dimensionless particle moment can be calculated by (9) and (10), respectively. As time advances, tends to zero due to coagulation; (11) can be simplified with some limit operation as

which is consistent with the asymptotic analysis shown by Xie et al. [18, 19].

3. The Choice of Expansion Point

In the TEMOM model, the choice of the expansion point at is not arbitrary. Some researchers think that the expansion for the characteristic size should take account of dispersion in the size spectrum and that is best done by using the well-known log-normal based expression in a second order closure formalism [6]. Using the expansion point, the corresponding moment equations based on TEMOM in the continuum regime are rewritten as

The detailed derivation is provided in the appendix. The comparison of numerical results between (6) and (13) is shown in Figure 1. And the initial conditions are selected as lognormal distribution: ; ; . The results show that the relative errors are small, and the two sets of equations are equivalent approximately. However, this selection results in more complicated moment equations. In fact, for any expression, when it is operated using the Taylor expansion technique, the selection of the expansion point is not unique. In mathematics, we only need to make sure in the targeted range than the Taylor-series expansion is convergent, and the final constructed moment equations are simple in the form. From this viewpoint, the selection of cannot be considered to be prior to the selection of .

(a)

(b)

(c)

(d)

Figure 1: The comparison of numerical results among (6), (13), (15), and (16) with the initial condition , , and ; (a) the numerical results of ; (b) the relative errors for ; (c) the numerical results of ; (d) the relative errors for .

4. The High-Order Moment Equations

The accuracy of the TEMOM model largely depends on the truncation errors of Taylor-series expansion. One method to determine the truncated errors is comparing the results of different TEMOM models, for example, the first three-moment model, the first four-moment model, the first five-moment model, and so forth. Similar to the derivation of the first three moment equations, is expanded at and truncated at the first four terms as

then the closed first four-moment equations can be obtained as follows:

with the same process, the closed first five-moment equations are

The comparison of numerical results among (5), (14), and (15) are also shown in Figure 1. The results show that the relative errors are small, and the three sets of equations can be considered equivalent nearly. Since the complexity of moment equations increases rapidly with the increasing number of reserved items but with little differences in the numerical results, we usually prefer using the first three-moment model proposed originally by Yu et al. [3].

5. Conclusion

Without a prior assumption for the shape of particle size distribution, the TEMOM has been considered as a promising method to model the aerosol population balance equation. In this study, the fundamental problems of the TEMOM model in the continuum regime due to Brownian coagulation are clarified, such as the choice of the expansion point and the error of the high-order moment equations. It benefits the understanding of PBE and TEMOM model.

Appendix

With the same process as Yu et al. [3] shown, we can get the first three-moment equations as follows:

Then using (3) and (5), the fractional moments can be approximated by the combination of the first three integral moments, for example, the is

Now replacing all the fractional moments by the integral moments (i.e., , , and ) in (A.2) and substituting and into (A.1), (13) is obtained.

Acknowledgments

The authors should appreciate the useful discussion and suggestion from Dr. Yu Mingzhou at China Jiliang University. This work is supported by the National Natural Science Foundation of China with Grant nos. 50806023 and 50721005, the Fundamental Research Funds for the Central Universities (Project no. 2013TS078), and the Program of Introducing Talents of Discipline to Universities (“111” Project no. B06019), China.

References

S. H. Park, K. W. Lee, E. Otto, and H. Fissan, “The log-normal size distribution theory of Brownian aerosol coagulation for the entire particle size range. Part I: analytical solution using the harmonic mean coagulation kernel,” Journal of Aerosol Science, vol. 30, no. 1, pp. 3–16, 1999.

S. K. Friedlander, Smoke, Dust, and Haze: Fundamentals of Aerosol Dynamics, Oxford University Press, London, UK, 2nd edition, 2000.

M. Z. Yu, J. Z. Lin, and T. L. Chan, “A new moment method for solving the coagulation equation for particles in Brownian motion,” Aerosol Science and Technology, vol. 42, no. 9, pp. 705–713, 2008.

S. S. Talukdar and M. T. Swihart, “Aerosol dynamics modeling of silicon nanoparticle formation during silane pyrolysis: a comparison of three solution methods,” Journal of Aerosol Science, vol. 35, no. 7, pp. 889–908, 2004.

N. M. Morgan, C. G. Wells, M. J. Goodson, M. Kraft, and W. Wagner, “A new numerical approach for the simulation of the growth of inorganic nanoparticles,” Journal of Computational Physics, vol. 211, no. 2, pp. 638–658, 2006.

S. E. Pratsinis, “Simultaneous nucleation, condensation, and coagulation in aerosol reactors,” Journal of Colloid And Interface Science, vol. 124, no. 2, pp. 416–427, 1988.

E. Otto, F. Stratmann, H. Fissan, S. Vemury, and S. E. Pratsinis, “39 P 17 Brownian coagulation in the transition regime II: a comparison of two modelling approaches,” Journal of Aerosol Science, vol. 24, no. 1, pp. S535–S536, 1993.

M.-Z. Yu, J.-Z. Lin, and T.-L. Chan, “Effect of precursor loading on non-spherical TiO₂ nanoparticle synthesis in a diffusion flame reactor,” Chemical Engineering Science, vol. 63, no. 9, pp. 2317–2329, 2008.

R. R. Upadhyay and O. A. Ezekoye, “Evaluation of the 1-point quadrature approximation in QMOM for combined aerosol growth laws,” Journal of Aerosol Science, vol. 34, no. 12, pp. 1665–1683, 2003.

R. McGraw, “Description of aerosol dynamics by the quadrature method of moments,” Aerosol Science and Technology, vol. 27, no. 2, pp. 255–265, 1997.

R. O. Fox, Computational Models for Turbulent Reacting Flows, Cambridge University Press, Cambridge, UK, 2003.

M. Z. Yu and J. Z. Lin, “Taylor-expansion moment method for agglomerate coagulation due to Brownian motion in the entire size regime,” Journal of Aerosol Science, vol. 40, no. 6, pp. 549–562, 2009.

M. Z. Yu and J. Z. Lin, “Solution of the agglomerate Brownian coagulation using Taylor-expansion moment method,” Journal of Colloid and Interface Science, vol. 336, no. 1, pp. 142–149, 2009.

M. Z. Yu and J. Z. Lin, “Binary homogeneous nucleation and growth of water-sulfuric acid nanoparticles using a TEMOM model,” International Journal of Heat and Mass Transfer, vol. 53, no. 4, pp. 635–644, 2010.

M. Z. Yu, J. Z. Lin, H. H. Jin, and Y. Jiang, “The verification of the Taylor-expansion moment method for the nanoparticle coagulation in the entire size regime due to Brownian motion,” Journal of Nanoparticle Research, vol. 13, no. 5, pp. 2007–2020, 2011.

M. L. Xie, M. Z. Yu, and L. P. Wang, “A TEMOM model to simulate nanoparticle growth in the temporal mixing layer due to Brownian coagulation,” Journal of Aerosol Science, vol. 54, pp. 32–48, 2012.

E. Otto, H. Fissan, S. H. Park, and K. W. Lee, “The log-normal size distribution theory of Brownian aerosol coagulation for the entire particle size range. Part II: analytical solution using Dahneke's coagulation kernel,” Journal of Aerosol Science, vol. 30, no. 1, pp. 17–34, 1999.

M. L. Xie and L. P. Wang, “Asymptotic solution of population balance equation based on TEMOM model,” Chemical Engineering Science, vol. 94, pp. 79–83, 2013.

M. L. Xie and Q. He, “Analytical solution of TEMOM model for particle coagulation due to Brownian motion,” Journal of Aerosol Science, vol. 66, pp. 24–30, 2013.