No Arabic abstract
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.
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.
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.
Spiral arms have been observed in nearly a dozen protoplanetary discs in near-infrared scattered light and recently also in the sub-millimetre continuum. While one of the most compelling explanations is that they are driven by planetary or stellar companions, in all but one cases such companions have not yet been detected and there is even ambiguity on whether the planet should be located inside or outside the spirals. Here we use 3D hydrodynamic simulations to study the morphology of spiral density waves launched by embedded planets taking into account the vertical temperature gradient, a natural consequence of stellar irradiation. Our simulations show that the pitch angle of the spirals in thermally stratified discs is the lowest in the disc mid-plane and increases towards the disc surface. We combine the hydrodynamic simulations with 3D radiative transfer calculations to predict that the pitch-angle of planetary spirals observed in the near-infrared is higher than in the sub-millimetre. We also find that in both cases the spirals converge towards the planet. This provides a new powerful observational method to determine if the perturbing planet is inside our outside the spirals, as well as map the thermal stratification of the disc.
The formation of protoplanetary discs during the collapse of molecular dense cores is significantly influenced by angular momentum transport, notably by the magnetic torque. In turn, the evolution of the magnetic field is determined by dynamical processes and non-ideal MHD effects such as ambipolar diffusion. Considering simple relations between various timescales characteristic of the magnetized collapse, we derive an expression for the early disc radius, $ r simeq 18 , {rm AU} , left({eta_{rm AD} / 0.1 , {rm s}} right)^{2/9} left({B_z / 0.1, {rm G}} right) ^{-4/9} left({M / 0.1 msol} right) ^{1/3},$ where $M$ is the total disc plus protostar mass, $eta_mathrm{AD}$ is the ambipolar diffusion coefficient and $B_z$ is the magnetic field in the inner part of the core. This is about significantly smaller than the discs that would form if angular momentum was conserved. The analytical predictions are confronted against a large sample of 3D, non-ideal MHD collapse calculations covering variations of a factor 100 in core mass, a factor 10 in the level of turbulence, a factor 5 in rotation, and magnetic mass-to-flux over critical mass-to-flux ratios 2 and 5. The disc radius estimates are found to agree with the numerical simulations within less than a factor 2. A striking prediction of our analysis is the weak dependence of circumstellar disc radii upon the various relevant quantities, suggesting weak variations among class-0 disc sizes. In some cases, we note the onset of large spiral arms beyond this radius.
We have performed three-dimensional numerical simulations of accretion discs in a close binary system using the Smoothed Particle Hydrodynamics method. Our result show that, contrary to previous claims, 3D discs do exist even when the specific heat ratio of the gas is as large as gamma=1.2. Although the disc is clearly more spread in the z-direction in this case than it is for the quasi-isothermal one, the disc height is compatible with the hydrostatic balance equation. Our numerical simulations with gamma=1.2 also demonstrate that spiral shocks exist in 3D discs. These results therefore confirm previous 2D simulations.