The application of a finite volume multigrid method to three-dimensional flow problems in a highly viscous fluid with a variable viscosity
- 1 December 1996
- journal article
- research article
- Published by Taylor & Francis in Geophysical & Astrophysical Fluid Dynamics
- Vol. 83 (3-4) , 261-291
- https://doi.org/10.1080/03091929608208968
Abstract
In this paper we discuss the application of a finite-volume multigrid method to solve three-dimensional thermally driven convection in a highly viscous, incompressible fluid with a variable viscosity. The conservation laws are solved in the primitive variable formulation, A second-order control volume method is used as discretization. Two schemes are used for time stepping, a second-order implicit-explicit scheme based on the Crank-Nicolson and Adams-Bashforth method, and a fully-implicit θ-method. The implicit system of nonlinear equations are solved using multigrid iteration with the SIMPLER method as smoother. In this paper, we describe the implemented multigrid method and investigate its efficiency and the robustness for different viscosity contrasts. Convergence tests showed that with a small modification of the SIMPLER method, the multigrid method exhibits a satisfactory convergence rate even for viscosity contrasts up to 109. Three cases of time-dependent thermally driven convection with viscosity contrasts up to 105 are considered and discussed and the multigrid method has demonstrated its robustness also for these cases. Further, we have also compared the computational efficiency of the two time stepping methods. It appeared that the fully-implicit scheme is a little more efficient than the implicit-explicit scheme for a constant viscosity and it was considerably more efficient for a viscosity contrast of 103Keywords
This publication has 37 references indexed in Scilit:
- Localization of toroidal motion and shear heating in 3‐D high Rayleigh number convection with temperature‐dependent viscosityGeophysical Research Letters, 1995
- Convergence Analysis of the Defect-Correction Iteration for Hyperbolic ProblemsSIAM Journal on Scientific Computing, 1995
- Three-Dimensional Fully Spectral Numerical Method for Mantle Convection with Depth-Dependent PropertiesJournal of Computational Physics, 1994
- Variable viscosity convection experiments with a stress-free upper boundary and implications for the heat transport in the Earth's mantlePhysics of the Earth and Planetary Interiors, 1993
- 3-D Convection With Variable ViscosityGeophysical Journal International, 1991
- A benchmark comparison for mantle convection codesGeophysical Journal International, 1989
- Two Multigrid Methods for Three-Dimensional Problems with Discontinuous and Anisotropic CoefficientsSIAM Journal on Scientific and Statistical Computing, 1987
- Heat transport by variable viscosity convection and implications for the Earth's thermal evolutionPhysics of the Earth and Planetary Interiors, 1984
- The Multi-Grid Method for the Diffusion Equation with Strongly Discontinuous CoefficientsSIAM Journal on Scientific and Statistical Computing, 1981
- Multi-level adaptive solutions to boundary-value problemsMathematics of Computation, 1977