There has been significant progress in the development of numerical geodyna
mo models over the last five years. Advances in computer technology have ma
de it possible to perform three-dimensional simulations, with thermal or co
mpositional convection as the driving mechanism. These numerical simulation
s give reasonable results for the morphology and strength of the field at t
he core-mantle boundary and the models are also capable of giving reversals
and excursions which can be compared with palaeomagnetic observations; the
y also predict differential rotation between the inner core and the mantle.
However, there are still a number of fundamental problems associated with t
he simulations, which are proving hard to overcome. Despite tl:le advances
in computing power, the models are still expensive and take a long time to
run. This problem may diminish as faster machines become available, and new
numerical methods exploit parallelization effectively, but currently there
ase no practical schemes available which work at low Ekman number.
Even with turbulent values of the diffusivities (and the question of whethe
r isotropic diffusivities are appropriate is still unresolved), the appropr
iate dynamical regime has not yet been reached. In consequence, modelling a
ssumptions about the nature of the flow near title boundaries have to be ma
de, and different choices can have profound effects on the dynamics. The na
ture of large-scale magnetoconvection at small E is still not well understo
od, and until we have more understanding of this issue, it will be difficul
t to have a great deal of confidence in the predictions of the numerical mo
dels.