多相空间数值模拟并行化的探究.pdfVIP

多相空间数值模拟并行化研究 陆林生”董超群”李志辉2’ ”(江南计算技术研究所无锡214083) 0008 2’(国家计算流体力学实验室北京1 3) 摘要 以三维绕流气体运动论数值算法为例，研究多相空间数值模拟并行计算。围绕区域分解方法、 向量归约技术和边界处理并行优化策略等方面进行了讨论，利用并行程序概念设计系统(PPCDS)完成了 HPF并行程序设计，在高性能计算机上进行了大规模并行计算，获得了较好的并行效率。 关键词多相空间，数据并行，区域分解，向量归约，边界处理，并行优化 1引言 三维绕流气体运动论数值模拟基于从稀薄流到连续流的气体运动论统一数值算法，是由中国空气 动力研究与发展中心张涵信院士指导下发展起来的一种新算法，其数学模型和计算方法在文献[1]中 进行了详细论述。本文以此为例着重讨论多相空间并行化技术。 三维绕流气体运动论数值模拟的计算空间是由离散速度坐标空间和位置空间组成的六维空间，可 形成二相空间。算法分成两个主要部分，一部分是在各个离散速度点处基于时间和位置空间求解关于 气体分子速度分布函数厂的具有非线性源项的双ttt守恒方程：另一部分是基于速度分布函数厂，使用 沧，下面给出模拟算法的主要计算式。 1．1 离散速度分布函数的控制方程和差分方程 设(X，Y，z)是直角坐标系的空间位置坐标，(t，n，‘)是变换坐标，J是坐标变换的雅可比行 列式，则控制方程为： 型+望+塑+塑：s at a{8q85 本课题得到国家自然科学基金(批准号资助。陆林生，男，1942年生，教授，博士生导师，主要研究领域为并行{5{刖与并行计算 冉法，并行训算环境等。 董超群，男，1972年生，博士，主要研究领域为并行计算环境、CAD／CASE。 李志辉，男，1968年生．博士 历，副研究员，主要研究领域为稀薄气体力学与计算流体力学。 了2 V=Zx。仉+0a仉+K一玑W=圪。六+％d乞+t一￡．■。，■一，屹一分别是离散速 度坐标O-，J，0处的气体分子速度分量；S=V(∥≥，口一厶正口)，V是气体分子碰撞频率。 并且， 心扩丽P eXp(一和+(1_Pr)e啄2c2肛5)／(5删2)] C=(Vx口一％)f+(矿，d—y。)，+(吃曰一U：)k 面=qxi+qyJ+qzk 用差分法求解方程(1．1)，在‘，n，0方向的网格点分别是i=l，2，．．，N． j=1，2，．．，N“k=l，2…，N。．在离散速度坐标o，6，0方向的网格点分别是O=1，…．N。： 6=1，2，．N s；o=1，2，．．，N。．利用非定常时间分裂法，(1．1)式的差分方程为： U“1 …(1．2) 2t(血／2)‘(出／2)厶(垃／2)乓(出)厶(＆／2)‘(血／2)t(m／2)u” 3)一(1．6)的差分算子： 其中，￡，，Lf，￡。，L。分别代表式(1 —OU—：5 ……(13) ’ Ot 丝+堡：o (1．4) 研 西 丝+塑：o (15)