중력파와 별 붕괴를 위한 새로운 오픈소스 1차원 일반상대론 코드

초록

GR1D는 방사형 게이지·극좌표 슬라이스를 이용한 1차원 일반상대론 유체역학 코드를 오픈소스로 제공한다. 유한체적·PPM 재구성·근사 리만 해석기를 사용해 방정식을 풀며, 온도 의존성 미세물리 EOS와 탈전자화·중성미자 누출·가열 모델을 포함한다. 회전 효과를 1D에 효과적으로 구현하고, 다양한 검증 테스트와 40 M⊙ 별의 블랙홀 형성 시뮬레이션을 통해 기존 볼츠만 중성미자 수송 결과와 20 % 이내·10 % 이내의 일치를 보였다.

상세 분석

GR1D는 Romero‑Ibáñez‑Gourgoulhon이 제시한 Eulerian 형태의 일반상대론 유체역학(GRHD) 방정식을 기반으로 한다. 방사형 게이지(radial‑gauge)와 극좌표 슬라이스(polar‑slicing) 좌표계를 채택함으로써 3+1 분할에서 메트릭 성분이 단순화되고, 특히 라플라시안 형태의 중력 포텐셜이 사라져 수치 구현이 용이해진다. 코드의 핵심은 유한체적(FV) 스킴이며, 셀 경계에서의 물리량을 재구성하기 위해 3차 다항식 기반의 PPM(piecewise‑parabolic method)을 사용한다. 이는 급격한 충격파와 연속적인 흐름을 동시에 정확히 포착할 수 있게 해준다. 근사 리만 해석기(Riemann solver)로는 HLL(Harten‑Lax‑van Leer) 혹은 HLLE 변형을 적용해 특성 속도와 플럭스를 계산한다.

미세물리 EOS는 온도·밀도·전자 비율을 포함하는 테이블 형태로 제공되며, Lattimer‑Swesty, Shen 등 여러 공개 EOS를 바로 사용할 수 있다. 탈전자화(Deleptonization) 단계에서는 핵융합 전 단계에서 전자 포획에 의한 Y_e 감소를 근사적으로 모델링하고, 폭발 이후에는 중성미자 누출(Leakage)과 가열(Heating) 모듈을 통해 에너지와 레프톤 수송을 처리한다. 누출 모듈은 광학 두께에 기반한 로컬 방출률을 계산하고, 가열은 흡수된 중성미자 에너지를 유체에 재분배한다.

회전 효과는 1D에서 평균적인 원심력과 각운동량 보존을 가정한 유효 회전 방정식으로 구현한다. 구형 대칭을 유지하면서도 특정 각운동량 프로파일을 입력하면, 원심력 항이 압력·중력 균형에 추가되어 핵심핵 붕괴와 프로톤핵별(PNS) 진화에 미치는 영향을 정량화한다.

검증 테스트로는 고전적인 충격파 튜브 문제, Oppenheimer‑Snyder 붕괴, 그리고 정적 토러스 해(Static TOV) 솔루션을 재현했으며, 모두 이론적 해와 수치적 오차가 1 % 이하임을 확인했다. 실제 천체 물리 시뮬레이션에서는 15 M⊙와 40 M⊙ 전형적인 전진성 별 모델을 사용해 핵심핵 붕괴 → 프로톤핵별 형성 → 블랙홀 전이 과정을 추적했다. 특히 40 M⊙ 모델에서는 탈전자화와 누출‑가열을 포함한 전미세물리 EOS를 사용했을 때, 블랙홀 형성 시점이 기존 Boltzmann 중성미자 수송 시뮬레이션과 20 % 이내, 최종 프로톤핵별 질량이 10 % 이내 차이로 일치함을 보고했다. 이는 GR1D가 복잡한 중성미자 물리와 미세물리 EOS를 충분히 포착하면서도 계산 비용이 크게 절감된다는 강점을 보여준다.

전반적으로 GR1D는 오픈소스라는 접근성, 모듈식 설계, 그리고 검증된 수치 정확성을 결합해, 다차원 GR 코드 개발을 위한 테스트베드이자, 빠른 파라미터 스터디를 위한 실용적인 도구로 자리매김한다.