A matrix-dependent transfer multigrid method for strongly variable viscosity infinite Prandtl number thermal convection

Abstract We apply a two-dimensional Cartesian finite element treatment to investigate infinite Prandtl number thermal convection with temperature, strain rate and yield stress dependent rheology using parameters in the range estimated for the mantles of the terrestrial planets. To handle the strong viscosity variations that arise from such nonlinear rheology in solving the momentum equation, we exploit a multigrid method based on matrix-dependent intergrid transfer and the Galerkin coarse grid approximation. We observe that the matrix-dependent transfer algorithm provides an exceptionally robust and efficient means for solving convection problems with extreme viscosity gradients. Our algorithm displays a convergence rate per multigrid cycle about five times better than what other published methods (e.g., CITCOM of Moresi and Solomatov, 1995) offer for cases with similar extreme viscosity variation. The algorithm is explained in detail in this paper. When this method is applied to problems with temperature and strain rate dependent rheologies, we obtain strongly time dependent solutions characterized by episodic avalanching of cold material from the upper boundary layer to the bottom of the convecting domain for a significantly broad range of parameter values. In particular, we observe this behavior for the relatively simple case of temperature dependent Newtonian rheology with a plastic yield stress. The intensity and temporal character of the episodic behavior depends sensitively on the yield stress value. The regions most strongly affected by the yield stress are thickened portions of the cold upper boundary layer which can suddenly become unstable and form downgoing diapirs. These computational results suggest that the finite yield properties of silicate rocks must play a vitally important role in planetary mantle dynamics. Although our example calculations were selected mainly to illustrate the power of our multigrid method, they suggest that many possible exotic behaviors in planetary mantles have yet to be discovered.

[1]  R. Trompert,et al.  Mantle convection simulations with rheologies that generate plate-like behaviour , 1998, Nature.

[2]  Ulrich Hansen,et al.  On the Rayleigh number dependence of convection with a strongly temperature-dependent viscosity , 1998 .

[3]  R. Kettler.  Analysis and comparison of relaxation schemes in robust multigrid and preconditioned conjugate gradient methods , 1982 .

[4]  Long Chen.  INTRODUCTION TO MULTIGRID METHODS , 2005 .

[5]  David A. Yuen,et al.  The effects of a composite non-Newtonian and Newtonian rheology on mantle convection , 1993 .

[6]  J. Weertman.  The creep strength of the Earth's mantle , 1970 .

[7]  D. Silvester,et al.  ITERATIVE METHODS FOR STABILIZED MIXED VELOCITY PRESSURE FINITE-ELEMENTS , 1992 .

[8]  Christian Wagner.  A Multigrid Method with Matrix-Dependent Transfer Operators Applied in Groundwater Science , 1994 .

[9]  John R. Baumgardner,et al.  Three-dimensional treatment of convective flow in the earth's mantle , 1985 .

[10]  S. Patankar.  Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.

[11]  A. Reusken,et al.  Multigrid with Matrix-dependent Transfer Operators for Convection-diffusion Problems , 1994 .

[12]  A. Wathen,et al.  Iterative solution techniques for the stokes and Navier‐Stokes equations , 1994 .

[13]  Louis Moresi,et al.  Numerical investigation of 2D convection with extremely large viscosity variations , 1995 .

[14]  B. Evans,et al.  Strength of the lithosphere: Constraints imposed by laboratory experiments , 1995 .

[15]  A. Brandt,et al.  The Multi-Grid Method for the Diffusion Equation with Strongly Discontinuous Coefficients , 1981 .

[16]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[17]  Patrick Wu,et al.  Rheology of the Upper Mantle: A Synthesis , 1993, Science.

[18]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[19]  J. T. Ratcliff,et al.  Transitions in thermal convection with strongly variable viscosity , 1997 .

[20]  Louis Moresi,et al.  Mantle convection with a brittle lithosphere: thoughts on the global tectonic styles of the Earth and Venus , 1998 .

[21]  Woo-Sun Yang.  Variable Viscosity Thermal Convection at Infinite Prandtl Number in a Thick Spherical Shell , 1997 .

[22]  G. Ranalli.  The Microphysical Approach to Mantle Rheology , 1991 .

[23]  M. F. Ashby,et al.  On the rheology of the upper mantle , 1973 .

[24]  A. J. Anderson,et al.  Convection in the Earth's mantle , 1984 .

[25]  D. Yuen,et al.  Rheological transition in mantle convection with a composite temperature-dependent, non-Newtonian and Newtonian rheology , 1995 .

[26]  J. Weertman,et al.  High Temperature Creep of Rock and Mantle Viscosity , 1975 .

[27]  R. Verfürth.  A combined conjugate gradient - multi-grid algorithm for the numerical solution of the Stokes problem , 1984 .

[28]  D. Yuen,et al.  Convectively induced transition in mantle rheological behavior , 1995 .

[29]  Wolfgang Hackbusch,et al.  Multi-grid convergence theory , 1982 .

[30]  G. Glatzmaier.  Numerical simulations of mantle convection: time-dependent, three-dimensional, compressible, spherical shell , 1988 .

[31]  P. M. deZeeuw.  Matrixdependent prolongations and restrictions in a blackbox multigrid solver , 1988 .

[32]  P. M. De Zeeuw,et al.  Matrix-dependent prolongations and restrictions in a blackbox multigrid solver , 1990 .

[33]  D. Turcotte,et al.  Studies of finite amplitude non‐Newtonian thermal convection with application to convection in the Earth's mantle , 1976 .

[34]  D. Yuen,et al.  Is the lower‐mantle rheology Newtonian today? , 1996 .

[35]  Louis Moresi,et al.  Heat transport efficiency for stagnant lid convection with dislocation viscosity: Application to Mars and Venus , 1998 .

[36]  Ulrich Hansen,et al.  The application of a finite volume multigrid method to three-dimensional flow problems in a highly viscous fluid with a variable viscosity , 1996 .

[37]  C. Doglioni,et al.  Differential Rotation Between Lithosphere and Mantle' A Consequence of Lateral Mantle Viscosity Variations , 1991 .

[38]  Thomas J. Ahrens,et al.  Rock physics & phase relations : a handbook of physical constants , 1995 .