Exact shearing box solutions of magnetohydrodynamic flows with resistivity, viscosity and cooling |
| |
Authors: | Pierre Lesaffre Steven A Balbus † |
| |
Institution: | Laboratoire de Radioastronomie, 24 rue Lhomond, 75231 Paris Cedex 05, France;Department of Astronomy, University of Virginia, Charlottesville, VA 22903, USA |
| |
Abstract: | Axisymmetric incompressible modes of the magnetorotational instability (MRI) with a vertical wavenumber are exact solutions of the non-linear local equations of motion for a disc (shearing box). They are referred to as 'channel solutions'. Here, we generalize a class of these solutions to include energy losses, viscous, and resistive effects. In the limit of zero shear, we recover the result that torsional Alfvén waves are exact solutions of the non-linear equations. Our method allows the extension of these solutions into the dissipative regime. These new solutions serve as benchmarks for simulations including dissipation and energy loss, and to calibrate numerical viscosity and resistivity in the zeus3d code. We quantify the anisotropy of numerical dissipation and compute its scaling with time and space resolution. We find a strong dependence of the dissipation on the mean magnetic field that may affect the saturation state of the MRI as computed with zeus3d . It is also shown that elongated grid cells generally preclude isotropic dissipation and that a Courant time-step smaller than that which is commonly used should be taken to avoid spurious anti-diffusion of magnetic field. |
| |
Keywords: | accretion accretion discs instabilities MHD waves methods: analytical methods: numerical |
|
|