求解考虑颗粒凝并的通用动力学方程的多重MonteCarlo算法

赵海波; 郑楚光; 徐明厚

求解考虑颗粒凝并的通用动力学方程的多重MonteCarlo算法

华中科技大学煤燃烧国家重点实验室, 武汉 430074

基金项目: 国家重点基础研究专项经费资助项目(2002CB211602);国家自然科学基金资助项目(重点)(90410017)

详细信息

作者简介:
赵海波(1977- ),男,湖南宁乡人,讲师,博士(联系人.Tel/Fax:+86-27-87545526;E-mail:klinsmannzhb@163.com).

中图分类号: O343.5
计量
- 文章访问数: 2546
- HTML全文浏览量: 109
- PDF下载量: 1371
- 被引次数: 0
出版历程
- 收稿日期: 2004-08-12
- 修回日期: 2005-03-08
- 刊出日期: 2005-07-15

Multi-Monte Carlo Method for General Dynamic Equation Considering Particle Coagulation

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

摘要

摘要: Monte Carlo(MC)方法被广泛用于通用动力学方程的求解,然而普通MC方法的计算代价较高而计算精度不稳定．提出一种新的多重Monte Carlo(MMC)算法来求解GDE,该算法同时具有基于时间驱动MC方法、常数目法和常体积法的特点．首先详细介绍了该算法,包括加权虚拟颗粒的引入,MMC算法的计算流程,时间步长的设置,颗粒是否发生凝并事件的判断,凝并伙伴的寻找,凝并事件的后果处理．然后利用MMC算法对存在理论分析解的5种特殊工况进行数值求解,模拟结果与理论解符合很好,证明MMC算法具有良好的计算精度和较低的计算代价．最后分析了不同类型的凝并核对于凝并过程的影响,常凝并核和连续区布朗凝并核对小颗粒影响大一些,而线性凝并核和二次方凝并核对大颗粒影响大一些．
- 数值算法 /
- 粒数衡算方程 /
- 虚拟颗粒 /
- 颗粒尺度分布 /
- 凝并核 /
- 计算代价 /
- 计算精度
Abstract: Monte Carlo (MC) method is widely adopted to take into account general dynamic equation (GDE) for particle coagulation, however popular MC method has high computation cost and statistical fatigue. A new Multi-Monte Carlo (MMC) method, which has characteristics of time-driven MC method, constant number method and constant volume method, was promoted to solve GDE for coagulation. Firstly MMC method was described in details, including the introduction of weighted fictitious particle, the scheme of MMC method, the setting of time step, the judgment of the occurrence of coagulation event, the choice of coagulation partner and the consequential treatment of coagulation event. Secondly MMC method was validated by five special coagulation cases in which analytical solutions exist. The good agreement between the simulation results of MMC method and analytical solutions shows MMC method conserves high computation precision and has low computation cost. Lastly the different influence of different kinds of coagulation kernel on the process of coagulation was analyzed: constant coagulation kernel and Brownian coagulation kernel in continuum regime affect small particles much more than linear and quadratic coagulation kernel, whereas affect big particles much less than linear and quadratic coagulation kernel.
- numerical solution /
- population balance equation /
- fictitious particle /
- particle size distribution /
- constant kernel /
- computation cost /
- computation precision

HTML全文

参考文献(13)

[1]	Tucker W G. An overview of PM 2.5 sources and control strategies[J].Fuel Processing Technology,2000,65(1):379—392. doi: 10.1016/S0378-3820(99)00105-8
[2]	Meng Z,Dabdub D,Seinfeld J H. Size-resolved and chemically resolved model of atmospheric aerosol dynamics[J].Journal of Geophysical Research,1998,103(3):3419—3435. doi: 10.1029/97JD02796
[3]	Debry E,Sportisse B,Jourdain B. A stochastic approach for the numerical simulation of the general dynamics equation for aerosols[J].Journal of Computational Physics,2003,184(2):649—669. doi: 10.1016/S0021-9991(02)00041-4
[4]	Liffman K. A direct simulation Monte Carlo method for cluster coagulation[J].Journal of Computational Physics,1992,100(1):116—127. doi: 10.1016/0021-9991(92)90314-O
[5]	Kruis F E,Maisels A,Fissan H. Direct simulation Monte Carlo method for particle coagulation and aggregation[J].AICHE Journal,2000,46(9):1735—1742. doi: 10.1002/aic.690460905
[6]	Smith M,Matsoukas T. Constant-number Monte Carlo simulation of population balances[J].Chemical Engineering Science,1998,53(9):1777—1786. doi: 10.1016/S0009-2509(98)00045-1
[7]	Lee K,Matsoukas T. Simultaneous coagulation and break-up using constant-N Monte Carlo[J].Powder Technology,2000,110(1/2):82—89. doi: 10.1016/S0032-5910(99)00270-3
[8]	Lin Y,Lee K,Matsoukas T. Solution of the population balance equation using constant-number Monte Carlo[J].Chemical Engineering Science,2002,57(12):2241—2252. doi: 10.1016/S0009-2509(02)00114-8
[9]	Nanbu K. Direct simulation scheme derived from the Boltzmann equation Ⅰ-Monocomponent gases[J].Journal of the Physical Society of Japan,1980,49(11):2042—2049. doi: 10.1143/JPSJ.49.2042
[10]	Nanbu K,Yonemura S. Weighted particles in coulomb collision simulations based on the theory of a cumulative scattering angle[J].Journal of Computational Physics,1998,145(2):639—654. doi: 10.1006/jcph.1998.6049
[11]	ZHAO Hai-bo,ZHENG Chu-guang,XU Ming-hou. Multi-Monte Carlo method for particle coagulation:description and validation[J].Applied Mathematics and Computation,2005,168(2):.
[12]	Friedlander S K,Wang C S. The self-preserving particle size distribution for coagulation by Brownian motion[J].Journal of Colloid and Interface Science,1966,22(1):126—132. doi: 10.1016/0021-9797(66)90073-7
[13]	Lee K W,Lee Y J,Han D S. The log-normal size distribution theory for Brownian coagulation in the low Knudsen number regime[J].Journal of Colloid and Interface Science,1997,188(2):486—492. doi: 10.1006/jcis.1997.4773

施引文献

资源附件(0)

访问统计

计量

文章访问数: 2546
HTML全文浏览量: 109
PDF下载量: 1371
被引次数: 0

姓名
邮箱
手机号码
标题
留言内容
验证码

留言板

求解考虑颗粒凝并的通用动力学方程的多重MonteCarlo算法

作者简介:
赵海波(1977- ),男,湖南宁乡人,讲师,博士(联系人.Tel/Fax:+86-27-87545526;E-mail:klinsmannzhb@163.com).

计量

Multi-Monte Carlo Method for General Dynamic Equation Considering Particle Coagulation

计量

目录

留言板

求解考虑颗粒凝并的通用动力学方程的多重MonteCarlo算法

作者简介: 赵海波(1977- ),男,湖南宁乡人,讲师,博士(联系人.Tel/Fax:+86-27-87545526;E-mail:klinsmannzhb@163.com).

计量

出版历程

Multi-Monte Carlo Method for General Dynamic Equation Considering Particle Coagulation

计量

出版历程

目录

作者简介:
赵海波(1977- ),男,湖南宁乡人,讲师,博士(联系人.Tel/Fax:+86-27-87545526;E-mail:klinsmannzhb@163.com).