Implementation of a multigrid solver on a GPU for Stokes equations with strongly variable viscosity based on Matlab and CUDA

Authors:
Liang Zheng;Huai Zhang;Taras Gerya;Matthew Knepley;David A Yuen;Yaolin Shi
Affiliations:
Key Laboratory of Computational Geodynamics, University of Chinese Academy of Sciences, China, College of Earth Science, University of Chinese Academy of Sciences, China, Minnesota Supercomputing ...;Key Laboratory of Computational Geodynamics, University of Chinese Academy of Sciences, China, College of Earth Science, University of Chinese Academy of Sciences, China;Institute of Geophysics, ETH-Zurich, Switzerland;Computational Institute, University of Chicago, Chicago, IL, USA;Minnesota Supercomputing Institute, University of Minnesota, MN, USA, School of Environmental Sciences, China University of Geosciences, Wuhan, China;Key Laboratory of Computational Geodynamics, University of Chinese Academy of Sciences, China, College of Earth Science, University of Chinese Academy of Sciences, China
Venue:
International Journal of High Performance Computing Applications
Year:
2014

Citing 3
Cited 0

Multigrid iterative algorithm using pseudo-compressibility for three-dimensional mantle convection with strongly variable viscosity

Journal of Computational Physics
Multigrid Methods for the Stokes System

Computing in Science and Engineering
Development of a Stokes flow solver robust to large viscosity jumps using a Schur complement approach with mixed precision arithmetic

Journal of Computational Physics

Quantified Score

Hi-index	0.00

Visualization

Abstract

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.