Simplified Turbulence in a 2D Planetary Core

The magnetic field of the Earth and planets is created by the movement of an electrically conducting fluid in the deep planetary interior. In the case of the Earth, this is liquid iron in the fluid outer core. We can learn about how planets generate this magnetic field by numerically simulating this fluid motion with supercomputers. Unfortunately, even with the largest clusters we are still very far from accurately simulating the dynamics of planetary cores.

One way around this is to make certain simplifications. For example, fluid motion in rotating fluid systems can be thought of as quasi-two dimensional. If we simulate a two-dimensional, instead of a three-dimensional fluid, we can get closer to the correct parameter regime. To this end, I wrote a numerical model in Python that simulates an incompressible fluid in a rotating, two-dimensional shell. This would correspond to the fluid outer core of a planet. Below is a video of a simulation.

This is a movie showing vorticity in a turbulent, 2D, rotating convective system. This is meant to approximate flow in the Earth's core, this represents a slice through the equator. The outer edge of the simulation corresponds to the core-mantle boundary and the white area in the middle is the solid inner core. In this plot, blue patches correspond to fluid parcels that are rotating counterclockwise and red patches correspond to fluid parcels that are rotating clockwise. The units are non-dimensional.