计算物理 ›› 2008, Vol. 25 ›› Issue (1): 65-74.

• 研究论文 • 上一篇    下一篇

求解Boltzmann模型方程的气体运动论大规模并行算法

李志辉1,2, 张涵信1,2   

  1. 1. 国家计算流体力学实验室, 北京 100083;
    2. 中国空气动力研究与发展中心超高速研究所, 四川 绵阳 621000
  • 收稿日期:2006-09-11 修回日期:2007-01-30 出版日期:2008-01-25 发布日期:2008-01-25
  • 作者简介:李志辉(1968-),男,四川眉山,研究员,博士后,从事稀薄气体动力学与计算流体力学研究.
  • 基金资助:
    国家自然科学基金(90205009,10321002)资助项目

Massive Parallelization of Gas-kinetic Algorithm for Boltzmann Model Equation

LI Zhihui1,2, ZHANG Hanxin1,2   

  1. 1. National Laboratory for Computational Fluid Dynamics, Beijing 100083, China;
    2. China Aerodynamics Research and Development Center, HAI, Mianyang 621000, China
  • Received:2006-09-11 Revised:2007-01-30 Online:2008-01-25 Published:2008-01-25

摘要: 研究各流域三维流动问题的Boltzmann模型方程计算方法,建立直接求解分子速度分布函数的气体运动论耦合迭代数值格式;基于变量依赖关系、数据通信与并行可扩展性分析,使用区域分解并行化方法,建立气体运动论数值算法并行方案,发展求解各流域三维绕流问题的气体运动论并行算法.拟定高低不同马赫数下来自不同流域的三维球体及返回舱绕流算例,进行高性能Fortran(HPF)大规模并行计算,将计算结果与有关实验数据、相关理论预测等进行比较分析,研究揭示不同流区复杂绕流现象及流动机理.研究表明,所发展的气体运动论并行算法具有很好的并行独立性,基本达到线性加速的并行效果,显示出良好的并行可扩展性.

关键词: Boltzmann模型方程, 速度分布函数, 离散速度坐标法, 有限差分方法, 并行计算

Abstract: In a numerical study of Bohzmann model equation, a gas-kinetic finite difference scheme with coupling and iteration is constructed to solve molecular velocity distribution function directly. The parallel strategy is established by using parallel technique of domain decomposition based on variable dependency relation, data communication and parallel expansibility. Gas-kinetic HPF(High Performance Fortran) parallel algorithm is developed to solve three-dimensional problems in various flow regimes. Hypersonic gas flows around a sphere and a spacecraft at various Knudsen numbers, Mach numbers and flying angles are computed at a high performance computer with massive scale HPF parallel. The computational results are in good agreement with experimental and theoretical ones. It is shown that the parallel speed-up increases approximately linearly with the numbers of processors. It indicates high parallel efficiency and expansibility with good load balance and data communication. It suggests that a gas-kinetic parallel algorithm on large scale can be used for three-dimensional complex hypersonic flow problems in various flow regimes.

Key words: Boltzmann-model equation, velocity distribution function, discrete velecity ordinate method, finite difference scheme, parallel computing

中图分类号: