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 . 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 . 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 . 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 . 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 . 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.

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is about the influence of numerical diffusion on the accretion disk dynamo effect in magnetohydrodynamics (MHD) simulations. Specifically, the paper focuses on the fact that in global MHD simulations of accretion disks, the numerical magnetic diffusion may be too strong, thus limiting the enhancement of the magnetic field and resulting in the magnetic field strength in the simulation results being far lower than the actually observed values. The paper points out that both the existing local simulations (such as the shearing box method) and global simulations have failed to generate sufficient magnetic field strength to explain the observational data, especially in the highly ionized accretion disks of dwarf novae. Observations show that the plasma β values in these disks are close to 1 (that is, the magnetic field pressure is comparable to the gas pressure), while the β values in the simulation results are much greater than 1. To explore this problem, the paper first reviews the basic physical properties of accretion disks and their observational characteristics, especially the hot, highly ionized accretion disks in dwarf novae. Then, the paper discusses the application of numerical simulations in studying the accretion disk dynamo effect and points out the limitations of local and global simulations in explaining observational data. Finally, the paper proposes a hypothesis that numerical magnetic diffusion may be an important reason for the discrepancy between the simulation results and the observations, and illustrates this through theoretical calculations. The paper calculates the time evolution of a weak magnetic field in an incompressible fluid in the presence of linear shear and magnetic diffusion and finds that the maximum increase factor of the initial magnetic field energy depends on \( R_m^{2/3} \) of the magnetic Reynolds number \( R_m \). Based on this, the paper estimates that the numerical magnetic Reynolds number of the current global numerical simulations is about six orders of magnitude lower than the physical value in the accretion disks of dwarf novae, which may cause the simulations to be unable to truly reproduce the dynamo effect in the observable accretion disks.

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

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

Helical and nonhelical large-scale dynamos in thin accretion discs

Gravito-turbulence and dynamo in poorly ionised protostellar discs. I. Zero-net-flux case

Gravitoturbulent dynamos in astrophysical discs

Accretion Disks and Dynamos: Toward a Unified Mean Field Theory

Vorticity and magnetic dynamo from subsonic expansion waves II: Dependence on magnetic Prandtl number, forcing scale, cooling time

Gravitoturbulent dynamo in global simulations of gaseous disks

Helical dynamo growth and saturation at modest versus extreme magnetic Reynolds numbers

Vorticity and magnetic dynamo from subsonic expansion waves. II. Dependence on the magnetic Prandtl number, forcing scale, and cooling time

The role of magnetic boundaries in kinematic and self-consistent magnetohydrodynamic simulations of precession-driven dynamo action in a closed cylinder

On the dynamo driven accretion disks

Global simulations of self-gravitating magnetized protoplanetary disks

Magnetically threaded accretion disks in resistive magnetohydrodynamic simulations and asymptotic expansion

Non-linear magnetic buoyancy instability and galactic dynamos

Large scale magnetic fields in viscous resistive accretion disks. I. Ejection from weakly magnetized disks

Extracting scaling laws from numerical dynamo models

Beyond the standard accretion disc model: coupled magnetic disc--corona solutions with a physically motivated viscosity law

Magnetic helicity dissipation and production in an ideal MHD code

Transition from small-scale to large-scale dynamo in a supernova-driven, multiphase medium