[논문 리뷰] The relativistic Burgers equation on a FLRW background and its finite volume approximation

Tuba Ceylan, Philippe G. LeFloch외 1명|arXiv (Cornell University)|2015. 12. 26.

Computational Fluid Dynamics and Aerodynamics참고 문헌 6인용 수 3

한 줄 요약

이 논문은 압력이 없는 유체 가정 하에 상대론적 유체역학 방정식에서 유도된 Friedmann–Lemaître–Robertson–Walker (FLRW) 시공간 배경에서 상대론적 버거스 방정식을 수립한다. 이는 공간적으로 균일한 해를 유지하고 다양한 곡률 매개변수(k = -1, 0, 1)에서 충격파와 희박화를 정확하게 포착하는 잘 균형 잡힌 유한체적 스킴을 제안한다. 수치 결과는 k = -1일 때 특성 속도가 증가하여 더 빠른 파동 전파를 보임을 보여준다.

ABSTRACT

A relativistic generalization of the inviscid Burgers equation was proposed by LeFloch, Makhlof, and Okutmustur and then investigated on a Schwarzschild background. Here, we extend their analysis to a Friedmann-Lemaitre-Robertson-Walker (FLRW) background. This problem is more challenging due to the existence of non-trivial spatially homogeneous solutions. First, we derive the relativistic Burgers model of interest and determine its spatially homogeneous solutions. Second, we design a numerical scheme based on the finite volume methodology, which is well-balanced in the sense that spatially homogeneous solutions are preserved at the discrete level of approximation. Numerical experiments demonstrate the efficiency of the proposed method for weak solutions containing shock waves.

연구 동기 및 목표

• 우주적 팽창과 공간 곡률을 고려한 FLRW 배경 시공간에서 버거스 방정식의 상대론적 일반화를 유도하는 것.
• 곡률이 있는 배경에서 비자명한 공간적으로 균일한 해를 식별하고 분석하는 것.
• 이산 수준에서 균일한 해를 유지하는 잘 균형 잡힌 유한체적 수치 스킴을 설계하는 것.
• 다양한 곡률 매개변수(k = -1, 0, 1)에 대해 충격파와 희박화를 포함한 약한 해에 대해 스킴을 수치적으로 검증하는 것.

제안 방법

• 압력이 0인 유체 가정 하에 곡률이 있는 FLRW 시공간에서 상대론적 오일러 방정식으로부터 상대론적 버거스 방정식을 유도함으로써 기하학적 보존 법칙을 도출한다.
• 스케일 인자 a(t)와 곡률 매개변수 k를 사용하여 모델을 보존형으로 표현함으로써, 소스 항을 포함한 쌍곡형 균형 법칙을 도출한다.
• 수치적 플럭스와 재구성 기법을 기반으로 일阶 및 이阶 잘 균형 잡힌 유한체적 고두노프 유형 스킴을 구성한다.
• 수치 플럭스와 소스 항이 공간적으로 균일한 해와 일치하도록 하여 잘 균형 잡힌 성질을 확보한다.
• 셀 경계에서 수치 플럭스를 계산하기 위해 Rusanov 유형의 리만 해석기를 사용하며, 플럭스와 소스 항은 셀 평균에서 평가한다.
• k = -1, 0, 1에 대해 충격파와 희박화를 포함한 시험 케이스에 스킴을 적용하여 수렴성과 안정성을 평가한다.

실험 결과

연구 질문

• RQ1 비자명한 공간 곡률(k = -1, 0, 1)을 가진 FLRW 시공간에서 상대론적 버거스 방정식은 어떻게 행동하는가?
• RQ2 FLRW 배경에서 상대론적 버거스 모델의 공간적으로 균일한 해는 어떤 성질을 갖는가?
• RQ3 공간적으로 균일한 해를 유지하면서 충격파와 희박화를 정확히 포착할 수 있는 유한체적 스킴을 설계할 수 있는가?
• RQ4 곡률 매개변수 k는 수치 해에서 충격파와 희박화의 전파 속도와 구조에 어떤 영향을 미치는가?
• RQ5 시간에 따라 변화하는 곡률이 있는 배경에서 불연속 해에 대해 제안된 잘 균형 잡힌 스킴의 수렴성과 안정성은 어떻게 되는가?

주요 결과

• 시간에 따라 변화하는 스케일 인자 a(t)와 곡률 k 덕분에 FLRW 시공간에서 상대론적 버거스 방정식은 비자명한 공간적으로 균일한 해를 갖는다.
• 모델의 특성 속도는 (1 - kr²)^{1/2}로 주어지며, k가 감소할수록 증가하므로 k = -1일 때 더 빠른 파동 전파가 발생함을 설명한다.
• 제안된 유한체적 스킴은 이산 수준에서 균일한 해를 유지하는 잘 균형 잡힌 성질을 지녀, 천체물리학적 배경에서 장기간 시뮬레이션에 매우 중요하다.
• 수치 실험 결과는 모든 세 곡률 케이스에서 충격파와 희박화를 정확히 포착하는 스킴의 안정성과 효율성을 확인한다.
• 고두노프 스킴을 통해 이계 정확도를 달성하였으며, 희박화 및 충격 해 모두에 대해 효과적인 수렴이 이루어진다.
• 그림 3과 4의 빨간 선(k = -1)은 초록선(k = 0)과 파랑선(k = 1)보다 더 빠른 수렴과 전파를 보이며, 특성 속도의 분석 예측을 확인한다.