[论文解读] Simulating bulk viscosity in neutron stars. I. Formalism

Giovanni Camelio, Lorenzo Gavassino+3 more|arXiv (Cornell University)|Apr 25, 2022

Pulsars and Gravitational Waves Research参考文献 63被引用 3

一句话总结

本文開發並比較了三種理論框架——多組分反應流體、Hiscock-Lindblom（HIS）二階Müller-Israel-Stewart理論，以及Maxwell-Cattaneo線性極限——用於在廣義相對論流體力學下模擬中子星中的體積黏度。推導了球對稱下的完整動力學方程組，建立了體積黏性流體與反應混合物之間的數學對偶性，並計算了Hiscock-Lindblom非線性區的等效聲速，為包含真實耗散效應的中子星動力學數值模擬奠定了基礎。

ABSTRACT

The faithful inclusion of the effects of bulk viscosity induced by the presence of chemical reactions is an important issue for simulations of core-collapse supernovae, binary neutron star mergers, and neutron star oscillations, where particle abundances are locally pushed out of chemical equilibrium by rarefaction and compression of the fluid elements. In this work, we discuss three different approaches that can be used to implement bulk viscosity in general relativistic hydrodynamic simulations of neutron stars: the exact multi-component reacting fluid, and two M\"uller-Israel-Stewart theories, namely the second order Hiscock-Lindblom model and its linear limit, the Maxwell-Cattaneo model. After discussing the theory behind the three approaches, we specialize their dynamics equations to spherical symmetry in the radial gauge-polar slicing (i.e., Schwarzschild) coordinates. We also discuss a particular choice for the equation of state of the fluid and the associated neutrino emission rates, which are used in a companion paper for the numerical comparison of the three frameworks, and we obtain the effective sound speed for the Hiscock-Lindblom theory in the non-linear regime.

研究动机与目标

• 開發一個一致的理論框架，用於在廣義相對論中子星模擬中納入化學反應引起的體積黏度。
• 比較三種不同的方法——多組分反應流體、Hiscock-Lindblom（HIS）理論與Maxwell-Cattaneo模型——用於模擬中子星中的體積黏度。
• 將這些模型的流體動力學方程專用于球對稱，使用史瓦西座標（徑向規範-極坐標切片）以實現一維數值實現。
• 建立體積黏性流體與反應混合物之間的數學對偶性，驗證多組分流體作為Müller-Israel-Stewart流體力學的數值可行替代方案。
• 計算Hiscock-Lindblom理論非線性區的等效聲速，這對數值穩定性與物理準確性至關重要。

提出的方法

• 推導多組分反應流體的完整相對論流體動力學方程組，包含反應過程中的中微子發射率與能量損失。
• 應用Müller-Israel-Stewart形式化方法來模擬體積黏度，特別是實現二階Hiscock-Lindblom理論及其線性化Maxwell-Cattaneo極限。
• 使用史瓦西座標系（徑向規範-極坐標切片）將所有三種模型專用于球對稱，從而實現一維數值實現。
• 引入特定的狀態方程與粒子反應速率（基於Urca過程），以供後續數值工作使用。
• 利用特徵線法計算Hiscock-Lindblom模型在非線性區的特徵信號傳播速度（等效聲速）。
• 推導熱力學一致性條件，包括矩陣 ∂²u/∂Ya∂Yb 的非負性，以確保反應流體框架中的熱力學穩定性與物理有效性。

实验结果

研究问题

• RQ1 多組分反應流體框架能否作為Müller-Israel-Stewart流體力學在模擬中子星體積黏度時的數值穩健替代方案？
• RQ2 Hiscock-Lindblom與Maxwell-Cattaneo模型在非線性區的等效聲速如何比較？它們與平衡態聲速有何關係？
• RQ3 在Müller-Israel-Stewart框架中，反應導致的能量損失的正確表達式是什麼？它如何影響流體動力學方程？
• RQ4 在球對稱條件下，反應混合物與體積黏性流體之間的數學對偶性如何體現在流體動力學方程中？
• RQ5 Hiscock-Lindblom理論非線性區的特徵信號傳播速度是什麼？它如何確保雙曲性與數值穩定性？

主要发现

• 非線性Hiscock-Lindblom理論中的等效聲速推導為 c²_s,uv = a1 − b1/τ + ˜ζ/(τh) + (n/h)(a2 − b2/τ)，在線性極限下簡化為Maxwell-Cattaneo結果。
• 當 χ⁻¹ = ρ³uab(∂Y_eq_a/∂ρ)ₛ(∂Y_eq_b/∂ρ)ₛ 時，多組分反應流體的平衡紫外聲速與體積黏性理論一致，確認了數學對偶性。
• 多組分流體模型已擴展以包含中微子發射率，從而實現中子星模擬中真實的能量損失建模。
• 所有三種模型的流體動力學方程均已完整推導，並專用于史瓦西座標下的球對稱形式，實現了直接的數值實現。
• 證明Hiscock-Lindblom非線性區的聲速下界由紅外聲速決定，確保了熱力學一致性。
• 推導結果確認，在平衡態下矩陣 ∂²u/∂Ya∂Yb 為非負定矩陣，這是反應流體框架中熱力學穩定性的必要條件。