[논문 리뷰] A Mathematical Description of the IDSA for Supernova Neutrino transport, its discretization and a comparison with a finite volume scheme for Boltzmann's Equation

Heiko Berninger, Emmanuel Frénod외 4명|Archive ouverte UNIGE (University of Geneva)|2012. 11. 29.

Numerical methods for differential equations참고 문헌 7인용 수 4

한 줄 요약

이 논문은 핵붕괴 초신성에서 중성미자 운반에 대한 등방성 확산 소스 근사법(IDSA)의 엄밀한 수학적 공식화를 제공하며, 시간 분할 및 유한체적 방법을 사용해 풀어낸 축소된 데볼트츠 방정식과 비교한다. 이는 낮은 이산화 오차에도 불구하고, 고착 및 자유확산 중성미자 영역 간의 전이 영역에서 주로 국소화된 약 20%의 상대적 $L^2$ 노름 오차를 포함하는 상당한 모델링 오차를 유발한다는 것을 입증한다.

ABSTRACT

In this paper we give an introduction to the Boltzmann equation for neutrino transport used in core collapse supernova models as well as a detailed mathematical description of the \emph{Isotropic Diffusion Source Approximation} (IDSA). Furthermore, we present a numerical treatment of a reduced Boltzmann model problem based on time splitting and finite volumes and revise the discretization of the IDSA for this problem. Discretization error studies carried out on the reduced Boltzmann model problem and on the IDSA show that the errors are of order one in both cases. By a numerical example, a detailed comparison of the reduced model and the IDSA is carried out and interpreted. For this example the IDSA modeling error with respect to the reduced Boltzmann model is numerically determined and localized.

연구 동기 및 목표

• 핵붕괴 초신성에서 중성미자 운반에 대한 등방성 확산 소스 근사법(IDSA)의 수학적으로 엄밀한 유도 및 기술 제공.
• 고정 배경 물질을 가진 축소된 데볼트츠 모델을 기반으로 한 시간 분할 및 유한체적 방법을 사용한 수치적 스킴 개발.
• 일致된 비교를 위해 IDSA 이산화를 축소 모델에 재표현 및 적응.
• 수치 실험을 통해 IDSA와 축소된 데볼트츠 모델 간의 모델링 오차를 정량화하고 국소화.
• 특히 고착 및 자유확산 행동 간 전이에서의 핵심 물리적 영역을 포괄하는 데 있어 IDSA의 타당성과 한계 평가.

제안 방법

• 에너지, 반경 및 각도 의존성을 포함한 구형 대칭 조건에서 중성미자 운반에 대한 $Ø(v/c)$ 데볼트츠 방정식 유도.
• 분포 함수를 고착 및 자유확산 성분으로 분해하고, 이를 확산 기반 소스 항을 통해 결합한 IDSA 도입.
• 반응 및 운반 연산자를 분리하기 위해 시간 분할을 적용하고, 반응 부분은 해석적으로, 운반 부분은 보존적인 유한체적 방법으로 풀이.
• 이전 연구에서 유도된 IDSA 이산화를 축소 모델에 적응시켜 공간 및 각도 격자에서의 일致성 확보.
• 이산화 오차가 1% 이하가 되도록 고해상도 격자($N_r=257$, $N_\mu=129$, $\Delta t=1.5625\cdot10^{-3}$) 사용.
• 모델링 오차를 IDSA와 데볼트츠 해의 $L^2$ 상대적 차이로 계산하고, 공간적으로 절대 차이를 맵핑하여 오차 국소화.

Figure 1: Flow map of the flux of neutrinos due to transport and without reactions. We can see that incoming neutrinos ( $\mu<0$ ) turn into outgoing neutrinos ( $\mu>0$ ) and, therefore, eventually leave the computational domain again.

실험 결과

연구 질문

• RQ1 IDSA 모델은 축소된 데볼트츠 방정식에 비해 해의 정확도 측면에서 어떻게 비교되는가?
• RQ2 특히 물리적 영역와 관련하여 IDSA의 모델링 오차는 공간적으로 어디에 가장 뚜렷하게 나타나는가?
• RQ3 IDSA의 모델링 오차 크기는 축소된 데볼트츠 모델에 비해 얼마이며, 전이 영역에 따라 어떻게 달라지는가?
• RQ4 두 모델의 이산화 오차는 비교의 신뢰성에 어떻게 영향을 미치는가?
• RQ5 왜 IDSA는 반응 지배 영역에서 데볼트츠 모델보다 더 빠르게 평형 상태에 도달하는가?

주요 결과

• IDSA와 축소된 데볼트츠 모델 간의 상대적 모델링 오차는 $L^2$ 노름 기준 약 20%이며, 오차 범위는 ±2%이다.
• 모델링 오차는 주로 반응 지배 영역과 자유확산 영역 간의 중간 전이 영역, 특히 공간 간격 $(3,6)$ 에 집중되어 있다.
• 축소된 데볼트츠 모델과 IDSA의 이산화 오차는 모두 1차 수준이지만, 고해상도 격자를 통해 1% 이하로 통제되었다.
• IDSA는 반응 영역과 운반 영역 모두에서 데볼트츠 해를 과소평가하며, 가장 큰 과소평가가 전이 영역에서 발생한다.
• IDSA는 정적 상태 근사에 의해 고착 및 확산 성분 간 강한 결합으로 인해 반응 지배 영역에서 더 빠르게 평형 상태에 도달한다.
• IDSA의 확산 성분은 기대와는 달리 반응 지배 영역에서도 0이 되지 않아 고착 성분이 감소하고, 이는 모델링 오차의 원인이 된다.

(a) Reduced Boltzmann model with $T=0.0005$ : The dependency of the error in the infinity norm on $\Omega$ on the spatial mesh size $h=\Delta r$ is displayed. The curve with round markers represents the relative error, the curve with the plus markers displays the absolute error. The straight line ha