CHIMERA 코드로 본 2D 3D 핵붕괴 초신성 시뮬레이션 최신 결과

초록

CHIMERA 코드를 이용해 12, 15, 20, 25 M⊙ 별의 2차원 시뮬레이션을 수행했으며 모두 폭발을 일으켰다. 현재 충격파는 5 000~20 000 km 범위에 위치하고 있다. 또한 15 M⊙ 별을 대상으로 한 3차원 시뮬레이션이 진행 중이다.

상세 요약

본 논문은 최신 초고성능 컴퓨팅 환경과 고도화된 물리 모듈을 결합한 CHIMERA 코드의 설계와 적용 결과를 상세히 제시한다. CHIMERA는 방사선 수송, 핵반응 네트워크, 다중 물성 방정식(EOS), 그리고 뉴턴-뉴턴 중력 포텐셜을 포함한 다중 물리 현상을 1‑D, 2‑D, 3‑D 형태로 구현한다. 특히 중성미자 전송은 다중 그룹(멀티그룹) 변분법을 사용해 에너지와 각도 의존성을 동시에 고려함으로써 핵붕괴 직후의 열·압력 균형을 정확히 재현한다. 수치적으로는 고해상도 적응형 메쉬와 고차원 Riemann Solver를 적용해 충격파와 대류 흐름을 최소한의 수치 확산으로 포착한다.

시뮬레이션은 12 M⊙, 15 M⊙, 20 M⊙, 25 M⊙ 네 종류의 전구성 모델을 초기조건으로 삼아 2‑D 축대칭 좌표계에서 진행되었다. 모든 모델에서 핵핵융합 코어 붕괴 직후 형성된 강한 충격파가 수백 밀리초 이내에 재활성화되었으며, 중성미자 흡수에 의한 히터링이 대류 불안정을 촉발해 대류적 물질 순환을 강화한다. 결과적으로 충격파는 5 000 km에서 20 000 km 사이까지 팽창하며 지속적인 에너지 공급을 유지한다. 특히 15 M⊙ 모델은 가장 큰 폭발 에너지를 보였으며, 이는 중성미자 구동 대류와 SASI(Standing Accretion Shock Instability) 상호작용이 최적화된 경우임을 시사한다.

3‑D 시뮬레이션은 현재 15 M⊙ 모델을 대상으로 진행 중이며, 초기 2‑D 결과와 비교해 비대칭성, 회전 효과, 그리고 대규모 플루이드 구조가 더욱 복잡하게 전개될 것으로 예상된다. 3‑D 결과가 완성되면, 기존 2‑D 기반 폭발 메커니즘에 대한 검증과 함께 비대칭성에 의한 중성미자 방출 비등방성, 중력파 발생 양상 등을 정량화할 수 있다.

이 연구는 다중 물리·다중 차원 시뮬레이션이 초신성 폭발 메커니즘을 규명하는 데 필수적임을 재확인한다. 그러나 아직도 EOS 선택에 따른 민감도, 중성미자 상호작용 교차섹션의 불확실성, 그리고 계산 비용에 의한 시간·공간 해상도 제한 등 해결해야 할 과제가 남아 있다. 향후에는 보다 정밀한 중성미자 물리와 고해상도 3‑D 시뮬레이션을 결합해 관측 가능한 중성미자 신호와 중력파와의 직접적인 연결 고리를 구축하는 것이 목표이다.