[논문 리뷰] Heat Flow and Boundary Heterogeneity in Rotating Convection

J. E. Mound, Christopher J. Davies|arXiv (Cornell University)|2017. 08. 31.

Geomagnetism and Paleomagnetism Studies참고 문헌 69인용 수 11

한 줄 요약

이 연구는 구형 기하학에서 수치 시뮬레이션을 통해 빠르게 회전하는 대류에서 비균일한 열경계 조건이 있는 경우의 열전달을 조사한다. 경계에 의해 유도된 축방향 속도와 온도 비정상 간의 상관관계가 균일한 경우에 비해 나우셀 수(Nusselt number)를 최대 25% 향상시키며, 강한 초임계성과 높은 경계 비균일성 조건에서 규모 지수(γ)는 약 20% 증가한다. 이는 행성 핵의 역학에 중대한 영향을 미친다.

ABSTRACT

Convection in the metallic cores of terrestrial planets is likely to be subjected to lateral variations in heat flux through the outer boundary imposed by creeping flow in the overlying silicate mantles. Boundary anomalies can significantly influence global diagnostics of core convection when the Rayleigh number, $Ra$, is weakly supercritical; however, little is known about the strongly supercritical regime appropriate for planets. We perform numerical simulations of rapidly rotating convection in a spherical shell geometry and impose two patterns of boundary heat flow heterogeneity: a hemispherical $Y_1^1$ spherical harmonic pattern; and one derived from seismic tomography of Earth's lower mantle. We consider Ekman numbers $10^{-4} \le E \le 10^{-6}$, flux-based Rayleigh numbers up to $\sim800$ times critical, and Prandtl number unity. The amplitude of the lateral variation in heat flux is characterised by $q_{ m L}^*=0, 2.3, 5.0$, the peak-to-peak amplitude of the outer boundary heat flux divided by its mean. We find that the Nusselt number, $Nu$, can be increased by up to $\sim25$\% relative to the equivalent homogeneous case due to boundary-induced correlations between the radial velocity and temperature anomalies near the top of the shell. The $Nu$ enhancement tends to become greater as the amplitude and length scale of the boundary heterogeneity are increased and as the system becomes more supercritical. This $Ra$ dependence can steepen the $Nu \propto Ra^γ$ scaling in the rotationally dominated regime, with $γ$ for our most extreme case approximately $20\%$ greater than the equivalent homogeneous scaling. Therefore, it may be important to consider boundary heterogeneity when extrapolating numerical results to planetary conditions.

연구 동기 및 목표

• 핵-매니플라 경계에서 열역학적 열류의 횡방향 변동이 행성 핵 내 열전달에 미치는 영향을 이해하기 위해.
• 지구형 행성에 관련된 강한 초임계성, 빠르게 회전하는 조건에서 경계 비균일성이 대류 열전달에 미치는 영향을 조사하기 위해.
• 열경계 비정상이 자전에 의해 지배되는 대류에서 나우셀 수 규모 지수 γ에 미치는 영향을 정량화하기 위해.
• 수치 모델에서 일반적으로 사용되는 균일한 경계 조건이 행성 내부의 열전달 효율성을 잘못 표현하는지 평가하기 위해.
• 경계 조건의 진폭과 길이 척도가 대류 역학과 열전달 규모 지수를 어떻게 변화시키는지 평가하기 위해.

제안 방법

• 에카먼 수(Ekman number)가 10⁻⁶에서 10⁻⁴ 사이인 구형 셸 기하학에서 빠르게 회전하는 대류의 수치 시뮬레이션.
• 경계 열류 비균일성에 대해 두 가지 패턴을 적용: 반구형 Y₁¹ 구면 조화함수 및 지구의 내부 맨틀 지질학적 탐사 데이터에서 유도된 패턴.
• 열류 기반 레이놀즈 수(Rayleigh number)를 약 800배 임계값까지 시뮬레이션하며, 프란틀 수(Prandtl number)는 일정하게 유지함.
• 횡방향 열류 변동의 진폭은 q*ₗ = 0, 2.3, 5.0로 표현되며, 평균 대비 최대-최소 변동 비율을 의미함.
• 열전달 및 유동 역학 분석을 위해 나우셀 수(Nu), 레이놀즈 수, 운동에너지 및 열에너지 소산률을 계산함.
• E, Ra 및 경계 비균일성의 변화에 따라 Nu ∝ Ra^γ 관계의 규모 관계를 분석하여 균일한 규모와의 편차 평가함.

실험 결과

연구 질문

• RQ1 횡방향 비균일성이 있는 열경계 조건이 빠르게 회전하는 대류에서 나우셀 수에 어떤 영향을 미치는가?
• RQ2 경계에 의해 유도된 축방향 속도와 온도 비정상 간의 상관관계가 열전달 효율을 얼마나 향상시키는가?
• RQ3 비균일한 경계 조건에서 나우셀 수 규모 관계 Nu ∝ Ra^γ의 규모 지수 γ는 균일한 경우에 비해 어떻게 변화하는가?
• RQ4 열전달 향상 정도가 경계 비균일성의 진폭과 길이 척도에 따라 달라지는가?
• RQ5 레이놀즈 수가 강한 초임계성 영역으로 증가함에 따라 시스템 행동은 어떻게 변화하는가?

주요 결과

• 경계에 의해 유도된 축방향 속도와 온도 비정상 간의 상관관계로 인해 상부 경계 근처에서 나우셀 수(Nu)가 균일한 경우에 비해 최대 25% 향상됨.
• 비균일성의 진폭(q*ₗ)과 길이 척도가 증가할수록 Nu 향상 효과가 더욱 두드러지며, 특히 고레이놀즈 수 조건에서 두드러짐.
• 가장 극단적인 경우(q*ₗ = 5.0)에서 Nu ∝ Ra^γ 관계의 규모 지수 γ는 균일한 경우에 비해 약 20% 증가하여 Ra에 대한 의존성이 더욱 뚜렷해짐.
• 이 효과는 특히 E ≤ 10⁻⁵인 자전에 의해 지배되는 영역에서 가장 두드러지며, 이 경우 시스템은 강력한 초임계성 상태를 유지하고 경계층 역학이 크게 변화함.
• 결과적으로 경계 비균일성이 열전달 규모 지수를 상당히 변화시킬 수 있으며, 수치 결과를 행성 조건으로 외삽할 때 균일한 경계 조건 가정의 타당성을 의심스럽게 만든다.
• 관측된 Nu 향상은 열비정상과 정렬된 지속적인 축방향 유동 구조와 관련이 있으며, 이는 상부 경계 영역에서 효율적인 열 이송을 촉진함.