【精选】离散相系统群体平衡模型的求解算法_苏军伟.pdfVIP

中国科学: 化学 2010 年 第40 卷 第2 期: 144 ~ 160 SCIENTIA SINICA Chimica SCIENCE CHINA PRESS 评述 离散相系统群体平衡模型的求解算法 ① ②* ③ 苏军伟 , 顾兆林 , XU X.Yun ① 西安交通大学机械工程学院机械工程及其自动化系, 西安 710049; ② 西部灾害与环境力学教育部重点实验室, 兰州 730000; 西安交通大学人居学院环境科学与技术系, 西安710049; ③ Department of Chemical Engineering, Imperial College London, London, SW7 2AZ * 通讯作者, E-mail: guzhaoln@ 收稿日期: 2008-12-30; 接受日期：2009-02-26 摘要 准确预测离散相系统中微观粒子的尺度分布演变对系统动态流动行为的准确确 关键词 定起关键性作用. 粒子的尺度分布演变以及引起尺度分布变化的离散相微观行为(聚并、 群体平衡方程 直接离散方法 破碎、长大等) 由群体平衡模型来描述. 该模型是关于数值密度函数的非线性双曲型方程, Monte Carlo 方法 数值求解为主要手段. 本文对群体平衡方程的直接离散方法、Monte Carlo 、矩方法从实现 矩方法 难易程度、计算机资源消耗、计算精度三方面进行了详细阐述, 并着重介绍了几种性能优 离散相系统 越的矩方法——矩积分方法(QMOM) 、矩直接积分方法(DQMOM) 、可调节矩积分方法 (M-QMOM)、自适应矩直接积分方法(ADQMOM)、定点矩积分方法(FPQMOM)、粒子游 动算法(MPEM)和局部定点矩积分方法(LFPQMOM). 最后根据算法的优缺点及其当前发 展状况对不同算法的未来发展做了预测. 1 引言 力和相间传质封闭至关重要, 这使得该方程成为研 究离散相体系中介观结构对宏观结构影响重要工具. 两相、多相共存的离散相系统广泛存在于工业和 群体平衡模型本身描述的离散相微观行为的微观频 自然界中. 双流体模型以其能够得到系统的宏观流 度项需要从微观进行封闭, 所以该方程为离散相系 动结构、计算量低等优点在科学研究和工业实际模拟 统架起了一个从微观到宏观的桥梁. 群体平衡方程 中广泛采用. 该模型是通过单相流动的N-S 方程通过 是关于数值密度函数的连续形式, 该方程为非线性 某种方式平均(体积平均、时间平均、总体均值或体 双曲型方程, 且通常求解域为半无穷域, 没有一般意 积时间平均等)得到的, 这些平均过程使双流体模型 义的分析解, 数值求解成为研究该方程的主要手段[1]. 存在着与生俱来的封闭问题——湍