개수로 유동

유체역학 전산유체역학 마지막 수정: 2026-08-08 04:51:29

1. 개요[편집]

개수로 유동 Open channel flow
구동력	중력 (관로유동의 압력구배에 대응)
특징	자유수면 존재 — 수심 자체가 미지수
지배 무차원수	프루드 수
등류 공식	매닝 공식 · 셰지 식
대표 불연속	도수 (hydraulic jump)
수치 모델	얕은 물 방정식 유한체적 · HEC-RAS

관로 해석은 경계가 정해져 있고 유동을 푼다. 개수로 해석은 경계를 유동과 함께 푼다.

개수로 유동(open channel flow)은 위쪽이 대기에 열린 자유수면을 가진 흐름으로, 압력구배가 아니라 중력이 유동을 구동하고 수심 자체가 풀어야 할 미지수인 유동이다. 하천, 수로, 하수관거의 부분 만관 흐름, 도로 배수, 심지어 홍수터 범람까지 전부 여기에 들어간다.

관유동과의 차이는 단순히 “뚜껑이 없다”가 아니다. 관로에서는 단면적이 주어지고 압력강하가 결과지만, 개수로에서는 수면 위치가 곧 해다. 자유수면은 압력이 대기압으로 고정된 경계인데 그 위치를 모르므로, 경계가 해와 함께 움직이는 문제가 된다. 그 결과 지배 방정식이 비선형이 되고, 해가 여러 개 존재하며(대응수심), 불연속(도수)이 자발적으로 생긴다. 유체역학에서 자유수면 유동이 늘 골치인 이유의 압축판이다.¹

2. 흐름 상태 — 사류와 상류[편집]

개수로에서는 정보가 흐름을 타고 하류로도 가지만, 중력파(표면파)를 타고 상류로도 거슬러 올라간다. 얕은 물에서 미소 진폭 표면파의 속도는 $\sqrt{gD}$ ( $D=A/T$ , 통수 단면적을 수면폭으로 나눈 수리수심)이므로, 유속과 이 파속의 비인 프루드 수 $\mathrm{Fr}=V/\sqrt{gD}$ 가 흐름의 성격을 가른다.

$\mathrm{Fr}<1$ 상류(常流, subcritical) — 파가 상류로 거슬러 올라갈 수 있다. 하류 조건이 상류 수면을 지배한다. 깊고 느린 하천.
$\mathrm{Fr}=1$ 한계류(critical) — 임계수심 $y_c$ . 직사각 수로에서 단위폭 유량 $q$ 에 대해 $y_c=(q^2/g)^{1/3}$ .
$\mathrm{Fr}>1$ 사류(射流, supercritical) — 파가 하류로만 간다. 상류 조건만이 해를 결정한다. 얕고 빠른 여수로.

이 구조는 압축성 유동의 아음속/초음속과 수학적으로 동형이고, 계산 관점에서 중요한 것은 경계조건을 어디에 걸어야 하는가가 상태에 따라 뒤집힌다는 점이다. 상류 구간에서는 하류 끝에서, 사류 구간에서는 상류 끝에서 조건을 준다. 이걸 반대로 걸면 솔버가 조용히 이상한 답을 낸다.

3. 비에너지와 대응수심[편집]

수로 바닥을 기준으로 잰 에너지를 비에너지(specific energy)라 한다.

E = y + \frac{V^2}{2g} = y + \frac{Q^2}{2gA(y)^2}

유량 $Q$ 를 고정하고 $E$ 를 $y$ 의 함수로 그리면 특유의 곡선이 나온다. $y$ 가 크면 두 번째 항이 사라져 $E \to y$ (45° 점근선), $y$ 가 작으면 속도수두가 폭발해 $E\to\infty$ . 따라서 곡선은 최솟값을 갖는 U자 형태이고, 하나의 $E$ 값에 대해 수심이 두 개 존재한다. 이 둘을 대응수심(alternate depths)이라 부르며, 하나는 상류(깊고 느림), 하나는 사류(얕고 빠름)다.

$dE/dy=0$ 을 풀면 정확히 $\mathrm{Fr}=1$ 이 나온다. 즉 최소 비에너지 상태가 임계 상태이며, 직사각 수로에서는 $E_{\min}=\tfrac{3}{2}y_c$ 다. 이 사실이 실무에서 강력한 이유는, 수로가 좁아지거나 바닥이 솟아 흐름이 임계에 도달하는 지점이 생기면 그 단면이 유량과 수심을 묶는 제어단면(control section) 이 되기 때문이다. 위어, 여수로 마루, 파샬 플룸 같은 계측 구조물이 전부 이 원리로 유량을 잰다 — 임계단면 하나만 만들어 놓으면 수심 하나로 유량이 결정된다.

4. 등류 — 매닝 공식[편집]

바닥 경사 $S_0$ 가 일정하고 단면이 변하지 않는 긴 수로에서는 중력의 구동과 벽면 마찰이 균형을 이뤄 수심이 일정해진다. 이때의 수심이 정상수심(normal depth) $y_n$ 이다. 셰지(1770년대)는 $V=C\sqrt{R_h S}$ 를 제시했고, 오늘날 현장을 지배하는 것은 매닝 공식이다(SI 단위).

V = \frac{1}{n} R_h^{2/3} S_0^{1/2}, \qquad R_h = \frac{A}{P}

$R_h$ 는 동수반경(통수 단면적 ÷ 윤변), $n$ 은 조도계수다. 셰지 계수와는 $C=R_h^{1/6}/n$ 로 연결된다.

매닝 공식의 실체는 완전난류 영역에서 마찰계수가 레이놀즈수에 무관해진다는 사실에 기댄 경험식이다. 그래서 $n$ 이 모든 것을 결정하는데, 이 값을 고르는 일이 곧 하천 해석의 예술 영역이다.

수로 상태	매닝 n
마감 콘크리트	0.011 ~ 0.014
흙수로, 잘 정비됨	0.020 ~ 0.025
자연 하천, 잡초·사행	0.033 ~ 0.050
산지 하천, 큰 호박돌	0.040 ~ 0.070

수치 모델의 보정(calibration)이 사실상 “관측 수위에 맞을 때까지 $n$ 을 조정하는 작업”으로 귀결되는 일이 흔하다. 그래서 개수로 수치해석에서는 격자 수렴보다 $n$ 불확실성이 훨씬 큰 오차원인 경우가 많다.²

5. 점변류와 수면형[편집]

실제 하천에서 수심은 좀처럼 $y_n$ 으로 가만히 있지 않는다. 댐이 막고, 낙차가 있고, 단면이 변한다. 수심이 거리에 따라 완만히 변하는 흐름을 점변류(gradually varied flow, GVF)라 하고, 에너지식을 미분해 정리하면 딱 한 줄의 상미분방정식이 나온다.

\frac{dy}{dx} = \frac{S_0 - S_f}{1 - \mathrm{Fr}^2}

$S_f$ 는 매닝 공식으로 계산한 마찰경사다. 이 한 식에 개수로 수면 거동이 전부 들어 있다.

분자 $S_0-S_f$ 의 부호는 $y$ 가 $y_n$ 보다 큰가 작은가로 결정된다.
분모 $1-\mathrm{Fr}^2$ 의 부호는 $y$ 가 $y_c$ 보다 큰가 작은가로 결정된다.
$y \to y_c$ 에서 분모가 0이 되어 $dy/dx \to \infty$ — 수면이 수직이 된다. GVF 가정 자체가 깨지는 지점이며, 여기서부터는 급변류로 넘어간다.

$y_n$ 과 $y_c$ 의 대소로 수로를 완경사 M(mild, $y_n>y_c$ ), 급경사 S(steep), 임계 C, 수평 H, 역경사 A로 나누고, 실제 수심이 $y_n$ · $y_c$ 가 만드는 세 구역 중 어디에 있는지로 번호를 붙인다. 그래서 M1, M2, M3, S1, S2, S3 … 같은 수면형 분류가 나온다.

M1 — 배수곡선(backwater). 완경사 하천에 댐이나 보를 세우면 상류로 수위 상승이 수 km씩 뻗는다. 홍수 위험도 평가의 핵심이며, “댐 하나 세우면 어디까지 잠기는가”가 바로 이 곡선이다.
M2 — 저하곡선(drawdown). 자유 낙차(free overfall) 직전에서 수면이 임계수심을 향해 내려간다.
M3 — 수문 아래로 빠져나온 사류가 하류 상류 흐름을 만나기까지의 구간. 이 구간의 끝에서 도수가 터진다.

이 ODE를 푸는 것이 하천 모델링의 본체다. 단면이 균일한 인공수로에서는 $y$ 를 등간격으로 잡고 $x$ 를 구하는 직접축차법(direct step method)이, 단면이 제각각인 자연하천에서는 $x$ 를 단면 위치로 고정하고 $y$ 를 반복 수렴시키는 표준축차법(standard step method)이 표준이다. HEC-RAS의 정상류 계산이 바로 후자다.

6. 도수 — 에너지로는 못 푸는 문제[편집]

사류가 상류로 전이할 때는 완만한 변화가 불가능하다. 수면이 급격히 솟구치며 거센 난류 롤러가 형성되는데, 이것이 도수(hydraulic jump)다.

여기서 초보자가 반드시 한 번 밟는 지뢰가 있다. 도수를 에너지 방정식으로 풀려고 하는 것. 도수 구간에서는 난류로 흩어지는 에너지 손실이 상당한데, 그 손실량을 미리 모르므로 에너지식은 미지수 두 개짜리 식 하나가 되어 닫히지 않는다. 대신 운동량을 쓴다. 짧은 구간이라 벽면 마찰력이 무시 가능하고, 그러면 운동량 플럭스와 정수압력의 합인 비력(specific force)이 보존된다.

M = \frac{Q^2}{gA} + \bar z A = \text{const}

직사각 수로에 적용하면 상류 프루드 수 $\mathrm{Fr}_1$ 만으로 도수 전후의 공액수심(conjugate depths) 비가 나온다.

\frac{y_2}{y_1} = \frac{1}{2}\left(\sqrt{1+8\,\mathrm{Fr}_1^2}-1\right)

그리고 에너지 손실은 계산 결과로 따라 나온다: $\Delta E = (y_2-y_1)^3/(4y_1y_2)$ . 항상 양수이며, $\mathrm{Fr}_1$ 이 커질수록 급격히 커진다. 여수로 끝에 감세공(stilling basin)을 두어 일부러 도수를 일으키는 것은 이 손실을 활용하는 설계다 — 하류 하상을 파먹을 에너지를 물속에서 미리 태워 버린다.³

구조적으로 이것은 충격파와 완전히 같다. 압축성 유동에서 랭킨-위고니오 관계를 질량·운동량·에너지 보존으로 유도하고 엔트로피 증가로 방향을 정하듯, 도수도 보존형에서 유도하고 에너지 감소로 방향을 정한다. 대응수심(에너지 곡선의 두 해)과 공액수심(운동량 곡선의 두 해)이 서로 다른 쌍이라는 점도 반드시 구분해야 한다.⁴

7. 수치해석[편집]

개수로 수치 모델의 표준 지배식은 얕은 물 방정식(생브낭 방정식)이다. 쌍곡형 보존법칙계이므로 유한체적법 + 리만 솔버(HLL, HLLC, 로) 조합이 정석이고, 도수는 별도 처리 없이 충격 포착으로 자연스럽게 잡힌다. 보존형으로 풀어야 공액수심이 정확히 나온다는 점도 충격파와 동일하다.

이 분야만의 고약한 문제 두 가지가 있다.

웰밸런스드 도식(well-balanced). 정지한 호수( $u=0$ , 수면 수평, 바닥은 울퉁불퉁)는 자명한 정상해여야 한다. 그런데 플럭스 항과 바닥경사 소스항을 따로 이산화하면 둘이 정확히 상쇄되지 않아 가만히 있는 물에서 가짜 유속이 생긴다. 지형 기복이 큰 실제 하천에서는 이 오차가 실제 신호보다 커질 수 있어서, 수면고도 기반 재구성(hydrostatic reconstruction) 같은 웰밸런스드 처리가 사실상 필수다.
마름/젖음(wet-dry). 홍수터 범람은 수심 0인 셀이 젖는 과정이다. $h\to 0$ 이면 $u=hu/h$ 가 0으로 나누기가 되고 파속 $\sqrt{gh}$ 로 잡은 시간 간격이 무너진다. 임계 수심 아래를 마른 셀로 처리하고 플럭스를 특수 처리하는 로직이 붙는데, 여기서 질량이 새는 버그가 이 바닥 단골 이슈다.

실무 코드는 목적에 따라 갈린다. 1차원 단면 기반 정상류·부정류 계산에는 HEC-RAS(미 육군공병단, 무료라 사실상 표준), 도시 우수관망과 지표 흐름 연동에는 SWMM, 상용 통합 모델로는 MIKE·Delft3D·TELEMAC이 쓰인다. 최근에는 2차원 침수 해석이 기본이 되어 GPU 가속 SWE 솔버가 빠르게 자리를 잡았다.

프루드 상사가 지배한다는 점에서 개수로는 자유수면 유동 일반과 한 가족이다. 선박 조파저항 모형시험, 탱크 내 슬로싱, 댐 붕괴 문제 모두 프루드 수를 맞춰 축소모형을 만들고 — 그 대가로 레이놀즈수는 포기한다.

8. 관련 문서[편집]

그래서 개수로 문제를 3차원 VOF 방법으로 통째로 푸는 일은 여수로 국부 해석 정도에서나 한다. 하천 수십 km를 그렇게 풀면 계산기가 먼저 항복한다. 연직 방향을 적분해 없앤 얕은 물 방정식이 압도적으로 실용적인 이유. ↩
조도계수 표는 채우(Ven Te Chow)의 1959년 교과서 표가 60년 넘게 인용되고 있다. 사진과 함께 “이런 하천은 n이 얼마”라고 적어 놓은 그 표 말이다. 수치해석 논문에서 5차 정확도 도식을 자랑해도, 입력 $n$ 의 불확실성이 ±30%면 그 정확도는 장식이 된다. ↩
도수의 형태도 $\mathrm{Fr}_1$ 에 따라 이름이 다르다. 1.0~~1.7 파상도수, 1.7~~2.5 약도수, 2.5~~4.5 진동도수, 4.5~~9.0 정상도수, 9.0 초과 강도수. 설계자가 가장 싫어하는 것은 진동도수인데, 에너지 소산은 시원찮으면서 주기적인 파를 하류로 계속 내보내 호안을 갉아먹는다. ↩
사류 구간에 다리 교각처럼 장애물을 두면 수면에 사선 형태의 사파(oblique standing wave) 가 생긴다. 초음속 유동의 사각충격파와 판박이라, 도수를 “수직 충격파”, 사파를 “사각 충격파”에 대응시켜 배우면 두 과목을 한 번에 정리할 수 있다. ↩