[논문 리뷰] Three-Dimensional Simulations of Classical Novae

A. Kercek, W. Hillebrandt외 1명|arXiv (Cornell University)|1998. 11. 17.

Astrophysical Phenomena and Observations참고 문헌 1인용 수 23

한 줄 요약

이 연구는 3D 유체역학 시뮬레이션을 통해 고전형 초신성 폭발의 첫 번째 결과를 제시한다. 3D 오일러형 코드를 사용하여 카르본+산소 백색왜성(1 M⊙)에서의 난류 수소 연소를 모의하기 위해 카르테시안 격자에 핵반응 네트워크를 적용하였다. 주요 발견은 폭발 기간 동안 축적된 대기의 탄소와 산소로의 자가풍부화가 너무 느리며, 심지어 태양 금속성의 5배일지라도 빠른 초신성의 형성을 유도하지 못한다는 것이다. 이는 빠른 초신성 형성에 있어 폭발 이전의 풍부화 메커니즘이 필수적임을 시사한다.

ABSTRACT

We present first results of three-dimensional (3D-) calculations of turbulent and degenerate hydrogen-burning on top of a C+O white dwarf of one solar mass. The simulations are carried out by means of a code which solves Euler's equation for an arbitrary equation of state together with a nuclear reaction network and the energy input from nuclear reactions on a Cartesian grid covering a fraction of the white dwarf's surface and accreted atmosphere. The flow patterns we obtain are very different from those of earlier 2D simulations using the same initial conditions and the same numerical resolution. The possibility of self-enrichment of the accreted hydrogen-rich atmosphere with carbon and oxygen from the surface layers of the white dwarf during the violent phase of the burning is investigated, and it is demonstrated that self-enrichment proceeds too slowly if the accreted gas has near-solar CNO-abundances at the onset of the thermonuclear runaway. As a result, we do not find a fast nova outburst. This conclusion remains valid if the initial metallicity of the accreted gas is raised by a factor of five. Therefore we conclude that fast nova outbursts indeed require huge enrichments of C and O, as postulated from spherically symmetric models, and that the mechanism which leads to such enhancements must operate prior to the outburst.

연구 동기 및 목표

• 고전형 초신성 폭발에서 3D 유체역학의 역할을 조사하며, 특히 열핵반응기 동안 난류 대류 및 혼합의 역학을 다루는 것.
• 백색왜성 표면의 탄소와 산소로 축적된 수소 풍부 대기가 폭발 기간 동안 자가풍부화되어 빠른 초신성 폭발을 유도할 수 있는지 테스트하는 것.
• 3D 결과를 이전의 2D 시뮬레이션과 비교하여 차원 수의 영향을 혼합, 에너지 생성 및 폭발 역학에 미치는 영향을 평가하는 것.
• 초기 금속성의 향상(최대 5배 태양 금속성)이 연소 과정을 충분히 가속화하여 빠른 초신성 특성을 유도할 수 있는지 평가하는 것.
• 3D 시뮬레이션에서 수치적 확산 또는 해상도가 부족한 소규모 운동이 2D에 비해 혼합을 더 억제하는지 평가하여 폭발적 폭발의 가능성에 영향을 미치는지 확인하는 것.

제안 방법

• 수치 시뮬레이션은 3D 카르테시안 격자에서 유체역학 방정식을 해석하는 PPM 유형의 수치 해법 코드인 PROMETHEUS를 사용하여 수행되었으며, 백색왜성 표면과 대기의 일부를 커버하였다.
• 유체역학 방정식은 12종의 입자 종류(1H, 4He, 12C, 13C, 13N, 14N, 15N, 14O, 15O, 16O, 17O, 17F)를 포함한 핵반응 네트워크와 결합하여 해석되었으며, 반응률은 Wallace & Woosley(1981) 및 Thielemann(비공개 통신)의 자료를 사용하였다.
• 에너지 생성은 네트워크 방정식과 동시에 해석하여 수치적 불안정성을 방지하였으며, 일관된 상태방정식과 암시적 시간 적분을 사용하였다.
• 초기 조건은 축적과 폭발 이전 진화를 모의한 1D 암시적 유체역학 모델(Glasner 등 1997)에서 유도되었으며, 1.0 M⊙의 C+O 백색왜성과 Z=0.02(일부 시뮬레이션에서는 5× 태양 금속성)의 축적 대기(2×10⁻⁵ M⊙)를 가진 상태였다.
• 3D 격자는 1D 모델에서 유도된 수압 평형 상태로 초기화되었으며, 수십 개의 동역학 시간 스케일 동안 안정화된 후, 대류와 연소를 유도하기 위해 흔들림을 가하였다.
• 비물리적 반사 현상을 방지하기 위해 주기적 경계 조건을 적용하였으며, 안정성과 적분 양의 수렴성을 확보하기 위해 중간 해상도(예: 220×100 격점)로 시뮬레이션를 수행하였다.

실험 결과

연구 질문

• RQ1 고전형 초신성 폭발 중 3D 유체역학 흐름은 2D 시뮬레이션과 어떻게 다를까? 특히 코리올리력에 의한 구조와 혼합 효율 측면에서 비교해보라.
• RQ2 백색왜성 표면의 탄소와 산소로 축적된 수소 풍부 대기가 열핵반응기 동안 충분히 빠르게 자가풍부화되어 빠른 초신성 폭발을 유도할 수 있는가?
• RQ3 축적 기체의 초기 금속성(최대 5× 태양 금속성)을 높이면 연소 과정이 충분히 가속화되어 빠른 초신성 폭발을 유도할 수 있는가?
• RQ4 3D 시뮬레이션은 2D 모델에 비해 대류의 과도 및 C, O의 끌어올림을 얼마나 억제하는가? 이는 에너지 방출과 폭발 역학에 어떤 영향을 미치는가?
• RQ5 3D 시뮬레이션에서의 수치적 확산이 혼합을 이끌어내기에 충분한가, 아니면 정확한 혼합을 포괄하기 위해 명시적 미세구역 난류 모델링이 필요한가?

주요 결과

• 3D 시뮬레이션은 2D 시뮬레이션과는 다소 다른 유동 패턴을 생성하며, 소규모 난류 운동의 에너지가 더 많아져 대류 과도 현상과 C, O의 혼합 억제를 유도한다.
• 폭발 기간 동안 축적된 대기의 탄소와 산소로의 자가풍부화가 너무 느리며, 심지어 초기 금속성이 5× 태양 금속성으로 증가하더라도 빠른 초신성을 유도하지 못한다. 이는 온도 상승 시간 스케일이 여전히 약 100초 수준으로 너무 길기 때문이다.
• 고금속성 시뮬레이션에서 최고 속도는 150초 후 약 2×10⁷ cm/s에 도달하였으며, 저금속성 케이스에 비해 약 두 배 수준이지만, 에너지 생성률과 최종 대기 금속성은 다소 높을 뿐이다.
• 고금속성 케이스는 400초에 도달했을 때 저금속성 케이스와 구분되지 않는 연소 상태에 도달하여, 초기 금속성의 증가가 근본적인 시간 스케일 문제를 해결하지 못함을 시사한다.
• 3D 시뮬레이션 결과, 혼합 감소로 인해 폭발적 수소 연소가 약화되어 2D 시뮬레이션보다 폭발이 덜 격렬하며, 따라서 자가풍부화에 의한 빠른 초신성 형성 가능성이 낮아진다.
• 결과적으로, 폭발 중 자가풍부화는 빠른 초신성 형성의 타당한 메커니즘이 아니며, 대신 큰 양의 C와 O 풍부화가 폭발 이전에 발생해야 하며, 이는 아마도 난류 불안정성 또는 축적 단계 동안의 장기 혼합에 의해 발생할 것이다.