On the Use of Finite Difference Matrix-vector Products in Newton-krylov Solvers for Implicit Climate Dynamics with Spectral Elements

Abstract Efficient solutions of global climate models require effectively handling disparate length and time scales. Implicit solution approaches allow time integration of the physical system with a step size governed by accuracy of the processes of interest rather than by stability of the fastest time scales present. Implicit approaches, however, require the solution of nonlinear systems within each time step. Usually, a Newton's method is applied to solve these systems. Each iteration of the Newton's method, in turn, requires the solution of a linear model of the nonlinear system. This model employs the Jacobian of the problem-defining nonlinear residual, but this Jacobian can be costly to form. If a Krylov linear solver is used for the solution of the linear system, the action of the Jacobian matrix on a given vector is required. In the case of spectral element methods, the Jacobian is not calculated but only implemented through matrix- vector products. The matrix-vector multiply can also be approximated by a finite difference approximation which may introduce inaccuracy in the overall nonlinear solver. In this paper, we review the advantages and disadvantages of finite difference approximations of these matrix- vector products for climate dynamics within the spectral element shallow water dynamical core of the Community Atmosphere Model (CAM).

[1]  Katherine J. Evans,et al.  Recent Progress in Nonlinear and Linear Solvers , 2011 .

[2]  Carol S. Woodward,et al.  User Documentation for KINSOL v2.2.0 , 2004 .

[3]  Paul Fischer,et al.  An Overlapping Schwarz Method for Spectral Element Solution of the Incompressible Navier-Stokes Equations , 1997 .

[4]  R. K. Scott,et al.  An initial-value problem for testing numerical models of the global shallow-water equations , 2004 .

[5]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

[6]  Mark A. Taylor,et al.  A Scalable and Adaptable Solution Framework within Components of the Community Climate System Model , 2009, ICCS.

[7]  P. Swarztrauber,et al.  A standard test set for numerical approximations to the shallow water equations in spherical geometry , 1992 .

[8]  T. A. Zang,et al.  Spectral Methods: Fundamentals in Single Domains , 2010 .

[9]  Yousef Saad,et al.  Hybrid Krylov Methods for Nonlinear Systems of Equations , 1990, SIAM J. Sci. Comput..

[10]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[11]  Homer F. Walker,et al.  On Using Approximate Finite Differences in Matrix-Free Newton-Krylov Methods , 2008, SIAM J. Numer. Anal..

[12]  Katherine J. Evans,et al.  AMIP Simulation with the CAM4 Spectral Element Dynamical Core , 2013 .

[13]  M. Taylor,et al.  Accuracy Analysis of a Spectral Element Atmospheric Model Using a Fully Implicit Solution Framework , 2010 .

[14]  Tao Feng,et al.  On finite difference approximation of a matrix-vector product in the Jacobian-free Newton-Krylov method , 2011, J. Comput. Appl. Math..

[15]  Yousef Saad,et al.  A Flexible Inner-Outer Preconditioned GMRES Algorithm , 1993, SIAM J. Sci. Comput..

[16]  Paul Fischer,et al.  Hybrid Schwarz-Multigrid Methods for the Spectral Element Method: Extensions to Navier-Stokes , 2005 .

[17]  R. Sadourny.  Conservative Finite-Difference Approximations of the Primitive Equations on Quasi-Uniform Spherical Grids , 1972 .

[18]  Paul A. Ullrich,et al.  The spectral element method (SEM) on variable-resolution grids: evaluating grid sensitivity and resolution-aware numerical viscosity , 2014 .

[19]  Katherine J. Evans,et al.  A Spectral Deferred Correction Method Applied to the Shallow Water Equations on a Sphere , 2013 .

[20]  Katherine J. Evans,et al.  Algorithmically scalable block preconditioner for fully implicit shallow-water equations in CAM-SE , 2014, Computational Geosciences.

[21]  P. Brown.  A local convergence theory for combined inexact-Newton/finite-difference projection methods , 1987 .

[22]  Mark A. Taylor,et al.  A compatible and conservative spectral element method on unstructured grids , 2010, J. Comput. Phys..