High-order adaptive discontinuous finite elements for the shallow water equations with sub-grid irregular bathymetry

기술 요약: 하부 격자 불규칙 해저 지형을 가진 얕은 물 방정식을 위한 고차 적응 불연속 유한 요소

문제 제기
얕은 물 방정식 (SWE) 을 이용한 연안 유동의 수치 모델링은 이용 가능한 해저 지형 데이터의 높은 해상도와 고차 유한 요소 방법에서 일반적으로 사용되는 더 거친 격자 사이의 불일치라는 중대한 도전에 직면해 있습니다. 표준 고차 접근법은 종종 해저 지형을 해 변수와 동일한 차수의 다항식으로 근사화해야 합니다. 불규칙한 해저 지형의 경우, 이는 깁스 현상을 피하기 위해 평활화 또는 총변동 감쇠 (TVD) 제한이 필요하며, 이는 귀중한 고해상도 데이터를 폐기하게 됩니다. furthermore, 이러한 격자에서 젖음과 마름 (홍수) 을 처리하면서도 질량 보존, 수심 양수성, 그리고 "호수 정지" 정상 상태 (C-속성) 를 유지하는 것은 수학적으로 복잡하며, 특히 격자가 해저 지형 기울기나 젖음 - 마름 전선과 정렬되지 않을 때 더욱 그렇습니다.

방법론
저자들은 정규성 가정을 하지 않고 실제 불규칙 해저 지형 데이터를 활용하도록 설계된 deal.II 라이브러리 내에서 구현된 고차 불연속 갤러킨 (DG) 방법을 제안합니다. 핵심 방법론은 다음과 같습니다:

• 예측 변수: 연속 방정식의 주요 예측 변수로 수심 ($h$) 대신 자유 수면 높이 ($\zeta$) 를 선택합니다. 이는 고차 유한 요소에 적합한 매끄러운 표현을 용이하게 합니다. 수심은 $h = \max(\zeta + z_b, 0)$를 통해 적분점에서 비선형적으로 복원되며, 여기서 $z_b$는 참조 데이터에서 직접 샘플링된 해저 지형입니다.
• 하부 격자 해저 지형 표현: 해저 지형은 다항식 근사나 후처리 없이 적분점에서 직접 평가됩니다. 이를 통해 거친 격자에서도 정확한 하부 격자 해저 지형 특징을 복원할 수 있습니다.
• 이산화 체계:
  • 공간: 노드 기반 DG 방법이 사용됩니다. 정수압 압력 구배의 비보존적 형태는 [67] 에서 영감을 받은 특정 공식과 요소 인터페이스에서의 정수압 재구성을 결합하여 이산화됩니다. 이는 모든 해저 지형 함수에 대해 체계가 균형을 이루도록 보장합니다.
  • 시간: 암시적 - 명시적 (IMEX) Runge-Kutta (RK) 체계가 사용됩니다. 마찰 항은 얕은 영역의 강성을 처리하기 위해 암시적 (TR-BDF2 사용) 으로 처리되는 반면, 대류와 압력 구배는 명시적 (3 단계 2 차 RK 체계 사용) 으로 처리됩니다.
  • 젖음과 마름: 차수 적응 전략이 사용됩니다. 수심이 0 이거나 0 에 가까운 요소 (젖음 - 마름 전선) 에서 다항식 차수가 국소적으로 0 ($r=0$) 으로 감소되어, 사실상 1 차 셀 중심 유한 체적 근사로 전환됩니다. 이는 해저 지형의 임의의 제한 없이 안정성과 양수성을 보장합니다.
  • 추적자 수송: 수동 추적자를 위한 일관된 보존 이산화가 제안됩니다. 깊이 평균 농도를 진화시키고 연속 방정식과 추적자 방정식의 이산 형태가 일관되도록 함으로써, 불규칙한 해저 지형 시나리오에서 종종 발생하는 인위적인 단조성 위반을 방지합니다.
• 적응 격자 세분화 (AMR): 이 방법은 동적 $h$- 및 $p$-적응을 위한 deal.II 의 기능을 활용합니다. 격자 세분화는 오차 지표 (예: 와도 또는 추적자 기울기) 에 의해 주도되는 반면, 다항식 차수 적응은 젖음/마름 영역을 처리합니다.

주요 기여

1. 고차 하부 격자 방법: 이 논문은 해저 지형의 정규성을 가정하거나 해저 지형 데이터의 평활화/제한을 요구하지 않으면서 SWE 에 대해 균형 유지, 질량 보존, 양수성 보존을 동시에 결합하는 최초의 고차 방법을 제시합니다.
2. 수학적 성질: 저자들은 제안된 체계가 다음과 같음을 증명합니다:
   • 균형 유지: 임의의 해저 지형, 불연속적인 것을 포함하여 호수 정지 상태를 정확히 보존합니다.
   • 양수성 보존: 적분점 해상도에 의존하는 온화한 CFL 조건 하에서 수심은 음수가 되지 않습니다.
   • 일관성: 추적자 수송의 이산화는 연속 방정식과 일관성 (CWC 속성) 을 가지며, 단조성 위반을 방지합니다.
3. 불규칙성에 대한 강건성: 이 방법은 불규칙한 해저 지형과 격자 규모에서 과소 해석된 특징이 존재할 때 강건함을 보여주며, 다항식 기반 해저 지형 근사에서 흔히 발생하는 진동과 불안정성을 피합니다.

결과
이 방법은 일련의 이상화된 및 실제 벤치마크를 통해 검증되었습니다:

• 호수 정지: 이 체계는 불연속 해저 지형을 가진 직교 및 왜곡된 격자 모두에서 젖음 - 마름 인터페이스가 있더라도 기계 정밀도까지 정수압 평형을 보존했습니다.
• 수로 유동: 해석되지 않은 해저 지형 장애물이 있는 테스트에서 제안된 하부 격자 접근법은 정확한 정상 상태 해를 산출했습니다. 반면, 해저 지형의 다항식 근사는 인위적인 진동을 보이거나 올바른 정상 상태로 수렴하지 못했습니다.
• 이동하는 젖음 - 마름 전선 (Thacker 테스트): 이 방법은 젖음 - 마름 전환을 가진 주기적 진동을 올바르게 처리하여 수치 반올림 오차 내에서 질량 보존과 추적자 일관성을 유지했습니다.
• 이동하는 와류: 적응 격자 세분화 (AMR) 시뮬레이션은 정적 고해상도 격자와 동일한 정확도를 달성했으나, 자유도 (최대 9 배 감소) 와 계산 시간 (5~7 배 빠름) 은 크게 줄였습니다.
• 베니스 석호 벤치마크: 베니스 석호의 조석 순환에 대한 실제 시뮬레이션은 이 방법이 복잡한 기하학과 하부 격자 수로를 해결할 수 있음을 보여주었습니다. 하부 격자 접근법은 제한된 조각 선형 해저 지형 근사에 비해 더 강한 유속과 더 정확한 홍수 패턴을 생성했습니다. 동적 AMR 은 작은 질량 보존 오차 ($\sim 10^{-5}$) 를 도입했지만, 이는 일반적인 경계 데이터 불확실성에 비해 수용 가능한 것으로 판단되었습니다.

의의 및 주장
이 논문은 이 작업이 고해상도 해저 지형 데이터 세트를 완전히 활용하는 강건하고 정확한 연안 모델링 프레임워크를 제공한다고 주장합니다. 해저 지형 표현을 해의 다항식 기저와 분리함으로써, 이 방법은 전통적인 고차 체계에 내재된 정확성과 안정성 사이의 트레이드오프를 피합니다. 저자들은 그들의 접근 방식이 데이터 정규화 없이 불규칙한 해저 지형과 젖음/마름 전선을 가진 복잡한 연안 환경의 시뮬레이션을 가능하게 한다고 강조합니다. deal.II 와의 통합은 효율적인 병렬 처리와 동적 적응을 가능하게 하여, 이 방법을 고품질이고 계산 효율적인 연안 유동 시뮬레이션에 적합하게 만듭니다. 이 작업은 전파 속도에 기반한 시간 단계 제한을 완화하기 위한 압력 구배의 암시적 처리와 같은 향후 개발의 기초 역할을 합니다.