两相流系统的离散玻尔兹曼建模.docx

两相流系统的离散玻尔兹曼建模+许爱国}}l，2，3，张广财l，3，甘延标4l北京应用物理与计算数学研究所国家级计算物理重点实验室，北京100088 2．北京大学应用物理与技术研究中心和高能量密度物理数值模拟教育部重点实验室，北京1008713、理论物理国家重点实验室(中国科学院理论物理研宄所)，北京1001904一b华航天工业学院，廊坊065000摘要：本文借助离散玻尔兹曼方程构建液气二相流系统的非平衡相变模型。通过该模型， 人们可以同时研究二相流系统的动力学非平衡与相变动理学过程中的热力学非平衡。讨论了测量系统热力学非平衡程度的方案以及研究结果对宏观流体建模的贡献。关键词：多相流、非平衡相变、可压流体、格子玻尔兹曼、离散玻尔兹曼 1．引言非平衡与多相复杂流体系统形式多样，广泛存在于自然界和工程物理等领域。复杂流体 系统的最显著特征之一就是系统中存在大量的界面。大量界面的存在与复杂演化必定伴随着 丰富、复杂的非平衡效应。这类系统的物理建模与算法设计均存在较强的挑战性。传统连续 流体模型建立的基础是观测所用的时间尺度远大于系统热力学弛豫所需的时间尺度，所以可 以成功地描述动力学非平衡过程，而对热力学非平衡关注较少；或者从空间尺度上来说，传 统连续流体模型建立的基础是所关心的特征结构的尺度远大于“微观粒子”之闽的平均间距 (1000倍以上)，所以它描述的主要是大量“微观粒子”的平均效应，而对相对于平均效应的涨落关注较少。系统的状态是客观的，而“系统处于平衡与非平衡”以及“非平衡程度”的认识却是依赖于我们观测所用的时间和空间尺度的。在等离子体系统、包含流体不稳定性或非平衡相变等效应的系统中，存在着大量的不同 时间．空间尺度的结构和动理学过程。在爆燃与爆轰等系统中，与化学反应相伴随的放热、压缩甚至冲击等过程的介入使得系统行为的复杂程度进一步剧增。这类系统需要我们从不同 的空间尺度、不同的时间尺度上去理解、分析和综合。随着对系统粗线条、主要特征的逐渐 掌握，随着对冲击、爆轰等复杂物理过程模拟精度要求的提高，系统或材料内的各种微介观 结构与非平衡效应日益引起人们的关注。在最近的二十几年间格子玻尔兹曼(LB，LatticeBoltzmann)方法及其在不同类型的复杂流体系统模拟研究方面取得了长足的进展11．31。根据离散格式的不同，LB方法可以粗略地分为标准LB方法【l】、有限差分LB方法【4卅、有限体积LB方法【7】和有限元LB方法【8．9]，等等滞l据目的和功能，LB方法可以粗略地划分为偏微分方程解法器和微介观动理学模型。 第一类LB模型最初作为Navier-Stokes等流体方程的一种全新的渐进式、逼近式数值解法；目前这种全新的求解思路己经被推广应用于一系列微分方程的求解。第二类LB方法将LB 看作是一种全新的物理建模；它要求LB建模包含所要研究的物理因素，它希望LB微介观 模型具有某些传统连续流体模型所不具有或不方便处理的功能；这些功能以描述或度量热力学非平衡程度为主【1．3】。作为偏微分方程解法器的LB模型，其构建规则相对灵活：只要最+本工作受国家自然科学基金(批准号、计算物理重点实验基金、爆炸科学与技术国家重点实 验室基金开放基金(批准号：KFJJl4．IM)和理论物理国家重点实验室(中国科学院理论物理研究所开放课题f编号：Y4KFl51CJI)资助。 ”许爱国，1970年生，男，研究员。主要研究方向：冲击与爆轰相关的微介观结构与非平衡效应研究。E—mail：苎!尘』g!垃f《j§匹地：4￡：鲤122终的数值结果与所求解的方程一致便可。这就意味着，计算过程中的某些处理无需过多地考 虑其物理对应。鉴于目前文献中的LB模型以第一类为主，且占了绝大多数，在本文中，那些作为复杂系统微介观、粗粒化物理建模出现的LB模型称为离散玻尔兹曼模型(DBM，DiscreteBoltzmannModel)f2-3|。离散Boltzmann方法(DBM．DiscreteBolt,maannMethod)的最大特点就在于它继承了Boltzmann方程描述的跨尺度特性。通过它，人们可以同时研究动力学非平衡和热力学非平衡。它不介入分子运动细节，但又与分子动力学描述类似，可以 直接从微介观角度考虑动力学和热力学过程的耦合，是用来研究这些微介观结构和非平衡效 应的有效工具之一。作为复杂系统微介观粗粒化物理模型的LB方法即DBM在近年来取得了快速的进展。 通过DBM同时研究流体系统的宏观动力学行为和热力学非平衡行为的思想最初出现在综述文章121r9，进而在随后的几个模型中得到了具体化和实际应用【10．14】。单松弛因子DBM的例子可参见【11．13】，多松弛因子的DBM参见1141。在文献110，13lqb，实现了爆燃和爆轰问题的DBM建模与模