1 개요 (Introduction & 패러다임 변화)
성능위주설계(Performance-Based Design, PBD)의 정착과 함께 화재안전성 평가의 핵심 도구로 미국 국립표준기술연구원(NIST)의 FDS(Fire Dynamics Simulator)가 표준 전산유체역학(CFD) 코드로 광범위하게 활용되고 있다. 과거 경험적 대공간 화재 수식이나 존 모델(Zone Model)의 공간적 균일성(Well-stirred assumption) 한계를 극복하고, 복잡한 3차원 공간 내 난류 연소 플룸, 연기 전파, 열전달 현상을 시간-공간적으로 정밀하게 포착할 수 있게 되었다.
그러나 FDS는 비압축성 저마하수(Low-Mach number) 근사 나비에-스톡스(Navier-Stokes) 방정식 기반의 대와류모사(Large Eddy Simulation, LES)를 채택하므로, 계산 격자 크기($\delta x$)가 해석 신뢰도와 계산 비용을 전적으로 지배한다. 화재특성직경($D^*$) 대비 격자 해상도($D^*/\delta x$)의 적정성 평가와 수치확산(Numerical Diffusion), 아격자(Sub-grid Scale, SGS) 연소모델의 한계를 엔지니어링적으로 검증하지 않은 시뮬레이션은 치명적인 과소/과대 설계를 유발하므로, 수치해석적 메커니즘과 타당성 검증 체계의 정립이 필수적이다.
2 핵심 메커니즘 및 지배방정식 유도
1. 저마하수(Low-Mach Number) 지배방정식 체계
화재 현상은 유속이 음속보다 극히 작은 저유속($Ma = u/c < 0.3$) 열유동이므로, 음파(Acoustic Wave)에 의한 고속 압력파동을 필터링하여 계산 시간 간격($\Delta t$)의 극단적 축소를 방지하는 저마하수 근사를 적용한다. 전체 압력($P$)은 열역학적 배경압력($p_0(z,t)$)과 유동을 구동하는 섭동압력($\tilde{p}(x,y,z,t)$)으로 분리된다($P = p_0 + \tilde{p}$).
(1) 질량 보존 방정식 (연속방정식):
$$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{u}) = 0$$여기서 $\rho$는 유체 밀도($kg/m^3$), $\mathbf{u} = (u,v,w)$는 속도 벡터($m/s$)이다.
(2) 운동량 보존 방정식 (Navier-Stokes):
$$\frac{\partial (\rho \mathbf{u})}{\partial t} + \nabla \cdot (\rho \mathbf{u} \mathbf{u}) = -\nabla \tilde{p} + \nabla \cdot \tau_{ij} + (\rho - \rho_\infty)\mathbf{g} + \mathbf{f}_b$$여기서 $\tau_{ij}$는 점성 응력 텐서, $(\rho - \rho_\infty)\mathbf{g}$는 부력항(Boussinesq 근사를 배제한 완전 밀도차 부력), $\mathbf{f}_b$는 스프링클러 액적 등 외력 항이다.
(3) 화학종 및 에너지 보존 방정식:
$$\frac{\partial (\rho Y_k)}{\partial t} + \nabla \cdot (\rho Y_k \mathbf{u}) = \nabla \cdot (\rho D_k \nabla Y_k) + \dot{m}'''_k$$$$\frac{\partial (\rho h_s)}{\partial t} + \nabla \cdot (\rho h_s \mathbf{u}) = \frac{dp_0}{dt} - \nabla \cdot \mathbf{q}'' + \dot{q}''' + \tau_{ij}:\nabla \mathbf{u}$$여기서 $Y_k$는 화학종 $k$의 질량분율, $h_s$는 현열(Sensible Enthalpy), $\mathbf{q}''$는 전도 및 복사 열유속 벡터, $\dot{q}'''$는 단위체적당 발열속도($kW/m^3$)이다.
2. 대와류모사(LES) 공간 필터링 및 격자해상도 지수 ($D^*/\delta x$)
FDS의 LES는 격자 크기($\delta x = (\Delta x \Delta y \Delta z)^{1/3}$)보다 큰 에너지 함유 와류(Large Eddies)는 직접 풀고, 격자 미만의 미세 와류는 Deardorff 또는 Smagorinsky 아격자(SGS) 모델로 와점성계수($\mu_{sgs} = \rho C_s \Delta^2 |\bar{S}|$)를 산정하여 소산 효과를 근사한다.
(1) 화재특성직경 (Characteristic Fire Diameter, $D^*$): 화원의 총 발열량($\dot{Q}$)과 대기 조건에 따른 열기둥(Plume)의 특성 치수를 무차원화한 기준 직경이다.
$$D^* = \left( \frac{\dot{Q}}{\rho_\infty c_p T_\infty \sqrt{g}} \right)^{2/5}$$여기서 $\dot{Q}$는 화원 열방출률($kW$), $\rho_\infty$는 주변 공기 밀도(약 $1.204\,kg/m^3$), $c_p$는 정압비열(약 $1.005\,kJ/kg\cdot K$), $T_\infty$는 절대온도(약 $293.15\,K$), $g$는 중력가속도($9.81\,m/s^2$)이다. 상온 공기 조건에서 $D^* \approx 0.0433 \times \dot{Q}^{2/5}$로 간략화된다.
(2) 무차원 격자 해상도 지수 ($R^* = D^*/\delta x$):
$$R^* = \frac{D^*}{\delta x}$$미국 원자력규제위원회(US NRC NUREG-1824) 및 NIST 표준에 따르면, 플룸 내 대류 열전달 및 연기 유동을 안정적으로 포착하기 위한 권장 격자 해상도는 $4 \le D^*/\delta x \le 16$ 범위를 만족해야 한다.
3 정량 비교 및 동역학 해석
| 구분 | DNS (직접수치모사) | LES (대와류모사 - FDS 기본) | RANS (시간평균모사) |
|---|---|---|---|
| 격자 스케일 | 콜모고로프 마이크로 스케일 ($\eta$) | 관성역 영역 격자 크기 ($\delta x$) | 기하형상 형상 기준 거시 격자 |
| 난류 해석 모델 | 전 영역 직접 해석 (모델 불필요) | 대형 와류 직접해석 + SGS 모델 | 전 와류 통계 모델링 ($k-\epsilon$, $k-\omega$) |
| 격자 민감도 ($D^*/\delta x$) | $D^*/\delta x > 100$ (초고격자) | $4 \le D^*/\delta x \le 16$ (NUREG-1824) | 격자 수렴 후 준평형 도출 |
| 계산 부하 | 극대 ($Re^3$ 비례, 실무 적용 불가) | 중등도 (병렬 클러스터 활용 실무화) | 낮음 (정상상태 엔지니어링 계산) |
| 비정상 와류 거동 | 완전 포착 | 플룸 주기적 퍼핑(Puffing) 완벽 재현 | 시간평균값으로 평활화(와류 상실) |
4 실무 실패 및 한계 메커니즘 (Engineering Failure Modes)
1. 수치 불안정성 및 CFL(Courant-Friedrichs-Lewy) 조건 위반
FDS의 양해법(Explicit scheme) 기반 시간 적분 시, 유체 이동 속도와 시간 간격의 비가 격자 크기를 초과하면 수치 발산(Numerical Blow-up)이 발생한다.
$$\text{CFL} = \frac{|\mathbf{u}| \Delta t}{\delta x} \le 1.0 \quad (\text{FDS 내부 가변 시간 제어: } 0.8 \sim 1.0)$$국소적으로 격자를 무리하게 세분화($\delta x \downarrow$)할 경우 $\Delta t$가 지나치게 작아져 전체 계산 시간이 비정상적으로 급증하거나, 반대로 강제 배연 디퓨저의 고속 토출부($|\mathbf{u}| \uparrow$)에서 CFL 조건 위반으로 연산 오류가 초래된다.
2. 아격자(Sub-grid Scale) 연소모델의 한계와 과대 환기 현상
실무에서 주로 사용되는 단순 연소 모델은 혼합분율(Mixture Fraction) 모델 또는 EDC(Eddy Dissipation Concept)이다. 이는 격자 내에서 연료 증기와 산소가 순간적으로 무한 반응속도로 반응한다고 가정(Infinitely fast chemistry: Mixed is burned)하므로, 화재실 구석이나 격자가 거친 영역에서 실제로는 소멸(Extinction)되거나 불완전 연소되어야 할 상황에서도 100% 완전 연소되는 것으로 왜곡되어 열방출률 및 독성가스 발생 농도가 부정확해진다.
3. 벽면 경계층 및 열전달 모델링 왜곡 (Wall Boundary Effect)
FDS는 고체 표면 근처의 난류 경계층을 완전 직접 해석할 수 없으므로 로그 벽함수(Log-law Wall Function)를 사용한다. 격자 크기가 $10\,cm$ 이상으로 조대화될 경우, 천장 제트 기류(Ceiling Jet)의 얇은 온도/속도 경계층(두께 수 $cm$)을 해상하지 못하여 천장면 전도 열손실이 과소평가되고, 스프링클러 헤드 감열체(RTI) 위치에서의 기류 속도가 왜곡되어 작동 시간이 늦어지는 설계 실패가 발생한다.
4. 멀티메시(Multi-Mesh) 경계 보간 오류
대규모 건축물 해석 시 다중 메시를 분할 적용할 때 메시 간 격자비(Mesh ratio)가 $1:2$를 초과하거나 정수배로 분할되지 않으면, 접합면(Mesh Boundary)에서 압력 및 속도 보간 불일치로 인한 인공적 와류(Spurious Eddy) 및 질량 누설(Mass Leakage)이 발생한다.
5 기술사적 개선 대책 및 제언 (Conclusion)
1. 격자 민감도 해석(Grid Sensitivity Analysis)의 법적 의무화
국내 소방청 성능위주설계 가이드라인에 단일 격자 수행 결과를 그대로 심의 도서로 제출하는 관행을 근절해야 한다. 동일 형상에 대해 최소 3단계 해상도($D^*/\delta x = 4, 8, 16$)의 격자 민감도 해석을 수행하여 리차드슨 외삽법(Richardson Extrapolation) 또는 격자 수렴 지수(GCI, Grid Convergence Index)를 산정하고 5% 이내의 수렴성을 입증하도록 제도화해야 한다.
2. 국소 영역 적응형 다중 격자(Nested Multi-Mesh) 설계 표준화
대공간 전체를 미세 격자로 해석하는 것은 하드웨어 한계상 불가능하므로, 플룸 발달 구역 및 감열체/연기감지기 직하부는 $D^*/\delta x = 10 \sim 16$의 미세 격자($5 \sim 10\,cm$)를 배치하고, 원거리 연기 하강 영역은 $D^*/\delta x = 4 \sim 6$의 조대 격자($20 \sim 30\,cm$)로 단계별 매칭($1:2$ 정수비 준수)하는 최적 도메인 분할 설계를 적용해야 한다.
3. V&V(Verification & Validation) 보고서 작성 및 PBD 심의 선진화
미국 NFPA 92(연기제어) 및 SFPE Handbook의 권고에 부합하도록, NUREG-1824 기준 실험 측정치와의 편차(Bias Factor)를 반영한 불확실성 마진(안전율 $1.2 \sim 1.5$)을 피난허용시간(ASET) 산정 시 필수적으로 반영하여 화재안전성 평가의 공학적 신뢰성을 확보해야 한다.