[논문 리뷰] Geostrophic convective turbulence: The effect of boundary layers

Rudie Kunnen, Rodolfo Ostilla–Mónico외 3명|arXiv (Cornell University)|2014. 09. 23.

Fluid Dynamics and Turbulent Flows참고 문헌 43인용 수 28

한 줄 요약

이 연구는 다양한 에크만 수(Ek = 1.3×10⁻⁷에서 2×10⁻⁶)와 레일라이흐 수(Ra = 10¹⁰ 및 5×10¹⁰)를 사용한 직접 수치 시뮬레이션을 통해 회전하는 레일라이흐–베나르 대류에서 지구풍성 대류로의 전이를 조사한다. 이는 슬립 경계 조건과 완전 슬립 경계 조건의 영향을 중심으로 한다. 두 조건 모두 유사한 에크만 수(~9×10⁻⁷ 및 ~3×10⁻⁷)에서 전이가 발생하지만, 유동 현상은 다르다: 완전 슬립 판은 반전 에너지 캐스케이드를 동반한 대규모 바로트로픽 소용돌이를 형성하는 반면, 슬립 판은 공명 구조를 억제하고 반전 캐스케이드를 제거한다. 이는 경계 조건이 지구풍성 난류 역학에 결정적인 영향을 미친다는 것을 시사한다.

ABSTRACT

Rayleigh--Bénard (RB) convection, the flow in a fluid layer heated from below and cooled from above, is used to analyze the transition to the geostrophic regime of thermal convection. In the geostrophic regime, which is of direct relevance to most geo- and astrophysical flows, the system is strongly rotated while maintaining a sufficiently large thermal driving to generate turbulence. We directly simulate the Navier--Stokes equations for two values of the thermal forcing, i.e. $Ra=10^{10}$ and $Ra=5\cdot10^{10}$, a constant Prandtl number~$Pr=1$, and vary the Ekman number in the range $Ek=1.3\cdot10^{-7}$ to $Ek=2\cdot10^{-6}$ which satisfies both requirements of super-criticality and strong rotation. We focus on the differences between the application of no-slip vs. stress-free boundary conditions on the horizontal plates. The transition is found at roughly the same parameter values for both boundary conditions, i.e. at~$Ek\approx 9 imes 10^{-7}$ for~$Ra=1 imes 10^{10}$ and at~$Ek\approx 3 imes 10^{-7}$ for~$Ra=5 imes 10^{10}$. However, the transition is gradual and it does not exactly coincide in~$Ek$ for different flow indicators. In particular, we report the characteristics of the transitions in the heat transfer scaling laws, the boundary-layer thicknesses, the bulk/boundary-layer distribution of dissipations and the mean temperature gradient in the bulk. The flow phenomenology in the geostrophic regime evolves differently for no-slip and stress-free plates. For stress-free conditions the formation of a large-scale barotropic vortex with associated inverse energy cascade is apparent. For no-slip plates, a turbulent state without large-scale coherent structures is found; the absence of large-scale structure formation is reflected in the energy transfer in the sense that the inverse cascade, present for stress-free boundary conditions, vanishes.

연구 동기 및 목표

• 비틀림 경계 조건과 완전 슬립 경계 조건 하에서의 회전하는 레일라이흐–베나르 대류에서 지구풍성 난류로의 전이를 조사한다.
• 슬립 및 완전 슬립 판에서 지구풍성 영역으로의 전이가 동일한 매개변수 값에서 발생하는지 확인한다.
• 경계 조건이 지구풍성 영역에서 열전달 척도 법칙, 경계층 두께, 소산 분포 및 유동 구조에 미치는 영향을 분석한다.
• 이론적 척도 법칙(예: 줄리엔 등(2012a) 및 킹 등)이 다양한 경계 조건 하에서 유효한지 평가한다.
• 특히 경계층과 부피 동역학의 역할을 중심으로 전이를 이끄는 물리적 메커니즘을 명확히 한다.

제안 방법

• 등점성 수 Pr = 1을 가진 비압축성 나비에-스토크스 방정식의 직접 수치 시뮬레이션을 수행하고, 두 레일라이흐 수(Ra = 10¹⁰ 및 5×10¹⁰)를 사용한다.
• 강한 회전과 초임계 대류를 달성하기 위해 에크만 수(Ek = 1.3×10⁻⁷에서 2×10⁻⁶)를 체계적으로 변화시킨다.
• 수평 판에서 슬립 조건과 완전 슬립 조건 하의 유동 행동을 비교한다.
• Ra와 Ek에 대한 Nusselt 수(Nu)를 통해 열전달 척도를 분석한다.
• 열 및 운동 에너지 경계층 두께를 측정하고, 에크만 수에 따른 척도 법칙을 비교한다.
• 난류 운동 에너지와 열분산 분산의 공간 분해를 부피와 경계층 간에 수행한다.

실험 결과

연구 질문

• RQ1 Ra = 10¹⁰ 및 5×10¹⁰ 조건에서 슬립 및 완전 슬립 경계 조건 하에서 지구풍성 대류로의 전이가 일어나는 에크만 수는 얼마인가?
• RQ2 지구풍성 영역에서 슬립 및 완전 슬립 판 간의 열전달 척도 법칙(Nu ∼ Ra^β Ek^γ)은 어떻게 다를까?
• RQ3 경계층, 특히 열 및 운동 에너지 경계층 두께 척도 법칙 측면에서 전이에 미치는 역할는 무엇인가?
• RQ4 전이 과정에서 부피와 경계층 간 운동 에너지 및 열분산 소산 분포는 어떻게 변화하는가?
• RQ5 왜 지구풍성 영역에서 완전 슬립 조건 하에서는 반전 에너지 캐스케이드가 지속되지만, 슬립 조건 하에서는 사라지는가?

주요 결과

• 지구풍성 대류로의 전이가 Ra = 1×10¹⁰일 경우 약 Ek ≈ 9×10⁻⁷, Ra = 5×10¹⁰일 경우 약 Ek ≈ 3×10⁻⁷에서 발생하며, 슬립 및 완전 슬립 경계 조건 모두에서 유사한 전이 임계값을 보인다.
• 완전 슬립 판의 경우 대규모 바로트로픽 소용돌이가 형성되며 반전 에너지 캐스케이드가 발생하지만, 슬립 판의 경우 대규모 공명 구조가 없는 난류 상태를 나타낸다.
• 완전 슬립 판의 Nusselt 수 척도 법칙은 줄리엔 등(2012a)의 이론적 척도 법칙(Nu ∼ Ra^{3/2} Ek²)과 가장 잘 일치하지만, 슬립 판은 기존 이론으로는 설명되지 않는 다른 척도 법칙 행동을 보인다.
• 열 경계층 두께 척도 법칙은 Ra = 1×10¹⁰일 경우 약 Ek ≈ 8×10⁻⁷, Ra = 5×10¹⁰일 경우 약 Ek ≈ 3×10⁻⁷에서 급격히 변화하며, 이는 경계층 내부의 구조적 전이를 시사한다.
• 열분산 소산(εᵤ)은 주로 부피에서 발생하며, 전이 지점에서 척도 법칙이 변화함에 따라 주요 소산 메커니즘이 변화함을 시사한다.
• 전이는 점진적이며 모든 유동 지표가 완전히 동기화되지는 않으며, 지구풍성 영역의 시작을 명확히 정의할 수 있는 단일 기준은 존재하지 않음을 시사한다.