[论文解读] Three-Dimensional Modeling of the Deflagration Stage of a Type Ia Supernova Explosion

A. M. Khokhlov|arXiv (Cornell University)|Aug 29, 2000

Gamma-ray bursts and supernovae被引用 10

一句话总结

本文首次对钱德拉塞卡质量的碳氧白矮星的爆燃阶段进行了三维、高分辨率的数值模拟，采用火焰捕捉方法并结合湍流火焰速度模型。结果表明，湍流火焰速度在膨胀导致大尺度湍流冻结之前增长至声速的约30%，随后下降，从而限制能量释放，导致未燃烧的CO层和中间质量元素的形成；这些挑战可通过延迟引爆模型解决。

ABSTRACT

The paper describes a physical model and numerical algorithm for modeling Type Ia supernova (SNIa) explosions in three dimensions and presents first results of modeling a deflagration explosion in a nonrotating, Chandrasekhar-mass carbon-oxygen (CO) white dwarf. Simulations show that the turbulent flame speed grows exponentially, reaches approximately 30% of the speed of sound, and then declines as the large-scale turbulence is frozen by expansion. The freezing of turbulent motions appears to be a crucial physical mechanism regulating the rate of deflagration in SNIa. The energy of the explosion is comparable to that of a typical SNIa. However, the presence of the outer layer of unburned CO and the formation of intermediate mass elements and pockets of unburned CO near the center pose problems for the modeling of SNIa spectra. Delayed detonation is a way to alleviate these problems and to produce consistent spectra.

研究动机与目标

• 开发并应用三维数值模型，模拟Ia型超新星的爆燃阶段，解决二维模拟的局限性。
• 研究三维爆燃模拟中湍流火焰速度的演化，特别是湍流能量级联与膨胀驱动的大尺度运动冻结之间的竞争。
• 评估纯爆燃模型是否能产生能量和成分与观测一致的爆炸，或是否需要延迟引爆。
• 探究爆炸结果对着火位置和初始碳含量的敏感性，识别SNIa多样性的潜在‘隐藏参数’。

提出的方法

• 使用三维欧拉型流体动力学求解器模拟无粘性、可压缩流体，包含重力和核能沉积。
• 通过反应-扩散方程捕捉火焰，引入标量火焰变量f，并设定预设的湍流火焰速度S_t。
• 对简并物质的方程状态进行建模，包含电子、正电子、离子和辐射，T和P由ρ、E_i、Y_e和Y_i计算得出。
• 采用亚网格尺度模型描述湍流火焰速度，考虑三维中的瑞利-泰勒不稳定性及能量级联。
• 利用球对称密度分布的角平均重力场，模拟引力场。
• 从点火开始模拟至约2秒，追踪能量释放、膨胀及成分演化。

实验结果

研究问题

• RQ1 在钱德拉塞卡质量的CO白矮星三维爆燃模拟中，湍流火焰速度如何演化？膨胀在其中起何种调控作用？
• RQ2 纯爆燃模型能否产生能量和成分与观测到的Ia型超新星一致的爆炸？
• RQ3 三维湍流和大尺度浮力上升流对未燃烧CO层及中间质量元素形成有何影响？
• RQ4 在此三维设置中，爆燃向引爆的转变（DDT）在物理上是否可行？何种条件可能触发它？
• RQ5 爆炸结果对初始条件（如点火位置和初始碳丰度）的敏感性如何？

主要发现

• 湍流火焰速度（S_t）呈指数增长，达到声速约30%后，由于膨胀引起的大尺度湍流运动冻结而开始下降。
• 膨胀导致大尺度湍流冻结被识别为限制SNIa中能量释放和爆燃速率的关键物理机制。
• 模拟产生的爆炸能量与典型Ia型超新星相当，但外层存在显著的未燃烧CO层，内区则含有中间质量元素。
• 内区存在未燃烧CO和成分不均一性，对重现观测到的SNIa光谱构成挑战，表明纯爆燃模型存在局限性。
• 大尺度浮力上升流由瑞利-泰勒不稳定性驱动，但膨胀最终克服浮力，导致外层快速膨胀。
• 延迟引爆被证明是克服上述问题的可行方案，因其能高效烧尽剩余未燃烧物质并实现成分均质化。