基于JASMIN的地下水流大规模并行数值模拟.pdfVIP

第30卷第3期 计 算 物 理 V01．30，No．3 CHINESE OFCOMPUTATIONALPHYSICS 2013年5月 JOURNAL May，2013 文章编号：1001—246X(2013J03-0317-09 基于JASMIN的地下水流大规模并行数值模拟 程汤培1’2， 莫则尧1’2， 邵景力3 (1．北京应用物理与计算数学研究所计算物理重点实验室，北京100088； 2．北京应用物理与计算数学研究所高性能计算中心，北京 100094； 3．中国地质大学(北京)水资源与环境学院，北京100083) 摘要：针对具有精细网格剖分、长时间跨度特征的地下水流模拟中计算时间长、存储开销大等瓶颈问题，基于 MODFLOW三维非稳定流计算方法，提出基于网格片的核心算法以及基于影像区的通信机制，并在JASMIN框架上 研制了大规模地下水流并行数值模拟程序JOGFLOW．通过河南郑州市中牟县雁鸣湖水源地地下水流的模拟，对程 序正确性和性能进行了验证；通过建立一个具有精细网格剖分的假想地下水概念模型对可扩展性进行测试．相对 于32核的并行程序，在512以及1024个处理机上的并行效率分别可达77．2％和67．5％．数值模拟结果表明， JOGFLOW具有较好的计算性能与可扩展性，能够有效使用数百上千计算核心，支持千万量级以上网格剖分的地下 水流模型的大规模并行计算． 关键词：高仿真地下水流模拟；大规模并行计算；JOGFLOW；MODFLOW；JASMIN框架 中图分类号：T301．6，P641．2 文献标识码：A 0 研究背景 随着地下水演变机制研究以及地下水资源调控的不断深入，对精细网格剖分、长时间跨度特征的高仿真 地下水模型的数值模拟提出了迫切需求．然而，由于传统串行数值模拟软件在处理器、内存以及数据通道上 的限制，在建立水文地质模型时往往只能采用不精确的粗粒度网格去逼近．这种近似方法，极大地降低了模 拟的精确度和可信度． 针对此类问题，国外的一些研究团体提出了并行计算方法并研制了相关的并行计算程序，如TOUGH2一 置、地下水修复、多孔以及裂隙介质中多向流、多组分运移以及热传导的模拟等领域中．国内的研究工作主 要集中于区域分解算法”9I、线性解法器并行化¨¨141、渗流场的有限元并行计算¨5。171以及水文地质参数并 行反求与优化¨“2叫几个方面．MODFLOW作为世界上求解地下水流相关问题使用最广泛的数值模拟软件 提高问题求解规模方面取得了一定的进展，但较少考虑负载均衡等制约并行可扩展性的关键因素，难于有效 利用更多的处理器进行求解．而面向MODFLOW地下水流模拟在数百上千计算核心上的大规模高效并行工 作，国内外研究中尚未见报导． 针对MODFLOW中所涉及的地下水流动数值模拟过程，本文基于并行自适应结构网格应用支撑框架 收稿日期：2012—08—24；修回日期：2012—11—27 发展计划(863)(2010AA012301)资助项目 万方数据 计 算 物 理 第30卷 structuredmeshes JASMIN(Jadaptive applications 行了代码重构，提出了相应的并行实现策略以及负载均衡方法，编制了面向大规模地下水流模型的高效并行 OrientedGroundwaterFLOW 计算程序JOGFLOW(JASMIN program)．针对两个典型数值实例，对程序的正确 性以及性能可扩展性进行了验证．结果表明：JOGFLOW显著地加快地下水流模型计算速度、提高问题求解 规模，并具有良好的并行可扩展性． 1 数学模型 不考虑水的密度变化，孔隙介质中地下水在三维空间的流动可用下面的偏微分方程描述心3|：