Chondrules are crystallised droplets of silicates formed by rapid heating to high temperatures (> 1800 K) of solid precursors followed by hours or days of cooling. Dating of chondrules is consistent with the formation timescale of Jupiter in the core-accretion model (1-4 Myrs). Here we investigate if the shocks generated by a massive planet could generate flash heating episodes necessary to form chondrules using high resolution 2D simulations with the multi-fluid code ROSSBI. We use different radiative cooling prescriptions, i.e. different cooling rates and models, and vary planet mass, orbit and disk models. For long disk cooling rates (> 1000 orbits) and a massive protoplanet (> 0.75 jupiter masses), we obtain hot enough temperatures for chondrule formation, while using more realistic thermodynamics is not successful in the Minimum Mass Solar Nebula (MMSN) model. However, sufficient flash heating is achieved by a Jupiter mass planet in a 5 times more massive disk, which is a conceivable scenario for the young solar nebula and exoplanetary systems. Interestingly, a gap-forming massive planet triggers vortices, which can trap dust, i.e. chondrule precursors, and generates a high pressure environment that is consistent with cosmochemical evidence from chondrules. A massive gas giant can thus, in principle, both stimulate the concentration of chondrule precursors in the vicinity of the shocking regions, and flash-heat them via the shocks.

We present a novel study of dust-vortex evolution in global two-fluid disk simulations to find out if evolution toward high dust-to-gas ratios can occur in a regime of well-coupled grains with low Stokes numbers ($St=10^{-3}-{4\times 10^{-2}}$). We design a new implicit scheme in the code RoSSBi, to overcome the short timesteps occurring for small grain sizes. We discover that the linear capture phase occurs self-similarly for all grain sizes, with an intrinsic timescale (characterizing the vortex lifetime) scaling as $1/St$. After vortex dissipation, the formation of a {{global active dust ring}} is a generic outcome confirming our previous results obtained for larger grains. We propose a scenario in which, irrespective of grain size, multiple pathways can lead to local dust-to-gas ratios of order unity and above on relatively short timescales, $< 10^5$ yr, in the presence of a vortex, even with $St=10^{-3}$. When $St>10^{-2}$, the vortex is quickly dissipated by two-fluid instabilities, and large dust density enhancements form in the global dust ring. When $St<10^{-2}$, the vortex is resistant to destabilization. As a result, dust concentrations occur locally due to turbulence developing inside the vortex. Whatever the Stokes number, dust-to-gas ratios in the range $1-10$, a necessary condition to trigger a subsequent streaming instability, or even a direct gravitational instability of the dust clumps, appears to be an inevitable outcome. Although quantitative connections with other instabilities still need to be made, we argue that our results support a new scenario of vortex-driven planetesimal formation.

We present the first 2D hydrodynamical finite volume simulations in which dust is fully coupled with the gas, including its back-reaction onto it, and at the same time the dust size is evolving according to coagulation and fragmentation based on a sub-grid model. The aim of this analysis is to present the differences occurring when dust evolution is included relative to simulations with fixed dust size, with and without an embedded Jupiter-mass planet that triggers gap formation. We use the two-fluid polar Godunov-type code RoSSBi developed by Surville et al. 2016 combined with a new local sub-grid method for dust evolution based on the model by Birnstiel et al. 2012. We find striking differences between simulations with variable and fixed dust sizes. The timescales for dust depletion differ significantly and yield a completely different evolution of the dust surface density. In general sharp features such as pile-ups of dust in the inner disk and near gap edges, when a massive planet is present, become much weaker. This has important implications on the interpretation of observed substructure in disks, suggesting that the presence of a massive planet does not necessarily cause sharp gaps and rings in the dust component. Also, particles with different dust sizes show a different distribution, pointing to the importance of multi-wavelength synthetic observations in order to compare with observations by ALMA and other instruments. We also find that simulations adopting fixed intermediate particle sizes, in the range $10^{-2} - 10^{-1}$ cm, best approximate the surface density evolution seen in simulations with dust evolution.

We study dust capture by vortices and its long-term consequences in global two-fluid inviscid disk simulations using a new polar grid code RoSSBi. We perform the longest integrations so far, several hundred disk orbits, at the highest resolution attainable in global simulations of disks with dust, namely 2048x4096 grid points. This allows to study the dust evolution well beyond vortex dissipation. We vary a wide range of parameters, most notably the dust-to-gas ratio in the initial setup varies in the range $10^{-3}$ to $0.1$. Irrespective of the initial dust-to-gas ratio we find rapid concentration of the dust inside vortices, reaching dust-to-gas ratios of order unity inside the vortex. We present an analytical model that describes very well the dust capture process inside vortices, finding consistent results for all dust-to-gas ratios. A vortex streaming instability develops which causes invariably vortex destruction. After vortex dissipation large-scale dust-rings encompassing a disk annulus form in most cases, which sustain very high dust concentration, approaching ratios of order unity. The rings are long lived lasting as long as the duration of the simulations. They also develop a streaming instability, which manifests itself in eddies at various scales within which the dust forms compact high density clumps. Such clumps would be unstable to gravitational collapse in absence of strong dissipation by viscous forces. When vortices are particularly long lived, rings do not form but dust clumps inside vortices become then long lived features and would likely undergo collapse by gravitational instability. Rings encompass almost an Earth mass of solid material, while even larger masses of dust do accumulate inside vortices in the earlier stage. We argue that rapid planetesimal formation would occur in the dust clumps inside the vortices as well as in the post-vortex ring.

We study the evolution of planet-induced vortices in radially stratified disks, with initial conditions allowing for radial buoyancy. For this purpose we run global two dimensional hydrodynamical simulations, using the PLUTO code. Planet-induced vortices are a product of the Rossby wave instability (RWI) triggered in the edges of a planetary gap. In this work we assess the influence of radial buoyancy for the development of the vortices. We found that radial buoyancy leads to smoother planetary gaps, which generates weaker vortices. This effect is less pronounced for locally isothermal and quasi-isothermal (very small cooling rate) disks. We observed the formation of two generations of vortices. The first generation of vortices is formed in the outer wall of the planetary gap. The merged primary vortex induces accretion, depleting the mass on its orbit. This process creates a surface density enhancement beyond the primary vortex position. The second generation of vortices arise in this surface density enhancement, indicating that the bump in this region is sufficient to trigger the RWI. The merged secondary vortex is a promising explanation for the location of the vortex in the Oph IRS48 system. Finally, we observed a nonmonotonic behavior for the vortex lifetimes as a function of the thermal relaxation timescale, agreeing with previous studies. The birth times of the secondary vortices also display a nonmonotonic behavior, which is correlated with the growth time of the primary vortex and its induced accretion.

This contribution describes the evolution of the protoplanetary disk using 2D numerical simulations. The 2D Euler equations are solved with the finite volume method. The numerical simulations are used to study the persistence and migration of anticyclonic vortices. Two cases are presented : (1) vortices produced by the Rossby Waves Instability, (2) a non-linear vortex model initially implemented into the disk. The migration of the vortices is due to spiral density waves excited by the vortex in the gas of the disk.

We exploit our formula for the gravitational potential of finite size, power-law disks to derive a general expression linking the mass of the black hole in active galactic nuclei (AGN), the mass of the surrounding disk, its surface density profile (through the power index s), and the differential rotation law. We find that the global rotation curve v(R) of the disk in centrifugal balance does not obey a power law of the cylindrical radius R (except in the confusing case s = -2 that mimics a Keplerian motion), and discuss the local velocity index. This formula can help to understand how, from position-velocity diagrams, mass is shared between the disk and the black hole. To this purpose, we have checked the idea by generating a sample of synthetic data with different levels of Gaussian noise, added in radius. It turns out that, when observations are spread over a large radial domain and exhibit low dispersion (standard deviation less than 10% typically), the disk properties (mass and s-parameter) and black hole mass can be deduced from a non linear fit of kinematic data plotted on a (R, Rv 2)-diagram. For a deviation higher than 10%, masses are estimated fairly well from a linear regression (corresponding to the zeroth-order treatment of the formula), but the power index s is no longer accessible. We have applied the model to 7 AGN disks whose rotation has already been probed through water maser emission. For NGC3393 and UGC3789, the masses seem well constrained through the linear approach. For IC1481, the power-law exponent s can even be deduced. Because the model is scale-free, it applies to any kind of star/disk system. Extension to disks around young stars showing deviation from Keplerian motion is thus straightforward.