Computer Science > Numerical Analysis
[Submitted on 25 Jun 2014 (this version), latest version 9 Jul 2015 (v3)]
Title:Solution of nonlinear Stokes equations discretized by high-order finite elements on nonconforming and anisotropic meshes, with application to ice sheet dynamics
View PDFAbstract:Motivated by the need for efficient and accurate simulation of the dynamics of the polar ice sheets, we design high-order finite element discretizations and scalable solvers for the solution of nonlinear incompressible Stokes equations. In particular, we focus on power-law, shear thinning rheologies commonly used in modelling ice dynamics and other geophysical flows. We use nonconforming hexahedral finite element meshes and the inf-sup stable finite element velocity-pressure pairings $\mathbb{Q}_k \times \mathbb{Q}_{k-2}^{\text{disc}}$ or $\mathbb{Q}_k \times \mathbb{P}_{k-1}^{\text{disc}}$, where $k \geq 2$ is the polynomial order of the velocity space. To solve the nonlinear equations, we propose a Newton-Krylov method with a block upper triangular preconditioner for the linearized Stokes systems. The diagonal blocks of this preconditioner are sparse approximations of the (1,1)-block and of its Schur complement. The (1,1)-block is approximated using linear finite elements based on the nodes of the high-order discretization, and is inverted using algebraic multigrid. This preconditioner is designed to be efficient on anisotropic meshes, which are necessary to match the high aspect ratio domains typical for ice sheets. Two different Schur complement approximations are compared, namely an inverse viscosity-weighted mass matrix and a least squares commutator approximation. In a comprehensive numerical study we find that our solver yields fast convergence that is independent of the mesh refinement, the occurrence of nonconforming element faces and the element aspect ratio, and depends only weakly on the polynomial finite element order. We simulate the ice flow in a realistic description of the Antarctic ice sheet derived from field data, and study the parallel scalability of our solver for problems with up to 615 million unknowns.
Submission history
From: Tobin Isaac [view email][v1] Wed, 25 Jun 2014 14:00:31 UTC (2,220 KB)
[v2] Fri, 13 Mar 2015 02:10:26 UTC (1,870 KB)
[v3] Thu, 9 Jul 2015 15:01:58 UTC (1,871 KB)
Current browse context:
math.NA
References & Citations
Bibliographic and Citation Tools
Bibliographic Explorer (What is the Explorer?)
Litmaps (What is Litmaps?)
scite Smart Citations (What are Smart Citations?)
Code, Data and Media Associated with this Article
CatalyzeX Code Finder for Papers (What is CatalyzeX?)
DagsHub (What is DagsHub?)
Gotit.pub (What is GotitPub?)
Papers with Code (What is Papers with Code?)
ScienceCast (What is ScienceCast?)
Demos
Recommenders and Search Tools
Influence Flower (What are Influence Flowers?)
Connected Papers (What is Connected Papers?)
CORE Recommender (What is CORE?)
arXivLabs: experimental projects with community collaborators
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.