Supersonic turbulence simulations with GPU-based high-order Discontinuous Galerkin hydrodynamics
arXiv:2401.06841 · doi:10.1093/mnras/stae2192
Abstract
We investigate the numerical performance of a Discontinuous Galerkin (DG) hydrodynamics implementation when applied to the problem of driven, isothermal supersonic turbulence. While the high-order element-based spectral approach of DG is known to efficiently produce accurate results for smooth problems (exponential convergence with expansion order), physical discontinuities in solutions, like shocks, prove challenging and may significantly diminish DG's applicability to practical astrophysical applications. We consider whether DG is able to retain its accuracy and stability for highly supersonic turbulence, characterized by a network of shocks. We find that our new implementation, which regularizes shocks at sub-cell resolution with artificial viscosity, still performs well compared to standard second-order schemes for moderately high Mach number turbulence, provided we also employ an additional projection of the primitive variables onto the polynomial basis to regularize the extrapolated values at cell interfaces. However, the accuracy advantage of DG diminishes significantly in the highly supersonic regime. Nevertheless, in turbulence simulations with a wide dynamic range that start with supersonic Mach numbers and can resolve the sonic point, the low numerical dissipation of DG schemes still proves advantageous in the subsonic regime. Our results thus support the practical applicability of DG schemes for demanding astrophysical problems that involve strong shocks and turbulence, such as star formation in the interstellar medium. We also discuss the substantial computational cost of DG when going to high order, which needs to be weighted against the resulting accuracy gain. For problems containing shocks, this favours the use of comparatively low DG order.
22 pages, 14 figures, accepted for publication in MNRAS
References in corpus (32)
- E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh
- GADGET: A code for collisionless and gasdynamical cosmological simulations
- Athena: A New Code for Astrophysical MHD
- Cosmological SPH simulations: The entropy equation
- Comparing the statistics of interstellar turbulence in simulations and observations: Solenoidal versus compressive turbulence forcing
- The Density Probability Distribution in Compressible Isothermal Turbulence: Solenoidal versus Compressive Forcing
- The Arepo public code release
- One-Point Probability Distribution Functions of Supersonic Turbulent Flows in Self-Gravitating Media
- On the universality of supersonic turbulence
- On the Origin of the Global Schmidt Law of Star Formation
- The Fractal Density Structure in Supersonic Isothermal Turbulence: Solenoidal versus Compressive Energy Injection
- Kinetic and Structural Evolution of Self-gravitating, Magnetized Clouds: 2.5-Dimensional Simulations of Decaying Turbulence
- The Formation of Stellar Clusters in Turbulent Molecular Clouds: Effects of the Equation of State
- Effects of the Equation of State on the Formation of Star Clusters
- Physical vs. Observational Properties of Clouds in Turbulent Molecular Cloud Models
- Subsonic turbulence in smoothed particle hydrodynamics and moving-mesh simulations
- The sonic scale revealed by the world's largest supersonic turbulence simulation
- A localised subgrid scale model for fluid dynamical simulations in astrophysics I: Theory and numerical tests
- The Distribution of Pressures in a Supernova-Driven Interstellar Medium. I. Magnetized Medium
- A comparison between grid and particle methods on the statistics of driven, supersonic, isothermal turbulence
- Astrophysical turbulence modeling
- From diffuse gas to dense molecular cloud cores
- Astrophysical hydrodynamics with a high-order discontinuous Galerkin scheme and adaptive mesh refinement
- Forming early-type galaxies without AGN feedback: a combination of merger-driven outflows and inefficient star formation
- Density profile evolution during prestellar core collapse: Collapse starts at the large scale
- Early Science from POSSUM: Shocks, turbulence, and a massive new reservoir of ionised gas in the Fornax cluster
- High-order Magnetohydrodynamics for Astrophysics with an Adaptive Mesh Refinement Discontinuous Galerkin Scheme
- ALMA observations of NGC 6334S I: Forming massive stars and cluster in subsonic and transonic filamentary clouds
- The role of the turbulence driving mode for the Initial Mass Function
- Planet-disc interactions with Discontinuous Galerkin Methods using GPUs
- High-order Discontinuous Galerkin hydrodynamics with sub-cell shock capturing on GPUs
- Density and Velocity Correlations in Isothermal Supersonic Turbulence