Hydrodynamic simulations of AGN jets the impact of Riemann solvers and spatial reconstruction schemes on jet evolution

Numerical simulations play an essential role in helping us to understand the physical processes behind relativistic jets in active galactic nuclei. The large number of hydrodynamic codes available today enables a variety of different numerical algorithms to be utilized when conducting the simulations. Since many of the simulations presented in the literature use different combinations of algorithms it is important to quantify the differences in jet evolution that can arise due to the precise numerical schemes used. We conduct a series of simulations using the FLASH (magneto-)hydrodynamics code in which we vary the Riemann solver and spatial reconstruction schemes to determine their impact on the evolution and dynamics of the jets. For highly refined grids the variation in the simulation results introduced by the different combinations of spatial reconstruction scheme and Riemann solver is typically small. A high level of convergence is found for simulations using third-order spatial reconstruction with the Harten–Lax–Van-Leer with contact and Hybrid Riemann solvers.