半隐式半Langrange+时间积分及其可扩展并行算法设计.pdfVIP

半隐式半 Langrange 时间积分及其可扩展并行算法设计 张卫民、曹小林 （国防科技大学计算机学院，wmzhang104@163.net ） 摘 要 目前半隐式半 Langrange 谱模式仍然是全球天气预报业务模式的主流。本文首先介绍已被 广泛应用到业务模式中的两个时间层半 Langrange 时间积分格式，然后研究了应用到半 Langrange 模式中的与三次插值方法具有相同精度的准三次插值方法，最后针对机群系统研究了两个时间层半 Langrange 时间积分格式的可扩展并行算法设计，该方法在由双 CPU 结点组成的 PC linux 机群环境 下取得了良好的运行效果。 关键字 半 Langrange ；可扩展并行算法；linux 机群系统 1 引言 所谓半隐式格式，是指在同一方程中，对激发快波的项用隐式表示，对描述慢波的项用显示表 示所构成的差分格式。由于激发快波的项有时不全是线性的，为了求解方便，只取其中的线性部分 用隐式表示。而半拉格朗日法是指在整个积分过程中不对同一气块沿其路径追踪，而是在同一时间 步长内追踪终点总是在网格点上的气块。半 Langrange 方法与 Euler 方法相比具有较好的计算稳定性 和较高的计算精度，半 Langrange 方法可以加长积分的时间步长，从而提高数值天气预报模式的积 分效率。如欧洲中期数值天气预报中心（ECMWF ）的T213L31 模式，若采用 Euler 方法求解，其时 间步长一般是取为 3 分钟，而若采用三个时间层的半 Langrange 格式则时间步长可取为 15 分钟，若 采用二个时间层的半 Langrange 格式则可取为 30 分钟，而两者进行 10 天的预报结果几乎是一致的。 在数值天气预报中的半 Langrange 方法的发展始于 80 年代，Robert （1981）首先将半Langrange 格式与半隐式格式结合起来，提出了三个时间层的半隐式半 Langrange 格式，并论证了这种方法的 时效优势。Temperton(1987)和 McDonald （1987）等提出了两个时间层的半Langrange 格式，进一步 提高了半 Langrange 格式的数值积分效率。ECMWF 于 1991 年在业务系统中实现了三个时间层的半 Langrange 谱模式，1996 年 ECMWF 进一步将两个时间层的半 Langrange 应用于谱模式中，目前时 间积分采用半隐式半Langrange 的谱模式仍然是各预报中心的全球中期预报业务模式的主流。 分布式并行计算的实质是，利用多个处理机通过数据通信协调完成单一计算任务。在并行计算 机上实现并行计算的关键是将计算任务划分为多个独立的子任务。半 Langrange 方法由于需要计算 各网格点在上一个时间步的起始点和时间步中间点处的各种量值，因此半 Langrange 谱模式相对于 基于 Euler 描述的谱模式在并行计算实现方面有更大的困难性。ECMWF 的 Dent 等利用区域分裂的 思想研究半隐式半 Langrange 谱模式的并行计算，并首先针对共享主存的 CRAY Y-MP 实现了多任 务并行的 T213L31 业务系统。之后 Dent 等通过引入虚拟边界区 Halo 的概念研究半 Langrange 谱模 式在分布并行计算上的实现方法，并于 1996 年在 Fujitsu 的 VPP700 向量并行计算机上实现了 T213L31 分布并行版本的业务化。由于一般的 Halo 并行方法的可扩展性差，半隐式半 Langrange 方 法的可扩展并行算法是一个有待进一步深入研究的问题。 2 二个时间层的半 Langrange 时间积分格式 描述大学运动的模式控制方程可写成如下一般的微分方程形式： dX R L +N dt 其中L 为R 的线性部分，N 为R 的非线性部分。二个时间层的半 Langrange 方法的关键思想是根据 t 和t ??t 时刻的已知风速外推出t +?t 2 时刻的