볼륨 항 적응을 이용한 고차 DG 스킴의 효율·안정성 향상

본 논문은 고차 불연속 Galerkin(DG) 방법을 이용한 시간 의존 편미분 방정식(PDE) 해석에서, 볼륨 항의 이산화 방식을 동적으로 전환하는 새로운 적응 기법인 “v‑adaptivity”(볼륨 항 적응)를 제안한다. 전통적인 DG 스킴은 두 가지 주요 구성 요소, 즉 면(term)과 볼륨(term)으로 이루어진다. 면 항은 수치적 리만 해석을 통해 처리되며, 다양한 수치 플럭스가 존재한다. 반면 볼륨 항은 설계 자유도가 큰데, 엔트로피 보존·안정을 위해 특별히 설계된 flux‑differencing(FD) 형태가 널리 사용된다. FD 형태는 두 점 플럭스를 이용해 각 노드 쌍 사이의 상호작용을 정확히 포착함으로써 엔트로피 보존을 보장하지만, 차수가 커질수록 연산량이 급격히 증가한다. 특히 텐서곱 구조를 갖지 않는 삼각형·프리즘·피라미드와 같은 요소에서는 모든 노드 쌍에 대해 두 점 플럭스를 계산해야 하므로 비용이 p²·d 수준으로 급증한다. 이에 저자들은 “v‑adaptivity”라는 개념을 도입한다. 핵심 아이디어는 매 러크‑쿠타(RK) 단계마다 각 요소에 대해 현재 해의 특성을 평가하고, 그 결과에 따라 볼륨 항을 약한 형태(Weak Form, WF)와 FD 형태 중 하나로 선택하는 것이다. WF는 전통적인 약한 형태로, 각 노드당 물리적 플럭스만 한 번 평가하면 되므로 연산 비용이 매우 낮다. 반면 FD는 엔트로피 보존을 위한 두 점 플럭스를 사용해 비용이 높지만, 강인성을 제공한다. 전환 기준은 엔트로피 생산량에 기반한다. 엔트로피 변수 w = ∂S/∂u(여기서 S는 수학적 엔트로피)와 시간 미분 ˙u를 이용해 각 형태에서 발생하는 엔트로피 생산 ˙S_WF와 ˙S_FD를 계산한다. 구체적으로, - WF에서의 엔트로피 생산: ˙S_WF = -J⁻¹ M⁻¹ Dᵀ M f, - FD에서의 엔트로피 생산: ˙S_FD = -J⁻¹ (2D - M⁻¹ B) f_EC, 여기서 f는 물리적 플럭스, f_EC는 엔트로피 보존 두 점 플럭스, M은 질량 행렬, D는 미분 행렬, B는 경계 행렬이다. 두 가지 전환 전략이 제시된다. 1) **강인성 중심 전략**: WF가 FD보다 더 많은 엔트로피를 생성하면(즉, ˙S_WF > ˙S_FD) 해당 요소를 FD 형태로 전환한다. 이는 엔트로피가 과다 생산되는 상황을 방지해 스킴을 엔트로피 안정적으로 만든다. 2) **효율성 중심 전략**: 엔트로피 증가가 사전에 정의된 허용 한계 ε 이하이면 WF를 유지한다. 이 경우, 엔트로피가 약간 증가하더라도 전체 시뮬레이션 비용을 크게 절감할 수 있다. 이론적 검증으로는 다음이 수행된다. - **정확도**: 두 형태 모두 동일한 차수(p)의 다항식 근사 공간을 사용하므로, 전환이 일어나도 전역 수렴 차수는 변하지 않는다. 수치 실험에서 L2 오차는 기대 차수와 일치한다. - **선형 안정성**: FD 형태는 SBP(Summation‑By‑Parts) 속성을 만족해 에너지(엔트로피) 경계 항을 정확히 보존한다. WF는 전환 기준에 따라 엔트로피 과다 생산을 억제하도록 설계돼, 전체 스킴은 선형 안정성을 유지한다. - **엔트로피 적합성**: 제안된 지표는 엔트로피 불안정성을 정량적으로 감지하고, 필요 시 FD로 전환해 엔트로피 불안정성을 방지한다. 수치 실험에서는 다양한 1D·2D·3D 압축성 흐름 문제에 적용하였다. - **1D 충격파 문제**: 급격한 압력·밀도 불연속이 발생하는 영역에서 FD가 자동으로 활성화되어 비물리적 진동을 억제한다. - **2D 복잡 난류 케이스(예: 라인즈톤 폭발, 켈리-헬름홀츠 난류)**: 고차(p≥5)와 비텐서곱형 메쉬에서도 FD 전환이 국소화되어 전체 플럭스 평가 횟수가 30~50% 감소한다. - **3D 압축성 Navier‑Stokes‑Fourier 흐름**: 고온 고압 플라즈마 시뮬레이션에서, 엔트로피 기반 전환이 고온 코어 영역에만 적용돼, 외부 평탄 영역은 저비용 WF를 유지한다. 성능 측면에서, FD 형태는 분석 플럭스 대비 약 1.4~1.6배 더 비싸다(표 1). 그러나 v‑adaptivity를 적용하면 전체 시뮬레이션에서 FD가 차지하는 비율이 10~20% 수준으로 감소한다. 결과적으로 전체 CPU 시간은 균일 FD 스킴 대비 2배 이상 절감되었다. 또한, 저자는 기존 MOOD(다중 차원 최적 차수 검출)와 같은 사후 검증 적응 기법과 비교했을 때, v‑adaptivity가 사전 지표(엔트로피 생산량)를 사용해 전진 단계에서 바로 전환을 수행하므로, 재계산 비용이 크게 감소한다는 장점을 강조한다. 결론적으로, 본 연구는 DG 스킴에서 볼륨 항 선택을 동적으로 최적화함으로써 고차 정확도와 계산 효율 사이의 전통적인 트레이드오프를 크게 완화한다. 제안된 엔트로피 기반 지표는 다양한 물리 모델(압축성 Euler, Navier‑Stokes‑Fourier)과 메쉬 유형(텐서곱·비텐서곱)에서 적용 가능하며, 향후 고차 보존법 전반에 걸친 적응 전략 개발에 중요한 토대를 제공한다.