Browse > Article
http://dx.doi.org/10.12989/csm.2012.1.4.325

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)
Publication Information
Coupled systems mechanics / v.1, no.4, 2012 , pp. 325-343 More about this Journal
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
partial differential equations; operator-splitting methods; laplacian operator, diagonal shift operator iterative methods; eigenvalue approach; ice streams; membrane stress; shallow shelf approximation;
Citations & Related Records
연도 인용수 순위
  • Reference
1 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.   DOI   ScienceOn
2 Geiser, J. (2011), Iterative splitting methods for differential equations, CRC Press, New York, USA.
3 Geiser, J. (2010), Iterative splitting methods for differential equations. series: numerical analysis and scientific computing, CRC Press, Taylor and Franics Group, London.
4 Greve, R. and Blatter, H. (2009), Dynamics of ice sheets and glaciers, Springer, Berlin, Germany.
5 Hackbusch, W. (1985), Multi-grid methods and applications, Springer, Berlin, Germany.
6 Hackbusch, W. (1986), Theorie und numerik elliptischer differentialgleichungen, Teubner, Stuttgart, Germany.
7 Hestenes, M. and Stiefel, E. (1952), "Methods of conjugate gradients for solving linear systems", J. Res. Nat. Bureau Standards, 49(6), 409-438.   DOI
8 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.
9 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.   DOI
10 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.   DOI   ScienceOn
11 LAPACK92. LAPACK (1992), Linear algebra package, Software Library, http://www.netlib.org/lapack/.
12 LIS11. LIS (2011), A library of iterative solvers for linear systems, Software Library, http://www.ssisc.org/.
13 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.   DOI   ScienceOn
14 MacAyeal, D.R. (1993), "A tutorial in the use of control methods in ice-sheet modeling", J. Glaciol., 39(131), 91-98.   DOI
15 MacAyeal, D.R. (1995), "Basal friction of Ice Stream E,West Antarctica", J. Glaciol., 41(138), 247-262.   DOI
16 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.
17 Saad, Y. (2003), Iterative methods for sparse linear systems,Society for Industrial Mathematics, Philiadelphia, USA.
18 Smith, B. Bjoerstad, P. and Gropp, W. (1999), Domain decomposition: parallel multilevel methods for elliptic partial differential equations, Cambridge Univ Press, New York, USA.
19 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.
20 Thomas, J.W. (1999), Numerical partial differential equations: conservation laws and elliptic equations, Springer, New York.
21 Toselli, A. and Widlund, O. (2005), Domain decomposition methods: algorithms and theory. Springer, Berlin, Germany.
22 Trottenberg, U., Oosterlee, C. and Schuller, A. (2001), Multigrid, Academic Press, San Diego, USA.
23 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.   DOI   ScienceOn
24 Zeidler, E. (1990), Nonlinear functional analysis and its applications II/B: nonlinear monotone operators, Springer, New York.
25 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.   DOI   ScienceOn
26 Arakawa, A. (1972), Design of the ucla general circulation model, Technical report no. 7, Department of Meteorology, University of California.
27 Barbu, V. (1998), Partial differential equations and boundary value problems, Kluwer Academic Publishers, Dordrecht.
28 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.   DOI
29 Geiser, J. (2008), "Operator-splitting methods respecting eigenvalue problems for convection diffusion and wave equations", Int. J. Appl. Math. Mech., 4(5), 93-104.