[논문 리뷰] Admissible states and physical constraints preserving numerical schemes for special relativistic magnetohydrodynamics

Kailiang Wu, Huazhong Tang|arXiv (Cornell University)|2016. 03. 22.

Computational Fluid Dynamics and Aerodynamics참고 문헌 31인용 수 7

한 줄 요약

이 논문은 특수 상대론적 자기유체역학(RMHD)에서 물리적 제약을 유지하는(PCP) 수치적 스킴의 수학적 기초를 구축하며, 허용 상태 집합 𝒢의 명시적 등가 형태를 유도한다. 𝒢의 볼록성을 증명하고, 이산적 발산 없는 자기장 조건 하에서 일반화된 라크스-프리드리히스 분할 성질을 도입함으로써, CFL 조건 하에서 증명 가능한 물리적 일致성을 갖는 1D 및 2D RMHD에서의 최초의 고차 정확도 PCP 스킴을 가능하게 한다.

ABSTRACT

The paper studies mathematical properties of the admissible state set $\mathcal G$ of relativistic magnetohydrodynamics(RMHD), paving a way for developing numerical schemes for RMHDs with the solutions in $\mathcal G$, referred as physical-constraints-preserving (PCP) schemes. To overcome the difficulties arising from the extremely strong nonlinearities and no explicit formulas of the primitive variables and the flux vectors with respect to the conservative vector, two equivalent forms of $\cal G$ with explicit constraints on the conservative vector are skillfully discovered. The first form is derived by analyzing polynomial roots and transferring successively, and used to prove the convexity of $\mathcal G$ with the aid of semi-positive definiteness of the second fundamental form of a hypersurface. While the second form is derived based on the convexity, and followed by the orthogonal invariance of $\cal G$. The Lax-Friedrichs (LxF) splitting property does not hold generally for the nonzero magnetic field, but by a constructive inequality and pivotal techniques, we discover the generalized LxF splitting properties, combining the convex combination of some LxF splitting terms with a discrete divergence-free condition of the magnetic field. Based on the above theoretical analyses, several one- and two-dimensional PCP schemes are studied for the first time. In the 1D case, a first-order accurate LxF type scheme is first proved to be PCP under the CFL condition, and then the high-order accurate PCP schemes are proposed with the PCP limiter. In the 2D case, the discrete divergence-free condition and PCP property are analyzed for a first-order accurate LxF type scheme, and followed by two sufficient conditions for high-order accurate PCP schemes. The current analyses reveal in theory for the first time that the discrete divergence-free condition is closely connected with the PCP property.

연구 동기 및 목표

• 특수 상대론적 자기유체역학(RMHD)에서 수치적 해의 물리적 타당성을 보장하는 허용 상태 집합 𝒢의 수학적 구조를 확립하는 것.
• RMHD에서 원시 변수와 유량에 대한 명시적 공식이 부족하고 극도로 비선형적인 문제를 해결하기 위해, 보존 변수에 대한 등가 명시적 제약 조건을 도출하는 것.
• 초등형 표면의 제2 기본형의 준정부정성에 기반하여, 허용 상태 집합 𝒢의 볼록성을 증명하는 것.
• 이산적 발산 없는 자기장 조건을 통합한 일반화된 라크스-프리드리히스 분할 성질을 개발하여, 안정적이고 물리적으로 일致한 수치적 스킴을 가능하게 하는 것.
• CFL 조건 하에서 증명 가능한 안정성과 물리적 일치성을 갖는 1D 및 2D RMHD에 대한 최초의 고차 정확도 물리적 제약을 유지하는(PCP) 스킴을 구축하는 것.

제안 방법

• 다항식 근의 분석과 초등형 표면 기하학을 활용하여, 허용 상태 집합 𝒢의 두 가지 등가 명시적 형태를 유도하고, 제2 기본형의 준정부정성에 기반한 볼록성 증명을 수행한다.
• 𝒢의 볼록성을 활용하여, 수직 불변성을 갖는 두 번째 등가 형태를 도출함으로써, 구조적 명확성과 수치적 설계의 향상을 도모한다.
• 이산적 발산 없는 자기장 조건을 구성하고, LxF 항의 볼록 조합과 결합하여 일반화된 라크스-프리드리히스 분할을 도입한다.
• 이산적 발산 없는 조건이 PCP 성질과 본질적으로 연결되어 있음을 증명하여, 수치적 스킴에서의 안정성과 물리적 일치성을 보장한다.
• CFL 조건 하에서 증명 가능한 PCP 성질을 갖는 1D 및 2D의 1차 정확도 LxF 유형 스킴을 설계하고, PCP 리미터를 활용하여 1D에서 고차 정확도로 확장한다.
• 이산적 발산 없는 제약 조건과 재구성 및 유량 제한 기법을 조합하여, 2D에서 고차 정확도 PCP 스킴을 위한 두 가지 충분 조건을 확립한다.

실험 결과

연구 질문

• RQ1 특수 상대론적 자기유체역학에서 허용 상태 집합 𝒢를 정의하는 명시적이고 수학적으로 엄밀한 제약 조건은 무엇인가?
• RQ2 원시 변수와 유량에 대한 명시적 공식이 없음에도 불구하고, 허용 상태 집합 𝒢의 볼록성을 어떻게 엄밀히 증명할 수 있는가?
• RQ3 이산적 발산 없는 자기장 조건이 수치적 스킴에서 물리적 제약 유지에 어떤 영향을 미치는가?
• RQ4 표준 LxF 분할이 실패하는 비영인 자기장을 가진 RMHD에 대해, 일반화된 라크스-프리드리히스 분할을 어떻게 구성할 수 있는가?
• RQ5 2차원 RMHD 스킴에서 고차 정확도를 확보하면서 물리적 제약을 유지하기 위한 충분 조건는 무엇인가?

주요 결과

• 강한 비선형성에도 불구하고, 수학적 구조 분석이 가능한 두 가지 등가 명시적 형태의 허용 상태 집합 𝒢가 도출되었다.
• 초등형 표면의 제2 기본형의 준정부정성에 기반하여, 허용 상태 집합 𝒢가 볼록임을 증명하였다.
• 이산적 발산 없는 자기장 조건과 LxF 항의 볼록 조합을 결합하여, 일반화된 라크스-프리드리히스 분할 성질을 구성하였다.
• 이산적 발산 없는 조건이 물리적 제약 유지 성질과 이론적으로 밀접하게 연결되어 있음을 입증하였으며, 이는 스킴 설계에 핵심적인 통찰을 제공한다.
• 1D에서 표준 CFL 조건 하에서 증명 가능한 PCP 성질을 갖는 1차 정확도 LxF 유형 스킴이 입증되었으며, 이는 PCP 리미터를 통해 고차 정확도로 확장 가능하다.
• 2D에서 고차 정확도 PCP 스킴을 위한 두 가지 충분 조건이 유도되었으며, 다차원에서 안정적이고 물리적으로 일치하는 수치적 방법을 설계하는 데 이론적 기반을 마련하였다.