Coupled Models and Parallel Simulations for Three-Dimensional Full-Stokes Ice Sheet Modeling
Coupled Models and Parallel Simulations for Three-Dimensional Full-Stokes Ice Sheet Modeling
复制标题
DOI：
10.4208/nmtma.2011.m1031
复制
发表时间：
2011
影响因子：
1.3
通讯作者：
Huai-xin Zhang;L. Ju
中科院分区：
文献类型：
--
作者：
Huai-xin Zhang;L. Ju
一个三维全斯托克斯计算模型被认为是确定的动力学，温度和厚度的冰盖。控制热力学方程由动量为非线性流变的三维全Stokes方程组、温度演化的对流扩散能量方程和冰厚变化的质量守恒方程组成。在这里，我们讨论的可变分辨率网格，有限元离散化，并采用模型组件的并行算法。求解器通过一个精心设计的耦合器进行集成，用于组件之间的参数数据交换。离散化利用高品质，可变分辨率的质心Voronoi Delaunay三角网格和现有的并行求解器。我们展示了网格技术，离散化方案，并行求解器的效率和可扩展性，通过计算实验，使用简化的几何基准测试问题和现实的格陵兰冰盖的几何形状。AMS科目分类：65M60、65M55、65Y05、65Z05、68U20、68W10
A three-dimensional full-Stokes computational model is considered for deter- mining the dynamics, temperature, and thickness of ice sheets. The governing thermo- mechanical equations consist of the three-dimensional full-Stokes system with nonlin- ear rheology for the momentum, an advective-diffusion energy equation for temper- ature evolution, and a mass conservation equation for ice-thickness changes. Here, we discuss the variable resolution meshes, the finite element discretizations, and the parallel algorithms employed by the model components. The solvers are integrated through a well-designed coupler for the exchange of parametric data between com- ponents. The discretization utilizes high-quality, variable-resolution centroidal Voronoi Delaunay triangulation meshing and existing parallel solvers. We demonstrate the grid- ding technology, discretization schemes, and the efficiency and scalability of the paral- lel solvers through computational experiments using both simplified geometries arising from benchmark test problems and a realistic Greenland ice sheet geometry. AMS subject classifications: 65M60, 65M55, 65Y05, 65Z05, 68U20, 68W10