적응형 격자와 플럭스 제한 확산을 이용한 방사선 유체역학 솔버 개발 및 원시성운 핵 붕괴 적용

초록

본 논문은 RAMSES 코드에 플럭스 제한 확산(FLD)과 회색 근사를 적용한 방사선-유체역학(RHD) 솔버를 구현하고, 이를 적응형 격자(AMR)와 결합하여 1 M⊙ 원시성운 핵의 3차원 붕괴 시뮬레이션에 적용한 결과를 제시한다. 명시적 고전적 유체 부분과 암시적 방사선 확산·물질-방사선 결합을 시간 분할 방식으로 처리하여 2차 정확도를 달성하였다.

상세 분석

본 연구는 방사선 전달이 별 형성 과정에 미치는 영향을 정량적으로 평가하기 위해, RAMSES 코드에 새로운 RHD 모듈을 삽입하였다. 방사선 양은 코무빙 프레임에서 회색(주파수 적분) 형태로 기술되며, 방사선 에너지 방정식은 플럭스 제한 확산(FLD) 근사를 통해 확산 방정식 형태로 변환된다. FLD에서 사용된 Minerbo 제한자는 광학 두께가 큰 영역에서는 λ≈1/3, 투명 영역에서는 λ≈1/R을 만족하여, 광학적 전이 영역에서도 안정적인 플럭스를 제공한다.

수치적 구현은 두 단계로 나뉜다. 첫 번째는 RAMSES의 기존 2차 MUSCL Godunov 스킴을 활용한 명시적 단계로, 물질의 질량·운동량·전체 에너지 보존 방정식을 방사선 압력(Pr=Er/3)까지 포함해 수정한다. 이때 고유 파동속도는 기존 음속에 방사선 압력 기여를 더해 λγP/ρ+4Er/9ρ 형태로 확장되며, CFL 조건도 이에 맞게 조정된다. 두 번째는 방사선 확산항과 물질-방사선 흡·복사 결합항을 암시적으로 풀기 위한 단계이다. 방정식의 비선형성(특히 T⁴ 항)을 선형화하기 위해 온도 변화를 작은 것으로 가정하고 1차 테일러 전개를 적용, 결과적으로 각 셀에서 온도와 방사선 에너지 사이의 선형 연립식을 얻는다. 이 연립식은 전역 행렬 A·x=b 형태로 구성되고, Conjugate Gradient 알고리즘을 이용해 효율적으로 해결한다.

AMR 구조와의 결합에서는, 각 정교화 레벨마다 동일한 암시적 스키마를 적용하되, 경계 셀에서 플럭스 제한자 λ와 Rosseland·Planck 불투명도를 현재 단계의 온도로 평가해 선형성을 유지한다. 또한, 방사선 압력 항을 명시적 단계에 포함시키면서도, 방사선 에너지 흐름은 암시적 단계에서만 업데이트함으로써 시간 스케일 차이를 효과적으로 조정한다.

검증 테스트로는 방사선 충격파, 방사선 전파, 그리고 방사선-유체 결합 문제 등을 수행했으며, 모두 기존 문헌 결과와 일치함을 확인했다. 최종 3차원 무회전 1 M⊙ 핵 붕괴 시뮬레이션에서는, 방사선에 의한 온도 상승과 첫 번째 코어 형성 시점의 질량·반지름이 이전 연구와 유사하게 재현되었다. 이는 제안된 RHD 솔버가 별 형성 시뮬레이션에 필요한 물리적 정밀도와 AMR 기반 고해상도 계산을 동시에 만족한다는 강력한 증거가 된다.