de.mpg.escidoc.pubman.appbase.FacesBean
English
 
Help Guide Disclaimer Contact us Login
  Advanced SearchBrowse

Item

ITEM ACTIONSEXPORT

Released

Journal Article

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

MPS-Authors

Radice,  D.
Astrophysical Relativity, AEI-Golm, MPI for Gravitational Physics, Max Planck Society;

http://pubman.mpdl.mpg.de/cone/persons/resource/persons20670

Rezzolla,  L.
Astrophysical Relativity, AEI-Golm, MPI for Gravitational Physics, Max Planck Society;

Locator
There are no locators available
Fulltext (public)

aa19735-12.pdf
(Any fulltext), 5MB

Supplementary Material (public)
There is no public supplementary material available
Citation

Radice, D., & Rezzolla, L. (2012). THC: a new high-order finite-difference high-resolution shock-capturing code for special-relativistic hydrodynamics. Astronomy and Astrophysics, 547: A26. doi:http://dx.doi.org/10.1051/0004-6361/201219735.


Cite as: http://hdl.handle.net/11858/00-001M-0000-000E-B072-5
Abstract
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.