[论文解读] Two-stage fourth-order accurate time discretizations for 1D and 2D special relativistic hydrodynamics

Yuhuan Yuan, Huazhong Tang|arXiv (Cornell University)|Dec 15, 2017

Advanced Numerical Methods in Computational Mathematics参考文献 33被引用 7

一句话总结

本文提出了一种用于一维和二维特殊相对论流体动力学的两阶段四阶时间格式，利用直接欧拉型广义黎曼问题（GRP）方法，并通过解析求解局部准一维GRP状态。该格式在捕捉激波和复杂波结构方面实现了高阶精度与鲁棒性，通过标准黎曼问题和光滑流动测试的数值实验得到验证。

ABSTRACT

This paper studies the two-stage fourth-order accurate time discretization \cite{LI-DU:2016} and applies it to special relativistic hydrodynamical equations. It is shown that new two-stage fourth-order accurate time discretizations can be proposed. With the aid of the direct Eulerian GRP (generalized Riemann problem) methods \cite{Yang-He-Tang:2011,Yang-Tang:2012} and the analytical resolution of the local "quasi 1D" GRP, the two-stage fourth-order accurate time discretizations are successfully implemented for the 1D and 2D special relativistic hydrodynamical equations. Several numerical experiments demonstrate the performance and accuracy as well as robustness of our schemes.

研究动机与目标

• 开发适用于特殊相对论流体动力学（RHD）的更高阶精度数值格式，同时保持激波捕捉的鲁棒性与精度。
• 将两阶段四阶时间格式框架扩展至一维和二维RHD方程，克服传统龙格-库塔方法计算成本过高的问题。
• 采用直接欧拉型GRP方法，通过解析求解局部准一维GRP状态，实现更高效率与更高精度。
• 通过在黎曼问题和光滑流动上的全面数值实验，验证该格式的性能。

提出的方法

• 采用基于GRP通量函数推导的两阶段四阶时间格式，避免了对完整阶段黎曼求解器调用的需求。
• 采用直接欧拉型GRP方法，通过解析求解局部黎曼问题，计算通量的时间导数至二阶。
• 利用局部准一维GRP的解析解，计算保守变量的高阶时间导数，从而实现时间上的四阶精度。
• 将该格式应用于一维和二维特殊相对论性欧拉方程的守恒形式，采用$̳$-律状态方程。
• 通过特征分析和黎曼不变量，推导出原始变量时间导数的演化方程，尤其针对声学与激波情形。
• 以一致方式处理洛伦兹因子与相对论性通量，确保物理一致性与数值稳定性。

实验结果

研究问题

• RQ1 能否有效将两阶段四阶时间格式适配于一维和二维特殊相对论流体动力学方程？
• RQ2 与传统龙格-库塔方法相比，该方法在RHD中的精度与效率如何？
• RQ3 采用解析GRP求解的直接欧拉型GRP方法，能否在不显著增加计算成本的前提下实现高阶时间精度？
• RQ4 该格式在解析强激波与稀疏波等相对论流场中的鲁棒性如何？
• RQ5 该格式在标准一维和二维黎曼问题及光滑流动基准测试中的表现如何？

主要发现

• 收敛性研究表明，所提出的两阶段四阶时间格式在时间和空间上均实现了对一维和二维特殊相对论流体动力学的四阶精度。
• 该格式能以高分辨率准确捕捉复杂波结构，包括激波与稀疏波，且振荡极小。
• 数值实验表明，该方法在极端相对论区域（如高洛伦兹因子与强激波）下仍保持高度鲁棒性。
• 与标准四阶段龙格-库塔方法相比，计算成本显著降低，因为每个阶段均避免了完整的黎曼求解器调用。
• 对于光滑解，该格式在L1与L∞范数下均表现出最优收敛率，证实了理论精度阶数。
• 局部准一维GRP的解析求解能够高效准确地计算至二阶的时间导数，构成高阶精度的核心机制。