A parallel Jacobian-free Newton-Krylov solver for a coupled sea ice-ocean model

Journal of Computational Physics

Volumn 257, Issue , 2014, Pages 901-911

  Losch, Martin ^a   Fuchs, Annika ^a   Lemieux, Jean François ^b   Vanselow, Anna ^c  

^a ALFRED WEGENER INSTITUTE FOR POLAR AND MARINE RESEARCH   (Germany)

^b Environnement et Changement Climatique   (Canada)

^c UNIVERSITY OF OLDENBURG   (Germany)

Author keywords

Jacobian free Newton Krylov solver; Numerical sea ice modeling; Parallel implementation; Preconditioning; Sea ice dynamics; Vector implementation

Indexed keywords

CLIMATE MODELS; ELASTICITY; NUMERICAL MODELS; PARALLEL ARCHITECTURES; PROGRAM PROCESSORS;

JACOBIAN-FREE NEWTON-KRYLOV; JACOBIAN-FREE NEWTON-KRYLOV SOLV; MOMENTUM EQUATION; NEWTON/KRYLOV SOLVERS; NUMERICAL SEA ICE MODELING; PARALLEL IMPLEMENTATIONS; PRECONDITIONING; SEA ICE DYNAMICS; SEA ICE MODELING; VECTOR IMPLEMENTATION;

SEA ICE;

EID: 84888022899     PISSN: 00219991     EISSN: 10902716     Source Type: Journal    
DOI: 10.1016/j.jcp.2013.09.026     Document Type: Article

References (39)

1. Proshutinsky A., Kowalik Z. Preface to special section on Arctic Ocean Model Intercomparison Project (AOMIP) studies and results. J. Geophys. Res. 2007, 112.
2. Rampal P., Weiss J., Dubois C., Campin J.-M. IPCC climate models do not capture Arctic sea ice drift acceleration: Consequences in terms of projected sea ice thinning and decline. J. Geophys. Res. 2011, 116.
3. Hibler W.D. A dynamic thermodynamic sea ice model. J. Phys. Oceanogr. 1979, 9:815-846.
4. Wilchinsky A.V., Feltham D.L. Modelling the rheology of sea ice as a collection of diamond-shaped floes. J. Non-Newton. Fluid Mech. 2006, 138:22-32.
5. Feltham D.L. Sea ice rheology. Annu. Rev. Fluid Mech. 2008, 40:91-112.
6. Tsamados M., Feltham D.L., Wilchinsky A.V. Impact of a new anisotropic rheology on simulations of Arctic sea ice. J. Geophys. Res. 2013, 118:91-107.
7. Wang K., Wang C. Modeling linear kinematic features in pack ice. J. Geophys. Res. 2009, 114.
8. Girard L., Bouillon S., Jérôme W., Amitrano D., Fichefet T., Legat V. A new modelling framework for sea ice models based on elasto-brittle rheology. Ann. Glaciol. 2011, 52:123-132.
9. Lemieux J.-F., Tremblay B. Numerical convergence of viscous-plastic sea ice models. J. Geophys. Res. 2009, 114.
10. Lemieux J.-F., Knoll D., Tremblay B., Holland D.M., Losch M. A comparison of the Jacobian-free Newton-Krylov method and the EVP model for solving the sea ice momentum equation with a viscous-plastic formulation: a serial algorithm study. J. Comput. Phys. 2012, 231:5926-5944.
11. Zhang J., Hibler W.D. On an efficient numerical method for modeling sea ice dynamics. J. Geophys. Res. 1997, 102:8691-8702.
12. Losch M., Menemenlis D., Campin J.-M., Heimbach P., Hill C. On the formulation of sea-ice models. Part 1: Effects of different solver implementations and parameterizations. Ocean Model. 2010, 33:129-144.
13. Losch M., Danilov S. On solving the momentum equations of dynamic sea ice models with implicit solvers and the elastic viscous-plastic technique. Ocean Model. 2012, 41:42-52.
14. Lemieux J.-F., Tremblay B., Sedláček J., Tupper P., Thomas S., Huard D., Auclair J.-P. Improving the numerical convergence of viscous-plastic sea ice models with the Jacobian-free Newton-Krylov method. J. Comput. Phys. 2010, 229:2840-2852.
15. Cai X.-C., Sarkis M. A restricted additive Schwarz preconditioner for general sparse linear systems. SIAM J. Sci. Comput. 1999, 21:792-797.
16. Hovland P.D., McInnes L.C. Parallel simulation of compressible flow using automatic differentiation and PETSc. Parallel Comput. 2001, 27:503-519.
17. Godoy W.F., Liu X. Parallel Jacobian-free Newton Krylov solution of the discrete ordinates method with flux limiters for 3D radiative transfer. J. Comput. Phys. 2012, 231:4257-4278.
18. Paniconi C., Putti M. A comparison of Picard and Newton iteration in the numerical solution of multidimensional variably saturated flow problems. Water Resour. Res. 1994, 30:3357-3374.
19. Marshall J., Adcroft A., Hill C., Perelman L., Heisey C. A finite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers. J. Geophys. Res. 1997, 102:5753-5766.
20. MITgcm User Manual MITgcm Group, Online documentation, MIT/EAPS, Cambridge, MA 02139, USA. http://mitgcm.org/public/r2_manual/latest/online_documents.
21. Hundsdorfer W., Trompert R.A. Method of lines and direct discretization: a comparison for linear advection. Appl. Numer. Math. 1994, 13:469-490.
22. Roe P. Some contributions to the modelling of discontinuous flows. Lectures in Applied Mathematics 1985, vol. 22:163-193. American Mathematical Society. B. Engquist, S. Osher, R. Somerville (Eds.).
23. Geiger C.A., Hibler W.D., Ackley S.F. Large-scale sea ice drift and deformation: Comparison between models and observations in the western Weddell Sea during 1992. J. Geophys. Res. 1998, 103:21893-21913.
24. Lemieux J.-F., Tremblay B., Thomas S., Sedláček J., Mysak L.A. Using the preconditioned Generalized Minimum RESidual (GMRES) method to solve the sea-ice momentum equation. J. Geophys. Res. 2008, 113.
25. Hunke E.C., Dukowicz J.K. The elastic-viscous-plastic sea ice dynamics model in general orthogonal curvilinear coordinates on a sphere-incorporation of metric terms. Mon. Weather Rev. 2002, 130:1847-1865.
26. Knoll D., Keyes D. Jacobian-free Newton-Krylov methods: a survey of approaches and applications. J. Comput. Phys. 2004, 193:357-397.
27. Giering R., Kaminski T. Recipes for adjoint code construction. ACM Trans. Math. Softw. 1998, 24:437-474.
28. Heimbach P., Hill C., Giering R. An efficient exact adjoint of the parallel MIT general circulation model, generated via automatic differentiation. Future Gener. Comput. Syst. 2005, 21:1356-1371.
29. Saad Y. A flexible inner-outer preconditioned GMRES method. SIAM J. Sci. Comput. 1993, 14:461-469.
30. Eisenstat S.C., Walker H.F. Choosing the forcing terms in an inexact Newton method. SIAM J. Sci. Comput. 1996, 17:16-32.
31. Cai X.-C., Gropp W.D., Keyes D.E., Tidriri M.D. Newton-Krylov-Schwarz methods in CFD. Proceedings of the International Workshop on the Navier-Stokes Equations 2013, 17-30. F.-K. Hebeker, R. Rannacher, G. Wittum (Eds.).
32. Cai X.-C., Gropp W.D., Keyes D.E., Melvin R.G., Young D.P. Parallel Newton-Krylov-Schwarz algorithms for the transonic full potential equation. SIAM J. Sci. Comput. 1998, 19:246-265.
33. Thomas L.H. Elliptic problems in linear differential equations over a network 1949, Watson Sci. Comput. Lab Report, Columbia University, New York.
34. Duff I.S., Meurant G.A. The effect of ordering on preconditioned conjugate gradients. BIT Numer. Math. 1989, 29:635-657.
35. Nguyen A.T., Kwok R., Menemenlis D. Source and pathway of the Western Arctic upper halocline in a data-constrained coupled ocean and sea ice model. J. Phys. Oceanogr. 2012, 43:80-823.
36. Rignot E., Fenty I., Menemenlis D., Xu Y. Spreading of warm ocean waters around Greenland as a possible cause for glacier acceleration. Ann. Glaciol. 2012, 53:257-266.
37. Karcher M., Beszczynska-Möller A., Kauker F., Gerdes R., Heyen S., Rudels B., Schauer U. Arctic ocean warming and its consequences for the Denmark Strait overflow. J. Geophys. Res. 2011, 116.
38. Hill C., Menemenlis D., Ciotti B., Henze C. Investigating solution convergence in a global ocean model using a 2048-processor cluster of distributed shared memory machines. Sci. Program. 2007, 12:107-115.
39. Schlichting J.J.F.M., van der Vorst H.A. Solving 3D block bidiagonal linear systems on vector computers. J. Comput. Appl. Math. 1989, 27:323-330.