[논문 리뷰] Monolithic Convex Limiting for Legendre-Gauss-Lobatto Discontinuous Galerkin Spectral Element Methods

Andrés M. Rueda-Ramírez, Benjamin Bolm외 2명|arXiv (Cornell University)|2023. 03. 01.

Advanced Numerical Methods in Computational Mathematics인용 수 4

한 줄 요약

이 논문은 레전드르–갈루아–로바트로프 비연속 갈레르킨 스펙트럴 요소 방법(Legendre–Gauss–Lobatto discontinuous Galerkin spectral element methods, LGL-DGSEM)을 위한 단일체 볼록 제한(MCL) 프레임워크를 소개한다. LGL 노드가 지닌 하위셀 유량 구조와 합산-통합-성질을 활용하여 불변 도메인 유지 및 엔트로피 안정성을 보장한다. 이 방법은 시간 스텝에 의존하지 않는 안정화를 달성하여, 인위적 난류에 의존하지 않고 강한 충격파, 급격한 기울기, 소용루어가 지배하는 영역을 포함한 압축성 유동의 견고한 시뮬레이션을 가능하게 한다.

ABSTRACT

We extend the monolithic convex limiting (MCL) methodology to nodal discontinuous Galerkin spectral element methods (DGSEM). The use of Legendre-Gauss-Lobatto (LGL) quadrature endows collocated DGSEM space discretizations of nonlinear hyperbolic problems with properties that greatly simplify the design of invariant domain preserving high-resolution schemes. Compared to many other continuous and discontinuous Galerkin method variants, a particular advantage of the LGL spectral operator is the availability of a natural decomposition into a compatible subcell flux discretization. Representing a high-order spatial semi-discretization in terms of intermediate states, we perform flux limiting in a manner that keeps these states and the results of Runge-Kutta stages in convex invariant domains. Additionally, local bounds may be imposed on scalar quantities of interest. In contrast to limiting approaches based on predictor-corrector algorithms, our MCL procedure for LGL-DGSEM yields nonlinear flux approximations that are independent of the time-step size and can be further modified to enforce entropy stability. To demonstrate the robustness of MCL/DGSEM schemes for the compressible Euler equations, we run simulations for challenging setups featuring strong shocks, steep density gradients and vortex dominated flows.

연구 동기 및 목표

• 레전드르–갈루아–로바트로프(Legendre–Gauss–Lobatto, LGL) 구분을 사용하는 노달 비연속 갈레르킨 스펙트럴 요소 방법에 단일체 볼록 제한(MCL)을 확장한다.
• LGL-DGSEM에서의 하위셀 유량 분해와 대각 행렬 질량을 활용하여 고해상도, 불변 도메인 유지 가능 시스템의 설계를 단순화한다.
• 런지–쿠타 단계 전반에 걸쳐 볼록 불변 도메인을 유지하고 스칼라 양의 국소적 범위를 설정할 수 있는 유량 제한 전략을 개발한다.
• 시간 스텝 크기와 무관한 MCL 기반 유량 수정을 통해 엔트로피 안정성을 확보한다.
• 강한 충격파, 밀도 기울기, 난류 유동을 포함한 도전적인 초구형 문제에서의 견고성을 입증한다.

제안 방법

• 중간 상태를 통해 고차원 공간 반연속 이산화를 표현하여, 볼록 불변 도메인을 유지하면서도 유량 제한을 가능하게 한다.
• LGL-DGSEM에 내재된 자연스러운 하위셀 유량 분해를 활용하여 인위적 재구성 없이 반-확산 보정을 정의한다.
• 예측-수정 시간 분할을 피하기 위해 제한된 유량을 잔여항에 직접 통합하는 단일체적 접근을 구현한다.
• 유량 합계가 아닌 개별 유량을 제약하여 불변 도메인 유지 조건을 강제함으로써 국소적 및 전역적 범위를 보장한다.
• 제한기 기반 수정을 통해 반연속 엔트로피 부등식을 유지함으로써 엔트로피 안정성을 확보한다.
• LGL-DGSEM에서 대각 행렬 질량과 합산-통합(SBP) 연산자를 활용하여 엔트로피 안정성과 고차원 정확도를 지원한다.

실험 결과

연구 질문

• RQ1 단일체 볼록 제한이 LGL-DGSEM에 효과적으로 확장되어 불변 도메인 유지 및 엔트로피 안정성을 달성할 수 있는가?
• RQ2 시간 스텝 의존성과 난류 측면에서 MCL은 전통적인 예측-수정 FCT 기법과 비교해 어떻게 다를까?
• RQ3 MCL/LGL-DGSEM의 성능은 강한 충격파, 급격한 기울기, 소용루어 지배 유동을 포함한 기준 문제에서 어떻게 나타나는가?
• RQ4 안정성에 영향을 주지 않으면서 MCL 프레임워크 내에서 스칼라 양의 국소적 범위를 강제로 설정할 수 있는가?
• RQ5 단일체 형식이 비선형 초구형 시스템에서 시간 스텝에 의존하지 않는 안정화를 보장하는가?

주요 결과

• MCL/LGL-DGSEM 프레임워크는 시간 스텝에 의존하지 않는 안정화를 달성하여, 수치적 안정화가 필요한 문제에서 기대하는 시간 수렴성을 확보한다.
• 이 방법은 임의의 진동이나 도메인 위반 없이 강한 충격파, 급격한 밀도 기울기, 소용루어 지배 유동을 성공적으로 포착한다.
• 예측-수정 FCT 기법과 비교해 MCL는 시간 스텝 크기 의존성 난류를 감소시켜 시간 수렴 특성을 향상시킨다.
• LGL 노드의 사용은 자연스러운 하위셀 유량 분해와 대각 행렬 질량을 가능하게 하여 볼록 제한의 구현을 단순화한다.
• 엔트로피 안정성은 MCL 기반 유량 수정을 통해 강제로 구현되며, 제한 전략을 통해 반연속 엔트로피 부등식이 유지된다.
• Trixi.jl에서의 구현은 재현 가능하며 공개되어 있어 코드와 결과의 완전한 복제를 지원한다.