[论文解读] Directed Relativistic Blast Wave

Andrei Gruzinov|ArXiv.org|Apr 23, 2007

Laser-Plasma Interactions and Diagnostics被引用 8

一句话总结

本文提出了一种自相似的、定向的超相对论性激波解，该解作为动量几乎等于能量的通用定向爆炸的中间渐近吸引子。通过相对论流体力学和数值模拟，表明此类爆炸会演化为一种普遍的、狭窄的抛物形激波结构，其张角约为 ∼Γ⁻¹，其中 Γ 为洛伦兹因子，且该吸引子由一个收敛缓慢的能量积分表征，表明 β ≈ 2.032。

ABSTRACT

A spherically symmetrical ultra-relativistic blast wave is not an attractor of a generic asymmetric explosion. Spherical symmetry is reached only by the time the blast wave slows down to non-relativistic velocities, when the Sedov-Taylor-von Neumann attractor solution sets in. We show however, that a directed relativistic explosion, with the explosion momentum close to the explosion energy, produces a blast wave with a universal intermediate asymptotic -- a selfsimilar directed ultra-relativistic blast wave. This universality might be of interest for the astrophysics of gamma-ray burst afterglows.

研究动机与目标

• 识别动量 P ≈ 能量 E 的定向超相对论性爆炸的中间渐近行为。
• 确定此类定向爆炸是否导致一种与初始非对称性无关的普遍自相似激波结构。
• 研究激波角结构的动力学起源，特别是其张角与洛伦兹因子的标度关系。
• 通过数值模拟评估自相似解的有效性及其收敛性质。
• 阐明因果性与能量积分收敛性在决定吸引子物理可实现性中的作用。

提出的方法

• 在圆柱坐标系中推导出相对论流体力学的自相似假设，假设具有 Γ ≫ 1 的超相对论性流动。
• 引入无量纲变量 ξ = e⁻²ᵗ⁄ᵀ(t−z)/τ 和 η = e⁻ᵗ⁄ᵀr/τ，其中 τ = C(E/ρ)¹ᐟ³，以将偏微分方程简化为自相似函数 P(ξ,η)、Q(ξ,η)、U(ξ,η) 的常微分方程组。
• 应用激波间断条件和相对论流体力学方程（在 φ=0 处 Tᵘᵛ∂ᵥφ = τᵘᵛ∂ᵥφ），推导激波面上的边界条件。
• 对下游距离较大的情况执行渐近分析，提出含参数 β 的标度假设，并求解得到 β ≈ 2.032 的常微分方程组。
• 在 Lax 有限差分格式中采用以光速移动的坐标网格（800×800 分辨率）模拟轴对称的相对论性流体动力学。
• 将模拟结果（能量密度、激波位置、速度分布）与理论预测进行比较，以验证吸引子解。

实验结果

研究问题

• RQ1 动量 P ≈ E 的定向超相对论性爆炸是否演化为一种普遍的自相似激波结构？
• RQ2 在该类定向爆炸中，激波张角与洛伦兹因子 Γ 的标度关系如何？
• RQ3 为何 Blandford-McKee 的自相似解无法作为非球形、非定向爆炸的吸引子？
• RQ4 能量积分的收敛速率如何影响向吸引子状态的趋近过程？
• RQ5 数值模拟能否证实存在 β ≈ 2.032 的自相似解，且能量积分呈现对数发散？

主要发现

• 动量 P ≈ E 的定向超相对论性爆炸会演化为一种普遍的自相似激波结构，其激波张角 ∼Γ⁻¹，且与初始非对称性无关。
• 该自相似解为中间渐近吸引子，数值模拟结果表明不同初始构型下演化行为一致，得到验证。
• 当 β = 2 时，能量积分对数发散；数值解仅在 2 < β < 2.032 范围内存在，表明 β ≈ 2.032 为临界值。
• 数值模拟显示最大能量密度 ε_max 的衰减略快于指数衰减，与能量积分的缓慢收敛一致。
• 激波位置 x_s 几乎呈指数增长，且 ε_max ∝ x_s⁻²，证实了理论标度关系 (23)–(24)。
• 从 ε_max 衰减推断的时间常数 τ 与 τ ≈ 8x_s q_max 的理论预测在 10% 以内一致，验证了理论假设。