THC: a new high-order finite-difference high-resolution shock-capturing code for special-relativistic hydrodynamics

By David Radice, Luciano Rezzolla. 2012

151 citations recorded by the source corpus.

We present THC: a new high-order flux-vector-splitting code for Newtonian and special-relativistic hydrodynamics designed for direct numerical simulations of turbulent flows. Our code implements a variety of different reconstruction algorithms, such as the popular weighted essentially non oscillatory and monotonicity-preserving schemes, or the more specialised bandwidth-optimised WENO scheme that has been specifically designed for the study of compressible turbulence. We show the first systematic comparison of these schemes in Newtonian physics as well as for special-relativistic flows. In particular we will present the results obtained in simulations of grid-aligned and oblique shock waves and nonlinear, large-amplitude, smooth adiabatic waves. We will also discuss the results obtained in classical benchmarks such as the double-Mach shock reflection test in Newtonian physics or the linear and nonlinear development of the relativistic Kelvin-Helmholtz instability in two and three dimensions. Finally, we study the turbulent flow induced by the Kelvin-Helmholtz instability and we show that our code is able to obtain well-converged velocity spectra, from which we benchmark the effective resolution of the different schemes.

Subfields: Astronomy and Astrophysics, Computational Mechanics

View Paper: https://www.semanticscholar.org/paper/0786aba27adcf942b12e38d728a4dee015502ca9

How researchers see this paper

The researchers whose trusted networks rank it highest, strongest advocate first. Every researcher gets a different order.

In Astronomy and Astrophysics

In Computational Mechanics

Citations in the corpus

Beyond second-order convergence in simulations of binary neutron stars in full general-relativity (2013, methodology): “It is smaller than the nominal one of the scheme (which is fourth since we use fourth-order finite-differencing for the spacetim ), but this is to be expected because HRSC methods typically reach their nominal convergence order only at very high resolutions (Shu 1997; Radice & Rezzolla 2012).”
GW170817: Joint Constraint on the Neutron Star Equation of State from Multimessenger Observations (2017, methodology): “We perform 29 merger simulations using the GR hydrodynamics code WhiskyTHC (Radice & Rezzolla 2012; Radice et al. 2014a, 2014b).”
Binary Neutron Star Mergers: Mass Ejection, Electromagnetic Counterparts, and Nucleosynthesis (2018, methodology): “We evolve the initial data using the WhiskyTHC code (Radice & Rezzolla 2012; Radice et al. 2014a, 2014b, 2015).”
The black hole accretion code (2016, background): “f In the general-relativistic hydrodynamic WhiskyTHC code (Radice and Rezzolla 2012; Radice et al. 2014), this desirable property allows to set floors on density close to the limit of floating point precision ∼ 10ρref .”