研究论文

有限系统中两组分颗粒凝并问题的异权值蒙特卡罗模拟

展开
  • 上海大学 上海市应用数学和力学研究所,上海 200072

收稿日期: 2017-04-11

  网络出版日期: 2018-12-21

基金资助

国家自然科学基金资助项目(11332006);国家自然科学基金资助项目(11272196)

Multi-Monte Carlo simulation for bicomponent coalescence in finite system

Expand
  • Shanghai Institute of Applied Mathematics & Mechanics, Shanghai University, Shanghai 200072, China

Received date: 2017-04-11

  Online published: 2018-12-21

摘要

对于两组分颗粒系统,某组分 A 在整个颗粒群中的混合均匀程度是一个非常重要的物理量。采用异权值蒙特卡罗方法研究有限系统中凝并核与组分相关的两组分颗粒凝并混合程度问题。首先,验证了核函数齐次性指数 lambda 对系统凝并混合的影响,并得到了有限系统中分散指数与 lambda 的变化规律。其次,发现并总结了有限系统中考虑不同组分间吸引或排斥作用时颗粒凝并混合规律:临界时间之前系统在吸引排斥因素以及核函数的共同作用下发生颗粒凝并,且排斥作用的影响远远大于吸引作用的影响;临界时间之后系统在核函数的单独作用下进行颗粒凝并,排斥或吸引作用的影响可以忽略不计。最后,拟合得到了临界时间与吸引排斥指数 alpha 的幂函数关系,以及临界时间与齐次性指数 $\lambda$ 的指数函数关系。研究结果对药物混合工业过程的设计具有一定的理论指导意义。

本文引用格式

左昊, 沈杰, 卢志明 . 有限系统中两组分颗粒凝并问题的异权值蒙特卡罗模拟[J]. 上海大学学报(自然科学版), 2020 , 26(4) : 617 -627 . DOI: 10.12066/j.issn.1007-2861.2055

Abstract

The mixing state of a bicomponent population of granules is characterized by the total variance of component A, which measures the deviation of the composition of each granule from the overall mean. By means of the multi-weighted Monte Carlo method, first the fundamental influence of the homogeneity index lambda of the kernel is verified, and the scaling rate of segregation index, which is a function of homogeneity index lambda, is obtained. Second, two stages in the evolution of mixing are identified. Before the critical time, both interaction between components and the kernel take important effects on the coalescence. And the repulsive interaction is much stronger than the attractive interaction. After the critical time, only the kernel itself affects the coalescence and the effects of interaction between components are ignored. Finally, a power-law function between the critical time and alpha is obtained, whereas an exponential function relationship between the critical time and $\lambda$ is identified by a curve fitting technique. The results are of great importance in the pharmaceutical engineering.

参考文献

[1] 赵海波, 郑楚光. 离散系统动力学演变过程的颗粒群平衡模拟 [M]. 北京: 科学出版社, 2008: 1-30.
[2] Kuang C, McMurry P H, McCormick A V. Determination of cloud condensation nuclei production from measured new particle formation events[J]. Geophysical Research Letters, 2009,36(9). DOI: 10.1029/2009GL037584.
[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.
[4] Krapivsky P L, Ben-Naim E. Aggregation with multiple conservation laws[J]. Physical Review E, 1996,53(1):291-298.
[5] 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(D9). DOI: 10.1029/2008JD011073.
[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.
[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.
[8] 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(11):3715-3727.
[9] 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.
[10] Lushnikov A A. Evolution of coagulating systems: Ⅲ. coagulating mixtures[J]. Journal of Colloid and Interface Science, 1976,54(1):94-101.
[11] 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.
[12] Vigil R D, Ziff R M. On the scaling theory of two-component aggregation[J]. Chemical engineering science, 1998,53(9):1725-1729.
[13] 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]. Europhysics Letters, 2007,78(5):56002.
[14] 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.
[15] Piskunov V N. Analytical solutions for coagulation and condensation kinetics of composite particles[J]. Physica D: Nonlinear Phenomena, 2013,249:38-45.
[16] 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.
[17] Barrasso D, Walia S, Ramachandran R. Multi-component population balance modeling of continuous granulation processes: a parametric study and comparison with experimental trends[J]. Powder Technology, 2013,241:85-97.
[18] Hashemian N, Ghanaatpishe M, Armaou A. Development of a reduced order model for bi-component granulation processes via laguerre polynomials[C]// American Control Conference, 2016. DOI: 10.1109/ACC.2016.7525483.
[19] Burgalat J, Rannou P. Brownian coagulation of a bi-modal distribution of both spherical and fractal aerosols[J]. Journal of Aerosol Science, 2107,105:151-165.
[20] Matsoukas T, Lee K, Kim T. Mixing of components in two component aggregation[J]. AIChE journal, 2006,52(9):3088-3099.
[21] Lee K, Kim T, Rajniak P, et al. Compositional distributions in multicomponent aggre-gation[J]. Chemical Engineering Science, 2008,63(5):1293-1303.
[22] Zhao H B, Kruis F E, Zheng C. Monte Carlo simulation for aggregative mixing of nanoparticles in two-component systems[J]. Industrial & Engineering Chemistry Research, 2011,50(13):10652-10664.
[23] Zhao H B, Einar Kruis F. Dependence of steady-state compositional mixing degree on feeding conditions in two-component aggregation[J]. Industrial & Engineering Chemistry Research, 2014,53(14):6047-6055.
[24] Zhao H B, Zheng C G, Xu M H. Multi-Monte Carlo approach for general dynamic equation considering simultaneous particle coagulation and breakage[J]. Powder Technology, 2005,154(2):164-178.
[25] Zhao H B, Kruis F E, Zheng C G. 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.
[26] Zhao H Q, Xu Z W, Zhao H B. Predictions on dynamic evolution of compositional mixing degree in two-component aggregation[J]. Journal of Aerosol Science, 2016,101:10-21.
[27] 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.
[28] Matsoukas T, Jr Marshall C L. Bicomponent aggregation in finite systems[J]. Europhysics Letters, 2010,92(4):46007.
[29] Jiang Z, Shen J, Lu Z M. Bivariate Taylor-series expansion method of moment for particle population balance equation in Brownian coagulation[J]. Journal of Aerosol Science, 2017,114:94-106.
文章导航

/