Abstract
We have used a numerical scheme based on higher-order finite differences to investigate effects of adiabatic heating and viscous dissipation on 3-D rapidly rotating thermal convection in a Cartesian box with an aspect-ratio of 2⊗2⊗1. Although we omitted coupling with the magnetic field, which can play a key role in the dynamics of the Earth's core, the understanding of non-linear rotating convection including realistic thermodynamic effects is a necessary prerequisite for understanding the full complexity of the Earth's core dynamics. The system of coupled partial differential equations has been solved in terms of the principal variables vorticity ω, vector potential A and temperature T. The use of the vector potential A allows the velocity field to be calculated with one spatial differentiation in contrast to the spheroidal and toroidal function approach. The temporal evolution is governed by a coupled time-dependent system consisting of ω and T. The equations are discretized in all directions by using an eighth-order, variable spaced scheme. Rayleigh number Ra of 106, Taylor number Ta of 108 and a Prandtl number Pr of 1 have been employed. The dissipation number of the outer core was taken to be 0.2. A stretched grid has been employed in the vertical direction for resolving the thin shear boundary layers at the top and bottom. This vertical resolution corresponds to around 240 regularly spaced points with an eighth-order accuracy. For the regime appropriate to the Earth's outer core, the dimensionless surface temperature T 0 takes a large value, around 4. This large value in the adiabatic heating/cooling term is found to cause stabilization of both the temperature and velocity fields.
Similar content being viewed by others
References
Bercovici D., Schubert G. and Glatzmaier G.A., 1992. Three-dimensional convection of an infinite Prandtl number compressible fluid in a basally heated spherical shell. J. Fluid Mech., 239, 683–719.
Boehler R., 1996. Melting temperature of the Earth's mantle and core: Earth's thermal structure. Annu. Rev. Earth Planet. Sci., 24, 15–40.
Braginsky S.I. and Roberts P.H., 1995. Equations governing convection in the Earth's core and the geodynamo. Geophys. Astrophys. Fluid Dynamics, 79, 1–97.
Breuer D. and Spohn T., 1993. Cooling of the Earth, Urey ratios, and the problem of potassium in the core. Geophys. Res. Lett., 20, 1655–1658.
Brummell N.H., Cattaneo F. and Toomre J., 1995. Turbulent dynamics in the solar convection zone. Science, 269, 1370–1379.
Brummell N.H., Hurlburt N.E. and Toomre J., 1998. Turbulent compressible convection with rotation. II. Mean-flows and differential rotation. Astrophys. J., 493, 955–969.
Busse F.H., Grote E. and Tilgner A., 1998. On convection driven dynamos in rotating spherical shells. Stud. geophys. geod., 42, 211–223.
Cardin P. and Olson P., 1994. Chaotic thermal convection in a rapidly rotating spherical shell: consequences for flow in the outer core. Phys. Earth Planet. Inter., 82, 235–259.
Christensen U., 1989. Mantle rheology, constitution and convection. In: W.R. Peltier (Ed.), Mantle Convection-Plate Tectonics and Global Dynamics. Gordon and Breach Science Publishers, New York, 595–655.
Christensen U., Olson P. and Glatzmaier G.A., 1998. A dynamo model interpretation of geomagnetic field structures. Geophys. Res. Lett., 25, 1565–1568.
Christensen U. and Yuen D.A., 1985. Layered convection induced by phase transitions. J. Geophys. Res., 90, 10291–10300.
Demircan A., Scheel S. and Seehafer N., 2000. Heteroclinic behavior in rotating Rayleigh-Benard convection. Europ. Phys. J. B, 13, 765–775.
Douglas J. Jr. and Rachford H.H., 1956. On the Numerical Solution of Heat Conduction Problems in Two and Three Space Variables. Trans. Amer. Math. Soc., 82, 421–430.
Fornberg B., 1992. Fast generation of weights in finite-difference formulas. In: G.D. Byrne and W.E. Schiesser (Eds.), Recent Developments in Numerical Methods and Software for ODEs/DAEs/PDEs. World Scientific Publ. Co., Singapore.
Fornberg B., 1995. A pseudospectral approach for polar and spherical geometries. SIAM J. Sci. Comput., 16, 1071–1081.
Fornberg B., 1996. A Practical Guide to Pseudospectral Methods. Cambridge Univ. Press, Cambridge.
Fornberg B. and Merrill D., 1997. Comparison of finite-difference and pseudospectral methods for convective flow over a sphere. Geophys. Res. Lett., 24, 3245–3248.
Glatzmaier G.A. and Roberts P.H., 1995. A three-dimensional convective dynamo solution with rotating and finitely conducting inner core and mantle. Phys. Earth Planet. Inter., 91, 63–75.
Glatzmaier G.A. and Roberts P.H., 1996. An inelastic evolutionary geodynamo simulation driven by compositional and thermal convection. Physica D, 97, 81–94.
Hansen U.P. and Yuen D.A., 1995. Different consequences from compressibility on thermal and thermal-chemical convection in the Earth's outer core. EOS, 76, 280–281.
Jarvis G.T. and McKenzie D.P., 1980. Convection in a compressible fluid with infinite Prandtl number. J. Fluid Mech., 96, 515–583.
Jarvis G.T. and Peltier W.R., 1989. Convection models and geophysical observations. In: W.R. Peltier (Ed.), Mantle Convection-Plate Tectonics and Global Dynamics, Gordon and Breach Science Publishers, New York, 479–593.
Kuang W. and Bloxham J., 1997. An Earth-like numerical dynamo model. Nature, 389, 371–374.
Kuang W. and Bloxham J., 1999. Numerical modeling of magnetohydrodynamic convection in a rapidly rotating spherical shell: weak and strong field dynamo action. J. Comput. Phys., 153, 51–81.
Larsen T.B., Yuen D.A., Moser J. and Fornberg B., 1997. A high-order finite-difference method applied to large Rayleigh number mantle convection. Geophys. Astrophys. Fluid Dynamics, 84, 53–83.
Lister J.R. and Buffet B.A., 1995. The strength and efficiency of thermal and compositional convection in the geodynamo. Phys. Earth Planet. Inter., 91, 17–30.
Malevsky A.V., 1996. Spline-characteristic method for simulation of convective turbulence. J. Comput. Phys., 123, 466–475.
Malevsky A.V. and Thomas S.J., 1997. Parallel algorithms for semi-Lagrangian advection. J. Numer. Methods Fluids, 25, 455–473.
Merilees P.E., 1973. The pseudospectral approximation applied to the shallow water equations on a sphere. Atmosphere, 11, 13–20.
Merilees P.E., 1974. Numerical experiments with the pseudospectral method in spherical coordinates. Atmosphere, 12, 77–96.
Moser J., 1994. Mantle Dynamics and Rotation of the Earth. Ph.D. Thesis, Faculty of Mathematics and Physics, Charles University, Prague.
Pedlosky J., 1987. Geophysical Fluid Dynamics, Springer-Verlag, pp. 43–50.
Roberts P.H., 1992. Dynamo theory. In: D.A. Yuen (Ed.), Chaotic Processes in the Geological Sciences, The IMA Volumes in Mathematics and its Applications, Vol. 41, Springer-Verlag, New York, pp. 237–280.
St. Pierre M.G., 1996. On the local nature of turbulence in Earth's outer core. Geophys. Astrophys. Fluid Dynamics, 83, 293–306.
Sun Z.P. and Schubert G., 1995. Numerical simulations of thermal convection in a rotating spherical fluid at high Taylor and Rayleigh numbers. Phys. Fluids, 7, 2686–2699.
Tritton D.J., 1988. Physical Fluid Dynamics, 2nd edn., Oxford Univ. Press, New York.
VelĂmskĂ˝ J. and Matyska C., 2000. The influence of adiabatic heating/cooling on magnetohydrodynamic systems. Phys. Earth Planet. Inter., 117, 197–207.
Vitásek E., 1987. Numerické Metody, SNTL, Praha, 136–139 (in Czech).
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Mistr, Z., Matyska, C. & Yuen, D.A. Modeling of Rapidly Rotating Thermal Convection Using Vorticity and Vector Potential. Studia Geophysica et Geodaetica 46, 59–81 (2002). https://doi.org/10.1023/A:1019837616461
Published:
Issue Date:
DOI: https://doi.org/10.1023/A:1019837616461