This paper focuses on the development of an efficient,three-dimensional,thermo-mechanical,nonlinear-Stokes flow computational model for ice sheet simulation.The model is based on the parallel finite element model deve...This paper focuses on the development of an efficient,three-dimensional,thermo-mechanical,nonlinear-Stokes flow computational model for ice sheet simulation.The model is based on the parallel finite element model developed in[14]which features high-order accurate finite element discretizations on variable resolution grids.Here,we add an improved iterative solution method for treating the nonlinearity of the Stokes problem,a new high-order accurate finite element solver for the temperature equation,and a new conservative finite volume solver for handling mass conservation.The result is an accurate and efficient numerical model for thermo-mechanical glacier and ice-sheet simulations.We demonstrate the improved efficiency of the Stokes solver using the ISMIP-HOM Benchmark experiments and a realistic test case for the Greenland ice-sheet.We also apply our model to the EISMINT-II benchmark experiments and demonstrate stable thermo-mechanical ice sheet evolution on both structured and unstructured meshes.Notably,we find no evidence for the“cold spoke”instabilities observed for these same experiments when using finite difference,shallow-ice approximation models on structured grids.展开更多
基金the U.S.Department of Energy,Office of Science,Advanced Scientific Computing Research and Biological and Environmental Research programs through the Scientific Discovery through Advanced Computing(SciDAC)project PISCEES,and by the US National Science Foundation under the grant number DMS-1215659'the National 863 Project of China under the grant number 2012AA01A309'the National Center for Mathematics and Interdisciplinary Sciences of the Chinese Academy of Sciences.
文摘This paper focuses on the development of an efficient,three-dimensional,thermo-mechanical,nonlinear-Stokes flow computational model for ice sheet simulation.The model is based on the parallel finite element model developed in[14]which features high-order accurate finite element discretizations on variable resolution grids.Here,we add an improved iterative solution method for treating the nonlinearity of the Stokes problem,a new high-order accurate finite element solver for the temperature equation,and a new conservative finite volume solver for handling mass conservation.The result is an accurate and efficient numerical model for thermo-mechanical glacier and ice-sheet simulations.We demonstrate the improved efficiency of the Stokes solver using the ISMIP-HOM Benchmark experiments and a realistic test case for the Greenland ice-sheet.We also apply our model to the EISMINT-II benchmark experiments and demonstrate stable thermo-mechanical ice sheet evolution on both structured and unstructured meshes.Notably,we find no evidence for the“cold spoke”instabilities observed for these same experiments when using finite difference,shallow-ice approximation models on structured grids.
