25점 논술 출제빈도 ★★★★★ Module 2. 수계소화설비 & 소방수리학

Hazen-Williams 공식과 Darcy-Weisbach 공식의 수리학적 유도, 차이점 및 배관망 해석

1 개요 (Introduction & 패러다임 변화)

수계소화설비(스프링클러, 옥내소화전, 물분무, 미분무설비 등)의 수리학적 설계(Hydraulic Design)에서 배관 마찰손실수두($h_L$)의 산정은 소방펌프의 설계 양정($H$), 배관 구경의 경제성, 그리고 최원격 헤드의 최소 방수압력 및 방수량 확보를 좌우하는 소방수리학의 핵심 근간이다.

배관 마찰손실을 계산하는 대표적인 수식으로는 차원 해석과 물리적 운동량 보존 법칙에 기반한 이론식인 Darcy-Weisbach 공식과, 상온의 물 유동에 대한 광범위한 실험 데이터를 바탕으로 정립된 경험적 공식인 Hazen-Williams 공식이 있다. 국내 화재안전성능기준(NFPC 102/103) 및 NFPA 13은 실무적 간편성을 이유로 일반 수계설비의 수리계산에 Hazen-Williams 공식을 주로 적용해 왔으나, 이 공식은 물 이외의 유체(부동액, 포소화약제 수용액), 극고온/극저온수, 미분무설비의 초고압·초고속 영역, 또는 층류·천이 영역에서는 점성 계수의 결여로 인해 심각한 계산 오차를 유발한다.

최근 소방방재 엔지니어링 패러다임은 과거의 단순 규약배관방식(Pipe Schedule System)에서 탈피하여 컴퓨터 프로그램을 활용한 전양정 완전 수리계산방식(Hydraulically Designed System)으로 전환되었으며, 초고층·대공간 및 복합건축물의 복잡한 루프(Loop) 및 그리드(Grid) 격자 배관망을 정확히 해석하기 위해 Hardy Cross 수렴해법 및 Darcy-Weisbach 보편 유동 모델의 융합 적용이 소방기술사에게 강력히 요구되고 있다.

2 핵심 메커니즘 및 수식 유도

[메커니즘 1] 운동량 보존 법칙에 기반한 Darcy-Weisbach 공식 완전 유도

직경 $D$, 길이 $L$인 수평 원관 내에서 비압축성 유체가 일정한 단면평균유속 $v$로 정상 등속 유동할 때, 길이 $L$ 구간의 검사체적(Control Volume)에 대한 힘의 평형(정적 정상 유동 시 운동량 변화율 $\Sigma F = 0$)으로부터 마찰손실 압력강하식($\Delta P$)을 유도한다.

1) 검사체적의 힘의 평형 방정식:
유동 방향으로 작용하는 상류 정압력에 의한 힘 $F_1 = P_1 \cdot A$, 하류 정압력에 의한 반작용력 $F_2 = P_2 \cdot A$, 그리고 관 내벽면 전체에서 유체에 작용하는 전단응력(마찰저항력) $F_f = \tau_0 \cdot (\pi D L)$의 평형을 수립한다.

$$ \Sigma F_x = P_1 A - P_2 A - \tau_0 (\pi D L) = 0 $$$$ \Delta P \cdot \left(\frac{\pi D^2}{4}\right) = \tau_0 \cdot (\pi D L) \implies \Delta P = \frac{4 L}{D} \tau_0 $$

2) 관벽 전단응력($\tau_0$)과 마찰계수($f$)의 무차원화 정의:
유체역학에서 벽면 전단응력 $\tau_0$는 유체의 동압(Dynamic Pressure, $\frac{1}{2}\rho v^2$)에 비례하며, Darcy 마찰계수(Darcy Friction Factor, $f$)는 다음과 같이 정의된다.

$$ \tau_0 = \frac{1}{8} f \rho v^2 $$

3) 압력손실($\Delta P$) 및 수두손실($h_L$)의 도출:
전단응력 $\tau_0$를 힘의 평형식에 대입하면 Darcy-Weisbach 압력강하식이 완성된다.

$$ \Delta P = \frac{4 L}{D} \left(\frac{1}{8} f \rho v^2\right) = f \cdot \frac{L}{D} \cdot \frac{\rho v^2}{2} \quad \left[\text{N/m}^2\right] $$

유체의 비중량 $\gamma = \rho g$로 나누어 수두($\text{m}$) 단위로 표현하고, 유량 $Q = A v = \frac{\pi D^2}{4} v \implies v = \frac{4 Q}{\pi D^2}$를 대입하여 유량 기준 수식으로 전개한다.

$$ h_L = \frac{\Delta P}{\rho g} = f \cdot \frac{L}{D} \cdot \frac{v^2}{2g} \quad \left[\text{m}\right] $$$$ h_L = f \cdot \frac{L}{D} \cdot \frac{1}{2g} \left(\frac{4Q}{\pi D^2}\right)^2 = \frac{8 f L Q^2}{\pi^2 g D^5} \approx 0.0826 \cdot f \cdot \frac{L \cdot Q^2}{D^5} $$

$\implies$ 손실수두는 유속의 제곱($v^2$) 또는 유량의 제곱($Q^2$)에 비례하고, 관경의 5승($D^5$)에 반비례하는 유체역학적 엄밀식이다.

[메커니즘 2] 마찰계수($f$) 결정 이론: Hagen-Poiseuille, Moody 선도 및 Colebrook-White 식

Darcy-Weisbach 식의 핵심 파라미터인 마찰계수 $f$는 유동의 레이놀즈수($Re = \frac{\rho v D}{\mu}$)와 상대조도($\epsilon / D$)에 따라 결정된다.

1) 층류 영역 ($Re < 2,100$):
점성력 지배 유동으로 표면 거칠기($\epsilon$)와 무관하며, Hagen-Poiseuille 적분해로부터 엄밀 유도된다.

$$ f = \frac{64}{Re} $$

2) 난류 천이 및 활면/조면 영역 ($Re > 4,000$):
Colebrook과 White(1939)가 정립한 난류 마찰계수 음함수(Implicit Equation)로, Moody Chart의 수치적 근간이다.

$$ \frac{1}{\sqrt{f}} = -2.0 \log_{10} \left( \frac{\epsilon / D}{3.7} + \frac{2.51}{Re \sqrt{f}} \right) $$

3) 완전 난류(완전 조면 난류, Completely Turbulent Flow):
고레이놀즈수 영역에서는 점성 저층(Viscous Sublayer)이 관벽의 거칠기 돌기($\epsilon$)보다 얇아져 파괴되므로, $f$는 레이놀즈수와 무관해지고 오직 상대조도($\epsilon / D$)만의 함수가 된다(von Kármán-Prandtl 식).

$$ \frac{1}{\sqrt{f}} = 1.14 - 2.0 \log_{10} \left( \frac{\epsilon}{D} \right) = -2.0 \log_{10} \left( \frac{\epsilon / D}{3.7} \right) $$
[메커니즘 3] Hazen-Williams 공식의 수리학적 체계 및 실무 단위 환산

Hazen과 Williams(1902)는 상온($15\sim 20^\circ\text{C}$)의 물이 흐르는 관로의 방대한 실험 데이터를 통계 분석하여 Chezy 공식($v = C_{\text{Ch}} \sqrt{R_h S}$)을 변형한 지수형 경험식을 제안하였다.

1) 기본 평균유속식 (SI 기본 단위: $v\,[\text{m/s}], R_h\,[\text{m}], S\,[-]$):

$$ v = 0.8492 \cdot C \cdot R_h^{0.63} \cdot S^{0.54} $$

여기서 $C$는 Hazen-Williams 조도계수(무차원, 관 내벽 평활도), $R_h$는 동수반경($R_h = A / P_w = D / 4$), $S$는 동수경사(단위길이당 손실수두, $S = h_L / L$)이다.

2) 단위길이당 마찰손실수두($h_L / L$) 및 압력손실($\Delta P_m$)의 공학 단위식 유도:
유속식을 동수경사 $S$에 대해 정리하고 원관 조건($R_h = D/4$)과 $v = \frac{4 Q}{\pi D^2}$를 대입한다.

$$ S = \frac{h_L}{L} = \left( \frac{v}{0.8492 \cdot C \cdot (D/4)^{0.63}} \right)^{\frac{1}{0.54}} = \left( \frac{4 Q}{\pi D^2 \cdot 0.8492 \cdot C \cdot 4^{-0.63} D^{0.63}} \right)^{1.852} $$$$ \therefore h_L = 10.67 \cdot C^{-1.852} \cdot D^{-4.87} \cdot Q^{1.852} \cdot L \quad \left[\text{m}\right] $$

3) 소방설비 실무 배관 수리계산 표준식 (NFPA 13 및 국내 실무 단위: $P\,[\text{bar}], Q\,[\text{L/min}], d\,[\text{mm}], L\,[\text{m}]$):

$$ \Delta P_m = 6.05 \times 10^4 \cdot \frac{Q_m^{1.85}}{C^{1.85} \cdot d_{mm}^{4.87}} \cdot L_m \quad \left[\text{bar}\right] $$

$\text{MPa}$ 단위로 변환 시($1\,\text{bar} = 0.1\,\text{MPa}$): $\Delta P_{\text{MPa}} = 6.05 \times 10^3 \cdot \frac{Q_m^{1.85}}{C^{1.85} \cdot d_{mm}^{4.87}} \cdot L_m$

[메커니즘 4] Hardy Cross 반복법을 이용한 폐루프(Loop) 배관망 해석 알고리즘

스프링클러 격자망이나 상수도 소화배관망과 같은 다중 폐루프 배관망에서는 Kirchhoff의 전기회로 법칙과 동일한 2대 수리학적 지배법칙이 성립해야 한다.

  • 제1법칙 (절점 연속방정식, Node Law): 임의의 절점(Node)에서 유입 유량의 합은 유출 유량의 합과 같다. (질량보존: $\sum Q_{\text{node}} = 0$)
  • 제2법칙 (폐루프 수두손실 총합, Loop Law): 임의의 폐회로(Loop)를 시계방향으로 일주할 때 마찰손실수두의 대수합은 0이어야 한다. (에너지보존: $\sum_{\text{loop}} h_L = 0$)

유량 보정계수($\Delta Q$) 수식 유도:
임의의 관로 손실수두를 일반화된 지수식 $h_L = K Q |Q|^{n-1} = K Q^n$ (H-W: $n=1.85$, D-W: $n=2.0$)으로 놓는다. 가정 유량 $Q_0$에 미소 보정량 $\Delta Q$를 더한 실제 유량 $Q = Q_0 + \Delta Q$에 대해 폐루프 손실합을 1차 테일러 전개한다.

$$ \sum_{\text{loop}} h_L = \sum K (Q_0 + \Delta Q)^n \approx \sum K Q_0^n + \sum \left( n K Q_0^{n-1} \Delta Q \right) = 0 $$$$ \therefore \Delta Q = - \frac{\sum K Q_0^n}{n \sum K Q_0^{n-1}} = - \frac{\sum h_L}{n \sum \left| \frac{h_L}{Q_0} \right|} $$

$\implies$ Hazen-Williams 공식 적용 시 분모의 계수 $n = 1.85$, Darcy-Weisbach 공식 적용 시 $n = 2.0$을 대입하여 $|\Delta Q| < \epsilon_{\text{tol}}$(허용오차 $10^{-4}$)이 될 때까지 컴퓨터 반복 연산을 수행한다.

3 정량 비교 및 동역학 해석

비교 항목Darcy-Weisbach 공식Hazen-Williams 공식소방기술 실무적 의미 및 기술사 평가
공식의 본질차원 해석 및 검사체적 운동량 보존에 기반한 이론적 보편식상온 수유동 실험 데이터에 기반한 경험적 통계식D-W는 물리 법칙 엄밀성, H-W는 실무 계산 편의성 제공
손실수두 수식$$h_L = f \frac{L}{D} \frac{v^2}{2g} = 0.0826 f \frac{L Q^2}{D^5}$$$$h_L = 10.67 C^{-1.85} D^{-4.87} Q^{1.85} L$$유량 지수: D-W는 $Q^{2.0}$, H-W는 $Q^{1.852}$로 고유속에서 편차 발생
조도 파라미터절대조도 $\epsilon\,[\text{mm}]$ (관벽 거칠기 높이의 물리량) $\to$ 상대조도 $\epsilon/D$무차원 C계수 (C값이 클수록 매끄러움: 신관강관 120, CPVC 150)$\epsilon$은 물리적 조도이며, C계수는 점성·유속 영향이 뭉뚱그려진 계수
적용 가능 유체모든 뉴턴 유체 (물, 부동액, 포수용액, 가스, 오일 등 전 유체)오직 상온($15\sim 25^\circ\text{C}$)의 청수(Water)에만 한정부동액(글리콜) 혼합 수계설비에 H-W 사용 시 치명적 양정 부족 유발
유동 영역 한계층류($Re < 2100$), 천이구역, 완전난류 전 영역 적용천이 난류 및 완전 난류 일부에 한정 (층류·극저유속 적용 불가)미분무수(초난류) 또는 고점도 부동액(층류~천이)에서 H-W 파탄
화재안전기준 규정NFPA 750(미분무), 특수 유체 설비 수리계산에 의무/권장NFPA 13, NFPC 103(스프링클러) 표준 수리계산식으로 명시현행 실무는 H-W 지배적이나 PBD 고도화로 D-W 전환 가속
📈 [차트 1] 유량 Q(L/min) 증가에 따른 배관(100A, L=100m) 마찰손실수두(m) 비교 해석마우스 오버 시 수치 확인 가능
[그림 1] 유량-마찰손실수두 특성 상관 곡선: 배관 100A(내경 102mm, 길이 100m) 기준. Hazen-Williams 공식($Q^{1.85}$, 신관 $C=120$)과 Darcy-Weisbach 공식($Q^{2.0}$, 상대조도 기준)은 일반 설계유속 범위($2\sim 4\,\text{m/s}$, $1,000\sim 2,000\,\text{L/min}$)에서 상호 근접하나, 노후 배관($C=90$, 주황색) 시 스케일링으로 손실수두가 약 1.7배 폭증함. 또한 초고속 유동($>3,000\,\text{L/min}$)에서는 지수 차이로 인해 Darcy-Weisbach의 마찰손실이 H-W($C=120$)를 초과하여 이론적 분기가 발생함.
[도해 A] Darcy-Weisbach 원관 운동량 검사체적P₁·AP₂·A관벽 전단응력 τ₀ (둘레 πD × L)배관 길이 L (내경 D)힘의 평형: (P₁ - P₂)·A = τ₀·(π D L)τ₀ = ⅛ f ρ v² 대입 ⟹ h_L = f (L/D) (v²/2g)층류 f=64/Re | 난류 Colebrook-White 식[도해 B] Hardy Cross 폐루프 배관망 해석 모델Node ANode BNode CNode DQ_inQ_outQ₁ (h_L1)Q₂Q₄Q₃ (h_L3)Loop 1루프 수두손실 총합: Σ h_L = Σ K Qⁿ = 0보정유량 ΔQ = - Σ h_L / (n Σ |h_L / Q|)(H-W 적용 시 n=1.85, D-W 적용 시 n=2.0)
[그림 2] 시험장 3분 컷 손작도 핵심 도해: (좌) Darcy-Weisbach 원관 내 마찰 검사체적의 유체 정압력($P_1 A, P_2 A$)과 관벽 전단응력($\tau_0 \pi D L$) 평형 유도, (우) 다중 폐루프 소방배관망에서 절점 유량 보존($\Sigma Q = 0$) 및 루프 손실 총합($\Sigma h_L = 0$) 수렴을 위한 Hardy Cross 유량 보정($\Delta Q$) 알고리즘.

4 실무 실패 및 한계 메커니즘

[실패 1] 배관 노후화 및 미생물부식(MIC) 시 C계수 과대평가에 따른 방수압 미달

소방설비 수리계산 시 실무자들이 설계 기준서의 초기 신관 C값($C=120$)을 준공 후 수십 년간 불변하는 상수로 오판하는 오류가 빈번하다. 습식 스프링클러 배관 내부에서는 용존산소 부식, 슬러지 퇴적, 미생물 유도 부식(MIC: Microbiologically Influenced Corrosion)으로 인해 관벽에 결절(Tubercles)과 스케일이 증식하여 실제 유효 조도계수가 $C=80\sim 90$ 이하로 급락한다. 마찰손실수두는 $h_L \propto (1/C)^{1.85}$의 관계를 가지므로, C값이 120에서 90으로 저하되면 관마찰 손실이 무려 70% ($1.333^{1.85} \approx 1.70$) 폭증한다. 이로 인해 화재 시 최원격 헤드 기준 방수압력($0.1\,\text{MPa}$)이 $0.05\,\text{MPa}$ 이하로 강하되어 소화에 완전 실패한다.

[실패 2] 동결방지 부동액(Antifreeze) 및 포수용액 배관에 Hazen-Williams 맹목 적용

한랭지 옥외 주차장이나 냉동창고의 건식/프리아액션 대체용 부동액(프로필렌글리콜 또는 에틸렌글리콜 수용액) 계통이나 포소화설비(AFFF 폼 수용액)의 수리계산 시, 설계자가 관행적으로 Hazen-Williams 공식을 그대로 적용하는 치명적 실수가 발생한다. H-W 공식에는 점성계수($\mu$) 항이 전혀 없으므로 점성이 물의 $3\sim 8$배에 달하고 유동이 층류~천이 영역으로 전락하는 글리콜 혼합액의 점성 마찰손실을 반영하지 못한다. Darcy-Weisbach 식을 적용하여 레이놀즈수 저하에 따른 마찰계수 급증($f = 64/Re$)을 계산하지 않으면, 실제 요구 마찰손실이 설계치 대비 200~300% 이상 과소평가되어 소방펌프 양정 부족으로 인한 방수 불능 참사가 발생한다.

[실패 3] 미분무수(Water Mist) 초고압·초고속 영역에서 Hazen-Williams 공식의 지수 파탄

미분무수 소화설비는 배관 사용압력이 $10\sim 14\,\text{MPa}$에 달하고, 관내 유속이 $6\sim 10\,\text{m/s}$ 이상으로 치솟아 레이놀즈수가 $Re > 10^6$인 완전 조면 난류(Rough Turbulent) 영역에서 작동한다. Hazen-Williams 공식의 유량 지수 $1.85$는 중저유속 영역의 평활관 실험식에 기반하므로, 초고속 난류 영역에서 실제 유동 저항($Q^{2.0}$에 비례)을 추종하지 못하고 마찰손실을 과소평가한다. 따라서 NFPA 750 표준은 미분무수 배관의 수리계산 시 Hazen-Williams 공식 사용을 전면 금지하고 Darcy-Weisbach 공식의 적용을 강제하고 있다.

[실패 4] 밸브류 및 관부속류 등가길이($L_e$) 산정 시 C계수 미보정에 따른 누적 오차

소방배관 실무에서 엘보, 티, 게이트밸브 등의 국부마찰손실을 등가길이표(Equivalent Length Table)에서 찾아 단순 합산할 때, 설계 도면의 배관 조도계수($C$)와 등가길이표의 기준 조도계수($C=120$) 간의 보정을 누락하는 경우가 다반사이다. 등가길이 환산은 $L_e \propto (C_{\text{actual}} / 120)^{1.85}$의 보정계수를 곱해야 함에도 이를 간과하여 CPVC 배관($C=150$)에서 부속류 손실을 과대평가하거나, 노후 강관($C=100$)에서 부속류 손실을 심각하게 과소평가하는 수리학적 오류가 반복되고 있다.

5 기술사적 개선 대책 및 제언 (Conclusion)

소방기술사 종합 고찰 및 실무 엔지니어링 대책

소방 배관망 수리학은 단순 규정 수치의 대입을 넘어 유체의 열유체역학적 물성과 배관의 시계열적 노후도를 정밀하게 반영하는 고도화된 엔지니어링 프로세스이어야 한다. 이에 기술사적 관점에서 [유체별 공식 분리 $\to$ 장기 조도 C값 안전율 확보 $\to$ 루프망 수치해석 고도화 $\to$ 성능위주설계(PBD) 감리 검증]의 5단계 혁신 대책을 제언한다.

  1. 유체 물성 및 설비 특성에 따른 수리계산 공식 법적 분리 규정화:
    국가화재안전성능기준(NFPC)을 개정하여 일반 상온 청수 스프링클러 설비는 Hazen-Williams 공식($C=100\sim 120$)을 허용하되, 1) 부동액 혼합 계통, 2) 포소화설비 혼합액 배관, 3) 고압 미분무설비(NFPA 750), 4) 유체 온도 $40^\circ\text{C}$ 이상의 특수 공정 소화배관에 대해서는 Colebrook-White 기반의 Darcy-Weisbach 공식 사용을 법적 의무화해야 한다.
  2. 배관 내구연한 및 부식 환경을 고려한 설계 C계수 하향 표준화:
    현재 실무에서 관행적으로 적용하는 습식 강관 $C=120$, 건식 강관 $C=100$은 20년 이상 경과 시 방수압 결손을 방어하지 못한다. 따라서 설계 기준을 개정하여 습식 배관 설계 시 $C=100$을 표준 적용하도록 하고, 건식/준비작동식 배관에는 질소 봉입(Nitrogen Inerting System, 순도 98% 이상)을 의무화하여 산소와 수분을 차단함으로써 30년 이상 $C=120$의 조도를 영구 보존하도록 유도해야 한다.
  3. 격자형(Grid) 및 폐루프(Loop) 배관망 설계와 Hardy Cross 전산 연계:
    말단 분기배관(Tree Type) 방식은 단일 지점 폐쇄 시 광역 소화 실패를 초래하므로, 대형 물류창고 및 지하주차장에는 루프/그리드 배관망 설계를 적극 도입한다. 이와 함께 Hardy Cross 알고리즘 기반의 다중 루프 수치해석 프로그램을 활용하여 절점별 잔차 유량($\Delta Q < 10^{-4}$)을 수렴시키는 전산 검증 체계를 표준화한다.
  4. 국부마찰손실 저항계수(K-factor / Velocity Head Method) 도입:
    기존의 부정확한 등가길이($L_e$) 합산 방식에서 탈피하여, 밸브 및 이음쇠류의 마찰손실을 국부저항계수 $K_L$에 기반한 동압 손실($h_m = K_L \frac{v^2}{2g}$)로 직접 산정하는 방식을 도입함으로써 고유속·특수 배관망에서의 수리학적 계산 정밀도를 선진국(NFPA/FM Global) 수준으로 격상시킨다.
  5. BIM 기반 3D 소방수리계산(CFD/HydraCAD) 및 준공 수리시험 실측 검증:
    성능위주설계(PBD) 심의 대상 건축물에 대해 3차원 BIM 배관 모델과 수리해석 엔진(HydraCAD, PIPENET 등)을 직접 연동하여 설계 단계에서 압력 및 유속 분포를 가시화하고, 준공 감리 시 최원격 시험밸브함(Inspector's Test Connection)에서 실측된 동압·유량 데이터와 해석 프로그램 산출값의 오차율을 $\pm 5\%$ 이내로 상호 검증하는 디지털 엔지니어링 감리제를 정착시켜야 한다.
공프리(gongfree.kr) 무료 엔지니어링 & 학습 도구 모음

소방기술사 수험생 및 엔지니어를 위해 별도 설치 없이 브라우저에서 즉시 작동하는 100% 무료 스마트 웹툴: