샬로우 워터 모멘트 방정식의 엔트로피 보존 및 엔트로피 안정 DG 방법

초록

본 논문은 샬로우 워터 모멘트 방정식(SWME)이 전체 에너지에 해당하는 엔트로피 보존 법칙을 만족함을 증명하고, 뉴턴식 슬립 마찰과 Manning 마찰이 엔트로피 소산 항임을 보인다. 이를 바탕으로 엔트로피 보존 수치 플럭스를 설계하고, 엔트로피 안정성과 웰밸런스를 갖는 노드형 DG 스펙트럴 요소법을 제안한다. 수치 실험을 통해 이론적 결과와 알고리즘의 효율성을 확인한다.

상세 분석

본 연구는 1차원 샬로우 워터 모멘트 방정식(SWME)을 대상으로 연속적인 엔트로피 구조를 체계적으로 분석한다. 먼저, 수직 좌표에 대한 다항식 전개를 이용해 평균 속도와 N개의 모멘트 α_i를 정의하고, 이들에 대한 보존형식(3a‑3c)을 제시한다. 기존 SWE와 달리 모멘트 간 비선형 결합항 A_i와 B_i가 존재하는데, 저자들은 이 항들을 정밀히 다루어 전체 에너지 식을 도출한다. Lemma 1‑3을 통해 평균 운동에너지, 포텐셜 에너지, 그리고 모멘트에 의해 발생하는 추가 운동에너지를 각각 보존형식으로 변환하고, 최종적으로 총 에너지 E와 엔트로피 플럭스 F를 (20)‑(21) 형태로 정의한다. 여기서 E는 h·u_m²/2 + Σ α_i²/(2i+1) + gh²/2 + ghb 로, 볼록성은 헤시안 행렬이 양정인 것으로 확인된다. 엔트로피 변수 w는 E를 보존 변수 u에 대해 미분한 결과이며, w₁ = −u_m²/2 − Σ α_i²/(2i+1) + g(h+b), w₂ = u_m, w_{i+2}=α_i/(2i+1) 로 얻어진다.

마찰 항에 대해서는 두 가지 대표 모델을 고려한다. (24)식의 뉴턴식 슬립 마찰은 바닥 마찰과 내부 점성 마찰을 포함하고, (28)식의 Manning 마찰은 물 깊이에 대한 비선형 저항을 반영한다. Lemma 4와 Lemma 5에서 각각 wᵀS ≤ 0 를 증명함으로써, 이들 마찰이 엔트로피를 소산시키는 항임을 보인다. 이는 물리적 에너지 손실과 수치적 안정성 사이의 일관성을 확보한다.

수치적 구현에서는 노드형 DG 스펙트럴 요소를 채택한다. 엔트로피 보존 플럭스는 Tadmor의 엔트로피 호환 조건을 만족하도록 설계되었으며, 인터페이스에서 적절한 수치 점성(예: 로컬 Lax‑Friedrichs) 을 추가해 엔트로피 안정성을 확보한다. 또한, 바닥 높이 b(x)와 정적 평형 상태를 정확히 보존하도록 웰밸런스 기법을 도입하였다. 이러한 설계는 SBP‑like 성질을 갖는 차분 연산자와 결합되어 반사 경계와 급격한 수위 변동에서도 비진동적인 해를 제공한다.

마지막으로, 1차 및 2차 모멘트 모델에 대한 다양한 테스트 케이스(정적 호수, 급격한 수위 파동, 마찰에 의한 감쇠)를 수행한다. 결과는 엔트로피 감소가 기대대로 관측되며, 기존 비안정적 DG 스키마에 비해 수치 진동이 현저히 감소함을 보여준다. 특히, Manning 마찰이 포함된 경우에도 엔트로피 감소가 유지되어 물리적 에너지 손실을 정확히 포착한다. 전체적으로, 연속 엔트로피 분석을 기반으로 한 DG 설계가 SWME의 복잡한 비선형 구조와 마찰 효과를 동시에 다룰 수 있음을 입증한다.