[논문 리뷰] A two-dimensional high-order well-balanced scheme for the shallow water equations with topography and Manning friction

Victor Michel-Dansac, Christophe Berthon외 2명|arXiv (Cornell University)|2020. 04. 09.

Computational Fluid Dynamics and Aerodynamics참고 문헌 77인용 수 7

한 줄 요약

이 논문은 지형과 맨니지 마찰을 고려한 얕은 수면 방정식에 대해 이차원 고차수 잘 균형 잡힌 유한체적 스킴을 제안한다. 고차수 다항식 재구성과 잘 균형 잡힌 일차 스킴을 안정 상태 감지기로 제어하는 볼록 조합을 통해 구현되며, 이는 매끄러운 영역에서는 고차수 정확도를 유지하면서도 이동하는 안정 상태를 정확히 유지하고, MOOD 절차를 통해 수치적 안정성을 확보한다. 실제 수심 데이터를 기반으로 한 2011년 일본 쓰나미의 정확한 시뮬레이션을 통해 검증되었다.

ABSTRACT

We develop a two-dimensional high-order numerical scheme that exactly preserves and captures the moving steady states of the shallow water equations with topography or Manning friction. The high-order accuracy relies on a suitable polynomial reconstruction, while the well-balancedness property is based on the first-order scheme from Michel-Dansac et. al., 2016 and Michel-Dansac et. al., 2017, extended to two space dimensions. To get both properties, we use a convex combination between the high-order scheme and the first-order well-balanced scheme. By adequately choosing the convex combination parameter following a very simple steady state detector, we ensure that the resulting scheme is both high-order accurate and well-balanced. The method is then supplemented with a MOOD procedure to eliminate the spurious oscillations coming from the high-order polynomial reconstruction and to guarantee the physical admissibility of the solution. Numerical experiments show that the scheme indeed possesses the claimed properties. The simulation of the 2011 Japan tsunami, on real data, further confirms the relevance of this technique.

연구 동기 및 목표

• 지형과 맨니지 마찰이 존재하는 이차원 얕은 수면 유동에서 이동하는 안정 상태를 정확히 유지하는 고차수 수치 스킴을 개발하는 것.
• 특히 비트레비얼한 안정 상태에 대해서도 잘 균형 잡힘 상태를 유지하면서 매끄러운 영역에서 고차수 정확도를 확보하는 것.
• 불필요한 진동과 수심의 음수 값을 방지함으로써 수치적 안정성을 확보하는 것.
• 실세계 데이터, 특히 2011년 일본 쓰나미를 포함하여 실제 수심 데이터에서의 검증을 통해 복잡한 지형에서의 정확성과 안정성을 입증하는 것.

제안 방법

• 안정 상태 감지기를 통해 동적으로 조정되는 조합 계수를 가진 고차수 재구성과 일차 수준의 잘 균형 잡힌 스킴 간의 볼록 조합을 사용한다.
• 매끄러운 영역에서 제2 및 제3차 정확도를 달성하기 위해 고차수 다항식 재구성을 적용한다.
• 고차수 재구성에서 발생하는 불필요한 진동을 탐지하고 제거하기 위해 MOOD(Multi-dimensional Optimal Order Detection) 절차를 적용한다.
• 지형과 맨니지 마찰 항을 잘 균형 잡힌 프레임워크 내에서 정확하게 처리하는 고두노프 유형의 유한체적 프레임워크를 구현한다.
• 지역 유동 조건에 따라 고차수 모드와 잘 균형 잡힌 모드를 전환하기 위한 새로운 안정 상태 감지기를 도입한다.
• Rusanov 유형의 수치적 플럭스와 소스항 균형 조정을 사용한 유한체적 이산화를 구현하여 잘 균형 잡힘 상태를 유지한다.

실험 결과

연구 질문

• RQ1 지형과 맨니지 마찰이 모두 존재하는 2차원 얕은 수면 유동에서 이동하는 안정 상태를 정확히 유지할 수 있는 고차수 유한체적 스킴을 설계할 수 있는가?
• RQ2 비트레비얼한 안정 상태, 특히 마찰이 존재하는 경우에도 잘 균형 잡힘 상태를 유지하면서 고차수 정확도를 유지하는 방법은 무엇인가?
• RQ3 고차수 재구성은 수치적 진동에 어떤 영향을 미치며, 정확도나 안정성을 희생시키지 않고 이를 제어하는 방법은 무엇인가?
• RQ4 제안된 스킴은 2011년 일본 쓰나미와 같이 극한의 지형 기울기를 가진 실세계 수심 데이터에서 어떻게 성능을 발휘하는가?
• RQ5 복잡한 메esh에서 고차수 스킴의 계산 비용은 실용적인 쓰나미 시뮬레이션에서 향상된 정확성으로 인해 정당화되는가?

주요 결과

• PWB1 스킴(제2차 수준 잘 균형 잡힘)은 수치적 확산이 현저히 감소하고 파동 구조의 해상도가 향상되어 P0 스킴(일차 수준)과 유사한 결과를 도출하였다.
• PWB3 스킴(제3차 수준 잘 균형 잡힘)은 극한의 지형 기울기로 인해 과도한 MOOD 안정화가 발생하여 기대한 정확도를 확보하지 못했으며, 결과적으로 일차 수준 스킴과 유사한 해를 보였다.
• PWB1 스킴은 정확한 쓰나미 전파 시간을 잘 캡처했고, 불필요한 진동이 없었으며, 효과적인 잘 균형 잡힘 상태와 안정성을 입증하였다.
• P0 및 PWB1 스킴 모두 물리적 센서 데이터와 잘 일치했으며, 특히 파동의 도착 시간과 진폭에서 PWB1 스킴이 P0 스킴보다 더 우수한 근사치를 제공하였다.
• 약 1,300만 개의 격자 요소를 사용한 실제 수심 데이터 시뮬레이션에서, 스킴은 특히 키테이르 해구 근처에서 큰 지형 기울기를 정확히 처리하는 데 성공했다.
• PWB1의 계산 비용은 P0 대비 약 2.5배 높았고, PWB3는 약 10배 높아, 복잡한 메쉬에서 정확도와 효율성 사이의 상충 관계를 보였다.