2007/03/31 by Peng Wang, Tom Abel, Weiqun Zhang · 1 citation
Engineering · Mathematics · Physics and Astronomy · #Adaptive mesh refinement #Aerodynamics and Acoustics in Jet Flows #Applied mathematics #Classical mechanics #Computational Fluid Dynamics and Aerodynamics #Computational science #Computer science #Fluid Dynamics and Turbulent Flows #Mathematics #Mechanics #Physics #Statistical physics #astro-ph
paper · pdf · doi:10.1086/529434
14 pages, 23 figures. A section on 3D GRB jet simulation added. Accepted by ApJS
arxiv created 2007/12/31 · openalex publication_date 2008/05/28 · arxiv updated 2009/12/01 · openalex created_date 2025/10/10 · openalex updated_date 2026/08/05
Astrophysical relativistic flow problems require high resolution three-dimensional numerical simulations. In this paper, we describe a new parallel three-dimensional code for simulations of special relativistic hydrodynamics (SRHD) using both spatially and temporally structured adaptive mesh refinement (AMR). We used the method of lines to discretize the SRHD equations spatially and a total variation diminishing (TVD) Runge-Kutta scheme for time integration. For spatial reconstruction, we have implemented piecewise linear method (PLM), piecewise parabolic method (PPM), third order convex essentially non-oscillatory (CENO) and third and fifth order weighted essentially non-oscillatory (WENO) schemes. Flux is computed using either direct flux reconstruction or approximate Riemann solvers including HLL, modified Marquina flux, local Lax-Friedrichs flux formulas and HLLC. The AMR part of the code is built on top of the cosmological Eulerian AMR code \sl enzo. We discuss the coupling of the AMR framework with the relativistic solvers. Via various test problems, we emphasize the importance of resolution studies in relativistic flow simulations because extremely high resolution is required especially when shear flows are present in the problem. We also present the results of two 3d simulations of astrophysical jets: AGN jets and GRB jets. Resolution study of those two cases further highlights the need of high resolutions to calculate accurately relativistic flow problems.