复杂流动与流体行为拟颗粒模拟.pdfVIP

复杂流动与流体行为的拟颗粒模拟 1 1,2 1 1,2 1,2 1,2 1 王利民 ，葛蔚 ，陈飞国 ，侯超峰 ，卢健新 ，张家元 1．中国科学院过程工程研究所多相复杂系统国家重点实验室，北京（100080 ） 2. 中国科学院研究生院，北京（100049 ） E-mail ：wge@ 摘 要：复杂流动和流体行为中的多尺度结构对化学工程的理论和应用研究都是巨大的挑 战。多尺度的粒子模拟是处理复杂界面、物性以及连续介质模型难以适用的纳微和极端条件 下流动的有力手段之一。本文将重点介绍本实验室近年来以拟颗粒模拟为基础的微观粒子模 拟研究及其在流态化和纳微流动中的应用，并在此基础上探索离散模拟通用化的可行性。 关键词：化学工程；多尺度结构；拟颗粒模拟；多相系统；纳微流动 中图分类号：TQ 028 文献标识码：A 自 然 的 方 法 ， 即 分 子 动 力 学 方 法 0. 引 言 [2,3] MD(Molecular Dynamics) 。而宏观上， 化工系统中普遍存在流动、传递与反 离散元方法DEM(Discrete Element Method) 应等多种过程的耦合、气液固等多相流动 [4]也是模拟颗粒物料的直观方法。但更多 的耦合以及高温高压条件等，使系统具有 的粒子方法可以理解为一种计算模型，使 显著的非平衡性，从而表现出典型的动态 用的是某种虚拟粒子，如微观的直接模拟 多尺度结构、多态性和突变现象、放大效 蒙特卡罗方法 DSMC(Direct Simulation 应等[1] 。要严格描述这些多尺度结构通常 Monte Carlo)[5,6]，介观的耗散粒子动力学 有两种途径：宏观结构的形成可用“ 自顶向 方 法 DPD(Dissipative Particle 下”(Top Down, TD)方法追溯到连续介质 Dynamics)[7] ，格子气自动机 LGA(Lattice 模型能够成立的最小尺度，即进行直接数 Gas Automata)[8] ，格子波尔兹曼方法 值模拟(DNS) 。传统上基于网格方法的 LBM(Lattice Boltzmann Method)[9]，以及宏 DNS 在处理运动颗粒或气泡等复杂边界 观的光滑粒子动力学方法 SPH(Smooth 时存在不少困难，比如规则网格中界面的 Particle Hydrodynamics)[10~12] 和流体粒子 描述问题和不规则网格中界面的追踪和网 方法 FPM(Fluid Particle Model)[13] 。本文重 格映射问题。而对介观结构的形成，宏观 点 介 绍 的 拟 颗 粒 模 拟 (Pseudo-particle 连续介质理论已无法适用，目前还缺乏完 modeling, PPM)[14,15] 中也构造了一种变径 善的理论，需要采取“ 由底而上”(Bottom 的微观虚拟粒子，即拟颗粒(PP) 。 Up, BU)方法, 通