[논문 리뷰] A robust and stable numerical scheme for a depth-averaged Euler system

Nora Aïssiouene, Marie-Odile Bristeau외 2명|arXiv (Cornell University)|2015. 06. 10.

Computational Fluid Dynamics and Aerodynamics참고 문헌 38인용 수 5

한 줄 요약

이 논문은 자유 표면을 가진 압축성 불가약 유동 방정식에서 유도된 비수압 균형 깊이 평균 유동 방정식을 위한 안정적이고 강건한 유한체적 방법을 제시한다. 운동량 기반의 투영-보정 방법을 사용함으로써, 잘 균형 잡힌 상태, 양성 보장, 그리고 완전 이산화된 엔트로피 부등식을 확보하며, 건조 상태에서도 안정성을 유지한다. 분석 기준 문제인 탄성 파라볼로이드 우물과 단일파에 대해 수렴성과 정확성이 검증되었다.

ABSTRACT

We propose an efficient numerical scheme for the resolution of a non-hydrostatic Saint-Venant type model. The model is a shallow water type approximation of the incompressbile Euler system with free surface and slightly differs from the Green-Naghdi model. The numerical approximation relies on a kinetic interpretation of the model and a projection-correction type scheme. The hyperbolic part of the system is approximated using a kinetic based finite volume solver and the correction step implies to solve an elliptic problem involving the non-hydrostatic part of the pressure. We prove the numerical scheme satisfies properties such as positivity, well-balancing and a fully discrete entropy inequality. The numerical scheme is confronted with various time-dependent analytical solutions. Notice that the numerical procedure remains stable when the water depth tends to zero.

연구 동기 및 목표

• 고전적 샹티에 방정식을 확장한 비수압 균형 깊이 평균 에일러 모델을 위한 안정적이고 정확한 수치적 방법을 개발하는 것.
• 해법이 잘 균형 잡힘, 수심의 양성, 완전 이산화된 엔트로피 부등식을 유지하도록 보장하는 것.
• 특히 수심이 0에 수렴할 때 발생하는 습식/건식 전면을 강력하게 처리하는 것.
• 단일파 및 탄성 파라볼로이드 우물 구조를 포함한 시간에 의존하는 해석적 해와의 검증.
• 보정 단계의 변분형식을 통해 2차원 비정렬 메esh로의 확장 기반 마련

제안 방법

• 해법은 분할 스텝 투영-보정 전략을 사용하여 시스템을 초구형 및 타원형 하위단계로 분리한다.
• 초구형 부분은 운동론 기반의 유한체적 방법으로 해결하여 일致성과 안정성을 확보한다.
• 보정 단계는 비수압 압력에 대한 이산 타원형 문제를 유한차분 이산화를 통해 해결한다.
• 이 방법은 완전 이산화된 엔트로피 부등식을 만족하고 이산 인프-서프 조건을 유지함을 증명하였다.
• 경계 조건은 유량 강제 방식으로 처리되며, 특히 단일파 유입과 같은 시간에 의존하는 유입 시나리오에 적합하다.
• 깊이 평균 에일러 시스템의 운동론적 해석을 통해 소스 항과 비수압 효과의 강력한 수치적 처리가 가능하다.

실험 결과

연구 질문

• RQ1 비수압 균형 깊이 평균 에일러 모델을 위한 강력한 수치적 해법을 설계할 수 있는가? 이는 잘 균형 잡힘과 양성 등의 핵심 물리적 성질을 유지할 수 있는가?
• RQ2 비수압 압력에 관련된 보정 단계를 유한차분을 사용해 안정적이고 정확하게 이산화할 수 있는가?
• RQ3 특히 습식/건식 인터페이스에서 수심이 0에 수렴할 때 해법이 안정성과 수렴성을 유지하는가?
• RQ4 단일파 및 탄성 파라볼로이드 유동과 같은 시간에 의존하는 해석적 해에 대해 해법의 수렴 속도는 어떻게 되는가?
• RQ5 보정 단계의 변분형식을 활용해 2차원 비정렬 메쉬로의 확장을 가능하게 할 수 있는가?

주요 결과

• 수심이 0에 수렴할 때조차도 해법이 안정성과 정확성을 유지하여 습식/건식 전면의 강력한 처리가 가능하다.
• 탄성 파라볼로이드 우물 시험 케이스에서, 수심의 L¹ 오차는 시간이 지나도 유한하게 유지되며, 일阶 및 이阶 해법 모두 이론적 예측 수렴 속도를 충족한다.
• 단일파 시험 케이스에서, 내부 및 경계 유입 조건 모두에서 이론적 순서와 일치하는 수렴 속도를 보이며, 최대 3200개의 격자에서 오차가 계산되었다.
• 두 기준 문제 모두에서 공간 및 시간에 대해 이차 수렴을 달성하여 정확성과 강력성을 확인하였다.
• 완전 이산화된 해법에서 이산 엔트로피 부등식이 만족됨을 확인하여 해법의 열역학적 일致성을 입증하였다.
• 다양한 메쉬 해상도와 경계 조건 조건에서도 해법이 강력함을 입증하였으며, 특히 경계에서 솔리톤이 들어오는 경우에도 안정성을 유지하였다.