光滑粒子流体动力学一种并行数值计算方案.PDFVIP

第 29 卷第 1 期 航 天 器 环 境 工 程 2012 年 2 月 SPACECRAFT ENVIRONMENT ENGINEERING 23 光滑粒子流体动力学的一种并行数值计算方案 周 浩，汤文辉，冉宪文，陈 华 （国防科学技术大学 理学院 工程物理研究所，长沙410073 ） 摘要：在三维超高速碰撞数值计算方面，针对三维光滑粒子动力学（SPH ）方法计算量大和耗时长的缺点， 文章提出了一种简单直接、易于编程实现的 SPH 并行计算方案，并简述了该方案的基本思想、任务划分、变量 存储、信息传递以及主要计算步骤。最后利用自编并行程序计算了两个超高速碰撞实例，结果表明：针对几百 万个粒子，在运算速度为每秒 5 万亿次的“银河”计算机上申请 23 个核并行计算，每步约需要 8 s ，加速比约 为 10，并行效率约为50% ，计算时间显著减少。 关键词：超高速碰撞；并行 SPH 算法；并行树搜索算法；加速比；并行效率 中图分类号：O246 文献标识码：A 文章编号：1673-1379(2012)01-0023-04 DOI ：10.3969/j.issn.1673-1379.2012.01.005 0 引言 质量、动量、能量三个守恒方程，利用附近粒子相 应物理量插值得到各个粒子物理量的随体导数，使 随着航天科技的迅速发展以及航天发射任务 拉氏流体力学偏微分方程组变成差分方程组，从而 的急剧增加，大量的空间碎片散布在近地轨道上， 易于数值求解。SPH 中的一个主要问题是计算量比 形成了复杂的空间碎片环境，对航天器的安全构成 较大，耗时较长，所以开展并行计算研究十分必要。 了极大的威胁。如何确保航天器的飞行安全，成为 当前航天器设计中一个十分重要的问题。对超高速 1 SPH 方法并行计算方案 碰撞问题进行理论研究非常困难，而现有的试验技 并行编程比较复杂，面临着执行结果的不确 术又很难达到超高速碰撞所需要的速度条件，因此 定、如何在各个进程间合理地划分任务、如何减少 数值模拟是一种重要的研究手段。有限元方法作为 通信，以及对于不同进程数（即不同数量的核）具 当前计算力学领域主要的数值计算方法，经过多年 有伸缩性等问题[4-6] 。本文从实际需要出发，采用 的发展，日渐成熟和完善，被广泛应用于解决工程 基于 MPI （Message Passing Interface ）的信息传递 问题。但是在面对超高速碰撞等涉及到大变形的问 并行程序设计平台和粒子分割算法。每个进程中都 题时，应用有限元方法经常会产生网格畸变，严重 有 array1 和 array2 两个数组，分别用来存储所有粒 影响计算精度，甚至导致计算终止。