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

LU并行求解算法的研究与实现 徐 金 秀 (江南计算技术研究所 无锡214083) corn) (jinxiu3186@sina successive over-relaxation对称逐次超松弛)求解方 仿真CFD中N．S方程的SSOR(symmetric 法，并针对该算法数据相关性的特点给出了并行求解的不同镱略及相应的分析，最后给出了 一些测试结果。 关键词并行算法、计算区域、数据分布、数据相关性、PE 1． 引言 从70年代开始世界上相继出现了多种不同类型的并行计算机，运算速度从当初的每秒 200万次浮点运算达到了今天的每秒十万亿次级浮点运算。这些超级并行计算机的出现不仅 为解决许多大型科学计算问题提供了硬件基础，而且对研究能充分发挥这些机器并行处理能 力及解决具有重大挑战意义的科学工程计算的并行算法的发展也起到了很大的促进作用。 三维Navier-Stokes(简称N．s)方程的求解在数值天气预报、计算流体力学、石油地震勘 探、航空航天等领域都有着广泛的应用，特别是由于超级计算机的出现，解决涉及数据量大、 运算量大、精度要求高等问题的计算已成为可能，目前是国际上研究的热点。本文以美国 NASA Ames研究中心提出的NASbenchmark中CFD(计算流体动力学)仿真应用Lu课题为蓝本， 通过对串行程序及数据相关性的分析，采取相应的数据分配策略，从计算过程中抽取高度的 并行机制，取得了较好的并行性能。 2．数学模型 2．1差分方程 可压缩流体与时间相关的近物面Navier—Stokes方程可表示如下： auOE(U)OF(U)go(u) af a善 a叩 af —。…一一 ．OT(U，Uf)．ov(u，％)．ow(u，％) 8{8q8‘ (1) +H(U，Uf，U口，Uf)， (f，{，■，f)∈D，xD 边界条件和初始条件分别满足下式 (r，f，刁，f)∈见xSD B(U，％，％，％)=U8(r，；，口，f)， U=Uo(孝，刁，f，)， (≠，野，f)∈D 当r=0 其中：U、E、F、G、T、V、w和H是5xl的向量函数；‘、n、‘是可计算区域的一般曲线 Or-／}。 坐标的三个方向，f是时间变量，D∈吼3是有界区域，aD是D的边界，口={『f ＼1 零 1殳用一个仞日_【ElJ琏垃歹刊基T伺呲斤升近似驯芏川尚戢拔uJ水侍力栈【1)剀尚靓逝1以 解。在差分过程中为了保证解的迭代收敛，要包含个辅助高阶耗散项，当仅对辅助高阶耗 13“A 、0E 8￡z 8n 8nz 8(8￡1“ ：△，f塑塑：．+塑±型+塑±型] 时，上 ‘ 。 a亭 a叩 af 述方程 出叫等+嘭等+酵等Ⅲ谢+ ∽ 差盏萎 其中H’为修改过的向量限制函数，上述二式左部(LHS)被称为impl