[논문 리뷰] Simulating bulk viscosity in neutron stars. I. Formalism

Giovanni Camelio, Lorenzo Gavassino외 3명|arXiv (Cornell University)|2022. 04. 25.

Pulsars and Gravitational Waves Research참고 문헌 63인용 수 3

한 줄 요약

이 논문은 일반 상대론적 유체역학 하에서 중성자별의 부피 점성도를 시뮬레이션하기 위해 다성분 반응 유체, Hiscock-Lindblom(HIS) 제2차 Müeller-Israel-Stewart 이론, 그리고 Maxwell-Cattaneo 선형 근사 모델의 세 가지 이론적 프레임워크를 개발하고 비교한다. 구형 대칭에서 동역학 방정식의 전반적인 집합을 유도하고, 부피 점성 유체와 반응 혼합물 간의 수학적 이중성(duaity)을 확립하며, 비선형 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 모델이라는 세 가지 다른 접근법을 사용하여 중성자별 내의 부피 점성도를 모델링하는 것.
• Schwarzschild 좌표계(반경 방향 게이지-극좌락)를 사용하여 이러한 모델의 유체역학 방정식을 구형 대칭으로 특수화하는 것.
• 부피 점성 유체와 반응 혼합물 간의 수학적 이중성을 확립하고, 다성분 유체 모델이 Müller-Israel-Stewart 유체역학에 대한 수치적으로 실현 가능한 대안임을 검증하는 것.
• 비선형 Hiscock-Lindblom 이론 영역에서의 효과적 음속을 계산하여 수치적 안정성과 물리적 정확성에 필수적인 요소를 확보하는 것.

제안 방법

• 다성분 반응 유체에 대한 일반 상대론적 유체역학 방정식 전체를 유도하며, 반응 중 중성자 방출률과 에너지 손실을 포함한다.
• 부피 점성도를 모델링하기 위해 Müller-Israel-Stewart 체계를 적용하며, 특히 제2차 Hiscock-Lindblom 이론과 그 선형 근사인 Maxwell-Cattaneo 근사를 구현한다.
• 모든 세 모델을 Schwarzschild 좌표계(반경 방향 게이지-극좌락)를 사용하여 구형 대칭으로 특수화하여 1차원 수치 구현을 가능하게 한다.
• 협업 수치 작업에 사용할 수 있도록 특정 상태 방정식과 입자 반응 속도(우르카 과정 기반)를 도입한다.
• 특성선 방법을 사용하여 비선형 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/∂ρ)ₛ 조건을 만족할 경우 부피 점성 이론과 일치함을 확인하여 수학적 이중성을 검증한다.
• 다성분 유체 모델은 중성자 방출률을 포함하도록 확장되어 중성자별 시뮬레이션에서 현실적인 에너지 손실 모델링을 가능하게 한다.
• 모든 세 모델의 유체역학 방정식이 Schwarzschild 좌표계에서 완전히 유도되고, 구형 대칭으로 특수화되어 직접적인 수치적 구현이 가능해졌다.
• 비선형 Hiscock-Lindblom 영역에서의 음속이 적외선 음속으로부터 아래bound되어 있음을 입증하여 열역학적 일관성을 보장한다.
• 유도 과정을 통해 평형 상태에서 행렬 ∂²u/∂Ya∂Yb가 음이 아닌 정의를 갖는 것으로 확인되어 반응 유체 프레임워크에서 열역학적 안정성에 필수적인 조건을 충족함을 확인한다.