Relativistic Hydrodynamics with Wavelets

Methods to solve the relativistic hydrodynamic equations are important in a large number of astrophysical simulations, which may be very dynamic and involve multiscale features. This requires computational methods that are highly adaptive and capable of automatically resolving numerous localized features and instabilities that emerge across the computational domain and over many temporal scales. While this has been historically accomplished with adaptive-mesh-refinement-based methods, alternatives using wavelet bases and the wavelet transformation have recently achieved significant success in adaptive representation for advanced engineering applications. The current work presents a new method, extending the wavelet adaptive multiresolution representation method, for the integration of the relativistic hydrodynamic equations using iterated interpolating wavelets and introduces a highly adaptive implementation for multidimensional simulation. The wavelet coefficients provide a direct measure of the local approximation error for the solution and place collocation points that naturally adapt to the fluid flow while providing good conservation of fluid quantities. The resulting implementation, oahu, is applied to a series of demanding 1D and 2D problems that explore high Lorentz factor outflows and the formation of several instabilities, including the Kelvin-Helmholtz instability and the Rayleigh-Taylor instability.