Geodynamo (ryadav)¶
Overview¶
Simulations for the geodynamo by Rakesh Yadav (Harvard University). Simulations are thermally driven magnetohydrodynamics (MHD) simulations in a rotating spherical shell modeled on the Earth’s outer core. In the present simulations, the ratio of the inner to outer radii is set to 0.35, and the spherical shell is filled by a electrically conductive fluid. The thermally driven convection is assumed with the Boussinesq approximation. The inner core is assumed to be corotating with the mantle (outside of the spherical shell).
The nonslip boundary condition is applied to the inner and outer boundaries for the velocity, and temperature is fixed to be 1.0 and 0.0 at the inner and outer boundaries, respectively. The inner core and mantle is assumed to be electrically insulated. Consequently, the magnetic field is connected to the potential magnetic fields at the inner and outer boundaries.
Spatial resolution, dimensionless number, and hyperdiffusivity settings for each case are listed in the following table as a parameter setting in ‘main_input’. The definition of each parameter is described in the Rayleigh document.
Publications¶
The data is used for the following paper:
Cao, H., and Yadav, R. K., and Aurnou, J. M., Geomagnetic polar minima do not arise from steady meridional circulation, Proceedings of the National Academy of Sciences, 115, 11186–11191, 2018.
List of parameters¶
Run_name 
\(N_{r}\) 
\(N_{\theta}\) 
\(E\) 
\(Pr\) 
\(Pm\) 
\(Ra\) 
Hyperdiffusion 
\(\alpha\) 
\(\beta\) 

512 
1536 
1.00E06 
1.00E+00 
5.00E02 
2.50E+10 
.TRUE. 
3.00E+00 
7.00E+00 

512 
768 
1.00E06 
1.00E+00 
1.00E01 
1.00E+10 
.TRUE. 
4.00E+00 
7.00E+00 

512 
768 
1.00E06 
1.00E+00 
2.00E01 
5.00E+09 
.FALSE. 
5.00E+00 
7.00E+00 

256 
768 
1.00E06 
1.00E01 
5.00E02 
2.50E+09 
.TRUE. 
4.00E+00 
7.00E+00 

256 
768 
1.00E06 
1.00E+00 
5.00E01 
3.00E+09 
.FALSE. 
0.00E+00 
0.00E+00 

256 
768 
1.00E06 
1.00E+00 
1.00E+00 
3.00E+09 
.TRUE. 
3.00E+00 
6.00E+00 

512 
1536 
1.00E07 
1.00E+00 
3.00E02 
3.00E+11 
.TRUE. 
8.00E+00 
5.00E+00 

512 
1536 
1.00E07 
1.00E+00 
5.00E02 
1.50E+11 
.TRUE. 
3.00E+00 
1.00E+01 

512 
1536 
1.00E07 
1.00E+00 
2.00E01 
6.00E+10 
.TRUE. 
2.00E+00 
9.00E+00 
Definition of parameters¶
Geometry and spatial resolution
\(N_{r}\): Number of radial grid points
\(N_{\theta}\): Number of meridional grid points
Dimensionless numbers
\(E = \frac{\nu}{\Omega L^2}\): Ekman number
\(Pr = \frac{\kappa}{\nu}\): Prandtl number
\(Pm = \frac{\eta}{\nu}\): Magnetic Prandtl number
\(Ra = \frac{\alpha g_{o} T L^3}{\kappa \nu}\): Rayleigh number
 Hyperdiffusivity parameters
Hyperdiffusion: Flag for hyperdiffusivity
 Set this to variable to .true. to enable hyperdiffusion. The default value is .false. When active, diffusivities \(\nu\), \(\kappa\), and \(\eta\), are multiplied by an additional factor as:
\(\{\nu,\kappa,\eta\}\rightarrow\{\nu,\kappa,\eta\}\left(1+\alpha\left(\frac{\ell1}{\ell_\mathrm{max}1}\right)^\beta\right)\)
\(\alpha\): The value of \(\alpha\) when hyper diffusion is active.
\(\beta\): The value of \(\beta\) when hyper diffusion is active.