[논문 리뷰] High order well-balanced finite volume methods for multi-dimensional systems of hyperbolic balance laws

Jonas P. Berberich, Praveen Chandrashekar외 1명|arXiv (Cornell University)|2019. 03. 12.

Computational Fluid Dynamics and Aerodynamics인용 수 4

한 줄 요약

이 논문은 다차원 초월적 균형 법칙에 대해 일반적인 고차 정확도를 갖는 균형 잡힌 유한 체적 프레임워크를 제시한다. 이 방법은 알려진 정적 또는 시간에 따라 변화하는 해(해석적 또는 이산적)를 정확히 유지하기 위해 목표 해에서의 편차를 진화시킴으로써 이를 실현한다. 이 방법은 고차 정확도를 유지하며, 코드 수정이 최소화되고, 고차 수치법을 사용할 경우에도 계산 비용이 단지 5–30% 정도 증가할 뿐이다.

ABSTRACT

We introduce a general framework for the construction of well-balanced finite volume methods for hyperbolic balance laws. We use the phrase well-balancing in a broader sense, since our proposed method can be applied to exactly follow any solution of any system of hyperbolic balance laws in multiple spatial dimensions and not only time independent solutions. The solution has to be known a priori, either as an analytical expression or as discrete data. The proposed framework modifies the standard finite volume approach such that the well-balancing property is obtained and in case the method is high order accurate, this is maintained under our modification. We present numerical tests for the compressible Euler equations with and without gravity source term and with different equations of state, and for the equations of compressible ideal magnetohydrodynamics.

연구 동기 및 목표

• 다양한 다차원 초월적 균형 법칙 시스템에 적용 가능한 일반적인 고차 정확도를 갖는 균형 잡힌 유한 체적 방법을 개발하는 것.
• 특정 상태 방정식이나 수압 균형 가정에 제한되지 않고, 정적, 정적, 또는 시간에 따라 변화하는 임의의 해를 정확히 유지할 수 있도록 하는 것.
• 기존의 유한 체적 솔버와의 호환성을 확보하기 위해 재구성 전략과 소스 항 이산화 방법만 수정하면서도 고차 정확도를 유지하는 것.
• 압축성 에일러 방정식에 중력이 작용하는 복잡한 시스템과 이상 자기유체역학(MHD)에 대해 안정성과 효율성을 입증하는 것. 이는 비제로 속도를 갖는 비정상적인 정적 상태를 포함한다.
• 다양한 정확도 차수와 해 유형에 대해 균형 잡힘 수정으로 인한 계산 비용 증가를 정량화하는 것.

제안 방법

• 이 방법은 해 자체가 아니라 알려진 목표 해에서의 편차를 진화시킴으로써 목표 해에서 정확히 유지되도록 유한 체적 스킴을 재구성한다.
• 수압 재구성 원리를 임의의 목표 해에 적용하여, 알려진 해의 구조에 기반한 재구성 전략을 사용한다.
• 소스 항은 이산 수준에서 유량과 소스 항이 정확히 균형을 이루도록 보장하는 균형 잡힌 적분 전략을 통해 수정된다.
• 이 프레임워크는 임의의 일致된 재구성, 수치적 유량, 적분 규칙, 시간 적분기와 호환되며, 임의의 고차 정확도를 달성할 수 있다.
• 이 방법은 정렬된 카르테시안 및 극좌표 격자에 적용 가능하며, 정적 또는 시간에 따라 변화하는 목표 해를 모두 지원한다.
• 이산 목표 해(예: 수치 시뮬레이션에서 유도된 해)의 경우, 정확도와 효율성을 유지하기 위해 사전 계산된 데이터 배열을 사용한다.

실험 결과

연구 질문

• RQ1 일반적인 유한 체적 프레임워크를 설계하여 초월적 균형 법칙의 임의의 알려진 해를 시스템이나 해 유형에 관계없이 정확히 유지할 수 있는가?
• RQ2 해가 단순한 정적 상태가 아닌 경우에도 고차 정확도를 유지할 수 있는 균형 잡힌 스킴를 어떻게 설계할 수 있는가?
• RQ3 시간에 따라 변화하는 해 또는 복잡한 정적 해를 균형 잡기 위해 계산 오버헤드는 얼마나 되는가?
• RQ4 복잡한 상태 방정식이나 비정상적인 속도장을 갖는 시스템, 예를 들어 비대칭으로 돌고 있는 항성과 같은 시스템에 이 방법을 적용할 수 있는가?
• RQ5 균형 잡힘 수정의 계산 비용은 정확도 차수 증가에 따라 어떻게 변화하는가?

주요 결과

• 압축성 에일러 방정식에 중력이 작용하는 경우와 이상 자기유체역학(MHD)의 정적 및 정적 해를 복잡한 상태 방정식이나 비제로 속도장을 포함하여 정확히 유지한다.
• 일차원 정적 해의 경우, 일차 및 이차 수준의 방법에서 계산 비용이 약 20% 증가하였고, 칠차 수준에서는 단지 5% 증가하였다.
• 극좌표 격자에서 케플러 운동 디스크(정적, 시간에 따라 변화하지 않는 해)를 고려한 경우, 일차 수준의 경우 비용 증가가 20% 미만, 이차 수준의 경우 10% 미만이었다.
• 시간에 따라 변화하는 해(편미분이 있는 에일러 파동)의 경우, 일차 수준에서 벽 시계 시간 비율이 최대 30%에 달했지만, 삼차 수준에서는 약 15%에 그쳤다.
• 모든 수치 실험에서 고차 정확도가 유지되었으며, 관측된 수렴률은 이론적 기대와 일치하였다.
• 기존의 유한 체적 코드에 최소한의 수정으로 쉽게 통합 가능하며, 주로 소스 항 및 재구성 전략 조정이 필요하다.