[논문 리뷰] Pencil-Code special edition: modelling supernova driven turbulence

Frederick A. Gent, Mordecai‐Mark Mac Low외 3명|arXiv (Cornell University)|2018. 06. 05.

Solar and Space Plasma Dynamics참고 문헌 24인용 수 1

한 줄 요약

이 논문은 인위적 점성, 열전도도, 질량 확산 및 일반화된 쿠란트 조건을 도입하여 펜슬 코드(Pencil Code)를 초신성에 의해 촉발되는 난류를 모의할 수 있도록 개선한다. 이는 고마하 수(shock)를 정확하게 포착하기 위해 필요하다. 이 방법은 1차원 및 3차원 충격관과 초신성 잔해의 진화를 성공적으로 시뮬레이션하며, 등온 및 냉각 조건에서 분석적 및 반분석적 해와 일치한다.

ABSTRACT

High Mach number shocks are ubiquitous in interstellar turbulence. The Pencil Code is particularly well suited to the study of magnetohydrodynamics in weakly compressible turbulence and the numerical investigation of dynamos because of its high-order advection and time evolution algorithms. However, the high-order algorithms and lack of Riemann solver to follow shocks make it less well suited to handling high Mach number shocks, such as those produced by supernovae (SNe). Here, we outline methods required to enable the code to efficiently and accurately model SNe, using parameters that allow stable simulation of SN-driven turbulence, in order to construct a physically realistic galactic dynamo model. These include the resolution of shocks with artificial viscosity, thermal conductivity, and mass diffusion; the correction of the mass diffusion terms; and generalization of the Courant condition to include all source terms in the momentum and energy equations. We test our methods with the numerical solution of the 1D Riemann shock tube (Sod 1978), also extended to a 1D adiabatic shock with parameters and Mach number relevant to SN shock evolution, including shocks with radiative losses. We extend our test with the 3D numerical simulation of individual SN remnant evolution for a range of ambient gas densities typical of the interstellar medium and relate these to the analytical solutions of Sedov-Taylor (adiabatic) and the snowplough and Cioffi et al. (1988) results incorporating cooling and heating processes.

연구 동기 및 목표

• 초신성(SN) 폭발에서 흔히 나타나는 고마하 수 충격을 정확하게 모의할 수 있도록 펜슬 코드를 확장하기 위해, 기존에 리만 해법기구와 고차수 스킴이 부족하여 어려움을 겪어온 문제를 해결한다.
• 약한 압축성의 자기유체역학 유동에서 강력한 충격을 안정적이고 정확하게 해석할 수 있는 수치적 방법을 개발한다.
• 실제로 초신성에 의해 촉발되는 난류를 모의하기 위해 인위적 점성, 열전도도 및 질량 확산을 포함한 물리적 과정을 통합한다.
• 모멘텀 및 에너지 방정식의 모든 소스 항을 고려하여 쿠란트 조건을 일반화함으로써 안정성과 정확도를 향상시킨다.
• 1차원 충격관과 다양한 간성간 매질 밀도에서의 3차원 초신성 잔해 진화에 대해 분석적 해와의 비교를 통해 개선된 코드의 정당성을 검증한다.

제안 방법

• 기존 고차수 스킴이 실패하는 강력한 충격을 해결하기 위해 인위적 점성을 도입한다.
• 충격면에서 에너지 이행과 혼합을 모의하기 위해 열전도도 및 질량 확산 항을 구현한다.
• 보존 법칙과 물리적 거동에 일관성을 확보하기 위해 질량 확산 항을 수정한다.
• 모멘텀 및 에너지 방정식의 모든 소스 항을 고려하여 쿠란트 조건을 일반화함으로써 시간 스텝 안정성을 향상시킨다.
• 1차원 리만 충격관(Sod 1978) 및 등온 충격 테스트를 통해 수치 방법의 정당성을 검증한다.
• 3차원 시뮬레이션을 통해 다양한 주변 기체 밀도에서 개별 초신성 잔해를 시뮬레이션하고, 세도프-테일러 및 Cioffi 등(1988) 모델과 비교한다.

실험 결과

연구 질문

• RQ1 인위적 점성과 수정된 수치 스킴을 사용하여 펜슬 코드를 초신성에서 발생하는 고마하 수 충격을 정확하게 모의할 수 있는가?
• RQ2 열전도도와 질량 확산은 펜슬 코드에서 충격 시뮬레이션의 안정성과 정확도에 어떤 영향을 미치는가?
• RQ3 일반화된 쿠란트 조건은 강력한 소스 항이 존재하는 시뮬레이션에서 수치적 안정성을 어느 정도 향상시키는가?
• RQ4 3차원 초신성 잔해 시뮬레이션은 등온 및 냉각 조건에서 분석적 해를 얼마나 잘 재현하는가?
• RQ5 개선된 펜슬 코드는 간성간 매질의 밀도 범위에 걸쳐 실제적인 초신성에 의해 촉발되는 난류를 시뮬레이션할 수 있는가?

주요 결과

• 개선된 펜슬 코드는 인위적 점성, 열전도도 및 수정된 질량 확산을 통해 고마하 수 충격을 성공적으로 해석하며, 초신성에 의해 촉발되는 난류의 안정적 시뮬레이션을 가능하게 한다.
• 1차원 충격관 테스트(Sod 1978)는 분석적 해와 양호한 일치를 보이며, 수치 구현의 정당성을 검증한다.
• 초신성 관련 파arameter를 사용한 1차원 등온 충격 테스트는 강력한 충격을 고정밀도로 모의할 수 있음을 확인한다.
• 1차원 테스트에 복사 손실을 포함시킴으로써 물리적 현실성이 향상되고 기대되는 충격 구조와 일치한다.
• 3차원 초신성 잔해 시뮬레이션은 등온 조건에서는 세도프-테일러 해를 재현하며, 냉각 및 가열을 고려한 경우 눈보물 모양(snowplough) 및 Cioffi 등(1988) 모델과 일치한다.
• 일반화된 쿠란트 조건은 다양한 소스 항 크기에서 시간 스텝의 안정성과 시뮬레이션의 강건성을 크게 향상시킨다.