一个三维多松弛时间全速域格子Boltzmann模型-计算物理.PDF

第３５卷 第４期 计　 　 算　 　 物　 　 理 Ｖｏｌ．３５，Ｎｏ．４ ２０１８年７月 　 ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＣＯＭＰＵＴＡＴＩＯＮＡＬ ＰＨＹＳＩＣＳ　 Ｊｕｌ． ，２０１８ 文章编号：１００１⁃２４６Ｘ（２０１８）０４⁃０３７９⁃０９ 一个三维多松弛时间全速域格子Ｂｏｌｔｚｍａｎｎ模型 １，∗ ２，３，∗ ２ １ 陈　 锋 ，　 许爱国 ，　 张广财 ，　 焦培刚 （１．山东交通学院，济南　 ２５０３５７；　 ２．北京应用物理与计算数学研究所计算物理国家重点实验室，北京　 １０００８８； ３．北京大学应用物理与技术研究中心和高能量密度物理数值模拟教育部重点实验室，北京　 １００８７１） 摘　 要：构建一个既适用于低速不可压流体又适用于高速可压缩流体的三维自由参数多松弛时间格子Ｂｏｌｔｚｍａｎｎ 模型．模型中，根据ＳＯ（３）群的不可约表述基函数构造转化矩阵，根据恢复可压Ｎａｖｉｅｒ⁃Ｓｔｏｋｅｓ方程的需要选取非守 恒矩平衡值．通过ｖｏｎ Ｎｅｕｍａｎｎ稳定性分析模型参数对数值稳定性的影响，并给出建议选择范围．模型经过基准问 题的验证，模拟结果与解析解及其它数值结果符合较好． 关键词：格子Ｂｏｌｔｚｍａｎｎ；可压缩流体；多松弛时间 中图分类号：Ｏ３５ 文献标识码：Ａ ＤＯＩ：１０．１９５９６／ ｊ．ｃｎｋｉ．１００１⁃２４６ｘ．７６８０ ０　 引言 高速可压流体普遍存在于自然界和工程技术的各个领域，如爆轰物理、航空航天、天体物理、生物物理 等．可压流体的数值算法在国防、航空航天、海洋、流体机械、化工及生命科学等高科技领域均有重要的理论 意义和广泛的应用前景．但由于其非线性和非平衡特性，对于可压流系统，特别是内含激波和接触间断系统 的数值模拟仍是一项有意义和富有挑战性的工作． ［１－６］ 格子Ｂｏｌｔｚｍａｎｎ （ＬＢ）理论 基于非平衡统计物理基本方程Ｂｏｌｔｚｍａｎｎ方程，适用于连续介质假设不再 成立的非平衡系统．但是作为计算流体领域的新方法，ＬＢ理论一直受到数值稳定性的约束．ＬＢ方法的稳定 ［７］ ［８］ 性问题已经受到许多学者的关注，也作出很多努力和尝试．Ａｌｅｘａｎｄｅｒ 、俞慧丹 等采用调节声速的办法减 ［９］ 轻小马赫数约束的影响，所建模型适用于近似等温的可压流体系统． 闫广武等 提出多速度多能级的二维 ［１０］ ［１１］ 可压模型构造方案．孙成海等 提出局部自适应模型，模型可以适用于较宽范围的马赫数．单肖文等 提 ［１２］ 出根据Ｈｅｒｍｉｔｅ 函数推导分布函数构造可压模型的理论．舒常等 利用圆函数代替麦克斯韦分布函数构造 ［１３］ 二维高速可压模型．何雅玲等 通过耦合双分布函数方法也构造了一系列高速可压ＬＢ模型．冉政课题组基 于冲击波界面的捕捉技术、Ｌｉｅ 对称守恒和熵等思想在可压流体系统的ＬＢ 方法方面取得了一系列进 展［１４－１５］．Ｔｓｕｔａｈａｒａ 等［１６－１８］提出并发展了有限差分可压ＬＢ 建模思路．我们课题组［１９－２２］通过补充物理粘性 在高马赫可压流体的ＬＢ 方法方面也取得了一系列进展．最近，我们通过矩关系反求离散速度空间的平衡态 分布函数，在不使用额外黏性项的条件下实现了高速可压流体的模型构建［２３－２４］．这些方案多采用单松弛因 子模型即ＢＧＫ模型．在原始的ＢＧＫ 模型中，所有动理学模式趋于热力学平衡的速率都用一个常系数来描 述，它不随系统中密度、压强