[논문 리뷰] Hydromagnetic waves in an expanding universe -- cosmological MHD code tests using analytic solutions

Thomas Berlok|arXiv (Cornell University)|2022. 03. 22.

Cosmology and Gravitation Theories참고 문헌 75인용 수 7

한 줄 요약

이 논문은 우주론적 MHD 코드를 테스트하기 위해 공동운동 좌표계에서 선형 수자외파의 해석적 및 수치적 해를 제시한다. 아인슈타인-데사이트르(EdS) 우주에서 알프레드파와 자기소음파에 대한 정확한 해를 유도하고, 이를 통해 arepo 코드를 검증하며, 충분한 해상도를 사용할 경우 뛰어난 일치를 보인다. 주요 기여는 척도 인자와 그 도함수와 관련된 미세한 코드 오류를 파악하기 위한 웨이브 감쇠 분 析 프레임워크를 제공하는 것이다.

ABSTRACT

We describe how analytic solutions for linear hydromagnetic waves can be used for testing cosmological magnetohydrodynamic (MHD) codes. We start from the comoving MHD equations and derive analytic solutions for the amplitude evolution of linear hydromagnetic waves in a matter-dominated, flat Einstein-de-Sitter (EdS) universe. The waves considered are comoving, linearly polarized Alfv\'en waves and comoving, magnetosonic (fast) waves modified by self-gravity. The solution for compressible waves is found for a general adiabatic index and we consider the limits of hydrodynamics without self-gravity in addition to the full solution. In addition to these analytic solutions, the linearized equations are solved numerically for a $\Lambda$CDM cosmology. We use the analytic and numeric solutions to compare with results obtained using the cosmological MHD code AREPO and find good agreement when using a sufficient number of grid points. We interpret the numerical damping clearly evident in simulations with few grid points by further deriving the Alfv\'en wave solution including physical Navier-Stokes viscosity. A comparison between Alfv\'en wave simulations and theory reveals that the dissipation can be described by a numerical viscosity coefficient $\eta_\mathrm{num} \propto a^{-5/2}$ where $a$ is the scale factor. We envision that our examples could be useful when developing a new cosmological MHD code or for regression testing of existing codes.

연구 동기 및 목표

• 척도 인자 $a$와 그 시간 도함수 $\dot{a}$에서 기인하는 오류를 특별히 타겟으로 삼는 우주론적 MHD 코드를 위한 철저한 테스트 프레임워크를 개발하기 위해.
• 기본적인 평탄한 아인슈타인-데사이트르(EdS) 우주에서 선형 유체역학적 자기파의 해석적 해를 제공하여 기준 기준선으로 삼기 위해.
• 더 넓은 적용성을 위해 $\Lambda$CDM 우주론에서 파동의 수치적 해를 포함하도록 테스트 세트를 확장하기 위해.
• 적은 격점 수를 가진 시뮬레이션에서의 수치적 감쇠가 $\eta_{\text{num}} \propto a^{-5/2}$와 비례하는 수치적 점성도로 해석될 수 있음을 보여주어 오류 진단이 가능하게 하기 위해.
• 이러한 파동 기반 테스트를 우주론적 MHD 코드의 자동 회귀 테스트 파이프라인에 통합할 것을 주장하기 위해.

제안 방법

• 선형화된 공동운동 MHD 방정식을 사용하여 아인슈타인-데사이트르(EdS) 우주론에서 공동운동 알프레드파와 중력 수정 자기소음파에 대한 해석적 해를 유도한다.
• 4차 룬게-쿠타 방법과 적응 시간 스텝을 사용하여 $\Lambda$CDM 우주론의 선형화된 방정식을 수치적으로 해결한다.
• 척도 인자 $a$와 그 도함수 $\dot{a}$를 도입하여 표준 MHD 방정식을 공동운동 좌표계로 변환하고, 운동량 및 유도 방정식을 적절히 수정한다.
• 일차원 주기적 영역에 파동 해를 구현하여 표준 하드웨어에서 빠르고 계산 비용이 낮은 테스트를 가능하게 한다.
• arepo 코드의 시뮬레이션 결과를 해석적 및 수치적 기준 해와 비교하여 정확도와 수렴성을 평가한다.
• 수치적 감쇠를 분석하기 위해 물리적 나비에-스토크스 점성도를 포함한 알프레드파 해를 유도하고, 이를 바탕으로 효과적 점성도 $\eta_{\text{num}} \propto a^{-5/2}$에 적합시킨다.

실험 결과

연구 질문

• RQ1 공동운동 좌표계에서 유체역학적 자기파의 해석적 해를 어떻게 도출할 수 있으며, 이를 우주론적 MHD 코드 테스트의 기준 해로 활용할 수 있는가?
• RQ2 아인슈타인-데사이트르 우주에서 알프레드파와 자기소음파의 진폭은 어떻게 행동하며, 척도 인자 $a$와 함께 어떻게 진화하는가?
• RQ3 이러한 파동의 수치적 시뮬레이션 결과가 $\Lambda$CDM 우주론에서 고정밀 수치 해와 어떻게 비교되는가?
• RQ4 저해상도 시뮬레이션에서의 수치적 감쇠는 양적 관계로 수치적 점성도 계수와 연결될 수 있는가? 그리고 이 점성도는 척도 인자에 어떻게 의존하는가?
• RQ5 이러한 파동 기반 테스트는 척도 인자 $a$와 $\dot{a}$를 다루는 코드 오류를 얼마나 잘 탐지할 수 있는가?

주요 결과

• 충분한 격점 해상도를 사용할 경우, EdS 우주에서의 알프레드파 및 자기소음파 해석적 해는 arepo 시뮬레이션과 뛰어난 일치를 보인다.
• 적은 격점 수를 가진 시뮬레이션에서의 수치적 감쇠는 $\eta_{\text{num}} \propto a^{-5/2}$와 비례하는 수치적 점성도로 잘 기술되며, 이는 유도된 물리적 스케일링과 일치한다.
• $\Lambda$CDM에 대한 arepo 시뮬레이션과 고정밀 수치 해의 비교에서, $a_i = 1/128$에서 $a = 1$까지의 우주 시간 범위가 더 길어지므로 파동의 진동 빈도가 EdS보다 더 높게 나타난다.
• arepo 코드에 의 intensionally 삽입된 버그(예: $\dot{a}$ 항의 잘못된 스케일링)는 해석적 해와의 비교를 통해 시뮬레이션 출력 결과를 분석함으로써, 명백한 수치적 실패가 없더라도 성공적으로 탐지되었다.
• 파동 기반 테스트 세트는 척도 좌표계에서 점성항의 잘못된 스케일링과 같은 미세한 코드 오류를 효과적으로 탐지할 수 있으며, 이는 일반적으로 발견되지 않을 수 있다.
• 저자들은 이러한 테스트가 우주론적 MHD 코드의 코드 신뢰도를 향상시키고 회귀를 방지하기 위해 자동 회귀 테스트 프레임워크에 통합되어야 한다고 결론 내린다.