다중군 방사선 유체역학을 위한 새로운 압축성 솔버 CASTRO

초록

CASTRO 코드에 다중군 방사선 유체역학을 O(v/c) 차수까지 정확히 구현하였다. 방사선‑유체 결합을 하이퍼볼릭, 주파수 공간 이동을 별도, 확산·소스‑싱크를 파라볼릭으로 분리해 각각 고차 Godunov과 후진 오일러 방식으로 풀며, AMR 기반 격자에서 중성미자와 광자를 모두 다룰 수 있다.

상세 요약

본 논문은 압축성 천체 물리 시뮬레이션 코드 CASTRO에 다중군 방사선 유체역학(Multi‑Group Radiation Hydrodynamics, MGRHD) 모듈을 추가한 방법론을 상세히 제시한다. 핵심은 코밋 프레임(comoving‑frame) 접근법을 채택해 방사선 전송 방정식을 물질 속도 v에 대한 1차 항까지 보존함으로써, 상대론적 효과를 O(v/c) 수준에서 정확히 기술한다는 점이다. 이를 위해 플럭스 제한 확산(Flux‑Limited Diffusion, FLD) 근사를 적용해 방사선 플럭스를 고전적인 확산 형태로 표현하면서도 광속에 근접한 자유 경로를 제한한다.

수치 알고리즘은 시스템을 세 부분으로 분할한다. 첫 번째는 방사선과 유체가 상호작용하는 하이퍼볼릭 서브시스템으로, 질량·운동량·에너지 보존 방정식에 방사선 압력·에너지 교환 항을 포함한다. 이 부분은 고차 정확도의 Godunov 스킴(재구성‑리컨스트럭션 + Riemann solver)으로 명시적으로 풀어, 충격파와 급격한 흐름 변화를 정확히 포착한다. 두 번째는 방사선 에너지 밀도가 주파수(에너지) 공간에서 이동하는 항을 담당하는 ‘주파수 공간 어드벡션’이다. 이 역시 고전적인 Godunov 방식으로 구현해, 그룹 간 에너지 전이를 안정적으로 계산한다. 세 번째는 확산·소스‑싱크 항을 포함하는 파라볼릭 서브시스템으로, 시간 스텝 제한이 심각한 확산 문제를 완화하기 위해 1차 후진 오일러(Backward Euler) 방법을 사용한 암시적 해법을 적용한다. 암시적 단계는 다중군 선형 시스템을 효율적으로 풀기 위해 멀티그리드(preconditioned Krylov) 기법을 결합한다.

CASTRO는 블록 구조 적응형 격자(Adaptive Mesh Refinement, AMR)를 이용해 공간·시간 해상도를 동시에 가변적으로 조정한다. 각 레벨의 격자는 논리적으로 직사각형이며, 시간 정밀도는 레벨마다 반비례하게 감소한다(시간 서브사이클링). 이 설계는 고해상도가 요구되는 핵심 영역(예: 초신성 코어, 충돌 충격면)에서 계산 비용을 크게 절감한다. 또한, 다중군 구현은 중성미자와 광자 모두에 적용 가능하도록 설계돼, 에너지 그룹별 흡수·산란 계수를 입력하면 자동으로 방사선‑물질 상호작용을 처리한다.

논문은 여러 검증 테스트를 제시한다. 첫째, 방사선‑유체 결합 파동 테스트에서 O(v/c) 항을 포함했을 때와 제외했을 때의 차이를 정량적으로 비교한다. 둘째, 다중군 확산 문제에서 그룹 간 에너지 교환이 정확히 보존되는지를 확인한다. 셋째, 실제 천체 물리 시나리오인 핵 붕괴 초신성 모델에 적용해, 중성미자 방출과 충격파 전파가 기존 단일군 모델보다 더 정밀하게 재현되는 것을 보여준다. 전체적으로, 제안된 알고리즘은 정확도, 안정성, 확장성 측면에서 기존 방법을 능가한다는 결론을 도출한다.