Robust and efficient primal-dual Newton-Krylov solvers for viscous-plastic sea-ice models

by   Yu-hsuan Shih, et al.

We present a Newton-Krylov solver for a viscous-plastic sea-ice model. This constitutive relation is commonly used in climate models to describe the large scale sea-ice motion. Due to the strong nonlinearity of the momentum equation, the development of fast, robust and scalable solvers is still a substantial challenge. We propose a novel primal-dual Newton linearization for the momentum equation. In contrast to existing methods, it converges faster and more robustly with respect to mesh refinement, and thus allows fully resolved sea-ice simulations. Combined with an algebraic multigrid-preconditioned Krylov method for the Newton linearized systems, which contain strongly varying coefficients, the resulting solver scales well and can be used in parallel. We present highly resolved benchmark solutions and solve problems with up to 8.4 million spatial unknowns.


page 9

page 10

page 11


An adaptive solution strategy for Richards' equation

Flow in variably saturated porous media is typically modelled by the Ric...

A uniform preconditioner for a Newton algorithm for total-variation minimization and minimum-surface problems

Solution methods for the nonlinear partial differential equation of the ...

Parallel matrix-free higher-order finite element solvers for phase-field fracture problems

Phase-field fracture models lead to variational problems that can be wri...

Robust and Efficient Multilevel-ILU Preconditioned Newton-GMRES for Incompressible Navier-Stokes

We introduce a new preconditioned Newton-GMRES method for solving the no...

Large-scale Binary Quadratic Optimization Using Semidefinite Relaxation and Applications

In computer vision, many problems such as image segmentation, pixel labe...

A modified combined active-set Newton method for solving phase-field fracture into the monolithic limit

In this work, we examine a numerical phase-field fracture framework in w...

A Primal-Dual Solver for Large-Scale Tracking-by-Assignment

We propose a fast approximate solver for the combinatorial problem known...

Please sign up or login with your details

Forgot password? Click here to reset