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

两相流系统的离散玻尔兹曼建模+ 许爱国}}l，2，3，张广财l，3，甘延标4 l北京应用物理与计算数学研究所国家级计算物理重点实验室， 北京100088 2．北京大学应用物理与技术研究中心和高能量密度物理数值模拟教育部重点实验室，北京100871 3、理论物理国家重点实验室(中国科学院理论物理研宄所)，北京100190 4一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， Boltzmann Discrete Bolt,maannMethod) Model)f2-3|。离散Boltzmann方法(DBM．Discrete 的最大特点就在于它继承了Boltzmann方程描述的跨尺度特性。通过它，人们可以同时研究 动力学非平衡和热力学非平衡。它不