Planet-disc interactions with discontinuous Galerkin methods using GPUs

David A.Velasco Romero, Maria Han Veiga, Romain Teyssier, Frédéric S. Masset

Research output: Contribution to journalArticlepeer-review

5 Scopus citations


We present a two-dimensional Cartesian code based on high-order discontinuous Galerkin methods, implemented to run in parallel over multiple graphics processing units. 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-8r2Ω, r and Ω being respectively the local orbital radius and frequency, for fifth-order schemes and resolution of ~10-2r. 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.

Original languageEnglish (US)
Pages (from-to)1855-1865
Number of pages11
JournalMonthly Notices of the Royal Astronomical Society
Issue number2
StatePublished - Aug 1 2018
Externally publishedYes

All Science Journal Classification (ASJC) codes

  • Astronomy and Astrophysics
  • Space and Planetary Science


  • Hydrodynamics
  • Methods: numerical
  • Planet-disc interactions
  • Protoplanetary discs


Dive into the research topics of 'Planet-disc interactions with discontinuous Galerkin methods using GPUs'. Together they form a unique fingerprint.

Cite this