INTRODUCTION

Numerical computations comprise an increasingly important tool in the understanding of the Earth's dynamo, and, with the increased accessibility of supercomputers, direct, realistic simulations of the geodynamo are not far off. Any such simulation must solve the equations governing a three dimensional, rapidly rotating, dynamically consistent dynamo with Lorentz force J × B present in the dominant balance of forces. The simplest dynamo with these characteristics, first proposed by Childress & Soward (1972), uses the convective motions of rapidly rotating Benard convection to drive a dynamically consistent MHD dynamo. Computationally, the Childress–Soward dynamo has the advantage of permitting the expansion of the unknown fields in Fourier series in all directions, allowing three dimensional fast Fourier transforms (FFT's) to be used in calculating the nonlinear terms. Since no fast Legendre transform exists at the moment, the resulting programs will be faster than more realistic spherical dynamo simulations, while at the same time reflecting the important features of these models.

In this paper, the strong field branch of the Childress–Soward dynamo is investigated using direct numerical simulations.