Large-scale vertical magnetic fields are believed to play a key role in the evolution of protoplanetary discs. Associated with non-ideal effects, such as ambipolar diffusion, they are known to launch a wind that could drive accretion in the outer part of the disc ($R> 1$ AU). They also potentially lead to self-organisation of the disc into large-scale axisymmetric structures, similar to the rings recently imaged by sub-millimetre or near-infrared instruments (ALMA and SPHERE). The aim of this paper is to investigate the mechanism behind the formation of these gaseous rings, but also to understand the dust dynamics and its emission in discs threaded by a large-scale magnetic field. To this end, we performed global magneto-hydrodynamics (MHD) axisymmetric simulations with ambipolar diffusion using a modified version of the PLUTO code. We explored different magnetisations with the midplane $\beta$ parameter ranging from $10^5$ to $10^3$ and included dust grains -- treated in the fluid approximation -- ranging from $100 \mu$m to 1 cm in size. We first show that the gaseous rings (associated with zonal flows) are tightly linked to the existence of MHD winds. Secondly, we find that millimetre-size dust is highly sedimented, with a typical scale height of 1 AU at $R=100$ AU for $\beta=10^4$, compatible with recent ALMA observations. We also show that these grains concentrate into pressure maxima associated with zonal flows, leading to the formation of dusty rings. Using the radiative transfer code MCFOST, we computed the dust emission and make predictions on the ring-gap contrast and the spectral index that one might observe with interferometers like ALMA.

Gravitational instability (GI) controls the dynamics of young massive protoplanetary discs. Apart from facilitating gas accretion on to the central protostar, it must also impact on the process of planet formation: directly through fragmentation, and indirectly through the turbulent concentration of small solids. To understand the latter process, it is essential to determine the dust dynamics in such a turbulent flow. For that purpose, we conduct a series of 3D shearing box simulations of coupled gas and dust, including the gas's self-gravity and scanning a range of Stokes numbers, from 0.001 to ~0.2. First, we show that the vertical settling of dust in the midplane is significantly impeded by gravitoturbulence, with the dust scale-height roughly 0.6 times the gas scale height for centimetre grains. This is a result of the strong vertical diffusion issuing from (a) small-scale inertial-wave turbulence feeding off the GI spiral waves and (b) the larger-scale vertical circulations that naturally accompany the spirals. Second, we show that at R=50 AU concentration events involving sub-metre particles and yielding order 1 dust to gas ratios are rare and last for less than an orbit. Moreover, dust concentration is less efficient in 3D than in 2D simulations. We conclude that GI is not especially prone to the turbulent accumulation of dust grains. Finally, the large dust scale-height measured in simulations could be, in the future, compared with that of edge-on discs seen by ALMA, thus aiding detection and characterisation of GI in real systems.

Semi-analytical models of disk outflows have successfully described magnetically-driven, self-confined super-Alfv\'enic jets from near Keplerian accretion disks. These Jet Emitting Disks are possible for high levels of disk magnetization $\mu$ defined as $\mu=2/\beta$ where beta is the usual plasma parameter. In near-equipartition JEDs, accretion is supersonic and jets carry away most of the disk angular momentum. However, these solutions prove difficult to compare with cutting edge numerical simulations, for the reason that numerical simulations show wind-like outflows but in the domain of small magnetization. In this work, we present for the first time self-similar isothermal solutions for accretion-ejection structures at small magnetization levels. We elucidate the role of MRI-like structures in the acceleration processes that drive this new class of solutions. The disk magnetization $\mu$ is the main control parameter: massive outflows driven by the pressure of the toroidal magnetic field are obtained up to $\mu \sim 10^{-2}$, while more tenuous centrifugally-driven outflows are obtained at larger $\mu$ values. The generalized parameter space and the astrophysical consequences are discussed. We believe that these new solutions could be a stepping stone in understanding the way astrophysical disks drive either winds or jets. Defining jets as self-confined outflows and winds as uncollimated outflows, we propose a simple analytical criterion based on the initial energy content of the outflow, to discriminate jets from winds. We show that jet solution are achieved at all magnetization level, while winds could be obtained only in weakly magnetized disks that feature heating.

Dwarf novae (DNe) and X-ray binaries exhibit outbursts thought to be due to a thermal-viscous instability in the accretion disk. The disk instability model (DIM) assumes that accretion is driven by turbulent transport, customarily attributed to the magneto-rotational instability (MRI). Recent results point out that MRI turbulence alone fails to reproduce the light curves of DNe. We aim to study the impact of wind-driven accretion on the light curves of DNe. Local and global simulations show that magneto-hydrodynamic winds are present when a magnetic field threads the disk, even for relatively high ratios of thermal pressure to magnetic pressure ($\beta \approx 10^{5}$). These winds are very efficient in removing angular momentum but do not heat the disk; they do not behave as MRI-driven turbulence. We add wind-driven transport in the angular momentum equation of the DIM, assuming a fixed magnetic configuration: dipolar or constant with radius. We use prescriptions for the wind torque and the turbulent torque derived from shearing box simulations. The wind torque enhances the accretion of matter, resulting in light curves that look like DNe outbursts when assuming a dipolar field with a moment $\mu\approx10^{30}\,\mathrm{G\,cm^{3}}$. In the region where the wind dominates, the disk is cold, optically thin and the accretion speed is sonic. This acts as if the inner disk was truncated, leading to higher quiescent X-ray luminosities from the white dwarf boundary layer than expected with the standard DIM. The disk is stabilized if the wind-dominated region is large enough, potentially leading to `dark' disks emitting little radiation. Wind-driven accretion can play a key role in shaping the light curves of DNe and X-ray binaries. Future studies will need to include the time evolution of the magnetic field threading the disk to fully assess its impact on the dynamics of the accretion flow.

Rings and gaps have been observed in a wide range of protoplanetary discs, from young systems like HLTau to older discs like TW Hydra. Recent disc simulations have shown that magnetohydrodynamic (MHD) turbulence (in the ideal or non-ideal regime) can lead to the formation of rings and be an alternative to the embedded planets scenario. In this paper, we investigate how these ring form in this context and seek a generic formation process, taking into account the various dissipative regimes and magnetizations probed by the past simulations. We identify the existence of a linear and secular instability, driven by MHD winds, and giving birth to rings of gas having a width larger than the disc scale height. We show that the linear theory is able to make reliable predictions regarding the growth rates, ring/gap contrast and spacing, by comparing these predictions to a series of 2D (axisymmetric) and 3D MHD numerical simulations. In addition, we demonstrate that these rings can act as dust traps provided that the disc is sufficiently magnetised, with plasma beta lower than $10^4$. Given its robustness, the process identified in this paper could have important implications, not only for protoplanetary discs but also for a wide range of accreting systems threaded by large-scale magnetic fields.

Context. Disks around pre-main-sequence stars evolve over time by turbulent viscous spreading. The main contender to explain the strength of the turbulence is the Magneto-Rotational-Instability (MRI) model, whose efficiency depends on the disk ionization fraction. Aims. We aim at computing self-consistently the chemistry including PAH charge chemistry, the grain charging and an estimate of an effective value of the turbulence alpha parameter in order to find observational signatures of disk turbulence. Methods. We introduced PAH and grain charging physics and their interplay with other gas-phase reactions in the physico-chemical code ProDiMo. Non-ideal magnetohydrodynamics effects such as Ohmic and ambipolar diffusion are parametrized to derive an effective value for the turbulent parameter alpha_eff . We explored the effects of turbulence heating and line broadening on CO isotopologue sub-millimeter lines. Results. The spatial distribution of alpha_eff depends on various unconstrained disk parameters such as the magnetic parameter beta_mag or the cosmic ray density distribution inside the protoplanetary disks. The inner disk midplane shows the presence of the so-called dead-zone where the turbulence is quasi-inexistent. The disk is heated mostly by thermal accommodation on dust grains in the dead-zone, by viscous heating outside the dead-zone up to a few hundred astronomical units, and by chemical heating in the outer disk. The CO rotational lines probe the warm molecular disk layers where the turbulence is at its maximum. However, the effect of turbulence on the CO line profiles is minimal and difficult to distinguish from the thermal broadening. Conclusions. Viscous heating of the gas in the disk midplane outside the dead-zone is efficient. The determination of alpha from CO rotational line observations alone is challenging.

Dwarf novae (DNe) are accreting white dwarfs that show eruptions due to a thermal-viscous instability in the accretion disk. The outburst timescales constrain $\alpha$, the ratio of the viscous stress to the thermal pressure, and so the mechanism of angular momentum transport. The eruptive state has $\alpha\approx0.1$ while the quiescent state has $\alpha\approx0.03$. Turbulent transport due to the magneto-rotational instability (MRI) is generally considered to be the source of angular momentum transport in DNe. Here, we perform 3D local magnetohydrodynamic (MHD) shearing box simulations including vertical stratification, radiative transfer and a net constant vertical magnetic flux to investigate how transport changes between the outburst and quiescent states of DNe. We find that a constant $B_z$ provides a higher $\alpha$ in quiescence than in outburst, in opposition to what is expected. Including resistivity quenches MRI turbulence in quiescence, suppressing transport, unless the magnetic field is high enough, which again leads to $\alpha\approx0.1$. A major difference between simulations with a net poloidal flux and simulations without is that angular momentum transport in the former is shared between turbulent and wind-driven transport. We find that wind-driven transport dominates in quiescence even for low magnetic fields $\sim 1$ G. This can have a major impact on observational signatures since wind-driven transport does not heat the disk. Furthermore, wind transport cannot be reduced to an $\alpha$ prescription. We provide fits for $\alpha$ and the wind torque with $\beta$, ratio of thermal to magnetic pressure. We conclude that the evolution of the thermal-viscous instability, and its consequences on the outburst cycles of CVs, needs to be seriously revised to take into account that most of the accretion energy may be carried away by a wind instead of being locally dissipated.

MHD turbulence plays a crucial role in the dust dynamics of protoplanetary discs. It affects planet formation, vertical settling and is one possible origin of the large scale axisymmetric structures, such as rings, recently imaged by ALMA and SPHERE. Among the variety of MHD processes, the magnetorotational instability (MRI) has raised particular interest since it provides a source of turbulence and potentially organizes the flow into radial structures. However, the weak ionization of discs prevents the MRI from being excited beyond 1 AU. The strong sedimentation of millimetre dust measured in T-Tauri discs is also in contradiction with predictions based on ideal MRI turbulence. In this paper, we study the effects of non-ideal MHD and winds on the dynamics and sedimentation of dust grains. We consider a weakly ionized plasma subject to ambipolar diffusion characterizing the disc outer regions (>1 AU). For that, we perform numerical MHD simulations in the stratified shearing box, using the PLUTO code. Our simulations show that the mm-cm dust is contained vertically in a very thin layer, with typical heightscale ~0.4 AU at 30 AU, compatible with recent ALMA observations. Horizontally, the grains are trapped within pressure maxima induced by ambipolar diffusion, leading to the formation of dust rings. For micrometer grains, dust and gas scaleheights are similar. In this regime, the settling cannot be explained by a simple 1D diffusion theory but results from a large scale 2D circulation induced by both MHD winds and zonal flows. Overall, our results show that non-ideal MHD effects and their related winds play a major role in shaping the radial and vertical distribution of dust in protoplanetary discs. Leading to substantial accretion efficiency, non-ideal effects also a promising avenue to reconcile the low turbulent activity measured in discs with their relatively high accretion rates.

The eruptive cycles of dwarf novae (DN) are thought to be due to a thermal-viscous instability in the accretion disk surrounding the white dwarf (WD). This model has long been known to imply a stress to pressure ratio \alpha ~0.1 in outburst compared to \alpha ~ 0.01 in quiescence. Such an enhancement in $\alpha$ has recently been observed in simulations of turbulent transport driven by the magneto-rotational instability (MRI) when convection is present, without requiring a net magnetic flux. We independently recover this result by carrying out PLUTO MHD simulations of vertically stratified, radiative, shearing boxes with the thermodynamics and opacities appropriate to DN. The results are robust against the choice of vertical boundary conditions. The thermal equilibrium solutions found by the simulations trace the well-known S-curve in the density-temperature plane. We confirm that the high values of \alpha ~ 0.1 occur near the tip of the hot branch of the S-curve, where convection is active. However, we also present thermally-stable simulations at lower temperatures that have standard values of \alpha ~ 0.03 despite the presence of vigorous convection. We find no simple relationship between \alpha and the strength of the convection, as measured by the ratio of convective to radiative flux. The cold branch is only very weakly ionized so, in the second part of this work, we studied the impact of non-ideal MHD effects on transport. We include resistivity in the simulations and find that the MRI-driven transport is quenched (\alpha ~ 0) below the critical density at which the magnetic Reynolds number R_m \leq 10^4. This is problematic as X-ray emission observed in quiescent systems requires ongoing accretion onto the WD. We verify that these X-rays cannot self-sustain MRI-driven turbulence by photo-ionizing the disk and discuss possible solutions to the issue of accretion in quiescence.

We examine the properties of strongly magnetized accretion discs in a global framework, with particular focus on the evolution of magnetohydrodynamic instabilities such as the magnetorotational instability (MRI). Work by Pessah and Psaltis showed that MRI is stabilized beyond a critical toroidal field in compressible, differentially rotating flows and, also, reported the appearance of two new instabilities beyond this field. Their results stemmed from considering geometric curvature effects due to the suprathermal background toroidal field, which had been previously ignored in weak-field studies. However, their calculations were performed under the local approximation, which poses the danger of introducing spurious behavior due to the introduction of global geometric terms in an otherwise local framework. In order to avoid this, we perform a global eigenvalue analysis of the linearized MHD equations in cylindrical geometry. We confirm that MRI indeed tends to be highly suppressed when the background toroidal field attains the Pessah-Psaltis limit. We also observe the appearance of two new instabilities that emerge in the presence of highly suprathermal toroidal fields. These results were additionally verified using numerical simulations in PLUTO. There are, however, certain differences between the the local and global results, especially in the vertical wavenumber occupancies of the various instabilities, which we discuss in detail. We also study the global eigenfunctions of the most unstable modes in the suprathermal regime, which are inaccessible in the local analysis. Overall, our findings emphasize the necessity of a global treatment for accurately modeling strongly magnetized accretion discs.

The radial transport of angular momentum in accretion disk is a fundamental process in the universe. It governs the dynamical evolution of accretion disks and has implications for various issues ranging from the formation of planets to the growth of supermassive black holes. While the importance of magnetic fields for this problem has long been demonstrated, the existence of a source of transport solely hydrodynamical in nature has proven more difficult to establish and to quantify. In recent years, a combination of results coming from experiments, theoretical work and numerical simulations has dramatically improved our understanding of hydrodynamically mediated angular momentum transport in accretion disk. Here, based on these recent developments, we review the hydrodynamical processes that might contribute to transporting angular momentum radially in accretion disks and highlight the many questions that are still to be answered.

The intermittent small-scale structure of turbulence governs energy dissipation in many astrophysical plasmas and is often believed to have universal properties for sufficiently large systems. In this work, we argue that small-scale turbulence in accretion disks is universal in the sense that it is insensitive to the magnetorotational instability (MRI) and background shear, and therefore indistinguishable from standard homogeneous magnetohydrodynamic (MHD) turbulence at small scales. We investigate the intermittency of current density, vorticity, and energy dissipation in numerical simulations of incompressible MHD turbulence driven by the MRI in a shearing box. We find that the simulations exhibit a similar degree of intermittency as in standard MHD turbulence. We perform a statistical analysis of intermittent dissipative structures and find that energy dissipation is concentrated in thin sheet-like structures that span a wide range of scales up to the box size. We show that these structures exhibit strikingly similar statistical properties to those in standard MHD turbulence. Additionally, the structures are oriented in the toroidal direction with a characteristic tilt of approximately 17.5 degrees, implying an effective guide field in that direction.

The planet-forming region of protoplanetary disks is cold, dense, and therefore weakly ionized. For this reason, magnetohydrodynamic (MHD) turbulence is thought to be mostly absent, and another mechanism has to be found to explain gas accretion. It has been proposed that magnetized winds, launched from the ionized disk surface, could drive accretion in the presence of a large-scale magnetic field. The efficiency and the impact of these surface winds on the disk structure is still highly uncertain. We present the first global simulations of a weakly ionized disk that exhibits large-scale magnetized winds. We also study the impact of self-organization, which was previously demonstrated only in non-stratified models. We perform numerical simulations of stratified disks with the PLUTO code. We compute the ionization fraction dynamically, and account for all three non-ideal MHD effects: ohmic and ambipolar diffusions, and the Hall drift. Simplified heating and cooling due to non-thermal radiation is also taken into account in the disk atmosphere. We find that disks can be accreting or not, depending on the configuration of the large-scale magnetic field. Magnetothermal winds, driven both by magnetic acceleration and heating of the atmosphere, are obtained in the accreting case. In some cases, these winds are asymmetric, ejecting predominantly on one side of the disk. The wind mass loss rate depends primarily on the average ratio of magnetic to thermal pressure in the disk midplane. The non-accreting case is characterized by a meridional circulation, with accretion layers at the disk surface and decretion in the midplane. Finally, we observe self-organization, resulting in axisymmetric rings of density and associated pressure "bumps". The underlying mechanism and its impact on observable structures are discussed.

Understanding how accretion proceeds in proto-planetary discs and more generally their dynamics is a crucial issue for explaining the conditions in which planets form. The role that accretion of gas from the surrounding molecular cloud onto the disc may have on its structure needs to be quantified. We perform tri-dimensional simulations using the Cartesian AMR code RAMSES of an accretion disc subject to infalling material. For the aspect ratio of $H/R \simeq 0.15$ and disk mass $M_d \simeq 10^{-2}$ M$_\odot$ used in our study, we find that for typical accretion rates on the order of a few 10$^{-7}$ M$_\odot$ yr$^{-1}$, values of the $\alpha$ parameter as high as a few 10$^{-3}$ are inferred. The mass that is accreted in the inner part of the disc is typically at least $50\%$ of the total mass that has been accreted onto the disc. Our results suggest that external accretion of gas at moderate values, onto circumstellar discs may trigger prominent spiral arms, reminiscent of recent observations made with various instruments, and lead to significant transport through the disc. If confirmed from observational studies, such accretion may therefore influence disc evolution.

Understanding the diversity of planets requires to study the morphology and the physical conditions in the protoplanetary disks in which they form. We observed and spatially resolved the disk around the ~10 Myr old protoplanetary disk HD 100453 in polarized scattered light with SPHERE/VLT at optical and near-infrared wavelengths, reaching an angular resolution of ~0.02", and an inner working angle of ~0.09". We detect polarized scattered light up to ~0.42" (~48 au) and detect a cavity, a rim with azimuthal brightness variations at an inclination of 38 degrees, two shadows and two symmetric spiral arms. The spiral arms originate near the location of the shadows, close to the semi major axis. We detect a faint spiral-like feature in the SW that can be interpreted as the scattering surface of the bottom side of the disk, if the disk is tidally truncated by the M-dwarf companion currently seen at a projected distance of ~119 au. We construct a radiative transfer model that accounts for the main characteristics of the features with an inner and outer disk misaligned by ~72 degrees. The azimuthal brightness variations along the rim are well reproduced with the scattering phase function of the model. While spirals can be triggered by the tidal interaction with the companion, the close proximity of the spirals to the shadows suggests that the shadows could also play a role. The change in stellar illumination along the rim, induces an azimuthal variation of the scale height that can contribute to the brightness variations. Dark regions in polarized images of transition disks are now detected in a handful of disks and often interpreted as shadows due to a misaligned inner disk. The origin of such a misalignment in HD 100453, and of the spirals, is unclear, and might be due to a yet-undetected massive companion inside the cavity, and on an inclined orbit.

In Keplerian accretion disks, turbulence and magnetic fields may be jointly excited through a subcritical dynamo process involving the magnetorotational instability (MRI). High-resolution simulations exhibit a tendency towards statistical self-organization of MRI dynamo turbulence into large-scale cyclic dynamics. Understanding the physical origin of these structures, and whether they can be sustained and transport angular momentum efficiently in astrophysical conditions, represents a significant theoretical challenge. The discovery of simple periodic nonlinear MRI dynamo solutions has recently proven useful in this respect, and has notably served to highlight the role of turbulent magnetic diffusion in the seeming decay of the dynamics at low magnetic Prandtl number Pm (magnetic diffusivity larger than viscosity), a common regime in accretion disks. The connection between these simple structures and the statistical organization reported in turbulent simulations remained elusive, though. Here, we report the numerical discovery in moderate aspect ratio Keplerian shearing boxes of new periodic, incompressible, three-dimensional nonlinear MRI dynamo solutions with a larger dynamical complexity reminiscent of such simulations. These "chimera" cycles are characterized by multiple MRI-unstable dynamical stages, but their basic physical principles of self-sustainment are nevertheless identical to those of simpler cycles found in azimuthally elongated boxes. In particular, we find that they are not sustained at low Pm either due to subcritical turbulent magnetic diffusion. These solutions offer a new perspective into the transition from laminar to turbulent instability-driven dynamos, and may prove useful to devise improved statistical models of turbulent accretion disk dynamos.

Recently it has been proposed that the zombie vortex instability (ZVI) could precipitate hydrodynamical activity and angular momentum transport in unmagnetised regions of protoplanetary discs, also known as "dead zones". In this letter we scrutinise, with high resolution 3D spectral simulations, the onset and survival of this instability in the presence of viscous and thermal physics. First, we find that the ZVI is strongly dependent on the nature of the viscous operator. Although the ZVI is easily obtained with hyper-diffusion, it is difficult to sustain with physical (second order) diffusion operators up to Reynolds numbers as high as 10^7. This sensitivity is probably due to the ZVI's reliance on critical layers, whose characteristic lengthscale, structure, and dynamics are controlled by viscous diffusion. Second, we observe that the ZVI is sensitive to radiative processes, and indeed only operates when the Peclet number is greater than a critical value ~10^4, or when the cooling time is longer than ~10 Omega^{-1}. As a consequence, the ZVI struggles to appear at R>0.3 AU in standard 0.01 M_sun T Tauri disc models, though younger more massive disks provide a more hospitable environment. Together these results question the prevalence of the ZVI in protoplanetary discs.

Recent observations revealed organised structures in protoplanetary disks, such as axisymmetric rings or horseshoe concen- trations evocative of large-scale vortices. These structures are often interpreted as the result of planet-disc interactions. However, these disks are also known to be unstable to the magneto-rotational instability (MRI) which is believed to be one of the dominant angular momentum transport mechanism in these objects. It is therefore natural to ask if the MRI itself could produce these structures without invoking planets. The nonlinear evolution of the MRI is strongly affected by the low ionisation fraction in protoplanetary disks. The Hall effect in particular, which is dominant in dense and weakly ionised parts of these objects, has been shown to spontaneously drive self- organising flows in shearing box simulations. Here, we investigate the behaviour of global MRI-unstable disc models dominated by the Hall effect and characterise their dynamics. We perform 3D unstratified Hall-MHD simulations of keplerian disks for a broad range of Hall, ohmic and ambipolar Elsasser numbers. We confirm the transition from a turbulent to an organised state as the intensity of the Hall effect is increased. We observe the formation of zonal flows, their number depending on the available magnetic flux and on the intensity of the Hall effect. For intermediate Hall intensity, the flow self-organises into long-lived magnetised vortices. Neither the addition of a toroidal field nor ohmic or ambipolar diffusion drastically change this picture in the range of parameters we have explored. The ability of these structures to trap dust particles in this configuration is demonstrated. We conclude that Hall-MRI driven organisation is a plausible scenario which could explain some of the structures found in recent observations.

Accretion discs are ubiquitous in the universe and it is a crucial issue to understand how angular momentum and mass are being radially transported in these objects. Here, we study the role played by non-linear spiral patterns within hydrodynamical and non self-gravitating accretion disc assuming that external disturbances such as infall onto the disc may trigger them. To do so, we computed self-similar solutions that describe discs in which a spiral wave propagates. Such solutions present both shocks and critical sonic points that we carefully analyze. For all allowed temperatures and for several spiral shocks, we calculated the wave structure. In particular we inferred the angle of the spiral patern, the stress it exerts on the disc as well as the associated flux of mass and angular momentum as a function of temperature. We quantified the rate of angular momentum transport by means of the dimensionless $\alpha$ parameter. For the thickest disc we considered (corresponding to $h/r$ values of about 1/3), we found values of $\alpha$ as high as $0.1$, and scaling with the temperature $T$ such that $\alpha \propto T^{3/2} \propto (h/r)^3$. The spiral angle scales with the temperature as $\arctan(r/h)$. The existence of these solutions suggests that perturbations occurring at disc outer boundaries, such as for example perturbations due to infall motions, can propagate deep inside the disc and therefore should not be ignored, even when considering small radii.

The local properties of turbulence driven by the magnetorotational instability (MRI) in rotating, shearing flows are studied in the framework of a shearing-box model. Based on numerical simulations, we propose that the MRI-driven turbulence comprises two components: the large-scale shear-aligned strong magnetic field and the small-scale fluctuations resembling magnetohydrodynamic (MHD) turbulence. The energy spectrum of the large-scale component is close to $k^{-2}$, whereas the spectrum of the small-scale component agrees with the spectrum of strong MHD turbulence $k^{-3/2}$. While the spectrum of the fluctuations is universal, the outer-scale characteristics of the turbulence are not; they depend on the parameters of the system, such as the net magnetic flux. However, there is remarkable universality among the allowed turbulent states -- their intensity $v_0$ and their outer scale $\lambda_0$ satisfy the balance condition $v_0/\lambda_0\sim \mathrm d\Omega/\mathrm d\ln r$, where $\mathrm d\Omega/\mathrm d\ln r$ is the local orbital shearing rate of the flow. Finally, we find no sustained dynamo action in the $\mathrm{Pm}=1$ zero net-flux case for Reynolds numbers as high as $45\,000$, casting doubts on the existence of an MRI dynamo in the $\mathrm{Pm}\leq 1$ regime.

We numerically investigate the dynamics of a 2D non-magnetised protoplanetary disc surrounded by an inflow coming from an external envelope. We find that the accretion shock between the disc and the inflow is unstable, leading to the generation of large-amplitude spiral density waves. These spiral waves propagate over long distances, down to radii at least ten times smaller than the accretion shock radius. We measure spiral-driven outward angular momentum transport with 1e-4 < alpha < 1e-2 for an inflow accretion rate Mout>1e-8 Msun/yr. We conclude that the interaction of the disc with its envelope leads to long-lived spiral density waves and radial angular momentum transport with rates that cannot be neglected in young non-magnetised protostellar discs.

We characterize magnetically driven accretion at radii between 1 au and 100 au in protoplanetary discs, using a series of local non-ideal magnetohydrodynamic (MHD) simulations. The simulations assume a Minimum Mass Solar Nebula (MMSN) disc that is threaded by a net vertical magnetic field of specified strength. Confirming previous results, we find that the Hall effect has only a modest impact on accretion at 30 au, and essentially none at 100 au. At 1-10 au the Hall effect introduces a pronounced bi-modality in the accretion process, with vertical magnetic fields aligned to the disc rotation supporting a strong laminar Maxwell stress that is absent if the field is anti-aligned. In the anti-aligned case, we instead find evidence for bursts of turbulent stress at 5-10 au, which we tentatively identify with the non-axisymmetric Hall-shear instability. The presence or absence of these bursts depends upon the details of the adopted chemical model, which suggests that appreciable regions of actual protoplanetary discs might lie close to the borderline between laminar and turbulent behaviour. Given the number of important control parameters that have already been identified in MHD models, quantitative predictions for disc structure in terms of only radius and accretion rate appear to be difficult. Instead, we identify robust qualitative tests of magnetically driven accretion. These include the presence of turbulence in the outer disc, independent of the orientation of the vertical magnetic fields, and a Hall-mediated bi-modality in turbulent properties extending from the region of thermal ionization to 10 au.

Angular momentum transport in accretion discs is often believed to be due to magnetohydrodynamic turbulence mediated by the magnetorotational instability. Despite an abundant literature on the MRI, the parameters governing the saturation amplitude of the turbulence are poorly understood and the existence of an asymptotic behavior in the Ohmic diffusion regime is not clearly established. We investigate the properties of the turbulent state in the small magnetic Prandtl number limit. Since this is extremely computationally expensive, we also study the relevance and range of applicability of the most common subgrid scale models for this problem. Unstratified shearing boxes simulations are performed both in the compressible and incompressible limits, with a resolution up to 800 cells per disc scale height. The latter constitutes the largest resolution ever attained for a simulation of MRI turbulence. In the presence of a mean magnetic field threading the domain, angular momentum transport converges to a finite value in the small Pm limit. When the mean vertical field amplitude is such that {\beta}, the ratio between the thermal and magnetic pressure, equals 1000, we find {\alpha}~0.032 when Pm approaches zero. In the case of a mean toroidal field for which {\beta}=100, we find {\alpha}~0.018 in the same limit. Both implicit LES and Chollet-Lesieur closure model reproduces these results for the {\alpha} parameter and the power spectra. A reduction in computational cost of a factor at least 16 (and up to 256) is achieved when using such methods. MRI turbulence operates efficiently in the small Pm limit provided there is a mean magnetic field. Implicit LES offers a practical and efficient mean of investigation of this regime but should be used with care, particularly in the case of a vertical field. Chollet-Lesieur closure model is perfectly suited for simulations done with a spectral code.

The magnetorotational (MRI) dynamo has long been considered one of the possible drivers of turbulent angular momentum transport in astrophysical accretion disks. However, various numerical results suggest that this dynamo may be difficult to excite in the astrophysically relevant regime of magnetic Prandtl number (Pm) significantly smaller than unity, for reasons currently not well understood. The aim of this article is to present the first results of an ongoing numerical investigation of the role of both linear and nonlinear dissipative effects in this problem. Combining a parametric exploration and an energy analysis of incompressible nonlinear MRI dynamo cycles representative of the transitional dynamics in large aspect ratio shearing boxes, we find that turbulent magnetic diffusion makes the excitation and sustainment of this dynamo at moderate magnetic Reynolds number (Rm) increasingly difficult for decreasing Pm. This results in an increase in the critical Rm of the dynamo for increasing kinematic Reynolds number (Re), in agreement with earlier numerical results. Given its very generic nature, we argue that turbulent magnetic diffusion could be an important determinant of MRI dynamo excitation in disks, and may also limit the efficiency of angular momentum transport by MRI turbulence in low Pm regimes.

Protoplanetary discs are poorly ionised due to their low temperatures and high column densities, and are therefore subject to three "non-ideal" magnetohydrodynamic effects: Ohmic dissipation, ambipolar diffusion, and the Hall effect. The existence of magnetically driven turbulence in these discs has been a central question since the discovery of the magnetorotational instability. Early models considered Ohmic diffusion only and led to a scenario of layered accretion, in which a magnetically "dead" zone in the disc midplane is embedded within magnetically "active" surface layers at distances ~1-10 au from the central protostellar object. Recent work has suggested that a combination of Ohmic dissipation and ambipolar diffusion can render both the midplane and surface layers of the disc inactive and that torques due to magnetically driven outflows are required to explain the observed accretion rates. We reassess this picture by performing three-dimensional numerical simulations that include, for the first time, all three non-ideal MHD effects. We find that the Hall effect can generically "revive" dead zones by producing a dominant azimuthal magnetic field and a large-scale Maxwell stress throughout the midplane, provided the angular velocity and magnetic field satisfy Omega.B > 0. The attendant large magnetic pressure modifies the vertical density profile and substantially increases the disc scale height beyond its hydrostatic value. Outflows are produced, but are not necessary to explain accretion rates <10^{-7} Msun/yr. The flow in the disc midplane is essentially laminar, suggesting that dust sedimentation may be efficient. These results demonstrate that, if the MRI is relevant for driving mass accretion in protoplanetary discs, one must include the Hall effect to obtain even qualitatively correct results.

Planets appear to form in environments shaped by the gas flowing through protostellar disks to the central young stars. The flows in turn are governed by orbital angular momentum transfer. In this chapter we summarize current understanding of the transfer processes best able to account for the flows, including magneto-rotational turbulence, magnetically-launched winds, self-gravitational instability and vortices driven by hydrodynamical instabilities. For each in turn we outline the major achievements of the past few years and the outstanding questions. We underscore the requirements for operation, especially ionization for the magnetic processes and heating and cooling for the others. We describe the distribution and strength of the resulting flows and compare with the long-used phenomenological $\alpha$-picture, highlighting issues where the fuller physical picture yields substantially different answers. We also discuss the links between magnetized turbulence and magnetically-launched outflows, and between magnetized turbulence and hydrodynamical vortices. We end with a summary of the status of efforts to detect specific signatures of the flows.

The magnetorotational instability (MRI) is the most promising mechanism by which angular momentum is efficiently transported outwards in astrophysical discs. However, its application to protoplanetary discs remains problematic. These discs are so poorly ionised that they may not support magnetorotational turbulence in regions referred to as `dead zones'. It has recently been suggested that the Hall effect, a non-ideal magnetohydrodynamic (MHD) effect, could revive these dead zones by enhancing the magnetically active column density by an order of magnitude or more. We investigate this idea by performing local, three-dimensional, resistive Hall-MHD simulations of the MRI in situations where the Hall effect dominates over Ohmic dissipation. As expected from linear stability analysis, we find an exponentially growing instability in regimes otherwise linearly stable in resistive MHD. However, instead of vigorous and sustained magnetorotational turbulence, we find that the MRI saturates by producing large-scale, long-lived, axisymmetric structures in the magnetic and velocity fields. We refer to these structures as zonal fields and zonal flows, respectively. Their emergence causes a steep reduction in turbulent transport by at least two orders of magnitude from extrapolations based upon resistive MHD, a result that calls into question contemporary models of layered accretion. We construct a rigorous mean-field theory to explain this new behaviour and to predict when it should occur. Implications for protoplanetary disc structure and evolution, as well as for theories of planet formation, are briefly discussed.

Magnetorotational dynamo action in Keplerian shear flow is a three-dimensional, nonlinear magnetohydrodynamic process whose study is relevant to the understanding of accretion and magnetic field generation in astrophysics. Transition to this form of dynamo is subcritical and shares many characteristics of transition to turbulence in non-rotating hydrodynamic shear flows. This suggests that these different fluid systems become active through similar generic bifurcation mechanisms, which in both cases have eluded detailed understanding so far. In this paper, we investigate numerically the bifurcation mechanisms at work in the incompressible Keplerian magnetorotational dynamo problem in the shearing box framework. Using numerical techniques imported from dynamical systems research, we show that the onset of chaotic dynamo action at magnetic Prandtl numbers larger than unity is primarily associated with global homoclinic and heteroclinic bifurcations of nonlinear magnetorotational dynamo cycles. These global bifurcations are supplemented by local bifurcations of cycles marking the beginning of period-doubling cascades. This suggests that nonlinear magnetorotational dynamo cycles provide the pathway to turbulent injection of both kinetic and magnetic energy in incompressible magnetohydrodynamic Keplerian shear flow in the absence of an externally imposed magnetic field. Studying the nonlinear physics and bifurcations of these cycles in different regimes and configurations may subsequently help to better understand the conditions of excitation of magnetohydrodynamic turbulence and instability-driven dynamos in various astrophysical systems and laboratory experiments. The detailed characterization of global bifurcations provided for this three-dimensional subcritical fluid dynamics problem may also prove useful for the problem of transition to turbulence in hydrodynamic shear flows.

The aim of this paper is to investigate the properties of accretion disks threaded by a weak vertical magnetic field, with a particular focus on the interplay between MHD turbulence driven by the magnetorotational instability (MRI) and outflows that might be launched from the disk. For that purpose, we use a set of numerical simulations performed with the MHD code RAMSES in the framework of the shearing box model. We concentrate on the case of a rather weak vertical magnetic field such that the initial ratio beta0 of the thermal and magnetic pressures in the disk midplane equals 10^4. As reported recently, we find that MHD turbulence drives an efficient outflow out of the computational box. We demonstrate a strong sensitivity of that result to the box size: enlargements in the radial and vertical directions lead to a reduction of up to an order of magnitude in the mass-loss rate. Such a dependence prevents any realistic estimates of disk mass-loss rates being derived using shearing-box simulations. We find however that the flow morphology is robust and independent of the numerical details of the simulations. Its properties display some features and approximate invariants that are reminiscent of the Blandford & Payne launching mechanism, but differences exist. For the magnetic field strength considered in this paper, we also find that angular momentum transport is most likely dominated by MHD turbulence, the saturation of which scales with the magnetic Prandtl number, the ratio of viscosity and resistivity, in a way that is in good agreement with expectations based on unstratified simulations. This paper thus demonstrates for the first time that accretion disks can simultaneously exhibit MRI-driven MHD turbulence along with magneto-centrifugally accelerated outflows.

We present high resolution simulations with RAMSES of supersonic colliding stellar winds. The collision results in a double shock structure which is subject to different instabilities. The Kelvin-Helmholtz instability (KHI) introduces some mixing and variability. For isothermal winds, the Non-linear Thin Shell Instability violently affects the interaction region. Properly modelling these instabilities requires a high enough resolution and an adapted numerical method, especially when one of the winds strongly dominates the other one. At large scale, orbital motion is expected to turn the shocked zone into a spiral but we find that in some configurations the KHI may disrupt the spiral. A colliding wind structure is also expected in gamma-ray binaries composed of a massive star and a young pulsar which emits a highly relativistic wind. Numerical simulations are necessary to understand the geometry of such systems and should take into account the relativistic nature of the pulsar wind. We implemented a second order Godunov method to solve the equations of relativistic hydrodynamics with RAMSES, including the possibility for Adaptive Mesh Refinement. After a brief overview of our numerical implementation, we will present preliminary simulation that contrast the structure of the flow in gamma-ray and stellar binaries.

Gamma-ray binaries are colliding wind binaries (CWB) composed of a massive star a non-accreting pulsar with a highly relativistic wind. Particle acceleration at the shocks results in emission going from extended radio emission to the gamma-ray band. The interaction region is expected to show common features with stellar CWB. Performing numerical simulations with the hydrodynamical code RAMSES, we focus on their structure and stability and find that the Kelvin-Helmholtz instability (KHI) can lead to important mixing between the winds and destroy the large scale spiral structure. To investigate the impact of the relativistic nature of the pulsar wind, we extend RAMSES to relativistic hydrodynamics (RHD). Preliminary simulations of the interaction between a pulsar wind and a stellar wind show important similarities with stellar colliding winds with small relativistic corrections.

Magnetorotational turbulence and magnetically driven disc winds are often considered as separate processes. However, realistic astrophysical discs are expected to be subject to both effects, although possibly at different times and locations. We investigate here the potential link between these two phenomena using a mixed numerical and analytical approach. We show in particular that large-scale MRI modes which dominate strongly magnetised discs (plasma beta~10) naturally produce magnetically driven outflows in the nonlinear regime. We show that these outflows share many similarities with local and global disc wind solutions found in the literature. We also investigate the 3D stability of these outflows and show that they are unstable on dynamical timescales. The implications of these results for the transition between a jet-emitting disc and a standard "viscous" disc are discussed.

The collision of winds from massive stars in binaries results in the formation of a double-shock structure with observed signatures from radio to X-rays. We study the structure and stability of the colliding wind region as it turns into a spiral due to orbital motion. We focus on adiabatic winds, where mixing between the two winds is expected to be restricted to the Kelvin-Helmholtz instability (KHI). Mixing of the Wolf-Rayet wind with hydrogen-rich material is important for dust formation in pinwheel nebulae such as WR 104, where the spiral structure has been resolved in infrared. We use the hydrodynamical code RAMSES with an adaptive grid. A wide range of binary systems with different wind velocities and mass loss rates are studied with 2D simulations. A specific 3D simulation is performed to model WR 104. Orbital motion leads to the formation of two distinct spiral arms where the KHI develops differently. We find that the spiral structure is destroyed when there is a large velocity gradient between the winds, unless the collimated wind is much faster. We argue that the KHI plays a major role in maintaining or not the structure. We discuss the consequences for various colliding wind binaries. When stable, there is no straightforward relationship between the spatial step of the spiral, the wind velocities, and the orbital period. Our 3D simulation of WR 104 indicates that the colder, well-mixed trailing arm has more favourable conditions for dust formation than the leading arm. The single-arm infrared spiral follows more closely the mixing map than the density map, suggesting the dust-to-gas ratio may vary between the leading and trailing density spirals. However, the density is much lower than what dust formation models require. Including radiative cooling would lead to higher densities, and also to thin shell instabilities whose impact on the large structure remains unknown.

We study the interaction between tides and convection in astrophysical bodies by analysing the effect of a homogeneous oscillatory shear on a fluid flow. This model can be taken to represent the interaction between a large-scale periodic tidal deformation and a smaller-scale convective motion. We first consider analytically the limit in which the shear is of low amplitude and the oscillation period is short compared to the timescales of the unperturbed flow. In this limit there is a viscoelastic response and we obtain expressions for the effective elastic modulus and viscosity coefficient. The effective viscosity is inversely proportional to the square of the oscillation frequency, with a coefficient that can be positive, negative or zero depending on the properties of the unperturbed flow. We also carry out direct numerical simulations of Boussinesq convection in an oscillatory shearing box and measure the time-dependent Reynolds stress. The results indicate that the effective viscosity of turbulent convection falls rapidly as the oscillation frequency is increased, attaining small negative values in the cases we have examined, although significant uncertainties remain because of the turbulent noise. We discuss the implications of this analysis for astrophysical tides.

The nature of dynamo action in shear flows prone to magnetohydrodynamic instabilities is investigated using the magnetorotational dynamo in Keplerian shear flow as a prototype problem. Using direct numerical simulations and Newton's method, we compute an exact time-periodic magnetorotational dynamo solution to the three-dimensional dissipative incompressible magnetohydrodynamic equations with rotation and shear. We discuss the physical mechanism behind the cycle and show that it results from a combination of linear and nonlinear interactions between a large-scale axisymmetric toroidal magnetic field and non-axisymmetric perturbations amplified by the magnetorotational instability. We demonstrate that this large scale dynamo mechanism is overall intrinsically nonlinear and not reducible to the standard mean-field dynamo formalism. Our results therefore provide clear evidence for a generic nonlinear generation mechanism of time-dependent coherent large-scale magnetic fields in shear flows and call for new theoretical dynamo models. These findings may offer important clues to understand the transitional and statistical properties of subcritical magnetorotational turbulence.

The magnetorotational instability (MRI) is believed to be responsible for most of the angular momentum transport in accretion discs. However, molecular dissipation processes may drastically change the efficiency of MRI turbulence in realistic astrophysical situations. The physical origin of this dependency is still poorly understood as linear and quasi linear theories fail to explain it. In this paper, we look for the link between molecular dissipation processes and MRI transport of angular momentum in non stratified shearing box simulations including a mean vertical field. We show that magnetic helicity is unimportant in the model we consider. We perform a spectral analysis on the simulations tracking energy exchanges in spectral space when turbulence is fully developed. We find that the energy exchanges are essentially direct (from large to small scale) whereas some non linear interactions appear to be non local in spectral space. We speculate that these non local interactions are responsible for the correlation between turbulent transport and molecular dissipation. We argue that this correlation should then disappear when a significant scale separation is achieved and we discuss several methods by which one can test this hypothesis.

We consider the radial migration of vortices in two-dimensional isothermal gaseous disks. We find that a vortex core, orbiting at the local gas velocity, induces velocity perturbations that propagate away from the vortex as density waves. The resulting spiral wave pattern is reminiscent of an embedded planet. There are two main causes for asymmetries in these wakes: geometrical effects tend to favor the outer wave, while a radial vortensity gradient leads to an asymmetric vortex core, which favors the wave at the side that has the lowest density. In the case of asymmetric waves, which we always find except for a disk of constant pressure, there is a net exchange of angular momentum between the vortex and the surrounding disk, which leads to orbital migration of the vortex. Numerical hydrodynamical simulations show that this migration can be very rapid, on a time scale of a few thousand orbits, for vortices with a size comparable to the scale height of the disk. We discuss the possible effects of vortex migration on planet formation scenarios.

A subcritical transition to turbulence in magnetized Keplerian shear flows is investigated using a statistical approach. Three-dimensional numerical simulations of the shearing box equations with zero net magnetic flux are employed to determine the transition from decaying to sustained turbulence as a function of the magnetic Reynolds number Rm. The results reveal no clear transition to sustained turbulence as the average lifetime of the transients grows as an exponential function of Rm, in accordance with a type-II supertransient law.

The Sun is a magnetic star whose cyclic activity is thought to be linked to internal dynamo mechanisms. A combination of numerical modelling with various levels of complexity is an efficient and accurate tool to investigate such intricate dynamical processes. We investigate the role of the magnetic buoyancy process in 2D Babcock-Leighton dynamo models, by modelling more accurately the surface source term for poloidal field. Methods. To do so, we reintroduce in mean-field models the results of full 3D MHD calculations of the non-linear evolution of a rising flux tube in a convective shell. More specifically, the Babcock-Leighton source term is modified to take into account the delay introduced by the rise time of the toroidal structures from the base of the convection zone to the solar surface. We find that the time delays introduced in the equations produce large temporal modulation of the cycle amplitude even when strong and thus rapidly rising flux tubes are considered. Aperiodic modulations of the solar cycle appear after a sequence of period doubling bifurcations typical of non-linear systems. The strong effects introduced even by small delays is found to be due to the dependence of the delays on the magnetic field strength at the base of the convection zone, the modulation being much less when time delays remain constant. We do not find any significant influence on the cycle period except when the delays are made artificially strong. A possible new origin of the solar cycle variability is here revealed. This modulated activity and the resulting butterfly diagram are then more compatible with observations than what the standard Babcock-Leighton model produces.

In the recent years, MRI-driven turbulent transport has been found to depend in a significant way on fluid viscosity $\nu$ and resistivity $\eta$ through the magnetic Prandtl number $Pm=\nu/\eta$. In particular, the transport decreases with decreasing $Pm$; if persistent at very large Reynolds numbers, this trend may lead to question the role of MRI-turbulence in YSO disks, whose Prandtl number is usually very small. In this context, the principle objective of the present investigation is to characterize in a refined way the role of dissipation. Another objective is to characterize the effect of linear (channel modes) and quasi-linear (parasitic modes) physics in the behavior of the transport. These objectives are addressed with the help of a number of incompressible numerical simulations. The horizontal extent of the box size has been increased in order to capture all relevant (fastest growing) linear and secondary parasitic unstable modes. The major results are the following: i- The increased accuracy in the computation of transport averages shows that the dependence of transport on physical dissipation exhibits two different regimes: for $Pm \lesssim 1$, the transport has a power-law dependence on the magnetic Reynolds number rather than on the Prandtl number; for $Pm > 1$, the data are consistent with a primary dependence on $Pm$ for large enough ($\sim 10^3$) Reynolds numbers. ii- The transport-dissipation correlation is not clearly or simply related to variations of the linear modes growth rates. iii- The existence of the transport-dissipation correlation depends neither on the number of linear modes captured in the simulations, nor on the effect of the parasitic modes on the saturation of the linear modes growth. iv- The transport is usually not dominated by axisymmetric (channel) modes

The mechanism of angular momentum transport in accretion discs has long been debated. Although the magnetorotational instability appears to be a promising process, poorly ionized regions of accretion discs may not undergo this instability. In this letter, we revisit the possibility of transporting angular momentum by turbulent thermal convection. Using high-resolution spectral methods, we show that strongly turbulent convection can drive outward angular momentum transport at a rate that is, under certain conditions, compatible with observations of discs. We find however that the angular momentum transport is always much weaker than the vertical heat transport. These results indicate that convection might be another way to explain global disc evolution, provided that a sufficiently unstable vertical temperature profile can be maintained.

(abridged) Aims: We present new results exhibiting a subcritical baroclinic instability (SBI) in local shearing box models. We describe the 2D and 3D behaviour of this instability using numerical simulations and we present a simple analytical model describing the underlying physical process. Results: A subcritical baroclinic instability is observed in flows stable for the Solberg-Hoiland criterion using local simulations. This instability is found to be a nonlinear (or subcritical) instability, which cannot be described by ordinary linear approaches. It requires a radial entropy gradient weakly unstable for the Schwartzchild criterion and a strong thermal diffusivity (or equivalently a short cooling time). In compressible simulations, the instability produces density waves which transport angular momentum outward with typically alpha<3e-3, the exact value depending on the background temperature profile. Finally, the instability survives in 3D, vortex cores becoming turbulent due to parametric instabilities. Conclusions: The subcritical baroclinic instability is a robust phenomenon, which can be captured using local simulations. The instability survives in 3D thanks to a balance between the 2D SBI and 3D parametric instabilities. Finally, this instability can lead to a weak outward transport of angular momentum, due to the generation of density waves by the vortices.

The formation mechanism of planetesimals in protoplanetary discs is hotly debated. Currently, the favoured model involves the accumulation of meter-sized objects within a turbulent disc, followed by a phase of gravitational instability. At best one can simulate a few million particles numerically as opposed to the several trillion meter-sized particles expected in a real protoplanetary disc. Therefore, single particles are often used as super-particles to represent a distribution of many smaller particles. It is assumed that small scale phenomena do not play a role and particle collisions are not modeled. The super-particle approximation can only be valid in a collisionless or strongly collisional system, however, in many recent numerical simulations this is not the case. In this work we present new results from numerical simulations of planetesimal formation via gravitational instability. A scaled system is studied that does not require the use of super-particles. We find that the scaled particles can be used to model the initial phases of clumping if the properties of the scaled particles are chosen such that all important timescales in the system are equivalent to what is expected in a real protoplanetary disc. Constraints are given for the number of particles needed in order to achieve numerical convergence. We compare this new method to the standard super-particle approach. We find that the super-particle approach produces unreliable results that depend on artifacts such as the gravitational softening in both the requirement for gravitational collapse and the resulting clump statistics. Our results show that short range interactions (collisions) have to be modelled properly.

(abriged) MRI turbulence is a leading mechanism for the generation of an efficient turbulent transport of angular momentum in an accretion disk through a turbulent viscosity effect. It is believed that the same process could also transport large-scale magnetic fields in disks, reshaping the magnetic structures in these objects. This process, known as turbulent resistivity, has been suggested and used in several accretion-ejection models and simulations to produce jets. Still, the efficiency of MRI-driven turbulence to transport large-scale magnetic fields is largely unknown. We investigate this problem both analytically and numerically. We introduce a linear calculation of the MRI in the presence of a spatially inhomogeneous mean magnetic field. We show that, in this configuration, MRI modes lead to an efficient magnetic field transport, on the order of the angular momentum transport. We next use fully non linear simulations of MRI turbulence to compute the turbulent resistivity in several magnetic configurations. We find that the turbulent resistivity is on the order of the turbulent viscosity in all our simulations, although somewhat lower. The turbulent resistivity tensor is found to be highly anisotropic with a diffusion coefficient 3 times greater in the radial direction than in the vertical direction. These results support the possibility of driving jets from turbulent disks; the resulting jets may not be steady.

(Abriged) The existence of large-scale and long-lived 2D vortices in accretion discs has been debated for more than a decade. They appear spontaneously in several 2D disc simulations and they are known to accelerate planetesimal formation through a dust trapping process. However, the issue of the stability of these structures to the imposition of 3D disturbances is still not fully understood, and it casts doubts on their long term survival. Aim: We present new results on the 3D stability of elliptical vortices embedded in accretion discs, based on a linear analysis and several non-linear simulations. Methods: We derive the linearised equations governing the 3D perturbations in the core of an elliptical vortex, and we show that they can be reduced to a Floquet problem. We solve this problem numerically in the astrophysical regime and we present several analytical limits for which the mechanism responsible for the instability can be explained. Finally, we compare the results of the linear analysis to some high resolution simulations. Results: We show that most anticyclonic vortices are unstable due to a resonance between the turnover time and the local epicyclic oscillation period. In addition, we demonstrate that a strong vertical stratification does not create any additional stable domain of aspect ratio, but it significantly reduces growth rates for relatively weak (and therefore elongated) vortices. Conclusions: Elliptical vortices are always unstable, whatever the horizontal or vertical aspect-ratio is. The instability can however be weak and is often found at small scales, making it difficult to detect in low-order finite-difference simulations.

(Abriged) The magnetorotational instability (MRI) is believed to be an efficient way to transport angular momentum in accretion discs. It has also been suggested as a way to amplify magnetic fields in discs, the instability acting as a nonlinear dynamo. Recent numerical work has shown that a large-scale magnetic field, which is predominantly azimuthal, can be sustained by motions driven by the MRI of this same field. Following this idea, we present an analytical calculation of the MRI in the presence of an azimuthal field with a non-trivial vertical structure. We find that the mean radial EMF associated to MRI modes tends to reduce the magnetic energy, acting like a turbulent resistivity by mixing the non-uniform azimuthal field. Meanwhile, the azimuthal EMF generates a radial field that, in combination with the Keplerian shear, tends to amplify the azimuthal field and can therefore assist in the dynamo process. This effect, however, is reversed for sufficiently strong azimuthal fields, naturally leading to a saturation of the dynamo and possibly to a cyclic behaviour of the magnetic field, as found in previous numerical works.

(abridged) MHD turbulence is known to exist in shearing boxes with either zero or nonzero net magnetic flux. However, the way turbulence survives in the zero-net-flux case is not explained by linear theory and appears as a purely numerical result. Aims: We look for a nonlinear mechanism able to explain the persistence of MHD turbulence in shearing boxes with zero net magnetic flux, and potentially leading to large-scale dynamo action. Method: Spectral nonlinear simulations of the magnetorotational instability are shown to exhibit a large-scale axisymmetric magnetic field, maintained for a few orbits. The generation process of this field is investigated using the results of the simulations and an inhomogeneous linear approach. Results: The mechanism by which turbulence is sustained in zero-net-flux shearing boxes is shown to be related to the existence of a large-scale azimuthal field, surviving for several orbits. In particular, it is shown that MHD turbulence in shearing boxes can be seen as a dynamo process coupled to a magnetorotational-type instability.

We study the influence of the choice of transport coefficients (viscosity and resistivity) on MHD turbulence driven by the magnetorotational instability (MRI) in accretion disks. We follow the methodology described in paper I: we adopt an unstratified shearing box model and focus on the case where the net vertical magnetic flux threading the box vanishes. For the most part we use the finite difference code ZEUS, including explicit transport coefficients in the calculations. However, we also compare our results with those obtained using other algorithms (NIRVANA, the PENCIL code and a spectral code) to demonstrate both the convergence of our results and their independence of the numerical scheme. We find that small scale dissipation affects the saturated state of MHD turbulence. In agreement with recent similar numerical simulations done in the presence of a net vertical magnetic flux, we find that turbulent activity (measured by the rate of angular momentum transport) is an increasing function of the magnetic Prandtl number Pm for all values of the Reynolds number Re that we investigated. We also found that turbulence disappears when the Prandtl number falls below a critical value Pm_c that is apparently a decreasing function of Re. For the limited region of parameter space that can be probed with current computational resources, we always obtained Pm_c>1. We conclude that the magnitudes of the transport coefficients are important in determining the properties of MHD turbulence in numerical simulations in the shearing box with zero net flux, at least for Reynolds numbers and magnetic Prandtl numbers that are such that transport is not dominated by numerical effects and thus can be probed using current computational resources.

The magneto-rotational instability is presently the most promising source of turbulent transport in accretion disks. However, some important issues still need to be addressed to quantify the role of MRI in disks; in particular no systematic investigation of the role of the physical dimensionless parameters of the problem on the dimensionless transport has been undertaken yet. First, we complete existing investigations on the field strength dependence by showing that the transport in high magnetic pressure disks close to marginal stability is highly time-dependent and surprisingly efficient. Second, we bring to light a significant dependence of the global transport on the magnetic Prandtl number, with $\alpha\propto Pm^\delta$ for the explored range: $0.12<Pm<8$ and $200<Re<6400$ ($\delta$ being in the range 0.25 to 0.5). We show that the dimensionless transport is not correlated to the dimensionless linear growth rate, contrarily to a largely held expectation. More generally, these results stress the need to control dissipation processes in astrophysical simulations.

Hydrodynamic unstratified keplerian flows are known to be linearly stable at all Reynolds numbers, but may nevertheless become turbulent through nonlinear mechanisms. However, in the last ten years, conflicting points of view have appeared on this issue. We have revisited the problem through numerical simulations in the shearing sheet limit. It turns out that the effect of the Coriolis force in stabilizing the flow depends on whether the flow is cyclonic (cooperating shear and rotation vorticities) or anticyclonic (competing shear and rotation vorticities); keplerian flows are anticyclonic. We have obtained the following results: i/ The Coriolis force does not quench turbulence in subcritical flows; ii/ The resolution demand, when moving away from the marginal stability boundary, is much more severe for anticyclonic flows than for cyclonic ones. Presently available computer resources do not allow numerical codes to reach the keplerian regime. iii/ The efficiency of turbulent transport is directly correlated to the Reynolds number of transition to turbulence $Rg$, in such a way that the Shakura-Sunyaev parameter $\alpha\sim 1/Rg$. iv/ Even the most optimistic extrapolations of our numerical data show that subcritical turbulent transport would be too inefficient in keplerian flows by several orders of magnitude for astrophysical purposes. v/ Our results suggest that the data obtained for keplerian-like flows in a Taylor-Couette settings are largely affected by secondary flows, such as Ekman circulation.