ANZIAM J. 46(E) pp.C854--C870, 2005.

An angular spectral method for solution of the heat equation in spheroidal geometries

D. J. Ivers
(Received 28 October 2004, revised 24 June 2005)

Abstract

A spectral numerical method is presented for solving the heat equation in oblate or prolate spheroids. Cartesian coordinates are scaled to transform the spheroidal geometry into a spherical geometry. The diffusion term in the transformed equation is anisotropic, being enhanced in the polar directions. The transformed equation is discretised in angle using a truncated spherical harmonic expansion of the temperature in transformed spherical polar coordinates. The anisotropic diffusion term is reduced to block tridiagonal form using recurrence relations for spherical harmonics. The radial coordinate is discretised using finite differences in scaled radius but other radial schemes are possible. Without heat sources and with a homogeneous Dirichlet boundary condition the problem reduces to an eigenproblem for the decay rate. The results are compared to the separated variables solution, which employs oblate or prolate spheroidal wave functions. The method directly extends to other scalar problems in spheroidal geometries, which have the highest derivatives in Laplacian form, including passive advection-diffusion of a scalar. The method may be extended, with difficulty, to problems with vector diffusion or ellipsoidal geometries.

Download to your computer

Authors

D. J. Ivers
School of Mathematics and Statistics, University of Sydney, Sydney, Australia. mailto:D.Ivers@maths.usyd.edu.au

Published September 1, 2005. ISSN 1446-8735

References

  1. S. I. Braginsky and V. P. Meytlis. Local turbulence in the Earth's core. Geophys. Astrophys. Fluid Dynam., 55:71--87, 1990.
  2. C. Flammer. Spheroidal Wave Functions, Stanford University Press, Stanford, California. 1957
  3. D. J. Ivers. Geodynamo models with anisotropic turbulent thermal diffusion. Abstracts of the 9th SEDI (Studies of the Earth's Deep Interior) Symposium, page 41 Garmisch--Partenkirchen, Germany. 2004
  4. D. J. Ivers. Thermal instability of an oblate spheroid. Abstracts of the 8th SEDI (Studies of the Earth's Deep Interior) Symposium, Granlibakken, USA. 2002
  5. C. Niven. On the conduction of heat in ellipsoids of revolution. Phil. Trans. Roy. Soc. Lond., 171:117--151, 1880.
  6. G. C. Phillips, and D. J. Ivers. Spherical anisotropic diffusion models for the Earth's core. Phys. Earth Planet. Int., 117:209--223, 2000.
  7. G. C. Phillips and D. J. Ivers. Spectral interactions of rapidly-rotating anisotropic turbulent viscous and thermal diffusion in the Earth's core. Phys. Earth Planet. Int., 128:93--107, 2001.
  8. G. C. Phillips and D. J. Ivers. Strong field anisotropic diffusion models for the Earth's core. Phys. Earth Planet. Int., 140:13--28, 2003.
Copyright Austral. Mathematical Soc.; prepared September 1, 2005, from Iver.tex with anziamje.sty (v2.3h)