-
A Fast second-order solver for stiff multifluid dust and gas hydrodynamics
Authors:
Leonardo Krapp,
Juan Garrido-Deutelmoser,
Pablo Benítez-Llambay,
Kaitlin M. Kratter
Abstract:
We present MDIRK: a Multifluid second-order Diagonally-Implicit Runge-Kutta method to study momentum transfer between gas and an arbitrary number ($N$) of dust species. The method integrates the equations of hydrodynamics with an Implicit Explicit (IMEX) scheme and solves the stiff source term in the momentum equation with a diagonally-implicit asymptotically stable Runge-Kutta method (DIRK). In p…
▽ More
We present MDIRK: a Multifluid second-order Diagonally-Implicit Runge-Kutta method to study momentum transfer between gas and an arbitrary number ($N$) of dust species. The method integrates the equations of hydrodynamics with an Implicit Explicit (IMEX) scheme and solves the stiff source term in the momentum equation with a diagonally-implicit asymptotically stable Runge-Kutta method (DIRK). In particular, DIRK admits a simple analytical solution that can be evaluated with $\mathcal{O}(N)$ operations, instead of standard matrix inversion, which is $\mathcal{O}(N)^3$. Therefore the analytical solution significantly reduces the computational cost of the multifluid method, making it suitable for studying the dynamics of systems with particle-size distributions. We demonstrate that the method conserves momentum to machine precision and converges to the correct equilibrium solution with constant external acceleration. To validate our numerical method we present a series of simple hydrodynamic tests, including damping of sound waves, dusty shocks, a multi-fluid dusty Jeans instability, and a steady-state gas-dust drift calculation. The simplicity of MDIRK lays the groundwork to build fast high-order asymptotically stable multifluid methods.
△ Less
Submitted 26 February, 2024; v1 submitted 29 September, 2023;
originally announced October 2023.
-
A gap-sharing planet pair shaping the crescent in HD 163296: a disk sculpted by a resonant chain
Authors:
Juan Garrido-Deutelmoser,
Cristobal Petrovich,
Carolina Charalambous,
Viviana V. Guzmán,
Ke Zhang
Abstract:
ALMA observations of the disk around HD 163296 have resolved a crescent-shape substructure at around 55 au, inside and off-center from a gap in the dust that extends from 38 au to 62 au. In this work we propose that both the crescent and the dust rings are caused by a compact pair (period ratio $\simeq 4:3$) of sub-Saturn-mass planets inside the gap, with the crescent corresponding to dust trapped…
▽ More
ALMA observations of the disk around HD 163296 have resolved a crescent-shape substructure at around 55 au, inside and off-center from a gap in the dust that extends from 38 au to 62 au. In this work we propose that both the crescent and the dust rings are caused by a compact pair (period ratio $\simeq 4:3$) of sub-Saturn-mass planets inside the gap, with the crescent corresponding to dust trapped at the $L_5$ Lagrange point of the outer planet. This interpretation also reproduces well the gap in the gas recently measured from the CO observations, which is shallower than what is expected in a model where the gap is carved by a single planet. Building on previous works arguing for outer planets at $\approx 86$ and $\approx 137$ au, we provide with a global model of the disk that best reproduces the data and show that all four planets may fall into a long resonant chain, with the outer three planets in a 1:2:4 Laplace resonance. We show that this configuration is not only an expected outcome from disk-planet interaction in this system, but it can also help constraining the radial and angular position of the planet candidates using three-body resonances.
△ Less
Submitted 30 January, 2023;
originally announced January 2023.
-
Substructures in protoplanetary disks imprinted by compact planetary systems
Authors:
Juan Garrido-Deutelmoser,
Cristobal Petrovich,
Leonardo Krapp,
Kaitlin M. Kratter,
Ruobing Dong
Abstract:
The substructures observed in protoplanetary disks may be the signposts of embedded planets carving gaps or creating vortices. The inferred masses of these planets often fall in the Jovian regime despite their low abundance compared to lower-mass planets, partly because previous works often assume that a single substructure (a gap or vortex) is caused by a single planet. In this work, we study the…
▽ More
The substructures observed in protoplanetary disks may be the signposts of embedded planets carving gaps or creating vortices. The inferred masses of these planets often fall in the Jovian regime despite their low abundance compared to lower-mass planets, partly because previous works often assume that a single substructure (a gap or vortex) is caused by a single planet. In this work, we study the possible imprints of compact systems composed of Neptune-like planets ($\sim10-30\;M_\oplus$) and show that long-standing vortices are a prevalent outcome when their inter-planetary separation ($Δa$) falls below $\sim8$ times $H_{\rm p}$ -- the average disk's scale height at the planets locations. In simulations where a single planet is unable to produce long-lived vortices, two-planet systems can preserve them for at least $5,000$ orbits in two regimes: i) fully-shared density gaps with elongated vortices around the stable Lagrange points $L_4$ and $L_5$ for the most compact planet pairs ($Δa \lesssim 4.6\; H_{\rm p}$); ii) partially-shared gaps for more widely spaced planets ($Δa \sim 4.6 - 8\;H_{\rm p}$) forming vortices in a density ring between the planets through the Rossby wave instability. The latter case can produce vortices with a wide range of aspect ratios down to $\sim3$ and can occur for planets captured into the 3:2 (2:1) mean-motion resonances for disk's aspects ratios of $h\gtrsim 0.033$ ($h\gtrsim 0.057$). We suggest that their long lifetimes are sustained by the interaction of spiral density waves launched by the neighboring planets. Overall, our results show that distinguishing imprint of compact systems with Neptune-mass planets are long-lived vortices inside the density gaps, which in turn are shallower than single-planet gaps for a fixed gap width.
△ Less
Submitted 19 April, 2022;
originally announced April 2022.
-
Dust trapping around Lagrangian points in protoplanetary disks
Authors:
Matias Montesinos,
Juan Garrido-Deutelmoser,
Johan Olofsson,
Cristian A. Giuppone,
Jorge Cuadra,
Amelia Bayo,
Mario Sucerquia,
Nicolas Cuello
Abstract:
Trojans are defined as objects that share the orbit of a planet at the stable Lagrangian points $L_4$ and $L_5$. In the Solar System, these bodies show a broad size distribution ranging from micrometer($μ$m) to centimeter(cm) particles (Trojan dust) and up to kilometer (km) rocks (Trojan asteroids). It has also been theorized that earth-like Trojans may be formed in extra-solar systems. The Trojan…
▽ More
Trojans are defined as objects that share the orbit of a planet at the stable Lagrangian points $L_4$ and $L_5$. In the Solar System, these bodies show a broad size distribution ranging from micrometer($μ$m) to centimeter(cm) particles (Trojan dust) and up to kilometer (km) rocks (Trojan asteroids). It has also been theorized that earth-like Trojans may be formed in extra-solar systems. The Trojan formation mechanism is still under debate, especially theories involving the effects of dissipative forces from a viscous gaseous environment. We perform hydro-simulations to follow the evolution of a protoplanetary disk with an embedded 1--10 Jupiter-mass planet. On top of the gaseous disk, we set a distribution of $μ$m--cm dust particles interacting with the gas. This allows us to follow dust dynamics as solids get trapped around the Lagrangian points of the planet. We show that large vortices generated at the Lagrangian points are responsible for dust accumulation, where the leading Lagrangian point $L_4$ traps a larger amount of submillimeter (submm) particles than the trailing $L_5$, which traps mostly mm--cm particles. However, the total bulk mass, with typical values of $\sim M_{\rm moon}$, is more significant in $L_5$ than in $L_4$, in contrast to what is observed in the current Solar System a few gigayears later. Furthermore, the migration of the planet does not seem to affect the reported asymmetry between $L_4$ and $L_5$. The main initial mass reservoir for Trojan dust lies in the same co-orbital path of the planet, while dust migrating from the outer region (due to drag) contributes very little to its final mass, imposing strong mass constraints for the in situ formation scenario of Trojan planets.
△ Less
Submitted 12 October, 2020; v1 submitted 22 September, 2020;
originally announced September 2020.
