내부 가열·냉각으로 구현한 회전 구형 껍질 대류의 확산‑자유 역학

초록

본 연구는 구형 껍질 내부에 가열·냉각을 직접 지정함으로써 회전 대류를 모사하고, 열전달·온도 차와 같은 거시적 특성이 확산계수에 의존하지 않는 ‘확산‑자유’ 거동을 보임을 확인하였다. 회전 제약이 강한 경우에도 전단 흐름(제트류)의 구조와 속도는 확산에 민감하지만, 충분히 높은 초임계성에서는 확산‑자유 상태에 근접한다.

상세 분석

이 논문은 Boussinesq 근사하에 회전 구형 껍질 내부에서 내부 가열(q)과 냉각을 동시에 적용하는 새로운 실험 설계를 제시한다. 기존의 경계‑구동 방식은 열전도층을 형성해 전체 열전달을 경계 전도에 의존하게 만들며, 이는 분자 확산계수(ν, κ)의 영향을 크게 받는다. 저자들은 가열·냉각 구역을 고정된 두께 δ로 두고, 자유‑슬립·불투과·절연 경계조건을 적용함으로써 경계 전도층을 완전히 대체하였다.

비차원화는 껍질 두께 d를 길이 스케일, 열 플럭스 기반 온도 스케일 Λ=Fd/κ, 자유‑낙하 시간 τ=p d/Λ를 시간 스케일로 사용한다. 주요 무차원 수는 플럭스 기반 레일리수 Ra_F = d⁴F/(ν κ²), 타일러수 Ta = 4Ω₀² d⁴/ν², 프란틀수 Pr=ν/κ(=1) 및 껍질 반경비 η= r_i/r_o(=0.8)이다. 특히, 내부 가열·냉각에 의해 정의된 ‘플럭스 기반 대류 로스비수’ Ro_cv,F = (Ra_F /Ta³)^{1/2} (dF)^{1/3}/(2Ω₀ d) 은 확산계수에 독립적이며, 실험 설계에서 Ro_cv,F 를 고정하거나 변동시키는 것이 핵심 제어 변수다.

시뮬레이션은 38개의 파라미터 조합으로 수행되었으며, Ro_cv,F 를 0.013.1, Ta 를 10⁵10⁹, Ra_F 를 3×10⁵~9.5×10⁹ 범위에서 탐색하였다. 결과는 세 가지 회전 제약 regime 으로 구분된다. (1) Ro_cv,F < 0.1: 전 스케일이 회전에 강하게 얽혀 ‘회전‑제한’ regime; (2) 0.1 ≤ Ro_cv,F < 1: 큰 스케일은 회전에 얽히지만 작은 스케일은 자유로운 ‘회전‑영향’ regime; (3) Ro_cv,F ≥ 1: 회전 효과가 거의 사라진 ‘비회전’ regime.

핵심 발견은 다음과 같다. 첫째, 방사형 온도 차 ΔT와 대류 열전달(Nu) 은 Ro_cv,F 에만 의존하고, 동일 Ro_cv,F 에서는 ν, κ 의 변화에도 데이터가 거의 겹쳐 ‘확산‑자유’ 스케일링을 보인다. 이는 ΔT ∝ Ro_cv,F^{-α} (α≈0.50.7) 형태로, 회전이 강할수록 열전달이 억제됨을 의미한다. 둘째, 전단 흐름(적도 대류대류 흐름)의 방향 전환은 Ro_cv,F ≈ 0.30.5 부근에서 일어나며, 이 전환 역시 확산계수와 무관하게 Rossby 수만으로 예측된다. 셋째, 평균 레이놀즈수 Re와 소규모 와류의 로스비수 Ro_ω 은 확산에 민감해 같은 Ro_cv,F 내에서도 산포가 존재한다. 이는 작은 스케일이 점성·열확산에 의해 제어되기 때문이다. 그러나 Ro_cv,F 가 충분히 크게(>1) 되면 Re 역시 ν, κ 에 독립적인 ‘극한’ 스케일에 접근한다.

또한, 고초임계(초임계성 Ra/Ra_c ≫ 1) 상황에서는 전단 흐름의 구조가 점점 ‘확산‑자유’ 형태로 수렴한다. 즉, 전단 흐름의 강도와 위상은 초기에는 확산에 의존하지만, 난류가 충분히 발달하면 내부 가열·냉각이 주도하는 대류 패턴이 지배하게 된다. 이는 기존 경계‑구동 모델에서 관찰된 ‘프로그레시브 → 레트로그레시브’ 전환이 확산에 의해 인위적으로 촉진되는 것이 아니라, 실제 천체 내부에서 기대되는 로스비수 기준에 의해 자연스럽게 발생한다는 점을 시사한다.

결론적으로, 내부 가열·냉각 구성을 통한 구형 껍질 대류 모델은 ‘확산‑자유’ 열전달과 로스비수 기반 전단 흐름 전환을 성공적으로 재현한다. 이는 천체 내부(태양, 거대 행성)에서의 대류 모델링에 있어, 비현실적으로 큰 점성·열확산을 사용해야 하는 전통적 접근법을 대체할 수 있는 실용적인 경로를 제공한다.