椭圆型缓坡方程的微机群集并行数值模拟.pdfVIP

水 利 学 报 2003 年12 月 SHUILI XUEBAO 第12 期 文章编号：0559-9350 (2003) 12-0021-06 椭圆型缓坡方程的微机群集并行数值模拟 1 1 2 郑永红 ,游亚戈 ，沈永明 (1.中国科学院 广州能源研究所，广东 广州 510070；2.大连理工大学 海岸和近海工程国家重点实验室，辽宁 大连 116023) 摘要：提出了椭圆型缓坡方程的一种并行数值模拟方法，并在微机群集并行系统上进行了数值实现。为了验证方法 的精度和效率问题，对两种典型缓坡地形上的波浪传播变形进行了数值模拟。与对应的串型算法得到的结果比较表 明，本文的并行数值方法不仅可以提高数值模拟效率，而且可以大大减少单机的存储容量，从而使该系统可以进行 更大区域的波浪变形的数值模拟。 关键词：椭圆型缓坡方程；并行数值模拟；微机群集系统；水波 中图分类号：O353.2 文献标识码：A Berkhoff ［1］首次提出的椭圆型缓坡方程是研究缓坡地形上波浪变形的一种有效数学模型，由于其处理 复杂波浪场的通用性以及其结果和实验值能较好地符合，在海岸工程领域得到了较广泛的应用。但椭圆型 缓坡方程的数值求解存在着一些困难［2］ ［2-6］ ，使得其计算效率比较低下，因此一些学者 尝试了若干不同的 数值方法对其进行模拟，取得了比较令人满意的效果。但迄今为止，这些方法仍然是基于单机器的串行算 法，当计算域较大时(如10km×10km 区域)，这些方法的计算效率仍然比较低下，并且单机内存开销很大， 因此过去的这些串行算法是不太适合于模拟大区域的波浪变形问题的。 并行计算方法是提高计算效率的一种有效方法，其基本概念在20 世纪60 年代初就已出现，最初将这 种方法用于科学计算是在大型计算机或巨型计算机上进行的，自20 世纪90 年代中期，特别是WINDOWS 操 作系统出现后，并行计算方法开始在由普通微机组成的微机群集系统上得到实现，目前国内外在这方面的 研究正开展得方兴未艾。基于 PC 的廉价并行计算方案在达到超级计算机同样要求的情况下，具有十分显 著的优点：价格低廉(仅及超级计算机的几百分之一到几十分之一)；并行支撑软件大多为开放的免费软件， 如NT-MPICH；规模可大可小，可根据实际需要和资金确定；系统升级方便等。 就作者所知，过去应用椭圆型缓坡方程模拟线性波浪变形问题都是在单个计算机上采用串行计算方法 进行的，如文献 ［3-6］。为了进一步提高计算效率，克服单机内存不足的问题，本文在自建的微机群集并 行系统和原有串型算法的基础上，提出了椭圆型缓坡方程的一种并行数值模拟方法。为了验证算法的正确 性和效率问题，对两种典型地形上的波浪变形问题进行了数值模拟。由于对本文问题，没有其它的并行算 法得到的结果可以用于比较，因此这里将并行算法的精度和效率与其对应的串型算法的精度和效率进行对 比。从计算结果、并行效率和并行加速比等数据来看，本文提出的椭圆型缓坡方程的并行数值模拟方法是 比较成功的。 1 数学模型及其离散格式 收稿日期：2002-06-25 基金项目：国家十五863 计划项目(2001AA516010)；国家杰出青年科学基金项目；广东省自然科学基金项目 (021503) 作者简介：郑永红(1970-)，男，安徽潜山人，研究员，博士，从事近海海洋环境水力学和海洋能源方面的研究。 21 水 利 学 报 2003 年12 月 SHUILI XUEBAO 第12 期 Berkhoff ［1］提出的椭圆型缓坡模型的基本方程为：