Variable Viscosity Thermal Convection at Infinite Prandtl Number in a Thick Spherical Shell
Yang, Woo-Sun
This item is only available for download by members of the University of Illinois community. Students, faculty, and staff at the U of I may log in with your NetID and password to view the item. If you are trying to access an Illinois-restricted dissertation or thesis, you can request a copy through your library's Inter-Library Loan office or purchase a copy directly from ProQuest.
Permalink
https://hdl.handle.net/2142/86549
Description
Title
Variable Viscosity Thermal Convection at Infinite Prandtl Number in a Thick Spherical Shell
Author(s)
Yang, Woo-Sun
Issue Date
1997
Doctoral Committee Chair(s)
Hsui, Albert T.
Department of Study
Geology
Discipline
Geology
Degree Granting Institution
University of Illinois at Urbana-Champaign
Degree Name
Ph.D.
Degree Level
Dissertation
Keyword(s)
Geophysics
Language
eng
Abstract
Variable viscosity thermal convection at infinite Prandtl number in a thick spherical shell is modeled using the finite element method. The discretized tensor Navier-Stokes equations are solved with the multigrid method applied to the spherical elements. Numerically convergent solutions can be obtained when specially tailored matrix dependent transfer is introduced in the multigrid method and line Jacobi relaxation is applied in the radial direction of the shell. The conjugate gradient method is used to correct both pressure and velocity fields simultaneously to satisfy the incompressibility condition. Numerical calculations are performed on a Cray T3D to investigate internally heated time dependent thermal convection for Rayleigh numbers between 10$\sp5$ and 10$\sp6$ with both Newtonian and non-Newtonian viscosity which can be dependent on temperature and pressure (depth). Linear downwellings and diffuse plume-like upwellings are dominant features in the numerical calculations. When Newtonian rheology is applied throughout the shell, depth- dependence of viscosity is the dominant rheological influence controlling the flow characteristics--lateral variation of viscosity plays a secondary role. Both a high viscosity lid and a viscosity increase below 660 km depth provide a general shift to lower harmonic degrees in both thermal anomalies and kinetic energy. With non-Newtonian rheology, narrow linear downwellings appear at the surface, and the upper boundary layer becomes very mobile. Thermal heterogeneities are greatly reduced throughout the shell. When non-Newtonian rheology is limited to the top 300 km as suggested by Karato and Wu (1993), the lateral flow increases because of the change in relative viscosity around downwelling regions at the transition depth. The downwelling network becomes more irregular and exhibits irregular plate-like characteristics. The planform of upwellings is also greatly affected by the large lateral flow and, in some case, shows narrow extended shapes. A significant increase in toroidal energy (toroidal-poloidal rms velocity ratio is 19-35% compared to less than 10% for Newtonian viscosity cases) is obtained with this hybrid rheology. The knowledge acquired from the simulations is applied to various depth regions in the global seismic tomography model proposed by Su et al. (1994) with a geodynamical interpretation.
Use this login method if you
don't
have an
@illinois.edu
email address.
(Oops, I do have one)
IDEALS migrated to a new platform on June 23, 2022. If you created
your account prior to this date, you will have to reset your password
using the forgot-password link below.
