[논문 리뷰] 3D hydrodynamic simulations of massive main-sequence stars. I. Dynamics and mixing of convection and internal gravity waves

Falk Herwig, Paul R. Woodward외 7명|arXiv (Cornell University)|2023. 03. 09.

Stellar, planetary, and galactic studies인용 수 4

한 줄 요약

이 연구는 25 M⊙ 주계열 별의 3차원 유체역학 시뮬레이션을 통해 핵-대두 경계에서 대류 및 내부 중력파(IGWs)에 의한 혼합을 조사한다. 내부 중력파 혼합 효율성에 대한 스케일링 관계 $ D_{\text{IGW-hydro}} \propto L^{4/3} $ 를 확립하였으며, 내부 중력파 혼합이 열적으로 강화된 낙차에 의해 구동될 경우, 정상 상태 가열 조건에서 혼합 계수는 $ 2{-}3 \times 10^4\ \text{cm}^2/\text{s} $ 로 제한되며, 이는 거대 별에서 표면 조성에 미치는 영향이 제한적임을 시사한다.

ABSTRACT

We performed 3D hydrodynamic simulations of the inner $\approx 50\%$ radial extent of a $25\ \mathrm{M_\odot}$ star in the early phase of the main sequence and investigate core convection and internal gravity waves in the core-envelope boundary region. Simulations for different grid resolutions and driving luminosities establish scaling relations to constrain models of mixing for 1D applications. As in previous works, the turbulent mass entrainment rate extrapolated to nominal heating is unrealistically high ($1.58 imes 10^{-4}\ \mathrm{M_\odot/yr}$), which is discussed in terms of the non-equilibrium response of the simulations to the initial stratification. We measure quantitatively the effect of mixing due to internal gravity waves excited by core convection interacting with the boundary in our simulations. The wave power spectral density as a function of frequency and wavelength agrees well with the GYRE eigenmode predictions based on the 1D spherically averaged radial profile. A diffusion coefficient profile that reproduces the spherically averaged abundance distribution evolution is determined for each simulation. Through a combination of eigenmode analysis and scaling relations it is shown that in the $N^2$-peak region, mixing is due to internal gravity waves and follows the scaling relation $D_\mathrm{IGW-hydro} \propto L^{4/3}$ over a $\gtrapprox 2\ \mathrm{dex}$ range of heating factors. Different extrapolations of the mixing efficiency down to nominal heating are discussed. If internal gravity wave mixing is due to thermally-enhanced shear mixing, an upper limit is $D_\mathrm{IGW} \lessapprox 2$ to $3 imes 10^4\ \mathrm{cm^2/s}$ at nominal heating in the $N^2$-peak region above the convective core.

연구 동기 및 목표

• 거대 주계열 별의 핵-대두 경계 근처에서 대류 및 내부 중력파(IGWs)와 같은 혼합 과정을 정량화하는 것.
• 3차원 유체역학 시뮬레이션을 통해 내부 중력파 혼합 효율성에 대한 스케일링 관계를 수립하여 1차원 별 진화 모델에 통합하는 것.
• 시뮬레이션에서 관측된 유동성 침입률의 현실성 평가 및 과도하게 높은 값이 발생하는 원인 규명.
• 내부 중력파 혼합이 낙차에 의해 구동되는지, 또는 유동성 확률적 혼합과 같은 다른 메커니즘에 의해 일어나는지 확인하는 것.
• 3차원 시뮬레이션 결과를 GYRE를 사용한 1차원 고유모드 예측과 비교하기 위해, 구형 평균화된 프로파일을 활용하는 것.

제안 방법

• 25 M⊙ 주계열 별의 내부 50% 반경 범위에 대해 적응 메쉬 강화와 은밀한 시간 적분을 사용한 3차원 유체역학 시뮬레이션을 수행하였다.
• 혼합 효율성에 대한 스케일링 법칙을 유도하기 위해 다양한 격자 해상도와 가열 빛의 강도로 시뮬레이션을 수행하였다.
• 3차원 시뮬레이션의 구형 평균화된 반경 프로파일을 GYRE에 입력하여 내부 중력파의 고유모드를 계산하였다.
• 파워 스펙트럼 밀도를 분석하고 GYRE 예측과 비교하여 파동의 자극 및 전파를 검증하였다.
• 각 시뮬레이션의 평균 농도 변화를 재현할 수 있도록 혼합 계수 프로파일을 유도하였다.
• 고차원 속도 통계량(왜도 및 첨도)을 사용하여 대류 및 파동 지배 영역 간의 동적 경계를 식별하였다.

실험 결과

연구 질문

• RQ1 3차원 시뮬레이션에서 거대 별의 가열 빛의 강도에 따라 내부 중력파 혼합 효율성은 어떻게 스케일링되는가?
• RQ2 왜 3차원 시뮬레이션은 현실적으로 비현실적으로 높은 질량 침입률을 보이며, 이를 1차원 모델과 어떻게 조율할 수 있는가?
• RQ3 N²-피크 영역 내부의 내부 중력파가 혼합에 얼마나 기여하는가? 이는 1차원 고유모드 예측과 비교해 볼 때 어떠한가?
• RQ4 내부 중력파 혼합이 열적으로 강화된 낙차에 의해 구동될 경우, 시뮬레이션 외삽을 통해 정상 상태 가열 조건에서 혼합 효율의 상한선은 무엇인가?
• RQ5 유동성 확률적 혼합과 같은 대체 혼합 메커니즘이 낙차 기반 모델보다 더 높은 내부 중력파 혼합 효율을 설명할 수 있는가?

주요 결과

• 내부 중력파에 의한 혼합 효율성은 가열 빛 강도 범위가 $ \geq 2\ \text{dex} $ 에 걸쳐 $ D_{\text{IGW-hydro}} \propto L^{4/3} $ 의 비율로 스케일링됨.
• 시뮬레이션에서 측정된 파워 스펙트럼 밀도는 구형 평균 1차원 프로파일 기반의 GYRE 고유모드 예측과 잘 일치함.
• 정상 상태 가열 조건으로 외삽한 난류 질량 침입률은 $ 1.58 \times 10^{-4}\ \text{M}_\odot\ \text{yr}^{-1} $ 로, 평형 상태가 아닌 초기 조건으로 인해 비현실적으로 높음.
• 내부 중력파 혼합이 열적으로 강화된 낙차에 의해 구동될 경우, N²-피크 영역에서의 혼합 계수 상한선은 정상 상태 가열 조건에서 $ 2{-}3 \times 10^4\ \text{cm}^2/\text{s} $ 로 제한됨.
• 반경 속도의 확률 밀도 함수(PDF)는 동적 경계에서 왜도와 초과 첨도의 곱이 최대값을 보이며, 이는 대류에서 파동 지배 흐름으로의 전이를 확인함.
• 시뮬레이션 결과는 내부 중력파 혼합이 낙차 기반일 경우 국소적 영향을 미치며, 25 M⊙ 별의 주계열 수명 동안 표면 조성에 큰 변화를 주지 못할 것임을 시사한다.