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

Authors: He Qing; Xie Mingliang
Journal: Abstract and Applied Analysis
 (2013)

---

## 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 pointu, 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.


---

## Body

## 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]
(1)∂n(υ,t)∂t=12∫0vβ(υ1,υ-υ1)n(υ1,t)n(υ-υ1,t)dυ1-∫0∞β(υ1,υ)n(υ,t)n(υ1,t)dυ1,
in which n(υ,t) is the number density concentrations of the particles with volume from υ to υ+dυ at time t, β(υ,υ1) is the collision frequency function between particles with volumes υ and υ1.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υk 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:
(2)dMkdt=12∬0∞[(υ+υ1)k-υk-υ1k]cc×β(υ,υ1)n(υ,t)n(υ1,t)dυdυ1,ccccccccccccccccc(k=0,1,2,…),
where the moment Mk is defined as
(3)Mk=∫0∞υkn(υ)dυ.
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 u is set to be M1/M0 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
(4)β=B2(1υ1/3+1υ11/3)(υ1/3+υ11/3),
where the constant B2=3kBT/2μ, kB is the Boltzmann’s constant, T 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 υ=u=M1/M0 for υk, the fractional moments can be approximated by the combination of integral moments as follows:
(5)υk=uk+uk-1k(υ-u)+uk-2k(k-1)2(υ-u)2+uk-3k(k-1)(k-2)3!(υ-u)3+⋯;
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 M0, M1, and M2, 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:
(6)dM0dt=B2(2MC2-13MC-151)M0281,dM1dt=0,dM2dt=-2B2(2MC2-13MC-151)M1281,
where the dimensionless moment MC is defined as
(7)MC=M0M2M12.
It is clear that M1 remains constant due to the rigorous mass conservation requirement, and its initial conditions can be noted as M00, M20, and MC0=M00M20/M12. 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:
(8)dM2dM0=-2M12M02.
This equation can be dissolved directly as follows:
(9)M2=2M12M0+C1M12,C1=M20M12-2M00,
where C1 is the integration constant, and the dimensionless moment MC can be expressed as
(10)MC=2+M0C1.
Then substituting (10) into the first formula in (6), the relationship between M0 and t can be obtained:
(11)5lnM0+169C1M0-52ln|169+5C1M0-2C12M02|-701917arctanh|4C1M0-5|917=134B234C1t+C2,C2=5lnM00+169C1M00-52ln|169+5C1M00-2C12M002|-701917arctanh|4C1M00-5|917,
in which C2 is the integration constant. Then the second moment M2 and dimensionless particle moment MC can be calculated by (9) and (10), respectively. As time advances, M0 tends to zero due to coagulation; (11) can be simplified with some limit operation as
(12)limM0→01M0=limt→∞169B281t,
which is consistent with the asymptotic analysis shown by Xie et al. [18, 19].


## 3. The Choice of Expansion Pointu

In the TEMOM model, the choice of the expansion point atu=M1/M0 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 u=M0-3/2M12M2-1/2 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
(13)dM0dt=-B2M02(-2MC4+17MC5/2-4MC2-35MC+35MC1/2+151)81,dM1dt=0,dM2dt=-2B2M12(2MC4-8MC5/2-5MC2-145MC-8MC1/2+2)81MC.
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: M00=1; M1=1; M20=4/3. 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 u=M0-3/2M12M2-1/2 cannot be considered to be prior to the selection of u=M1/M0.The comparison of numerical results among (6), (13), (15), and (16) with the initial condition M00=1, M1=1, and M20=4/3; (a) the numerical results of M0; (b) the relative errors for M0; (c) the numerical results of M2; (d) the relative errors for M2.
(a)
(b)
(c)
(d)


## 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,υk is expanded at υ=u and truncated at the first four terms as
(14)υk=uk+uk-1k(υ-u)+uk-2k(k-1)2(υ-u)2+uk-3k(k-1)(k-2)3!(υ-u)3+⋯;
then the closed first four-moment equations can be obtained as follows:
(15)dM0dt=B26561(70M06M32M16-636M05M3M2M15cccc+1440M04M22M14+1225M04M3M13cccc-5160M03M2M12-10061M02),dM1dt=0,dM2dt=B26561(-32M04M32M14+408M03M3M2M13ccccc-1260M02M22M12-344M02M3M1ccccc+3570M0M2+23902M12-32M04M32M14+408M03M3M2M13),dM3dt=B22187(76M03M32M13-348M02M3M2M12cccc+925M3M0+36M0M22M1cccc+25791M1M2-236M13M0);
with the same process, the closed first five-moment equations are
(16)dM0dt=B2M0259049M18×(-59049M18+8520M04M4M2M12-11550M03M4M14-45474M03M3M2M13+60060M02M3M15-128700M0M2M16-3745M05M4M3M1+51480M02M22M14+350M06M42+10010M04M32M12),dM1dt=0,dM2dt=-2B259049M16×(-525M03M4M14-8160M03M3M2M13+4080M02M3M15-22050M0M2M16-106299M18+12600M02M22M14+1280M04M32M12+35M06M42+1365M04M4M2M12-424M05M4M3M1),dM3dt=B219683M15×(-1105M04M4M3M1+915M03M4M2M12-930M02M4M14-9870M02M3M2M13+13584M0M3M15+226458M2M16+3150M0M22M14+3920M03M32M12+74M05M42),dM4dt=-2B259049M14M02×(695M06M42-2905M05M4M3M1+8340M04M4M2M12-16545M03M4M14+14430M03M3M2M13+22770M0M2M16-481344M02M3M15-360297M02M22M14-7020M18-4810M04M32M12).
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 pointu and the error of the high-order moment equations. It benefits the understanding of PBE and TEMOM model.



---

*Source: 490123-2013-12-19.xml*