[논문 리뷰] Event and Apparent Horizon Finders for 3+1 Numerical Relativity

Jonathan Thornburg|ePrints Soton (University of Southampton)|2005. 12. 30.

Astrophysical Phenomena and Observations참고 문헌 67인용 수 5

한 줄 요약

이 논문은 3+1 수치相对론 시뮬레이션에서 사라짐과 실질적 사건의 사건의 경계를 찾는 데 사용되는 수치 알고리즘과 소프트웨어 도구에 대한 종합적인 리뷰를 제공한다. 이들은 사건의 경계를 영향을 미치는 빛의 경로를 뒤로 적분하거나 표면을 통해 찾는 방법과, 경계가 외부로 트랩된 표면을 찾는 비선형 타원형 편미분방정식을 푸는 방법을 통해 실질적 사건의 경계를 찾는 방법을 상세히 설명하며, 축대칭 및 일반적인 시공간 모두에서 효율성, 정확성, 강건성을 강조한다.

ABSTRACT

Event and apparent horizons are key diagnostics for the presence and properties of black holes. In this article I review numerical algorithms and codes for finding event and apparent horizons in numerically-computed spacetimes, focusing on calculations done using the 3+1 ADM formalism. There are 3 basic algorithms for finding event horizons, based respectively on integrating null geodesics \emph{forwards} in time, integrating null geodesics \emph{backwards} in time, and integrating null \emph{surfaces} backwards in time. The last of these is generally the most efficient and accurate. There are a large number of apparent-horizon finding algorithms, with differing trade-offs between speed, robustness, accuracy, and ease of programming. In axisymmetry, shooting algorithms work well and are fairly easy to program. In slices with no continuous symmetries, Nakamura et al.'s algorithm and elliptic-PDE algorithms are fast and accurate, but require good initial guesses to converge. In many cases Schnetter's "pretracking" algorithm can greatly improve an elliptic-PDE algorithm's robustness. Flow algorithms are generally quite slow, but can be very robust in their convergence.

연구 동기 및 목표

• 3+1 ADM 형식을 통해 계산된 시공간에서 사건의 경계와 실질적 사건의 경계를 탐지하기 위한 수치 알고리즘을 체계화하고 리뷰하는 것.
• 속도, 정확성, 강건성, 구현 용이성 측면에서 다양한 사건의 경계 탐지 기법의 강점과 약점을 식별하고 비교하는 것.
• 시공간 대칭성(예: 축대칭 대비 일반적)과 계산 제약 조건에 따라 적절한 알고리즘을 선택하는 데 도움을 주는 것.
• 공개된 수렴 가능한 소프트웨어를 목록화하고 수치 상대론 소프트웨어 엔지니어링의 최선의 실천 방법을 강조하여 코드 재사용과 공유를 촉진하는 것.
• 새로운 방법이 등장함에 따라 공동체 기여를 유도하여 계속 업데이트되고 확장되는 살아있는 참고 자료로 기능하는 것.

제안 방법

• 사건의 경계는 후처리 단계에서 세 가지 주요 알고리즘을 사용하여 찾는다: 빛의 경로를 앞으로 적분하는 방법, 빛의 경로를 뒤로 적분하는 방법, 그리고 빛의 표면을 뒤로 적분하는 방법—이 중 후자가 가장 효율적이고 정확하다.
• 실질적 사건의 경계는 경계가 외부로 트랩된 표면(MOTS) 조건을 만족하는 비선형 타원형 편미분방정식을 풀어 찾는다. 이 방정식은 3-메트릭, 그 공간 도함수, 그리고 외부 곡률을 포함한다.
• 축대칭 상황에서는 射격 알고리즘과 영점 찾기 방법이 효과적이며 쉽게 구현할 수 있다.
• 일반적인 시공간에서는 나카무라 등이 제안한 스펙트럼 적분-반복 방법과 타원형 편미분방정식 해법기가 빠르고 정확하지만, 좋은 초기 추정치가 필요하다.
• 사전 추적(pretracking)은 타원형 편미분방정식 해법기의 수렴을 향상시키기 위해 시험적인 사건의 경계를 진짜 해로 향해 진화시켜 사용된다.
• 유연성(강건성)을 우선시할 경우 유동 알고리즘(flow algorithms)이 사용되며, 이는 역동적 진화를 통해 MOTS에 수렴하도록 한다.

실험 결과

연구 질문

• RQ1 수치적으로 계산된 시공간에서 사건의 경계를 찾는 데 가장 효율적이고 정확한 수치 방법은 무엇인가?
• RQ2 다양한 시공간 대칭성에서 수렴성, 강건성, 계산 비용 측면에서 실질적 사건의 경계 탐지 알고리즘 간의 비교는 어떻게 이루어지는가?
• RQ3 타원형 편미분방정식 기반 실질적 사건의 경계 탐지기에서 초기 추정치의 선택이 수렴에 어떤 영향을 미치며, 이를 어떻게 개선할 수 있는가?
• RQ4 소프트웨어 엔지니어링 실천 방법은 수치 상대론에서 사건의 경계 탐지 코드의 신뢰성, 재사용성, 유지보수성에 어떻게 기여할 수 있는가?
• RQ5 공개된 수치 라이브러리와 해법기(예: 비선형 방정식 및 상미분방정식용)를 어떻게 활용하여 강력한 사건의 경계 탐지 도구를 구축할 수 있는가?

주요 결과

• 빛의 표면을 뒤로 적분하는 방법이 사건의 경계 탐지에서 가장 효율적이고 정확하며, 앞서 적분하거나 경로 기반 접근보다 뛰어나다.
• 타원형 편미분방정식 해법기와 나카무라 등이 제안한 스펙트럼 방법은 일반적인 시공간에서 실질적 사건의 경계 탐지에 높은 정확도와 빠른 속도를 제공하지만, 수렴을 위해 좋은 초기 추정치가 필요하다.
• 샨너터의 사전 추적 알고리즘이 타원형 편미분방정식 기반의 사건의 경계 탐지기의 강건성을 크게 향상시켜 해가 수렴하도록 이끈다.
• 유동 알고리즘은 수렴성에서 매우 강건하지만 계산 비용이 높아, 속도보다 신뢰성이 더 중요한 경우에 적합하다.
• 비선형 대수 방정식을 해결할 때는 속도와 신뢰성의 균형을 잘 맞추는 데서 유리한 점을 고려해 밴 위인지렌-데커-브렌트 알고리즘을 권장한다.
• 유연한 상미분방정식 해법기인 LSODE와 RKF45는 사건의 경계 탐지 파이프라인의 시간 적분 작업에 필수적이며, netlib를 통해 무료로 이용 가능하여 신뢰성 있고 효율적인 구현이 가능하다.