[논문 리뷰] Thermonuclear explosions of Chandrasekhar-mass C+O white dwarfs

M. Reinecke, W. Hillebrandt외 1명|ArXiv.org|1998. 12. 07.

Gamma-ray bursts and supernovae참고 문헌 7인용 수 7

한 줄 요약

이 논문은 수치적 난류 연소 및 레벨 세트 방법을 통해 난류flame 전파와 날카운 연소 front를 모델링하는 새로운 유한체적 코드를 사용하여 찬드라세카르 질량의 C+O 백색왜성에 대한 2차원 유체역학 시뮬레이션을 제시한다. 연소 구조와 난류를 정확하게 처리함에도 불구하고 중심점화 시뮬레이션은 여전히 유한한 상태에 머물러 있으며, 비중앙점화는 일반적인 Ia 초신성 에너지에 못 미치는 약한 폭발을 일으킨다.

ABSTRACT

First results of simulations are presented which compute the dynamical evolution of a Chandrasekhar-mass white dwarf, consisting of equal amounts of carbon and oxygen, from the onset of violent thermonuclear burning, by means of a new two-dimensional numerical code. Since in the interior of such a massive white dwarf nuclear burning progresses on microscopic scales as a sharp discontinuity, a so-called flamelet, which cannot be resolved by any numerical scheme, and since on macroscopic scales the burning front propagates due to turbulence, we make an attempt to model both effects explicitly in the framework of a finite-volume hydrodynamics code. Turbulence is included by a sub-grid model, following the spirit of large eddy simulations, and the well-localized burning front is treated by means of a level set, which allows us to compute the geometrical structure of the front more accurately than with previous methods. The only free parameters of our simulations are the location and the amount of nuclear fuel that is ignited as an initial perturbation. We find that models in which explosive carbon burning is ignited at the center remain bound by the time the front reaches low densities, where we stopped the computations because our description of combustion is no longer applicable. In contrast, off-center ignition models give rise to explosions which, however, are still too weak for typical Type Ia supernovae. Possible reasons for this rather disappointing result are discussed.

연구 동기 및 목표

• 폭발적인 열핵반응의 시작부터 찬드라세카르 질량의 C+O 백색왜성의 역학적 진화를 시뮬레이션하는 것.
• 수치적 방법을 사용하여 거대한 백색왜성에서 미세구조의flame 전면과 난류 연소를 해결하는 데 도전하는 것.
• 중심점화 또는 비중앙점화가 Ia 초신성과 일치하는 폭발적 붕괴를 유도하는지 판단하는 것.
• 난류와flame 전면 기하학이 폭발 에너지의 결정에 미치는 역할을 평가하는 것.

제안 방법

• 점화 이후 백색왜성의 전체 역학적 진화를 모델링하기 위해 새로운 2차원 유한체적 유체역학 코드를 사용한다.
• 이전 방법보다 더 높은 기하학적 정확도로 날카우며 잘 국소화된 연소 전면을 추적하기 위해 레벨 세트 방법을 사용한다.
• 해상도가 떨어지는 난류를 고려하기 위해 대규모 난류 시뮬레이션(Large Eddy Simulation, LES) 원리에 기반한 미세구조 스케일 모델을 구현한다.
• 초기 조건은 핵분열 연료의 점화 위치와 양만 지정되며, 연소 파rameter에 대한 추가 조정은 이루어지지 않는다.
• 연소 전면이 낮은 밀도에 도달할 경우 모델이 적용 불가가 되므로 시뮬레이션은 그 지점에서 중단된다.
• 결과를 알려진 물리적 제약과 비교하고 질량 및 에너지 보존을 확보함으로써 코드를 검증한다.

실험 결과

연구 질문

• RQ1 2차원 유체역학 시뮬레이션은 찬드라세카르 질량의 C+O 백색왜성에서 연소에서 폭발로의 전이를 정확하게 모델링할 수 있는가?
• RQ2 초기 점화 위치(중심 대비 비중앙)는 폭발 에너지와 잔여물의 유한성에 어떤 영향을 미치는가?
• RQ3 백색왜성 폭발의 맥락에서 미세구조 및 레벨 세트 기법을 사용하여 난류와flame 전면 기하학을 어느 정도 정확하게 모델링할 수 있는가?
• RQ4 중심점화 시뮬레이션은 폭발적 연소가 빈약한데도 불구하고 왜 유한한 잔여물을 남기나?
• RQ5 현재 수치 모델은 관측된 Ia 초신성의 에너지 특성을 재현하는 데 어떤 한계를 지니는가?

주요 결과

• 중심점화 시뮬레이션은 모델의 저밀도 절단 시점에 도달하기 전에 별을 해체할 수 있는 충분한 에너지를 얻지 못하여 유한한 잔여물을 낳는다.
• 비중앙점화는 폭발적 붕괴를 일으키지만, 결과 폭발 에너지는 일반적인 Ia 초신성 관측치에 못 미친다.
• 레벨 세트 방법은 이전 방법보다flame 전면의 기하학적 구조를 더 정확하게 표현하며, 난류 연소에서 전면 추적 성능을 향상시킨다.
• 미세구조 난류 모델은flame 전파에 영향을 미치는 대규모 난류 효과를 포착하지만, 전체 난류의 단순화된 근사치에 불과하다.
• 시뮬레이션은 여전히 유한하거나 약하게 유한한 상태를 유지하며, 현재 모델링 접근법이 필요한 폭발 에너지를 확보하지 못하고 있음을 시사한다.
• 저자들은 낮은 밀도에서의 해상도가 떨어지는flame 역학과 불완전한 연소 모델링이 폭발 강도가 과소평가되는 주요 원인로 지목한다.