Development of a Stokes flow solver robust to large viscosity jumps using a Schur complement approach with mixed precision arithmetic

Journal of Computational Physics

Volumn 230, Issue 24, 2011, Pages 8835-8851

  Furuichi, Mikito ^a   May, Dave A ^b   Tackley, Paul J ^b  

^a JAPAN AGENCY FOR MARINE EARTH SCIENCE AND TECHNOLOGY   (Japan)

^b ETH ZURICH   (Switzerland)

Author keywords

Double double; High precision; Mantle convection; Multi grid; Preconditioner; Stokes flow; Viscosity jump

Indexed keywords

CONVERGENCE OF NUMERICAL METHODS; ECONOMIC AND SOCIAL EFFECTS; GEODYNAMICS; ITERATIVE METHODS; PROBLEM SOLVING;

DOUBLE-DOUBLE; HIGH-PRECISION; MANTLE CONVECTION; MIXED PRECISION; MULTIGRIDS; PRECISION ARITHMETIC; PRECONDITIONERS; SCHUR COMPLEMENT; STOKES FLOWS; VISCOSITY JUMP;

VISCOSITY;

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

References (40)

1. Schubert G., Turcotte D.L., Olson P. Mantle Convection in the Earth and planets 2001, Cambridge University Press, pp. 940.
2. R.G. Gordon, Diffuse oceanic plate boundaries: Strain rates, vertically averaged rheology, and comparisons with narrow plate boundaries and stable plate interiors, in: M.A. Richards, R.G. Gordon, R.D. Van der Hilst (Eds.), History and Dynamics of Global Plate Motions, Geophysical Monograph 121, American Geophysical Union, Washington, DC, 2000, pp. 143-159.
3. Schmeling H., Babeyko A., Enns A., Faccenna C., Funiciello F., Gerya T., Golabek G., Grigull S., Kaus B.J.P., Morra G., Schmalholz S., van Hunen J. A benchmark comparison of spontaneous subduction models - towards a free surface. Phys. Earth Planet. Interiors 2008, 171:198-223.
4. Furuichi M., Kameyama M., Kageyama A. Three-Dimensional Eulerian method for large deformation of viscoelastic fluid: toward plate-mantle simulation. J. Comput. Phys. 2008, 227:4977-4997.
5. Furuichi M., Kameyama M., Kageyama A. Validity test of a Stokes flow solver by fluid rope coiling: toward plate-mantle simulation. Phys. Earth Planet. Interiors 2009, 176:44-53.
6. Stadler G., Gurnis M., Burstedde C., Wilcox L.C., Alisic L., Ghattas O. The dynamics of plate tectonics and mantle flow: from local to global scales. Science 2010, 329(5995):1033-1038.
7. Jadamec M.A., Billen M.I. Reconciling rapid 3-D mantle flow and surface plate motions near the Eastern Alaska slab edge. Nature 2010, 465:338-341.
8. Tackley P.J. 1993, Effects of strongly temperature-dependent viscosity on time-dependent, three-dimensional models of mantle convection. Geophys. Res. Lett. 2009, 20(20):2187-2190.
9. Tackley P.J., Xie S. Stag3D: a code for modeling thermo-chemical multiphase convection in Earth's mantle. Proceedings of the Second MIT Conference on Computational Fluid and Solid Mechanics 2003, 1524-1527. Elsevier B.V, Amsterdam. K.J. Bathe (Ed.).
10. Tackley P.J. Modelling compressible mantle convection with large viscosity contrasts in a three-dimensional spherical shell using the Yin-Yang grid. Phys. Earth Planet. Interiors. 2008, 171:7-18.
11. Moresi L., Solomatov V. Numerical investigation of 2D convection with extremely large viscosity variations. Phys. Fluids 1995, 7:2154-2162.
12. Trompert R.A., Hansen U. The application of a finite volume multigrid method to three-dimensional flow problems in a highly viscous fluid with a variable viscosity. Geophys. Astrophys. Fluid Dyn. 1996, 83(3):261-291.
13. Auth C., Harder H. Multigrid solution of convection problems with strongly variable viscosity. Geophys. Res. Lett. 1999, 137(3):793-804.
14. Albers M. A local mesh refinement multigrid method for 3-D convection problems with strongly variable viscosity. J. Comput. Phys. 2000, 160(1):126-150.
15. Kameyama M., Kageyama A., Sato T. Multigrid iterative algorithm using pseudo-compressibility for three-dimensional mantle convection with strongly variable viscosity. J. Comput. Phys. 2005, 206:162-181.
16. Kameyama M. A multigrid-based mantle convection simulation code and its optimization to the Earth Simulator. J. Earth Simul. 2005, 4:2-10.
17. Zhong S., McNamara A., Tan E., Moresi L., Gurnis M. A benchmark study on mantle convection in a 3-D spherical shell using CitcomS. Geochem. Geo-phys. Geosyst. 2008, 9:Q10017.
18. Burstedde C., Ghattas O., Stadler G., Tu T., Wilcox L.C. Parallel scalable adjoint-based adaptive solution for variable-viscosity Stokes flows. Comput. Methods Appl. Mech. Eng. 2009, 198:1691-1700.
19. C. Burstedde, O. Ghattas, G. Stadler, E. Tan, T. Tu, L.C. Wilcox, S. Zhong, Scalable adaptive mantle convection simulation on petascale supercomputers, in: Proceedings of ACM/IEEE SC08, 2008.
20. Patankar S.V. Numerical Heat Transfer and Fluid Flow 1980, McGraw-Hill, New York.
21. May D.A., Moresi L. Preconditioned iterative methods for Stokes flow problems arising in computational geodynamics. Phys. Earth Planet Interiors. 2008, 171:33-47.
22. Elman H.C. Preconditioning for the steady-state Navier-Stokes equations with low viscosity. SIAM J. Sci. Comput. 1996, 20:1299-1316.
23. Elman H.C., Silvester D.J., Wathen A.J. Finite Elements and Fast Iterative Solvers 2005, Oxford University Press, Oxford, p. 353.
24. Elman H.C., Howle V.E., Shadid J., Shuttleworth R., Tuminaro R. Block preconditioners based on approximate commutators. SIAM J. Sci. Comput. 2006, 27(5):1651-1668.
25. Silvester D., Elman H., Kay D., Wathen A. Efficient preconditioning of the linearized Navier-Stokes equations for incompressible flow. J. Comput. Appl. Math. 2001, 128:261-279.
26. Grinevich P.P., Olishanskii M.A. An iterative method for the Stokes-type problem with variable viscosity. SIAM J. Sci. Comput 2009, 31(5):3959-3978.
27. Geenen T., ur Rehman M., MacLachlan S.P., Segal G., Vuik C., van den Berg A.P., Spakman W. Scalable robust solvers for unstructured FE geodynamic modeling applications: solving the Stokes equation for models with large localized viscosity contrasts. Geochem. Geophys. Geosyst. 2009, 10:Q09002. 10.1029/2009GC002526.
28. D.H. Bailey, High-Precision Software Directory. http://www.crd.lbl.gov/~dhbailey/mpdist/.
29. Y. Hida, X. Li, D. Bailey, Algorithm for quad-double precision floating point arithmetic, in: Proceedings of the 15th IEEE Symposium on Computer Architecture, 2001, pp. 287-302.
30. Uzawa H. Iterative methods for concave programming. Studies in Linear and Nonlinear Programming 1958, 54-165. Stanford University Press, Stanford, CA. K.J. Arrow, L. Hurwicz, H. Uzawa (Eds.).
31. Benzi M., Golub G.H., Liesen J. Numerical solutions of saddle point problems. Acta Numer. 2005, 1-137.
32. Ramage A., Wathen A.J. Iterative solution techniques for the stokes and Navier-Stokes equations. Int. J. Numer. Methods Fluids 1994, 19(1):67-83.
33. Trottenberg U., Oosterlee C., Schuller A. Multigrid 2001, Academic Press, New York, p. 316.
34. Vanka S.P. Block-implicit multigrid solutions of Navier-Stokes equations in primitive variables. J. Comput. Phys. 1986, 65:138-158.
35. Wathen A.J., Silvester D.J. Fast iterative solution of stabilised Stokes systems I: using simple diagonal preconditioners. SIAM J. Numer. Anal. 1993, 30:630-649.
36. Sleijpen G.L.G., vander Vorst H.A., Modersitzki J. Differences in the effects of rounding errors in Krylov solvers for symmetric indefinite linear systems. SIAM J. Matrix Anal. Appl. 2001, 22 3:726-751.
37. H. Kotakemori, A. Fujii, H. Hasegawa, A. Nishida, Implementation of fast quad precision operation and acceleration with SSE2 for iterative solver library, IPSJ Trans. Adv. Comput. Syst. 1(1) (2008) 75-84 (in Japanese).
38. Shewchuk J.R. Adaptive precision floating-point arithmetic and fast robust geometric predicates. Discrete Comput. Geom. 1997, 18:305-363.
39. Dekker T. A floating-point technique for extending the available precision. Numer. Math. 1971, 18:224-242.
40. Knuth D. The Art of Computer Programming. Seminumerical Algorithms 1998, vol. 2. Addison Wesley, Reading, Massachusetts. first ed.