[논문 리뷰] A large time-step and well-balanced Lagrange-Projection type scheme for the shallow-water equations

Christophe Chalons, Pierre Kestener외 2명|arXiv (Cornell University)|2016. 04. 06.

Computational Fluid Dynamics and AerodynamicsEngineering참고 문헌 31인용 수 22

한 줄 요약

이 논문은 음파를 암시적으로, 물질파를 명시적으로 처리하는 음압-투과(implicit-explicit, IMEX) 접근 방식을 사용하여, 난류수방정식에 대한 큰 시간스텝, 잘 균형 잡힌 라그랑주-투과 방법을 제시한다. 이 방법은 음속에 기반한 CFL 제약 조건을 제거하여 안정적인 큰 시간스텝을 허용하며, 정체된 물의 상태를 유지하면서도 에너지 안정성과 수심의 양수성을 유지한다.

ABSTRACT

This work focuses on the numerical approximation of the Shallow Water Equations (SWE) using a Lagrange-Projection type approach. We propose to extend to this context recent implicit-explicit schemes developed in the framework of compressibleflows, with or without stiff source terms. These methods enable the use of time steps that are no longer constrained by the sound velocity thanks to an implicit treatment of the acoustic waves, and maintain accuracy in the subsonic regime thanks to an explicit treatment of the material waves. In the present setting, a particular attention will be also given to the discretization of the non-conservative terms in SWE and more specifically to the well-known well-balanced property. We prove that the proposed numerical strategy enjoys important non linear stability properties and we illustrate its behaviour past several relevant test cases.

연구 동기 및 목표

• 기존의 고두노프 유형의 수치해법에서 음파에 기반한 제한적인 CFL 조건을 극복하는 난류수방정식을 위한 수치해법을 개발하는 것.
• 정체된 물의 상태를 유지하는 잘 균형 잡힌 성질을 확보하여, 정적 수류 해법을 정확히 유지하는 것.
• 저속 및 초음속 흐름에서의 안정성을 확보하기 위해 수심의 양수성과 이산화된 엔트로피 부등식을 유지하는 것.
• 최근의 암시-명시적(Lagrange-Projection) 전략을 압축성 유동에서 난류수계로 확장하여, 정확성과 안정성을 희생시키지 않고 큰 시간스텝을 허용하는 것.
• 유량, 초음속, 비유일성 Riemann 해를 포함한 벤치마크 테스트 케이스에서 해법을 검증하는 것.

제안 방법

• 계산체계를 음파와 물질파 성분으로 분리하기 위해 라그랑주-투과 분해를 사용한다.
• 음파는 반암시적 시간 이산화를 통해 암시적으로 처리되어 음속에 의존하지 않는 시간스텝을 허용한다.
• 물질파는 명시적으로 처리되어 느린, 초음속 흐름에서의 정확성을 유지한다.
• 비균일 메esh에서도 정체된 물의 상태를 정확히 포착하기 위해 잘 균형 잡힌 재구성 기법을 사용한다.
• 이산 해법에서 엔트로피 안정성을 보장하기 위해 보존형 엔트로피 공식을 적용한다.
• 유한체적 프레임워크 내에서 라그랑주-투과 분해에서 유도된 수치적 플럭스를 사용하여 해법을 구현한다.

실험 결과

연구 질문

• RQ1 IMEX 라그랑주-투과 방법이 난류수 시뮬레이션에서 음속 기반의 안정성 또는 정확성 제약 없이 큰 시간스텝을 달성할 수 있는가?
• RQ2 제안된 해법이 지형과 비균일 메쉬 조건에서도 정체된 물의 상태를 유지하는 잘 균형 잡힌 성질을 보존하는가?
• RQ3 기존 방법들과 비교해 비유일 엔트로피 해를 갖는 Riemann 문제의 해를 어떻게 정확하게 포착하는가?
• RQ4 음파의 암시적 처리가 CFL 조건과 계산 효율성에 어떤 영향을 미치는가?
• RQ5 복잡한 흐름 조건에서 수심의 양수성과 이산화된 엔트로피 부등식을 유지할 수 있는가?

주요 결과

• 유량 영역 테스트 케이스에서, 기존의 명시적 CFL 한계보다 10배 큰 시간스텝을 사용하여 안정적인 시뮬레이션을 구현하였다.
• 라그랑주-투과 해법은 수압 재구성 기반 HRHLL 해법과 매우 유사한 해를 생성하였으며, ACU 해법보다 더 높은 정확도를 보였다.
• 충격파가 존재하는 초음속 영역에서, 정확한 충격 위치와 파동 구조를 잘 포착하였으며, 기존의 확립된 해법과 일치하는 결과를 보였다.
• 비유일 Riemann 문제에서, 수압 재구성 방법과 동일한 해로 수렴하였고, HLLACU는 다른, 덜 정확한 해를 생성하였다.
• 음파의 암시적 처리로 인해 음속 기반의 CFL 제약 조건이 제거되어 큰 시간스텝을 허용하면서도 잘 균형 잡힘과 양수성을 유지하였다.
• 모든 테스트 케이스에서 엔트로피 안정성과 수심의 양수성을 유지하였으며, 강한 지형 조건과 근접 정체 조건에서도 안정하였다.