1 개요 (Introduction)
초고층 및 대규모 복합건축물 화재 시 연기와 유독가스의 거동은 열부력(Stack Effect), 외기 풍압(Wind Effect), 공조 및 급기가압 제연설비(HVAC/Smoke Control System)의 복합적 유체역학 상호작용에 의해 지배된다. 전산유체역학(CFD, 예: FDS)은 화재실 국소 플룸 해석에 정밀하나, 수십 층 규모의 건물 전체에 대한 장시간 비정상상태(Transient) 해석 시 격자수 폭증과 연산 부하로 실무 적용이 불가능하다.
이에 NIST(미국표준기술연구원)에서 개발한 CONTAM(CONTAMW) 기반의 망해석(Network / Multi-zone Model) 기법은 각 구획(Room, Corridor, Shaft)을 균일한 압력·온도를 갖는 '노드(Node)'로 모델링하고, 출입문·틈새(Crack)·창문·덕트를 비선형 오리피스 '유량 경로(Flow Path)'로 연결하여 베르누이 방정식과 질량보존 법칙으로 건물 전체의 차압 분포, 기류 이동량, 연돌효과 및 피난안전구역/부속실 차압 성능을 거시적·정량적으로 신속히 해석하는 성능위주소방설계(PBD)의 핵심 기법이다.
2 핵심 메커니즘 및 수식 유도
1. Multi-zone 모델의 기본 가정
각 구획(Zone $i$) 내 공기는 완전 혼합(Well-mixed) 상태로 가정하며, 구획 전체에 걸쳐 단일 압력 $P_i$와 절대온도 $T_i$를 가진다. 층고 $z$에 따른 수직 정수압 구배(Hydrostatic Pressure Gradient)는 밀도 $ ho_i$를 적용하여 선형 분포로 보정한다.
2. 유량 경로(Flow Path)의 오리피스 유량 방정식 유도
인접한 두 구획 $i$와 $j$ 사이의 개구부 및 간극을 통한 체적유량 $Q_{ij}$ 및 질량유량 $\dot{m}_{ij}$는 베르누이 정리와 멱법칙(Power Law)에 의해 다음과 같이 유도된다:
$$\Delta P_{ij} = P_i(z) - P_j(z)$$개구부를 통과하는 이론 유속 $v = \sqrt{ rac{2 |\Delta P_{ij}|}{ ho}}$로부터 오리피스 축류 및 마찰 손실을 고려한 유량계수 $C_d$를 적용하면 질량유량 $\dot{m}_{ij}$는 다음과 같다:
$$\dot{m}_{ij} = ho \cdot Q_{ij} = C_d A_{ij} \sqrt{2 ho_0 |\Delta P_{ij}|} \cdot ext{sgn}(\Delta P_{ij}) \quad ext{또는} \quad \dot{m}_{ij} = C_F (\Delta P_{ij})^n$$여기서 $A_{ij}$는 유효 누설면적($\mathrm{m^2}$), $ ho_0$는 유입 유체의 밀도($\mathrm{kg/m^3}$), $C_F$는 유량계수($\mathrm{kg/(s\cdot Pa^n)}$), $n$은 유동지수(Flow Exponent: 완전 난류 오리피스는 $0.5$, 층류 간극은 $1.0$, 실무 균열은 통상 $0.6 \sim 0.65$)이다.
3. 연돌효과(Stack Effect) 수직 압력 구배 방정식
수직 계단실/승강기로(Shaft, 절대온도 $T_s$)와 건물 외기(Outdoor, 절대온도 $T_o$) 사이의 임의 높이 $z$에서 발생하는 차압 $\Delta P_{stack}(z)$는 중력가속도 $g$와 이상기체 상태방정식($ ho = rac{P_{atm}}{R T}$)을 결합하여 도출된다:
$$\Delta P_{stack}(z) = (P_{s,0} - P_{o,0}) - g P_{atm} \left( rac{1}{R T_s} - rac{1}{R T_o} ight) z = \Delta P_0 - ho_0 g \left( rac{T_s - T_o}{T_s} ight) z$$중성대(Neutral Pressure Level, NPL) $z_{npl}$에서는 $\Delta P(z_{npl}) = 0$이 되며, 동절기($T_s > T_o$)에는 하부에서 외기 유입($\Delta P < 0$), 상부에서 샤프트 배출($\Delta P > 0$) 현상이 급격화된다.
4. 노드 질량보존 연속방정식 및 연립방정식 해법
정상상태(Steady-state)에서 임의 구획 $i$로 유입/유출되는 순 질량유량의 합과 공조/기계제연 급배기량($\dot{m}_{mech, i}$)의 합은 $0$이어야 한다:
$$\sum_{j=1}^{N_j} \dot{m}_{ji}(P_i, P_j) + \dot{m}_{mech, i} = 0$$건물 전체 $N$개 노드에 대해 비선형 대수방정식 군(System of Non-linear Algebraic Equations)이 형성되며, CONTAM은 이를 Newton-Raphson 기법과 스카이라인 역행렬 해석기(Skyline Matrix Solver)를 통해 노드 압력 벡터 $\mathbf{P} = [P_1, P_2, \dots, P_N]^T$를 잔여오차 허용치($\epsilon \le 10^{-5}$)까지 반복 수렴하여 해를 결정한다.
3 정량 비교 및 동역학 해석
1. 화재/연기 거동 수치해석 툴 3종 비교 분석
| 구분 | 망해석 모델 (CONTAMW) | 존 모델 (CFAST) | 전산유체역학 (FDS/CFD) |
|---|---|---|---|
| 지배 메커니즘 | 베르누이 + 1차원 유량경로 질량보존 | 2층대(Upper/Lower Zone) 질량·에너지 보존 | Navier-Stokes (LES/DNS) 난류해석 |
| 해석 공간 범위 | 건물 전체 (100층 이상 초고층 가능) | 수 개 인접 구획 (1~10개 실) | 단일 실 또는 국소 아트리움 한정 |
| 연산 소요 시간 | 수 초 ~ 수 분 (실시간 최적화 가능) | 수 분 ~ 수십 분 | 수십 시간 ~ 수일/수주 소요 |
| 주요 출력 결과 | 층별 차압분포, 도어 개방력, 계단실 누설량 | 연기층 높이, 상/하층 온도, 복사열 | 3차원 속도장, 와류, 가시거리, CO 농도 |
| 제연 실무 적용 | 부속실 가압설계, 연돌효과 대책, 댐퍼 밸런싱 | 피난개시시간 산정, 플래시오버 판정 | 제연경계벽, 배기구 유효흡입, ASET 검증 |
각 구획(외기, 거실, 부속실, 샤프트)을 개별 노드로 치환하고, 방화문·새시 간극·급기 댐퍼를 비선형 멱법칙 유량경로(Path)로 결합하여 매트릭스 수렴 연산 수행
4 실무 실패 및 한계 메커니즘
1. 누설면적(Effective Leakage Area) 가정 오류에 따른 차압 실패
- 현상: 실무 설계 시 ASHRAE 핸드북의 기본 간극 데이터(Tight, Average, Loose)를 무비판적으로 차용하여 시뮬레이션을 수행함으로써 준공 후 차압 시험(TAB) 시 기준 미달(40Pa 미만) 또는 과압(60Pa 초과) 발생.
- 메커니즘: 층간 슬래브 조인트, 승강기로 피트, 커튼월 스팬드럴 패널 틈새 등 비구획 누설경로를 모델링에서 누락할 경우 실제 계단실/부속실 가압 공기가 거실로 빠져나가지 못하고 갇히거나(과압 형성), 반대로 샤프트 외부로 과다 누설되어 차압 형성이 붕괴됨.
2. 1차원 Well-mixed 가정의 공간적 한계 (대공간 성층화 무시)
- 현상: 로비, 아트리움 등 수직 대공간 화재 시 연기층의 하강 및 확산 예측 실패.
- 메커니즘: CONTAM은 노드 내부의 온도와 오염물질이 순간적으로 완전 혼합된다고 가정하므로, 화재 플룸의 부력 상승, 천장 제트 기류(Ceiling Jet), 수직 온도 성층화(Thermal Stratification)로 인한 연기층 정체 및 재침강 현상을 모사할 수 없음.
3. 동적 피난 시 도어 개방(Door Opening Dynamics) 시계열 미반영
- 현상: 피난 개시 후 재실자의 동시다발적 방화문 개방 시 급기가압 계통의 압력 서지(Surge) 및 차압 급락 발생.
- 메커니즘: 정적(Steady-state) 모델링만 수행할 경우, 출입문 개방 순간의 일시적 차압 붕괴(0~10Pa)와 개방 후 닫힘 시 덕트 관성력에 의한 과압 피크를 제어하지 못해 피난 장애(도어 개방력 110N 초과)를 초래함.
4. HVAC 댐퍼 및 팬 P-Q 곡선의 비연계 (정풍량 가정의 오류)
- 현상: 덕트 계통의 마찰 압력 손실 및 송풍기 동특성을 이상적인 정풍량원(Constant Flow Source)으로 단순화하여 계통 불안정 유발.
- 메커니즘: 차압 댐퍼가 닫히거나 풍도 저항이 급증할 때 송풍기 작동점이 실속(Stall) 또는 서징 영역으로 이동하는 현상을 반영하지 못해 실제 풍량이 급감하는 실패가 빈번함.
5 기술사적 개선 대책 및 제언 (Conclusion)
1. CONTAM-CFD(FDS) 하이브리드 커플링(Coupling) 해석 기법 정립
건물 전체의 거시적 압력장 및 연돌효과는 CONTAMW로 해석하여 경계조건(Boundary Condition: 각 개구부 압력 및 질량유량)을 도출하고, 화재가 발생한 기준층 및 아트리움 대공간은 FDS(CFD)의 계산 영역으로 연계 입력하는 다중스케일 결합 해석(Multi-scale Hybrid Simulation)을 PBD 표준으로 제도화해야 한다.
2. 현장 기밀도 측정(Blower Door Test) 기반의 모델 보정(Calibration)
설계 초기에는 보수적 틈새 면적을 적용하되, 골조 및 창호 시공 단계에서 가압 감압법(ASTM E779 / ISO 9972)을 통한 층별 기밀도 실측값을 CONTAM Flow Path 계수에 피드백 입력하는 Dynamic Calibration 프로세스를 감리/TAB 기준에 의무 반영해야 한다.
3. 초고층 복합 연돌효과 저감을 위한 능동형 차압 제어 아키텍처 구축
동절기 연돌효과 극복을 위해 계단실을 15~20개 층 단위로 수직 구획(Vertical Compartment)하고, CONTAM 과도해석을 통해 도출된 각 존별 차압 변동성에 대응할 수 있도록 인버터(VFD) 차압 연동 급기팬 및 고속 구동 모터식 차압제어댐퍼(MVD)를 설치하여 도어 개폐 시 3초 이내 응답 특성을 확보해야 한다.
4. 소방시설 성능위주설계 평가 가이드라인의 정량적 평가 기준 혁신
현행 NFPC 501A의 40~60Pa 정적 기준을 넘어, CONTAM 시뮬레이션 시 외기 최저온도 조건(-15℃), 10분 이상 피난 시나리오에 따른 도어 개방률(최대 3개 층 동시 개방), 송풍기 P-Q 동특성 반영을 필수 심의 항목으로 강제 규정함으로써, 서류상 허수를 배제하고 화재 시 실질적 피난 안전성을 담보하는 엔지니어링 패러다임 전환이 시급하다.