[논문 리뷰] Reliability of astrophysical jet simulations in 2D: On inter-code reliability and numerical convergence

Martin Krause, M. Camenzind|ArXiv.org|2001. 10. 22.

Computational Fluid Dynamics and Aerodynamics참고 문헌 12인용 수 19

한 줄 요약

이 연구는 NIRVANA MHD 코드를 사용하여 2차원 천체물리적 제트 시뮬레이션의 수치 수렴성과 코드 간 신뢰성을 평가한다. Lind et al. (1989) 및 Kössl & Müller (1988)의 두 개의 기존 코드와 결과를 비교한다. 수치적으로 수렴하는 경향이 있는 수소역학적 경우, 전역 제트 특성은 약 100 ppb에서 수렴하지만, 마그네틱 제트 시뮬레이션은 움직이는 막스 디스크로 인해 수렴하지 못하며, 고해상도 시뮬레이션은 질량 유인을 증가시키는 소규모 켈빈-헬름홀츠 및 레일리-테일러 모드와 같은 새로운 불안정성을 드러낸다.

ABSTRACT

In the present paper, we examine the convergence behavior and inter-code reliability of astrophysical jet simulations in axial symmetry. We consider both, pure hydrodynamic jets and jets with a dynamically significant magnetic field. The setups were chosen to match the setups of two other publications, and recomputed with the MHD code NIRVANA. We show that NIRVANA and the two other codes give comparable, but not identical results. We find that some global properties of a hydrodynamical jet simulation, like e.g. the bow shock velocity, converge at 100 points per beam radius (ppb) with NIRVANA. The situation is quite different after switching on the toroidal magnetic field: In this case, global properties converge even at 10 ppb. In both cases, details of the inner jet structure and especially the terminal shock region are still insufficiently resolved, even at our highest resolution of 70 ppb in the magnetized case and 400 ppb for the pure hydrodynamic jet. In the case of our highest resolution simulation, we can report two new features: First, small scale Kelvin-Helmholtz instabilities are excited at the contact discontinuity next to the jet head. This slows down the development of the long wavelength Kelvin-Helmholtz instability and its turbulent cascade to smaller wavelengths. Second, the jet head develops Rayleigh-Taylor instabilities which manage to entrain an increasing amount of mass from the ambient medium with resolution. This region extends in our highest resolution simulation over 2 jet radii in the axial direction.

연구 동기 및 목표

• NIRVANA MHD 코드를 사용하여 2차원 천체물리적 제트 시뮬레이션의 수치 수렴성을 평가하기 위해.
• NIRVANA, Lind et al. (1989), Kössl & Müller (1988)의 다수 코드 간 결과 비교를 통해 코드 간 신뢰성 평가를 위해.
• 인위적 점성과 수치 해상도가 제트 구조 및 충격 동역학 시뮬레이션에 미치는 영향을 조사하기 위해.
• 전역 제트 매개변수와 내부 구조(예: 종단 충격파 및 접촉 불연속면)의 수렴을 위한 해상도 요구 조건을 규명하기 위해.
• 수해상도에 따라 변화하는 조건에서 수소역학적 및 자기장이 있는 제트 시나리오에서 켈빈-헬름홀츠 및 레일리-테일러 불안정성의 발생 여부를 조사하기 위해.

제안 방법

• NIRVANA MHD 코드를 사용하여 다양한 공간 해상도(10에서 400점 per beam radius)로 벤치마크 제트 시뮬레이션(수소역학적 및 자기장이 있는 경우)을 재계산하였다.
• Lind et al. (1989) 및 Kössl & Müller (1988)의 이전 연구와 동일한 조건을 만족시키기 위해 축대칭(2D)을 적용하여 천체물리적 제트를 모델링하였다.
• 충격과 불연속면을 처리하기 위해 Rusanov 유량과 인위적 점성을 포함한 유한체적 MHD 해법기를 사용하였다.
• 다양한 해상도 수준에서 부두 충격파 속도 및 막스 디스크 위치와 같은 전역 특성의 수렴성을 추적하였다.
• 시간에 따른 시리즈 데이터와 등고선도를 분석하여 제트 머리와 코코넛 내에서의 불안정성 발생 및 구조적 진화를 파악하였다.
• 고해상도 시뮬레이션(수소역학적 경우 최대 400 ppb, 자기장이 있는 경우 최대 70 ppb)을 수행하여 해상도에 따라 변하는 특성, 예를 들어 소규모 난류 및 질량 유인과 같은 특징을 탐지하였다.

실험 결과

연구 질문

• RQ1 다양한 MHD 코드가 동일한 천체물리적 제트 설정에서 유사한 결과를 도출하는가?
• RQ2 수치 해상도가 부두 충격파 속도 및 막스 디스크 위치와 같은 전역 제트 매개변수의 수렴성에 어떤 영향을 미치는가?
• RQ3 특히 충격파 및 접촉 불연속면의 해상도에서 인위적 점성이 코드 간 출력 차이에 기여하는 역할은 무엇인가?
• RQ4 고해상도 시뮬레이션은 낮은 해상도에서 간과된 새로운 물리적 특성—소규모 불안정성 또는 증가한 질량 유인—을 드러내는가?
• RQ5 NIRVANA 코드를 사용한 자기장이 있는 제트 시뮬레이션에서 전역 매개변수의 수렴성은 보이지만, 움직이는 막스 디스크가 나타나는 이유는 무엇인가?

주요 결과

• NIRVANA 코드에서 수소역학적 제트 시뮬레이션의 전역 특성은 약 100점 per beam radius(ppb)에서 수렴한다.
• 자기장이 있는 제트 시뮬레이션은 전역 특성(예: 부두 충격파 형태, 속도)이 최소 10 ppb에서도 수렴하지만, 막스 디스크가 유입 경계 쪽으로 이동하여 수렴하지 않음을 나타낸다.
• 최고 해상도(수소역학적 경우 400 ppb, 자기장이 있는 경우 70 ppb)에서 소규모 켈빈-헬름홀츠 불안정성이 장파수 모드 이전에 나타나며, 이는 그 성장 억제와 난류 전이 지연를 초래한다.
• 최고 해상도 자기장이 있는 시뮬레이션에서 레일리-테일러 불안정성이 발생하여 주변 매질로부터 점점 더 많은 질량을 유인하며 축 방향으로 약 2개의 제트 반지름에 걸쳐 확장된다.
• 고해상도에도 불구하고 내부 제트 구조와 종단 충격 영역은 여전히 잘 해결되지 않으며, 수소역학적 경우 400 ppb에서도 난류 진화가 지속된다.
• 자기장 영향 하에서 특히 코드 간 수렴 행동이 다를 수 있음을 시사하는 FLOW와 NIRVANA 간의 격차—특히 NIRVANA에서 수렴하지 않는 막스 디스크—는 이전의 시뮬레이션 결과에 도전하며, 자기장 존재 하에서 코드에 따라 수렴 특성이 다를 수 있음을 시사한다.