[논문 리뷰] Bound-preserving and entropy-stable algebraic flux correction schemes for the shallow water equations with topography

Hennes Hajduk, Dmitri Kuzmin|arXiv (Cornell University)|2022. 07. 15.

Advanced Numerical Methods in Computational Mathematics인용 수 5

한 줄 요약

이 논문은 고도를 고려한 얕은 수면 방정식에 대해 경계를 유지하고 엔트로피 안정적인 대칭적 유량 보정(AFC) 방법을 제안하며, 정지 상태에서의 균형, 수심의 비음성, 엔트로피 안정성을 확보하기 위해 단일화된 볼록 제한 전략을 사용한다. 이 방법은 기계 정밀도 수준의 정착 상태에서 강수량을 유지하고, 마찰 기반 알고리즘을 통해 강건한 젖고 마른 영역 시뮬레이션을 실현하며, 표준 CFL 조건 하에서 높은 정확도와 안정성으로 1차원 기준 시험을 통과하였다.

ABSTRACT

A well-designed numerical method for the shallow water equations (SWE) should ensure well-balancedness, nonnegativity of water heights, and entropy stability. For a continuous finite element discretization of a nonlinear hyperbolic system without source terms, positivity preservation and entropy stability can be enforced using the framework of algebraic flux correction (AFC). In this work, we develop a well-balanced AFC scheme for the SWE system including a topography source term. Our method preserves the lake at rest equilibrium up to machine precision. The low-order version represents a generalization of existing finite volume approaches to the finite element setting. The high-order extension is equipped with a property-preserving flux limiter. Nonnegativity of water heights is guaranteed under a standard CFL condition. Moreover, the flux-corrected space discretization satisfies a semi-discrete entropy inequality. New algorithms are proposed for realistic simulation of wetting and drying processes. Numerical examples for well-known benchmarks are presented to evaluate the performance of the scheme.

연구 동기 및 목표

• 수심의 비음성 및 엔트로피 안정성과 같은 핵심 물리적 성질을 유지하는 수치적 방법을 얕은 수면 방정식에 대해 개발하는 것.
• 소스 항이 있는 비균일 초월 방정식에 대해 대칭적 유량 보정(AFC) 방법을 확장하여, 특히 유한요소 해석에 적용하는 것.
• 비정렬 메esh에서 정착 상태 균형에 대한 균형 잡힌 성질을 확보하는 것 — 이는 유한요소 방법에서 오랫동안 해결되지 않은 과제이다.
• 불필요한 진동을 피하고 건조 영역에서의 안정성을 유지하는 강건한 젖고 마른 영역 알고리즘을 설계하는 것.
• 기존의 유한체적 접근법을 연속 유한요소 방법으로 일반화하면서 경계 및 엔트로피 성질을 유지하는 것.

제안 방법

• 기본 해석으로 연속적인 조각별 선형 유한요소를 사용하며, 성질 유지 유량 제한기를 통해 고차원으로 확장한다.
• 수심과 속도의 경계를 강제하기 위해 단일화된 볼록 제한(MCL) 전략을 적용하여 표준 CFL 조건 하에서 비음성 보장을 한다.
• 유량 보정된 공간 해석은 半-연속 엔트로피 부등식을 만족하여 엔트로피 안정성을 보장한다.
• 정착 상태에서의 균형 잡힘은 고도 소스 항을 정확히 상쇄시키는 수치적 유량 수정을 통해 달성되며, 기계 정밀도 수준까지 균형 상태를 유지한다.
• 새로운 마찰 기반 젖고 마른 영역 알고리즘이 도입되어 인위적 점성 또는 복잡한 재구성 기법에 의존하지 않는다.
• CFL 매개수 ν=0.05를 사용한 적응형 시간 스텝은 추가적인 확산이나 단계 반복 없이 안정성을 확보한다.

실험 결과

연구 질문

• RQ1 소스 항이 있는 얕은 수면 방정식에 대해 대칭적 유량 보정(AFC)을 확장하면서 수심의 비음성 및 엔트로피 안정성을 유지할 수 있는가?
• RQ2 소스 항이 있는 연속 유한요소 프레임워크에서 정착 상태 균형을 어떻게 달성할 수 있는가?
• RQ3 불필요한 진동을 피하고 건조 영역에서 안정성을 유지하는 강건한 젖고 마른 영역 알고리즘을 설계할 수 있는가?
• RQ4 제안된 AFC 방법이 이동하는 파면과 건조 상태를 포함한 기준 문제에서 정확성과 충격 포착 능력을 유지하는가?
• RQ5 이 방법을 고차원 유한요소로 일반화하고 코리olis 힘 또는 저면 마찰과 같은 추가 소스 항을 포함할 수 있는가?

주요 결과

• 제안된 AFC 방법은 기계 정밀도 수준까지 정착 상태 균형을 유지하며, 비정렬 메쉬에서 정확한 균형 잡힘을 보여준다.
• 표준 CFL 조건 하에서 수심의 비음성이 보장되며, 인위적 점성 또는 유량 보정 프레임워크 외의 추가 제한 없이도 가능하다.
• 유량 보정된 공간 해석은 半-연속 엔트로피 부등식을 만족하여 엔트로피 안정성을 보장한다.
• 댐 붕괴 문제의 수치 결과는 기준 솔루션과 뛰어난 일치를 보이며, 속도 피크가 이전 방법보다 약간 낮게 나타나 향후 정밀도 향상 가능성을 시사한다.
• 진동하는 포물선 형태의 호수 시험에서 MCL 및 MCL-SDE 방법은 정확한 주기적 솔루션을 잘 추적하며, 저차원 프로파일은 정확한 솔루션을 약간 뒤따른다.
• 마찰 기반 젖고 마른 영역 알고리즘은 ν=0.05에서 안정적이고 효과적이며, 시뮬레이션 붕괴나 불안정성을 유발하는 대체 접근법보다 뛰어난 성능을 보였다.