[논문 리뷰] A strongly hyperbolic first-order CCZ4 formulation of the Einstein equations and its solution with discontinuous Galerkin schemes

Michael Dumbser, Federico Maria Guercilena외 3명|arXiv (Cornell University)|2017. 07. 31.

Numerical methods for differential equations인용 수 8

한 줄 요약

이 논문은 제약 조건 감쇠를 포함한 CCZ4 체계에서 유도된 강력 타원형(first-order) 형식인 FO-CCZ4를 제시한다. 이는 안정적인 고차수 불연속 갈레르킨(DG) 시뮬레이션을 가능하게 하며, 첫 번째로 고차수 3차원 이동 펄서(binary black hole) 시스템에 대한 DG 시뮬레이션을 성공적으로 수행하였다. 이는 다차원 진공 테스트를 통한 검증과 함께, 적응형 메쉬 구조 및 서브셀 리미터(subcell limiter)를 통해 향상되었다.

ABSTRACT

We present a strongly hyperbolic first-order formulation of the Einstein equations based on the conformal and covariant Z4 system (CCZ4) with constraint-violation damping, which we refer to as FO-CCZ4. As CCZ4, this formulation combines the advantages of a conformal and traceless formulation, with the suppression of constraint violations given by the damping terms, but being first order in time and space, it is particularly suited for a discontinuous Galerkin (DG) implementation. The strongly hyperbolic first-order formulation has been obtained by making careful use of first and second-order ordering constraints. A proof of strong hyperbolicity is given for a selected choice of gauges via an analytical computation of the entire eigenstructure of the FO-CCZ4 system. The resulting governing partial differential equations system is written in non-conservative form and requires the evolution of 58 unknowns. A key feature of our formulation is that the first-order CCZ4 system decouples into a set of pure ordinary differential equations and a reduced hyperbolic system of partial differential equations that contains only linearly degenerate fields. We implement FO-CCZ4 in a high-order path-conservative arbitrary-high-order-method-using-derivatives (ADER)-DG scheme with adaptive mesh refinement and local time-stepping, supplemented with a third-order ADER-WENO subcell finite-volume limiter in order to deal with singularities arising with black holes. We validate the correctness of the formulation through a series of standard tests in vacuum, performed in one, two and three spatial dimensions, and also present preliminary results on the evolution of binary black-hole systems. To the best of our knowledge, these are the first successful three-dimensional simulations of moving punctures carried out with high-order DG schemes using a first-order formulation of the Einstein equations.

연구 동기 및 목표

• 고차수 불연속 갈레르킨(DG) 방법에 적합한 일阶, 강력 타원형인 아인슈타인 방정식의 형식을 개발하기 위해.
• 공형 및 트레이스리스 형식의 장점을 결합하고 제약 위반 감쇠를 통해 수치적 안정성을 향상시키기 위해.
• 특히 이원자 시스템을 포함한 블랙홀 시공간의 고차수 안정적 시뮬레이션을 DG 방법을 통해 가능하게 하기 위해.
• 블랙홀 진화 과정에서 발생하는 특이점 등의 수치적 과제를 고도화된 리미터 및 적응형 메쉬 구조로 해결하기 위해.

제안 방법

• 강력 타원성을 보장하기 위해 일阶 및 이阶 순서 제약 조건을 신중히 적용하여 CCZ4 체계의 일阶 형식을 유도하기 위해.
• 58개의 진화 변수를 가진 비보존형 타원형 PDE 체계를 구성하여 순수 미분방정식과 선형으로 탈퇴한 필드를 포함하는 감소된 타원형 체계로 분해하기 위해.
• 국부적 시간 스텝핑과 적응형 메쉬 구조를 통합한 고차수 임의의 고차수 방법을 사용하는 ADER-DG 스킴에 시스템을 구현하기 위해.
• 블랙홀 시뮬레이션에서 발생하는 특이점을 다루기 위해 제3차 ADER-WENO 서브셀 유한체적 리미터를 통합하기 위해.
• 표준 진공 테스트를 통해 형식을 검증하기 위해(일차원, 이차원, 삼차원 모두).
• 이중 블랙홀 시스템의 진화에 대해 해당 스킴을 적용하여, 첫 번째로 고차수 DG와 일阶 형식을 사용한 3차원 이동 펄서 시뮬레이션을 달성하기 위해.

실험 결과

연구 질문

• RQ1 아인슈타인 방정식의 일阶, 강력 타원형 형식을 구성할 수 있는가? 이는 고차수 DG 이산화를 지원할 수 있는가?
• RQ2 표준 게이지 선택 조건 하에서 FO-CCZ4 형식이 강력 타원성을 어떻게 확보하는가? 고유구조(eigenstructure)는 어떠한가?
• RQ3 적응형 메쉬 구조와 국부적 시간 스텝핑을 통합한 고차수 DG 방법이 이동 펄서를 포함한 이중 블랙홀 시스템을 성공적으로 진화시킬 수 있는가?
• RQ4 ADER-WENO 서브셀 리미터는 블랙홀 특이점 근처에서 시뮬레이션을 안정화시키는 데 얼마나 효과적인가?
• RQ5 FO-CCZ4 형식은 다차원 진공 테스트와 이중 블랙홀 진화에서 성능과 정확도가 어느 정도인가?

주요 결과

• 선택된 게이지 조건 하에서 분석적 계산을 통해 전체 고유구조를 확인함으로써, FO-CCZ4 형식이 강력 타원성을 갖는 것으로 확인되었다.
• 체계는 순수 상미분방정식과 선형으로 탈퇴한 필드만을 포함하는 감소된 타원형 체계로 분해되며, 이는 수치적 처리를 단순화시킨다.
• 이 형식은 고차수 불연속 갈레르킨 방법을 사용한 이동 펄서 이중 블랙홀 시스템에 대한 첫 번째 성공적인 3차원 시뮬레이션을 가능하게 하였다.
• 일차원, 이차원, 삼차원 진공 테스트 전역에서 높은 정확도와 안정성을 확보하였다.
• 적응형 메쉬 구조, 국부적 시간 스텝핑, ADER-WENO 서브셀 리미터의 조합은 블랙홀 특이점 근처의 수치적 불안정성을 효과적으로 관리하였다.
• 이중 블랙홀 진화에 대한 초보적 결과는 일阶 형식을 사용한 고차수 DG 방법이 복잡한 시공간 역학을 다룰 수 있음을 시사한다.