[论文解读] Very Light Jets I. Axisymmetric Parameter Study and Analytic Approximation

Martin Krause|ArXiv.org|Nov 20, 2002

Astrophysics and Cosmic Phenomena参考文献 28被引用 51

一句话总结

本研究对极轻喷流（密度对比度 η = 10⁻⁵ 至 10⁻²）进行了系统的轴对称流体动力学模拟，马赫数范围为 2.6 至 26，重点研究弓形激波和包层演化。研究推导了球形爆炸波阶段的解析解，表明喷流可将高达 10⁶⁰ erg 的能量传递至环境介质，在 r₁ 处，结论为紧凑对称源（CSOs）必须是年轻的，而非受致密介质约束，因为其头部传播速度极快。

ABSTRACT

The propagation of extragalactic jets is studied by a series of twelve axisymmetric hydrodynamic simulations. Motivated by observational constraints, but unlike most previous simulations, the regime of jet to external medium density (eta) from 10^-5 to 10^-2 is explored, for Mach numbers (M) between 2.6 and 26. The computational domain contained the bow shocks for the whole simulation time. The bow shocks are found to be spherical at source sizes below a critical value r1 (blastwave phase), which can reach up to 10 jet radii. [...] The numerical work is complemented by an analytic approach for the spherical phase. Extending previous work, the radial force balance could be integrated for arbitrary background density and energy input, which results in a global solution. The analytic results are shown to be consistent with the numerical work, and a lower limit to r1 can be calculated, which falls below the numerical results by a few jet radii. It is shown explicitely how a King density distribution changes the discussed aspects of the bow shock propagation. Because the jet head propagates very fast in the blastwave phase, it turns out that it is not possible to ``frustrate'' a jet by a high density environment. This is very important for the class of small radio galaxies (compact symmetric objects / GHz peaked sources): They have to be young. During its blastwave phase, a powerful jet can transfer typically 10^60 erg to the environmental gas. This is enough to balance the radiative losses in a cooling flow, if one of the cluster galaxies harbours a powerful jet every Gigayear.

研究动机与目标

• 研究极轻喷流（η < 10⁻³）在低喷流-环境密度对比度范围内的动力学行为，该领域在以往模拟中尚未被充分探索。
• 通过使用封闭的喷流注入边界和足够大的计算域，实现对整个模拟时间内弓形激波和包层演化的解析。
• 发展并验证喷流驱动弓形激波球形爆炸波阶段的解析模型，扩展先前工作，适用于任意背景密度和能量输入。
• 将来自天鹅座 A 的观测约束（长宽比 ≈1.2，包层-弓形激波宽度比 ≈1）与模拟喷流参数（特别是 η 和 M）进行对比。
• 确定致密环境是否可能阻碍喷流传播，尤其是对紧凑对称源（CSOs）和吉赫兹峰值源（GPS）的影响。

提出的方法

• 使用有限体积欧拉格式代码，开展 12 次轴对称流体动力学模拟，喷流注入侧采用封闭边界，计算域足够大以容纳整个模拟期间的弓形激波。
• 在 η 从 10⁻⁵ 至 10⁻² 和内部马赫数 M 从 2.6 至 26 的范围内变化，探索极轻喷流区域。
• 在球形阶段应用径向力平衡，对任意背景密度和能量输入积分运动方程，获得全局解析解。
• 采用金氏密度分布描述外部介质，并评估其对弓形激波和包层形态的影响。
• 将数值结果与解析解对比，验证模型，并推导出从球形到拉长弓形激波过渡半径 r₁ 的下限。
• 利用推导出的解析与数值结果，计算弓形激波在 r₁ 处传递给周围气体的能量。

实验结果

研究问题

• RQ1 在极轻喷流（η < 10⁻³）中，弓形激波的长宽比如何随时间演化？从球形向拉长形态的转变由什么决定？
• RQ2 致密外部介质是否可能阻碍喷流传播，还是在爆炸波阶段喷流头部始终快速前进？
• RQ3 当喷流到达半径 r₁ 时，其向周围介质传递的总能量是多少？该能量是否足以平衡星系团中的辐射冷却？
• RQ4 模拟得到的包层和弓形激波尺寸与钱德拉 X 射线望远镜对天鹅座 A 的观测相比如何，特别是长宽比和包层-弓形激波宽度比？
• RQ5 外部密度分布的选择（如金氏分布）在多大程度上影响球形阶段弓形激波和包层的动力学行为？

主要发现

• 在早期阶段（r < r₁），弓形激波保持球形，仅当激波变为亚音速时才发生向拉长形态的转变，此时 r₁ 可达喷流半径的 10 倍。
• 基于径向力平衡并考虑任意背景密度和能量输入的球形阶段解析解，与数值模拟一致，并为 r₁ 提供了一个下限，该下限比数值结果低约几个喷流半径。
• 强大喷流可在 r₁ 处将高达 10⁶⁰ erg 的能量传递给周围气体，计算公式为 E_gas = 10⁶⁰ erg × (L/10⁴⁶ erg/s)^1.5 × (v_j/c/3)^-2.5 × (η/10⁻⁴)^-1.25 × (n₀/0.1 cm⁻³)^-0.5。
• 该能量传递足以平衡星系团中的辐射冷却，意味着每 10⁹ 年一次此类喷流事件即可维持星系团中观测到的 keV 量级温度。
• 包层-外部介质界面处的开尔文-赫姆霍兹不稳定性普遍存在，可能负责卷入并加速受激波加热的气体，从而可能驱动 X 射线或发射线特征。
• 紧凑对称源（CSOs）和吉赫兹峰值源必须是年轻的，而非受致密环境约束，因为其喷流头部传播速度过快（例如，在致密星际介质中穿越 1 kpc 仅需约 0.3 Myr），即使在极端情况下，其年龄也难以超过约 0.3 Myr。