18 research outputs found
FARGO3D: A new GPU-oriented MHD code
We present the FARGO3D code, recently publicly released. It is a magnetohydrodynamics code developed with special emphasis on the physics of protoplanetary disks and planet-disk interactions, and parallelized with MPI. The hydrodynamics algorithms are based on finite-difference upwind, dimensionally split methods. The magnetohydrodynamics algorithms consist of the constrained transport method to preserve the divergence-free property of the magnetic field to machine accuracy, coupled to a method of characteristics for the evaluation of electromotive forces and Lorentz forces. Orbital advection is implemented, and an N-body solver is included to simulate planets or stars interacting with the gas. We present our implementation in detail and present a number of widely known tests for comparison purposes. One strength of FARGO3D is that it can run on either graphical processing units (GPUs) or central processing units (CPUs), achieving large speed-up with respect to CPU cores. We describe our implementation choices, which allow a user with no prior knowledge of GPU programming to develop new routines for CPUs, and have them translated automatically for GPUs.Fil: Benítez Llambay, Pablo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Astronomía Teórica y Experimental. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba. Instituto de Astronomía Teórica y Experimental; ArgentinaFil: Masset, Frédéric S.. Universidad Nacional Autónoma de México; Méxic
Horseshoe drag in three-dimensional globally isothermal disks
We study the horseshoe dynamics of a low-mass planet in a three-dimensional, globally isothermal, inviscid disk. We find, as reported in previous work, that the boundaries of the horseshoe region (separatrix sheets) have cylindrical symmetry about the disk´s rotation axis. We interpret this feature as arising from the fact that the whole separatrix sheets have a unique value of Bernoulli´s constant, and that this constant does not depend on altitude, but only on the cylindrical radius, in barotropic disks. We next derive an expression for the torque exerted by the horseshoe region on the planet, or horseshoe drag. Potential vorticity is not materially conserved as in two-dimensional flows, but it obeys a slightly more general conservation law (Ertel´s theorem) that allows an expression for the horseshoe drag identical to the expression in a two-dimensional disk to be obtained. Our results are illustrated and validated by three-dimensional numerical simulations. The horseshoe region is found to be slightly narrower than previously extrapolated from two-dimensional analyses with a suitable softening length of the potential. We discuss the implications of our results for the saturation of the corotation torque, and the possible connection to the flow at the Bondi scale, which the present analysis does not resolve.Fil: Masset, F. S.. Universidad Nacional Autónoma de México; MéxicoFil: Benítez Llambay, Pablo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Astronomía Teórica y Experimental. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba. Instituto de Astronomía Teórica y Experimental; Argentin
Cálculo de las alturas del relieve lunar
Con el presente trabajo se pretende dar un método sencillo y exacto para calcular la altura de algún accidente lunar. Como nos ha sido imposible encontrar un trabajo que resuma en una ecuación las variables más importantes a tener en cuenta en este proceso, creímos necesario desarrollar un método lo suficientemente sencillo que permita a cualquier persona calcular la altura de los accidentes lunares a partir de simples mediciones ópticas, aún careciendo del instrumental adecuado. Es importante destacar el hecho de que contando con una ecuación que exprese la altura de cualquier relieve en función de otras variables (que pueden ser medidas) es posible hacer uso de la computación para el desarrollo de un software que facilite dicho cálculo. El trabajo toma como punto de partida conocimientos de trigonometría muy sencillos. Sabemos que a partir de una relación trigonométrica es posible calcular la altura de un objeto conociendo la longitud de su sombra y el ángulo que los rayos del sol forman con la horizontal en ese momento
Formación y evolución de exoplanetas
Tesis (Lic. en Astronomía)--Universidad Nacional de Córdoba. Facultad de Matemática, Astronomía y Física, 2011.En este trabajo se estudian principalmente los procesos de formación y evolución planetaria en diferentes etapas. En la etapa temprana, se estudia el proceso desde la formación del disco protoplanetario hasta la formación del embrión planetario. Se estudia la interacción planeta-disco y las consecuencias de ésta. Se estudia el concepto de trampas planetarias y su aplicación a problemas concretos. En la etapa tardía se estudian las fuerzas de marea generadas por la interacción gravitatoria planeta-estrella, y en particular se aplican estos resultados a planetas de corto período. Como aplicación se realiza un estudio de la distribución masa-período de los exoplanetas de corto período conocidos y se propone una explicación para la forma de tal distribución. Por último, se aplican los estudios realizados a un caso exoplanetario particular: CoRoT-7.Pablo Benítez Llamba
Planet heating prevents inward migration of planetary cores
Planetary systems are born in the disks of gas, dust and rocky fragments that
surround newly formed stars. Solid content assembles into ever-larger rocky
fragments that eventually become planetary embryos. These then continue their
growth by accreting leftover material in the disc. Concurrently, tidal effects
in the disc cause a radial drift in the embryo orbits, a process known as
migration. Fast inward migration is predicted by theory for embryos smaller
than three to five Earth masses. With only inward migration, these embryos can
only rarely become giant planets located at Earth's distance from the Sun and
beyond, in contrast with observations. Here we report that asymmetries in the
temperature rise associated with accreting infalling material produce a force
(which gives rise to an effect that we call "heating torque") that counteracts
inward migration. This provides a channel for the formation of giant planets
and also explains the strong planet-metallicity correlation found between the
incidence of giant planets and the heavy-element abundance of the host stars.Comment: 19 pages, 4 figure
Long-term and large-scale hydrodynamical simulations of migrating planets
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.Comment: Accepted for publication in ApJ. For a movie describing the method,
see https://youtu.be/66o0Z2lX8N
Observability of planet-disc interactions in CO kinematics
Empirical evidence of planets in gas-rich circumstellar discs is required to
constrain giant planet formation theories. Here we study the kinematic patterns
which arise from planet-disc interactions and their observability in CO
rotational emission lines. We perform three-dimensional hydrodynamical
simulations of single giant planets, and predict the emergent intensity field
with radiative transfer. Pressure gradients at planet-carved gaps, spiral wakes
and vortices bear strong kinematic counterparts. The iso-velocity contours in
the CO(2-1) line centroids reveal large-scale perturbations,
corresponding to abrupt transitions from below sub-Keplerian to super-Keplerian
rotation along with radial and vertical flows. The increase in line optical
depth at the edge of the gap also modulates , but this is a mild
effect compared to the dynamical imprint of the planet-disc interaction. The
large-scale deviations from the Keplerian rotation thus allow the planets to be
indirectly detected via the first moment maps of molecular gas tracers, at ALMA
angular resolutions. The strength of these deviations depends on the mass of
the perturber. This initial study paves the way to eventually determine the
mass of the planet by comparison with more detailed models.Comment: 6 pages, 3 color figures. 1 animation (Figure 3, Adobe Reader
recommended). Accepted for publication in MNRAS Letter
Dust Settling Instability in Protoplanetary Discs
The streaming instability (SI) has been extensively studied in the linear and
non-linear regimes as a mechanism to concentrate solids and trigger
planetesimal formation in the midplane of protoplanetary discs. A related dust
settling instability (DSI) applies to particles while settling towards the
midplane. The DSI has previously been studied in the linear regime, with
predictions that it could trigger particle clumping away from the midplane.
This work presents a range of linear calculations and non-linear simulations,
performed with FARGO3D, to assess conditions for DSI growth. We expand on
previous linear analyses by including particle size distributions and
performing a detailed study of the amount of background turbulence needed to
stabilize the DSI. When including binned size distributions, the DSI often
produces converged growth rates with fewer bins than the standard SI. With
background turbulence, we find that the most favorable conditions for DSI
growth are weak turbulence, characterized by with
intermediate-sized grains that settle from one gas scale-height. These
conditions could arise during a sudden decrease in disc turbulence following an
accretion outburst. Ignoring background turbulence, we performed a parameter
survey of local 2D DSI simulations. Particle clumping was either weak or
occurred slower than particles settle. Clumping was reduced by a factor of two
in a comparison 3D simulation. Overall, our results strongly disfavor the
hypothesis that the DSI significantly promotes planetesimal formation.
Non-linear simulations of the DSI with different numerical methods could
support or challenge these findings
RAM: Rapid Advection Algorithm on Arbitrary Meshes
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.Comment: 15 pages, 7 figures. Submitted to ApJ. Comments are welcom