On the role of numerical diffusivity in MHD simulations of global accretion disc dynamos

Abstract

Observations, mainly of outbursts in dwarf novae, imply that the anomalous viscosity in highly ionized accretion discs is magnetic in origin and requires that the plasma ${\beta \sim 1}$ . Until now, most simulations of the magnetic dynamo in accretion discs have used a local approximation (known as the shearing box). While these simulations demonstrate the possibility of a self-sustaining dynamo, the magnetic activity generated in these models saturates at $\beta \gg 1$ . This long-standing discrepancy has previously been attributed to the local approximation itself. There have been recent attempts at simulating magnetic activity in global accretion discs with parameters relevant to the dwarf novae. These too find values of $\beta \gg 1$ . We speculate that the tension between these models and the observations may be caused by numerical magnetic diffusivity. As a pedagogical example, we present exact time-dependent solutions for the evolution of weak magnetic fields in an incompressible fluid subject to linear shear and magnetic diffusivity. We find that the maximum factor by which the initial magnetic energy can be increased depends on the magnetic Reynolds number as ${\mathcal {R}}_{m}^{2/3}$ . We estimate that current global numerical simulations of dwarf nova discs have numerical magnetic Reynolds numbers around six orders of magnitude less than the physical value found in dwarf nova discs of ${\mathcal {R}}_{m} \sim 10^{10}$ . We suggest that, given the current limitations on computing power, expecting to be able to compute realistic dynamo action in observable accretion discs using numerical MHD is, for the time being, a step too far.