[논문 리뷰] Radiation hydrodynamics with neutrinos: Variable Eddington factor method for core-collapse supernova simulations

Markus Rampp, Hans‐Thomas Janka|ArXiv.org|2002. 03. 07.

Neutrino Physics Research참고 문헌 86인용 수 213

한 줄 요약

이 논문은 핵붕괴 초신성에서 중성미자 운반을 시뮬레이션하기 위한 새로운 수치 코드를 제시한다. 시간 및 에너지 의존성 버틀로프 방정식을 해결하기 위해 가변 에딩턴 인자 방법을 사용하며, 시간 암시적 통합을 이용한 레이 격자 상에서의 모멘트 방정식과 모델 버틀로프 방정식을 반복적으로 풀어, 중성미자 에너지 및 렙톤 수의 교환을 정확하게 모델링할 수 있다. 모듈러 설계 덕분에 다차원 및 상대론적 응용에 적합하다.

ABSTRACT

Neutrino transport and neutrino interactions in dense matter play a crucial role in stellar core collapse, supernova explosions and neutron star formation. Here we present a detailed description of a new numerical code for treating the time and energy dependent neutrino transport in hydrodynamical simulations of such events. The code is based on a variable Eddington factor method to deal with the integro-differential character of the Boltzmann equation. The moments of the neutrino distribution function and the energy and lepton number exchange with the stellar medium are determined by iteratively solving the zeroth and first order moment equations in combination with a model Boltzmann equation. The latter is discretized on a grid of tangent rays. The integration of the transport equations and the neutrino source terms is performed in a time-implicit way. In the present version of the program, the transport part is coupled to an explicit hydrodynamics code which follows the evolution of the stellar plasma by a finite-volume method with piecewise parabolic interpolation, using a Riemann solver for calculating the hydrodynamic states. The neutrino source terms are implemented in an operator-split step. Neutrino transport and hydrodynamics can be calculated with different spatial grids and different time steppings. The structure of the described code is modular and offers a high degree of flexibility for an application to relativistic and multi-dimensional problems at different levels of refinement and accuracy. We critically evaluate results for a number of test cases, including neutrino transport in rapidly moving stellar media and approximate relativistic core collapse, and suggest a path for generalizing the code to be used in multi-dimensional simulations of convection in neutron stars and supernovae.

연구 동기 및 목표

• 핵붕괴 초신성에서 중성미자 운반을 고정밀도로 시뮬레이션할 수 있는 견고하고 영리한 수치 프레임워크를 개발하는 것.
• 다차원, 시간 의존성 유체역학 환경에서 적분미분 버틀로프 방정식을 해결하는 데 발생하는 계산적 과제를 해결하는 것.
• 핵붕괴 및 충격 복구 기간 동안 밀도 높은 항성 물질과의 중성미자 에너지 및 렙톤 수 교환을 정확하게 모델링할 수 있도록 하는 것.
• 미래의 다차원 및 상대론적 초신성 시뮬레이션에 적합한 모듈러이고 확장 가능한 코드 구조를 제공하는 것.
• 빠른 유동과 근사 상대론적 붕괴를 포함한 테스트 케이스를 통해 방법의 신뢰성을 검증하는 것.

제안 방법

• 중성미자 분포 함수의 각도 의존성을 근사하기 위해 가변 에딩턴 인자 방법을 사용하여 전체 버틀로프 방정식을 모멘트 방정식의 집합으로 축소한다.
• 영차 및 제1차 모멘트 방정식을 반복적으로 풀며, 탄젠트 레이 격자에 이산화된 모델 버틀로프 방정식과 결합한다.
• 강성 있고 급격히 변화하는 시스템에서 안정성을 확보하기 위해 운반 및 소스 항정 방정식에 시간 암시적 통합을 적용한다.
• 조각별 포arus 보간과 리만 해법을 사용한 유체역학 상태를 활용해 중성미자 운반을 명시적 유한부피 유체역학 코드와 결합한다.
• 운반 소스 항정을 연산자 분할을 통해 구현하여 유체역학과 중성미자 운반에 대해 별개의 공간 격자와 시간 간격을 허용한다.
• 핵 구성 전이를 밀도-온도 영역 3개로 나누어 모델링하며, 해리/재결합 요소 $f_{\mathrm{I}}, f_{\mathrm{II}}, f_{\mathrm{III}}$ 로 연속적인 전이를 제어한다.

실험 결과

연구 질문

• RQ1 핵붕괴 초신성 시뮬레이션에서 시간 및 에너지 의존성의 적분미분 버틀로프 방정식을 효율적이고 정확하게 해결할 수 있는가?
• RQ2 반복적인 모멘트 해법과 레이 기반 이산화를 갖는 가변 에딩턴 인자 방법이 중성미자 히팅 및 에너지 침착을 모델링하는 데 충분한 정확도를 확보할 수 있는가?
• RQ3 코드의 모듈러 설계는 향후 다차원 및 상대론적 시뮬레이션으로의 확장에 어떻게 기여하는가?
• RQ4 예를 들어 $^{56}$Ni 형성, $\alpha$-입자 재결합, 무거운 핵종 해리 등 핵 구성 전이가 붕괴 기간 동안 에너지 및 렙톤 운반에 어떤 영향을 미치는가?
• RQ5 빠르게 움직이는 매체 내 중성미자 운반 또는 근사 상대론적 붕괴와 같은 알려진 테스트 케이스와의 비교에서 수치 결과는 어떻게 되는가?

주요 결과

• 가변 에딩턴 인자 방법은 핵붕괴 초신성에서 중성미자 운반의 핵심 물리적 메커니즘, 즉 중성미자 에너지 및 렙톤 수의 항성 매질과의 교환을 성공적으로 포착한다.
• 시간 암시적 통합 기법은 핵붕괴 시점과 충격 형성과 같은 극도로 동적인 강성 환경에서도 수치적 안정성을 보장한다.
• 코드의 모듈러 아키텍처는 유체역학과 중성미자 운반에 대해 별개의 공간 및 시간 이산화를 가능하게 하여 계산의 영리성과 효율성을 향상시킨다.
• 빠르게 움직이는 매체 내 중성미자 운반과 근사 상대론적 붕괴를 위한 시뮬레이션 테스트 결과는 기대되는 물리적 거동과 양호하게 일치하며, 방법의 정확성을 검증한다.
• 핵 구성 전이를 $f_{\mathrm{I}}, f_{\mathrm{II}}, f_{\mathrm{III}}$ 요소를 통해 처리함으로써 플라즈마의 상 변화 기간 동안 에너지 보존과 열역학적 일관성이 유지된다.
• 이 방법은 중성별 및 초신성에서 대류와 중성미자 구동 폭발 메커니즘의 다차원 시뮬레이션에 향후 적용하기에 매우 적합하다.