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 ofmore »
- Authors:
-
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
- Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
- Publication Date:
- Research Org.:
- Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
- Sponsoring Org.:
- USDOE
- OSTI Identifier:
- 1201547
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Procedia Computer Science
- Additional Journal Information:
- Journal Volume: 51; Journal Issue: C; Journal ID: ISSN 1877-0509
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; 54 ENVIRONMENTAL SCIENCES; matrix-vector multiply; spectral element solvers; Newton's method
Citation Formats
Woodward, Carol S., Gardner, David J., and Evans, Katherine J. On the use of finite difference matrix-vector products in Newton-Krylov solvers for implicit climate dynamics with spectral elements. United States: N. p., 2015.
Web. doi:10.1016/j.procs.2015.05.468.
Woodward, Carol S., Gardner, David J., & Evans, Katherine J. On the use of finite difference matrix-vector products in Newton-Krylov solvers for implicit climate dynamics with spectral elements. United States. https://doi.org/10.1016/j.procs.2015.05.468
Woodward, Carol S., Gardner, David J., and Evans, Katherine J. Thu .
"On the use of finite difference matrix-vector products in Newton-Krylov solvers for implicit climate dynamics with spectral elements". United States. https://doi.org/10.1016/j.procs.2015.05.468. https://www.osti.gov/servlets/purl/1201547.
@article{osti_1201547,
title = {On the use of finite difference matrix-vector products in Newton-Krylov solvers for implicit climate dynamics with spectral elements},
author = {Woodward, Carol S. and Gardner, David J. and Evans, Katherine J.},
abstractNote = {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.},
doi = {10.1016/j.procs.2015.05.468},
journal = {Procedia Computer Science},
number = C,
volume = 51,
place = {United States},
year = {Thu Jan 01 00:00:00 EST 2015},
month = {Thu Jan 01 00:00:00 EST 2015}
}
Web of Science
Works referenced in this record:
Hybrid Schwarz-Multigrid Methods for the Spectral Element Method: Extensions to Navier-Stokes
book, January 2005
- Fischer, Paul F.; Lottes, James W.
- Lecture Notes in Computational Science and Engineering
An initial-value problem for testing numerical models of the global shallow-water equations
journal, January 2004
- Galewsky, Joseph; Scott, Richard K.; Polvani, Lorenzo M.
- Tellus A: Dynamic Meteorology and Oceanography, Vol. 56, Issue 5
On finite difference approximation of a matrix-vector product in the Jacobian-free Newton–Krylov method
journal, October 2011
- An, Heng-Bin; Wen, Ju; Feng, Tao
- Journal of Computational and Applied Mathematics, Vol. 236, Issue 6
A Local Convergence Theory for Combined Inexact-Newton/Finite-Difference Projection Methods
journal, April 1987
- Brown, Peter N.
- SIAM Journal on Numerical Analysis, Vol. 24, Issue 2
Hybrid Krylov Methods for Nonlinear Systems of Equations
journal, May 1990
- Brown, Peter N.; Saad, Youcef
- SIAM Journal on Scientific and Statistical Computing, Vol. 11, Issue 3
On Using Approximate Finite Differences in Matrix-Free Newton–Krylov Methods
journal, January 2008
- Brown, Peter N.; Walker, Homer F.; Wasyk, Rebecca
- SIAM Journal on Numerical Analysis, Vol. 46, Issue 4
AMIP Simulation with the CAM4 Spectral Element Dynamical Core
journal, February 2013
- Evans, K. J.; Lauritzen, P. H.; Mishra, S. K.
- Journal of Climate, Vol. 26, Issue 3
A Scalable and Adaptable Solution Framework within Components of the Community Climate System Model
book, January 2009
- Evans, Katherine J.; Rouson, Damian W. I.; Salinger, Andrew G.
- Lecture Notes in Computer Science
Accuracy Analysis of a Spectral Element Atmospheric Model Using a Fully Implicit Solution Framework
journal, August 2010
- Evans, Katherine J.; Taylor, Mark A.; Drake, John B.
- Monthly Weather Review, Vol. 138, Issue 8
Hybrid Schwarz-Multigrid Methods for the Spectral Element Method: Extensions to Navier-Stokes
book, January 2005
- Fischer, Paul F.; Lottes, James W.
- Lecture Notes in Computational Science and Engineering
An Overlapping Schwarz Method for Spectral Element Solution of the Incompressible Navier–Stokes Equations
journal, May 1997
- Fischer, Paul F.
- Journal of Computational Physics, Vol. 133, Issue 1
An initial-value problem for testing numerical models of the global shallow-water equations
journal, January 2004
- Galewsky, Joseph; Scott, Richard K.; Polvani, Lorenzo M.
- Tellus A: Dynamic Meteorology and Oceanography, Vol. 56, Issue 5
A Spectral Deferred Correction Method Applied to the Shallow Water Equations on a Sphere
journal, October 2013
- Jia, Jun; Hill, Judith C.; Evans, Katherine J.
- Monthly Weather Review, Vol. 141, Issue 10
Jacobian-free Newton–Krylov methods: a survey of approaches and applications
journal, January 2004
- Knoll, D. A.; Keyes, D. E.
- Journal of Computational Physics, Vol. 193, Issue 2
Algorithmically scalable block preconditioner for fully implicit shallow-water equations in CAM-SE
journal, October 2014
- Lott, P. Aaron; Woodward, Carol S.; Evans, Katherine J.
- Computational Geosciences, Vol. 19, Issue 1
A Flexible Inner-Outer Preconditioned GMRES Algorithm
journal, March 1993
- Saad, Youcef
- SIAM Journal on Scientific Computing, Vol. 14, Issue 2
Conservative Finite-Difference Approximations of the Primitive Equations on Quasi-Uniform Spherical Grids
journal, February 1972
- Sadourny, Robert
- Monthly Weather Review, Vol. 100, Issue 2
A compatible and conservative spectral element method on unstructured grids
journal, August 2010
- Taylor, Mark A.; Fournier, Aimé
- Journal of Computational Physics, Vol. 229, Issue 17
A standard test set for numerical approximations to the shallow water equations in spherical geometry
journal, September 1992
- Williamson, David L.; Drake, John B.; Hack, James J.
- Journal of Computational Physics, Vol. 102, Issue 1