[논문 리뷰] Hydrodynamic simulations of H entrainment at the top of He-shell flash convection

Paul R. Woodward, Falk Herwig외 1명|arXiv (Cornell University)|2013. 07. 15.

Stellar, planetary, and galactic studies참고 문헌 2인용 수 3

한 줄 요약

이 연구는 AGB 별의 헬륨 쉘 플래시 대류 영역 상단에서 수소가 얽히는 현상을 $4\pi$ 기하학에서 3차원, 완전 압축성 유체역학 시뮬레이션으로 처음으로 수행하였으며, 고해상도이고 안정적인 시뮬레이션을 달성하기 위해 PPM+PPB 복합유체 운반 기법을 사용하였다. 주요 결과는 $1536^3$ 해상도에서 수렴된 양피성율 $4.38 \pm 1.48 \times 10^{-13}\,M_{\odot}/\text{s}$를 도출하였으며, 이는 항성 내부의 딱딱한 대류 경계에서의 수치적 정밀도를 입증한다.

ABSTRACT

We present the first 3-dimensional, fully compressible gas-dynamics simulations in $4π$ geometry of He-shell flash convection with proton-rich fuel entrainment at the upper boundary. This work is motivated by the insufficiently understood observed consequences of the H-ingestion flash in post-AGB stars (Sakurai's object) and metal-poor AGB stars. Our investigation is focused on the entrainment process at the top convection boundary and on the subsequent advection of H-rich material into deeper layers, and we therefore ignore the burning of the proton-rich fuel in this study. We find that, for our deep convection zone, coherent convective motions of near global scale appear to dominate the flow. At the top boundary convective shear flows are stable against Kelvin-Helmholtz instabilities. However, such shear instabilities are induced by the boundary-layer separation in large-scale, opposing flows. This links the global nature of thick shell convection with the entrainment process. We establish the quantitative dependence of the entrainment rate on grid resolution. With our numerical technique simulations with $1024^3$ cells or more are required to reach a numerical fidelity appropriate for this problem. However, only the result from the $1536^3$ simulation provides a clear indication that we approach convergence with regard to the entrainment rate. Our results demonstrate that our method, which is described in detail, can provide quantitative results related to entrainment and convective boundary mixing in deep stellar interior environments with veryvstiff convective boundaries. For the representative case we study in detail, we find an entrainment rate of $4.38 \pm 1.48 imes 10^{-13}M_\odot \mathrm{/s}$.

연구 동기 및 목표

• AGB 별의 헬륨 쉘 플래시 대류 영역 상단에서 양성자-rich 연료가 얽히는 현상을 연구하여, 사쿠라이의 물체와 같은 후-AGB 별에서의 수소 흡입 플래시와 연관된 과정을 밝히는 것.
• 표준 혼합길이 모델이 실패하는 깊고 딱딱한 대류 경계 환경에서의 정량적 양피성율을 결정하는 것.
• 대류 경계 혼합을 위한 3차원 유체역학 시뮬레이션의 수치적 정밀도를 검증하여 해상도 수렴성과 안정성에 중점을 두는 것.
• PPM+PPB 복합유체 운반 기법이 고해상도에서 압축성이고 난류적인 흐름에서 날카로운 유체 경계를 정확하고 낮은 확산으로 추적할 수 있음을 입증하는 것.

제안 방법

• 헬륨 쉘 플래시 대류에 수소-rich 연료가 얽히는 현상을 모의하기 위해 $4\pi$ 구형 기하학에서 3차원, 완전 압축성 유체역학 시뮬레이션을 수행한다.
• 혼합 분율 $f$를 사용하는 복합유체 접근법을 적용하여 수소와 헬륨을 추적하며, 각 셀 내에서 압력과 온도 평형을 가정하여 각 유체의 별도 밀도 저장을 피한다.
• 혼합 분율 $f$를 세포 경계를 넘어서 고차수이고 안정적인 운반을 달성하기 위해 PPM+PPB(조각별 파라볼라 방법 + 파라볼라 프로파일 구속) 복원 기법을 적용한다.
• 유체 밀도를 혼합 밀도와 $f$에서 평균 분자량 비율을 사용해 보간하여 엄격한 질량 보존을 확보하며, 개별 유체 밀도의 직접 저장을 피한다.
• 해상도 수렴성과 수치적 정밀도의 임계 해상도를 평가하기 위해 여러 해상도($1024^3$에서 $1536^3$)로 시뮬레이션을 수행한다.
• 주요 결과로 고해상도 $1536^3$ 시뮬레이션을 사용하며, 그리드 정밀도 향상과 양피성율의 안정성에 기반한 수렴 분석을 수행한다.

실험 결과

연구 질문

• RQ1 헬륨 쉘 플래시 대류 영역 상단에서 수소-rich 물질의 양피성율은 얼마이며, 이는 수치적 해상도에 어떻게 의존하는가?
• RQ2 대규모 대류 운동과 유도 흐름이 딱딱하고 깊은 항성 대류 영역의 대류 경계의 안정성과 역학에 어떤 영향을 미치는가?
• RQ3 PPM+PPB 복합유체 운반 기법이 고해상도에서 압축성이고 난류적인 항성 대류에서 날카로운 유체 경계를 안정적이고 낮은 확산으로 시뮬레이션할 수 있는가?
• RQ4 두꺼운 셸 대류에서 전역적, 거의 전역적 수준의 대류 흐름의 구조가 양피성 과정을 얼마나 지배하는가?
• RQ5 해상도가 증가함에 따라 수치 시뮬레이션이 안정된 양피성율로 수렴하는가? 그리고 정량적 신뢰성에 필요한 해상도는 무엇인가?

주요 결과

• $1536^3$ 시뮬레이션에서 유일하게 수렴의 명확한 징후를 보이며, 양피성율이 정량적으로 $4.38 \pm 1.48 \times 10^{-13}\,M_{\odot}/\text{s}$로 결정되었다.
• $1024^3$ 셀 이하의 시뮬레이션은 충분한 수치적 정밀도를 확보하지 못하며, 오직 $1536^3$에서만 양피성율의 신뢰할 수 있는 추정치를 제공한다.
• 일관되고 거의 전역 수준의 대류 운동이 흐름을 지배하며, 경계층 분리 없이 경계에서의 셰어 흐름은 켈빈-헬름홀츠 불안정성에 대해 안정되어 있다.
• 대규모 상반된 흐름에서의 경계층 분리가 켈빈-헬름홀츠 불안정성을 유도하며, 이는 전역 대류 역학이 양피성 과정과 연결됨을 시사한다.
• PPM+PPB 방법은 유체 경계를 통해 단지 2~3개의 격자 셀 두께로도 매끄럽고 고해상도의 혼합 분율 $f$ 추적을 가능하게 하며, 수치 진동과 확산을 피한다.
• 국소적 압력 및 온도 평형 가정 덕분에 혼합 밀도와 $f$에서 정확한 밀도 재구성을 가능하게 하여 계산 비용을 크게 감소시키면서도 정확도를 유지한다.