거친 공이 더 멀리 날아가는 이유 — 항력 위기와 경계층 박리

거칠기 0.85는 비행점에서 Cd를 0.116으로 내리지만, 1.0까지 올리면 같은 70 m/s에서 0.124로 되오른다.

같은 속도면 매끈한 공이 공기를 덜 가른다고 생각하기 쉽다. 골프공은 그 생각을 부순다. 표면에 딤플(dimple, 작은 홈)이 수백 개인데도 매끈한 공보다 거의 두 배 멀리 간다. 항력계수 $C_D$ 가 어느 Reynolds 수에서 절벽처럼 떨어지는 구간, 항력 위기(drag crisis)가 그 차이다.

경계층 버튼을 층류에서 난류로 넘기면 박리점이 뒤로 밀리고, 그늘진 후류의 폭이 줄어든다.

Toggle the boundary layer. The turbulent one carries momentum closer to the wall, separates later, and leaves a thinner wake — the shaded region behind the cylinder shrinks and the pressure drag drops.

난류 쪽에서 후류가 얼마나 좁아지는지 보자. 폭이 줄면 압력 항력은 절반 아래로 떨어진다. 더 어지러운 경계층이 항력을 깎는 장면이다.

압력 항력과 마찰 항력#

공이 받는 힘은 두 갈래다. 표면을 스치는 점성이 만드는 마찰 항력(skin friction). 앞뒤 압력차가 만드는 압력 항력(form drag).

공처럼 뭉툭한 물체(bluff body)에서는 압력 항력이 거의 전부다. 앞면은 공기가 부딪혀 고압이다. 뒷면은 흐름이 떨어진 자리라 저압이 남는다. 그 차이가 공을 뒤로 당긴다.

전체를 계수 하나로 묶으면 이렇다.

$F_D = \tfrac{1}{2}\, \rho\, U^2\, A\, C_D$

$\rho$ 는 밀도, $U$ 는 상대속도, $A$ 는 정면 투영 면적, $C_D$ 는 항력계수다. 골프공이 건드리는 것은 $C_D$ 다. 그 값은 Reynolds 수에 따라 평탄하다가 갑자기 떨어진다.

$Re = \frac{\rho\, U\, D}{\mu}$

$D$ 는 지름, $\mu$ 는 점성계수다. $Re$ 는 관성력과 점성력의 비다.

경계층이 벽을 떠나는 순간#

압력 항력의 크기는 흐름이 언제 벽을 떠나느냐에 달렸다. 경계층 박리(separation)다.

공의 앞쪽 절반에서 공기는 빨라지고 압력은 내려간다. 뒤쪽 절반은 반대다. 느려지면서 압력이 다시 오른다. 역압력구배(adverse pressure gradient, 흐름 방향으로 압력이 커지는 상태)다.

벽 근처 공기는 이미 점성으로 운동량을 많이 잃었다. 이 맞바람을 끝까지 버티지 못한다. 멈추고, 거꾸로 밀린다. 그 점에서 경계층이 떨어지고, 뒤에는 넓은 저압 후류(wake)가 남는다.

박리가 이를수록 후류가 넓다. 후류가 넓을수록 압력 항력이 크다.

라미나르 경계층과 난류 경계층의 분업#

박리 각은 경계층이 층류인지 난류인지로 갈린다.

층류 경계층은 층끼리 섞이지 않는다. 벽 근처는 느리고 운동량이 빈약해서, 역압력구배를 만나면 일찍 무릎을 꿇는다. 앞쪽 정체점에서 약 82°다.

난류 경계층은 바깥의 빠른 공기를 벽 쪽으로 끌어온다. 벽 근처가 더 단단하다. 같은 맞바람을 더 오래 버틴다. 박리는 약 120°까지 밀리고 후류는 좁아진다. 위 버튼에서 본 폭의 차이가 이 각도 차이다. 벽 근처가 왜 버티는지는 관의 층이 갈라지는 구조를 같이 봐야 숫자가 된다.

항력 위기 — $C_D$ 가 절벽처럼 떨어진다#

매끈한 구의 $C_D$ 는 $Re \sim 10^3$ 에서 $3{\times}10^5$ 까지 약 0.4–0.5에 머문다.

$Re \approx 3{\times}10^5$ 부근에서 0.1 아래로 떨어진다. 흐름 방향이 바뀐 게 아니다. 경계층이 층류에서 난류로 천이(transition)하면서 박리점이 뒤로 밀린 것이다. 좁아진 후류가 계수를 끌어내린다. 이 급락이 항력 위기다.

거칠수록 천이가 더 낮은 $Re$ 에서 시작된다. 딤플은 일부러 거칠게 만든 표면이다. 골프공의 비행 $Re$ (약 $2{\times}10^5$ )는 매끈한 공이라면 아직 위기 앞, $C_D$ 가 높은 쪽이다. 딤플은 그 절벽을 앞당겨, 비행 속도에서 이미 낮은 $C_D$ 를 쓰게 만든다. 곡선을 손으로 옮기는 그림은 숫자 출력 바로 다음에 둔다.

Python — 골프공의 항력을 위기 앞뒤로#

매끈한 공과 딤플 공을 같은 속도에 넣는다. 지름 0.0427 m, 속도 70 m/s.

import math
 
def cd_subcritical(Re):
    # Clift-Gauvin 상관식 (Re < 3e5에서 유효)
    return (24.0 / Re) * (1 + 0.15 * Re**0.687) + 0.42 / (1 + 42500 * Re**-1.16)
 
def smoothstep(x, a, b):
    t = max(0.0, min(1.0, (x - a) / (b - a)))
    return t * t * (3 - 2 * t)
 
def drag_coefficient(Re, roughness):
    # roughness 0(매끈)~1(딤플): 임계 Re를 왼쪽으로 옮긴다
    sub = cd_subcritical(Re)
    log_re = math.log10(Re)
    log_crit = 5.55 - 1.35 * roughness
    t = smoothstep(log_re, log_crit - 0.18, log_crit + 0.18)
    sup = 0.08 + 0.13 * smoothstep(log_re, log_crit, 7.0)  # 위기 이후 가지
    return sub * (1 - t) + sup * t
 
def drag_force(U, D, roughness, rho=1.2, mu=1.8e-5):
    Re = rho * U * D / mu
    Cd = drag_coefficient(Re, roughness)
    A = math.pi * (D / 2) ** 2          # 정면 투영 면적
    Fd = 0.5 * rho * U**2 * A * Cd
    return Re, Cd, Fd
 
D = 0.0427   # 골프공 지름 [m]
U = 70.0     # 드라이버 직후 속도 [m/s]
 
for label, rough in [("매끈한 공", 0.0), ("딤플 공", 0.85)]:
    Re, Cd, Fd = drag_force(U, D, rough)
    print(f"{label}: Re={Re:.2e}  C_D={Cd:.3f}  F_D={Fd:.3f} N")

출력:

매끈한 공: Re=1.99e+05  C_D=0.487  F_D=2.049 N
딤플 공: Re=1.99e+05  C_D=0.116  F_D=0.488 N

roughness = 0.00Re marker = 1e5.0

cyan: sphere C_D · green: Stokes 24/Re. Raise roughness and the drag-crisis cliff slides toward lower Re — a dimpled ball trips the boundary layer early.

roughness 슬라이더를 올리면 절벽이 왼쪽으로 이동한다. logRe 마커를 5.3 근처에 두면 매끈한 곡선은 아직 평지고, 거친 곡선은 이미 바닥이다. 출력의 0.487과 0.116이 그 두 점이다.

출력의 2.049 N과 0.488 N은 같은 동압에서 갈라진다. $0.5 \times 1.2 \times 70^{2} = 2940$ Pa다. 지름 0.0427 m의 투영 면적은 약 $1.432 \times 10^{-3}$ m²이고, 동압과 곱하면 4.210 N이다. 0.487을 곱하면 2.049 N, 0.116을 곱하면 0.488 N이라 비는 4.20이다.

임계 로그는 $5.55 - 1.35 \times$ roughness다. 0이면 5.55, $Re$ 는 약 $3.55 \times 10^{5}$ 다. 0.85이면 4.40, $Re$ 는 약 $2.53 \times 10^{4}$ 다. 비행 $Re = 1.99 \times 10^{5}$ 의 로그는 5.299다. 매끈한 공의 천이 창은 5.37에서 열리므로 5.299는 창 밖이고, 아임계 가지가 0.487을 낸다. 딤플은 창을 지나 초임계 가지에 있다. 그 가지는 0.08에서 시작해 로그 7에서 0.21이 된다. roughness 0.85에서 smoothstep은 0.275이고, $0.08 + 0.13 \times 0.275 = 0.116$ 으로 출력과 맞는다.

roughness를 1까지 밀면 같은 5.299에서 smoothstep이 0.341이 되고 $C_D$ 는 0.124, 힘은 약 0.52 N으로 되오른다. 딤플 임계의 한가운데, $Re$ 약 $2.53 \times 10^{4}$ 는 8.9 m/s다. 70 m/s는 그 위에 있다. 공이 거기까지 느려지면 거친 표면도 높은 $C_D$ 쪽으로 돌아간다.

딤플, 그리고 다른 무대들#

같은 각도의 문제가 다른 장비에도 있다. 자동차는 후미를 다듬어 박리를 늦추고 후류를 줄인다. 크리켓은 공의 한쪽만 닦아 양쪽 박리점을 어긋나게 만든다. 실린더에 트립 와이어(trip wire, 경계층을 난류로 만드는 가는 철사)를 감으면 위기가 앞당겨진다.

공통점은 마찰이 직접 깎는 양보다, 흐름이 언제 벽을 떠나느냐가 항력을 정한다는 것이다.

거칠기를 더 올리면 항력이 되오르는 지점#

70 m/s는 이미 위기 뒤다. 거기서 roughness를 0.85보다 더 주면 $C_D$ 는 0.116에서 0.124로 되오른다. 출력의 4.20배 간격은 슬라이더의 끝이 아니다.

딤플을 열이 드나드는 면에 붙이면 벽 마찰도 같이 커진다. 그 마찰이 열전달로 어떻게 환산되는지는 항력만 재고 열을 맞히는 상사가 묻는다. 압력 항력을 줄이려고 천이를 앞당긴 표면에, 열 청구서가 따로 도착하는가?