[논문 리뷰] Investigating prescriptions for artificial resistivity in smoothed particle magnetohydrodynamics

James Wurster, Matthew R. Bate외 2명|arXiv (Cornell University)|2017. 06. 23.

Fluid Dynamics Simulations and Interactions참고 문헌 31인용 수 8

한 줄 요약

이 연구는 스무스드 파article 자기유체역학(SPMHD)에서 Price 등(2017), Tricco & Price(2013)의 항 평균화 및 가변 평균화 방식의 세 가지 인공 점성도 규정을 Orszag-Tang 소용파와 고해상도 별 형성 시뮬레이션을 통해 평가한다. Price 등 점성도는 최소한의 분산을 동반하면서 자기 구조를 가장 잘 해석하는 것으로 나타났으며, 별 형성 시뮬레이션에서는 기체 기포, 밀도 높은 디스크 또는 이동하는 싱크홀 입자 등 서로 다른 비수렴 결과를 유도하여, 인공 점성도가 결과에 상당한 편향을 주며 물리적 점성도가 수렴성을 확보하는 데 필수적임을 시사한다.

ABSTRACT

In numerical simulations, artificial terms are applied to the evolution equations for stability. To prove their validity, these terms are thoroughly tested in test problems where the results are well known. However, they are seldom tested in production-quality simulations at high resolution where they interact with a plethora of physical and numerical algorithms. We test three artificial resistivities in both the Orszag-Tang vortex and in a star formation simulation. From the Orszag-Tang vortex, the Price et. al. (2017) artificial resistivity is the least dissipative thus captures the density and magnetic features; in the star formation algorithm, each artificial resistivity algorithm interacts differently with the sink particle to produce various results, including gas bubbles, dense discs, and migrating sink particles. The star formation simulations suggest that it is important to rely upon physical resistivity rather than artificial resistivity for convergence.

연구 동기 및 목표

• 다양한 인공 점성도 규정이 SPMHD 시뮬레이션의 수치적 안정성과 정확성에 미치는 영향을 평가하기 위해.
• 표준 테스트 문제를 초월한 고해상도, 생산 수준의 시뮬레이션에서 인공 점성도 모델이 수렴 가능한 결과를 도출하는지 테스트하기 위해.
• 별 형성 시뮬레이션에서 점성도 선택에 따라 기체 기포, 밀도 높은 디스크 또는 이동하는 싱크홀 입자 등 서로 다른 결과가 발생하는 이유를 조사하기 위해.
• 점성도 규정을 수정함으로써(예: 속도 대신 운동량을 사용함) 기체 기포와 같은 비물리적 아티팩트를 억제할 수 있는지 평가하기 위해.
• 해상도에 독립적이고 물리적으로 타당한 결과를 확보하기 위해 향후 시뮬레이션에 물리적 점성도를 통합할 것을 주장하기 위해.

제안 방법

• 연구는 Phantom SPMHD 코드에 Price 등(2017), Tricco & Price(2013)의 항 평균화 및 가변 평균화 점성도 규정을 구현한다.
• Orszag-Tang 소용파 테스트를 통해 고해상도에서 자기 도시의 해석 능력과 분산 최소화 능력을 평가한다.
• 싱크홀 입자 동역학과 자기장 진화를 포함한 장기적 거동 평가를 위해 고해상도 별 형성 시뮬레이션을 수행한다.
• 신호 속도 계산에서 속도 대신 운동량을 사용함으로써 점성도 규정을 수정하여 수치적 아티팩트를 테스트한다.
• 기체 기포 형성 여부를 시뮬레이션 전반에서 모니터링하여, 아티팩트의 근본 원인이 인공 점성도인지 수치적 해상도인지 규명한다.
• 향후 시뮬레이션에서 해상도에 독립적일 뿐 아니라 물리적으로 타당한 대안으로 물리적 점성도를 권장한다.

실험 결과

연구 질문

• RQ1 Orszag-Tang 소용파 테스트에서 서로 다른 인공 점성도 규정은 자기 구조의 해상도에 어떤 영향을 미치는가?
• RQ2 고해상도 별 형성 시뮬레이션에서 인공 점성도 모델은 수렴 가능한 결과를 도출하는가, 아니면 상반된 물리적 결과를 초래하는가?
• RQ3 점성도 규정을 수정함으로써(예: 속도 대신 운동량 사용) 비물리적 기체 기포 형성을 억제할 수 있는가?
• RQ4 왜 서로 다른 점성도 모델이 별 형성 시뮬레이션에서 기체 기포, 밀도 높은 디스크 또는 이동하는 싱크홀 입자 등 서로 다른 결과를 초래하는가?
• RQ5 SPMHD 시뮬레이션에서 수렴성과 물리적으로 의미 있는 결과를 달성하기 위해 물리적 점성도가 필수적인가?

주요 결과

• Price 등(2017)의 인공 점성도는 가장 낮은 분산을 보이며, 테스트된 해상도에서 자기 도시를 성공적으로 해석한다.
• 항 평균화 및 가변 평균화 점성도 모델의 총 자기 에너지는 Price 등 모델 대비 최대 16% 높아, 자기 구조의 해상도가 열악함을 시사한다.
• 별 형성 시뮬레이션에서 Price 등 점성도는 싱크홀 입자 근처에 기체 기포를 생성하고, 항 평균화 모델은 큰 밀도 높은 디스크를 형성하며, 가변 평균화 모델은 싱크홀 입자의 수직 이동을 유도한다.
• 신호 속도 계산에서 속도 대신 운동량을 사용하도록 수정하면 기체 기포 형성이 약 0.02 자유낙하 시간 지연되지만, Price 등 점성도가 지배적일 경우 이를 완전히 방지하지 못한다.
• 모든 테스트된 인공 점성도 모델은 별 형성 시뮬레이션에서 수렴하지 않는 결과를 도출하여, 수치적 아티팩트를 유발하고 시뮬레이션 결과에 편향을 주는 것으로 나타났다.
• 연구는 물리적 점성도가 해상도에 독립적이고 물리적으로 타당하므로, 향후 시뮬레이션에서 인공 점성도보다 우선시되어야 한다고 결론 내린다.