[论文解读] Three-Dimensional Simulations of Classical Novae

A. Kercek, W. Hillebrandt+1 more|arXiv (Cornell University)|Nov 17, 1998

Astrophysical Phenomena and Observations参考文献 1被引用 23

一句话总结

本研究首次对经典新星爆发进行了3D流体动力学模拟，采用基于笛卡尔网格的3D欧拉型代码与核反应网络，模拟了1 M⊙ C+O白矮星上湍流氢燃烧过程。关键发现为：在爆发期间，吸积包层自增丰碳和氧的过程过于缓慢，即使在五倍太阳金属度下也无法产生快速新星，这意味着爆发前的增丰机制对快速新星形成至关重要。

ABSTRACT

We present first results of three-dimensional (3D-) calculations of turbulent and degenerate hydrogen-burning on top of a C+O white dwarf of one solar mass. The simulations are carried out by means of a code which solves Euler's equation for an arbitrary equation of state together with a nuclear reaction network and the energy input from nuclear reactions on a Cartesian grid covering a fraction of the white dwarf's surface and accreted atmosphere. The flow patterns we obtain are very different from those of earlier 2D simulations using the same initial conditions and the same numerical resolution. The possibility of self-enrichment of the accreted hydrogen-rich atmosphere with carbon and oxygen from the surface layers of the white dwarf during the violent phase of the burning is investigated, and it is demonstrated that self-enrichment proceeds too slowly if the accreted gas has near-solar CNO-abundances at the onset of the thermonuclear runaway. As a result, we do not find a fast nova outburst. This conclusion remains valid if the initial metallicity of the accreted gas is raised by a factor of five. Therefore we conclude that fast nova outbursts indeed require huge enrichments of C and O, as postulated from spherically symmetric models, and that the mechanism which leads to such enhancements must operate prior to the outburst.

研究动机与目标

• 研究3D流体动力学在经典新星爆发中的作用，特别是热核 runaway 期间湍流对流与混合的动力学过程。
• 检验白矮星表面的碳和氧是否能通过吸积的富氢包层自增丰，从而驱动快速新星爆发。
• 将3D模拟结果与以往2D模拟进行比较，评估维度变化对混合、能量生成及爆发动力学的影响。
• 评估初始金属度（最高至5倍太阳金属度）是否能显著加速燃烧过程，从而产生快速新星特征。
• 确定3D模拟中数值扩散或未解析的小尺度运动是否比2D模拟更强烈地抑制混合，从而影响爆发的可能。

提出的方法

• 采用PROMETHEUS代码进行数值模拟，该代码为PPM型流体动力学代码，基于3D笛卡尔网格对欧拉方程进行隐式积分，覆盖白矮星表面和大气层的一部分区域。
• 求解与12种核素反应网络（1H, 4He, 12C, 13C, 13N, 14N, 15N, 14O, 15O, 16O, 17O, 17F）耦合的流体动力学方程，反应速率采用Wallace & Woosley（1981）及Thielemann（私人通信）的数据。
• 同时求解核反应产生的能量生成与网络方程，以避免数值不稳定性，采用一致的状态方程和隐式时间积分方法。
• 初始条件基于Glasner等（1997）的1D隐式流体动力学模型，模拟吸积与爆发前演化过程，白矮星质量为1.0 M⊙，吸积包层质量为2×10⁻⁵ M⊙，金属度Z=0.02（部分模拟中为5倍太阳金属度）。
• 3D网格初始配置基于1D模型的静力平衡结构，经过多个动力时序的松弛后，施加扰动以触发对流与燃烧。
• 采用周期性边界条件以避免非物理反射，模拟在中等分辨率下进行（例如220×100网格点），以确保积分量的稳定性和收敛性。

实验结果

研究问题

• RQ1 经典新星爆发期间的3D流体动力学流动与2D模拟相比有何不同，特别是在涡度结构与混合效率方面？
• RQ2 在热核 runaway 期间，吸积的富氢包层能否足够快地从白矮星表面获得碳和氧的自增丰，从而产生快速新星？
• RQ3 提高吸积气体的初始金属度（最高至5倍太阳金属度）是否能显著加速燃烧过程并导致快速新星爆发？
• RQ4 3D模拟在多大程度上抑制了对流越界与C、O的上翻过程，与2D模型相比，这如何影响能量释放与爆发动力学？
• RQ5 3D模拟中的数值扩散是否足以驱动混合，还是必须采用显式亚网格尺度湍流建模才能准确捕捉新星爆发中的混合过程？

主要发现

• 3D模拟产生的流动模式与2D模拟显著不同，表现出更强的小尺度湍流运动，从而减少了对流越界并抑制了C和O向燃烧层的混合。
• 即使初始金属度提高至5倍太阳金属度，爆发期间吸积大气层的C和O自增丰过程仍过于缓慢，无法触发快速新星，因为温度上升 timescale 仍过长（约100秒量级）。
• 在高金属度模拟中，150秒后峰值速度达到约2×10⁷ cm/s，约为低金属度情况的两倍，但能量生成率与最终大气金属度仅 moderately 提高。
• 在400秒时，高金属度情况下的系统与低金属度情况下的系统在燃烧状态上已无法区分，表明提高初始金属度无法克服根本的时间尺度问题。
• 模拟结果表明，由于混合减少，3D中爆炸性氢燃烧被削弱，导致爆发比2D中更不剧烈，因此通过自增丰产生快速新星的可能性更低。
• 结果排除了爆发期间自增丰作为快速新星形成的可行机制，意味着大质量C和O增丰必须在爆发前发生，可能通过剪切不稳定性或吸积阶段的长期混合实现。