PIERNIK is a freely available, massively parallel, multi-fluid, multi-grid, magnetohydrodynamical (MHD) code. We have implemented the multi-state Harten-Lax-van Leer approximate Riemann solver for ideal magnetohydrodynamics - HLLD (”D” denotes Discontinuities) in PIERNIK. Our goal is to perform numerical simulations of hypersonic astrophysical flows in 1, 2 and 3 spatial dimensions, particularly in the presence of discontinuities, using the HLLD Riemann solver. The HLLD Riemann solver is more robust than the linearized Riemann solver, capable of resolving isolated discontinuities formed in the MHD system (hence named as HLLD) and corresponds to the HLLC (”C” denotes Contact) Riemann solver when the magnetic field vanishes. We demonstrate the robustness and accuracy of the HLLD Riemann solver implemented in PIERNIK through numerical tests and several interesting features available in the PIERNIK

Keywords: Computational methods, Magnetohydrodynamics, Turbulence