从稀薄流到连续流的统一数值算法的研究.pdfVIP

第十届全囝计算流体力学会议论文集 从稀薄流到连续流的统一数值算法研究‘ 李志辉张涵信 中国空气动力研究与发展中心．四川绵阳621000 捅矍从描述各个流域微观分子输运现象的Bolt=nama(玻尔兹曼1方程出发，利用气件分子对流运动、碰 撞松弛趋于平衡态的基车特性，建立能描述各个流域气体流动特性的气体分子速度分布函数的统一的简化控制 方程．发展和应用高教速度坐标法于速度空问，利用一套在物理空间和时间上连续而在速度空间离敞的分布函 数来代替原分布函数对速度空间的连续依赖性．基于非定常时间分裂数值计算方法和无波动、无自由参数的 NND耗散差分格式．建立直接模拟气体分子速度分布函数的有限差分方法．推广应用Gauss．Hcvmite无穷积分 法和华罗庚一王元提出的蛆单和逼近重积分的鼓论积分方法等，对离散速度空间进行宏观取矩获取物理空间 各点的气体流动参数，由此发展—套从稀薄流到连续流各流域统一的数值算法．通过对不同Knudscn(克努森) 数下一维激波管问西、二维四柱绕流和三维球件绕流数值实验表明本文发展的数值算法是可行的． 0引言 在自由分子流到连续流各流域空气动力学问题计算研究中，传统做法是稀薄流用稀薄流 的一套算法，如Monte Carlo方法；而连续流有连续流的一套方法，如数值求解N—s方程等： 两者相差较大，获得的气动数值彼此独立、难以光滑连接。本文研究的目的在于基于描述各 个流域微观分子输运现象的Boltzmann模型方程，发展一种连接连续流到稀薄流能有效模拟 各个流域气体流动特性的统一算法。 本文从分子运动论的基本方程——B01tzmsl]n方程”’出发，利用气体分子对流运动、碰撞 松弛趋于平衡态的基本特性．建立含流态控制参数的可描述各流域流动特性的统一的气体分 子速度分布函数的简化控制方程。推广和应用离散速度坐标法口J，去掉分布函数对于速度分量 K、蚱，雎的连续依赖性；研究并发展可用于离散速度空间宏观取矩的Gauss—Hermite无穷积 分法和以单和逼近重积分的华罗庚——王元黄金分割积分法闭等；将速度分布函数方程化为在 时间和位置空间连续而速度空间离散的具有非线性源项的双曲守恒方程。基于气体运动论关 于微观分子动力学的研究思想，吸收Monte Carlo方法关于分子运动和碰撞的去耦技术14l，使 用非定常时间分裂数值计算方法，借助无波动、无自由参数的NND耗散差分格式％建立一套 直接求解从稀薄流到连续流各流域统一的气体分子速度分布函数方程数值算法。为了证实本 文数值方法在求解稀薄流到连续流气体动力学问题方面的可靠性，本文已将该数值算法用于 不同Knudsen数下的一维激波管问题、二维园柱绕流和三维球体绕流的数值模拟。 1 各流域统一的气体分子速度分布函数简化控制方程的建立 从描述各个流域微观分子输运现象的B01％zmann方程【11出发．基于分子运动论的碰撞间隔 理论，利用气体分子对流运动、碰撞松弛趋于平衡态的基本特性16]，以Maxwell(麦克斯威尔) 分布函数厶作为权函数的Hermite(埃尔米特)多项式的渐近展开而得到气体分子当地平衡态 分布函数，”，井将各流域气体分子输运特性、熟力学效应及气体分子相互作用规则与分子碰 撞频率v和当地平衡态分布函数相联系m，由此建立含流态控制参数的各流域均适用的气体分 子速度分布函数的无量纲简化控制方程。 ·本工作得到国家自然科学基金项目的瓷助。 265一 ● ” ’4 笪+F．霎；w一一门 8t 舞 ，”：丘．I+O—片声．i∞：／T一班(5删2I (2) fu：一∽y，：唧}-c：川 (3) ，幽r¨／6石点。) ，墨：屯肚 (4)