[论文解读] A Simulation Study of Ultra-Relativistic Jets -- I. A New Code for Relativistic Hydrodynamics

Jeongbhin Seo, Hyesung Kang+3 more|arXiv (Cornell University)|Jun 8, 2021

Computational Fluid Dynamics and Aerodynamics参考文献 59被引用 14

一句话总结

本文提出了一种基于五阶有限差分WENO格式的新高阶精度相对论性流体动力学（RHD）代码，具备多种先进功能，包括SSPRK时间积分、高阶横向通量平均化，以及通过特征值修正抑制卡布恩克勒不稳定性。该代码专为在星系团介质（ICM）中模拟超相对论喷流而优化，WENO-Z、SSPRK和RC状态方程被选定为基准配置，以实现对激波、湍流及喷流动力学的稳健且精确模拟。

ABSTRACT

In an attempt to investigate the structures of ultra-relativistic jets injected into the intracluster medium (ICM) and the associated flow dynamics, such as shocks, velocity shear, and turbulence, we have developed a new special relativistic hydrodynamic (RHD) code in the Cartesian coordinates, based on the weighted essentially non-oscillatory (WENO) scheme. It is a finite difference scheme of high spatial accuracy, which has been widely employed for solving hyperbolic systems of conservation equations. The code is equipped with different WENO versions, such as the 5th-order accurate WENO-JS (Jiang & Shu 1996), WENO-Z, and WENO-ZA, and different time integration methods, such as the 4th-order accurate Runge-Kutta (RK4) and strong stability preserving RK (SSPRK), as well as the implementation of the equations of state (EOSs) that closely approximate the EOS of the single-component perfect gas in relativistic regime. In addition, it incorporates a high-order accurate averaging of fluxes along the transverse directions to enhance the accuracy of multi-dimensional problems, and a modification of eigenvalues for the acoustic modes to effectively control the carbuncle instability. Through extensive numerical tests, we assess the accuracy and robustness of the code, and choose WENO-Z, SSPRK, and the EOS suggested in Ryu et al. (2006) as the fiducial setup for simulations of ultra-relativistic jets. The results of our study of ultra-relativistic jets using the code is reported in an accompanying paper (Seo et al. 2021, Paper II).

研究动机与目标

• 开发一种高阶、稳健的RHD代码，用于模拟星系团介质（ICM）中的超相对论喷流。
• 解决多维流动中卡布恩克勒不稳定性与横向耗散等数值挑战。
• 识别数值格式与状态方程的最佳组合，以精确模拟喷流结构与动力学。

提出的方法

• 实现五阶有限差分WENO格式：WENO-JS、WENO-Z和WENO-ZA，用于高分辨率重构。
• 采用四阶精度的SSPRK与RK4时间积分方法，确保在激波主导流动中的稳定性与精度。
• 引入高阶横向通量平均化方法，提升多维问题的精度。
• 使用修正后的特征值处理声学模态，以抑制卡布恩克勒不稳定性，同时不损失分辨率。
• 采用RC状态方程（EOS）准确模拟喷流中γ ≈ 1.45至激波区域γ ≈ 4/3与5/3的绝热指数转变。
• 在1D、2D和3D问题中进行广泛测试，以验证精度、稳定性和性能。

实验结果

研究问题

• RQ1 哪种WENO重构格式在超相对论喷流模拟中能实现精度与鲁棒性的最佳平衡？
• RQ2 不同时间积分格式（SSPRK与RK4）如何影响具有大横向速度的激波模拟？
• RQ3 高阶横向通量平均化能否提升多维相对论流体动力学的精度？
• RQ4 特征值修正在相对论喷流弓形激波中在多大程度上能抑制卡布恩克勒不稳定性？
• RQ5 哪种状态方程（ID、TM、RC）最能准确捕捉喷流与ICM流体在激波加热与混合过程中的热力学演化？

主要发现

• WENO-Z格式在分辨率与鲁棒性之间实现了最佳平衡，优于WENO-JS（过于耗散）和WENO-ZA（在极端条件下过于敏感）。
• SSPRK时间积分在具有大横向速度的激波主导流动中表现优于RK4，确保了更高的稳定性并减少了振荡。
• 高阶横向通量平均化显著提升了多维模拟的精度，尤其在喷流包层与回流界面等复杂流动区域表现突出。
• 通过φ ≈ 5–10的特征值修正可有效抑制弓形激波表面的卡布恩克勒不稳定性，其中高功率喷流（Qjet ≳ 10⁴⁶ erg s⁻¹）采用φ ≈ 10，低功率喷流则采用φ ≈ 5。
• RC状态方程的模拟结果与TM相当，且比ID更准确逼近相对论性理想气体EOS，因此被选为准确模拟绝热指数从γ ≈ 1.45过渡至γ ≈ 5/3的默认配置。
• 采用基准配置（WENO-Z、SSPRK、RC EOS、特征值修正）的模拟成功捕捉了复杂的喷流形态，包括剪切不稳定性（KH）、弓形激波及湍流包层结构。