[논문 리뷰] Efficient Implementation of ADER Schemes for Euler and Magnetohydrodynamical Flows on Structured Meshes -- Comparison with Runge-Kutta Methods

Dinshaw S. Balsara, Chad Meyer외 3명|arXiv (Cornell University)|2010. 06. 10.

Computational Fluid Dynamics and Aerodynamics참고 문헌 47인용 수 5

한 줄 요약

이 논문은 구조적 메esh에서 Euler 및 자기유체역학(MHD) 방정식을 해결하기 위한 ADER(임의의 도함수) 방법의 실용적이고 효율적인 구현을 제시한다. 이는 전통적인 강안정성 보존 Runge-Kutta 스킴보다 뛰어난 성능을 보이는 단일 단계 시간 통합 방법을 제공한다. 저자들은 노드 공간과 모달 공간에서의 세부 수식을 제공하고, 발산 자유 MHD를 위한 수치적 플럭스와 전기장 계산을 최적화하며, 모든 시험된 정밀도 순서에서 Runge-Kutta 방법보다 거의 두 배 빠른 성능을 보임을 입증한다.

ABSTRACT

ADER (Arbitrary DERivative in space and time) methods for the time-evolution of hyperbolic conservation laws have recently generated a fair bit of interest. The ADER time update can be carried out in a single step, which is desirable in many applications. However, prior papers have focused on the theory while downplaying implementation details. The purpose of the present paper is to make ADER schemes accessible by providing two useful formulations of the method as well as their implementation details on three-dimensional structured meshes. We therefore provide a detailed formulation of ADER schemes for conservation laws with non-stiff source terms in nodal as well as modal space along with useful implementation-related detail. We also provide details for the efficient use of ADER schemes in obtaining the numerical flux for conservation laws as well as electric fields for divergence-free magnetohydrodynamics. An efficient WENO-based strategy for obtaining zone-averaged magnetic fields from face-centered magnetic fields in MHD is also presented. The schemes catalogued here have been implemented in the first author's RIEMANN code. The speed of ADER schemes is shown to be almost twice as fast as that of strong stability preserving Runge-Kutta time stepping schemes for all the orders of accuracy that we tested.

연구 동기 및 목표

• 이론적 ADER 스킴과 계산 유체역학에서의 실용적 구현 사이의 격차를 메우기 위해.
• 구조적 메쉬에서 노드 공간과 모달 공간에서의 ADER 방법에 대한 세부적이고 활용 가능한 수식을 제공하기 위해.
• 발산 자유 MHD 시뮬레이션을 위한 수치적 플럭스와 전기장 계산을 최적화하기 위해.
• 면 중심 데이터에서 영역 평균 자기장에 대한 WENO 기반 전략을 제시하기 위해.
• 강안정성 보존 Runge-Kutta 스킴에 비해 속도와 정밀도 측면에서 ADER의 계산적 우수성을 입증하기 위해.

제안 방법

• ADER 스킴은 Runge-Kutta 방법의 단계별 진행 방식을 피하기 위해 단일 단계 시간 갱신을 사용하여 구현된다.
• 이 방법은 3차원 구조적 메쉬에서 비강성 소스 항이 있는 쌍곡 보존 법칙을 해결하기 위해 노드 및 모달 수식을 활용한다.
• 수치적 플럭스는 보존 법칙에 맞게 조정된 고정밀도 라이만 해법을 사용해 효율적으로 계산된다.
• MHD에서의 전기장은 발산 자유 자기장 유지를 위해 플럭스 적분에서 직접 계산된다.
• 면 중심 자기장 데이터에서 영역 평균 자기장을 계산하기 위해 WENO 기반 재구성 전략이 사용된다.
• 모든 스킴은 RIEMANN 코드에 구현되어 고정밀도 및 성능 평가가 가능하도록 되어 있다.

실험 결과

연구 질문

• RQ1 계산 유체역학에서의 구조적 메쉬에 대해 ADER 스킴을 실용적으로 효율적으로 구현할 수 있는 방법은 무엇인가?
• RQ2 Euler 및 MHD 방정식에 대한 고성능 ADER 스킴을 실현하기 위한 핵심 구현 세부 사항은 무엇인가?
• RQ3 다양한 정밀도 순서에서 ADER의 계산 비용은 강안정성 보존 Runge-Kutta 스킴과 비교해 어떻게 되는가?
• RQ4 ADER 기반 MHD 시뮬레이션에서 자기장의 발산 자유 조건을 확보하기 위한 전략은 무엇인가?
• RQ5 ADER 스킴은 고정밀도를 유지하면서도 Runge-Kutta 방법 대비 거의 두 배의 계산 속도를 달성할 수 있는가?

주요 결과

• ADER 스킴은 모든 시험된 정밀도 순서에서 강안정성 보존 Runge-Kutta 스킴보다 거의 두 배 빠른 계산 속도를 달성한다.
• 노드 공간과 모달 공간에서의 ADER 구현은 구조적 메쉬에서 효율적이고 고정밀도 해를 제공한다.
• WENO 기반 재구성 전략의 사용은 면 중심 데이터에서 정확하고 안정적인 영역 평균 자기장을 보장한다.
• 수치적 플럭스와 전기장의 효율적 계산은 MHD에서 정확도 및 발산 자유 제약 조건을 유지하는 데 핵심적이다.
• ADER 방법의 단일 단계 시간 갱신은 다단계 Runge-Kutta 방법에 비해 뚜렷한 성능 우위를 제공한다.
• RIEMANN 코드 구현은 제안된 ADER 수식이 복잡한 유체역학 문제에 대해 강력하고 확장 가능하다는 것을 확인한다.