자연대류

시뮬레이션 전산유체역학 유체역학 마지막 수정: 2026-07-09 04:52:19

자연대류 Natural Convection
분야	유체역학 × 열전달 × 수치해석
구동력	부력 (밀도 차)
핵심 무차원수	레일리 수, 그라스호프 수, 프란틀 수
대표 근사	부시네스크 근사

1. 개요[편집]

자연대류(Natural Convection, 자유대류)는 팬이나 펌프 같은 외부 강제력 없이, 온도 차로 생긴 유체의 밀도 차와 중력이 만들어내는 부력(buoyancy)만으로 유동이 발생하는 열전달 현상이다. 따뜻해진 유체는 가벼워져 위로 뜨고 차가운 유체는 무거워져 가라앉으며, 이 순환이 스스로 흐름을 만들어 열을 실어 나른다.¹ 라디에이터 위로 아지랑이가 피어오르고, 뜨거운 커피잔에서 김이 하늘거리고, 대륙 규모의 대기 순환과 지구 맨틀의 대류까지 — 데워진 것이 뜨는 이 단순한 원리가 온갖 규모에서 반복된다.

강제대류(forced convection)와 달리 유동장과 온도장이 서로를 강하게 되먹임한다는 점이 자연대류의 본질이자 난점이다. 온도가 유동을 만들고(부력), 그 유동이 다시 온도를 재분포시킨다. 이 강한 연성 때문에 자연대류는 열전달 해석에서 가장 다루기 까다로운 축에 든다. 그럼에도 팬 없이 조용하고 전력 소모 없이 열을 뺄 수 있다는 미덕 때문에, 전자기기 냉각 설계자들에게는 놓칠 수 없는 카드다.

2. 지배하는 무차원수[편집]

자연대류의 거동을 결정하는 무차원수는 셋이다. 가장 중요한 것이 레일리 수(Rayleigh number) $Ra$ 로, 부력이 점성과 열확산의 저항을 얼마나 압도하는지를 나타낸다.

Ra = \frac{g \beta \Delta T\, L^3}{\nu \alpha}

여기서 $\beta$ 는 열팽창계수, $\Delta T$ 는 온도차, $L$ 은 특성 길이, $\nu$ 는 동점성계수, $\alpha$ 는 열확산계수다. $Ra$ 가 어떤 임계값(예: 수평 유체층에서 약 1708)을 넘으면 정지해 있던 유체가 갑자기 규칙적인 대류 셀을 이루며 움직이기 시작하는데, 이 유명한 불안정성이 레일리-베나르 대류(Rayleigh-Bénard convection)다. $Ra$ 가 더 커지면 유동은 층류에서 난류로 넘어간다.

레일리 수는 사실 두 무차원수의 곱으로 분해된다. 부력과 점성의 비를 나타내는 그라스호프 수(Grashof number) $Gr$ 과, 운동량 확산과 열확산의 비인 프란틀 수(Prandtl number) $Pr$ 이다.

Ra = Gr \cdot Pr, \qquad Gr = \frac{g\beta \Delta T\, L^3}{\nu^2}, \qquad Pr = \frac{\nu}{\alpha}

$Gr$ 은 자연대류에서 강제대류의 레이놀즈수에 해당하는 위치를 차지한다. 강제대류와 자연대류가 비등하게 섞인 혼합대류인지 판별할 때는 $Gr/Re^2$ 비를 본다 — 이 값이 1 근처면 어느 쪽도 무시할 수 없다는 뜻이다.

3. 부시네스크 근사[편집]

자연대류를 엄밀히 풀려면 밀도가 온도에 따라 변하는 압축성 문제를 다뤄야 하는데, 이건 압축성 유동만큼이나 골치 아프다. 다행히 대부분의 자연대류에서 온도차가 그리 크지 않다는 점을 이용한 우아한 지름길이 있다. 부시네스크 근사(Boussinesq approximation)다.

핵심 아이디어는 이렇다. 밀도 변화가 워낙 작으니 관성항·연속방정식 등 모든 곳에서는 밀도를 상수로 취급하되, 오직 중력항에서만 밀도 변화를 살려 부력으로 남긴다.² 그러면 부력항은 다음처럼 온도차에 선형으로 비례하는 깔끔한 형태가 된다.

\rho \mathbf{g} \approx \rho_0\big(1 - \beta (T - T_0)\big)\mathbf{g}

이 근사 덕분에 유동을 비압축성 나비에-스토크스 방정식의 틀 안에서 풀 수 있게 되고, 압축성의 지옥을 피해 갈 수 있다. 다만 온도차가 크면(대략 $\beta \Delta T \gtrsim 0.1$ ) 이 근사는 무너지므로, 그런 경우엔 밀도를 온도의 함수로 직접 다루는 완전 가변밀도 해석으로 넘어가야 한다.

4. 경계층과 열전달[편집]

자연대류에서도 경계층은 주인공이다. 가열된 수직 벽면을 생각해보자. 벽에 붙은 유체는 데워져 위로 상승하며 얇은 부력 경계층을 형성한다. 이 경계층 안에서 속도는 벽에서 0이었다가 부력으로 증가했다가 다시 주위 정지 유체로 접근하며 0이 되는 특유의 종 모양 분포를 그린다 — 강제대류 경계층에는 없는 자연대류만의 특징이다.

열전달의 세기는 누셀 수(Nusselt number) $Nu$ 로 표현되며, 자연대류에서는 $Nu = C\, Ra^n$ 꼴의 상관식이 널리 쓰인다. 지수 $n$ 은 층류에서 대략 $1/4$ , 난류에서 $1/3$ 부근이다.³ $n=1/3$ 이라는 값은 의미심장한데, 이 경우 $Nu \propto Ra^{1/3}$ 이 되어 특성 길이 $L$ 이 소거되므로 열전달 계수가 판의 크기와 무관해진다. 자연이 던지는 이런 우연한 선물은 설계 공학자에게 늘 반가운 법이다.

5. CFD 모델링과 응용[편집]

자연대류를 전산유체역학으로 풀 때의 실무적 어려움은 수렴성이다. 부력-속도-온도의 삼각 연성이 강해서 발산하기 쉽고, 특히 $Ra$ 가 높은 난류 영역에서는 난류 모델링의 부력 생성항을 제대로 켜야 한다. 부시네스크 근사를 켠 비압축성 솔버에 에너지 방정식을 연성하고, SIMPLE 알고리즘 계열로 압력-속도를 풀되 완화 계수를 낮춰 얌전히 달래는 것이 국룰이다. 정상상태 해가 잘 안 잡히면 아예 비정상(transient)으로 시간을 전진시키며 정상상태로 유도하는 꼼수도 흔하다.

응용 분야는 넓다. 전자기기 냉각에서는 팬 없는 무소음 방열판(heat sink) 설계의 핵심 도구이며, LED 조명, 전력 반도체, 노트북 저부하 냉각 등에서 자연대류만으로 얼마나 열을 뺄 수 있는지 예측한다. 건물 분야에서는 실내 온열 쾌적성 평가, 이중 외피(double skin façade)의 굴뚝 효과, 자연 환기 설계에 쓰인다. 그 밖에 원자로 사고 시 자연순환에 의한 잔열 제거, 태양열 집열기, 데이터센터 공조까지 — 전기 한 방울 안 쓰고 열을 옮기는 이 현상은 에너지 효율이 곧 미덕인 시대에 갈수록 몸값이 오르고 있다.

6. 관련 문서[편집]

“뜨거운 것이 위로 뜬다”는 게 워낙 당연해 보이지만, 무중력에서는 이 순환이 안 일어난다. 국제우주정거장에서 촛불이 강제 환기 없이는 산소 부족으로 꺼지는 이유가 바로 자연대류가 없어서다. 부력은 중력의 자식이다. ↩
밀도 변화를 어떤 곳에선 무시하고 어떤 곳에선 살린다니 얼핏 반칙 같지만, 척도 분석(scale analysis)으로 정당화된다. 관성항에서 밀도 변화의 기여는 무시할 만큼 작지만, 중력항에서는 그 작은 변화가 곧 유동을 낳는 유일한 구동력이라 절대 버리면 안 된다. 작은 것도 자리를 잘 만나면 주인공이 된다. ↩
상관식 계수 $C$ 와 지수 $n$ 은 형상(수직판, 수평판, 실린더, 구)과 $Ra$ 범위마다 다르며, 처칠-추(Churchill-Chu) 같은 실험 상관식이 핸드북에 빼곡히 정리돼 있다. CFD가 없던 시절 엔지니어들이 이 표 하나로 방열 설계를 다 해냈다는 걸 생각하면 절로 겸손해진다. ↩