Positivity-Preserving and Entropy-Stable Oscillation-Eliminating DGSEM for the Compressible Euler Equations on Curvilinear Meshes with Adaptive Mesh Refinement

논문 요약: 비정합 격자 (Nonconforming Meshes) 를 가진 곡선 격자에서의 적응적 격자 세분화 (AMR) 를 위한 엔트로피 안정성 및 양성 유지 DGSEM

1. 연구 배경 및 문제 제기

• 배경: 고차 엔트로피 안정성 (Entropy-Stable) 불연속 갤러킨 (DG) 방법은 복잡한 기하학적 구조를 가진 압축성 유동 시뮬레이션에 강력한 도구입니다. 이전 연구 [1] 에서 저자들은 정합 (conforming) 곡선 격자에서 엔트로피 안정성과 진동 제거 (Oscillation-Eliminating, OE) 를 결합한 DGSEM 을 개발했습니다.
• 문제점: 적응적 격자 세분화 (AMR) 를 도입하면 격자 크기나 차수가 다른 요소들이 만나는 비정합 (nonconforming) 인터페이스가 발생합니다. 기존 엔트로피 안정성 DGSEM 을 AMR 환경으로 확장할 때, 특히 양성 유지 (positivity preservation, 밀도와 압력이 음수가 되지 않음) 와 고차 정확도를 동시에 보장하는 것은 여전히 해결되지 않은 난제였습니다.
• 목표: 곡선 격자 기반 AMR 환경에서 엔트로피 안정성, 양성 유지, 진동 제거 기능을 모두 갖춘 DGSEM 프레임워크를 확장하고, 비정합 인터페이스에서의 수치 플럭스 (numerical flux) 설계 문제를 해결하는 것입니다.

2. 주요 방법론 및 제안된 기법

이 논문은 2:1 격자 세분화 제약 조건 하에서 2 차원 곡선 사각형 격자를 대상으로 하며, 다음과 같은 핵심 기법들을 제안합니다.

가. 비정합 인터페이스를 위한 두 가지 수치 플럭스 설계
저자들은 비정합 인터페이스 처리를 위해 두 가지 상반된 접근법을 제안하고 그 특성을 분석했습니다.

1. 엔트로피 안정성 플럭스 (Entropy-Stable Flux):

   • 구조: 인터페이스에서 엔트로피 불등식을 만족하고 전역 보존성을 보장하도록 설계되었습니다.
   • 한계: 다항식 차수 $N \ge 2$일 때, 비정합 보간 연산자 (interpolation operators) 에 음수 성분이 존재하게 됩니다. 엔트로피 안정성과 양성 유지를 증명하기 위해 이 음수 성분의 부호를 보정 (sign correction) 해야 하지만, 이로 인해 형식적 고차 정확도 (formal high-order consistency) 가 손실됩니다. 즉, $N \ge 2$에서는 플럭스가 최대 1 차 정확도까지만 보장됩니다.
   • 특징: $N=1$ (2 차 정확도) 에서는 음수 성분이 없어 정확도와 안정성을 동시에 달성합니다.

2. 모르타르 기반 플럭스 (Mortar-Based Flux):

   • 구조: 고해상도 (fine) 측의 모르타르 (mortar) 노드에서 해를 보간한 후, 표준적인 2 점 플럭스 (local Lax-Friedrichs) 를 계산하는 고전적인 접근법입니다.
   • 장점: 엔트로피 안정성 증명에는 실패하지만, 고차 정확도를 유지하며 양성 유지가 가능합니다.
   • 절충 (Trade-off): 엔트로피 안정성 (Entropy Stability) 과 고차 정확도 (High-order Accuracy) 는 비정합 인터페이스에서 동시에 증명하기 어렵다는 근본적인 트레이드오프를 드러냈습니다.

나. 곡선 AMR 격자를 위한 양성 유지 프레임워크 확장

• Zhang-Shu Limiter: Zhang-Shu 의 양성 유지 스케일링 리미터를 곡선 비정합 격자에 엄밀하게 확장했습니다.
• 이론적 증명: 전진 오일러 시간 적분과 적절한 CFL 조건 하에서, 제안된 두 가지 플럭스 (엔트로피 안정성 또는 모르타르 기반) 모두 셀 평균 밀도와 압력의 양성을 보장함을 증명했습니다.
• 데이터 전송: AMR 격자 간 (세분화 및 조대화) 데이터 전송 시에도 양성을 유지하면서 보존성을 확보하는 절차를 개발했습니다.

다. 진동 제거 (Oscillation-Eliminating, OE) 및 AMR 전략

• OE 절차: [14] 에서 제안된 진동 제거 절차를 비정합 격자에 적용했습니다. 비정합 인터페이스에서는 양쪽의 해를 공통의 노드 집합으로 보간한 후 점프 (jump) 를 계산하여 감쇠 계수를 결정합니다.
• AMR 전략: 충격 지표 (shock indicator) 를 기반으로 격자를 세분화하거나 조대화합니다. 세분화 시에는 다항식 제한 (restriction) 을, 조대화 시에는 $L^2$ 투영을 사용하여 데이터를 전송하며, 각 단계마다 Zhang-Shu 리미터를 적용하여 해의 물리적 허용 가능성 (admissibility) 을 유지합니다.

3. 주요 결과 및 수치 실험

MFEM 프레임워크를 사용하여 다양한 수치 실험을 수행하여 제안된 방법의 성능을 검증했습니다.

• 등엔트로피 와류 문제 (Isentropic Vortex):
  • 비정합 곡선 격자에서의 수렴성을 검증했습니다.
  • 결과: 엔트로피 안정성 플럭스는 $N \ge 2$에서 정확도가 저하되는 것을 확인했으며, 모르타르 기반 플럭스는 설계된 고차 정확도 ($N$ 차) 를 유지함을 보였습니다.
• 더블 마하 반사 (Double Mach Reflection):
  • 강한 충격파와 복잡한 유동 구조를 가진 문제입니다.
  • 결과: 모르타르 기반 플럭스 ($N=2$) 를 사용한 AMR 시뮬레이션은 정합 격자 시뮬레이션보다 더 풍부한 미세 구조 (슬립 라인, 하류 영역) 를 포착했으며, 엔트로피 안정성 플럭스 ($N=1$) 도 안정적으로 작동했습니다.
• 고마하 천체 제트 (High Mach Astrophysical Jet):
  • 밀도와 압력이 0 에 매우 근접하는 극한 조건에서의 양성 유지 능력을 검증했습니다.
  • 결과: 모르타르 기반 플럭스와 Zhang-Shu 리미터의 조합이 AMR 환경에서도 밀도와 압력이 음수가 되지 않도록 성공적으로 유지하여 계산이 붕괴되지 않음을 보였습니다.
• 원통 주위의 초음속 유동 (Supersonic Flow past Cylinder):
  • 곡선 격자에서의 AMR 효율성을 검증했습니다.
  • 결과: 켈빈 - 헬름홀츠 (Kelvin-Helmholtz) 불안정성과 같은 미세한 와류 구조를 AMR 을 통해 효율적으로 해상도하며, 격자 수가 적어도 정확한 결과를 얻었습니다.

4. 핵심 기여 및 의의

1. 이론적 통찰: 비정합 인터페이스에서 엔트로피 안정성, 양성 유지, 고차 정확도라는 세 가지 목표를 동시에 달성하는 것이 $N \ge 2$에서 불가능할 수 있음을 이론적으로 증명하고 분석했습니다. 이는 향후 연구 방향을 설정하는 중요한 통찰을 제공합니다.
2. 실용적 알고리즘: 곡선 격자 AMR 환경에서 양성 유지가 보장되는 완전한 DGSEM 알고리즘을 개발했습니다. 이는 충격파가 포함된 복잡한 압축성 유동 시뮬레이션에서 계산의 견고성 (robustness) 을 크게 향상시킵니다.
3. 모르타르 플럭스의 유효성: 엔트로피 안정성 증명 없이도 고차 정확도와 양성 유지를 동시에 달성할 수 있는 모르타르 기반 플럭스의 실용적 가치를 입증했습니다.
4. 확장성: 진동 제거 (OE), 충격 기반 AMR, 그리고 데이터 전송 절차를 통합하여 실제 공학적 문제에 적용 가능한 강력한 프레임워크를 제시했습니다.

5. 결론

이 논문은 곡선 격자 기반 AMR 환경에서 DGSEM 의 적용 가능성을 크게 확장했습니다. 엔트로피 안정성과 고차 정확도 사이의 본질적인 트레이드오프를 인정하면서도, 양성 유지가 보장되는 모르타르 기반 플럭스를 통해 고차 정확도와 견고성을 동시에 확보하는 실용적인 해결책을 제시했습니다. 이는 복잡한 기하학적 구조와 강한 충격파가 공존하는 압축성 유동 해석에 있어 매우 중요한 발전입니다.