TML

냉동공학

레일리수 Ra의 의미: 자연대류의 발생·발달과 Nu의 관계

레일리수(Rayleigh Number)는 자연대류를 해석할 때 매우 중요한 무차원수이다. 앞에서 배운 그라쇼프수 Gr는 온도차에 의해 발생하는 부력과 점성의 상대적인 크기를 나타냈고, 프란틀수 Pr는 운동량확산과 열확산의 상대적인 크기를 나타냈다.

레일리수는 이 두 무차원수를 결합한다.

\[
\boxed{
Ra=GrPr
}
\]

따라서 자연대류에서 부력에 의해 유체가 움직이려는 효과뿐 아니라 점성과 열확산의 영향까지 함께 고려할 수 있다.

목차25개 항목

레일리수는 왜 필요한가

자연대류의 출발점은 표면과 주변유체 사이의 온도차이다.

\[
\boxed{
\text{온도차}
\rightarrow
\text{밀도차}
\rightarrow
\text{부력}
\rightarrow
\text{유체운동}
}
\]

앞에서 그라쇼프수를

\[
Gr
=
\frac{
g\beta\Delta TL^3
}{
\nu^2
}
\]

로 정의했다.

그리고 물리적으로

\[
Gr
\sim
\frac{
\text{부력효과}
}{
\text{점성효과}
}
\]

라고 이해했다.

그러나 자연대류에는 한 가지 요소가 더 중요하다. 바로 열이 유체 내부에서 얼마나 빠르게 확산되는가이다.

자연대류의 원인이 온도차이므로 열이 주변으로 너무 빠르게 확산되면 온도차와 밀도차가 빠르게 완화될 수 있다.

따라서 자연대류를 종합적으로 이해하려면

\[
\boxed{
\text{부력}
+
\text{점성}
+
\text{열확산}
}
\]

을 함께 살펴볼 필요가 있다.

이를 나타내는 대표적인 무차원수가 레일리수이다.

레일리수의 정의

레일리수는 그라쇼프수와 프란틀수의 곱으로 정의한다.

\[
\boxed{
Ra
=
GrPr
}
\]

그라쇼프수는

\[
Gr
=
\frac{
g\beta\Delta TL^3
}{
\nu^2
}
\]

이고 프란틀수는

\[
Pr
=
\frac{\nu}{\alpha}
\]

이므로 두 식을 곱하면

\[
Ra
=
\frac{
g\beta\Delta TL^3
}{
\nu^2
}
\frac{\nu}{\alpha}
\]

가 된다.

동점성계수 ν 하나가 약분되므로

\[
\boxed{
Ra
=
\frac{
g\beta\Delta TL^3
}{
\nu\alpha
}
}
\]

이다.

따라서 다음 두 식을 함께 기억하면 좋다.

\[
\boxed{
Ra
=
GrPr
=
\frac{
g\beta\Delta TL^3
}{
\nu\alpha
}
}
\]

기호의미SI 단위
g중력가속도m/s²
β체적팽창계수1/K
ΔT표면과 주변유체의 온도차K
L특성길이m
ν동점성계수m²/s
α열확산계수m²/s

레일리수의 물리적 의미

레이놀즈수 Re는 관성효과와 점성효과를 비교했다.

\[
Re
\sim
\frac{
\text{관성효과}
}{
\text{점성효과}
}
\]

그라쇼프수 Gr는 부력효과와 점성효과를 비교했다.

\[
Gr
\sim
\frac{
\text{부력효과}
}{
\text{점성효과}
}
\]

레일리수는 여기에서 열확산까지 고려한다.

개념적으로는

\[
\boxed{
Ra
\sim
\frac{
\text{부력에 의한 자연대류 발생효과}
}{
\text{점성과 열확산에 의한 억제효과}
}
}
\]

라고 이해할 수 있다.

즉 Ra가 커진다는 것은 부력이 점성과 열확산의 영향을 이겨내고 자연대류가 발달하기 좋은 조건이 된다는 의미이다.

그라쇼프수와 프란틀수가 레일리수로 결합되고 누셀수와 대류열전달계수로 연결되는 자연대류 계산 흐름

왜 열확산이 중요한가

뜨거운 벽 가까이의 유체가 가열된다고 생각해보자.

유체온도가 상승하면 일반적으로 밀도가 감소하고 부력이 발생한다.

\[
T\uparrow
\rightarrow
\rho\downarrow
\rightarrow
\text{부력 발생}
\]

그러나 열확산계수 α가 매우 크다면 벽 근처의 열이 주변유체로 빠르게 퍼진다.

그러면 국부적인 온도차가 완화되고 자연대류를 일으키는 밀도차도 작아질 수 있다.

레일리수에서

\[
Ra
=
\frac{
g\beta\Delta TL^3
}{
\nu\alpha
}
\]

이므로 다른 조건이 같다면

\[
\boxed{
\alpha\uparrow
\quad\Rightarrow\quad
Ra\downarrow
}
\]

한다.

점성은 자연대류를 어떻게 억제하는가

점성은 유체의 움직임을 억제하는 방향으로 작용한다.

레일리수의 분모에는 동점성계수 ν가 포함되어 있다.

\[
\nu\uparrow
\quad\Rightarrow\quad
Ra\downarrow
\]

즉 점성이 큰 유체에서는 온도차로 부력이 발생하더라도 유체가 쉽게 움직이지 않는다.

따라서 레일리수의 분모

\[
\boxed{
\nu\alpha
}
\]

는 각각 유체운동을 억제하는 점성의 영향과 온도차를 완화하는 열확산의 영향을 나타낸다고 이해할 수 있다.

Ra 식을 변수별로 해석해보자

레일리수는

\[
Ra
=
\frac{
g\beta\Delta TL^3
}{
\nu\alpha
}
\]

이다.

온도차가 커지면

\[
\Delta T\uparrow
\quad\Rightarrow\quad
Ra\uparrow
\]

한다.

온도차가 커지면 밀도차와 부력이 커지기 때문이다.

체적팽창계수가 커져도

\[
\beta\uparrow
\quad\Rightarrow\quad
Ra\uparrow
\]

한다.

이는 같은 온도차에서도 밀도가 더 민감하게 변화한다는 뜻이다.

특성길이가 증가하면

\[
L\uparrow
\quad\Rightarrow\quad
Ra\uparrow
\]

하며 특히

\[
\boxed{
Ra\propto L^3
}
\]

이다.

따라서 자연대류에서는 장치의 크기와 형상이 매우 중요한 영향을 미친다.

Ra가 작다는 것은 무엇인가

Ra가 작으면 부력에 비해 점성과 열확산의 영향이 상대적으로 크다.

따라서 유체의 움직임이 약하고 전도적인 열전달의 비중이 상대적으로 커진다.

\[
\boxed{
Ra\downarrow
\rightarrow
\text{자연대류 약화}
}
\]

특정 밀폐 유체층과 같은 조건에서는 Ra가 충분히 작으면 유체운동이 발생하지 않고 거의 전도에 의해 열이 전달되는 경우도 있다.

Ra가 크다는 것은 무엇인가

Ra가 커지면 부력이 점성과 열확산의 영향을 극복하기 쉬워진다.

따라서 자연대류 유동이 보다 강하게 발달하는 경향을 보인다.

\[
\boxed{
Ra\uparrow
\rightarrow
\text{자연대류 발달 증가}
}
\]

많은 자연대류 상관식에서는 Ra가 증가하면 Nu도 증가한다.

\[
Ra\uparrow
\rightarrow
Nu\uparrow
\rightarrow
h\uparrow
\]

따라서 같은 면적과 온도차에서는 열전달량도 증가하는 방향으로 작용할 수 있다.

자연대류의 전체 관계가 완성된다

앞에서 자연대류의 누셀수를

\[
Nu
=
f(Gr,Pr)
\]

형태로 표현했다.

그런데

\[
Ra
=
GrPr
\]

이므로 많은 자연대류 상관식은

\[
\boxed{
Nu
=
f(Ra)
}
\]

형태로 표현할 수 있다.

경우에 따라서는 Pr의 영향을 별도로 고려하여

\[
\boxed{
Nu
=
f(Ra,Pr)
}
\]

형태를 사용하기도 한다.

단순화된 경험 상관식에서는 흔히

\[
\boxed{
Nu
=
CRa^n
}
\]

형태를 볼 수 있다.

계수 C와 지수 n은 형상, 유동영역과 적용범위에 따라 달라진다.

Gr와 Ra의 차이를 확실히 하자

구분GrRa
정의\(g\beta\Delta TL^3/\nu^2\)\(GrPr\)
핵심 영향부력 ↔ 점성부력·점성·열확산
Pr 포함아니오예
주요 역할자연대류 유동의 구동특성자연대류 발달과 열전달 상관식
Nu와의 관계\(Nu=f(Gr,Pr)\)\(Nu=f(Ra)\) 또는 \(f(Ra,Pr)\)

따라서

\[
\boxed{
Gr
=
\text{부력과 점성의 경쟁}
}
\]

이라고 본다면

\[
\boxed{
Ra
=
\text{열확산까지 포함한 자연대류 발달 지표}
}
\]

라고 이해할 수 있다.

왜 Ra = GrPr인가

각 무차원수의 의미를 다시 살펴보자.

\[
Gr
\sim
\frac{
\text{부력}
}{
\text{점성}
}
\]

이고

\[
Pr
=
\frac{
\text{운동량확산}
}{
\text{열확산}
}
\]

이다.

이를 곱하면 자연대류를 만드는 부력의 영향과 점성·열확산의 특성을 함께 포함할 수 있다.

수식으로는

\[
\boxed{
Ra
=
\frac{
g\beta\Delta TL^3
}{
\nu\alpha
}
}
\]

이며, 분자는 자연대류를 발생시키는 방향의 효과, 분모는 유동과 온도차를 완화하는 효과로 이해할 수 있다.

자연대류의 층류와 난류

강제대류에서는 Re가 층류와 난류 판단의 중요한 기준이었다.

자연대류에서는 Ra가 중요한 판단지표가 된다.

예를 들어 수직평판 외부 자연대류에서는 대표적인 공학적 기준으로

\[
Ra_L
\lesssim
10^9
\]

범위를 층류 자연대류로 다루는 경우가 많고,

\[
Ra_L
\gtrsim
10^9
\]

부근에서는 천이와 난류의 영향을 중요하게 고려한다.

그러나

\[
\boxed{
Ra=10^9
}
\]

을 모든 자연대류 문제의 절대적인 임계값으로 사용해서는 안 된다.

형상, 중력방향과 표면방향의 관계, 경계조건, 표면상태 등에 따라 자연대류의 발달특성이 달라지기 때문이다.

왜 형상에 따라 Ra 기준이 달라지는가

자연대류는 중력방향과 표면방향의 상대적인 관계에 크게 영향을 받는다.

예를 들어 다음 형상은 각각 자연대류 유동구조가 다르다.

  • 수직평판
  • 위쪽으로 가열되는 수평평판
  • 아래쪽으로 가열되는 수평평판
  • 수평원통
  • 구
  • 밀폐공간

따라서

\[
Nu
=
f(Ra)
\]

라는 기본적인 구조는 같더라도 사용하는 상관식의 C와 n, 특성길이와 적용 가능한 Ra 범위가 달라진다.

\[
\boxed{
\text{자연대류 상관식은 형상과 적용범위를 함께 확인한다}
}
\]

Ra의 의미를 잘 보여주는 Rayleigh-Bénard 대류

레일리수의 물리적 의미를 이해하기 좋은 대표적인 예가 아래쪽에서 가열되는 수평 유체층이다.

       차가운 면
────────────────
   ↓      ↓
 ↙   ↑  ↑   ↘
    대류셀
 ↘   ↓  ↓   ↙
────────────────
       뜨거운 면

온도차가 매우 작으면 부력이 충분하지 않아 유체가 거의 움직이지 않고 전도에 의해 열이 전달될 수 있다.

그러나 Ra가 일정한 임계값 이상으로 커지면 부력이 점성과 열확산의 억제효과를 극복하면서 대류셀이 발생한다.

이상화된 대표적인 Rayleigh-Bénard 조건에서는 대류가 발생하기 시작하는 임계 레일리수가 약

\[
\boxed{
Ra_{cr}
\approx
1708
}
\]

로 알려져 있다.

하지만 이 값은 특정한 수평 유체층의 안정성 문제에 해당한다.

따라서 수직평판에서 사용하는 \(Ra_L\sim10^9\) 수준의 천이 기준과 비교하여 같은 기준이라고 생각해서는 안 된다.

작은 Ra와 큰 Ra의 자연대류 비교

자연대류 Nu 상관식의 대표적인 형태

수직평판의 평균 자연대류를 단순한 형태로 설명할 때 특정 적용범위에서는

층류영역에 대해

\[
\overline{Nu}_L
\approx
0.59Ra_L^{1/4}
\]

형태의 경험식을 볼 수 있다.

또한 특정 난류영역에서는

\[
\overline{Nu}_L
\approx
0.10Ra_L^{1/3}
\]

형태가 사용되기도 한다.

이 식들의 계수와 적용범위는 특정 조건에 대한 경험적 상관식이다.

따라서 기초단계에서는 숫자 자체를 무조건 암기하기보다

\[
\boxed{
Nu
=
CRa^n
}
\]

이라는 구조와 상관식마다 적용조건이 존재한다는 점을 이해하는 것이 우선이다.

왜 Ra가 증가하면 Nu가 증가하는가

누셀수는

\[
Nu
=
\frac{hL}{k}
\]

로 정의한다.

Ra가 증가하면 부력에 의한 유체운동이 강해지는 경향이 있다.

유체운동이 활발해지면 벽면 근처에서 열이 효과적으로 운반되고 열경계층의 전열도 강화된다.

따라서 일반적인 자연대류 상관식에서는

\[
\boxed{
Ra\uparrow
\rightarrow
Nu\uparrow
\rightarrow
h\uparrow
\rightarrow
\dot Q\uparrow
}
\]

하는 경향을 볼 수 있다.

다만 실제 증가율은 형상과 해당 상관식의 지수에 따라 달라진다.

앞에서 계산한 수직평판 예제를 이어가보자

앞의 그라쇼프수 단원에서 높이 0.5 m인 수직평판에 대해 다음 값을 구했다.

\[
Gr_L
\approx
3.2\times10^8
\]

공기의 프란틀수를

\[
Pr
\approx
0.71
\]

이라고 하자.

사용 식

\[
Ra_L
=
Gr_LPr
\]

대입

\[
Ra_L
=
3.2\times10^8
\times
0.71
\]

계산 결과

\[
\boxed{
Ra_L
\approx
2.27\times10^8
}
\]

대표적인 수직평판 기준에서

\[
Ra_L<10^9 \]

이므로 이 예제는 층류 자연대류 영역으로 다룰 수 있다.

이제 Nu를 계산해보자

예시를 위해 해당 범위에 적용 가능한 단순한 층류 수직평판 상관식을 사용한다고 하자.

\[
\overline{Nu}_L
=
0.59Ra_L^{1/4}
\]

근사적으로

\[
Ra_L
\approx
2.3\times10^8
\]

을 사용하면

\[
\overline{Nu}_L
=
0.59
\left(
2.3\times10^8
\right)^{1/4}
\]

따라서

\[
\boxed{
\overline{Nu}_L
\approx
72.7
}
\]

정도가 된다.

대류열전달계수 h까지 구해보자

누셀수 정의에서

\[
\overline{Nu}_L
=
\frac{
\bar hL
}{
k
}
\]

이므로

\[
\boxed{
\bar h
=
\frac{
\overline{Nu}_Lk
}{
L
}
}
\]

이다.

평판 높이가

\[
L
=
0.5\mathrm{m}
\]

이고 공기의 열전도율을 예시로

\[
k
=
0.026\mathrm{W/(m\cdot K)}
\]

라고 하면

\[
\bar h
=
\frac{
72.7\times0.026
}{
0.5
}
\]

따라서

\[
\boxed{
\bar h
\approx
3.78\mathrm{W/(m^2\cdot K)}
}
\]

이다.

즉 처음에는 자연대류의 h를 알 수 없었지만

\[
\boxed{
Gr
\rightarrow
Ra
\rightarrow
Nu
\rightarrow
h
}
\]

의 순서로 실제 대류열전달계수를 계산할 수 있었다.

자연대류 문제의 전체 계산 순서

표면온도 Ts
주변온도 T∞
      ↓
ΔT = |Ts - T∞|
      ↓
막온도 Tf
      ↓
물성값 결정
β, ν, α, Pr, k
      ↓
Gr = gβΔTL³/ν²
      ↓
Ra = Gr·Pr
      ↓
형상과 Ra 범위 확인
      ↓
Nu 상관식 선택
      ↓
Nu
      ↓
h = Nu·k/L
      ↓
Q̇ = hAΔT

자연대류 문제에서는 이 순서를 먼저 잡아두는 것이 중요하다.

강제대류와 자연대류를 최종 비교하면

강제대류에서는 외부에서 유속이 주어지고 Re가 출발점이 된다.

\[
\boxed{
Re+Pr
\rightarrow
Nu
\rightarrow
h
}
\]

자연대류에서는 온도차와 부력으로부터 Gr를 구하고 Pr와 결합하여 Ra를 만든다.

\[
\boxed{
Gr\times Pr
=
Ra
\rightarrow
Nu
\rightarrow
h
}
\]

구분강제대류자연대류
유동 원인팬·펌프 등온도차에 의한 부력
대표 유동 무차원수ReGr
열확산 특성PrPr
종합 자연대류 지표–Ra = GrPr
열전달 결과NuNu

Gr·Pr·Ra·Nu를 하나로 연결해서 기억하자

네 개의 무차원수를 각각 독립적으로 외우면 관계가 헷갈리기 쉽다.

그라쇼프수는

\[
\boxed{
Gr
=
\text{부력과 점성의 관계}
}
\]

프란틀수는

\[
\boxed{
Pr
=
\text{운동량확산과 열확산의 관계}
}
\]

이다.

이 둘을 결합한 레일리수는

\[
\boxed{
Ra
=
GrPr
=
\text{자연대류의 발달을 나타내는 종합 지표}
}
\]

이다.

그리고 자연대류에 의한 전열성능을 누셀수로 나타낸다.

\[
\boxed{
Nu
=
\frac{hL}{k}
}
\]

따라서 자연대류의 핵심 연결은

\[
\boxed{
Gr+Pr
\rightarrow
Ra
\rightarrow
Nu
\rightarrow
h
}
\]

이다.

기술사 기출에서는 어떻게 연결되는가

현재 확인한 기출 범위에서는 레일리수 Ra만을 직접 명시하여 설명하도록 한 문제는 확인되지 않았다.

그러나 제131회에서는 자연대류에 이용되는 중요한 무차원수를 설명하도록 출제됐다.

또한 제139회에서는 Grashof 수, Prandtl 수, Nusselt 수를 함께 설명하도록 요구했다.

제133회에서도 자연대류열전달의 누셀수를 무차원수 관계식으로 표현하고 각 무차원수를 정의하도록 출제됐다.

따라서 자연대류에 관한 기술사 답안을 구성할 때는

\[
Ra
=
GrPr
\]

와

\[
Nu
=
f(Ra)
\]

의 연결을 이해해 두면 Gr·Pr·Nu를 하나의 자연대류 체계로 설명하는 데 도움이 된다.

기술사 답안에서 중요한 포인트

레일리수를 설명할 때는 먼저 자연대류에서 부력·점성·열확산의 영향을 종합적으로 나타내는 무차원수라는 점을 제시하는 것이 좋다.

그다음

\[
Ra
=
GrPr
=
\frac{
g\beta\Delta TL^3
}{
\nu\alpha
}
\]

를 제시하고 각 변수의 의미를 설명한다.

이어 Ra가 작으면 점성과 열확산의 영향이 상대적으로 크고, Ra가 증가하면 자연대류가 발달하는 경향이 있다는 물리적 의미를 설명한다.

마지막으로

\[
Nu=f(Ra)
\]

와

\[
h
=
\frac{Nu\,k}{L}
\]

를 연결하면 실제 자연대류 열전달까지 이어지는 답안이 된다.

층류·난류의 임계값을 제시할 때는 반드시 적용되는 형상과 조건을 함께 명시하는 것이 좋다.

시험 직전 핵심 정리

  • 레일리수는 그라쇼프수와 프란틀수의 곱이다.
  • \(Ra=GrPr=g\beta\Delta TL^3/(\nu\alpha)\)이다.
  • Gr가 부력과 점성을 비교한다면 Ra는 열확산까지 함께 고려한다.
  • Ra가 작으면 점성과 열확산의 영향이 상대적으로 크고 자연대류가 약한 경향이 있다.
  • Ra가 증가하면 자연대류가 발달하고 일반적으로 Nu와 h도 증가하는 방향으로 작용한다.
  • 자연대류 상관식은 \(Nu=f(Ra)\), \(Nu=f(Ra,Pr)\), 또는 \(Nu=CRa^n\) 형태로 나타낼 수 있다.
  • 수직평판에서 \(Ra_L\sim10^9\) 전후가 천이의 대표적 기준으로 사용되지만 모든 자연대류 형상에 적용되는 절대값은 아니다.
  • Rayleigh-Bénard 문제의 \(Ra_{cr}\approx1708\)은 특정 수평 유체층 조건의 대류 발생 임계값이다.
  • 자연대류 계산은 \(Gr+Pr\rightarrow Ra\rightarrow Nu\rightarrow h\rightarrow\dot Q\) 순으로 연결한다.

다음 단계에서는 ⑨-7 비오수 \(Bi\)를 공부한다. 지금까지의 Re·Pr·Nu·Gr·Ra가 주로 유체의 대류현상을 해석하는 무차원수였다면, Bi부터는 고체 내부의 온도분포와 비정상 열전도로 관심이 이동한다. 특히 \(Bi=hL_c/k_s\)가 왜 고체 내부 전도열저항과 표면 대류열저항의 비가 되는지, 그리고 \(Bi<0.1\) 조건이 집중정수법과 어떻게 연결되는지를 살펴본다.