ＧＰＵによって加速化された氷床モデリング用プログラムの開発と応用

GPU加速冰盖建模程序开发与应用

• 批准号: 15F15902
• 负责人: Greve Ralf
• 金额: $ 1.47万
• 依托单位: Hokkaido University
• 依托单位国家: 日本
• 项目类别: Grant-in-Aid for JSPS Fellows
• 财政年份: 2015
• 资助国家: 日本
• 起止时间: 2015-04-24 至 2017-03-31
• 项目状态: 已结题
• 来源: https://kaken.nii.ac.jp/grant/KAKENHI-PROJECT-15F15902/
• 关键词: Ice sheet; Dynamics; Full Stokes; Numerical model; GPU computing; Finite element method; Matrix assembly

Computer simulations have become established as an important technique for the understanding of ice sheet dynamics. In particular, they are relevant for predicting the evolution of ice sheets in a future changing climate. The goal of this project was to develop an operational ice sheet model that employs the technique of GPU- (graphics processing unit) accelerated computing.This model, termed "Sainou", is an Elmer/Ice (elmerice.elmerfem.org) derivative that solves the full-Stokes equations (no stresses neglected) of ice dynamics with the finite element method. In the second project year (FY 2016), the focus was on testing Sainou, using an established benchmark problem, namely the ISMIP-HOM experiment A (Pattyn et al., 2008, Cryosphere 2, 95-108). This experiment consists of an inclined slab of ice flowing over a bumpy bed.The computing time required by Sainou for the assembly of the finite element matrix of the full-Stokes system of partial differential equations was compared to that required by Elmer/Ice. Four different numbers of elements were tested for the ISMIP-HOM experiment A, and the simulations were carried out on a Mac Pro workstation, using only a single CPU. It was found that, for double-precision accuracy with 8000, 87500, 252000 and 391500 brick elements, respectively, Sainou is 7x, 17x, 23x and 25x faster than Elmer/Ice, and for single-precision accuracy it is even 9x, 38x, 45x and 49x faster. It could therefore be demonstrated successfully that GPU-accelerated computing is a useful method for fast and efficient ice sheet modelling.

计算机模拟已经成为理解冰盖动力学的一项重要技术。特别是，它们与预测未来气候变化中冰盖的演变有关。该项目的目标是开发一个可操作的冰盖模型，该模型采用GPU（图形处理单元）加速计算技术。该模型被称为“Sainou”，是Elmer/Ice（elmerice.elmerfem.org）的衍生物，该衍生物用有限元法求解冰动力学的全斯托克斯方程（不忽略应力）。在第二个项目年（2016财年），重点是使用已建立的基准问题，即ISMIP-HOM实验A（Pattyn等人，2008，Cryosphere 2，95-108）。本实验由一个倾斜的平板冰流过一个凹凸不平的床。所需的计算时间进行了比较，Sainou的全斯托克斯系统的偏微分方程的有限元矩阵的组装所需的Elmer/冰。ISMIP-HOM实验A测试了四种不同数量的元素，并且模拟在Mac Pro工作站上进行，仅使用单个CPU。结果发现，对于8000，87500，252000和391500砖元素的双精度精度，Sainou分别比Elmer/Ice快7倍，17倍，23倍和25倍，对于单精度，它甚至快9倍，38倍，45倍和49倍。因此，可以成功地证明，GPU加速计算是一个有用的方法，快速和有效的冰盖建模。

项目成果

Towards a GPU-accelerated ice sheet model: design philosophy and implementation

走向 GPU 加速的冰盖模型：设计理念和实现

• 发表时间: 2015
• 作者: Seddik;H. and R. Greve
• 通讯作者: H. and R. Greve

Full Stokes finite-element modeling of ice sheets using a graphics processing unit

使用图形处理单元对冰盖进行全斯托克斯有限元建模

• 发表时间: 2016
• 作者: Seddik;H. and R. Greve
• 通讯作者: H. and R. Greve