CASTRO의 고성능 회색 복사 유체역학 솔버: 하이퍼볼릭·패러볼릭 분할 구현

초록

CASTRO 코드에 플럭스 제한 회색 복사 전이(FRD) 모듈을 도입하여, 복사를 하이퍼볼릭(명시적 고든법)과 패러볼릭(후방 오일러 암시적) 부분으로 분리해 AMR 격자에서 효율적으로 계산한다. 혼합 프레임 접근을 사용해 에너지 보존을 유지하고, 방사선‑유체 결합 속도와 압력 변화를 정확히 포착한다.

상세 분석

본 논문은 CASTRO(Compressible Astrophysical Solver) 코드에 회색 복사 전이(FRD) 모듈을 추가한 방법론을 상세히 제시한다. 핵심은 복사‑유체 방정식을 혼합 프레임(mixed‑frame)으로 기술하고, 이를 하이퍼볼릭·패러볼릭 두 서브시스템으로 시간 분할(time‑splitting)한다는 점이다. 하이퍼볼릭 부분은 방사선 압력 구배 λ∇Er 를 포함한 수정된 음속 cs 를 이용해 고든법(Godunov) 기반의 고차원 비압축형 Riemann 솔버로 명시적 적분한다. 이때 고유값은 u±cs, u, u 로 단순화되며, λ와 Eddington factor f 사이의 관계식(3−f²=λ+1)을 가정해 근사한다. 패러볼릭 부분은 방사선 확산과 흡수·방출(source‑sink) 항을 포함하며, 강직성(stiffness) 때문에 후방 오일러(backward‑Euler) 방식으로 암시적 해결한다. 플럭스 제한자 λ는 Levermore‑Pomraning 형태를 채택해 광학 얇음·두꺼움 양극을 매끄럽게 연결한다. 또한, 복사량을 라보 프레임에서, 흡수·산란 계수를 코무빙 프레임에서 취함으로써 전체 에너지 보존을 보장한다. AMR 구현에서는 각 레벨에서 위 두 단계가 순차적으로 수행되며, 레벨 간 시간 동기화는 서브사이클링을 통해 이루어진다. 코드 스케일링 테스트는 수천 코어까지 거의 선형 효율을 보이며, 방사선‑유체 결합이 강한 초신성 폭발, 별 형성, 원반 물질 흐름 등 다양한 천체 물리 현상에 적용 가능함을 시연한다. 전체 알고리즘은 기존 CASTRO 수소·중력·핵반응 모듈과 완전히 호환되며, 복사 전이와 유체 역학을 동시에 다루는 고성능 시뮬레이션에 필요한 정확도와 안정성을 제공한다.