Guangqi: A two-dimensional radiation hydrodynamic code with realistic equation of states

TL;DR

Guangqi is a new two-dimensional radiation hydrodynamics code that accurately models complex equations of state and radiation transport for astrophysical simulations, with high scalability and verified performance.

Contribution

The paper introduces Guangqi, a novel 2D radiation hydrodynamics code with advanced features like realistic EoS support, angular momentum conservation in spherical coordinates, and scalable parallelization.

Findings

Successfully verified against standard benchmarks.

Supports complex equations of state for hydrogen and helium.

Employs innovative algorithms for angular momentum conservation.

Abstract

We present Guangqi, a new two-dimensional, finite-volume radiation hydrodynamics code designed for high-performance astrophysical simulations. The code simultaneously resolves the hydrodynamic equations for complex equations of state (EoS) and implicit radiation transport under the flux-limited diffusion approximation. Written in Fortran and parallelized via the Message Passing Interface. Guangqi supports analytic hydrogen and helium EoS under the assumption of local thermal and chemical equilibrium. The framework is compatible with both Cartesian and spherical-polar geometries -- utilizing non-uniform grid spacing -- and incorporates static (SMR) and adaptive mesh refinement to optimize computational efficiency. To address the inherent challenges of angular momentum conservation in spherical-polar coordinates, we implement a robust and consistent "passive scalar angular momentum algorithm" (PSAMA). Domain decomposition is managed through both Z-order and Hilbert space-filling curves to ensure scalability. The code has been rigorously verified against a suite of standard benchmarks and newly designed test cases specifically intended to diagnose the non-linear coupling between gas dynamics, intricate EoS, radiation transport, and angular momentum conservation.