[论文解读] Conservative Moment Equations for Neutrino Radiation Transport with Limited Relativity

Eirik Endeve, Christian Y. Cardall+1 more|arXiv (Cornell University)|Dec 17, 2012

Gamma-ray bursts and supernovae参考文献 12被引用 6

一句话总结

本文在有限相对论近似下，推导了核心坍缩超新星中中微子输运的保守、多维、能量依赖的矩方程，重点在于四动量和轻子数的一致守恒。文章提出了两种自洽的$$\mathcal{O}(v)$$形式——'加'与'减'——确保能量、动量和轻子数在流体速度的指定阶次下守恒，从而实现更精确的中微子驱动超新星爆炸模拟。

ABSTRACT

We derive conservative, multidimensional, energy-dependent moment equations for neutrino transport in core-collapse supernovae and related astrophysical systems, with particular attention to the consistency of conservative four-momentum and lepton number transport equations. After taking angular moments of conservative formulations of the general relativistic Boltzmann equation, we specialize to a conformally flat spacetime, which also serves as the basis for four further limits. Two of these---the multidimensional special relativistic case, and a conformally flat formulation of the spherically symmetric general relativistic case---are given in appendices for the sake of comparison with extant literature. The third limit is a weak-field, `pseudo-Newtonian' approach \citep{kim_etal_2009,kim_etal_2012} in which the source of the gravitational potential includes the trace of the stress-energy tensor (rather than just the mass density), and all orders in fluid velocity $v$ are retained. Our primary interest here is in the fourth limit: `$\mathcal{O}(v)$' moment equations for use in conjunction with Newtonian self-gravitating hydrodynamics. We show that the concept of `$\mathcal{O}(v)$' transport requires care when dealing with both conservative four-momentum and conservative lepton number transport, and present two self-consistent options: `$\mathcal{O}(v)$-plus' transport, in which an $\mathcal{O}(v^2)$ energy equation combines with an $\mathcal{O}(v)$ momentum equation to give an $\mathcal{O}(v^2)$ number equation; and `$\mathcal{O}(v)$-minus' transport, in which an $\mathcal{O}(v)$ energy equation combines with an $\mathcal{O}(1)$ momentum equation to give an $\mathcal{O}(v)$ number equation.

研究动机与目标

• 在有限相对论近似下，为核心坍缩超新星中的中微子输运发展保守、多维、能量依赖的矩方程。
• 在数值模拟中确保四动量和轻子数的一致守恒，这对精确建模中微子驱动爆炸至关重要。
• 解决现有形式中非保守能量和动量方程无法保持总能量和轻子数的问题。
• 提出两种自洽的$$\mathcal{O}(v)$$输运方案——'加'与'减'——确保能量、动量和数量方程均保持守恒。
• 为牛顿流体动力学中包含至$$\mathcal{O}(v)$$阶相对论修正的中微子输运提供稳健、高保真的模拟基础。

提出的方法

• 通过取共形平坦时空下保守的一般相对论玻尔兹曼方程的角矩，推导矩方程。
• 将一般形式特化为四种极限：多维狭义相对论、球对称广义相对论、弱场'伪牛顿'方法，以及主要关注的：牛顿流体动力学下的$$\mathcal{O}(v)$$输运。
• 提出两种自洽形式：'加'（$$\mathcal{O}(v^2)$$能量 + $$\mathcal{O}(v)$$动量 → $$\mathcal{O}(v^2)$$数量）与'减'（$$\mathcal{O}(v)$$能量 + $$\mathcal{O}(1)$$动量 → $$\mathcal{O}(v)$$数量）。
• 通过证明时间、空间和几何导数项的源项在方程间相互抵消，确保能量、动量和数量通量之间的一致性，从而实现守恒。
• 利用相同框架推导中微子数密度的保守演化方程，确保在离散数值方案中轻子数守恒。
• 通过与现有形式（如MJD10）比较验证该方法，表明非保守形式因速度依赖的源项而无法保持能量和轻子数守恒。

实验结果

研究问题

• RQ1 如何在多维、能量依赖的模拟中推导出中微子输运的保守矩方程，同时确保四动量和轻子数的一致守恒？
• RQ2 在中微子输运中使用非保守形式的能量和动量方程会产生何种影响，特别是对能量和轻子数守恒的影响？
• RQ3 如何将$$\mathcal{O}(v)$$阶相对论修正一致地纳入牛顿流体动力学中的中微子输运，而不违反守恒定律？
• RQ4 两种自洽的$$\mathcal{O}(v)$$形式——'加'与'减'——分别是什么？它们在能量、动量和数量方程中的精度阶次有何不同？
• RQ5 在离散形式中需满足何种条件，才能确保数值模拟中能量和轻子数守恒？

主要发现

• 本文在共形平坦时空中推导出中微子输运的保守矩方程，特别聚焦于适用于牛顿流体动力学的$$\mathcal{O}(v)$$极限。
• 提出两种自洽形式：'加'输运（能量方程为$$\mathcal{O}(v^2)$$，数量方程为$$\mathcal{O}(v^2)$$）与'减'输运（能量和数量方程均为$$\mathcal{O}(v)$$）。
• 保守能量和动量方程（方程B2和B3）通过时间、空间和几何源项的导数项相互抵消，与保守数量方程（方程B5）保持一致。
• 非保守形式（如MJD10中所用）因速度依赖的源项在无引力或相互作用时无法消去，导致能量和轻子数不守恒。
• 能量与数量方程之间的关系简单明了：$$\mathcal{J} = \epsilon \mathcal{J}_{\scriptscriptstyle{\mathcal{N}}}$$，当同时需要能量和数量守恒时，该关系可简化数值实现。
• 结果为开发在多维中微子输运模拟中同时守恒能量和轻子数的数值方案提供了严格基础，这对精确建模核心坍缩超新星至关重要。