빗방울은 왜 계속 빨라지지 않을까?
중력은 낙하하는 내내 같은 힘으로 끌어당긴다 — 질량 1kg당 약 , 1초가 지나든 1분이 지나든 변함없이. 진공이라면 자유낙하는 끝없이 빨라진다 — 포물선 운동에서 다룬 저항 없는 풍경. 그런데 공기 중에서 빗방울은 약 9 m/s, 낙하산을 편 사람은 약 5 m/s, 먼지 알갱이는 거의 가만히 떠 있다. 누구도 뉴턴 법칙을 어기지 않는다 — 두 힘이 정확히 상쇄되는
종단속도는 최대값이 아니다 — 평형이다. 중력과 저항이 상쇄되는 점근선.
속도는 미분, 가속도는 또 한 번의 미분
낙하하는 물체에는 위치 가 있다. 미분 모듈에서 본 대로,
두 힘 — 중력은 상수, 저항은 아니다
낙하하는 빗방울에 두 힘이 작용한다. 중력은 의 일정한 크기로 아래로 당긴다. 단위 질량당 그저 , 빠르기와 무관하다. **
로 들어가는 단위 질량당 알짜힘은 그 차이:
두 항. 첫 항은 상수. 둘째 항은 아니다. 그게 식 전부다.
dv/dt = g − kv가 실제로 말하는 것
왼쪽에서 오른쪽으로 읽자: 어느 순간이든, 속도가 변하는 비율 = 중력 − 저항. 재미있는 점은 저항이 v 자체에 달려 있다는 것 — v가 어떻게 변할지 결정하는 우변에 v가 또 들어 있다. 그래서 이 식은 미분방정식이다. 미지의 함수와 그 미분이 한 식에 같이 등장한다.
식만 보고도 세 영역을 읽을 수 있다:
- 가 작을 때 (정지에서 막 떨어진 직후): 가 작으니 — 진공에서의 자유낙하처럼 가속이 그대로 살아 있다.
- 가 중간일 때: 가 무시할 수 없는 크기. 는 여전히 양수지만 작아진다. 속도는 자라되 점점 느리게.
- 일 때: . . 속도가 더 이상 변하지 않는다. 이게
.종단속도 - v > g/k (예: 아래로 세게 던진 경우): kv > g. dv/dt < 0. 속도가 줄어든다 — 저항이 이긴다. 위쪽에서 로 내려올 때까지.
위젯의 파란 곡선은 점선 초록선을 아래에서부터 따라붙는다. 어느 쪽에서 출발하든 가 종착점.
import math
# Per-unit-mass equation of motion: dv/dt = g - k*v.
# Two competing forces, both per unit mass:
# gravity → +g (constant, pulls v upward in magnitude)
# drag → -k*v (proportional to current speed, opposes motion)
g = 9.8
def dvdt(v, k):
return g - k * v
# At v = 0, drag is zero; gravity is unopposed; full acceleration.
# As v grows, drag grows proportionally; net force shrinks; acceleration
# shrinks. The shrink is the whole point — the system has a built-in
# governor.
[(v, round(dvdt(v, k=0.5), 2)) for v in (0, 5, 10, 15, 19.6, 25)]
# → [(0, 9.8), free fall — nothing opposing
# (5, 7.3), drag has eaten 2.5 of g
# (10, 4.8),
# (15, 2.3),
# (19.6, 0.0), EQUILIBRIUM — v_t exactly
# (25, -2.7)] above v_t: drag wins, deceleration종단속도 = dv/dt = 0 풀기
“종착점”은 그저 이 되도록 만드는 값이다. 식을 세우고 풀면:
대수 한 줄. 저항이 클수록 (큰 ) 종단속도는 작아진다 — 깃털은 작은 속도에서 일찍 안정되고, 볼링공은 공기 중에서 에 닿는 데 영원에 가까운 시간이 걸린다 (실제로는 그 전에 땅에 부딪힌다). 이것이
식 는 깊은 정리가 아니다 — 고정점을 그대로 읽었을 뿐. 어떤 양에 그 양 자체에 따라 변하는 두 영향이 맞서고 있을 때면 어디서나 똑같은 모양이 나타난다. 시원한 방에서 식어가는 커피는 방 온도에 가까워지고, 저항을 거쳐 충전되는 축전기는 공급 전압에 가까워지고, 크기에 비례하는 출생률·사망률을 가진 인구는 환경수용력 (carrying capacity) 에 가까워진다. 물리는 달라도, 찾는 건 같은 고정점이다.
# Solve dv/dt = 0 for terminal velocity. One line of algebra.
# g - k*v_t = 0 → v_t = g/k.
def terminal_velocity(k, g=9.8):
return g / k
# Linear-drag terminal velocities for a few familiar k values.
[(name, k, round(terminal_velocity(k), 1))
for name, k in (("feather", 6.0), ("raindrop", 1.6),
("skydiver", 0.4), ("bowling ball", 0.05))]
# → [('feather', 6.0, 1.6), reaches ~1.6 m/s in well under a second
# ('raindrop', 1.6, 6.1),
# ('skydiver', 0.4, 24.5), classic ~50 mph free-fall figure
# ('bowling ball', 0.05, 196.0)] approaches very slowly, in practice
# it'd hit the ground long before
# Real raindrops sit closer to 9 m/s because actual drag is closer to
# v² than v; the linear case is a teaching simplification.다가감은 점근적 — 그리고 지수적
정지에서의 닫힌 형식 해는
종단속도와의 오차 — 격차 — 는 로 지수적으로 줄어든다. 초가 지나면 오차가 절반. 초가 지나면 오차는 의 1% 아래로 떨어진다. 위젯에서 그대로 보인다. 저항 슬라이더를 위로 (큰 ) 올리면 곡선이 1초도 안 돼 점근선에 들러붙고, 아래로 내리면 몇 초에 걸쳐 천천히 올라간다.
왜 지수일까? 오차 를 미분해 보면 . 오차가 변하는 비율이 오차 자체에 비례하고, 부호는 음수다. 지수 감쇠의 전형적인 모양. 자기 자신에 비례하는 비율로 변하는 양은 결국 모두 지수가 된다 — 로그 모듈에서는 같은 사실을 반대편에서 부른다. 선형 저항 낙하, 축전기 충전, 방사능 붕괴, 뉴턴 냉각 법칙 — 같은 식, 같은 답.
# Closed-form solution from rest, by separation of variables.
# v(t) = v_t * (1 - exp(-k*t))
# At t = 0, v = 0 (rest). As t → ∞, v → v_t. Half of v_t is reached at
# t = ln(2)/k ≈ 0.693/k. After 5 'time-constants' (5/k seconds), v is
# within 1% of v_t.
def v_of_t(t, k, g=9.8):
return (g / k) * (1 - math.exp(-k * t))
[(t, round(v_of_t(t, k=0.4), 2)) for t in (0, 1, 3, 5, 10, 20)]
# → [(0, 0.0),
# (1, 8.1), a third of v_t
# (3, 17.1), two-thirds
# (5, 21.2), ~87% of v_t = 24.5
# (10, 24.0), ~98%
# (20, 24.5)] within 0.001
# Same shape as a charging capacitor, a heating cup of coffee, a leaky
# bucket — every first-order linear ODE with a stable fixed point traces
# this curve. The terminal value is the fixed point; k is the rate at
# which deviation from the fixed point decays.종단속도는 결국 얼마나 빠르게 떨어지는지를 알려줄 뿐이다. 다음 질문은 자연스럽게 t초 동안 얼마나 떨어졌는가? — 이건 적분이 답한다. 를 알면 거리는 속도가 쌓인 양이 된다: . 0부터 까지 속도 곡선 아래의 넓이가 그때까지 떨어진 거리다. 를 대입하면 닫힌 형식 — — 이 나온다. 식에서 두 가지 영역이 그대로 읽힌다. 작은 에서는 곡선이 직선 아래에 놓인다 (처음부터 종단속도로 떨어진 게 아니니까). 큰 에서는 격차가 더 자라지 않고 가 에 가까워진다 — 종단속도 × 시간에서 “종단속도까지 가속하는 데 든 시간”을 빼는 상수 오프셋이 남는 셈.
실제 빗방울은 이차 저항을 따른다. : 힘이 속도의 1차가 아니라 제곱에 비례해 자란다. 종단속도를 구하는 논리는 그대로다 — 을 풀어 , — 다만 다가가는 모양은 더 이상 깔끔한 지수가 아니고 수치도 달라진다. 선형 저항의 는 저속·소형 (낮은 레이놀즈 수) 에 맞춘 교육용 단순화일 뿐. 답의 형태 (힘 균형 자리의 점근선) 는 어떤 저항 법칙에서도 똑같지만, 닫힌 형식은 달라진다.
종단속도는 최대값이 아니다 — 평형이다. 중력은 일정한 힘, 저항은 속도에 따라 자라는 힘. 둘이 균형을 이루는 자리에서 , 속도는 더 이상 변하지 않는다. 식 가 “낙하는 왜 평탄해질까?”라는 물음을 ” 의 고정점은 어디인가?”로 바꿔 놓고, 답은 대수 한 줄로 끝난다.
위젯에서 빗방울 프리셋 () 을 눌러 보자. 점선 점근선에서 계산 없이 종단속도를 눈으로 읽어 본 뒤, 와 로 식에서 다시 확인한다.
, 일 때 종단속도는? 이번에는 를 두 배 () 로 키우면 는 어떻게 바뀌나? 절반 () 으로 줄이면?
, 일 때 종단속도는 . 세 가지 속도, , , 에서 를 계산하라. 각각의 부호는 무엇을 말해 주는가?
한 학생이 종단속도로 떨어지는 빗방울을 보고 이런 결론을 내린다: “방울이 가속하지 않으니까, 중력은 더 이상 작용하지 않는 거다.” 한 문장으로 반박하라. 올바른 그림은 무엇인가?
길이가 같은 두 구간, 둘 다 다 — 하나는 낙하 초반 (에서 ), 다른 하나는 후반 (에서 , 그쯤이면 는 사실상 ). , . 각 구간에서 떨어진 거리를 계산하라. 같은 시간인데 후반 구간이 더 많이 떨어지는 이유는?
입문 물리는 다루기 쉬운 자유낙하 (저항 없음) 만 가르친 뒤, “실제 낙하에는 공기 저항이 있어서 복잡하다”라고 손사래를 친다. 그 복잡함은 식 하나에 항 하나가 더 붙는 것뿐이다: . 이 한 줄에서 종단속도는 대수 한 줄 () 로 나오고, 다가가는 모양은 적분 한 줄 () 로, 그리고 충전되는 축전기·식어가는 커피와의 연결은 한 문단으로 끝난다. Lemma는 식을 독자 눈앞에 그대로 둔다 — 물리는 거기서 진짜로 일어나기 때문이다. 교과서 사이드바를 채우는 “흔한 종단속도 표” 같은 데가 아니라. 종단속도는 빗방울에 관한 사실이 아니다. 안정한 고정점을 가진 모든 1계 시스템에 관한 사실이다.