LU+并行求解算法的研究与实现研究.pdfVIP

LU 并行求解算法的研究与实现 徐 金 秀 （江南计算技术研究所 无锡 214083） (jinxiu3186@) 摘 要 N-S 方程在 CFD （computational fluid dynamics ）领域有着广泛的应用。本文介绍了 仿真 CFD 中N-S 方程的 SSOR(symmetric successive over-relaxation 对称逐次超松弛)求解方 法，并针对该算法数据相关性的特点给出了并行求解的不同策略及相应的分析，最后给出了 一些测试结果。 关键词 并行算法、计算区域、数据分布、数据相关性、PE 1． 引言 从 70 年代开始世界上相继出现了多种不同类型的并行计算机，运算速度从当初的每秒 200 万次浮点运算达到了今天的每秒十万亿次级浮点运算。这些超级并行计算机的出现不仅 为解决许多大型科学计算问题提供了硬件基础，而且对研究能充分发挥这些机器并行处理能 力及解决具有重大挑战意义的科学工程计算的并行算法的发展也起到了很大的促进作用。 三维 Navier-Stokes(简称 N-S)方程的求解在数值天气预报、计算流体力学、石油地震勘 探、航空航天等领域都有着广泛的应用，特别是由于超级计算机的出现，解决涉及数据量大、 运算量大、精度要求高等问题的计算已成为可能，目前是国际上研究的热点。本文以美国 NASA Ames 研究中心提出的NAS benchmark 中CFD(计算流体动力学)仿真应用 LU 课题为蓝本， 通过对串行程序及数据相关性的分析，采取相应的数据分配策略，从计算过程中抽取高度的 并行机制，取得了较好的并行性能。 2． 数学模型 2.1 差分方程 可压缩流体与时间相关的近物面 Navier-Stokes 方程可表示如下： U E U F U G U ∂ ∂ ( ) ∂ ( ) ∂ ( ) + + ∂τ ∂ξ ∂η ∂ς ∂T U U ∂V U U ∂W U U ( , ξ ) ( , η ) ( , ζ ) + + + ∂ξ ∂η ∂ζ + ∈ × H U U U U D D ( , ξ , η , ζ ), (τ,ξ,η,ζ) τ (1) 边界条件和初始条件分别满足下式 B (U,U ,U ,U ) UB (τ,ξ,η,ς), (τ,ξ,η,ζ) ∈D ×∂D ξ η ς τ 0 U U (ξ,η,ζ,), (ξ,η,ζ) ∈D 当τ 0 其中：U、E、F、G、T、V、W 和 H 是 5×1 的向量函数； ξ、η、 ζ是可计算区域的一般曲线 坐标的三个方向，τ是时间变量，D ∈ℜ3 是有界区域，∂D 是 D 的边界，D {τ | 0≤τ≤T}。 τ