DOI QR코드

DOI QR Code

Operator-splitting methods respecting eigenvalue problems for shallow shelf equations with basal drag

  • Geiser, Jurgen (Department of Physics, Ernst-Moritz-Arndt University of Greifswald) ;
  • Calov, Reinhard (Potsdam Institute for Climate Impact Research)
  • Received : 2012.06.25
  • Accepted : 2012.11.14
  • Published : 2012.12.25

Abstract

We present different numerical methods for solving the shallow shelf equations with basal drag (SSAB). An alternative approach of splitting the SSAB equation into a Laplacian and diagonal shift operator is discussed with respect to the underlying eigenvalue problem. First, we solve the equations using standard methods. Then, the coupled equations are decomposed into operators for membranes stresses, basal shear stress and driving stress. Applying reasonable parameter values, we demonstrate that the operator of the membrane stresses is much stiffer than the operator of the basal shear stress. Here, we could apply a new splitting method, which alternates between the iteration on the membrane-stress operator and the basal-shear operator, with a more frequent iteration on the operator of the membrane stresses. We show that this splitting accelerates and stabilize the computational performance of the numerical method, although an appropriate choice of the standard method used to solve for all operators in one step speeds up the scheme as well.

Keywords

References

  1. Arakawa, A. (1972), Design of the ucla general circulation model, Technical report no. 7, Department of Meteorology, University of California.
  2. Barbu, V. (1998), Partial differential equations and boundary value problems, Kluwer Academic Publishers, Dordrecht.
  3. Bueler, E. and Brown, J. (2009), "Shallow shelf approximation as a "sliding law" in a thermomechanically coupled ice sheet model", J. Geophys. Res., 114(3), F03008. https://doi.org/10.1029/2008JF001179
  4. Bueler, E., Lingle, C.S., Kallen-Brown, J.A., Covey, D.N. and Bowman, L.N. (2005), "Exact solutions and verification of numerical models for isothermal ice sheets", J. Glaciol., 51(173), 291-306. https://doi.org/10.3189/172756505781829449
  5. Geiser, J. (2008), "Operator-splitting methods respecting eigenvalue problems for convection diffusion and wave equations", Int. J. Appl. Math. Mech., 4(5), 93-104.
  6. Geiser, J. (2009), "Operator-splitting methods in respect of eigenvalue problems for nonlinear equations and applications to burgers equations", J. Comput. Appl. Math., 231(2), 815-827. https://doi.org/10.1016/j.cam.2009.05.009
  7. Geiser, J. (2011), Iterative splitting methods for differential equations, CRC Press, New York, USA.
  8. Geiser, J. (2010), Iterative splitting methods for differential equations. series: numerical analysis and scientific computing, CRC Press, Taylor and Franics Group, London.
  9. Greve, R. and Blatter, H. (2009), Dynamics of ice sheets and glaciers, Springer, Berlin, Germany.
  10. Hackbusch, W. (1985), Multi-grid methods and applications, Springer, Berlin, Germany.
  11. Hackbusch, W. (1986), Theorie und numerik elliptischer differentialgleichungen, Teubner, Stuttgart, Germany.
  12. Hestenes, M. and Stiefel, E. (1952), "Methods of conjugate gradients for solving linear systems", J. Res. Nat. Bureau Standards, 49(6), 409-438. https://doi.org/10.6028/jres.049.044
  13. Jenson, J.W. (1994), A nonlinear numerical model of the Lake Michigan lobe, Laurentide ice sheet. Doctoral thesis, Department of Geology, Oregon State University, Corvallis, OR, USA.
  14. Joughin, I., Fahnestock, M. MacAyeal, D.R., Bamber, J.L. and Gogenini, S.P. (2001), "Observation and analysis of ice flow in the largest Greenland ice stream", J. Geophys. Res., 106(24), 34021-34034. https://doi.org/10.1029/2001JD900087
  15. Kanney, J., Miller, C. and Kelley, C. (2003), "Convergence of iterative split-operator approaches for approximating nonlinear reactive transport problems", Adv. Water. Resour., 26(3), 247-261. https://doi.org/10.1016/S0309-1708(02)00162-8
  16. LAPACK92. LAPACK (1992), Linear algebra package, Software Library, http://www.netlib.org/lapack/.
  17. LIS11. LIS (2011), A library of iterative solvers for linear systems, Software Library, http://www.ssisc.org/.
  18. Storti, M.M., Nigro, N. and Idelsohn, S. (1991), "Multigrid methods and adaptive refinement techniques in elliptic problems by finite element methods", Comput. Method. Appl. M., 93(1),13-30. https://doi.org/10.1016/0045-7825(91)90113-K
  19. MacAyeal, D.R. (1993), "A tutorial in the use of control methods in ice-sheet modeling", J. Glaciol., 39(131), 91-98. https://doi.org/10.1017/S0022143000015744
  20. MacAyeal, D.R. (1995), "Basal friction of Ice Stream E,West Antarctica", J. Glaciol., 41(138), 247-262. https://doi.org/10.1017/S0022143000016154
  21. Pollard, D. and DeConto, R.M. (2007), A coupled ice-sheet/ice-shelf/sediment model applied to a marine-margin flowline: forced and unforced variations, (Eds., Hambrey, M.J., Christoffersen, P., Glasser, N.F. and Hubbard, B.), Glacial Sedimentary Processes and Products, International Association of Sedimentologists Special Publication No. 39, 37-52. Blackwell Publ.
  22. Saad, Y. (2003), Iterative methods for sparse linear systems,Society for Industrial Mathematics, Philiadelphia, USA.
  23. Smith, B. Bjoerstad, P. and Gropp, W. (1999), Domain decomposition: parallel multilevel methods for elliptic partial differential equations, Cambridge Univ Press, New York, USA.
  24. Smith, B., Boyle, J., Dongarra, J., Garbow, B., Ikebe, Y., Klema, V. and Moler, C. (1976), Matrix eigensystem routines, EISPACK Guide. Lecture Notes in Computer Science, Springer, New York.
  25. Thomas, J.W. (1999), Numerical partial differential equations: conservation laws and elliptic equations, Springer, New York.
  26. Toselli, A. and Widlund, O. (2005), Domain decomposition methods: algorithms and theory. Springer, Berlin, Germany.
  27. Trottenberg, U., Oosterlee, C. and Schuller, A. (2001), Multigrid, Academic Press, San Diego, USA.
  28. van Gijzen, M.B., Erlangga, Y.A. and VuikKanney, C. (2007), "Spectral analysis of the discrete helmholtz operator preconditioned with a shifted laplacian", SIAM J. Sci. Comput., 29(5),1942-1958. https://doi.org/10.1137/060661491
  29. Zeidler, E. (1990), Nonlinear functional analysis and its applications II/B: nonlinear monotone operators, Springer, New York.