Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userZonal Flows and Long-Lived Axisymmetric Pressure Bumps in Magnetorotational Turbulence
112
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We study the behavior of magnetorotational turbulence in shearing box simulations with a radial and azimuthal extent up to ten scale heights. Maxwell and Reynolds stresses are found to increase by more than a factor two when increasing the box size beyond two scale heights in the radial direction. Further increase of the box size has little or no effect on the statistical properties of the turbulence. An inverse cascade excites magnetic field structures at the largest scales of the box. The corresponding 10% variation in the Maxwell stress launches a zonal flow of alternating sub- and super-Keplerian velocity. This in turn generates a banded density structure in geostrophic balance between pressure and Coriolis forces. We present a simplified model for the appearance of zonal flows, in which stochastic forcing by the magnetic tension on short time-scales creates zonal flow structures with life-times of several tens of orbits. We experiment with various improved shearing box algorithms to reduce the numerical diffusivity introduced by the supersonic shear flow. While a standard finite difference advection scheme shows signs of a suppression of turbulent activity near the edges of the box, this problem is eliminated by a new method where the Keplerian shear advection is advanced in time by interpolation in Fourier space.
rate research
Read More
Accretion disks are likely threaded by external vertical magnetic flux, which enhances the level of turbulence via the magnetorotational instability (MRI). Using shearing-box simulations, we find that such external magnetic flux also strongly enhances the amplitude of banded radial density variations known as zonal flows. Moreover, we report that vertical magnetic flux is strongly concentrated toward low-density regions of the zonal flow. Mean vertical magnetic field can be more than doubled in low-density regions, and reduced to nearly zero in high density regions in some cases. In ideal MHD, the scale on which magnetic flux concentrates can reach a few disk scale heights. In the non-ideal MHD regime with strong ambipolar diffusion, magnetic flux is concentrated into thin axisymmetric shells at some enhanced level, whose size is typically less than half a scale height. We show that magnetic flux concentration is closely related to the fact that the magnetic diffusivity of the MRI turbulence is anisotropic. In addition to a conventional Ohmic-like turbulent resistivity, we find that there is a correlation between the vertical velocity and horizontal magnetic field fluctuations that produces a mean electric field that acts to anti-diffuse the vertical magnetic flux. The anisotropic turbulent diffusivity has analogies to the Hall effect, and may have important implications for magnetic flux transport in accretion disks. The physical origin of magnetic flux concentration may be related to the development of channel flows followed by magnetic reconnection, which acts to decrease the mass-to-flux ratio in localized regions. The association of enhanced zonal flows with magnetic flux concentration may lead to global pressure bumps in protoplanetary disks that helps trap dust particles and facilitates planet formation.
We reveal and investigate a new type of linear axisymmetric helical magnetorotational instability which is capable of destabilizing viscous and resistive rotational flows with radially increasing angular velocity, or positive shear. This instability is double-diffusive by nature and is different from the more familiar helical magnetorotational instability, operating at positive shear above the Liu limit, in that it works instead for a wide range of the positive shear when ${rm (i)}$ a combination of axial/poloidal and azimuthal/toroidal magnetic fields is applied and ${rm (ii)}$ the magnetic Prandtl number is not too close to unity. We study this instability first with radially local WKB analysis and then confirm its existence using a global stability analysis of the magnetized flow between two rotating cylinders with conducting or insulating boundaries. From an experimental point of view, we also demonstrate the presence of the new instability in a magnetized viscous and resistive Taylor-Couette flow with positive shear for such values of the flow parameters, which can be realized in upcoming experiments at the DRESDYN facility. Finally, this instability might have implications for the dynamics of the equatorial parts of the solar tachocline and dynamo action there, since the above two necessary conditions for the instability to take place are satisfied in this region. Our global stability calculations for the tachocline-like configuration, representing a thin rotating cylindrical layer with the appropriate boundary conditions -- conducting inner and insulating outer cylinders -- and the values of the flow parameters, indicate that it can indeed arise in this case with a characteristic growth time comparable to the solar cycle period.
Axisymmetric magnetorotational instability (MRI) in viscous accretion disks is investigated by linear analysis and two-dimensional nonlinear simulations. The linear growth of the viscous MRI is characterized by the Reynolds number defined as $R_{rm MRI} equiv v_A^2/ uOmega $, where $v_A$ is the Alfv{e}n velocity, $ u$ is the kinematic viscosity, and $Omega$ is the angular velocity of the disk. Although the linear growth rate is suppressed considerably as the Reynolds number decreases, the nonlinear behavior is found to be almost independent of $R_{rm MRI}$. At the nonlinear evolutionary stage, a two-channel flow continues growing and the Maxwell stress increases until the end of calculations even though the Reynolds number is much smaller than unity. A large portion of the injected energy to the system is converted to the magnetic energy. The gain rate of the thermal energy, on the other hand, is found to be much larger than the viscous heating rate. Nonlinear behavior of the MRI in the viscous regime and its difference from that in the highly resistive regime can be explained schematically by using the characteristics of the linear dispersion relation. Applying our results to the case with both the viscosity and resistivity, it is anticipated that the critical value of the Lundquist number $S_{rm MRI} equiv v_A^2/etaOmega$ for active turbulence depends on the magnetic Prandtl number $S_{{rm MRI},c} propto Pm^{1/2}$ in the regime of $Pm gg 1$ and remains constant when $Pm ll 1$, where $Pm equiv S_{rm MRI}/R_{rm MRI} = u/eta$ and $eta$ is the magnetic diffusivity.
The search for young planets had its first breakthrough with the detection of the accreting planet PDS70b. In this study, we aim to broaden our understanding towards the formation of multi-planet systems such as HR8799 or the Solar System. Our previous study on HD169142, one of the closest Herbig stars, points towards a shadow-casting protoplanetary candidate. Here, we present follow-up observations to test our previously proposed hypothesis. We set our new data into context with previous observations to follow structural changes in the disk over the course of 6 years. We find spatially resolved systematic changes in the position of the previously described surface brightness dip in the inner ring. We further find changes in the brightness structure in azimuthal direction along the ring. And finally, a comparison of our SPHERE data with recent ALMA observations reveals a wavelength dependent radial profile of the bright ring. The time-scale on which the changes in the rings surface brightness occur suggest that they are caused by a shadow cast by a 1-10Mj planet surrounded by dust, an orbit comparable to those of the giant planets in our own Solar System. Additionally, we find the first indications for temperature-induced instabilities in the ring. And finally, we trace a pressure maxima, for the first time spatially resolved, with a width of 4.5au. The density distribution of the ring at mm wavelengths around the pressure maxima could further indicate effects from snow lines or even the dynamics and feedback of the larger grains.
We present results from the first 3D kinetic numerical simulation of magnetorotational turbulence and dynamo, using the local shearing-box model of a collisionless accretion disc. The kinetic magnetorotational instability grows from a subthermal magnetic field having zero net flux over the computational domain to generate self-sustained turbulence and outward angular-momentum transport. Significant Maxwell and Reynolds stresses are accompanied by comparable viscous stresses produced by field-aligned ion pressure anisotropy, which is regulated primarily by the mirror and ion-cyclotron instabilities through particle trapping and pitch-angle scattering. The latter endow the plasma with an effective viscosity that is biased with respect to the magnetic-field direction and spatio-temporally variable. Energy spectra suggest an Alfven-wave cascade at large scales and a kinetic-Alfven-wave cascade at small scales, with strong small-scale density fluctuations and weak non-axisymmetric density waves. Ions undergo non-thermal particle acceleration, their distribution accurately described by a kappa distribution. These results have implications for the properties of low-collisionality accretion flows, such as that near the black hole at the Galactic center.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
