-
RAM: Rapid Advection Algorithm on Arbitrary Meshes
Authors:
Pablo Benítez-Llambay,
Leonardo Krapp,
Ximena S. Ramos,
Kaitlin M. Kratter
Abstract:
The study of many astrophysical flows requires computational algorithms that can capture high Mach number flows, while resolving a large dynamic range in spatial and density scales. In this paper we present a novel method, RAM: Rapid Advection Algorithm on Arbitrary Meshes. RAM is a time-explicit method to solve the advection equation in problems with large bulk velocity on arbitrary computational…
▽ More
The study of many astrophysical flows requires computational algorithms that can capture high Mach number flows, while resolving a large dynamic range in spatial and density scales. In this paper we present a novel method, RAM: Rapid Advection Algorithm on Arbitrary Meshes. RAM is a time-explicit method to solve the advection equation in problems with large bulk velocity on arbitrary computational grids. In comparison with standard up-wind algorithms, RAM enables advection with larger time steps and lower truncation errors. Our method is based on the operator splitting technique and conservative interpolation. Depending on the bulk velocity and resolution, RAM can decrease the numerical cost of hydrodynamics by more than one order of magnitude. To quantify the truncation errors and speed-up with RAM, we perform one and two-dimensional hydrodynamics tests. We find that the order of our method is given by the order of the conservative interpolation and that the effective speed up is in agreement with the relative increment in time step. RAM will be especially useful for numerical studies of disk-satellite interaction, characterized by high bulk orbital velocities, and non-trivial geometries. Our method dramatically lowers the computational cost of simulations that simultaneously resolve the global disk and well inside the Hill radius of the secondary companion.
△ Less
Submitted 9 May, 2023;
originally announced May 2023.
-
Migration of Jupiter mass planets in discs with laminar accretion flows
Authors:
E. Lega,
A. Morbidelli,
R. P. Nelson,
X. S. Ramos,
A. Crida,
W. Bethune,
K. Batygin
Abstract:
Migration of giant planets in discs with low viscosity has been studied recently. The proportionality between migration speed and the disc's viscosity is broken by the presence of vortices that appear at the edges of the planet-induced gap. Our goal is to investigate vortex-driven migration in low-viscosity discs in the presence of radial advection of gas, as expected from angular momentum removal…
▽ More
Migration of giant planets in discs with low viscosity has been studied recently. The proportionality between migration speed and the disc's viscosity is broken by the presence of vortices that appear at the edges of the planet-induced gap. Our goal is to investigate vortex-driven migration in low-viscosity discs in the presence of radial advection of gas, as expected from angular momentum removal by magnetised disc winds. We performed three dimensional simulations using the grid-based code FARGOCA. We mimicked the effects of a disc wind by applying a synthetic torque on a surface layer of the disc characterised by a prescribed column density Sigma_A so that it results in a disc accretion rate of 10^-8 Solar masses per year. Discs with this structure are called 'layered' and the layer where the torque is applied is denoted as 'active'. We also consider the case of accretion focussed near the disc midplane to mimic transport properties induced by a large Hall effect or by weak Ohmic diffusion. We observe two migration phases: in the first phase, the migration of the planet is driven by the vortex and is directed inwards. This phase ends when the vortex disappear. Migration depends on the ability of the torque from the planet to block the accretion flow. When the flow is fast and unimpeded (small Sigma_A) migration is very slow. When the accretion flow is completely blocked, migration is faster and the speed is controlled by the rate at which the accretion flow refills the gap behind the migrating planet. The migration speed of a giant planet in a layered protoplanetary disc depends on the thickness of the accreting layer. The lack of large-scale migration apparently experienced by the majority of giant exoplanets can be explained if the accreting layer is sufficiently thin to allow unimpeded accretion through the disc.
△ Less
Submitted 8 November, 2021;
originally announced November 2021.
-
The dynamical evolution of escaped Jupiter Trojan asteroids, link to other minor body populations
Authors:
Romina P. Di Sisto,
Ximena S. Ramos,
Tabaré Gallardo
Abstract:
Jupiter Trojan asteroids are located around L4 and L5 Lagrangian points on relatively stable orbits, in 1:1 MMR with Jupiter. However, not all of them lie in orbits that remain stable over the age of the Solar System. Unstable zones allow some Trojans to escape in time scales shorter than the Solar System age. This may contribute to populate other small body populations. In this paper, we study th…
▽ More
Jupiter Trojan asteroids are located around L4 and L5 Lagrangian points on relatively stable orbits, in 1:1 MMR with Jupiter. However, not all of them lie in orbits that remain stable over the age of the Solar System. Unstable zones allow some Trojans to escape in time scales shorter than the Solar System age. This may contribute to populate other small body populations. In this paper, we study this process by performing long-term numerical simulations of the observed Trojans, focusing on the trajectories of those that leave the resonance. The orbits of current Trojans are taken as initial conditions and their evolution is followed under the gravitational action of the Sun and the planets. We find the rate of escape of Trojans from L5, ~1.1 times greater than from L4. The majority of escaped Trojans have encounters with Jupiter although they have encounters with the other planets too. Almost all escaped Trojans reach the comet zone, ~90% cross the Centaur zone and only L4 Trojans reach the transneptunian zone. Considering the real asymmetry between L4 and L5, we show that 18 L4 Trojans and 14 L5 Trojans with diameter D > 1 km are ejected from the resonance every Myr. The contribution of the escaped Trojans to other minor body populations would be negligible, being the contribution from L4 and L5 to JFCs and no-JFCs almost the same, and the L4 contribution to Centaurs and TNOs, orders of magnitude greater than that of L5. Considering the collisional removal, besides the dynamical one, and assuming that Trojans that escape due to collisions follow the same dynamical behavior that the ones removed by dynamics, we would have a minor contribution of Trojans to comets and Centaurs. However, there would be some specific regions were escaped Trojans could be important such as ACOs, Encke-type comets, S-L 9-type impacts on Jupiter and NEOs.
△ Less
Submitted 1 November, 2018;
originally announced November 2018.
-
Orbital migration and Resonance Offset of the Kepler-25 and K2-24 systems
Authors:
C. Charalambous,
X. S. Ramos,
P. Benítez-Llambay,
C. Beaugé
Abstract:
Based on the model described in Ramos et al., 2017, we present an analytical+numerical study of the resonance capture under Type-I migration for the Kepler-25 (Marcy et al., 2014) and K2-24 (Petigura et al., 2016) Kepler systems, both close to a 2/1 mean-motion resonance. We find that, depending on the flare index and the proximity to the central star, the average value of the period-ratio between…
▽ More
Based on the model described in Ramos et al., 2017, we present an analytical+numerical study of the resonance capture under Type-I migration for the Kepler-25 (Marcy et al., 2014) and K2-24 (Petigura et al., 2016) Kepler systems, both close to a 2/1 mean-motion resonance. We find that, depending on the flare index and the proximity to the central star, the average value of the period-ratio between two consecutive planets show a significant deviation with respect to the resonant nominal value, up to values well in agreement with the observations.
△ Less
Submitted 15 March, 2018;
originally announced March 2018.
-
Resonance Capture and Dynamics of 3-Planet Systems
Authors:
C. Charalambous,
J. G. Martí,
C. Beaugé,
X. S. Ramos
Abstract:
We present a series of dynamical maps for fictitious 3-planets systems in initially circular coplanar orbits. These maps have unveiled a rich resonant structure involving two or three planets, as well as indicating possible migration routes from secular to double resonances or pure 3-planet commensurabilities. These structures are then compared to the present-day orbital architecture of observed r…
▽ More
We present a series of dynamical maps for fictitious 3-planets systems in initially circular coplanar orbits. These maps have unveiled a rich resonant structure involving two or three planets, as well as indicating possible migration routes from secular to double resonances or pure 3-planet commensurabilities. These structures are then compared to the present-day orbital architecture of observed resonant chains. In a second part of the paper we describe N-body simulations of type-I migration. Depending on the orbital decay timescale, we show that 3-planet systems may be trapped in different combinations of independent commensurabilities: (i) double resonances, (ii) intersection between a 2-planet and a first-order 3-planet resonance, and (iii) simultaneous libration in two first-order 3-planet resonances. These latter outcomes are found for slow migrations, while double resonances are almost always the final outcome in high-density disks. Finally, we discuss an application to the TRAPPIST-1 system. We find that, for low migration rates and planetary masses of the order of the estimated values, most 3-planet sub-systems are able to reach the observed double resonances after following evolutionary routes defined by pure 3-planet resonances. The final orbital configuration shows resonance offsets comparable with present-day values without the need of tidal dissipation. For the 8/5 resonance proposed to dominate the dynamics of the two inner planets, we find little evidence of its dynamical significance; instead, we propose that this relation between mean motions could be a consequence of the interaction between a pure 3-planet resonance and a 2-planet commensurability between planets c and d.
△ Less
Submitted 14 March, 2018;
originally announced March 2018.
-
Planetary migration and the origin of the 2:1 and 3:2 (near)-resonant population of close-in exoplanets
Authors:
X. S. Ramos,
C. Charalambous,
P. Benítez-Llambay,
C. Beaugé
Abstract:
We present an analytical and numerical study of the orbital migration and resonance capture of fictitious two-planet systems with masses in the super-Earth range undergoing Type-I migration. We find that, depending on the flare index and proximity to the central star, the average value of the period ratio, $P_2/P_1$, between both planets may show a significant deviation with respect to the nominal…
▽ More
We present an analytical and numerical study of the orbital migration and resonance capture of fictitious two-planet systems with masses in the super-Earth range undergoing Type-I migration. We find that, depending on the flare index and proximity to the central star, the average value of the period ratio, $P_2/P_1$, between both planets may show a significant deviation with respect to the nominal value. For planets trapped in the 2:1 commensurability, offsets may reach values on the order of $0.1$ for orbital periods on the order of $1$ day, while systems in the 3:2 mean-motion resonance (MMR) show much smaller offsets for all values of the semimajor axis. These properties are in good agreement with the observed distribution of near-resonant exoplanets, independent of their detection method. We show that 2:1-resonant systems far from the star, such as HD82943 and HR8799, are characterized by very small resonant offsets, while higher values are typical of systems discovered by Kepler with orbital periods approximately a few days. Conversely, planetary systems in the vicinity of the 3:2 MMR show little offset with no significant dependence on the orbital distance. In conclusion, our results indicate that the distribution of Kepler planetary systems around the 2:1 and 3:2 MMR are consistent with resonant configurations obtained as a consequence of a smooth migration in a laminar flared disk, and no external forces are required to induce the observed offset or its dependence with the commensurability or orbital distance from the star.
△ Less
Submitted 21 April, 2017;
originally announced April 2017.
-
Long-term and large-scale hydrodynamical simulations of migrating planets
Authors:
Pablo Benítez-Llambay,
Ximena S. Ramos,
Cristian Beaugé,
Frédéric Masset
Abstract:
We present a new method that allows long-term and large-scale hydrodynamical simulations of migrating planets over a grid-based Eulerian code. This technique, which consists in a remapping of the disk by tracking the planetary migration, enables runs of migrating planets over a time comparable to the age of protoplanetary disks. This method also has the potential to address efficiently problems re…
▽ More
We present a new method that allows long-term and large-scale hydrodynamical simulations of migrating planets over a grid-based Eulerian code. This technique, which consists in a remapping of the disk by tracking the planetary migration, enables runs of migrating planets over a time comparable to the age of protoplanetary disks. This method also has the potential to address efficiently problems related with migration of multi-planet systems in gaseous disks, and to improve current results of migration of massive planets by including global viscous evolution as well as detailed studies of the co-orbital region during migration. We perform different tests using the public code FARGO3D to validate this method and compare its results with those obtained using a classical fixed grid.
△ Less
Submitted 5 May, 2016;
originally announced May 2016.
-
Giga-Year Evolution of Jupiter Trojans and the Asymmetry Problem
Authors:
Romina P. Di Sisto,
Ximena S. Ramos,
Cristián Beaugé
Abstract:
We present a series of numerical integrations of observed and fictitious Jupiter Trojan asteroids, under the gravitational effects of the four outer planets, for time-spans comparable with the age of the Solar System. From these results we calculate the escape rate from each Lagrange point, and construct dynamical maps of "permanence" time in different regions of the phase space. Fictitious astero…
▽ More
We present a series of numerical integrations of observed and fictitious Jupiter Trojan asteroids, under the gravitational effects of the four outer planets, for time-spans comparable with the age of the Solar System. From these results we calculate the escape rate from each Lagrange point, and construct dynamical maps of "permanence" time in different regions of the phase space. Fictitious asteroids in L4 and L5 show no significant difference, showing almost identical dynamical maps and escape rates. For real Trojans, however, we found that approximately 23% o f the members of the leading swarm escaped after 4.5 Gyrs, while this number increased to 28.3% for L5. This implies that the asymmetry between the two populations increases with time, indicating that it may have been smaller at the time of formation/capture of these asteroids. Nevertheless, the difference in chaotic diffusion cannot, in itself, account for the current observed asymmetry (~40%), and must be primarily primordial and characteristic of the capture mechanism of the Trojans. Finally, we calculate new proper elements for all the numbered Trojans using the semi-analytical approach of Beaugé and Roig (2001), and compare the results with the numerical estimations by Brož and Rosehnal (2011). For asteroids that were already numbered in 2011, both methods yield very similar results, while significant differences were found for those bodies that became numbered after 2011.
△ Less
Submitted 18 April, 2016;
originally announced April 2016.
-
The Resonance Overlap and Hill Stability Criteria Revisited
Authors:
X. S. Ramos,
J. A. Correa-Otto,
C. Beaugé
Abstract:
We review the orbital stability of the planar circular restricted three-body problem, in the case of massless particles initially located between both massive bodies. We present new estimates of the resonance overlap criterion and the Hill stability limit, and compare their predictions with detailed dynamical maps constructed with N-body simulations. We show that the boundary between (Hill) stable…
▽ More
We review the orbital stability of the planar circular restricted three-body problem, in the case of massless particles initially located between both massive bodies. We present new estimates of the resonance overlap criterion and the Hill stability limit, and compare their predictions with detailed dynamical maps constructed with N-body simulations. We show that the boundary between (Hill) stable and unstable orbits is not smooth but characterized by a rich structure generated by the superposition of different mean-motion resonances which does not allow for a simple global expression for stability.
We propose that, for a given perturbing mass $m_1$ and initial eccentricity $e$, there are actually two critical values of the semimajor axis. All values $a < a_{\rm Hill}$ are Hill-stable, while all values $a > a_{\rm unstable}$ are unstable in the Hill sense. The first limit is given by the Hill-stability criterion and is a function of the eccentricity. The second limit is virtually insensitive to the initial eccentricity, and closely resembles a new resonance overlap condition (for circular orbits) developed in terms of the intersection between first and second-order mean-motion resonances.
△ Less
Submitted 11 September, 2015;
originally announced September 2015.
