[논문 리뷰] New numerical framework for the generalized Baumgarte-Shapiro-Shibata-Nakamura formulation: The vacuum case for spherical symmetry

M. A. Alcoforado, R. F. Aranha외 2명|arXiv (Cornell University)|2021. 05. 19.

Pulsars and Gravitational Waves Research참고 문헌 62인용 수 7

한 줄 요약

이 논문은 구형 대칭에서 일반화된 BSSN 표현식을 해결하기 위한 새로운 오픈소스 스펙트럼 코드인 RIO를 제시한다. 갈레르킨-콜로케이션 방법을 기반으로 하며, 원점과 무한원에서의 정칙성을 내재한 기저 함수를 사용함으로써 명시적 정규화 없이 약속 조건에서 지수 수렴을 달성하며, 진공 구형 대칭 시험에서 해밀토니안 및 운동량 약속 조건에서 거의 머신 정밀도 수준의 정확도를 입증한다.

ABSTRACT

Here we report a developed high performance and simplified version of the code denominated RIO, which can be easily extended, for the generalized BSSN formulation. We implement a code which is regular at the center of symmetry, without use a special procedure for regularization, as usual. We get exponential convergence for constraints. The numerical algorithm is based on the Galerkin-Collocation method developed successfully for diverse physical scenarios by the Numerical Relativity Group at UERJ. For the sake of clarity in presentation, we consider here the most simple case to display the most salient features of the procedure. Thus, we focus on the definite tests of the new numerical framework. The timing and performance of the code show that we can reach a better accuracy close to the machine precision, for the Hamiltonian and momentum constraints. RIO will be an open source code; currently is under continuous development to consider more general and realistic problems.

연구 동기 및 목표

• 구형 대칭에서 일반화된 BSSN 표현식을 위한 고성능이며 확장 가능한 수치 프레임워크를 개발하기 위해.
• 원점과 극축에서의 보조 정규화 절차가 필요 없도록 정칙성을 스펙트럼 기저 함수에 내재시켜 제거하기 위해.
• 갈레르킨-콜로케이션 스펙트럼 방법을 사용하여 약속 위반에서 지수 수렴을 달성하기 위해.
• 향후 2차원 및 3차원 문제, 예를 들어 중력 수축과 브릴 웨이브 수축 문제로의 확장 기반을 마련하기 위해.
• 우에르주(UERJ)에서 수치相对론 연구를 위한 통합적이고 유지보수 가능한 코드베이스(RIO)를 구축하기 위해.

제안 방법

• 갈레르킨-콜로케이션(GC) 스펙트럼 방법을 적용하여 갈레르킨 투영과 가우스-레지온드르-로바트 점에서의 콜로케이션을 결합한다.
• L2[0, ∞) 힐베르트 공간에서 기저 함수를 선택하여 r=0 및 무한원에서의 내재적 정칙성을 확보한다.
• 이 기저 함수에 대한 스펙트럼 전개를 통해 편미분방정식 시스템을 자동화된 동역학 시스템으로 변환한다.
• 시간 진화는 스펙트럼 계수에 적용된 고정단계 4차 룬게-쿠타 적분기법을 사용한다.
• 결함이 발생할 수 있는 점들을 방지하기 위해 콜로케이션 점들을 분포시켰으며, 경계 조건은 기저 함수 선택을 통해 자연스럽게 처리된다.
• 이 방법은 알려진 초기 조건과 약속 조건을 가진 진공 및 구형 대칭 케이스에서 검증되었다.

실험 결과

연구 질문

• RQ1 일반화된 BSSN 표현식의 구형 대칭에서 스펙트럼 갈레르킨-콜로케이션 방법이 해밀토니안 및 운동량 약속 조건에서 지수 수렴을 달성할 수 있는가?
• RQ2 기저 함수 선택을 통해 원점과 무한원에서의 정칙성을 내재적으로 달성할 수 있으며, 특수한 정규화 절차가 필요 없게 할 수 있는가?
• RQ3 이 프레임워크의 성능 및 정확도는 기존 BSSN 구현체와 비교해 약속 위반과 수렴 속도 측면에서 어떻게 다른가?
• RQ4 이 프레임워크는 중력 수축 및 원통좌표계에서의 2차원 문제와 같은 더 복잡한 시나리오로 확장될 수 있는가?
• RQ5 RIO 코드는 진공 구형 대칭 시험에서 달성 가능한 정확도와 계산 효율성이 어느 정도인가?

주요 결과

• RIO 코드는 해밀토니안 및 운동량 약속 조건에서 지수 수렴을 달성하여 스펙트럼 정밀도를 나타낸다.
• 약속 위반은 거의 머신 정밀도 수준으로 감소하여 높은 수치 정밀도를 입증한다.
• 기저 함수 선택을 통해 원점과 무한원에서 정칙성이 자연스럽게 강제되며, 특수한 정규화 기법이 필요 없어진다.
• 고정단계 4차 룬게-쿠타 적분기법을 사용한 시간 진화에서 코드는 안정적이고 양호한 성능을 보인다.
• 이 프레임워크는 확장 가능하며, 향후 중력 수축 및 브릴 웨이브 수축 문제를 포함한 2차원 및 3차원 문제로의 응용을 위해 활발히 개발 중이다.
• RIO 코드는 오픈소스이며, 수치相对론 연구에서의 지속적 개발과 재사용을 위해 설계되었다.