[논문 리뷰] A Simple and Effective High-Order Shock-Capturing Limiter for Discontinuous Galerkin Methods

Scott Moe, James A. Rossmanith외 1명|arXiv (Cornell University)|2015. 07. 10.

Computational Fluid Dynamics and Aerodynamics참고 문헌 7인용 수 8

한 줄 요약

이 논문은 불연속 갈레르킨(DG) 방법을 위한 단순하고 고차수의 충격 포착 리미터를 소개한다. 이 리미터는 Barth-Jespersen 및 Zhang-Shu의 최대 원칙 유지(MPP) 프레임워크의 아이디어를 융합한 것으로, 국소적으로 인접한 요소의 데이터를 사용해 해 변수를 제한하고 고차수 보정을 최소한으로 감쇠시킴으로써 충격 근처의 진동을 효과적으로 억제하면서도 스무스 영역에서는 고차수 정확도를 유지한다. 이는 구조적 및 비구조적 메esh에서의 1D 및 2D 벤치마크를 통해 입증되었다.

ABSTRACT

The discontinuous Galerkin (DG) finite element method when applied to hyperbolic conservation laws requires the use of shock-capturing limiters in order to suppress unphysical oscillations near large solution gradients. In this work we develop a novel shock-capturing limiter that combines key ideas from the limiter of Barth and Jespersen [AIAA-89-0366 (1989)] and the maximum principle preserving (MPP) framework of Zhang and Shu [Proc. R. Soc. A, 467 (2011), pp. 2752--2776]. The limiting strategy is based on traversing the mesh element-by-element in order to (1) find local upper and lower bounds on user-defined variables by sampling these variables on neighboring elements, and (2) to then enforce these local bounds by minimally damping the high-order corrections. The main advantages of this limiting strategy is that it is simple to implement, effective at shock capturing, and retains high-order accuracy of the solution in smooth regimes. The resulting numerical scheme is applied to several standard numerical tests in both one and two-dimensions and on both Cartesian and unstructured grids. These tests are used as benchmarks to verify and assess the accuracy and robustness of the method.

연구 동기 및 목표

• 불연속성 근처에서 물리적으로 비합리적인 진동을 피하면서도 강건하고 고차수의 충격 포착 리미터를 개발하는 것.
• 충격 근처에서 소산을 최소화하면서도 스무스 영역에서 고차수 정확도를 유지하는 것.
• 리미터 과정에 탐지 기능을 내장함으로써 명시적 충격 탐지 단계가 필요 없도록 하는 것.
• 구조적 및 비구조적 메쉬에서 모두 쉽게 구현되고 효과적인 방법을 확보하는 것.
• MPP 프레임워크의 적용 범위를 정규성 유지 초과에서 강력한 충격을 다룰 수 있도록 확장하는 것.

제안 방법

• 리미터는 요소 단위로 작동하며, 인접한 요소의 점별 값에서 사용자 정의 변수에 대한 국소적 상한 및 하한을 계산한다.
• 해의 품질을 평가하기 위해 부드러움 지표를 사용하며, 기울기가 α에 비례하는 h^1.5에 기반한 임계값을 초과할 경우에만 리미팅을 시행한다.
• 국소적 경계는 고차수 다항식 보정의 감쇠를 최소화함으로써 강제하며, 셀 평균을 유지하고 고차수 정확도를 유지한다.
• 엔트로피 레이어와 같은 복잡한 유동 특징이 있는 영역에서는 원시 변수 대신 보존 변수의 셀 평균을 사용하여 경계 추정을 수행한다.
• 감도가 높은 영역(예: 전방면 스텝 위쪽)에서는 보다 낮은 잔류 엔트로피 레이어를 억제하기 위해 수정된 리미터 파라미터 α = 20h^1.5를 사용한다.
• SSPRK4 시간 적분과 LLF 수치적 플럭스를 사용하며, 직각좌표계 및 비구조적 메쉬에서 P3 요소를 적용하여 구현한다.

실험 결과

연구 질문

• RQ1 명시적 충격 탐지 없이도 스무스 영역에서 정확도를 유지하면서도 단순하고 고차수의 DG 리미터를 설계할 수 있는가?
• RQ2 Barth-Jespersen 및 MPP 원칙의 조합은 정확도의 순서를 떨어뜨리지 않고 충격 근처의 진동을 얼마나 효과적으로 억제하는가?
• RQ3 이 리미터는 수정이 거의 필요 없이 구조적 및 비구조적 메쉬에서 안정적으로 적용될 수 있는가?
• RQ4 엔트로피 레이어가 존재하는 전방면 스텝 문제와 같은 복잡한 충격 구조를 어떻게 잘 포착하는가?
• RQ5 고기울기 영역에서 경계 추정에 보존 변수를 원시 변수 대신 사용할 경우의 영향은 무엇인가?

주요 결과

• 1D 및 2D 유동 방정식에서 충격 근처의 물리적으로 비합리적인 진동을 효과적으로 억제하였으며, 특히 도전적인 전방면 스텝 문제에서도 성공적으로 작동하였다.
• 전방면 스텝 문제에서는 특히 스텝 상단 영역에서 α = 20h^1.5로 조정된 정밀한 리미터 파라미터를 사용함으로써 임의의 엔트로피 레이어를 효과적으로 완화하였다.
• 직각좌표계 메쉬 시뮬레이션에서는 뚜렷한 엔트로피 레이어가 관측되었지만, 이는 이전 연구에서 관측된 일반적인 범위 내에 있었으며, 리미터의 강건성을 시사한다.
• 비구조적 메쉬 결과는 직각좌표계 결과와 뛰어난 일치를 보였으며, 이는 리미터가 비구조적 메쉬에서도 효과적임을 확인한다.
• 접촉 불연속의 날카로운 해상도와 스무스 영역의 매끄러운 프로파일을 통해 리미터가 스무스 영역에서 고차수 정확도를 유지하고 있음을 입증하였다.
• 리미터는 추가적인 충격 탐지 단계가 필요 없었으며, 리미팅 결정이 국소적 경계 추정 과정에 내재되어 있기에 자연스럽게 작동하였다.