곡률 메쉬에서 오일러 방정식을 위한 엔트로피 안정 진동 억제 DGSEM

초록

**
본 논문은 2차원 압축성 오일러 방정식을 일반 곡률 메쉬에 적용할 수 있는 엔트로피 안정 고차 정확도 DGSEM을 제안한다. SBP 구조와 이산 기하학적 항등식을 이용해 반정밀 엔트로피 불평등을 보장하고, 진동 억제 DG(OEDG) 기법을 결합해 강한 충격파 근처의 Gibbs 현상을 효과적으로 억제한다. 새로운 충격 지표는 OEDG의 0차 감쇠계수를 활용해 셀 선택적으로 적용함으로써 계산 비용을 크게 절감한다.

**

상세 분석

**
이 연구는 두 가지 핵심 기술을 결합한다. 첫 번째는 Gauss‑Legendre‑Lobatto(GLL) 기반의 nodal DGSEM이 SBP(summation‑by‑parts) 성질을 만족하도록 설계된 점이다. 곡률 좌표 변환 시 발생하는 기하학적 항등식(metric identities)을 이산화 과정에서 정확히 보존하도록 메트릭 항을 재구성하고, 이를 통해 변환된 시스템이 자유 흐름(free‑stream) 보존과 동시에 엔트로피 안정성을 유지한다. 저자들은 엔트로피 변수와 포텐셜 함수를 이용해 엔트로피 보존 수치 플럭스를 정의하고, 요소 경계에서 이러한 플럭스를 사용함으로써 전역적인 엔트로피 불평등을 유도한다. 두 번째는 기존 OEDG 방법을 일반 곡률 메쉬에 적용 가능하도록 일반화한 것이다. 원래 OEDG는 직교 모달 기반의 삼각형 메쉬에만 적용 가능했으나, 여기서는 투영 연산자(projection operator)를 도입해 임의의 곡률 사각형 메쉬에서도 동일한 진동 억제 메커니즘을 구현한다. 특히, 0차 감쇠계수를 충격 지표로 재해석함으로써, 충격이 감지된 셀에만 OE 절차를 적용한다. 이는 전통적인 KXRCF 지표와 수학적으로 유사하면서도 구현이 간단하고 비용 효율적이다. 또한, 감쇠 연산을 의사시간(pseudo‑time) 문제로 풀어 정확히 적분함으로써 OFDG에서 발생하던 강성(stiffness) 문제를 회피한다. 전체 스키마는 다음과 같다: (1) SBP‑DGSEM으로 기본 고차 정확도와 엔트로피 안정성 확보, (2) 메트릭 항등식 보존을 위한 기하학적 처리, (3) OE 절차를 선택적 적용해 충격 영역에서만 인공 점성(damping) 추가, (4) 전역 보존성과 엔트로피 불평등 유지. 수치 실험에서는 정방형 및 비정형 곡률 메쉬에서의 소음 파동, 복잡한 충격 상호작용, 고마하톤 흐름 등을 테스트했으며, 모든 경우에서 기존 DGSEM 대비 진동이 현저히 감소하고, 엔트로피 감소가 이론적 한계 내에 머무르는 것을 확인했다. 이와 같이 본 논문은 고차 DG 방법의 두 가지 장벽—엔트로피 안정성 유지와 비물리적 진동 억제—을 동시에 해결한 최초의 곡률 메쉬용 프레임워크를 제공한다.

**