회전 레이리 벤드 대류의 전이와 경계조건에 따른 흐름 구분 기준

초록

본 연구는 회전 레이리-벤드 대류에서 무슬립·프리슬립 경계조건을 모두 고려한 직접 수치 시뮬레이션을 수행한다. 열 전달 효율(Nusselt 수)과 속도장 헬리시티를 각각 지표로 삼아, 코리올리 힘이 지배하는 흐름과 회전에 거의 영향을 받지 않는 흐름을 구분하는 기준을 탐색한다. 열 전달 기준에서는 두 경계조건 모두 동일한 조합의 레이놀즈·프루드·에크만 수값에서 전이가 일어나지만, 헬리시티 기준에서는 전이 위치가 서로 다르게 나타난다.

상세 요약

이 논문은 회전 레이리-벤드(RRB) 대류의 전이 현상을 정량적으로 규명하기 위해, 3차원 직접 수치 시뮬레이션(Direct Numerical Simulation, DNS)을 수행하였다. 시뮬레이션은 가로·세로 길이가 L인 정육면체 셀을 사용하고, 상하 경계는 온도 차이 ΔT를 유지하도록 설정하였다. 회전축은 중력과 일치하도록 수직으로 가정하고, 회전 속도 Ω에 따라 에크만 수 Ek = ν/(2ΩL²) 를 정의한다. 물성은 프루드 수 Pr = ν/κ 로 나타내며, 레이놀즈 수는 전형적인 대류 속도 U를 이용해 Re = UL/ν 로 계산한다.

두 종류의 경계조건—무슬립(no‑slip)과 자유슬립(free‑slip)—에 대해 각각 Ra = gαΔTL³/(νκ) 를 10⁶–10⁹ 범위, Pr = 0.7–7, Ek = 10⁻⁶–10⁻⁴ 범위에서 변화를 주었다. 주요 관측량은 열 전달 효율을 나타내는 Nusselt 수 Nu와, 속도장의 헬리시티 H = ⟨u·(∇×u)⟩/⟨|u|²⟩ 로 정의한다.

핵심 결과는 다음과 같다. 첫째, Nu가 회전에 무감각한(즉, 비회전 대류와 거의 동일한) 값을 보이는 전이점은 Re·Pr·Ek^{1/2} 혹은 Re·Pr^{1/2}·Ek^{1/2} 와 같은 특정 조합이 일정한 임계값을 초과할 때 일어난다. 이 임계값은 무슬립과 자유슬립 경계조건 모두에서 거의 동일하게 나타나며, 따라서 열 전달 기반 전이 기준은 경계조건에 독립적이다. 둘째, 헬리시티를 기준으로 전이를 정의하면, 전이 임계값이 경계조건에 따라 크게 달라진다. 무슬립 경우에는 헬리시티가 급격히 감소하는 임계점이 Ek·Re² ≈ const 로 나타나는 반면, 자유슬립에서는 훨씬 높은 Re·Ek 비율에서 헬리시티가 감소한다. 이는 경계층에서 발생하는 코리올리 효과와 경계조건에 따른 전단 응력 차이가 헬리시티 생성에 중요한 역할을 함을 시사한다.

또한, 전이 전후의 흐름 구조를 시각화한 결과, 전이 전에는 강한 코리올리 억제 하에 열 플럼프(plume)가 얇은 코리올리 층을 따라 회전하며, 전이 후에는 플럼프가 보다 자유롭게 상승·하강하면서 전통적인 비회전 대류와 유사한 대규모 롤 구조가 형성된다. 헬리시티 기반 전이에서는 특히 상하 경계층에서의 스위스 롤(swiss roll) 형태의 소용돌이가 사라지는 것이 관찰되었다.

이러한 결과는 지구 외핵과 같은 고속 회전 유체에서 열전달과 자기장 생성(다이너모) 메커니즘을 구분하는 데 중요한 함의를 가진다. 열 전달은 경계조건에 크게 좌우되지 않으나, 헬리시티—즉, 흐름의 비틀림—는 경계조건에 민감하므로, 실제 지구 외핵과 같은 경우에는 경계층 물성(예: 고체-액체 접촉면)의 특성을 정확히 반영해야 한다는 점을 강조한다.