Header

UZH-Logo

Maintenance Infos

Planet–disc interactions with discontinuous Galerkin methods using GPUs


Velasco Romero, David A; Han Veiga, Maria; Teyssier, Romain; Masset, Frédéric S (2018). Planet–disc interactions with discontinuous Galerkin methods using GPUs. Monthly Notices of the Royal Astronomical Society, 478(2):1855-1865.

Abstract

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.

Abstract

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.

Statistics

Citations

Dimensions.ai Metrics
2 citations in Web of Science®
1 citation in Scopus®
Google Scholar™

Altmetrics

Downloads

5 downloads since deposited on 01 Mar 2019
5 downloads since 12 months
Detailed statistics

Additional indexing

Item Type:Journal Article, refereed, original work
Communities & Collections:07 Faculty of Science > Institute for Computational Science
Dewey Decimal Classification:530 Physics
Uncontrolled Keywords:Space and Planetary Science, Astronomy and Astrophysics
Language:English
Date:2018
Deposited On:01 Mar 2019 13:36
Last Modified:17 Sep 2019 19:25
Publisher:Oxford University Press
ISSN:0035-8711
OA Status:Green
Free access at:Publisher DOI. An embargo period may apply.
Publisher DOI:https://doi.org/10.1093/mnras/sty1192

Download

Download PDF  'Planet–disc interactions with discontinuous Galerkin methods using GPUs'.
Preview
Content: Published Version
Filetype: PDF
Size: 1MB
View at publisher