[论文解读] Time dependent cosmic-ray shock acceleration with self-consistent injection

U. D. J. Gieseler, T. W. Jones+1 more|arXiv (Cornell University)|Nov 2, 2000

Solar and Space Plasma Dynamics被引用 13

一句话总结

本文提出了一种自洽的时变宇宙射线激波加速数值模型，基于流体物理的注入机制，利用波粒相互作用调节注入过程。结果表明，注入率稳定在热粒子的约10⁻³，强激波在达到稳态前的加速效率可达约20%，展示了宇宙射线激波特性的有效自调节。

ABSTRACT

One of the key questions to understanding the efficiency of diffusive shock acceleration of the cosmic rays (CRs) is the injection process from thermal particles. A self-consistent injection model based on the interactions of the suprathermal particles with self-generated magneto-hydrodynamic waves has been developed recently by Malkov (1998, Phys. Rev. E 58, 4911). By adopting this analytic solution, a numerical treatment of the plasma-physical injection model at a strong quasi-parallel shock has been devised and incorporated into the combined gas dynamics and the CR diffusion-convection code. In order to investigate self-consistently the injection and acceleration efficiencies, we have applied this code to the CR modified shocks of both high and low Mach numbers (M=30 and M=2.24) with a Bohm type diffusion model. We find the injection process is self-regulated in such a way that the injection rate reaches and stays at a nearly stable value after quick initial adjustment. For both shocks about 0.1% of the incoming thermal particles are injected into the CRs. For the weak shock, the shock has reached a steady state within our integration time and about 10% of the total available shock energy is transfered into the CR energy density. The strong shock has achieved a higher acceleration efficiency of about 20% by the end of our simulation, but has not yet reached a steady-state.

研究动机与目标

• 基于流体物理而非人为假设，发展一种自洽的激波宇宙射线注入模型。
• 研究波粒相互作用如何调节热粒子注入宇宙射线群体的过程。
• 模拟具有真实注入和扩散（Bohm型）的时变激波特性和加速过程，重点关注强激波和弱激波。
• 评估能量从激波向宇宙射线的转移效率以及对激波特性的反馈作用。
• 探索长期演化至宇宙射线主导状态的过程，以及当前计算方法的局限性。

提出的方法

• 采用Malkov的解析解，描述通过自激发阿尔芬波实现的超热粒子注入。
• 将注入模型整合进结合气体动力学与宇宙射线扩散-对流的代码中。
• 采用Bohm型扩散模型，空间网格分辨率足够精细，以捕捉近热动量扩散（与p²成正比）。
• 对准平行激波进行时变模拟，马赫数分别为M = 30（强激波）和M = 2.24（弱激波）。
• 跟踪相空间分布函数g(x,p)随时间的变化，直至最大动量p_max/m_p c ≈ 1。
• 计划通过自适应网格细化（AMR）扩展模拟，以研究更高宇宙射线动量和更长 timescales。

实验结果

研究问题

• RQ1 自洽的波粒相互作用注入如何影响激波处宇宙射线加速的时变演化？
• RQ2 其注入效率如何？是否随时间趋于稳定？
• RQ3 从激波向宇宙射线的能量转移效率如何随激波特性的强弱（马赫数）而变化？
• RQ4 宇宙射线反馈和波的生成在多大程度上调节注入过程？
• RQ5 与纯运动学注入模型相比，结果在倾斜激波中表现如何？

主要发现

• 注入过程具有自调节特性，在初始瞬态阶段后达到并维持近似稳定的速率。
• 对于强激波和弱激波，均有约10⁻³的入射热粒子被注入宇宙射线群体。
• 弱激波（M = 2.24）在模拟时间尺度内达到稳态，约10%的激波可用能量被转移至宇宙射线。
• 强激波（M = 30）实现了约20%的加速效率，但尚未达到稳态，表明宇宙射线压力仍在持续增长。
• 在粒子处于弱相对论性阶段的早期，宇宙射线压力保持亚主导，无法实现激波的完全修改。
• 该模型表明，未来使用自适应网格细化（AMR）的模拟对于解析向宇宙射线主导状态的长期演化至关重要。