[논문 리뷰] High-Order Discontinuous Galerkin Finite Element Methods with Globally Divergence-Free Constrained Transport for Ideal MHD

James A. Rossmanith|arXiv (Cornell University)|2013. 10. 16.

Computational Fluid Dynamics and Aerodynamics참고 문헌 30인용 수 3

한 줄 요약

이 논문은 이방성 자기유체역학(MHD)을 위한 고차수 불연속 갈레르킨 유한요소법을 제시하며, 새로운 제약 운반 기반 접근을 통해 전역적으로 수렴하지 않는 자기장 조건을 강제한다. 요소 중심의 자기벡터포ential을 도입하고, 이를 DG 스킴을 통해 갱신함으로써, 고차수 정확도로 모서리/면에 수직인 자기장 성분을 재구성함으로써 국소적 수렴하지 않는 조건과 전역적인 법선 자기장 투과성의 연속성을 모두 보장한다. 이로 인해 타원형 해법이나 자유 매개변수 없이 안정적이고 고차수 정확도를 갖는 시뮬레이션을 달성한다.

ABSTRACT

The modification of the celebrated Yee scheme from Maxwell equations to magnetohydrodynamics is often referred to as the constrained transport approach. Constrained transport can be viewed as a sort of predictor-corrector method for updating the magnetic field, where a magnetic field value is first predicted by a method that does not preserve the divergence-free condition on the magnetic field, followed by a correction step that aims to control these divergence errors. This strategy has been successfully used in conjunction with a variety of shock-capturing methods including WENO, central, and wave propagation schemes. In this work we show how to extend the basic CT framework to the discontinuous Galerkin finite element method on both 2D and 3D Cartesian grids. We first review the entropy-stability theory for semi-discrete DG discretizations of ideal MHD, which rigorously establishes the need for a magnetic field that satisfies the following conditions: (1) the divergence of the magnetic field is zero on each element, and (2) the normal components of the magnetic field are continuous across element edges/faces. In order to achieve such a globally divergence-free magnetic field, we introduce a novel CT scheme that is based on two ingredients: (1) we introduce an element-centered magnetic vector potential that is updated via a discontinuous Galerkin scheme on the induction equation; and (2) we define a mapping that takes element-centered magnetic field values and element-centered magnetic vector potential values and creates on each edge /face a representation of the normal component of the magnetic field; this representation is then mapped back to the elements to create a globally divergence-free element-centered representation of the magnetic field. For problems with shock waves, we make use of so-called moment-based limiters to control oscillations in the conserved quantities.

연구 동기 및 목표

• 이론적 MHD를 위한 고차수 불연속 갈레르킨 유한요소법을 개발하여 자기장의 수렴하지 않는 조건을 엄격히 강제한다.
• 2차원 및 3차원 카르테시안 격자에서 고차수 DG 방법에 대해 제약 운반(CT) 프레임워크를 확장하면서 고차수 정확도를 유지한다.
• 각 요소 내에서 국소적으로 수렴하지 않는 자기장과 요소 간면에서 연속적인 법선 자기장 성분을 모두 확보한다.
• 자기장 수렴 오차로 인한 수치적 불안정성, 즉 물리적으로 비현실적인 음수 압력 또는 밀도를 방지한다.
• 충격이 존재하는 상황에서 진동을 제어하기 위해 모멘트 기반 리미터를 통합함으로써 비선형 MHD 유동에서의 강건성을 향상시킨다.

제안 방법

• 유도 방정식을 푸는 불연속 갈레르킨 스킴을 통해 갱신되는 요소 중심의 자기벡터포텐셜을 도입한다.
• 벡터포텐셜 값들을 사용하여 요소의 모서리(2차원) 또는 면(3차원)에서 자기장 성분의 고차수 표현을 구성한다.
• 재구성된 모서리/면에 수직인 자기장 값을 다시 요소로 매핑하여 전역적으로 수렴하지 않는 요소 중심 자기장 표현을 생성한다.
• 예측-수정 전략을 활용한다: 먼저 표준 DG 진화를 통해 자기장을 예측하고, 이후 벡터포텐셜 기반 재구성으로 ∇·B = 0을 전역적으로 강제한다.
• 충격 근처에서 보조 양의 진동을 제어하기 위해 모멘트 기반 리미터를 사용하여 안정성을 유지한다.
• 각 요소에서의 수렴이 0이고, 요소 간면에서 법선 자기장 투과성이 연속되도록 하여 엔트로피 안정성을 확보한다.

실험 결과

연구 질문

• RQ1 제약 운반 프레임워크는 카르테시안 격자에서 고차수 불연속 갈레르킨 방법에 성공적으로 확장될 수 있는가?
• RQ2 스태그드 격자나 타원형 해법 없이 고차수 DG 맥락에서 전역적으로 수렴하지 않는 자기장을 어떻게 구성할 수 있는가?
• RQ3 자기벡터포텐셜은 고차수, 전역적으로 수렴하지 않는 자기장 재구성에 어떻게 기여하는가?
• RQ4 충격과 복잡한 유동 구조를 포함한 표준 MHD 테스트 케이스에서 제안된 방법의 성능은 어떠한가?
• RQ5 불연속성이 존재하는 상황에서 진동을 제어하면서도 고차수 정확도를 유지할 수 있는가?

주요 결과

• 각 요소 내에서 국소적으로 수렴하지 않는 조건과 요소 간면에서 법선 자기장 성분의 연속성을 강제함으로써 전역적으로 수렴하지 않는 자기장을 달성한다.
• 벡터포텐셜에서 재구성된 모서리/면에 수직인 자기장 성분은 고차수 다항식 근사가 사용되더라도 공간에서 2차 정확도를 보장한다.
• 부드러운 영역에서는 고차수 정확도를 유지하고, 충격 근처에서는 모멘트 기반 리미터를 통해 진동을 효과적으로 제어한다.
• 표준 MHD 테스트 케이스에 대한 수치 결과는 안정성과 강건성을 입증하며, 수렴 오차로 인한 음수 압력 또는 밀도가 관측되지 않았다.
• 프로젝션 방법에서 요구하는 타원형 해법이나 초월적 수치 해법에서의 자유 매개변수 없이 안정적이고 매개변수 없는 대안을 제공한다.
• 이 프레임워크는 2차원 및 3차원 카르테시안 메esh에 적용 가능하며, 모서리 및 면에서 고차수 자기장 성분을 재구성하기 위한 명시적 공식을 제공한다.