从稀薄流到连续流的HPF并行算法的研究.pdfVIP

2003空气动力学前沿研究论文集 从稀薄流到连续流的HPF并行算法研究。 牵志辉张涵信 (国家计算流体力学实验室，北京100083 中国空气动力研究与发展中心，四川绵阳621000) 摘要本文从修正的BGK．Boltzmann模型方程出发，引入离散速度坐标技术对气体 分子速度分量进行离散降维，基于非定常时间分裂数值计算方法和无波动、无自由参 数的NND耗散差分格式，发展直接求解气体分子速度分布函数的气体运动论有限 差分数值格式，提出一套能有效模拟各流域三维绕流问题的气体运动论统一算法。 在分析研究气体运动论数值算法内在并行度的基础上，开展各流域三维绕流问题统 一算法的HPF(高性能FORTRAN)并行化程序设计，建立一套能有效模拟各流域复 杂外形体绕流的唧F并行算法软件。通过对不同Knudsen(克努森)数下的三维球 体绕流及类“神舟号”返回舱外形体绕流的初步数值实验看出，计算结果与过渡区有 关实验数据及各流域的气体绕流现象一致性很好，证实了本文发展的气体运动论 HPF并行算法在求解稀薄流到连续流不同流域复杂绕流问题方面的可行性。 关键词 气体运动论 Boltzmann模型方程离散速度坐标法 NND差分格式 HPF并行计算 1 引言 航天飞行器再人过程先后经历自由分子流区、过渡流区、滑流和连续介质流区，所遇气动 环境十分复杂。在稀薄流到连续流各流域空气动力学问题计算研究中，传统的做法是稀薄流 用稀薄流的一套算法，如基于微观分子气体动力学将分子运动、碰撞耦合模拟，依托蒙特一卡 罗直接模拟(DSMC)方法求解稀薄气体流动问题-lt2J，并获得了广泛的应用；而连续流有连续 流的一套算法，如从宏观气体流动角度出发，以气体密度、温度、流动速度作为独立变量，通过 数值求解欧拉(Euler)方程、纳维尔一斯托克斯(N—s)方程发展了许多成熟的数值计算方 法13．4 o；但是两种方法相差较大，且计算结果彼此独立，很难随高度光滑连接。是否能发展一 种连接连续流到稀薄流能有效模拟不同流域气体流动问题的统一算法，这是近年来学术界所 关心的研究课题。 根据国际上关于气体运动论模型方程的研究现状L51J，结合DSMC方法和计算流体力学 有限差分方法的研究基础，文献[8，9]基于修正的BGK．Boltzmann模型方程，发展和应用离散 速度坐标技术于气体分子速度空间，利用一套在物理空间和时间上连续而速度空间离散的分 布函数来代替原分布函数对速度空间的连续依赖性，将各流域统一的气体分子速度分布函数 简化控制方程化为在各个离散速度点处基于时间和位置空间具有非线性源项的双曲守恒方 。国亲自然科学基金重大研究计划赍助项目 从稀薄流到连续流的HPF并行算法研究 程。研究并发展可用于离散速度空间宏观取矩的不同积分方法，并应用于不同马赫数绕流的 离散速度数值积分。基于非定常时间分裂数值计算方法，应用并推广NND格式，发展直接求 解气体分子速度分布函数的气体运动论有限差分格式，已经尝试建立起能有效模拟各流域一 维、二维气体流动问题的简化的气体运动论统一算法。本文将在上述气体运动论数值算法研 究的基础上，吸收气体运动论关于微观分子动力学的研究思想，利用已经发展起来的离散速度 坐标法和非定常时间分裂耦合迭代数值计算技术，依托二阶Runge—Kutta方法和原始变量的 NND一4(a)差分格式【1…，拟发展一套基于气体分子速度分布函数方程能有效模拟各流域三维 绕流问题的统一算法；并依靠我国自行研制的超级并行计算机系统，从分析基于Bo|tzmann模 型方程的气体运动论数值算法内在并行度出发，研究气体运动论数值算法所适合的并行方案， 起来的气体运动论统一算法进行HPF并行化程序设计，通过对不同流域的圆球绕流和类“神 舟号”返回舱外形体绕流的并行计算，以期建立一套基于Boltzrnann模型方程能有效模拟各流 域气体流动问题的可供实际应用的HPF并行算法软件。 2控制方程和数值方法 2．1控制方程 基于修正的BGK—Boltzmann模型方程，确定含流态控制参数的可描述不同流域微观分子 输运现象的气体分子速度分布函数简化控制方程[7，1…，其无量纲形式为： ¨ = 也 口