Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSpiral-driven accretion in protoplanetary discs - III tri-dimensional simulations
75
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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.
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.
We study the gravitational instability and fragmentation of primordial protostellar discs by using high-resolution cosmological hydrodynamics simulations. We follow the formation and evolution of spiral arms in protostellar discs, examine the dynamical stability, and identify a physical mechanism of secondary protostar formation. We use linear perturbation theory based on the spiral-arm instability (SAI) analysis in our previous studies. We improve the analysis by incorporating the effects of finite thickness and shearing motion of arms, and derive the physical conditions for SAI in protostellar discs. Our analysis predicts accurately the stability and the onset of arm fragmentation that is determined by the balance between self-gravity and gas pressure plus the Coriolis force. Formation of secondary and multiple protostars in the discs is explained by the SAI, which is driven by self-gravity and thus can operate without rapid gas cooling. We can also predict the typical mass of the fragments, which is found to be in good agreement with the actual masses of secondary protostars formed in the simulation.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
