Hostname: page-component-745bb68f8f-l4dxg Total loading time: 0 Render date: 2025-02-11T12:16:01.508Z Has data issue: false hasContentIssue false
  • English
  • Français

Numerical simulations of three-dimensional thermal convection in a fluid with strongly temperature-dependent viscosity

Published online by Cambridge University Press:  26 April 2006

Masaki Ogawa ,
Gerald Schubert  and
Abdelfattah Zebib
Masaki Ogawa
Affiliation:
Department of Earth Sciences, Ehime University, Ehime, 790, Japan
Gerald Schubert
Affiliation:
Department of Earth and Space Sciences and Institute of Geophysics and Planetary Physics, University of California, Los Angeles, CA 90024, USA
Abdelfattah Zebib
Affiliation:
Department of Mechanical and Aerospace Engineering, Rutgers University, New Brunswick, NJ 08903, USA
Get access Rights & Permissions [Opens in a new window]

Abstract

Numerical calculations are presented for the steady three-dimensional structure of thermal convection of a fluid with strongly temperature-dependent viscosity in a bottom-heated rectangular box. Viscosity is assumed to depend on temperature T as exp (− ET), where E is a constant; viscosity variations across the box r (= exp (E)) as large as 105 are considered. A stagnant layer or lid of highly viscous fluid develops in the uppermost coldest part of the top cold thermal boundary layer when r > rc1, where r = rc1 ≡ 1.18 × 103Rt0.309 and Rt is the Rayleigh number based on the viscosity at the top boundary. Three-dimensional convection occurs in a rectangular pattern beneath this stagnant lid. The planform consists of hot upwelling plumes at or near the centre of a rectangle, sheets of cold sinking fluid on the four sides, and cold sinking plume concentrations immersed in the sheets. A stagnant lid does not develop, i.e. convection involves all of the fluid in the box when r < rc1. The whole-layer mode of convection occurs in a three-dimensional bimodal pattern when r > rc2 = 3.84 × 106Rt−1.35. The planform of the convection is rectangular with the coldest parts of the sinking fluid and the hottest part of the upwelling fluid occurring as plumes at the four corners and at the centre of the rectangle, respectively. Both hot uprising plumes and cold sinking plumes have sheet-like extensions, which become more well-developed as r increases. The whole-layer mode of convection occurs as two-dimensional rolls when r < min (rc1, rc2). The Nusselt number Nu depends on the viscosity at the top surface more strongly in the regime of whole-layer convection than in the regime of stagnant-lid convection. In the whole-layer convective regime, Nu depends more strongly on the viscosity at the top surface than on the viscosity at the bottom boundary.

Type
Research Article
Information
Journal of Fluid Mechanics , Volume 233 , December 1991 , pp. 299 - 328
Copyright
© 1991 Cambridge University Press

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

Afrid, M. 1990 Three-dimensional natural convection. PhD dissertation, Rutgers University, The State University of New Jersey
Anderson, D. L. 1981 Plate tectonics on Venus. Geophys. Res. Lett. 8, 309–311.Google Scholar
Baumgardner, J. R. 1985 Three-dimensional treatment of convective flow in the Earth's mantle. J. Stat. Phys. 39, 501–511.Google Scholar
Bercovici, D., Schubert, G. & Glatzmaier, G. A. 1989a Three-dimensional spherical models of convection in the Earth's mantle. Science 244, 893–1016.Google Scholar
Bercovici, D., Schubert, G. & Glatzmaier, G. A. 1989b Influence of heating mode on three-dimensional mantle convection. Geophys. Res. Lett. 16, 617–620.Google Scholar
Bercovici, D., Schubert, G., Glatzmaier, G. A. & Zebib, A. 1989c Three-dimensional thermal convection in a spherical shell. J. Fluid Mech. 206, 75–104.Google Scholar
Bindschadler, D. L., Schubert, G. & Kaula, W. M. 1990 Mantle flow tectonics and the origin of Ishtar terra, Venus. Geophys. Res. Lett. 17, 1345–1348.Google Scholar
Blankenbach, B., Busse, F., Christensen, U., Cserepes, L., Gunkel, D., et al. 1989 A benchmark comparison for mantle convection codes. Geophys. J. Intl 98, 23–38.Google Scholar
Bolton, E. W. & Busse, F. H. 1985 Stability of convection rolls in a layer with stress-free boundaries. J. Fluid Mech. 150, 487–498.Google Scholar
Busse, F. H. 1967 The stability of finite amplitude cellular convection and its relation to an extremum principle. J. Fluid Mech. 30, 625–649.Google Scholar
Busse, F. H. & Frick, H. 1985 Square pattern convection in fluids with strongly temperature-dependent viscosity. J. Fluid Mech. 150, 451–465.Google Scholar
Chandrasekhar, S. 1961 Hydrodynamic and Hydromagnetic Stability. Clarendon.
Christensen, U. R. 1984a Heat transport by variable viscosity convection and implications for the Earth's thermal evolution. Phys. Earth Planet. Inter. 35, 264–282.Google Scholar
Christensen, U. R. 1984b Convection with pressure and temperature dependent non-Newtonian rheology. Geophys. J. R. Astron. Soc. 77, 343–384.Google Scholar
Christensen, U. & Harder, H. 1991 Three-dimensional convection with variable viscosity. Geophysical J. Intl 104, 213–226.Google Scholar
Cserepes, L., Rabinowicz, M. & Rosemberg-Borot, C. 1988 Three-dimensional convection in a two-layer mantle. J. Geophys. Res. 93, 12009–12025.Google Scholar
Glatzmaier, G. A. 1988 Numerical simulation of mantle convection: Time-dependent, three-dimensional, compressible, spherical shell. Geophys. Astrophys. Fluid Dyn. 43, 223–264.Google Scholar
Hager, B. H. & Gurnis, M. 1987 Mantle convection and the state of the Earth's interior. Rev. Geophys. 25, 1277–1285.Google Scholar
Houseman, G. 1988 The dependence of convection planform on mode of heating. Nature 332, 346–349.Google Scholar
Kaula, W. M. 1980 Material properties for mantle convection consistent with observed surface fields. J. Geophys. Res. 85, 7031–7044.Google Scholar
Machetel, P., Rabinowicz, M. & Bernardet, P. 1986 Three-dimensional convection in spherical shells. Geophys. Astrophys. Fluid Dyn. 37, 57–84.Google Scholar
Parmentier, E. M. & Turcotte, D. L. 1978 Two-dimensional mantle flow beneath a rigid accreting lithosphere. Phys. Earth Planet. Inter. 17, 281–289.Google Scholar
Patankar, S. V. 1980 Numerical heat transfer and fluid flow. Hemisphere, 197 pp.
Richter, F. M. 1978 Experiments on the stability of convection rolls in fluids whose viscosity depends on temperature. J. Fluid Mech. 89, 553–560.Google Scholar
Schnaubelt, M. & Busse, F. H. 1989 On the stability of two-dimensional convection rolls in an infinite Prandtl number fluid with stress-free boundaries. Z. angew. Math. Phys. 40, 153–162.Google Scholar
Schubert, G., Cassen, P. & Young, R. E. 1979 Subsolidus convective cooling histories of terrestrial planets. Icarus 38, 192–211.Google Scholar
Stengel, K. C., Oliver, D. S. & Booker, J. R. 1982 Onset of convection in a variable viscosity fluid. J. Fluid Mech. 120, 411–431.Google Scholar
Travis, B., Olson, P. & Schubert, G. 1990a The transition from two-dimensional to three-dimensional planforms in infinite-Prandtl-number thermal convection. J. Fluid Mech. 216, 71–91.Google Scholar
Travis, B., Weinstein, S. A. & Olson, P. 1990b Three-dimensional convection planforms with internal heat generation. Geophys. Res. Lett. 17, 243–246.Google Scholar
Weertman, J. 1970 The creep strength of the Earth's mantle. Rev. Geophys. Space Phys. 8, 145–168.Google Scholar
White, D. B. 1988 The planforms and onset of convection with a temperature-dependent viscosity. J. Fluid Mech. 191, 247–286.Google Scholar
Zebib, A., Goyal, K. & Schubert, G. 1985 Convective motion in a spherical shell. J. Fluid Mech. 152, 39–48.Google Scholar