천문·천체물리
SPHINCS_BSSN 코드, 표 형식 상태방정식용 중성자별 병합 원시변수 복원 알고리즘 3종 제시
중성자별 병합 시뮬레이션은 고밀도 상태방정식의 영향을 크게 받는데, 정교한 모델은 보통 표(table) 형태로 제공된다. 라그랑지안 수치상대론 코드 SPHINCS_BSSN은 시공간을 적응격자에서, 유체는 자유롭게 움직이는 입자로 진화시키기 때문에 기존 수치상대론과 다른 보존변수-원시변수(conservative-to-primitive) 변환 알고리즘이 필요하다. 연구진은 3D·2D 뉴턴-랩슨 방법과 리더스(Ridders) 방법 기반 1D 근 찾기 알고리즘 등 세 가지를 제시했다. DD2 상태방정식을 쓴 실제 병합 시뮬레이션에서 3D 방법은 빠르고 안정적이며 평균 실패율이 1%를 크게 밑돌았고, 1D 리더스 방법은 느리지만 사실상 실패하지 않았다. 이에 3D 뉴턴-랩슨을 기본으로 삼고 1D 리더스 방법을 안전장치로 사용하기로 했다.
arXiv astro-ph.HE (고에너지)