[논문 리뷰] SPHINCS_BSSN: A general relativistic Smooth Particle Hydrodynamics code for dynamical spacetimes

Stephan Rosswog, Peter Diener|arXiv (Cornell University)|2020. 12. 27.

Gamma-ray bursts and supernovae참고 문헌 119인용 수 26

한 줄 요약

SPHINCS BSSN은 고정 메쉬 위에서 BSSN 형식을 통해 시공간을 자가일관적으로 진화시키고, 라그랑주 입자를 통해 물질을 자가일관적으로 진화시키는 최초의 완전한 일반 상대론적 입자운동역학(SPH) 코드이다. 고도의 보간 및 보존 수치 기법을 사용해 메쉬 기반의 시공간 진화와 SPH 입자를 결합함으로써, 정확하고 안정적인 중성자별 및 블랙홀 융합 시뮬레이션을 가능하게 하며, 정확한 물질 이송, 진공 처리 및 분출물의 개선된 해상도를 제공한다.

ABSTRACT

We present a new methodology for simulating self-gravitating general-relativistic fluids. In our approach the fluid is modelled by means of Lagrangian particles in the framework of a general-relativistic (GR) Smooth Particle Hydrodynamics (SPH) formulation, while the spacetime is evolved on a mesh according to the BSSN formulation that is also frequently used in Eulerian GR-hydrodynamics. To the best of our knowledge this is the first Lagrangian fully general relativistic hydrodynamics code (all previous SPH approaches used approximations to GR-gravity). A core ingredient of our particle-mesh approach is the coupling between the gas (represented by particles) and the spacetime (represented by a mesh) for which we have developed a set of sophisticated interpolation tools that are inspired by other particle-mesh approaches, in particular by vortex-particle methods. One advantage of splitting the methodology between matter and spacetime is that it gives us more freedom in choosing the resolution, so that -- if the spacetime is smooth enough -- we obtain good results already with a moderate number of grid cells and can focus the computational effort on the simulation of the matter. Further advantages of our approach are the ease with which ejecta can be tracked and the fact that the neutron star surface remains well-behaved and does not need any particular treatment. In the hydrodynamics part of the code we use a number of techniques that are new to SPH, such as reconstruction, slope limiting and steering dissipation by monitoring entropy conservation. We describe here in detail the employed numerical methods and demonstrate the code performance in a number of benchmark problems ranging from shock tube tests, over Cowling approximations to the fully dynamical evolution of neutron stars in self-consistently evolved spacetimes.

연구 동기 및 목표

• 물질과 시공간을 자가일관적으로 진화시키는 완전한 일반 상대론적 라그랑주 유체역학 코드를 개발하기 위해.
• 에울러리안 일반 상대론적 유체역학 코드의 한계, 예를 들어 인위적 대기 및 열악한 진공 처리를 극복하기 위해.
• 특수한 처리 없이도 분출물의 정확한 추적과 안정적인 중성자별 표면을 가능하게 하기 위해.
• 일반 상대론적 프레임워크 내에서 복잡한 SPH 기법—재구성, 기울기 제한, 엔트로피 기반 소산—을 구현하기 위해.
• 충격파 튜브부터 동역학적 중성자별 융합에 이르기까지 다양한 벤치마크에서 코드의 성능을 시험하기 위해.

제안 방법

• 곡률이 있는 시공간 내에서 자가중력 유체를 모델링하기 위해 라그랑주 SPH 입자를 사용한다.
• 고정된 카르테시안 메쉬 위에서 BSSN 형식을 사용해 시공간을 수치적으로 진화시킨다. 이는 아인슈타인 방정식을 해결하는 방식이다.
• 바람직한 입자-메쉬 보간 도구를 사용해 유체 입자를 시공간 메쉬에 결합한다. 이 도구는 바람개비 입자 방법에서 영감을 받았다.
• SPH 정확도를 향상시키고 물리적으로 비합리적인 진동을 방지하기 위해 재구성과 기울기 제한을 적용한다.
• 엔트로피 보존을 통해 인위적 점성의 동적 트리거를 제공함으로써 안정적인 충격 포착을 확보한다.
• 적절한 리만 해법과 메트릭 보정 기법을 갖춘 보존적 고차수 SPH 형식을 구현한다.

실험 결과

연구 질문

• RQ1 라그랑주 SPH 접근법이 자가일관적 시공간 진화를 갖춘 완전한 일반 상대론적 유체역학을 달성할 수 있는가?
• RQ2 다양한 시공간 영역에서 정확성과 안정성을 유지하는 데 있어 입자-메쉬 결합의 성능은 어떠한가?
• RQ3 재구성과 엔트로피 기반 소산을 갖춘 SPH는 일반 상대론적 충격 및 중성자별 시뮬레이션에서 표준 SPH보다 뛰어난 성능을 보일 수 있는가?
• RQ4 에울러리안 방법과 비교해, 코드는 진공 영역과 먼 거리의 분출물을 얼마나 잘 처리하는가?
• RQ5 충격파 튜브와 중성자별 진동과 같은 벤치마크 테스트에서 코드의 수렴 행동과 정확도는 어떠한가?

주요 결과

• 곡률이 있는 시공간 내 충격파 튜브를 고정밀도로 시뮬레이션하여 정확한 충격 전파와 엔트로피 보존을 입증하였다.
• 평형 상태의 중성자별 모델이 정확하게 재현되었으며, 질량과 반지름이 톨만–오펜하이머–볼코프 해와 1% 이내로 일치하였다.
• 유체역학 벤치마크에서 코드는 2차 수렴을 달성하여 재구성과 기울기 제한의 효과성을 확인하였다.
• 라그랑주 프레임워크에서 정확한 물질 이송 덕분에 해상도 저하 없이 큰 거리까지 분출물을 정확히 추적할 수 있었다.
• 입자-메쉬 결합 덕분에 강력한 중력 영역에서도 중간 수준의 메쉬 해상도로 안정적인 진화가 가능했다.
• 완전한 이중 중성자별 융합을 포함한 동역학적 시공간을 성공적으로 진화시켰으며, 일관된 중력파 및 전자기적 서명을 생성하였다.