一种新高效并行算法.pdfVIP

282 包 芸等：一种新的高效并行算法 一种新的高效并行算法 包 芸 中山大学应用力学与工程系，510275 廖新浩 李力刚 中国科学院国家天文台，200030 摘 要 一种构造求解 3D 偏微分方程的高效并行计算算法的新思路。从差分格式上寻找切缝 分块进行无信息交换的并行差分计算，而后在处理器间同时进行切面上的信息交换。基于此思 想构造的半隐式切瓜并行算法，对天体力学地球坐标系中的单一热扩散方程行试验计算，实际 并行效率高达 95% 以上。 关键词 并行计算；数值模拟；天体力学 中图分类号 P14 随着科学发展到二十一世纪的今天，大型科学数值模拟计算在海洋、天文、环境等多个 领域中占有越来越重要的分量和位置。计算机的飞速发展也为解决这些领域中描述各种问题 的非线性、非定常偏微分方程组的数值求解提供了很好的硬件工具，使得各种大尺度、长周 期的复杂问题有了数值求解的可能。然而超快速运算的计算都是以并行计算来提高速度的。 如何进行并行计算，寻找高效的并行算法，成为了阻碍大尺度、长周期复杂问题数值计算发 展和应用的“瓶颈”，也是目前发展高性能计算方法的一个重要课题。以往的并行算法大都是 由计算力学的专家们选好差分格式、计算模式后，由专门的并行专家对已有计算方法进行并 行化处理，并行效率往往不尽人意[1~3] 。我们不妨换一个角度去考虑并行计算问题。由计算力 学专家在考虑差分格式、计算模式的同时就把并行计算作为一个目标，构造出几乎完全并行 的算法。基于这种思路，我们在研究地球外核流体力学的长周期数值模拟计算中，构造出地 球坐标系中的半隐式切瓜并行算法，并以单一热扩散方程数值求解对其并行效率进行验证。 1 基本方程 地球外核动力学研究区域为一有内核存在的球壳体 （如图1）。其中内半径为r ，外半径r i 0 且： 包 芸等：一种新的高效并行算法 283 η 1 r , r i 0 1−η 1−η η 为一参数，本文中η =0.4 。 r rii r0 取球体坐标系(r,θ ,φ ),r 为半径方向，θ 为北极到极半径与极轴的 r o 夹角，φ 为绕极轴的夹角。 地球外核热对流问题是由一组较复杂的偏微分方程组[4]来描述的。 图 1 研究区域 为了检验半隐式切瓜并行算法的并行效率，我们选取地球坐标系中较 简单的单一热扩散方程为计算算例。 设扩散方程为 ∂H 2 =∇ +