收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

GRAPES高分辨率气象数值预报模式并行计算关键技术研究

伍湘君  
【摘要】:建立高分辨率精细数值天气预报模式是大气科学研究和数值天气预报模式业务发展的主流方向。随着高性能计算机计算能力的提高,高性能计算机规模不断扩大。面对峰值性能达到每秒千万亿次的高性能计算机系统,数值天气预报模式能否充分利用计算平台所提供的计算能力,解决其庞杂的科学计算问题和海量数据处理问题,取决于数值模式的并行计算方案和并行实现方法,特别是其核心算法的计算效率。 以我国自主研发的新一代研究/业务数值天气预报系统GRAPES(Global/Regional Assimilation and PrEdiction System)的模式系统为基础,本文深入分析了GRAPES模式的科学计算原理,探讨了影响GRAPES模式并行计算效率的主要因素,设计了GRAPES分层软件框架结构和并行计算方案,建立了GRAPES模式的并行计算系统,并针对其中影响并行计算效率的关键问题提出了优化方案,重点研究了针对拉格朗日插值计算的并行实现方法,研究了GRAPES模式关键计算Helmholtz方程的高效求解算法,最后在天河-1A超级计算机上对GRAPES模式并行计算系统进行了一系列的测试分析。本文主要研究成果包括: 1、深入分析数值天气预报模式科学计算原理,揭示数据并行是数值天气预报模式并行计算的基本策略。由于模式数据流计算是随时间依次完成的,不同时间积分步之间数据流前后依赖,数值模式并行计算只能在一个积分步内进行。而同样由于数据相关性,数值模式并行计算通常采用数据并行方式,即采用将预报区域按照计算核数划分成块进行计算。 2、针对日益复杂的数值模式系统,分析指出数值天气预报软件系统必须采用软件工程方法进行组织、管理,软件开发必须遵循软件规范要求,并结合高性能计算机体系结构特点,设计了GRAPES分层软件框架结构,建立了符合软件工程规范的并行编程接口函数库PPI,实现了GRAPES模式并行版本基本软件架构。 3、针对拉格朗日插值计算,分析了影响GRAPES全球模式并行计算实现的极地区域网格变量聚集问题,提出了以―供方‖为中心的并行计算方案(put-scheme),实现了以―需方‖为中心的并行计算方案(get-scheme),改进了任务分配算法。测试结果表明,两种方案均有效解决了极地网格聚集对拉格朗日上游点插值并行计算的影响。但从计算效率而言,―get-scheme‖方案更具优势:1)减少了极地区域的大内存需求;2)减少了极地区域数据通讯的盲目性;3)增加了低纬度地区上游点位移量的允许范围。因此,―get-scheme‖方案具有更好负载平衡性和并行可扩展性能。 4、针对GRAPES模式中占主要计算开销的Helmholtz方程的求解,实现了基于科学计算可移植扩展工具包(Portable Extensible Toolkit for Scientific Computation,即PETSc)和高层并行预条件函数库(high performance preconditions,即Hypre)的广义极小残量法(GMRES)求解算法。与目前GRAPES模式版本中使用的广义共轭余差法(GCR)方法相比,GMRES方法具有收敛速度快、迭代次数少、求解精度高、并行可扩展性能好等特点。对于高分辨率精细模式,采用GRMES方法求解Helmholtz方程大大减少了GRAPES模式的计算开销,显著提高了GRAPES模式在大规模并行计算机上的运行效率。 5、通过对不同收敛精度Helmholtz方程求解的理想试验、实际资料绝热模式以及全物理过程多方面测试,揭示了GRAPES模式动力框架计算精度被物理过程计算精度所掩盖的问题。模式计算精度是整个动力框架计算和物理过程计算的综合结果,积分计算中每个过程的计算偏差都会在一定程度上反映到预报结果的偏差上,因此提高模式计算精度必须从模式计算的多个方面入手。 6、针对GRAPES核心算法,建立了IBM-cluster1600计算机上拉格朗日插值并行计算两种方法的并行通讯分析模型,以及Helmholtz方程两种求解方法的并行计算时间模型;通过IBM-cluster1600计算机上固定规模GRAPES模式可扩展性能测试,验证了GRAPES模式并行计算系统具有良好的并行可扩展性能。GRAPES模式在天河-1A超级计算机上的测试分析表明:1)GRAPES模式积分计算部分并行效率较高,并行计算保持了高可扩展性能,10天预报在2048个计算核上的计算效率接近90%;2)目前影响GRAPES模式并行可扩展性能的瓶颈一个是I/O操作,另一个就是如何将GRAPES模式更好应用于分层设计的计算机体系结构。 通过本文的研究,实现了具有良好可扩展性能的GRAPES模式并行计算系统。目前GRAPES并行模式系统已在国家气象中心业务运行(GRAPES区域模式业务运行,全球模式准业务运行),计算正确稳定,满足实时性业务要求。GRAPES模式并行计算系统的建立为GRAPES资料同化系统并行积累了经验,为GRAPES数值预报系统发展奠定了基础。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 王福军,程建钢,姚振汉;结构非线性动力分析显式积分并行算法[J];清华大学学报(自然科学版);2002年04期
2 冯诗齐;用并行计算构建一个虚拟太空[J];世界科学;2002年10期
3 苑明哲,于海斌,周悦;流水线型并行控制器[J];控制理论与应用;2004年04期
4 于方;;基于LINUX平台的PC机群并行环境的构建与应用[J];阴山学刊(自然科学版);2007年04期
5 刘巧英;乔玉兰;孙玉强;;并行技术领域中调度算法研究[J];中国科技信息;2009年05期
6 徐晓飞;曹祥玉;姚旭;陈盼;;一种基于Doolittle LU分解的线性方程组并行求解方法[J];电子与信息学报;2010年08期
7 费景高;并行显式Runge-Kutta公式的实现[J];计算机工程与设计;1994年05期
8 李强,邹经湘,黄文虎;结构动响应的并行算法[J];工程力学;1999年03期
9 胖质玲,赵芳;电力系统并行计算的实现[J];山东科技大学学报(自然科学版);1999年04期
10 朱小谦,张卫民,宋君强;天气预报有限差分格点模式并行计算研究[J];计算机工程与科学;2004年04期
11 何素贞,李书平,吴晨旭;高性能计算集群系统的设计和实现[J];厦门大学学报(自然科学版);2004年06期
12 王开健,刘西拉,顾雷;基于MPI机群环境下的广义逆力法并行化初探[J];岩石力学与工程学报;2005年01期
13 刘胜利,张晔,许世菊,马富明;四阶抛物方程的一个并行有限差分格式[J];吉林大学学报(理学版);2005年06期
14 王丹;陈硕;尚智;易杰;范西俊;;方形通道内流体流动的耗散粒子动力学并行计算[J];上海交通大学学报;2007年06期
15 刘玉雪;李文敬;杨立志;张海志;;Credit-Recovery分布式并行终止检测算法的改进[J];广西师范学院学报(自然科学版);2009年02期
16 李文强;刘晓;;追赶法并行求解循环三对角方程组[J];科技导报;2009年18期
17 郭永辉;王宏亮;田宙;;堆积颗粒消波效应模型及其并行计算[J];计算机工程与科学;2009年S1期
18 蒋伯诚,成克懋;快速傅里叶变换的并行计算[J];计算物理;1989年01期
19 陆益君;;超立方机上一阶线性递推问题的并行求解(英文)[J];应用数学;1991年04期
20 张学惠,韩永飞;计算并行性的研究和描述[J];吉林大学学报(理学版);1992年03期
中国重要会议论文全文数据库 前10条
1 范晓樯;李桦;田正雨;;超声速/高超声速飞行器复杂流场大规模并行数值仿真[A];计算流体力学研究进展——第十二届全国计算流体力学会议论文集[C];2004年
2 齐进;叶文华;;三维激光烧蚀瑞利-泰勒不稳定性并行计算[A];中国空气动力学学会第十届物理气体动力学专业委员会会议论文集[C];2001年
3 张望;王辉;;个性化服务中的并行K-Means聚类算法[A];2007年全国开放式分布与并行计算机学术会议论文集(下册)[C];2007年
4 丛鹏;;MPI并行计算实现工业CT图像重建[A];2004年CT和三维成像学术年会论文集[C];2004年
5 丁国昊;罗凯;李伟;李桦;;乘波飞行器气动特性数值模拟与并行计算[A];第三届高超声速科技学术会议会议文集[C];2010年
6 罗文彩;陈小前;;并行计算的多方法优化协作[A];第二十四届中国控制会议论文集(上册)[C];2005年
7 耿江东;薛正辉;高本庆;;应用并行GTD算法计算阵列天线近场受扰[A];第17届全国电磁兼容学术会议论文集[C];2007年
8 唐维军;张景琳;蔚喜军;;三维流体界面不稳定性的并行计算[A];中国工程物理研究院科技年报(2000)[C];2000年
9 左风丽;莫则尧;叶文华;;计算流体三维分裂格式的高效并行计算[A];中国工程物理研究院科技年报(2003)[C];2003年
10 杜志文;曾文华;;网格计算在文本分类中的应用[A];2006年全国开放式分布与并行计算机学术会议论文集(三)[C];2006年
中国博士学位论文全文数据库 前10条
1 陈军;分布式存储环境下并行计算可扩展性的研究与应用[D];中国人民解放军国防科学技术大学;2000年
2 尹欣;三维弹性问题边界元法并行计算及其工程应用[D];清华大学;2000年
3 陈晓春;基于并行计算的大涡模拟方法及其工程应用基础研究[D];西安建筑科技大学;2004年
4 王开健;基于特大增量步算法的网络并行计算[D];清华大学;2005年
5 张理论;面向气象预报数值模式的高效并行计算研究[D];中国人民解放军国防科学技术大学;2002年
6 寇哲君;可扩展冲击—接触并行计算及其在汽车碰撞模拟中的应用[D];清华大学;2003年
7 刘丽;人工免疫网络研究及应用[D];江南大学;2008年
8 蔡勇;基于GPU的车身结构接触碰撞过程并行计算方法[D];湖南大学;2013年
9 张永彬;岩石破裂过程分析并行计算方法研究[D];东北大学;2007年
10 姚路;图分割算法及其在大规模数值并行计算中的应用研究[D];国防科学技术大学;2013年
中国硕士学位论文全文数据库 前10条
1 刘战合;基于复合形方法的翼型优化设计及并行计算研究[D];西北工业大学;2004年
2 金光浩;有向图并行计算中的多目标剖分算法[D];中国工程物理研究院;2005年
3 岳常智;135四气门直喷式柴油机工作过程数值计算[D];大连理工大学;2006年
4 曹大勇;O_2在MgO(100)表面吸附的第一原理分子动力学研究[D];大连理工大学;2004年
5 叶敏娇;Manticore体系结构设计[D];浙江大学;2006年
6 赵玉林;H.26L视频编码器并行性算法研究与实现[D];西南交通大学;2003年
7 吴颖;一种基于LINUX的集群计算系统的设计与实现[D];四川大学;2005年
8 彭虹;分形视频图像压缩并行算法设计[D];福州大学;2006年
9 万婷婷;并行计算的网格应用研究[D];广东工业大学;2006年
10 朱丽莎;基于GPU的一维热传导算法研究[D];暨南大学;2011年
中国重要报纸全文数据库 前10条
1 轶嘉;英特尔全球首个并行计算中心落户无锡[N];人民邮电;2009年
2 江锡民;英特尔并行计算中心落户无锡[N];新华日报;2009年
3 刘琦;伯克利专家展望未来并行计算[N];中国计算机报;2008年
4 均儿;通用计算核动力[N];电脑报;2009年
5 本报记者 马文方;英特尔为何要牵头并行计算[N];中国计算机报;2009年
6 英特尔并行计算实验室研究员 TimothyMattson;并行计算:减少串行软件[N];中国计算机报;2007年
7 英特尔 赵军(Jun Zhao);PC机并行计算革命尚未成功[N];中国计算机报;2009年
8 ;并行计算成PC产业发展瓶颈[N];人民邮电;2008年
9 刘霞;计算能力的提升需要一场革命[N];科技日报;2010年
10 张云泉;并行计算:迎接多核时代的挑战[N];计算机世界;2006年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978