[논문 리뷰] High-order schemes for the Euler equations in cylindrical/spherical coordinates

Sheng Wang, Eric Johnsen|arXiv (Cornell University)|2017. 01. 17.

Computational Fluid Dynamics and Aerodynamics참고 문헌 23인용 수 8

한 줄 요약

이 논문은 원통좌표계 및 구좌표계에서의 오일러 방정식에 대해 고차 정확도와 보존성이라는 두 성질을 동시에 만족하는 새로운 고차 유한차 WENO 스킴을 제안한다. 이전까지는 이 두 성질을 동시에 확보하는 것이 어려웠다. 새로운 방법은 수렴성 연구와 충격파관통 문제/Sedov 폭발 문제를 통해 검증되었으며, 원통좌표계나 구좌표계의 축방향만을 고려하는 경우 기존 방법들보다 정확도와 보존성을 동시에 유지함으로써 뛰어난 성능을 보였다.

ABSTRACT

We consider implementations of high-order finite difference Weighted Essentially Non-Oscillatory (WENO) schemes for the Euler equations in cylindrical and spherical coordinate systems with radial dependence only. The main concern of this work lies in ensuring both high-order accuracy and conservation. Three different spatial discretizations are assessed: one that is shown to be high-order accurate but not conservative, one conservative but not high-order accurate, and a new approach that is both high-order accurate and conservative. For cylindrical and spherical coordinates, we present convergence results for the advection equation and the Euler equations with an acoustics problem; we then use the Sod shock tube and the Sedov point-blast problems in cylindrical coordinates to verify our analysis and implementations.

연구 동기 및 목표

• 원통좌표계 및 구좌표계에서 오일러 방정식에 대해 고차 정확도와 보존성을 동시에 확보하는 데 오랫동안 해결되지 않은 과제를 해결하기 위해.
• 세 가지 공간 이산화 방법을 평가하고 비교하기 위해: 하나는 정확하지만 비보존적인 방법, 하나는 보존적이지만 고차 정확도가 아닌 방법, 그리고 새로운 방법은 두 조건을 모두 만족한다.
• 제안된 방법의 성능을 평활 문제(진동, 음향)와 충격이 지배하는 유동(소드 충격파관통, Sedov 점폭발)에 대해 검증하기 위해, 특히 축방향 기하구조에서 수행한다.
• 향후 이산 갈레르킨 방법 및 확산 효과로의 확장을 위한 기반을 마련하기 위해, 특히 캐비테이션 기포의 형성과 붕괴를 모델링하기 위해.

제안 방법

• 축방향 의존성만을 갖는 원통좌표계 및 구좌표계에서의 오일러 방정식을 유량차분 형태로 기술하며 기하학적 소스항을 포함한다.
• 세포 경계에서의 수치 유량을 재구성하기 위해 비선형 가중치를 사용하는 오차가 제5차 정확도인 유한차 WENO 스킴을 적용하여 불연속성의 고해상도 캡처를 보장한다.
• 로컬 라크-프리드리히스 유량 분할을 사용하여 수치 유량을 계산하며, 안정성과 정확도를 확보하기 위해 점별 고유값 기반 파동 속도 추정을 수행한다.
• 시간 정확도 확보를 위해 제4차 정확도의 명시적 룬게-쿠타 시간 적분 스킴을 구현한다.
• 새로운 공간 이산화 방법(제3의 방법)을 도입하여 유량항을 재구성함으로써 보존성을 유지하면서도, 축방향 메트릭 항의 신중한 처리를 통해 고차 정확도를 유지한다.
• 부드러움 지표와 비선형 가중치를 사용하여 충격 근처에서 본질적으로 진동하지 않는 행동을 유지하며, 분모가 0이 되는 것을 방지하기 위해 작은 매개변수 ε = 10⁻⁶를 사용한다.

실험 결과

연구 질문

• RQ1 원통좌표계 및 구좌표계에서 오일러 방정식에 대해 고차 정확도와 보존성을 동시에 확보하는 유한차 WENO 스킴을 구성할 수 있는가?
• RQ2 곡선좌표계에서 축방향 유동에 적용했을 때, 다양한 공간 이산화 방법들이 정확도와 보존성 측면에서 어떻게 비교되는가?
• RQ3 오차가 제5차 정확도인 고차 WENO 스킴이 평활 문제인 진동 방정식과 음향 방정식에 대해 원통좌표계 및 구좌표계에서 수렴 성질을 어떻게 보이는가?
• RQ4 제안된 스킴이 소드 충격파관통 문제 및 Sedov 점폭발 문제와 같은 벤치마크 문제에서 충격파와 불연속성을 얼마나 잘 캡처하는가?
• RQ5 강한 불연속성과 기하학적 소스항이 존재하는 문제에서, 새로운 방법이 정확도를 희생시키지 않고도 보존 성질을 유지할 수 있는가?

주요 결과

• 제안된 제3의 방법은 고차 정확도와 보존성을 동시에 확보하지만, 나머지 두 방법은 이 두 성질을 동시에 만족하지 못한다.
• 수렴성 연구 결과, 새로운 방법은 평활 문제, 즉 진동 방정식과 음향 문제에서 최적의 제5차 정확도를 보였다.
• 소드 충격파관통 문제를 통해 새로운 방법이 충격 위치를 정확히 캡처하고 질량, 동량, 에너지의 보존성을 유지함을 확인하였다.
• Sedov 점폭발 문제에서는 새로운 방법이 자가유사 폭발파 구조를 정확히 재현하며 총 에너지 오차를 최소화하여 보존성을 유지함을 보였다.
• 연구 결과, 곡선좌표계에서 표준 WENO 구현은 특별히 수정되지 않은 한 물리적 양의 보존에 실패하는 경우가 많음을 확인하였으며, 이는 제안된 재구성의 필요성을 강조한다.
• 결과적으로 새로운 방법이 충격과 계면을 포함한 구형 또는 원통형 대칭 유동을 시뮬레이션하는 데 적합하다는 것이 밝혀졌으며, 이는 관성적 압축융합 또는 캐비테이션 기포의 동역학 모델링에 유용할 것이다.