The Stokes equations are frequently used to simulate geodynamic processes, including mantle convection, lithospheric dynamics, lava flow, and among others. In this study, the multigrid (MG) method is adopted to solve Stokes and continuity equations with strongly temperature-dependent viscosity. By taking advantage of the enhanced computing power of graphics processing units (GPUs) and the new version of Matlab 2010b, MG codes are optimized through Compute Unified Device Architecture (CUDA). Herein, we illustrate the approach that implements a GPU-based MG solver with Red–Black Gauss–Seidel (RBGS) smoother for the three-dimensional Stokes and continuity equations, in a hope that it helps solve the synthetic incompressible sinking problem in a cubic domain with strongly variable viscosity, and finally analyze our solver’s efficiency on a GPU.
AuthCHarderH (1999) Multigrid solution of convection problems with strongly variable viscosit. Geophysical Journal International137: 793–804.
2.
BachvalovN (1966) On the convergence of a relaxation method with natural constraints on the elliptic operator. Computational Mathematics and Mathematical Physics6: 101–135.
DeubelbeissYKausB (2008) Comparison of Eulerian and Lagrangian numerical techniques for the Stokes equations in the presence of strongly varying viscosity. Physics of the Earth and Planetary Interiors171: 92–111.
6.
ElmanHSilvesterDWathenA (2005) Finite elements and fast iteratice solvers with applications to incompressible fluid dynamics. Oxford University Press.
7.
FedorenkoR (1964) The speed of convergence of one iterative process. Computational Mathematics and Mathematical Physics4: 227–235.
8.
FuruichiMMayDATackleyPJ (2011) Development of a Stokes flow solver robust to large viscosity jumps using a Schur complement approach with mixed precision arithmetir. Journal of Computational Physics230: 8835–8851.
9.
GeryaT (2010) Introduction to numerical Geodynamic modelling. Cambridge: Cambridge University Press.
10.
HackbuschW (1977) On the Convergence of a Multi-grid Iteration Applied to Finite Element Equations. Report, Universitat Koln, pp. 76–12.
11.
HackbuschW (1978) On the multi-grid method applied to difference equations. Computing20: 291–306.
12.
KameyamaMKageyamaASatoT (2005) Multigrid iterative algorighm using pseudo-compressibility for three-dimensional mantle convection with strongly variable viscosity. Journal of Computational Physics206: 162–181.
13.
KlöcknerAPintoNLeeYCatanzaroBIvanovPFasihA (2009) PyCUDA: GPU Run-time Code Generation for High-performance Computing. Technical Report 2009-40, Scientific Computing Group, Brown University, Providence, RI, USA.
LynchD (2005) Numerical Partial Differential Equations for Environmental Scientists and Engineers: A Practical First Course. Berlin: Springer-Verlag.
16.
MayDMoresiL (2008) Preconditioned iterative methods for Stokes flow problems arising in computational geodynamics. Physics of the Earth and Planetary Interiors171: 33–47.
17.
MoresiLZhongSGurnisM (1996) The accuracy of finite element solutions of Stokes’ flow with strongly varying viscosity. Physics of the Earth and Planetary Interiors97: 83–94.
18.
NVIDIA (2011) NVIDIA CUDA C Programming Guide. NVIDIA Corporation.
19.
OosterleeCLorenzF (2006) Multigrid methods for the Stokes system. Computing in Science and Enginering8: 34–43.
20.
PatankarS (1980) Numerical Heat Transfer and Fluid Flow. New York: McGraw-Hill.
21.
PayneLPellW (1960) The Stokes flow problem for a class of axially symmetric bodies. Journal of Fluid Mechanics7: 529–549.
22.
SchubertGTurcotteDLOlsonP (2001) Mantle Convection in the Earth and Planets. Cambridge: Cambridge University Press.
23.
TackleyP (2009) Effects of strongly temperature-dependent viscosity on time-dependent, three-dimensional models of mantle convection. Geophysical Research Letters124: 18–28.
24.
TurcotteDLSchubertG (2002) Geodynamics: Applications of Continuum Physics to Geological Problems, 2nd edn. Cambridge: Cambridge University Press.
25.
WesselingP (1991) An Introduction to Multigrid Methods. New York: John Wiley and Sons.
26.
ZhengLGeryaTKnepleyMYuenDAZhangHShiY (2013) GPU implementation of multigrid solver for Stokes equation with strongly variable viscosity. In: YuenDAWangLChiXJohnssonLGeWShiY (eds.), GPU Solutions to Multi-scale Problems in Science and Engineering. New York: Springer-Verlag, pp. 293–304.
