Chapter 8

외부 유동과 경계층

Blasius 경계층 해, 층류-난류 전이, 그리고 평판 마찰 항력을 대화형으로 탐구합니다.

Blasius 해 경계층 두께 층류/난류 전이 마찰 계수

경계층의 개념

1904년 Prandtl은 높은 레이놀즈 수 유동에서 점성 효과가 물체 표면 근처의 얇은 층에 집중된다는 경계층(boundary layer) 개념을 제시했습니다. 경계층 바깥은 비점성 유동으로 취급할 수 있어, Euler 방정식과 Navier-Stokes 방정식을 각 영역에 맞게 적용하는 것이 가능해졌습니다.

경계층 두께 $\delta$는 유속이 자유류 속도 $U_\infty$의 99%에 도달하는 벽면으로부터의 거리로 정의됩니다. 배제 두께(displacement thickness) $\delta^*$와 운동량 두께(momentum thickness) $\theta$는 경계층의 질량 및 운동량 결손을 정량적으로 나타냅니다.

$$\delta^* = \int_0^\infty \left(1 - \frac{u}{U_\infty}\right) dy$$
$$\theta = \int_0^\infty \frac{u}{U_\infty}\left(1 - \frac{u}{U_\infty}\right) dy$$

Blasius 경계층 해

Blasius는 평판 위 층류 경계층에 대한 상사 해(similarity solution)를 도출했습니다. 상사 변수 $\eta = y\sqrt{U/(nu x)}$를 도입하면, 편미분 방정식이 상미분 방정식(Blasius 방정식) $f''' + \frac{1}{2}f f'' = 0$으로 환원됩니다. 무차원 속도 분포 $u/U = f'(\eta)$는 $\eta$만의 함수가 되어, 모든 $x$ 위치에서 동일한 형태를 보입니다.

이 해로부터 경계층 두께 $\delta \approx 5.0\sqrt{\nu x / U}$, 벽면 전단 응력 $\tau_w = 0.332 \,\rho\, U^2 / \sqrt{Re_x}$가 도출됩니다. $\delta$는 $\sqrt{x}$에 비례하므로 경계층은 하류로 갈수록 천천히 성장합니다.

$$\delta(x) = \frac{5.0\,x}{\sqrt{Re_x}}, \qquad Re_x = \frac{U\,x}{\nu}$$
$$\delta^*(x) = \frac{1.72\,x}{\sqrt{Re_x}}, \qquad C_f = \frac{0.664}{\sqrt{Re_x}}$$
Figure 8.1 Blasius 경계층 속도 분포
위치 x 0.20 m
자유류 U 3.0 m/s
δ(x)
δ*(x)
Rex
상태
Q. 평판 위 층류 경계층 두께 $\delta$가 $x$의 $\sqrt{x}$에 비례하는 물리적 이유는?
경계층의 성장은 벽면에서의 점성 확산(viscous diffusion)에 의해 결정됩니다. 확산 거리는 $\sqrt{\nu t}$에 비례하며, 유체 입자가 선단(leading edge)에서 $x$ 위치까지 이동하는 데 걸리는 시간은 $t \sim x/U$입니다. 따라서 $\delta \sim \sqrt{\nu x / U} = x / \sqrt{Re_x}$가 됩니다.

층류-난류 전이

평판 경계층은 선단에서 층류로 시작하여 하류로 가면서 임계 레이놀즈 수 $Re_{x,\text{crit}} \approx 5 \times 10^5$에 도달하면 난류로 전이됩니다. 전이 위치는 자유류 난류 강도, 표면 거칠기, 압력 구배 등에 의해 앞당겨지거나 지연될 수 있습니다.

난류 경계층은 층류보다 훨씬 두껍고 벽면 전단 응력이 높지만, 유동 분리에 대한 저항성이 강합니다. 난류 경계층 두께는 1/7 멱법칙에 의해 $\delta \approx 0.37 x / Re_x^{1/5}$로 근사됩니다.

$$\delta_{\text{turb}} \approx \frac{0.37\,x}{Re_x^{1/5}}, \qquad C_{f,\text{turb}} = \frac{0.0592}{Re_x^{1/5}}$$

평판 마찰 항력

평판의 총 마찰 항력은 국소 마찰 계수 $C_f(x)$를 평판 길이에 걸쳐 적분하여 구합니다. 평균 마찰 계수 $\bar{C}_f$는 층류 영역(Blasius)에서 $\bar{C}_f = 1.328 / \sqrt{Re_L}$, 난류 영역(1/7 멱법칙)에서 $\bar{C}_f = 0.074 / Re_L^{1/5}$입니다.

혼합 경계층(선단 부근 층류 + 하류 난류)의 경우, 전이 위치를 고려한 복합 상관식을 사용합니다. 전이 $Re_{x,\text{crit}} = 5 \times 10^5$를 가정하면 $\bar{C}_f = 0.074 Re_L^{-1/5} - 1742 Re_L^{-1}$이 대표적입니다.

$$\bar{C}_{f,\text{lam}} = \frac{1.328}{\sqrt{Re_L}}, \qquad \bar{C}_{f,\text{turb}} = \frac{0.074}{Re_L^{1/5}}$$
Figure 8.2 평판 마찰 계수 분포
평판 길이 L 1.0 m
U 5.0 m/s
ReL
평균 Cf (층류)
평균 Cf (난류)
전이 위치
Q. 동일 Re에서 난류 경계층의 벽면 전단 응력이 층류보다 큰 이유는?
난류 경계층에서는 격렬한 혼합(레이놀즈 응력)에 의해 높은 운동량의 외부 유체가 벽면 근처로 운반됩니다. 이로 인해 벽면 바로 위에 매우 얇은 점성 저층(viscous sublayer)이 형성되고, 여기서의 속도 구배 $du/dy$가 층류보다 훨씬 가파릅니다. $\tau_w = \mu (du/dy)_{y=0}$이므로 전단 응력이 커집니다.
Q. Blasius 해가 성립하기 위한 핵심 가정이 아닌 것은?
Blasius 해는 비압축성, 정상 상태, 제로 압력 구배(평판), 2차원 층류 경계층에 대한 해석해입니다. 난류 경계층에는 적용할 수 없으며, 난류 영역에서는 1/7 멱법칙 등의 경험적 모델을 사용합니다.