물줄기가 흐트러지는 순간 — 레이놀즈 임계수와 난류 경계층의 다층 구조

첫 격자점을 로그층에 두면 점성저층이 계산에서 빠지고, 벽 마찰은 B 한 숫자에 기댄다.

수도꼭지를 천천히 열면 물줄기는 매끈한 유리 막대처럼 떨어진다. 그런데 손잡이를 조금 더 돌리는 순간, 매끈한 줄기가 잘게 부서지며 어수선해진다. 같은 물, 같은 노즐, 바뀐 건 속도뿐이다. 무엇이 그 임계점을 결정하는가? 1883년 맨체스터의 레이놀즈가 색소 한 줄기로 답한 이 질문은 결국 한 무차원수에 압축된다. 임계가 왜 Re ≈ 2300 근처에 놓이는지, 그 너머의 경계층이 어떻게 점성저층·완화층·로그 영역으로 갈리는지가 이 노트의 범위다. Re 슬라이더를 1500에 두면 색실은 곧고, 8000 쪽으로 밀면 관 전체가 색으로 퍼진다.

Reynolds number Re = U·d/ν: 1500Laminar (Re_crit ≈ 2300)

2300을 지나는 동안 끝만 먼저 흔들린다. 4000을 넘기 전에는 와가 구간구간 나타났다 사라진다. Pause로 그 순간을 붙잡으면, 천이가 한 프레임에 끝나지 않는 게 보인다.

매끈한 색실이 갑자기 흐트러지는 임계점#

레이놀즈는 유리관 한가운데로 색소를 주입했다. 관 직경 $d$ , 평균 유속 $u$ , 동점성계수 $\nu$ 를 다양하게 바꿔가며 한 가지 질문에 답하려 했다. 색실은 언제 흐트러지는가?

답은 한 무차원수의 값이었다. 점도, 유속, 관 직경을 어떻게 조합해 바꾸어도 흐트러짐은 거의 같은 값에서 시작했다.

\mathrm{Re} = \frac{u\, d}{\nu}

여기서 $u$ 는 평균 유속, $d$ 는 특성 길이(관에서는 직경), $\nu = \mu / \rho$ 는 동점성계수(점성계수/밀도). Re가 임계값을 넘으면 색실은 평형을 잃는다.

Re = 2300의 의미 — 관성과 점성의 비#

Re는 단순한 실험 상수가 아니다. 운동량 방정식 양쪽을 무차원화하면 점성항 앞에 정확히 $1/\mathrm{Re}$ 가 붙는다. 즉 Re는 관성력 대 점성력의 비다.

Re가 작다 → 점성이 우세, 작은 교란을 점성이 빠르게 균질화 → 층류가 살아남는다.
Re가 크다 → 관성이 우세, 작은 교란이 비선형 항을 통해 자기증폭 → 난류로 천이.

원관 내부에서 실험적으로 알려진 임계 Re는 약 2300이다. 단, 이 값은 외란 상태에 따라 흔들린다. 매끈하게 가공된 관, 진동 없는 받침대, 깨끗한 입구 조건에서는 Re = 5만 까지도 층류를 유지시킨 보고가 있다. 따라서 "임계 Re ≈ 2300"은 공학 실용 기준이고, 천이의 본질은 외란을 키울 수 있는 에너지 통로가 열렸느냐다. 같은 무차원수가 곤충의 헤엄과 배의 항적을 가르는 이야기는 Reynolds 수가 유동의 스케일을 바꾸는 방식에 따로 있다.

평판 위 경계층(boundary layer, 벽 근처에 점성이 지배하는 얇은 층)이라면 척도가 다르다. 전연(leading edge)으로부터의 거리 $x$ 로 정의한 $\mathrm{Re}_x = U_\infty x / \nu$ 의 임계값은 약 $3 \times 10^5$ , 경계층 두께 $\delta$ 로 정의한 $\mathrm{Re}_\delta$ 는 약 2700에서 흐트러진다. 같은 천이 현상이 형상에 따라 다른 무차원수로 다시 나타나는 것이다.

Prandtl이 본 경계층 — 외층과 내층#

레이놀즈가 임계 조건을 알린 1년 뒤(1904년), 프란틀은 다른 질문을 던졌다. 점성이 작아도 벽 가까이서는 왜 마찰력이 무시할 수 없는가?

답은 점성이 효과를 발휘하는 영역을 명시적으로 분리하는 것이었다. 그는 벽에 점착(no-slip)된 유체로부터 자유류 속도까지 가속되는 얇은 층을 경계층으로 명명하고, 그 바깥은 포텐셜 유동으로 다뤘다. 이 분할 덕분에 항공기 날개 위 2.5–25 mm 두께의 얇은 층만 보면 점성의 영향을 추적할 수 있게 됐다.

난류 경계층은 한 층이 아니다. 벽으로부터 거리 $y$ 를 적절히 무차원화하면 보편 구조가 드러난다.

y^+ = \frac{y\, u_\tau}{\nu}, \qquad u_\tau = \sqrt{\tau_w / \rho}

$u_\tau$ 는 벽면 마찰속도(wall friction velocity, 벽 전단응력 $\tau_w$ 를 밀도로 나눈 후 제곱근). $\nu / u_\tau$ 는 점성이 운동량을 운반하는 거리, 곧 점성길이다. 즉 $y^+$ 는 "점성길이로 잰 벽으로부터의 거리"다.

벽법칙 — 점성저층, 완화층, 로그 영역#

평판 난류 경계층 내부의 속도 분포를 $u^+ = u / u_\tau$ 로 무차원화하면 세 영역이 나타난다.

영역	범위	속도 분포
점성저층 (viscous sublayer)	$0 < y^+ < 5$	$u^+ = y^+$
완화층 (buffer layer)	$5 < y^+ < 30$	두 식의 부드러운 천이
로그 영역 (log-law / inertial sublayer)	$30 < y^+ < 0.2\, \delta^+$	$u^+ = \frac{1}{\kappa} \ln y^+ + B$

여기서 $\kappa \approx 0.41$ 은 폰 카르만 상수, $B \approx 5.0$ 은 매끈한 벽에 대한 상수.

벽 바로 옆에서는 와괴의 운동이 벽에 의해 억제되어 분자확산이 지배한다. 따라서 전단응력이 거의 일정한 영역에서 $\tau_w = \mu \, du/dy$ 를 적분하면 정확히 $u^+ = y^+$ 가 나온다. 두께는 풍속 20 m/s, 전연으로부터 2 m 떨어진 위치에서 $y^+ = 1$ 이 약 0.02 mm — 점성저층은 머리카락 굵기의 1/5다.

로그 영역은 프란틀의 혼합거리 이론(mixing length, 와괴가 주변과 동화될 때까지 이동하는 평균 거리)에서 따라 나온다. 벽 근방에서 혼합거리가 거리 $y$ 에 비례한다고 두면 $\ell = \kappa y$ , 일정 전단응력 가정과 결합해 적분하면 위 로그식이 자연스럽게 떠오른다.

아래 그래프에서 로그-법 상수 $B$ 를 바꿔보자.

Log-law constant B: 5.00 smooth buffer blend shade zones

$B$ 를 5.0에서 5.5로 올리면 로그 직선이 위로 평행 이동한다 — 매끄럽지 않은 벽에서는 $B$ 가 작아지고, 흡입(suction)이 있으면 커진다. 즉 로그 직선의 절편 한 숫자에 표면 상태가 압축돼 있다.

y+=30에 로그식만 대면 u+는 (1/0.41)×ln(30)+5.0이다. ln(30)은 약 3.40이고, 0.41로 나누면 약 8.30, B를 더해 약 13.3이다. 같은 높이의 선형식은 u+=30이라 둘 사이가 17 가까이 벌어진다. y+=5에서 선형식은 5다. 로그식을 그 자리에 억지로 넣으면 ln(5)가 약 1.61이라 u+는 약 8.9가 된다. 5와 8.9 사이의 빈칸이 완화층이다. 한 식으로 이으면 그 구간 속도가 이중으로 계산된다.

풍속 20 m/s, 전연에서 2 m인 조건에서 y+=1이 0.02 mm였으니 y+=5는 0.10 mm다. 선형 법칙이 사는 두께가 그 0.10 mm다. 벽함수가 쓰는 30에서 300까지는 y+가 열 배 늘 때 ln이 약 2.30 늘고, u+는 0.41로 나눠 약 5.6만 오른다. B를 5.0에서 5.5로 옮기면 로그층 전체가 0.5만큼 평행 이동한다. 그 0.5는 한 디케이드 상승 5.6의 약 9%다. 절편이 거칠기 때문에 0.5만 어긋나도, 벽함수가 내놓는 마찰속도가 같이 밀린다.

혼합거리 이론 — 와괴를 분자처럼#

층류에서 운동량을 운반하는 주체는 분자다. 뉴턴의 점성법칙 $\tau = \mu\, du/dy$ 가 그 결과다. 난류에서는 와괴(eddy, 다양한 크기의 조직화된 유체 덩어리)가 그 자리를 차지한다.

프란틀은 기체분자운동론에서 평균자유행정이 분자 마찰을 결정한다는 사실을 빌려와, 와괴의 평균 이동 거리를 혼합거리 $\ell$ 이라 명명했다. 변동속도 $u'$ 는 거리 $\ell$ 을 이동하면서 만난 평균속도 차로 근사된다.

\overline{u' v'} \approx -\ell^2 \left(\frac{d\bar{u}}{dy}\right)^2

이로부터 와점성(eddy viscosity) $\nu_t$ 가 정의된다.

\nu_t = \ell^2 \left|\frac{d\bar{u}}{dy}\right|

이 양은 벽 근방에서 분자 점성보다 $10^4$ 배 클 수 있다. 공기의 분자 점성 $\nu \approx 1.5 \times 10^{-5}$ m²/s에 비해 지표면 부근의 와점성은 0.01–0.1 m²/s, 대규모 대기 현상에서는 1–100 m²/s에 이른다.

혼합거리 이론은 가설이지만, 0차 RANS 모델인 혼합거리 모델, k-ε, k-ω의 사고 골격이 모두 여기에 닿아 있다. 와동을 분자처럼 다루겠다는 단순화가 1925년에 만들어졌고, 100년 후에도 산업용 CFD 솔버의 기본 옵션으로 살아 있다.

NumPy로 본 층류 vs 난류 속도분포#

원관 내 층류는 정확한 해가 있다 — 회전 포물면(Hagen–Poiseuille). 난류는 거듭제곱 법칙(power law)으로 근사하면 같은 무차원 단면에서 한참 더 평탄하다.

import numpy as np
import matplotlib.pyplot as plt
 
def laminar_poiseuille(r, R, U_max):
    """Hagen–Poiseuille: 회전 포물면 속도 분포."""
    return U_max * (1.0 - (r / R) ** 2)
 
def turbulent_power_law(r, R, U_centerline, n=7):
    """1/n 거듭제곱 난류 근사 (n=7은 Re ~ 1e5에서 잘 맞음)."""
    return U_centerline * (1.0 - np.abs(r) / R) ** (1.0 / n)
 
def law_of_the_wall(y_plus, kappa=0.41, B=5.0):
    """세 영역 결합: 점성저층 / 부드러운 완화층 / 로그 영역."""
    y_plus = np.asarray(y_plus, dtype=float)
    u_visc = y_plus
    u_log = (1.0 / kappa) * np.log(np.maximum(y_plus, 1e-12)) + B
    in_buffer = (y_plus > 5) & (y_plus < 30)
    t = np.clip((y_plus - 5) / 25.0, 0.0, 1.0)
    u_blend = u_visc * (1 - t) + u_log * t
    return np.where(y_plus <= 5, u_visc,
           np.where(y_plus >= 30, u_log, u_blend))
 
R = 1.0
r = np.linspace(-R, R, 201)
u_lam = laminar_poiseuille(r, R, U_max=1.0)
u_tur = turbulent_power_law(r, R, U_centerline=1.0, n=7)
 
bulk_lam = np.trapz(u_lam * np.abs(r), r) / np.trapz(np.abs(r), r)
bulk_tur = np.trapz(u_tur * np.abs(r), r) / np.trapz(np.abs(r), r)
print(f"laminar  bulk/centerline = {bulk_lam:.3f}")
print(f"turbulent bulk/centerline = {bulk_tur:.3f}")
# laminar  bulk/centerline = 0.500
# turbulent bulk/centerline = 0.817
 
y_plus = np.logspace(-1, 3, 200)
u_plus = law_of_the_wall(y_plus)
 
fig, axes = plt.subplots(1, 2, figsize=(10, 4))
axes[0].plot(r, u_lam, label="laminar")
axes[0].plot(r, u_tur, label="turbulent (n=7)")
axes[0].set(xlabel="r/R", ylabel="u / U_max", title="pipe profile")
axes[0].legend()
axes[1].semilogx(y_plus, u_plus)
axes[1].set(xlabel="y+", ylabel="u+", title="law of the wall")
plt.tight_layout()

핵심 결과는 두 줄이다. 층류의 평균유속은 중심선 속도의 절반 — 회전 포물면 부피와 외접 원기둥 부피의 비. 난류의 평균유속은 중심선 속도의 81%. 같은 평균유속을 얻으려면 난류는 더 낮은 중심선 속도로 충분하다는 뜻이고, 거꾸로 말하면 난류는 단면 전체에 운동량을 더 균질하게 퍼뜨린다. 이것이 와점성이 분자 점성보다 만 배 크다는 진술의 가시적 표현이다.

천이가 산업에서 의미하는 것#

층류와 난류의 갈림은 단순한 이론 문제가 아니다.

항공기 날개에서 층류 영역을 길게 끌면 마찰항력이 50% 가까이 줄어든다 — NLF(natural laminar flow) 설계의 핵심 동기.
열교환기는 난류를 의도적으로 일으켜 열전달계수를 높인다. 난류의 운동량 균질화가 곧 열의 균질화이기도 하기 때문.
파이프 송유는 임계 Re 부근에서 펌프 동력이 갑자기 커진다. 천이 영역을 피해 작동점을 잡는 것이 운영 비용을 좌우한다.
LES(Large Eddy Simulation)나 RANS 모델에서 첫 격자점이 $y^+ \approx 1$ 에 들어가야 점성저층을 해상하고, $30 < y^+ < 300$ 이면 벽법칙으로 보강한다. 한 자리 숫자가 격자 비용을 좌우한다.

첫 격자점을 로그층에 두면 점성저층이 빠진다#

첫 셀이 y+ 30 근처에 있으면 u+=y+인 구간은 계산에 없다. 벽 전단은 로그식의 절편 B에 기댄다. 매끈한 벽의 B는 5.0이다. 거친 벽에서 그 숫자를 그대로 쓰면, 위에서 본 0.5의 평행 이동이 마찰속도 쪽으로 넘어간다. 코드의 0.817은 n=7, Re 10만 근처의 매끈한 관 근사다. 같은 중심 속도면 난류 유량은 층류 0.500 대비 1.63배다. 절편이 틀린 관에 그 비를 붙이면 유량이 어긋난다.

그 비는 속도분포가 이미 발달했다는 전제에서만 성립한다. 입구부터 포물선, 혹은 1/7 거듭제곱이 되기까지 관이 얼마나 필요한지는 발달 길이와 Hagen–Poiseuille이 따로 센다. 입구 길이를 빼먹고 출구 단면만 보면, 층류 0.500과 난류 0.817의 비교가 아직 존재하지 않는 단면에 적용된다.