[논문 리뷰] An entropy-stable and fully well-balanced scheme for the Euler equations with gravity

Christophe Berthon, Victor Michel-Dansac외 1명|arXiv (Cornell University)|2024. 06. 21.

Navier-Stokes equation solutions인용 수 1

한 줄 요약

이 논문은 중력이 포함된 에일러 방정식에 대해 완전히 잘 균형 잡힌, 엔트로피 안정성과 양성 보존을 갖춘 신규한 고두노프 유형의 유한체적 스킴을 제시한다. 하르텐, 랙스, 반 룰의 HLL 리만 해법에 적분 일致성과 엔트로피 안정성 기반의 소스 항 처리를 확장함으로써, 이 스킴은 기계 정밀도 수준에서 모든 움직이는 상태와 정적 등온 평형 상태를 정확히 유지하면서 이산 엔트로피 부등식을 만족시키고 밀도와 압력을 양성으로 유지한다.

ABSTRACT

The present work concerns the derivation of a numerical scheme to approximate weak solutions of the Euler equations with a gravitational source term. The designed scheme is proved to be fully well-balanced since it is able to exactly preserve all moving equilibrium solutions, as well as the corresponding steady solutions at rest obtained when the velocity vanishes. Moreover, the proposed scheme is entropy-preserving since it satisfies all fully discrete entropy inequalities. In addition, in order to satisfy the required admissibility of the approximate solutions, the positivity of both approximate density and pressure is established. Several numerical experiments attest the relevance of the developed numerical method. An extension to two-dimensional problems is given, applying the one-dimensional framework direction by direction on Cartesian grids.

연구 동기 및 목표

• 중력이 포함된 에일러 방정식의 모든 등온 평형 해—움직이는 상태이든 정적 상태이든—을 정확히 유지하는 수치 스킴을 개발하기.
• 이 스킴이 이산 엔트로피 부등식을 만족시켜 약한 해의 물리적 적합성을 보장하기.
• 모든 이산 해에서 밀도와 압력의 엄격한 양성을 유지하기, 이는 견고성에 필수적인 조건이다.
• 일차 스킴을 고차 정확도로 확장하면서도 완전히 잘 균형 잡힌 성질을 유지하기.
• 이중 차원에서 복잡한 평형 상태와 중력 작용이 있는 동적 유동을 정확히 포착할 수 있는 능력을 입증하기.

제안 방법

• 소스 항에 대해 하르텐, 랙스, 반 룰의 HLL 적분 일치 조건을 만족하는 두 개의 중간 상태를 갖는 근사 리만 해법을 설계한다.
• 모든 평형 해, 비트리비얼한 움직이는 평형 상태까지 정확히 잘 균형 잡힌 해법이 되도록 소스 항 이산화를 통합한다.
• 엔트로피 적분 일치 조건을 강제 적용하여 이산 엔트로피 부등식이 만족됨을 보장하며, HLL 엔트로피의 감소 원리를 활용한다.
• 수정된 리만 해법을 기반으로 한 고두노프 유형의 유한체적 스킴을 구성하여 일致성, 양성 보존, 엔트로피 안정성을 확보한다.
• MUSCL 유형의 고차원 재구성 기법을 적용하여 일차 스킴을 고차 정확도로 업그레이드하면서도 완전히 잘 균형 잡힌 성질을 유지한다.
• 직각 좌표계에서 차원 분리 기법을 통해 일차원 프레임워크를 이중 차원으로 확장하고, 각 방향에 대해 순차적으로 잘 균형 잡힌 절차를 적용한다.

실험 결과

연구 질문

• RQ1 중력이 포함된 에일러 방정식의 모든 움직이는 상태와 정적 등온 평형 해를 정확히 유지하는 수치 스킴을 구성할 수 있는가?
• RQ2 중력 소스 항이 존재하는 상황에서 일치성과 양성을 유지하면서도 이산 엔트로피 안정성을 어떻게 확보할 수 있는가?
• RQ3 중력 하에서 잘 균형 잡힘과 엔트로피 안정성을 동시에 달성하기 위해 HLL 리만 해법에 어떤 수정이 필요한가?
• RQ4 일차 스킴을 고차 정확도로 확장할 수 있는 정도는 어느 정도이며, 이 과정에서 완전히 잘 균형 잡힌 성질을 잃지 않는가?
• RQ5 중력이 한 방향에서만 작용하는 이중 차원 평형 상태를 정확히 포착할 수 있는가?

주요 결과

• 제안된 스킴은 1차원 및 2차원 수치 시험을 통해 모든 등온 평형 해—움직이는 상태이든 정적 상태이든—을 기계 정밀도 수준에서 정확히 유지한다.
• 소스 항 처리를 수정한 HLL 기반 리만 해법은 엔트로피 감소 원리를 만족하여 이산 엔트로피 부등식 만족의 증명이 가능하다.
• 2차원 움직이는 평형 상태 시험에서 잘 균형 잡힌 스킴은 밀도와 압력에 대해 상대 L2 오차가 10−16 이하로 유지되었고, 잘 균형 잡히지 않은 HLL 스킴은 10−3에서 10−2 수준의 오차를 보였다.
• 2차원 움직이는 평형 상태의 미세한 변화 시험에서 비물리적 파동이나 잡음이 발생하지 않아, 소용돌이 상태를 작은 외란 하에서도 유지할 수 있음을 확인했다.
• 가우시안 중력 포텐셜을 가진 압축 시험에서 반지름 방향의 압축과 충격 형성 현상을 정확히 포착하였고, 비물리적 진동이 발생하지 않았다.
• 모든 시뮬레이션 동안 밀도와 압력의 양성 보존이 유지되었으며, 모든 테스트 케이스에서 엔트로피 부등식이 만족됨을 확인했다.