A New Moment Method for Solving the Coagulation Equation for Particles in Brownian Motion
|
|
- Garey Leonard
- 6 years ago
- Views:
Transcription
1 Aerosol Science and Technology ISSN: (Prin (Online Journal homepage: A New Moment Method for Solving the Coagulation Equation for Particles in Brownian Motion Mingzhou Yu, Jianzhong Lin & Tatleung Chan To cite this article: Mingzhou Yu, Jianzhong Lin & Tatleung Chan (8 A New Moment Method for Solving the Coagulation Equation for Particles in Brownian Motion, Aerosol Science and Technology, 4:9, 75-73, DOI:.8/ To link to this article: Published online: Aug 8. Submit your article to this journal Article views: View related articles Citing articles: 83 View citing articles Full Terms & Conditions of access and use can be found at
2 Aerosol Science and Technology, 4:75 73, 8 Copyright c American Association for Aerosol Research ISSN: print / online DOI:.8/ A New Moment Method for Solving the Coagulation Equation for Particles in Brownian Motion Mingzhou Yu,, Jianzhong Lin,, and Tatleung Chan 3 The State Key Laboratory of Fluid Power Transmission and Control, Zhejiang University, Hangzhou, China China Jiliang University, Hangzhou, China 3 Department of Mechanical Engineering, the Hong Kong Polytechnic University, Hung Hom, Kowloon, Hong Kong A new numerical approach for solving coagulation equation, model, is first presented. In this model, the closure of the moment equations is approached using the Taylor-series expansion technique. Through constructing a system of three first-order ordinary differential equations, the most important indexes for describing aerosol dynamics, including particle number density, particle mass and geometric standard deviation, are easily obtained. This approach has no prior requirement for particle size spectrum, and the limitation existing in the log-normal distribution theory automatically disappears. This new approach is tested by comparing it with known accurate solutions both in the free molecular and the continuum regime. The results show that this new approach can be used to solve the particle general dynamic equation undergoing Brownian coagulation with sufficient accuracy, while less computational cost is needed.. INTRODUCTION In theory, the evolution of aerosol size distribution in the flow field can be highly traced through solving the particle general dynamic equation (PGDE. This PGDE equation has an ability to describe the combined physical processes for very fine particles, such as the internal processes of nucleation, condensation, growth, and coagulation as well as the external transport processes of diffusion, convection, and thermophoresis. Unfortunately, the general dynamic equation is none other than Boltzmann s transport equation, in which only a limited number of known analytical solutions exist due to its own non-linear integro-differential structure. Hence, an alternative method, the numerical technique, has to be used to obtain approximate solution for this problem. However, the numerical calculations often become impractical, even with a modern super-computer, due to Received November 7; accepted 7 May 8. We would like to acknowledge the financial supported by the Major Program of the National Natural Science Foundation of China (Grant No Address correspondence to Jianzhong Lin, China Jiliang University, Hangzhou, 38, China. mecjzlin@public.zju.edu.cn the requirement of large computational cost. In order to break the limitation in computational cost, three prominent methods were used, i.e., the moment method (MM (Pratsinis 988; Hulbert and Katz 994, Chan et al. 6, Lin et al. 7, the sectional method (SM (Gelbard and Seinfeld 98; Talukdar and Swihart 4, and the stochastic particle method (SPM (Wells and Kraft 5; Morgan et al. 6. These methods have both advantages and disadvantages in accuracy and efficiency, and now they are used in different fields in terms of the particular requirement. Because of the relative simplicity of implementation and low computational cost, the MM has been extensively used to solve most particulate problems, and has become a powerful tool for investigating aerosol microphysical processes in most cases (Wright et al. ; Yu et al. 8b. As a key issue for the MM, the closure of the moment evolution equations has been achieved by making a prior assumption for the shape of the aerosol size distribution (Pratsinis 988, approximating the integral moment by an n-point Gaussian quadrature (McGraw 997 and assuming the pth-order polynomial form for the moments (Barrett and Jheeta 996. Typically, the moment method proposed by Pratsinis ( 988 and the quadrature method of moments (QMOM proposed by McGraw (997 are two widely used moment methods now. The is the early form of the MM, and has been developed to obtain the representative properties of the aerosol size distribution such as total particle concentration, mean size, and polydispersity only through solving the first three moments. However, the introduction of prior assumption with log-normal distribution in the will lead to the method invalid at early stage of particle formation (Yu et al. 8a, 8b. The QMOM, to some extent, is an improved form of the and has been proved to be applicable to any form of the growth laws and coagulation kernels (Upadhyay and Ezekoye 3. When the population balance is written in terms of one internal coordinate (e.g., particle length or size, the closure problem can be successfully solved with the use of a quadrature approximation without prior assumption. However, 75
3 76 M. YU ET AL. the weights and abscissas of the quadrature approximation must be additionally obtained using the product-difference algorithm (Gordon 968, which greatly increases the computational cost in contrast to the, especially in the case where the PGDE is required to be incorporated into the computational fluid dynamics (CFD codes. In recent years, the particle formation occurring in a chemically reacting flow has been increasingly received attention, and accordingly it is possible to simultaneously capture the details of the fluid flow and transport, the evolution of the particle size distribution and complex chemical kinetics (Yu et al. 8a, 8b. For these particulate systems, consuming large computational cost will be inevitable for further investigations. Although some researchers such as Friedlander ( and Frenklach ( have shown the moment method without requiring a priori knowledge of the size spectrum is not new, it is still necessary to construct a new approach with respect to moment equation, which is easy to implement with low computational cost like the and has not the prior requirement for particle size spectrum like the QMOM, to adapt to the requirement of modern complicated particulate industries. In this article we present a new approach to solve the PGDE undergoing Brownian coagulation. The underlying idea of the approach is that the closure of the moment evolution equations is approached using the Taylor-series expansion technique. The Taylor-series expansion method of moments ( has no prior requirement for particle size spectrum, and the number of moment equations needed is equal to the order of the Taylorseries expansion. Theoretically, the precision of solution and the computational cost increase with increasing the terms of Taylorseries expansion. Based on the trial numerical simulation, keeping 3 terms of Taylor-series is preferable if both precision and efficiency are simultaneously considered. The present investigation is concerned with the Brownian coagulation both in the free molecule and the continuum regime. Additionally, it should be pointed out that the is easy to be extended to investigate other aerosol dynamics such as coalescence, surface growth and nucleation, and the derivation of Bivariate extension of the simultaneously tracking the particle size and surface is also easy to be obtained, which is similar to the QMOM method (Wright et al... THEORY The scientific description for ultrafine particle dynamics starts from Schmoluchowski s discrete coagulation equation (Smoluchowski 97 and follows the integro-differential equation governing the continuous size distribution for the number concentration function developed by Muller (98. It is believed that the Muller s equation remains the absolute predominated status for investigating aerosol Brownian coagulation problems since its appearance in 98, which is written as: n(v, t = v β(v, v v n(v, n(v v, dv n(v, β(v, vn(v, dv, [] where n(v,dv is the number of particles whose volume is between v and v+dv at time t, and β(v, v isthe collision kernel for two particles of volumes v and v.inprinciple, depending on the functional form of the coagulation kernel and on the method of solution of the Smoluchowski equation, there can be various analytical solutions for Equation (. The general disposition for this problem is to transform Equation ( into an ordinary differential Equation with respect to the moment m k. The moment transformation involves multiplying Equation ( by v k and then integrating over the entire size distribution, and finally the transformed moment Equation based on the size distribution are obtained (Upadhyay and Ezekoye 3: m k t = [ (v + v k v k v k ] β(v, v n(v, n(v, dvdv, (k =,,, [] where the moment m k is defined by: m k = v k n(vdv. [3] In the past, some efforts have been made to achieve the closure of Equation (. Three prominent methods, i.e., making a prior assumption for the shape of the aerosol size distribution (Pratsinis 988; Yu et al. 6, approximating the integral moment by an n-point Gaussian quadrature (McGraw 997; Yu et al. 7, and assuming the pth-order polynomial form for the moments (Barrett and Jheeta 996, have been evaluated and compared (Barrett and Webb 998. In this article, we develop a new approach, Taylor expansion technique, to solve this problem... Brownian Coagulation in the Free Molecule Regime In this size regime the collision frequency function β FM is (Pratsinis 988: β FM = B (/v + /v /( v /3 + v /3, [4] where B = (3/4π /6 (6k b T/ρ /, k b is the Boltzmann constant, T is the gas temperature and ρ is the mass density of the particles. Let us first concentrate on (/v + /v / in Equation (4 and then rewrite it in the form: C = (/v + /v /. [5] In principle, as Equation (4 is introduced into Equation (, at the moment Equation ( cannot yet be closed because of the presence of the quadratic form of the variable C. Usually, Equation ( is brought to an integrable form using the approximation (Pratsinis and Kim 989; Lee and Chen 984: C = b ( /v / + /v /, [6]
4 ANEW MOMENT METHOD FOR BROWNIAN COAGULATION 77 where b is a strong function of the spread of the aerosol size spectrum. In this study, we expand Equation (5 in Taylor-series, then substitute it into Equation (, which is subsequently solved. In general, it is more convenient to expand (v + v / in binary Taylor-series than (/v + /v /. So, we are more willing to rewrite Equation (4 in the following particular form: β FM = B (v + v /( v /6 v / + v /6 v /6 + v / v /6. [7] Here, there is only necessary to expand (v + v / for Equation (7 in a Taylor-series about point (v = u, v = u. Using the l Hôpital s rule, it can be verified that the Taylor-series expansion converges in the interval (, u. In this study, we define the expansion point u as the mean particle size v. Atthe case, the convergence of the Smoluchowski Equation is highly dependent on a number of factors, including the size range of the particles, and on the time differencing scheme employed in the numerical simulation. In this study, we use fully implicit time differencing scheme to solve this stiff system. Under the condition of Taylorseries expansion, (v + v / should be: (v + v / = (v u u + 4 (v u + u 4 u (v u (v u(v u 3u 3/ 6u 3/ (v u + [8] 3u 3/ For Equation (8, we take the first three terms of Taylor-series expansion, which is consistent with the following treatments for the fractional moment. Disposing Equation (7 with Equation (8 and substituting it into Equation ( results in the following complicated equation with respect to the first three moments m, m, and m : where dm = B [ξ φ + ξ φ + ξ 3 φ 3 ]n(v, n(v, dvdv [9.] dm = [9.] dm = B [ς φ + ς φ + ς 3 φ 3 ]n(v, n(v, dvdv [9.3] ξ = v/6 v + v /6 v /6 ξ = v7/6 v + v5/6 v /6 + v/6, v + v5/6 v7/6 + /6 v + vv /6 v + v v /6, ς 3 = 4v 5/6 v 5/6 + v 7/6 3/ v + v v /6 + v 3/ v /6 + v v 7/6 + v3/6 + v/6 v v /6 + v3/6 + v/6 v v, /6 ς = 4v 5/6 v 5/6 + v 7/6 v + v 7/6 v, ς = v v 3/6 + 4v 5/6 v /6 + 4v 5/6 v /6 + v 3/6 7/6 v + v v 3/ + v 7/6 v 3/, ς 3 = 8v /6 v /6 + 4v 3/6 v 3/ + v 5/ v 7/6 + v 5/ v 7/6 + 4v 3/6 v 3/ + v v 9/6 + 4v 5/6 v 7/6, + vv 9/6 + 4v 5/6 v 7/6 φ = 3 u 8,φ = 3 8 u and φ 3 = 3u 3/. Using the transformed Equation (3, the right terms of Equation (9 can be integrated. At the case, their detailed expressions take the following form: where dm = B [ξ φ + ξ φ + ξ3 φ 3] [.] dm = [.] dm = B [ς φ + ς φ + ς3 φ 3] [.3] ξ = m /6m / + m /6 m /6 + m /6 m /, ξ = m 7/6m / + m 5/6 m /6 + m 5/6 m /6 + m 7/6 m / + m / m /6 + m / m /6, ξ = 4m 5/6m 5/6 + m 7/6 m / + m 3/ m /6 + m 3/ m /6 + m / m 7/6 + m 3/6 m /, +m /6 m /6 + m 3/6 m / + m /6 m /6, ς = 4m 5/6m 5/6 + m 7/6 m / + m 7/6 m /, ς = m /6m 3/6 + 4m 5/6 m /6 + 4m 5/6 m /6 + m 3/6 m / + m 7/6 m 3/ + m 7/6 m 3/, ς3 = 8m /6m /6 + 4m 3/6 m 3/ + m 5/ m 7/6 + m 5/ m 7/6 + 4m 3/6 m 3/ + m / m 9/6 + 4m 5/6 m 7/6 + m / m 9/6 + 4m 5/6 m 7/6 For all the coefficients from ξ to ς 3,itisclear that all the same terms are not grouped and all the terms completely correspond to the coefficients from ξ to ς 3 in Equation (9. In particular, Equation (9 denoted by these coefficients is absolutely composed with fractional-order moments, not integral moments. In order to further close the moment equation, it is necessary to dispose these fractional moments. Therefore, we continue to use the Taylor-series expansion technique to dispose Equation (3 with respect to fractional kth moment. In the Equation (3, v k can be expanded with Taylor-series about point v = u (this is
5 78 M. YU ET AL. consistent with the expansion point u in (v + v / : v k = u k +u k k(v u+ uk k(k (v u +. [] Similar to the disposition for (v + v /,wetake the first three terms of Taylor-series. In order to meet the requirement of moment transformation, it is necessary to group the terms of truncated Equation (: ( u v k k k = + uk k uk k v + ( u k k + u k kv + u k 3uk k. [] Substituting Equation ( into Equation (3, we have ( u m k = v k k k n(vdv = uk k + u k km + (u k + uk k 3uk k m + ( u k k m. [3] It is obvious that m k in Equation (3 can be considered as a function of the first three moments m, m, and m.inprinciple, the number of moments in Equation (3 is consistent with the reserved terms of Taylor-series in Equation (. For example, there is a need to take the first four terms of Taylor-series if m k is required to be a function of four moments m, m, m, and m 3.As the Equation (3 is substituted into Equation (, the moment equation merely comprised by integral moment variables can be easily obtained: dm = ( B 7m m /6 4m m m5/6 6859m 4 m /6 59m /6 [5.3] It is obvious that Equation (5 is a system of first-order ordinary differential equations and all the right terms are denoted by the first three moments m, m, and m, and thus this system can be automatically closed. Under these conditions, the first three moments, which are also the three predominate parameters for describing aerosol dynamics, are obtained through solving this first-order ordinary differential system. Here, it should be pointed out that the whole derivation of the equations does not involve any assumptions for particle size spectrum, while the final mathematical form is much simpler than the model... Brownian Coagulation in the Continuum Regime The disposition for aerosol Brownian coagulation equation in the free molecule regime can be also expanded to be in the continuum regime. In this regime the collision frequency function is (Barrett and Webb 998: ( β C = B v + (v /3 /3 v /3 + v /3, [6] where B = k bt, µ is the gas viscosity. Similar to the solution for Brownian coagulation in the free molecule regime, we µ substitute Equation (6 into Equation ( and have: dm dm dm = =, = B ( 4m m u 4388m u + 44m m u + 669m u4 65m + 69m m u 3 584u 3/6, [4.] [4.] B ( 6748m u + 7m 34m m u 76m m u 3 376m m u + 65m u4 59u /6 [4.3] Here, the expansion point of Taylor-series expansion is denoted by the mean particle size v(= m /m. At the case, Equation (4 can be written in the following form: dm = dm ( B 65m m 3/6 m m m7/6 93m 4 m/6 584m 3/6, [5.] =, [5.] dm dm dm = B (v /3 v /3 + + v /3 v /3 n(v, n(v, dvdv [7.] =, [7.] = B + v /3 v 4/3 ( v 4/3 v /3 + 4vv n(v, n(v, dvdv, [7.3]
6 ANEW MOMENT METHOD FOR BROWNIAN COAGULATION 79 and so on. After performing the moment transformation with Equation (3 together with the 3rd-order Taylor-series expansion for fractional moment, we can obtain the final form for aerosol Brownian coagulation in the continuum regime: dm = B ( 5m 4 + m m 3m m m m 8m 4, [8.] dm =, [8.] dm = B ( m m 3m m m 5m 4 8m [8.3] 3. COMPUTATIONS The numerical computations are all performed on an Intel (R Pentium 4 CPU 3. GHz computer with memory GB. The 4th-order Runge-Kutta method with fixed time step is used to solve the system of ordinary differential equations, and the Jacobi method is applied to obtain the eigenvalues and the eigenvectors of Jacobi matrix, which occurs in the computation of the QMOM. In the computations, all the numerical simulations are based on dimensionless time, and the time step is supposed to be small enough,., in order to allow an accurate result. The detailed dimensionless time is defined by τ = B t in the free molecular regime and τ = B t in the continuum regime. All the programs are written using the C Programming language and are performed on Microsoft Visual C++ 6. compiler. 4. RESULTS AND DISCUSSION We now test our new approach by comparing it with known accurate models (Barrett and Jheeta 996; Barrett and Webb 998 both in the free molecular and the continuum regime. In order to make better comparison, we have written the codes based on the QMOM model, the model and the model. For QMOM, five different cases are included, i.e., quadrature approximation with, 3, 4, 5, and 6 nodes. The used here is based on the log-normal distribution approximation developed by Pratsinis (988. The initial distribution of all cases is assumed to be log-normal, i.e., n(v, = N exp{ (ln v ln v g /w g }/{ πvw g } with N =, v g = 3/ and wg = ln(4/3, which is consistent with the initial distribution conducted by Barrett and Jheeta (996, Barrett and Webb (998, and Williams (986. Under the case with log-normal initial distribution, the initial kth moments m k are obtained through solving m k = Nv k g ek w g /. The calculation of relative error for any variables follows the definition: %Error = φ φ reference φ reference %, [9] where φ is the arbitrary variable and φ reference is the referenced variable. 4.. Brownian Coagulation in the Free Molecular Regime In this regime the computation time is scaled by the definition τ = B t, which is consistent with the disposition by Barrett and Webb (998. Here, we follow the study of Barrett and Webb (998 in which they compared various models at several special time points using tabulation approach. Table shows the values of zero moment m and second moment m given by various methods (the first moment m is a constant when only Brownian coagulation is considered at time τ =., 5., and.. In this table, the values are all from Table 3ofBarrett and Webb (998 except for the and the. For QMOM,, and models, it should be noted that our own codes give the same values with those provided by Barrett and Webb (998. The last column in Table shows the consumed CPU time up to τ =. Since some studies have shown that the QMOM model with or 3 nodes produces the accurate result, and the precision can be increased with increasing the number of nodes (McGraw 997; Upadhyay and Ezekoye 3, it is feasible to take the as a reference to TABLE Moments m and m in the free molecular regime at different dimensionless time and the consumed dimensionless time for various methods τ = τ = 5 τ = CPU time (B Method m m m m m m (τ = Laguerre quad QMOM points QMOM 3 points QMOM 6 points (Pratsinis Gamma (Williams (presen
7 7 M. YU ET AL.. m QMOM Error% 5 QMOM τ(b (a τ(b (b FIG.. The zero moment m and the relative error for various methods in the free molecular regime. The relative error denotes the ratio of m from various methods to the QMOM with 6 points. validate other models. In fact, all the approximate methods in Table give reasonable agreement with the reference results. For the zero moment, the results obtained by are all identical to those by the except at τ =.. At τ =., the result of is closer to than Gamma model, QMOM, and models. Although different models provide different results for the second moment, the model produces the same precision as and model. In addition, it can be easily seen that the consumed CPU time for the is the shortest among all the investigated models. In order to further validate this new method, it is necessary to have the comparisons between our model and other known models over large evolution time. The evolutions of m with scaled time τ together with the relative error in various methods are shown in Figure. The relative error denotes the ratio of m from the various methods to the. In fact, all the errors mentioned in this paper, including the zero and second moment, are all the relative values of,, or QMOM- to. In Figure a, all the curves overlap each other, and it is difficult to distinguish from one to the other. Hence, the differences between them must be evaluated using the relative error, Error%, which is shown in Figure b. From this figure, it can be found that the error for the is always less than % from τ = toτ =. In addition, it is clear that the curves of error exhibit the same trend for all the models, and more importantly, the gives the least relative error beyond τ = 6. For all the investigated models, the relative error reaches constant values at τ>6, indicating that these models can give the reasonable results for zero moment. The second moment m as well as the relative error are shown in Figure. Like the zero moment, the relative error is still denoted by the ratio of various methods to the. Similar to Figure a, all the curves overlap in Figure a. In Figure b, the curves of relative error are very similar for the and, except that values are negative for the former but positive for the latter. For the, the error initially increases in the negative direction and then holds nearly m 3 QMOM Error% QMOM τ(b (a τ(b (b FIG.. The second moment m and the relative error for various methods in the free molecular regime. The relative error denotes the ratio of m from various methods to the QMOM with 6 points.
8 ANEW MOMENT METHOD FOR BROWNIAN COAGULATION 7 σ g QMOM τ(b FIG. 3. The geometric standard deviation σ g in the free molecular regime for various methods. constant at.65%. For the, the error starts to show its largest value at τ =., and then decreases until it attains its nearly constant value at about τ =. However, the absolute values of relative error of the are less than that of the in the entire time range. This indicates the can achieve the same precision as the for solving the second moment, although the precision is little less than that produced by the. There exists an asymptotic solution for the distribution of reduced particle size, the so-called self-preserving size distribution (SPSD (Frenklach 985, which has become an important tool to explore the aerosol coagulation mechanisms. The selfpreserving form is usually approximated by a lognormal distribution with a geometric standard deviation σ g which is obtained by solving the following equation (Pratsinis 988: ln σ g = 9 ln ( m m m. [] Although Equation ( was originally derivated from lognormal theory, it may be also applied in QMOM and models. Figure 3 shows the variations of σ g with time in the free molecular regime for various methods. We can see that the SPSD attains for all methods at about τ = 5. For the QMOM, the asymptotic values of σ g decrease from.378 to.346 with an increase of nodes from to 6. The asymptotic value of.345 for the is the closest to the among all the investigated models. In addition, the asymptotic value of.355 for the is identical to the value given by Lee et al. (984 who used log-normal functions for particle size distribution. The above comparisons suggest that the has a higher precision in describing the asymptotic aerosol size distribution than other methods when compared to. 4.. Brownian Coagulation in the Continuum Regime Similar to the treatment for Brownian coagulation in the free molecular regime, we now compare various methods in the continuum regime. The results from the are considered to be exact and are used as references. The computational time in Equation (8 is scaled by τ = B t. The variations of the zero moment m and the relative error for various methods in the continuum regime are shown in Figure 4. In Figure 4a, all the curves overlap with each other and thus we cannot distinguish one from another. Hence, it is still necessary to use the relative error to distinguish the precisions produced by all the investigated models. In the whole range up to τ =, Figure 4b shows the model nearly exhibits the same curve with the model, and in particular, the results of are closer to the. Moreover, like the and QMOM, the relative error produced by converges to a nearly constant value beyond τ =, and more importantly, this error is always below.5%. The results clearly indicate that the is capable of handling the evolution of zero moment with high precision in the continuum regime. The variations of second moment and relative error for various methods are shown in Figure 5. Similarly, all the curves. 4 m QMOM τ(b (a Error% τ(b (b QMOM FIG. 4. The zero moment m and the relative error for various methods in the continuum regime. The relative error denotes the ratio of m from various methods to the QMOM with 6 points.
9 7 M. YU ET AL. 3. m QMOM τ(b (a Error% τ(b (b QMOM FIG. 5. The second moment m and the relative error for various methods in the continuum regime. The relative error denotes the ratio of m from various methods to the QMOM with 6 points. exhibit the same trend and there are no visible differences with each other in Figure 5a. In the entire evolution time up to τ =, it is clear that in all the investigated models, the QMOM model produces the largest error, but its value is always below.6%. These suggest all the investigated models can achieve very high precision for representing second moment in this regime. Although the error produced by the is slightly larger than the, their curves almost show the same trend. Like the zero moment, the error produced by the converges to a nearly constant value of.4% beyond τ =. Typically, the curves of other models show the same trend except that there are negative values below τ = 5 for the. From the above comparisons, for the second moment, it can be easily found that the gives relative low error close to the. In the continuum regime, the particle size distribution can also achieve the self-preserving size spectrum. All the studies (Lee et al. 997; Pratsinis 988; Lee 983; Lee et al. 984 showed that the asymptotic value based on the geometric standard deviation σ g is near to.3 if the log-normal distribution theory is used. In Figure 6, the values of the geometric standard σ g QMOM τ(b FIG. 6. The geometric standard deviation σ g in the continuum regime for various methods. deviation for various methods are presented and compared. It is obvious that all curves show the same trend and attain each asymptotic value at about τ =. However, like in the free molecular regime, the diversity for these asymptotic values still exists in this regime. For the QMOM, the asymptotic value of σ g decreases from.35 to.35 with an increase of nodes from to 6. In this figure, the curves for are not shown because their asymptotic values are always between QMOM and. Here, it is worth pointing out that the gives the same curve as the model, and their asymptotic values are the same value of.39, which is consistent with the asymptotic value of.3 reported in the literatures (Lee et al. 997; Lee 983; Lee et al This suggests that the has an ability to handle the evolution of particle size spectrum with time based on geometric standard deviation, in particular with the same precision as the model. 5. CONCLUSIONS In this study we used the Taylor-series expansion technique to dispose the collision terms and the fractional moments to obtain a new form for the moment equations. This approach requires no prior requirements for particle size spectrum, and the number of moment equations to be solved simultaneously is equal to the reserved terms of Taylor-series. The present study was focused on Brownian coagulation in the free molecular regime as well as in the continuum regime, and the new approach was validated by comparing it with known methods. In order to validate model in accuracy and efficiency, three important indexes of aerosol dynamics, i.e., the zero moment, the second moment and the geometric standard deviation were compared. The results show that the new approach can be used to solve particle general dynamic equation undergoing Brownian coagulation with sufficient accuracy, while less computational cost is needed.
10 ANEW MOMENT METHOD FOR BROWNIAN COAGULATION 73 REFERENCES Barrett, J. C., and Jheeta, J. S. (996. Improving the Accuracy of the Moments Method for Solving the Aerosol General Dynamic Equation, J. Aerosol Sci. 7:35 4. Barrett, J. C., and Webb, N. A. (998. A Comparison of Some Approximate Methods for Solving the Aerosol General Dynamic Equation, J. Aerosol Sci. 9:3 39. Chan, T. L., Lin, J. Z., Zhou, K., and Chan, C. K. (6. Simultaneous Numerical Simulation of Nano and Fine Particle Coagulation and Dispersion in a Round Jet, J. Aerosol Sci. 37: Frenklach, M. (985. Dynamics of Discrete Distribution for Smoluchowski Coagulation Model, J. Coll. Interface Sci. 8:37 4. Frenklach, M. (. Method of Moments with Interpolative Closure, Chem. Engineer. Sci. 57:9 39. Frenklach, M., and Harris, S. J. (987. Aerosol Dynamics Modeling Using the Method of Moments, J. Coll. Interface Sci. 8:5 6. Friedlander, S. K. (. Smoke, Dust, and Haze. Fundamentals of Aerosol Dynamics. Oxford University Press. Gelbard, F., and Seinfeld, J. H. (98. Simulation of Multicomponent Aerosol Dynamics, J. Coll. Interface Sci. 78: Gordon, R. G. (968. Error Bounds in Equilibrium Statistical Mechanics, Journal Mathematical Physics 9: Hulbert, H. M., and Katz, S. (964. Some Problems in Particle Technology: A Statistical Mechanical Formulation, Chem. Engineer. Sci. 9: Lee, K. W. (983. Change of Particle Size Distribution During Brownian Coagulation, J. Coll. Interface Sci. 9: Lee, K. W., Chen, H., and Gieseke, J. A. (984. Log-Normally Preserving Size Distribution for Brownian Coagulation in the Free-Molecule Regime, Aerosol Sci. Technol. 3:53 6. Lee, K. W., and Chen, H. (984. Coagulation Rate of Polydisperse Particles, Aerosol Sci. Technol. 3: Lee, K. W., Lee, Y. J., and Han, D. S. (997. The Log-Normal Size Distribution Theory for Brownian Coagulation in the Low Knudsen Number Regime, J. Coll. Interface Sci. 88: Lin, J. Z., Chan, T. L., Lin, S., Zhou, K., Zhou, Y., and Lee, S. C. (7. Effects of Coherent Structures on Nanoparticle Coagulation and Dispersion in a Round Jet, Intl. J. Nonlinear Sci. Numer. Simul. 8: McGraw, R. (997. Description of Aerosol Dynamics by the Quadrature Method of Moments, Aerosol Sci. Technol. 7: Morgan, N. M., Wells, C. G., Goodson, M. J., Kraft, M., and Wagner, W. (6. A New Numerical Approach for the Simulation of the Growth of Inorganic Nanoparticles, J. Computational Physics : Muller, H. (98. Zur Allgemeinen Theorie Der Raschen Koagulation, Kolloideihefte 7:3 5. Pratsinis, S. E. (988. Simultaneous Nucleation, Condensation, and Coagulation in Aerosol Reactor, J. Coll. Interface Sci. 4: Pratsinis, S. E., and Kim, K. S. (989. Particle Coagulation, Diffusion and Thermophoresis in Laminar Tube Flows, J. Aerosol Sci. :. Settumba, N., and Garrick, S. C. (3. Direct Numerical Simulation of Nanoparticle Coagulation in a Temporal Mixing Layer via a Moment Method, J. Aerosol Sci. 34: Smoluchowski, V. (97. Versuch Einer Mathematischen Theorie Der Koagulationskinetik Kollider Losungen, Zeitschrift für Physikalische Chemie 9:9 68. Talukdar, S. S., and Swihart, M. T. (4. Aerosol Dynamics Modeling of Silicon Nanoparticle Formation During Silane Pyrolysis: A Comparison of Three Solution Methods, J. Aerosol Sci. 35: Upadhyay, R. R., and Ezekoye, O. A. (3. Evaluation of the -Point Quadrature Approximation in QMOM for Combined Aerosol Growth Laws, J. Aerosol Sci. 34: Wells, C. G., and Kraft, M. (5. Direct Simulation and Mass Flow Stochastic Algorithms to Solve a Sintering-Coagulation Equation, Monte Carlo Methods and Applications : Williams, M. M. R. (986. Some Topics in Nuclear Aerosol Dynamics, Progress in Nuclear Energy 7: 5. Wright, D. L., McGraw, R. and Rosner, D. E. (. Bivariate Extension of the Quadrature Method of Moments for Modeling Simultaneous Coagulation and Sintering of Particle Populations, J. Coll. Interface Sci. 36:4 5. Yu,M.Z., Lin, J. Z., and Chan, T. L. (6. Large Eddy Simulation of a Planar Jet Flow with Nanoparticle Coagulation, Acta Mechanica Sinica :93 3. Yu,M.Z., Lin, J. Z., and Chen L. H. (7. Nanoparticle Coagulation in a Planar Jet Via Moment Method, Appl. Math. Mech. 8: Yu, M. Z., Lin, J. Z., and Chan, T. L. (8a. Numerical Simulation of Nanoparticle Synthesis in Diffusion Flame Reactor, Powder Technology 8:9. Yu,M.Z., Lin, J. Z., and Chan, T. L. (8b. Effect of Precursor Loading on Non Spherical TiO Nanoparticle Synthesis in a Diffusion Flame Reactor, Chem. Engineer. Sci. 63:37 39.
THE VERIFICATION OF THE TAYLOR-EXPANSION MOMENT METHOD IN SOLVING AEROSOL BREAKAGE
144 THERMAL SCIENCE, Year 1, Vol. 16, No. 5, pp. 144-148 THE VERIFICATION OF THE TAYLOR-EXPANSION MOMENT METHOD IN SOLVING AEROSOL BREAKAGE by Ming-Zhou YU a,b * and Kai ZHANG a a College of Science, China
More informationNANOPARTICLE COAGULATION AND DISPERSION IN A TURBULENT PLANAR JET WITH CONSTRAINTS
THERMAL SCIENCE, Year 2012, Vol. 16, No. 5, pp. 1497-1501 1497 NANOPARTICLE COAGULATION AND DISPERSION IN A TURBULENT PLANAR JET WITH CONSTRAINTS by Cheng-Xu TU a,b * and Song LIU a a Department of Mechanics,
More informationNUMERICAL MODELING OF FINE PARTICLE FRACTAL AGGREGATES IN TURBULENT FLOW
THERMAL SCIENCE, Year 2015, Vol. 19, No. 4, pp. 1189-1193 1189 NUMERICAL MODELING OF FINE PARTICLE FRACTAL AGGREGATES IN TURBULENT FLOW by Feifeng CAO a, Zhanhong WAN b,c*, Minmin WANG b, Zhenjiang YOU
More informationFor one such collision, a new particle, "k", is formed of volume v k = v i + v j. The rate of formation of "k" particles is, N ij
Coagulation Theory and The Smoluchowski Equation: Coagulation is defined as growth of particles by collisions among particles. It is usually associated with dense, 3-d growth, in contrast to aggregation
More informationINSTANTANEOUS AEROSOL DYNAMICS IN A TURBULENT FLOW
1492 THERMAL SCIENCE, Year 2012, Vol. 16, No. 5, pp. 1492-1496 INSTANTANEOUS AEROSOL DYNAMICS IN A TURBULENT FLOW by Kun ZHOU * Clean Combustion Research Center, King Abdullah University of Science and
More informationJournal of Aerosol Science
Journal of Aerosol Science 5 (1) 3 8 Contents lists available at SciVerse ScienceDirect Journal of Aerosol Science journal homepage: www.elsevier.com/locate/jaerosci A TEMOM model to simulate nanoparticle
More informationMODELING ON THE BREAKUP OF VISCO-ELASTIC LIQUID FOR EFFERVESCENT ATOMIZATION
1446 THERMAL SCIENCE, Year 2012, Vol. 16, No. 5, pp. 1446-1450 MODELING ON THE BREAKUP OF VISCO-ELASTIC LIQUID FOR EFFERVESCENT ATOMIZATION by Li-Juan QIAN * China Jiliang University, Hangzhou, China Short
More informationMonte Carlo Simulation of Multicomponent Aerosols Undergoing Simultaneous Coagulation and Condensation
Aerosol Science and Technology ISSN: 0278-6826 (Print) 1521-7388 (Online) Journal homepage: http://www.tandfonline.com/loi/uast20 Monte Carlo Simulation of Multicomponent Aerosols Undergoing Simultaneous
More informationA new differentially weighted operator splitting Monte Carlo method for aerosol dynamics
This paper is part of the Proceedings of the 24 International Conference th on Modelling, Monitoring and Management of Air Pollution (AIR 2016) www.witconferences.com A new differentially weighted operator
More informationA Simple Numerical Algorithm and Software for Solution of Nucleation, Surface Growth, and Coagulation Problems
Aerosol Science and Technology, 37:892 898, 2003 Copyright c American Association for Aerosol Research ISSN: 0278-6826 print / 1521-7388 online DOI: 10.1080/02786820390225826 A Simple Numerical Algorithm
More informationOPTIMAL WAVELENGTH SELECTION ALGORITHM OF NON-SPHERICAL PARTICLE SIZE DISTRIBUTION BASED ON THE LIGHT EXTINCTION DATA
THERMAL SCIENCE, Year 2012, Vol. 16, No. 5, pp. 1353-1357 1353 OPTIMAL WAVELENGTH SELECTION ALGORITHM OF NON-SPHERICAL PARTICLE SIZE ISTRIBUTION BASE ON THE LIGHT EXTINCTION ATA by Hong TANG * College
More informationParticle Dynamics: Brownian Diffusion
Particle Dynamics: Brownian Diffusion Prof. Sotiris E. Pratsinis Particle Technology Laboratory Department of Mechanical and Process Engineering, ETH Zürich, Switzerland www.ptl.ethz.ch 1 or or or or nucleation
More informationflame synthesis of inorganic nanoparticles
flame synthesis of inorganic nanoparticles Shraddha Shekar and Jethro Akroyd Markus KRAFT Karlsruhe 12 July 2010 Part 1 Introduction Precursor (TES) Mesoporous silica nanoparticles Aim: To answer the following
More informationOPTIMAL WAVELENGTH SELECTION ALGORITHM OF NON-SPHERICAL PARTICLE SIZE DISTRIBUTION BASED ON THE LIGHT EXTINCTION DATA
THERMAL SCIENCE, Year 202, Vol. 6, No. 5, pp. 353-357 353 OPTIMAL WAVELENGTH SELECTION ALGORITHM OF NON-SPHERICAL PARTICLE SIZE ISTRIBUTION BASE ON THE LIGHT EXTINCTION ATA y Hong TANG * College of Metrology
More informationA Fast Monte Carlo Method Based on an Acceptance-Rejection Scheme for Particle Coagulation
Aerosol and Air Quality Research, 3: 73 8, 03 Copyright Taiwan Association for Aerosol Research ISSN: 680-8584 print / 07-409 online doi: 0.409/aaqr.0..0369 A Fast Monte Carlo Method Based on an Acceptance-Rejection
More informationCHARACTERISTIC DISTRIBUTION OF SUBMICRON AND NANO-PARTICLES LADEN FLOW AROUND CIRCULAR CYLINDER
1386 THERMAL SCIENCE, Year 2012, Vol. 16, No. 5, pp. 1386-1390 CHARACTERISTIC DISTRIBUTION OF SUBMICRON AND NANO-PARTICLES LADEN FLOW AROUND CIRCULAR CYLINDER by Cheng-Xu TU a,b* and Jian ZHANG b a Department
More informationDYNAMIC STABILITY OF NON-DILUTE FIBER SHEAR SUSPENSIONS
THERMAL SCIENCE, Year 2012, Vol. 16, No. 5, pp. 1551-1555 1551 DYNAMIC STABILITY OF NON-DILUTE FIBER SHEAR SUSPENSIONS by Zhan-Hong WAN a*, Zhen-Jiang YOU b, and Chang-Bin WANG c a Department of Ocean
More informationModels for Nucleation and Condensation on Nanoparticles
PFAU X Models for Nucleation and Condensation on Nanoparticles Stefan Radl, 1 Mario Schriefl 2,3 1 TU Graz, Institute of Process and Particle Engineering 2 TU Graz, Institute of Electronics Tristan Reinisch,
More informationKinetic Monte Carlo simulation of the effect of coalescence energy release on the size and shape evolution of nanoparticles grown as an aerosol
JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 6 8 AUGUST 2003 Kinetic Monte Carlo simulation of the effect of coalescence energy release on the size and shape evolution of nanoparticles grown as an aerosol
More informationModeling and simulation of multi-component aerosol dynamics
Modeling and simulation of multi-component aerosol dynamics Y. Efendiev March 1, 23 Abstract In the paper we consider modeling of heterogeneous aerosol coagulation where the heterogeneous aerosol particles
More informationAggregate Growth: R =αn 1/ d f
Aggregate Growth: Mass-ractal aggregates are partly described by the mass-ractal dimension, d, that deines the relationship between size and mass, R =αn 1/ d where α is the lacunarity constant, R is the
More informationImplementation of a discrete nodal model to probe the effect of size-dependent surface tension on nanoparticle formation and growth
Aerosol Science 37 (2006) 1388 1399 www.elsevier.com/locate/jaerosci Implementation of a discrete nodal model to probe the effect of size-dependent surface tension on nanoparticle formation and growth
More informationEvaluation of diffusion models for airborne nanoparticles transport and dispersion
Advances in Fluid Mechanics VII 111 Evaluation of diffusion models for airborne nanoparticles transport and dispersion F. Morency 1, S. Hallé 1, L. Dufresne 1 & C. Émond 2 1 École de technologie supérieure,
More informationPLEASE SCROLL DOWN FOR ARTICLE. Full terms and conditions of use:
This article was downloaded by: On: 3 June 2010 Access details: Access Details: Free Access Publisher Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
More informationarxiv: v1 [physics.comp-ph] 22 Jul 2010
Gaussian integration with rescaling of abscissas and weights arxiv:007.38v [physics.comp-ph] 22 Jul 200 A. Odrzywolek M. Smoluchowski Institute of Physics, Jagiellonian University, Cracov, Poland Abstract
More informationABSTRACT INTRODUCTION
MODELING OF ALUMINUM NANOPARTICLE FORMATION R. Schefflan D. Kalyon S. Kovenklioglu Stevens Institute of Technology Highly Filled Materials Institute Castle Point Station Hoboken, N.J. 73 ABSTRACT Picatinny
More informationA Unified Gas-kinetic Scheme for Continuum and Rarefied Flows
A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows K. Xu and J.C. Huang Mathematics Department, Hong Kong University of Science and Technology, Hong Kong Department of Merchant Marine, National
More informationNanoparticle coagulation via a Navier Stokes/nodal methodology: Evolution of the particle field
Aerosol Science 37 (2006) 555 576 www.elsevier.com/locate/jaerosci Nanoparticle coagulation via a Navier Stokes/nodal methodology: Evolution of the particle field S.C. Garrick a,, K.E.J. Lehtinen b, M.R.
More informationPrediction of Minimum Fluidisation Velocity Using a CFD-PBM Coupled Model in an Industrial Gas Phase Polymerisation Reactor
Journal of Engineering Science, Vol. 10, 95 105, 2014 Prediction of Minimum Fluidisation Velocity Using a CFD-PBM Coupled Model in an Industrial Gas Phase Polymerisation Reactor Vahid Akbari and Mohd.
More informationResearch of Micro-Rectangular-Channel Flow Based on Lattice Boltzmann Method
Research Journal of Applied Sciences, Engineering and Technology 6(14): 50-55, 013 ISSN: 040-7459; e-issn: 040-7467 Maxwell Scientific Organization, 013 Submitted: November 08, 01 Accepted: December 8,
More informationAnalysis of Particle Contamination in Plasma Reactor by 2-Sized Particle Growth Model
Korean J. Chem. Eng., 20(2), 392-398 (2003) Analysis of Particle Contamination in Plasma Reactor by 2-Sized Particle Growth Model Dong-Joo Kim, Pil Jo Lyoo* and Kyo-Seon Kim Department of Chemical Engineering,
More informationStudy of Numerical Diffusion in a Discrete- Sectional Model and Its Application to Aerosol Dynamics Simulation
Aerosol Science and Technology ISSN: 0278-6826 (Print) 1521-7388 (Online) Journal homepage: http://www.tandfonline.com/loi/uast20 Study of Numerical Diffusion in a Discrete- Sectional Model and Its Application
More informationAerosol nucleation and growth in a turbulent jet using the Stochastic Fields method
Loughborough University Institutional Repository Aerosol nucleation and growth in a turbulent jet using the Stochastic Fields method This item was submitted to Loughborough University's Institutional Repository
More informationOn the detailed modelling of high temperature nanoparticles synthesis
On the detailed modelling of high temperature nanoparticles synthesis Prof. http://como.cheng.cam.ac.uk/ The 9 th International Conference on Chemical Kinetics (ICCK) Department of Chemical Engineering
More informationEVOLUTION OF PARTICLE SIZE DISTRIBUTION IN AIR IN THE RAINFALL PROCESS VIA THE MOMENT METHOD
137 THERMAL SCIENCE, Year 1, Vol. 16, No. 5, pp. 137-1376 EVOLUTION OF PARTICLE SIZE DISTRIBUTION IN AIR IN THE RAINFALL PROCESS VIA THE MOMENT METHOD by Fu-Jun GAN a an Jian-Zhong LIN a,b * a Department
More informationSummary We develop an unconditionally stable explicit particle CFD scheme: Boltzmann Particle Hydrodynamics (BPH)
Summary Steps in BPH Space is divided into Cartesian cells Finite number of particles ~10 6-10 7 Particles fly freely between t n and t n+1 Mass, momentum and total energy are conserved Relax into a (LTE)
More informationPhysical Modeling of Multiphase flow. Boltzmann method
with lattice Boltzmann method Exa Corp., Burlington, MA, USA Feburary, 2011 Scope Scope Re-examine the non-ideal gas model in [Shan & Chen, Phys. Rev. E, (1993)] from the perspective of kinetic theory
More informationLow Variance Particle Simulations of the Boltzmann Transport Equation for the Variable Hard Sphere Collision Model
Low Variance Particle Simulations of the Boltzmann Transport Equation for the Variable Hard Sphere Collision Model G. A. Radtke, N. G. Hadjiconstantinou and W. Wagner Massachusetts Institute of Technology,
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationMultiphase Flow Simulations in Inclined Tubes with Lattice Boltzmann Method on GPU
Multiphase Flow Simulations in Inclined Tubes with Lattice Boltzmann Method on GPU Khramtsov D.P., Nekrasov D.A., Pokusaev B.G. Department of Thermodynamics, Thermal Engineering and Energy Saving Technologies,
More informationMonte Carlo simulation for simultaneous particle coagulation and deposition
222 Science in China: Series E Technological Sciences 2006 Vol.49 No.2 222 237 DOI: 10.1007/s11431-006-0222-3 Monte Carlo simulation for simultaneous particle coagulation and deposition ZHAO Haibo & ZHENG
More informationResearch Article Parallel Monte Carlo Simulation of Aerosol Dynamics
Advances in Mechanical Engineering, Article ID 435936, 11 pages http://dx.doi.org/10.1155/2014/435936 Research Article Parallel Monte Carlo Simulation of Aerosol Dynamics Kun Zhou, 1,2 Zhu He, 1 Ming Xiao,
More informationNUMERICAL PREDICTIONS OF DEPOSTION WITH A PARTICLE CLOUD TRACKING TECHNIQUE
Committed Individuals Solving Challenging Problems NUMERICAL PREDICTIONS OF DEPOSTION WITH A PARTICLE CLOUD TRACKING TECHNIQUE by James R. Valentine Reaction Engineering International Philip J. Smith Department
More informationLattice Bhatnagar Gross Krook model for the Lorenz attractor
Physica D 154 (2001) 43 50 Lattice Bhatnagar Gross Krook model for the Lorenz attractor Guangwu Yan a,b,,liyuan a a LSEC, Institute of Computational Mathematics, Academy of Mathematics and System Sciences,
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations Simo Särkkä Aalto University, Finland November 18, 2014 Simo Särkkä (Aalto) Lecture 4: Numerical Solution of SDEs November
More informationauthor's personal copy
Mechanics Research Communications 40 (2012) 46 51 Contents lists available at SciVerse ScienceDirect Mechanics Research Communications jo ur nal homep age : www.elsevier.com/locate/mechrescom Verifying
More informationAggregation and Breakage of Nanoparticle Dispersions in Heterogeneous Turbulent Flows
ggregation and reakage of Nanoparticle Dispersions in Heterogeneous Turbulent Flows M. Soos,, D. Marchisio*, J. Sefcik, M. Morbidelli Swiss Federal Institute of Technology Zürich, Z CH *Politecnico di
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationThe Brownian Oscillator
Chapter The Brownian Oscillator Contents.One-DimensionalDiffusioninaHarmonicPotential...38 The one-dimensional Smoluchowski equation, in case that a stationary flux-free equilibrium state p o (x) exists,
More informationTime-Dependent Aerosol Models and Homogeneous Nucleation Rates
Aerosol Science and Technology ISSN: 0278-6826 (Print) 1521-7388 (Online) Journal homepage: https:www.tandfonline.comloiuast20 Time-Dependent Aerosol Models and Homogeneous Nucleation Rates S. L. Girshick,
More informationPopulation balance modeling -an application in particle technology
Population balance modeling -an application in particle technology Sheryl Ehrman IITB and University of Maryland Link to Matlab: www.glue.umd.edu/~sehrman/popbal.htm Outline Aerosol reactors in industry
More informationSYNTHESIS OF INORGANIC NANOPARTICLES IN SUPERCRITICAL FLUIDS. A SIMPLE MODEL TO PREDICT PARTICLE SIZE
SYNTHESIS OF INORGANIC NANOPARTICLES IN SUPERCRITICAL FLUIDS. A SIMPLE MODEL TO PREDICT PARTICLE SIZE A. Erriguible 1, F. Marias 2, F. Cansell 3, C. Aymonier* 3 1 Laboratoire TREFLE (UMR C.N.R.S. 8508),
More informationComparison of maximum entropy and quadrature-based moment closures for shock transitions prediction in one-dimensional gaskinetic theory
Comparison of maximum entropy and quadrature-based moment closures for shock transitions prediction in one-dimensional gaskinetic theory Jérémie Laplante and Clinton P. T. Groth Citation: AIP Conference
More informationCrossover from ballistic to Epstein diffusion in the freemolecular
This is the author s final, peer-reviewed manuscript as accepted for publication. The publisher-formatted version may be available through the publisher s web site or your institution s library. Crossover
More informationKinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles
Kinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3. The 8th International Conference on Computational
More informationPredictive multiscale computational modeling of nanoparticles in flame reactors
Graduate Theses and Dissertations Graduate College 2012 Predictive multiscale computational modeling of nanoparticles in flame reactors Maulik Mehta Iowa State University Follow this and additional works
More informationSimplified Hyperbolic Moment Equations
Simplified Hyperbolic Moment Equations Julian Koellermeier and Manuel Torrilhon Abstract Hyperbolicity is a necessary property of model equations for the solution of the BGK equation to achieve stable
More informationCost-accuracy analysis of a variational nodal 2D/1D approach to pin resolved neutron transport
Cost-accuracy analysis of a variational nodal 2D/1D approach to pin resolved neutron transport ZHANG Tengfei 1, WU Hongchun 1, CAO Liangzhi 1, LEWIS Elmer-E. 2, SMITH Micheal-A. 3, and YANG Won-sik 4 1.
More informationMehrphasensimulationen zur Charakterisierung einer Aerosol-Teilprobenahme
Mehrphasensimulationen zur Charakterisierung einer Aerosol-Teilprobenahme G. Lindner 8.41, Modellierung und Simulation February 29, 2016 Review CFD Study for the Design of the Dilution Section The mixing
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel
More informationChapter 9 Generation of (Nano)Particles by Growth
Chapter 9 Generation of (Nano)Particles by Growth 9.1 Nucleation (1) Supersaturation Thermodynamics assumes a phase change takes place when there reaches Saturation of vapor in a gas, Saturation of solute
More informationOn the Boltzmann equation: global solutions in one spatial dimension
On the Boltzmann equation: global solutions in one spatial dimension Department of Mathematics & Statistics Colloque de mathématiques de Montréal Centre de Recherches Mathématiques November 11, 2005 Collaborators
More informationX i t react. ~min i max i. R ij smallest. X j. Physical processes by characteristic timescale. largest. t diff ~ L2 D. t sound. ~ L a. t flow.
Physical processes by characteristic timescale Diffusive timescale t diff ~ L2 D largest Sound crossing timescale t sound ~ L a Flow timescale t flow ~ L u Free fall timescale Cooling timescale Reaction
More informationA stencil of the finite-difference method for the 2D convection diffusion equation and its new iterative scheme
International Journal of Computer Mathematics Vol. 87, No. 11, September 2010, 2588 2600 A stencil of the finite-difference method for the 2D convection diffusion equation and its new iterative scheme
More informationThe Measure of Diffusion Skewness and Kurtosis in Magnetic Resonance Imaging
The Measure of Diffusion Skewness and Kurtosis in Magnetic Resonance Imaging Xinzhen Zhang, Chen Ling, Liqun Qi, and Ed Xuekui Wu April 19, 008, Revised on September 4, 008 This paper is dedicated to the
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More informationDiffusion and Adsorption in porous media. Ali Ahmadpour Chemical Eng. Dept. Ferdowsi University of Mashhad
Diffusion and Adsorption in porous media Ali Ahmadpour Chemical Eng. Dept. Ferdowsi University of Mashhad Contents Introduction Devices used to Measure Diffusion in Porous Solids Modes of transport in
More informationDirect observation of aggregative nanoparticle growth: Kinetic modeling of the size distribution and growth rate
Direct observation of aggregative nanoparticle growth: Kinetic modeling of the size distribution and growth rate Taylor J. Woehl, * Chiwoo Park, James E. Evans, Ilke Arslan, William D. Ristenpart, and
More informationModelling and numerical methods for the diffusion of impurities in a gas
INERNAIONAL JOURNAL FOR NUMERICAL MEHODS IN FLUIDS Int. J. Numer. Meth. Fluids 6; : 6 [Version: /9/8 v.] Modelling and numerical methods for the diffusion of impurities in a gas E. Ferrari, L. Pareschi
More informationStatistical Mechanics for the Truncated Quasi-Geostrophic Equations
Statistical Mechanics for the Truncated Quasi-Geostrophic Equations Di Qi, and Andrew J. Majda Courant Institute of Mathematical Sciences Fall 6 Advanced Topics in Applied Math Di Qi, and Andrew J. Majda
More informationPrinciples of Convective Heat Transfer
Massoud Kaviany Principles of Convective Heat Transfer Second Edition With 378 Figures Springer Contents Series Preface Preface to the Second Edition Preface to the First Edition Acknowledgments vii ix
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations Simo Särkkä Aalto University Tampere University of Technology Lappeenranta University of Technology Finland November
More informationarxiv: v1 [physics.ao-ph] 23 Feb 2015
Approximate analytical solutions to the condensation-coagulation equation of aerosols arxiv:5.6399v [physics.ao-ph] 3 Feb 5 Naftali Smith, Nir J. Shaviv,, Henrik Svensmark 3 Racah Institute of Physics,
More informationExponential moments for the homogeneous Kac equation
Exponential moments for the homogeneous Kac equation Maja Tasković University of Pennsylvania Young Researchers Workshop: Stochastic and deterministic methods in kinetic theory, Duke University, November
More informationPolitecnico di Torino. Porto Institutional Repository
Politecnico di Torino Porto Institutional Repository [Proceeding] Simulation of macromolecule self-assembly in solution: multiscale approach a Original Citation: Lavino, A; Di Pasquale, N; Carbone, P;
More informationA Taylor series solution of the reactor point kinetics equations
1 A Taylor series solution of the reactor point kinetics equations David McMahon and Adam Pierson Department of Nuclear Safety Analysis,Sandia National Laboratories, Albuquerque, NM 87185-1141 E-mail address:
More informationOverview of Turbulent Reacting Flows
Overview of Turbulent Reacting Flows Outline Various Applications Overview of available reacting flow models LES Latest additions Example Cases Summary Reacting Flows Applications in STAR-CCM+ Ever-Expanding
More informationLecture 3: From Random Walks to Continuum Diffusion
Lecture 3: From Random Walks to Continuum Diffusion Martin Z. Bazant Department of Mathematics, MIT February 3, 6 Overview In the previous lecture (by Prof. Yip), we discussed how individual particles
More informationDirect Modeling for Computational Fluid Dynamics
Direct Modeling for Computational Fluid Dynamics Kun Xu February 20, 2013 Computational fluid dynamics (CFD) is new emerging scientific discipline, and targets to simulate fluid motion in different scales.
More information19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007
19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, -7 SEPTEMBER 007 INVESTIGATION OF AMPLITUDE DEPENDENCE ON NONLINEAR ACOUSTICS USING THE DIRECT SIMULATION MONTE CARLO METHOD PACS: 43.5.Ed Hanford, Amanda
More informationStatistical mechanics of random billiard systems
Statistical mechanics of random billiard systems Renato Feres Washington University, St. Louis Banff, August 2014 1 / 39 Acknowledgements Collaborators: Timothy Chumley, U. of Iowa Scott Cook, Swarthmore
More informationMultiscale Modeling of TiO2 Nanoparticle Production in Flame Reactors: Effect of Chemical Mechanism
Chemical and Biological Engineering Publications Chemical and Biological Engineering 2010 Multiscale Modeling of TiO2 Nanoparticle Production in Flame Reactors: Effect of Chemical Mechanism Maulik Mehta
More informationFast Probability Generating Function Method for Stochastic Chemical Reaction Networks
MATCH Communications in Mathematical and in Computer Chemistry MATCH Commun. Math. Comput. Chem. 71 (2014) 57-69 ISSN 0340-6253 Fast Probability Generating Function Method for Stochastic Chemical Reaction
More informationA Study On The Heat Transfer of Nanofluids in Pipes
Project Report 2014 MVK160 Heat and Mass Transport May 15, 2014, Lund, Sweden A Study On The Heat Transfer of Nanofluids in Pipes Koh Kai Liang Peter Dept. of Energy Sciences, Faculty of Engineering, Lund
More informationA Gauss Lobatto quadrature method for solving optimal control problems
ANZIAM J. 47 (EMAC2005) pp.c101 C115, 2006 C101 A Gauss Lobatto quadrature method for solving optimal control problems P. Williams (Received 29 August 2005; revised 13 July 2006) Abstract This paper proposes
More informationA Sectional Model for Investigating Microcontamination in a Rotating Disk CVD Reactor
Aerosol Science and Technology, 38:1161 1170, 2004 Copyright c American Association for Aerosol Research ISSN: 0278-6826 print / 1521-7388 online DOI: 10.1080/027868290896799 A Sectional Model for Investigating
More informationStefan Stefanov Bulgarian Academy of Science, Bulgaria Ali Amiri-Jaghargh Ehsan Roohi Hamid Niazmand Ferdowsi University of Mashhad, Iran
Stefan Stefanov Bulgarian Academy of Science, Bulgaria Ali Amiri-Jaghargh Ehsan Roohi Hamid Niazmand Ferdowsi University of Mashhad, Iran Outlines: Introduction DSMC Collision Schemes Results Conclusion
More informationAssignment 3: The Atomic Nature of Matter
Introduction Assignment 3: The Atomic Nature of Matter Science One CS 2015-2016 Instructor: Mike Gelbart Last Updated Sunday, Jan 24, 2016 around 7:45pm Part One due Saturday, Jan 30, 2016 at 5:00pm Part
More informationNumerical Methods. King Saud University
Numerical Methods King Saud University Aims In this lecture, we will... find the approximate solutions of derivative (first- and second-order) and antiderivative (definite integral only). Numerical Differentiation
More informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationNUMERICAL MODELING OF THE GAS-PARTICLE FLUID FLOW AND HEAT TRANSFER IN THE SLIP REGIME
Proceedings of the Asian Conference on Thermal Sciences 2017, 1st ACTS March 26-30, 2017, Jeju Island, Korea ACTS-P00394 NUMERICAL MODELING OF THE GAS-PARTICLE FLUID FLOW AND HEAT TRANSFER IN THE SLIP
More informationTHREE-DIMENSIONAL HAUSDORFF DERIVATIVE DIFFUSION MODEL FOR ISOTROPIC/ANISOTROPIC FRACTAL POROUS MEDIA
Cai, W., et al.: Three-Dimensional Hausdorff Derivative Diffusion Model... THERMAL SCIENCE: Year 08, Vol., Suppl., pp. S-S6 S THREE-DIMENSIONAL HAUSDORFF DERIVATIVE DIFFUSION MODEL FOR ISOTROPIC/ANISOTROPIC
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationMonte Carlo simulations of dense gas flow and heat transfer in micro- and nano-channels
Science in China Ser. E Engineering & Materials Science 2005 Vol.48 No.3 317 325 317 Monte Carlo simulations of dense gas flow and heat transfer in micro- and nano-channels WANG Moran & LI Zhixin Department
More informationInvestigation of soot morphology and particle size distribution in a turbulent nonpremixed flame via Monte Carlo simulations
Investigation of soot morphology and particle size distribution in a turbulent nonpremixed flame via Monte Carlo simulations Ahmed Abdelgadir* a, Marco Lucchesi b, Antonio Attili a, Fabrizio Bisetti a
More informationTime-Dependent Statistical Mechanics 1. Introduction
Time-Dependent Statistical Mechanics 1. Introduction c Hans C. Andersen Announcements September 24, 2009 Lecture 1 9/22/09 1 Topics of concern in the course We shall be concerned with the time dependent
More informationComputational Aspects of Data Assimilation for Aerosol Dynamics
Computational Aspects of Data Assimilation for Aerosol Dynamics A. Sandu 1, W. Liao 1, G.R. Carmichael 2, D. Henze 3, J.H. Seinfeld 3, T. Chai 2, and D. Daescu 4 1 Department of Computer Science Virginia
More informationResearch Article A Novel Differential Evolution Invasive Weed Optimization Algorithm for Solving Nonlinear Equations Systems
Journal of Applied Mathematics Volume 2013, Article ID 757391, 18 pages http://dx.doi.org/10.1155/2013/757391 Research Article A Novel Differential Evolution Invasive Weed Optimization for Solving Nonlinear
More informationA Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional Hyperbolic Telegraph Equation
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.13(01) No.3,pp.59-66 A Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional
More information