风暴潮数值模式的并行化.doc

风暴潮数值模式的并行化 ComputerEngineeringandApplications计算机工程与应用 风暴潮数值模式的并行化 苗春葆,赵鹏,沈飙,刘永玲 MIAOChunbao,ZHAOPeng,SHENBiao,LIUYongling 中国海洋大学海洋环境学院,山东青岛266100 CollegeofPhysicatandEnvironmentalOceanography,OceanUniversityofChina,Qingdao,Shandong266100,China MIAOChunbao,ZHAOPeng,SHENBiao,eta1.Parallelizationofnumericalstormsurgemode1.ComputerEngineeringand Applications,2012,48(2):39-42. Abstract:AnumericalstormsurgemodelisparallelizedusingMessagePassingInterface(MPI).Anewparallelalgorithmforsolving traditionallinearequationswithhighefficiencyandsimplicityisintroduced.Loadbalanceshouldbeprimarilyconsideredduringopti- mizationofaparallelprogram.Thecomputationloadisperfectlybalancedbyreasonabledecompositionofsimulatedregionaccording totheamountofwatergrids,andresultsinasignificantperformanceenhancement.TheaccelerationrateonSMPplatformsusing8 CPUsCanreachto7.0,and6.5onclusterplatforms. Keywords:parallelization;stormsurge;tridiagonallinearequation;numericalmodel ● 摘要:实现了风暴潮数值模式基于MPI的并行化;根据该模式数值计算的特点提出了一种并行求解三对角方程组的新方法,相 对于传统算法编程简单而且并行效率更高;负载平衡是并行程序性能优化首先要解决的问题,以水格点的个数作为任务分解的 标准,实现了较好的负载平衡,相比水陆格点不作区分的分解方法性能有明显的提高;在SMP平台上使用8个CPU时加速比可以 达到70,在集群平台上为6.5. 关键词:并行;风暴潮;三对角方程组;数值模式 DOI:10.3778~.issn.1002.8331.2012.02.012文章编号:1002.833l(2012)02.0039.04文献标识码:A中图分类号:TP3l1.11 l引言 风暴潮是指由于强烈的大气扰动所导致的海面异常升高 的现象.如果风暴潮发生时I合遇上天文大潮的高潮阶段,则 往往会使其影响所及的海域水位暴涨,给沿海人民生命和财 产带来重大损失.因此,风暴潮是海洋工程中必须考虑的一 个重要因素. 2O世纪80年代以来,随着计算机技术的发展,风暴潮数值 模拟及预报也取得了较大的进展.近些年来,风暴潮也受到 了广大海洋研究者的高度重视u,是当前的研究热点之一. 从数值模拟的角度来说,网格分辨率越高,时间步长越小,数 值模拟的精度亦会越高.此外,风暴潮过程受地形影响很大, 只有高分辨率的计算网格才能实现风暴潮的精细预报;但增 加网格分辨率时,计算量也会显着增大.而风暴潮必须在尽 量短的时间内做出预报,因此如何在当前的计算条件下快速 准确地对风暴潮做出预报成为各地海洋预报部门关注的问 题.高性能计算机的应用及并行计算的引入,使更高网格分 辨率的风暴潮数值模拟及陕速预报成为现实.本文正是从这 一 点出发,研究了风暴潮数值模式的并行算法,并基于MPI实 现了其并行化. 2并行算法设计 2.1数值模型 模型的控制方程采用以下二维深度平均流模型: 筹+[++_-- ++=一g 善一鲁+ 百OV++=筹一+ 其中,为相对于平均水面的水位高度,即风暴潮增水水位,h 为静水水深,厂为柯氏参数,g为重力加速度,和分别为 X和Y方向的深度平均流,P为海平面气压,P=1kedm为海 水密度., .,. 和r, , , 分别为x和Y方向的海面风应 力,海底摩擦力. 风应力:Ta=CDP,其中:为风速,c_D是拖曳系数, 由经验公式来给定.虽然采用椭圆形风场模型的模拟效果 会更好一些,但为简便起见,这里仍采用圆形风场模型.海 底摩擦力:”ga