Simulations of MHD Instabilities in Protoplanetary Disks

We have developed a 3D spectral, anelastic, magnetohydrodynamic code for rapidly rotating, strongly sheared, stratified protoplanetary disks. The Cartesian domain is resolved with a Fourier-Fourier-Chebyshev basis, and uses horizontal coordinates that advect with the background shear. With this tool, we simulate the magnetorotational instability at high resolution and investigate the development of MHD turbulence.