两组分颗粒系统中凝并混合程度问题的多重蒙特卡罗模拟

doi:10.12066/j.issn.1007-2861.1972

摘要/Abstract

摘要：

多组分颗粒凝并问题有重要应用背景,如制药过程中有效成分分散、雨滴形成和晶体结晶过程等,因而对其深入研究有重要理论和实际意义.采用多重蒙特卡罗方法研究两组分颗粒系统中布朗凝并问题,其中凝并核与组分有关. 数值结果表明:两组分颗粒系统的矩统计量与单组分颗粒系统中的矩统计量有相同的渐近行为.当系统混合度达到稳定值并且分散指数以$1/\overline{v}$下降时,系统趋于自保持状态. 另外,还得到了不同初始条件下某组分总偏差的渐近稳定值与"凝并效率"参数的函数关系式.

关键词: 两组分, 凝并, 组分相关核, 混合度

Abstract:

Multi-component aggregation has important applications in areas such as dispersion of active ingredients in a pharmaceutical process, formation of raindroplets, and crystallization. In this paper, a population balance equation with composition dependent kernel is numerically simulated using multiple Monte Carlo methods. The results show that the moments of a two-component system have a scaling behavior similar to a single component system regardless of the kernel type. While the degree of mixing reaches a constant value and the segregation index decreases by $1/\overline {v}$, the system approaches a self-preserving state. A fitted function of the stable value of the power density of excess component A on collision efficiency is derived.

Key words: two-component, aggregation, composition dependent kernel, mixing degree

中图分类号:

O302
O359

姜志, 沈杰, 卢志明. 两组分颗粒系统中凝并混合程度问题的多重蒙特卡罗模拟[J]. 上海大学学报(自然科学版), 2019, 25(2): 245-255.

JIANG Zhi, SHEN Jie, LU Zhiming. Monte Carlo simulation for mixing degree in two-component system[J]. Journal of Shanghai University（Natural Science Edition）, 2019, 25(2): 245-255.

图/表 10

图1

图2

图3

图4

图5

图6

图7

图8

参考文献 28

[1]	赵海波, 郑楚光 . 离散系统动力学演变过程的颗粒群平衡模拟[M]. 北京: 科学出版社, 2008: 1-30.
[2]	Kuang C ,Mc Murry P H,Mc Cormick A V. Determination of cloud condensation nuclei production from measured new particle formation events[J]. Geophysical Research Letters, 2009. DOI: 10.1029/2009GL037584. doi: 10.1029/2019GL084173 pmid: 31762521
[3]	Senior C L, Zeng T, Che J , et al. Distribution of trace elements in selected pulverized coals as a function of particle size and density[J]. Fuel Processing Technology, 2000,63(2):215-241. doi: 10.1016/S0378-3820(99)00098-3
[4]	Riemer N, West M, Zaveri R A , et al. Simulating the evolution of soot mixing state with a particle-resolved aerosol model[J]. Journal of Geophysical Research: Atmospheres, 2009,114:D09202.
[5]	Liu Y, Jiang Y, Zhang X , et al. Structural and magnetic properties of the ordered FePt$_3$, FePt and Fe$_3$Pt nanoparticles[J]. Journal of Solid State Chemistry, 2014,209:69-73. doi: 10.1016/j.jssc.2013.10.027
[6]	Singh J ,Patterson R I A,Kraft M, et al. Numerical simulation and sensitivity analysis of detailed soot particle size distribution in laminar premixed ethylene flames[J]. Combustion and Flame, 2006,145(1):117-127. doi: 10.1016/j.combustflame.2005.11.003
[7]	Rajniak P, Mancinelli C, Chern R T , et al. Experimental study of wet granulation in fluidized bed: impact of the binder properties on the granule morphology[J]. International Journal of Pharmaceutics, 2007,334(1):92-102. doi: 10.1016/j.ijpharm.2006.10.040 pmid: 17207948
[8]	Marshall C L, Rajniak P, Matsoukas T . Numerical simulations of two-component granulation: comparison of three methods[J]. Chemical Engineering Research and Design, 2011,89(5):545-552. doi: 10.1016/j.cherd.2010.06.003
[9]	Puel F ,Févotte G,Klein J P. Simulation and analysis of industrial crystallization processes through multidimensional population balance equations. Part 1: a resolution algorithm based on the method of classes[J]. Chemical Engineering Science, 2003,58(16):3715-3727. doi: 10.1016/S0009-2509(03)00254-9
[10]	Krapivsky P L, Ben-Naim E . Aggregation with multiple conservation laws[J]. Physical Review E, 1996,53(1):291-298. doi: 10.1103/physreve.53.291 pmid: 9964259
[11]	Lushnikov A A . Evolution of coagulating systems: Ⅲ. Coagulating mixtures[J]. Journal of Colloid and Interface Science, 1976,54(1):94-101. doi: 10.1016/0021-9797(76)90288-5
[12]	Gelbard F M, Seinfeld J H . Coagulation and growth of a multicomponent aerosol[J]. Journal of Colloid and Interface Science, 1978,63(3):472-479. doi: 10.1016/S0021-9797(78)80008-3
[13]	Vigil R D, Ziff R M . On the scaling theory of two-component aggregation[J]. Chemical Engineering Science, 1998,53(9):1725-1729. doi: 10.1016/S0009-2509(98)00016-5
[14]	FernÁndez-DÍaz J M, Gómez-GarcÍa G J . Exact solution of Smoluchowski's continuous multi-component equation with an additive kernel[J]. EPL (Europhysics Letters), 2007,78(5):56002. doi: 10.1209/0295-5075/78/56002
[15]	FernÁndez-Díaz J M, Gómez-García G J . Exact solution of a coagulation equation with a product kernel in the multicomponent case[J]. Physica D: Nonlinear Phenomena, 2010,239(5):279-290. doi: 10.1016/j.physd.2009.11.010
[16]	Piskunov V N . Analytical solutions for coagulation and condensation kinetics of composite particles[J]. Physica D: Nonlinear Phenomena, 2013,249:38-45. doi: 10.1016/j.physd.2013.01.008
[17]	Yu M, Lin J, Chan T . A new moment method for solving the coagulation equation for particles in Brownian motion[J]. Aerosol Science and Technology, 2008,42(9):705-713. doi: 10.1080/02786820802232972
[18]	Matsoukas T, Lee K, Kim T . Mixing of components in two-component aggregation[J]. AIChE Journal, 2006,52(9):3088-3099. doi: 10.1016/j.colsurfb.2017.03.013 pmid: 28347946
[19]	Lee K, Kim T, Rajniak P , et al. Compositional distributions in multicomponentaggregation[J]. Chemical Engineering Science, 2008,63(5):1293-1303. doi: 10.1007/s11356-019-06784-6 pmid: 31797274
[20]	Zhao H, Kruis F E, Zheng C . Monte Carlo simulation for aggregative mixing of nanoparticles in two-component systems[J]. Industrial & Engineering Chemistry Research, 2011,50(18):10652-10664. doi: 10.6133/apjcn.201912_28(4).0017 pmid: 31826378
[21]	Zhao H B, Kruis F E . Dependence of steady-state compositional mixing degree on feeding conditions in two-component aggregation[J]. Industrial & Engineering Chemistry Research, 2014,53(14):6047-6055. doi: 10.6133/apjcn.201912_28(4).0017 pmid: 31826378
[22]	Haibo Z, Chuguang Z, Minghou X . Multi-Monte Carlo approach for general dynamic equation considering simultaneous particle coagulation and breakage[J]. Powder Technology, 2005,154(2):164-178. doi: 10.1016/j.powtec.2005.04.042
[23]	Zhao H, Kruis F E, Zheng C . Reducing statistical noise and extending the size spectrum by applying weighted simulation particles in Monte Carlo simulation of coagulation[J]. Aerosol Science and Technology, 2009,43(8):781-793. doi: 10.1080/02786820902939708
[24]	Matsoukas T, Kim T, Lee K . Bicomponent aggregation with composition-dependent rates and the approach to well-mixed state[J]. Chemical Engineering Science, 2009,64(4):787-799. doi: 10.1016/j.ces.2008.04.060
[25]	Mller H . Zur allgemeinen theorie ser raschen koagulation[J]. Fortschrittsberichte ber Kolloide und Polymere, 1928,27(6):223-250.
[26]	Wright D L ,McGraw R,Rosner D E.Bivariate extension of the quadrature method of moments for modeling simultaneous coagulation and sintering of particle populations[J]. Journal of Colloid and Interface Science, 2001,236(2):242-251. doi: 10.1006/jcis.2000.7409 pmid: 11401370
[27]	Yu M, Zhang X, Jin G , et al. A new analytical solution for solving the population balance equation in the continuum-slip regime[J]. Journal of Aerosol Science, 2015,80:1-10. doi: 10.1016/j.jaerosci.2014.10.007
[28]	Xie M L, Wang L P . Asymptotic solution of population balance equation based on TEMOM model[J]. Chemical Engineering Science, 2013,94:79-83. doi: 10.1016/j.ces.2013.02.025