[논문 리뷰] A 3+1 Computational Scheme for Dynamic Spherically Symmetric Black Hole Spacetimes -- II: Time Evolution

Jonathan Thornburg|ArXiv.org|1999. 06. 06.

Black Holes and Theoretical Physics참고 문헌 86인용 수 7

한 줄 요약

이 논문은 블랙홀 배제 기법을 사용하여 구형 대칭 블랙홀 시공간과 스칼라 장을 결합한 3+1 수치적 방법을 제안한다. 자유 진화 방식으로 에딩턴-핑켈스타인 유사 좌표계를 사용하고, 공간 및 시간에서 네째계수 유한차분법과 룬게-쿠타 적분기를 적용하여, t = 4000m를 초과하는 안정적인 장기 수치 시뮬레이션을 실현하였으며, 계량장 성분과 제약 조건에서 오차가 10^-5 이하로 유지된다.

ABSTRACT

This is the second in a series of papers describing a 3+1 computational scheme for the numerical simulation of dynamic black hole spacetimes. We discuss the numerical time-evolution of a given black-hole-containing initial data slice in spherical symmetry. We avoid singularities via the "black-hole exclusion" or "horizon boundary condition" technique, where the slices meet the black hole's singularity, but on each slice a spatial neighbourhood of the singularity is excluded from the domain of the numerical computations. After first discussing some of the key design choices which arise with the black hole exclusion technique, we then give a detailed description of our numerical evolution scheme for spherically symmetric scalar field evolution, assuming that a black hole is already present on the initial slice. We use a free evolution, with Eddington-Finkelstein-like coordinates and the inner boundary placed at a fixed coordinate radius well inside the horizon. Our numerical scheme is based on the method of lines (MOL), where spacetime PDEs are first finite differenced in space only, yielding a system of coupled ODEs for the time evolution of the field variables along the spatial-grid-point world lines. These ODEs are then time-integrated by standard methods. We use 4th order finite differencing in both space and time, with 5 and/or 6 point spatial molecules (off-centered near the grid boundaries), and a Runge-Kutta time integrator. The spatial grid is smoothly nonuniform, but not adaptive. We present numerical black hole + scalar field evolutions showing that this scheme is stable, can evolve "forever" (we have gone to t > 4000m), and is very accurate. At a resolution Delta_r/r = 3% near the horizon, typical errors in g_ij(K_ij) at t=100m are <= 1e-5(3e-7), and the energy constraint is < 3e-5.

연구 동기 및 목표

• 구형 대칭에서 스칼라 장을 가진 동적 블랙홀 시공간을 시뮬레이션하기 위한 강력한 3+1 수치적 프레임워크를 개발하는 것.
• 계산 도메인에서 블랙홀 특이점 주변의 공간 영역을 제거함으로써 수치적 특이점을 해결하는 것.
• 고차수 유한차분법과 표준 ODE 적분기를 사용하여 높은 정확도와 장기적 안정성을 확보하는 것.
• 수렴성과 제약 조건 유지 여부를 수치적 진화를 통해 검증하는 것.
• 통제 가능한 수치 오차로 구형 대칭에서 블랙홀 형성과 진화를 시뮬레이션할 수 있는 기초를 마련하는 것.

제안 방법

• 블랙홀 특이점 주변의 공간 영역을 계산 도메인에서 제거함으로써 블랙홀 배제 기법을 적용하여 좌표 특이점을 피한다.
• 시공간의 수평면에서 정규성을 확보하고 안정적인 진화를 가능하게 하기 위해 에딩턴-핑켈스타인 유사 좌표계를 사용한다.
• 선형 방법(MOL)을 적용하여 시공간 편미분방정식(PDE)을 시간에 대해 이산화하지 않고 공간에 대해서만 이산화함으로써, 세계선 沿해 필드 진화를 위한 연립 상미분방정식(ODE) 시스템을 도출한다.
• 공간에 대해 네째계수 유한차분법를 적용하며, 5점 또는 6점 스테이션을 사용하고 경계 근처에서는 중심에서 약간 이격된 방식으로 정확도를 유지한다.
• 시간에 대해 네째계수 룬게-쿠타 적분기를 사용하여 ODE 시스템을 시간에 따라 진행시킨다.
• 공간 격자는 부드럽게 비균일하지만 적응형이 아니며, 시공간의 수렴성에 따라 Δr/r ≈ 3% 정도의 해상도를 갖는다.

실험 결과

연구 질문

• RQ1 3+1 수치적 방법이 스칼라 장과 결합된 구형 대칭 블랙홀 시공간의 안정적이고 장기적인 진화를 달성할 수 있는가?
• RQ2 블랙홀 배제 기법은 수치 정확도를 유지하면서 좌표 특이점을 효과적으로 피할 수 있는가?
• RQ3 이 맥락에서 네째계수 유한차분법과 룬게-쿠타 시간 적분기를 사용할 경우 도달할 수 있는 정확도 수준은 무엇인가?
• RQ4 에너지 제약 조건과 같은 제약 조건이 장기적인 진화 동안 얼마나 잘 유지되는가?
• RQ5 블랙홀 내부 경계를 시공간의 수렴성 이하의 좌표 반경에 고정하더라도, 시스템이 임의로 오랜 시간 동안 붕괴 없이 정확하게 진화할 수 있는가?

주요 결과

• t = 4000m를 초과하는 시간 동안 안정적인 진화를 달성하여 장기적인 수치적 강건성을 입증하였다.
• 시공간의 수렴성에 따라 Δr/r ≈ 3% 수준에서, t = 100m 시점에 계량장 성분 g_ij의 오차는 ≤ 1×10^-5 이하, 외적 곡률 K_ij의 오차는 ≤ 3×10^-7 이하로 유지되었다.
• 에너지 제약 조건은 동일한 시간 간격 동안 3×10^-5 이하로 유지되어 강력한 제약 조건 충족을 나타내었다.
• 네째계수 공간 및 시간 이산화 방식을 통해 높은 정확도와 수렴성을 확보하였다.
• 내부 경계를 시공간의 수렴성 이하의 좌표 반경에 고정하더라도, 시스템은 안정적이고 정확하게 유지되었다.
• 결과적으로 블랙홀 배제 기법과 에딩턴-핑켈스타인 유사 좌표계의 조합이 좌표 특이점 없이 신뢰할 수 있는 블랙홀 시공간 시뮬레이션을 가능하게 한다는 것이 확인되었다.