A hydrostatic MHD code for modeling stellar interiors

¹ Argelander Institut für Astronomie, University of Bonn, Auf dem Hügel 71, 53121 Bonn, Germany
e-mail: jmitchell@astro.uni-bonn.de
² Instituto de Astrofísica, Facultad de Física, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, 7820436 Macul, Santiago, Chile

Received: 29 May 2015
Accepted: 14 October 2015

Abstract

In gravitationally stratified fluids, the length scales of fluid motion are often much greater in the horizontal direction than in the vertical. When modeling these fluids, it can be advantageous to use the hydrostatic approximation often used in atmospheric physics, which filters out vertically propagating sound waves and thus allows a longer time step. We describe a finite-difference numerical scheme that uses this approximation. This code is suitable for modeling (magneto)hydrodynamic processes in radiative stellar interiors, such as the baroclinic instability and stratified turbulence, or the Tayler-Spruit dynamo.