Solenoidal force balances in numerical dynamos

Numerical simulations of the geodynamo (and other planetary dynamos) have made significant progress in recent years. As computing power has advanced, some new models claim to be ever more appropriate for understanding Earth's core dynamics. One measure of the success of such models is the ability to replicate the expected balance between forces operating within the core; Coriolis and Lorentz forces are predicted to be most important. The picture is complicated for an incompressible flow by the existence of the pressure gradient force which renders the gradient parts of all other forces dynamically unimportant. This can confuse the situation, especially when the scale dependence of forces is considered. In this work we investigate force balances through the alternative approach of eliminating gradient parts of each force to form ‘solenoidal force balances’. We perform a length-scale-dependent analysis for several spherical simulations and find that removal of gradient parts offers an alternative picture of the force balance compared with looking at traditional forces alone. Solenoidal force balances provide some agreement with the results of previous studies but also significant differences. They offer a cleaner overall picture of the dynamics and introduce differences at smaller scales. This has implications for geodynamo models purporting to have reached Earth-like regimes: in order to achieve a meaningful comparison of forces, only the solenoidal part of forces should be considered.