Planet-disc interactions with Discontinuous Galerkin Methods using GPUs

Abstract

We present a two-dimensional Cartesian code based on high order discontinuous Galerkin methods, implemented to run in parallel over multiple GPUs. A simple planet-disc setup is used to compare the behaviour of our code against the behaviour found using the FARGO3D code with a polar mesh. We make use of the time dependence of the torque exerted by the disc on the planet as a mean to quantify the numerical viscosity of the code. We find that the numerical viscosity of the Keplerian flow can be as low as a few $10^{-8}{r}^2\Omega$, $r$ and $\Omega$ being respectively the local orbital radius and frequency, for fifth order schemes and resolution of $\sim 10^{-2}r$. Although for a single disc problem a solution of low numerical viscosity can be obtained at lower computational cost with FARGO3D (which is nearly an order of magnitude faster than a fifth order method), discontinuous Galerkin methods appear promising to obtain solutions of low numerical viscosity in more complex situations where the flow cannot be captured on a polar or spherical mesh concentric with the disc.

arXiv commentAccepted for publication in MNRAS

References in corpus (26)

• E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh
• On the width and shape of gaps in protoplanetary disks
• FARGO: a fast eulerian transport algorithm for differentially rotating disks
• A torque formula for non-isothermal Type I planetary migration - II. Effects of diffusion
• Wind-driven Accretion in Protoplanetary Disks. I: Suppression of the Magnetorotational Instability and Launching of the Magnetocentrifugal Wind
• A comparative study of disc-planet interaction
• Smoothed Particle Hydrodynamics in Astrophysics
• Thanatology in Protoplanetary Discs: the combined influence of Ohmic, Hall, and ambipolar diffusion on dead zones
• Orbital Migration and Mass Accretion of Protoplanets in 3D Global Computations with Nested Grids
• FARGO3D: A new GPU-oriented MHD code
• How Empty are Disk Gaps Opened by Giant Planets?
• Turbulence in the TW Hya Disk
• Standing on the shoulders of giants: Trojan Earths and vortex trapping in low mass self-gravitating protoplanetary disks of gas and solids
• Particle Concentration At Planet Induced Gap Edges and Vortices: I. Inviscid 3-D Hydro Disks
• On the coorbital corotation torque in a viscous disk and its impact on planetary migration
• Improved torque formula for low and intermediate mass planetary migration
• Saturated torque formula for planetary migration in viscous disks with thermal diffusion: recipe for protoplanet population synthesis
• A conservative orbital advection scheme for simulations of magnetized shear flows with the PLUTO code
• Protoplanetary Disks as (Possibly) Viscous Disks
• Gas and dust hydrodynamical simulations of massive lopsided transition discs - I. Gas distribution
• Astrophysical hydrodynamics with a high-order discontinuous Galerkin scheme and adaptive mesh refinement
• DISCO: A 3D Moving-Mesh Magnetohydrodynamics Code Designed for the Study of Astrophysical Disks
• Numerical simulations of the type III migration:I. Disc model and convergence tests
• RODEO: a new method for planet-disk interaction
• Spiral-driven accretion in protoplanetary discs - III tri-dimensional simulations
• Modeling gravitational instabilities in self-gravitating protoplanetary disks with adaptive mesh refinement techniques

Cited by in corpus (7)

• High-order Magnetohydrodynamics for Astrophysics with an Adaptive Mesh Refinement Discontinuous Galerkin Scheme
• High-order Discontinuous Galerkin hydrodynamics with sub-cell shock capturing on GPUs
• Spectral Difference method with a posteriori limiting: Application to the Euler equations in one and two space dimensions
• WENO-Wombat: Scalable Fifth-Order Constrained-Transport Magnetohydrodynamics for Astrophysical Applications
• Supersonic turbulence simulations with GPU-based high-order Discontinuous Galerkin hydrodynamics
• Spectral Difference method with a posteriori limiting: II- Application to low Mach number flows
• Capturing near-equilibrium solutions: a comparison between high-order discontinuous Galerkin methods and well-balanced schemes