Euler方程的双分布函数格子.pdfVIP

第 15 卷 第 2 期 　 计 　　算 　　物 　　理 　　 Vol . 15 ,No . 2 1998 年 3 月 　 CH IN ESE J OU RNAL OF COM PU TA T IONAL PH YSICS 　　 　Mar . , 1998 　 Euler 方程的双分布函数格子 Boltz mann Godunov 方法 阎广武 ( 中国科学院力学研究所 LNM 开放实验室 ,北京 100080) 胡守信 　施卫平 ( 吉林大学数学系力学教研室 ,长春 130023) 摘 　要 　应用双分布函数系统 ,通过 Godunov 分解 ,构造了一维 Euler 方程的格子 Boltzmann 算法 。 解决了传统格子气固有的 GC 问题与能量方程之间的矛盾 ,实现了分布函数与宏观物理量之间的 一一对应 。 关键词 　Euler 方程 　格子 Boltzmann 方法 　双分布函数 　Godunov 分解 中图分类号 　O24 1 　O354 0 　引　言 格子Boltzmann 方法是一种用虚拟气体的运动模拟流体力学方程的方法[ 1～3 ] , 虽然其 自 ( ) 身带有压缩性质 本身是气体 ,但却被限制在低 Mach 数流动的范围内。最近几年 ,人们开始 ( ) [4～6 ] 研究可压缩流动 特别是含有激波和接触间断问题 的格子 Boltzmann 方法 。然而用传统 格子 Boltzmann 方法模拟可压缩的气体流动存在一些困难 : ( ) D , 虽解除了限制 ,但是对激波的捕获能力较低 , 即 1 低速度限制 ,可以消除压缩因子 ij 此类模型仅适用于无激波的可压缩流动[7 ] 。 ( ) ( ) [ 8 ] γ 2 所谓的 IC Ideal Case 问题限制 ,二阶矩 自身的要求 ,导致了比热比在一维时 = 3 (二维时 γ ) 的理想情况 。通常我们研究的气体的比热比为 14 ,无 IC 限制是构造可压缩 = 2 模型的一个基本要求 。 (3) GC ( General Case) 问题的引入是通过定义第三个维度上的振动能量 ε [ 8 ] 来实现的 ,然 而却得到一个附加的关于该振动能的守恒方程 。目前构造高精度高分辨的模型尚在研究之 中。 (4) 从传统的格子 Boltzmann 方程到 目前的几个可压缩模型 ,均存在一个共同的问题 , 即 (Δ 2 2) 模型的精度与稳定性 。从模型的构造上我们发现 ,一般的精度是 O x , K ,这里的 K 是 Knudsen 数 ,其分析的基础是建立在下面的假设上 : (a) 小的 Knudsen 数 , (b) 低速度 。在流体 ( ) 力学中我们采用如下的表达式 K = M / R e , M 为 Mach 数 , R e 为 Reynolds 数 ,当 R e 数不能