최근 건축물의 성능위주설계 및 화재안전성평가 시 화재시뮬레이션의 활용이 증가하고 있으며, 시간에 따른 3차원 공간의 화재현상의 예측이 가능한 대표적인 Field 모델인 Fire dynamics simulator (FDS)[1]가 가장 널리 활용되고 있다. FDS의 예측결과는 사용자에 의해 입력되는 물리 및 수치적 입력인자에 따라 상당한 차이를 가질 수 있기 때문에, 적용되는 입력인자의 정확성은 화재시뮬레이션 결과의 신뢰성에 핵심적이라 할 수 있다. 특히 열발생률(Heat release rate)은 유효연소열과의 관계를 통해 연료공급량으로도 해석될 수 있어, 가장 중요한 물리량으로 인식된다[2].

열발생률은 건축물의 화재안전성평가 시 화원정보로 입력되며, 최대 열발생률과 화재성장률(Fire growth rate)이 고 려된 화재성장곡선(Fire growth curve)의 형태로 표현된다. 추가로 화원의 크기 역시 단위 면적당 열발생률(Heat release rate per unit area, HRRPUA) 또는 연료의 공급유량과 직접적으로 연관되는 중요한 화원정보라 할 수 있다. 일반적으로 화원의 면적은 실규모 화재실험에 적용된 가연물의 표면적이 주로 인용된다. 그러나 실물 화재의 유효 연소면적은 가연물의 형상 및 배치, 화재확산 방향 및 속도에 따라 크게 변화될 수 있다. 그 결과 화재시뮬레이션에서 모델화원의 면적을 정의하는 것은 본질적 한계를 갖고 있다.

구체적으로 동일한 열발생률의 조건에서 화원면적의 변화는 화염길이의 변화에 의해 구획 내부의 화재특성에 상당한 차이를 야기할 수 있으며[3], 결과적으로 건축물의 화재위험성이 과대 또는 과소평가될 수 있다[4]. 이에 적절한 화원면적의 설정을 위한 다양한 연구들이 이루어지고 있다. 가장 기초적인 접근 방법으로서 다양한 화재조건에 대한 DB 확보를 위해 성능위주설계에서 주로 고려되는 공간에 대한 실규모 화재실험[5-7]이 수행되었다. 또한 가연물의 형상 및 배치에 따라 화재현상을 지배하는 화원면적의 도출[8], 건물용도에 따른 HRRPUA의 통계적 분석[9]이 이루어졌다. 수치해석적 연구로는 뉴질랜드의 건축법 C/VM2[10] 및 Fleischmann[11]의 연구를 통해 엔지니어가 활용할 수 있는 HRRPUA의 입력가이드가 제공되고 있다. 그럼에도 불구하고 화원면적의 선정은 매우 다양한 화재 조건에 따른 높은 불확실성으로 인하여 주로 사용자의 주관적 선택에 의존하고 있는 실정이다.

또한 화원면적의 변화는 수치해석의 수행 시에 격자 설정의 어려움을 동반 할 수 있다. 일반적으로 격자 설정은 특성화재직경(Characteristic fire diameter, D*)내에 삽입되는 격자의 수(D*/∆, 이하 N_grid가 4∼16개[12,13]일 때 적절한 예측결과를 제시한다는 경험적 제안에 근거하여 이루어지고 있다. 그러나 특성화재직경에는 화원의 물리적 크기에 대한 정보가 고려되지 않기 때문에, 특성화재직경 내 삽입되는 격자의 수는 화원면적 내 설정되는 격자의 수를 의미하지 않는다. 또한 특성화재직경과 화원의 물리적 크기가 상이할 때, 격자크기의 변화는 연소 반응대의 공간 해상도를 저하시킬 수 있다. 그 결과 성능위주설계에서 중요시되는 허용피난시간(Available safe egress time, ASET)과 연관된 주요 물리량들의 예측오차를 동반하게 된다[14]. Large eddy simulation (LES)의 기법이 적용된 FDS에서 격자의 의미는 난류 유동 및 연소모델의 의존도와 직접적으로 연관됨을 고려할 때, 화원의 면적변화에 따른 격자크기의 설정은 더욱 중요시 될 수 있다.

이러한 배경 하에 본 연구에서는 동일한 최대 열발생률을 갖는 과환기조건의 ISO 9705 표준화재실을 대상으로 수치해석적 연구가 수행되었다. 구체적으로 화원의 면적이 변화할 때 N_grid에 대한 온도, 화학종농도 및 열발생률 등 ASET과 관련된 물리량의 의존도가 검토되었다. 또한 격자크기의 영향으로부터 독립적인 조건에서 화원면적의 변화에 따른 구획 내부 및 구획 외부에서의 주요 물리량의 정량적 차이를 확인하였다. 본 연구결과는 성능위주설계를 위한 화재시뮬레이션의 예측 신뢰성을 향상시키는데 유용한 정보를 제공할 것으로 기대된다.

화원면적의 변화에 따른 특성화재직경 기반의 격자크기 설정 및 화재특성의 변화를 검토하기 위하여, 과환기 조건의 ISO 9705 표준 화재실이 수치해석 대상으로 선정되었다. Figure 1은 수치해석에 적용된 구획의 해석영역을 나타낸다. 경계조건이 수치해석에 미치는 영향을 최소화하기 위하여 구획 외부에 복도 형태의 공간을 배치하여 해석영역의 크기는 2.4m(x)×7.2m(y)×2.4m(z)로 설정되었다. 그림에 도시되지 않았으나, 구획 외부에서 수치해석 조건에 따른 물리량의 변화를 검토하기 위해 구획 외부의 바닥 및 천장, 좌·우면에 벽면조건이 부여되었다. 개구부 크기는 표준형상과 같이 0.8m(W)×2.0m(H)로 설정되었다. 화원면적 및 N_grid의 변화가 ASET 관련 물리량의 예측에 미치는 영향을 검토하기 위하여 구획 내·외부의 다양한 위치에서 물리량의 예측이 이루어졌다. 수치해석 조건에 따른 구획 내 열적환경의 검토를 위하여 구획 폭 방향의 중심선(x=1.2m)을 따라 y=0.6m, 1.8m, 3.0m의 위치에 열전대 트리가 3개 배치되었다. 각 열전대 트리를 통해 구획 내 0.2 m∼2.2 m의 높이에 대한 온도가 0.2 m 간격으로 측정되었다. 구획 내 바닥 및 천장의 앞부분(y=03.0m)과 뒷부분(y=0.6m)에서는 열유속의 예측이 이루어졌다. 수치해석 조건에 따른 ASET 관련 물리량의 변화를 검토하기 위하여 x=1.2m의 중심선을 따른 6개 위치에서 물리량의 예측이 이루어졌다. 개구부 두께(0.3m)의 중심부(y=3.75m)에서는 수치해석 조건에 따른 구획 내부의 물리량이 예측되었다. 구획 외부에서는 개구

부를 기준으로 0.6 m 떨어진 지점으로부터 3.0 m 떨어진 지점까지 0.6 m 간격을 갖는 5개 위치에서물리량이 예측되었다. 고려된 물리량은 O₂, CO₂, CO 및 온도, 가시도이며, 인명안전 검토 기준인 z=1.8m의 동일한 높이에서 예측되었다. 정사각형 화원면적은 세 개 조건이 고려되었으며, 한 변의 길이가 0.6, 1.2 및 1.8 m로 변화되었다.

Figure 2는 시간에 따른 화원의 공급 열량을 도시한 것이다. 화재지속시간은 건축물의 화재안전성평가에서 고려되는 600 s가 적용되었다. 본 연구의 목적과 연료종은 큰 연관성을 가지지 않으므로, 액체연료 헵탄(C₇ H₁₆)이 적용되었으며, 이를 고려하여 화재성장률은 Ultra-fast로 설정되었다. 최대 열발생률(1.04 MW)은 개구부를 통한 유입기와 화원에서 공급되는 연료 질량유량을 통해 얻어지는 총괄당량비(Global equivalence ratio)가 0.3이 되도록 설정되었다. 총괄당량비가 0.3인 과환기화재조건을 고려한 이유로 화재초기에 대한 안전성을 평가하는 성능위주설계의 목적과 환기부족화재에 대한 FDS 예측성능의 한계를 들 수 있다.

Table 1은 본 연구에서 고려된 수치해석 조건을 제시하고 있다. 앞서 언급된 것과 같이 최대 열발생률은 1.04 MW, 화재성장률은 Ultra-fast이며, 화원면적은 0.36, 1.44 및 3.24 m²의 세 가지 조건이 고려되었다. 동일한 최대 열발생률의 조건에서 화원면적의 감소는 연료공급유속을 증가시키게 되므로, 연료흐름의 모멘텀에 의해 난류부력화염의 거동을 상실할 수 있다. 이에 Froude number의 제곱근으로 식 (1)과 같이 정의되는 무차원 열발생률(Dimensionless heat release rate, Q*이 고려되었다. 일반적인 화재의 Q*는 2.5 미만이며, 이를 초과할 때 제트화염의 거동을 나타내는 것으로 알려져 있다[15]. 검토 결과 본 연구에서 고려된 화원면적에 따른 최대 Q*값은 2.48이며, 모든 조건에서 난류부력의 지배를 받는 거동을 나타낼 것으로 확인되었다. 참고로 Q_max, ρ_a, C_a, T_a, D및 g는 각각 화원의 최대 열발생률, 공기의 밀도 및 비열, 절대온도, 화원직경과 중력가속도를 의미한다.

격자크기에 따른 N_grid의 범위는 선행연구[12,13]에서 제시된 기준을 따라 선정되었다. 격자의 조밀도에 따른 4개(Coarse), 10개(Moderate), 16개(Fine)를 기준으로 6개, 12개의 격자가 삽입되는 조건을 추가하여 5개의 격자 크기가 고려되었다. 화원의 최대 열발생률에 따른 특성화재직경 D*는 식 (2)를 통해 도출되며, 여기서 Q_max, ρ_a, C_a, T_a, g는 각각 화원의 최대 열발생률, 공기의 밀도 및 비열, 절대온도, 그리고 중력가속도를 의미한다. 본 연구에서 적용된 화원의 최대 열발생률로부터 도출되는 D*는 0.97 m이며, 고려된 N_grid값에 따라 정육면체 형상을 갖는 격자의 한 변이 0.20 m∼0.06 m로 변화되었다.

각 조건은 N_grid값에 따라 N₄∼N₁₆으로 명명되었다. 이와 같은 격자계가 적용될 때 가장 적은 수의 격자는 5,184개(N₄)이며, 가장 많은 수의 격자는 192,000개(N₁₆)이다. 위와 같은 조건의 수치해석이 세 개의 화원면적에 대해 동일하게 고려되어, 총 15개 조건이 수행되었다. 이를 통해 동일한 최대 열발생률을 갖는 과환기조건의 구획화재를 대상으로 N_grid 및 화원면적의 변화가 ASET과 관련된 주요 물리량의 예측에 미치는 영향이 평가되었다.

FDS는 사용자가 입력하는 화원정보에 따라 화재현상을 예측하게 되므로 입력된 열발생률이 정확히 구현되어야만 N_grid 및 화원면적의 영향이 독립적으로 검토되었다 할 수 있다. 이에 수치해석 조건에 따른 열발생률의 예측결과가 비교되었으며, 그 결과가 Figure 3에 제시되었다. Figure 3(a)는 A_f = 1.44 m²에서 N_grid에 따른 열발생률의 예측결과와 입력된 화재성장곡선을 시간에 대한 함수형태로 비교하여 나타낸다. 그림을 살펴보면 시간에 따른 열발생률은 Kim[16]의 연구결과와 동일하게 N_grid의 변화와 관계없이 입력된 화 원정보에 대해 유사한 예측결과를 제시한다. Figure 3(b)는 열발생률에 대한 화원면적의 영향을 검토하기 위해 정상상태에서의 평균값을 비교한 결과이다. 정상상태는 Figure 3(a)에 음영으로 제시된 300∼600 s 구간이 고려되었다. 그림의 실선은 화원정보의 최대 열발생률(1.04 MW)을 나타낸다. 결과를 살펴보면 모든 조건이 N_grid 및 화원면적의 변화와 관계없이 입력값과 큰 차이를 나타내지 않음을 확인할 수 있다. 각 조건의 상대오차 중 최대값은 5.1%, 열발생률의 관점에서는 53 kW로 입력값에 비해 매우 작은 값을 가진다. 모든 조건의 평균 상대 오차는 2%로 N_grid 및 화원면적의 변화가 열발생률의 구현에 큰 영향을 미치지 않음을 알 수 있다. 결론적으로 다양한 물리량들의 비교를 통해 N_grid 및 화원면적의 영향을 독립적으로 검토하는 것의 타당성을 확보할 수 있다.

FDS를 통해 예측되는 가장 기본적인 물리량은 온도, 열유속 등 열적 물리량을 들 수 있다. 수치해석 조건의 변화가 열적 물리량에 미치는 영향을 검토하기 위해 구획 내(x=1.2, y=3.0)에서 예측된 온도(z=2.2m) 및 열유속(z=2.4m)의 평균값이 비교되었다. 평균값은 300∼600 s 의 정상상태 구간에 대해 도출되었다. Figure 4(a)는 온도의 비교를 나타내며, 동일한 N_grid에서 화원면적이 감소할수록 높은 값을 나타낸다. 구체적으로 화원면적에 따른 온도의 차이는 최대 216 ℃의 값을 나타낸다. 이러한 경향은 열유속에 대한 비교에서도 동일하게 확인되며, N₄ 조건에서 화원면적에 따른 예측값은 50 kW/m² 이상의 차이를 나타낸다. 이를 통해 동일한 최대 열발생률의 조건에서 화원면적의 변화는 내부 열적 물리량의 예측값에 상당한 차이를 야기함을 알 수 있다. 그러나 동일한 화원면적의 조건에서 N_grid에 따른 온도와 열유속의 비교를 살펴보면, N₄ 조건에서 상대적으로 높은 값을 나타내는 것을 제외하면 N₆∼N₁₆의 조건에서는 비교적 일정한 값을 제시하고 있다. 그러나 실무에서는 N₄의 기준 역시 활용되고 있으므로[17], 동일한 화원면적의 조건에서 N₄가 상대적으로 높은 열적 물리량의 분포를 나타내는 원인에 대한 검토가 요구된다.

Figure 5는 동일한 화원면적의 조건에서 N_grid가 변화할 때, N₄에서 높은 열적물리량을 나타내는 원인 분석을 위한 N_grid에 따른 정상상태 개구부 환기량의 비교를 나타낸다. 환기량은 300∼600 s 구간의 예측값에 대한 평균을 통해 도출되었다. 비교 결과, N₄ 조건은 모든 화원면적에서 N₆∼N₁₆에 비해 낮은 환기량의 예측값을 제시하고 있다. 정량적으로 N₄의 평균 환기량(0.94 kg/s)은 모든 N₆∼N₁₆ 조건의 평균값(1.03 kg/s)에 대해 약 9.4%의 상대오차를 가진다. 이는 격자크기 증가에 따른 개구부 유동에 대한 분해능 저하로 나타난 것으로 해석될 수 있다. 즉, 환기량의 과소예측에 의해 구획 내 연소환경이 상대적으로 과농해지며, 이로 인^해 내부 열적물리량이 증가하게 된다. 이는 N₄∼N₁₆의 범위 에서 적절한 화재현상의 예측이 가능하다는 연구결과[12,13]와는 상이한 경향을 보인다. 이러한 상반된 결과의 원인은 기존의 N_grid 범위가 개방된 공간에서의 실험에 대한 검증을 통해 제시되었기 때문으로 판단된다. Mun 등[18]은 환기부족화재에 대한 수치해석 시 격자의 크기는 개구부유동에 대한 해상도가 확보되어야 함을 보고한 바 있다. 그러나 Figure 4 및 5를 통해 과환기조건에서도 개구부 영역 내 격자 크기가 예측결과에 영향을 미치게 되며,N_grid 뿐만 아니라 개구부 유동의 해상도가 추가적으로 고려되어야만 예측 신뢰성이 확보될 것임을 알 수 있다. 추가적으로, 환기량이 적절히 예측될 때, 동일한 화원면적에 대한 N_grid의 변화는 예측결과에 큰 영향을 미치지 않음을 확인하였다.

Figure 6는 동일한 N_grid 조건에서 화원면적에 따른 구획 내 열적물리량의 변화 원인을 검토하기 위한 화원면적별 단면도의 비교를 나타낸다. 구체적으로 정상상태에서 화원면적에 따른 이론 화염면의 위치와 단면 온도분포 및 유동장을 나타낸다. 앞선 N₆∼N₁₆의 조건에서 N_grid의 영향이 작은 것으로 확인되었으므로, 상세한 분석을 위해 가장 조밀한 N₁₆의 조건이 비교되었다. 화염면 위치의 비교에서 화원면적의 감소에 따라 화염면이 개구부 방향으로 이동하는 것을 확인할 수 있다. 이는 연료공급의 흐름이 갖는 운동에너지에 의한 것으로 해석될 수 있다. 동일한 최대 열발생률, 즉 동일한 연료공급량을 갖는 화원의 면적 변화는 연료공급유속의 변화를 가져온다. 따라서 화원면적이 증가할 때 연료유동은 개구부로부터 유입되는 외기에 밀려 후방으로 이동하게 된다. 반면 화원면적이 좁아 연료흐름의 유속이 증가하는 경우 화염은 후방으로 이동하지 않는다. 이 경우 유입기는 화염에 의해 열적팽창을 거친 뒤 개구부로 이동하며, 이 유동에 의해 화염이 개구부 방향으로 이동하게 된다. 화원면적에 따른 유입기의 경로는 유선(Streamline)의 비교를 통해 확인될 수 있다. 이와 같이 화원면적의 변화는 내부유동 및 화염면 위치에 차이를 야기한다. 따라서 Figure 4에서 확인된 물리량의 편차는 측정위치와 화염면의 상대적 거리, 고온기체의 유동경로 변화에 의한 것으로 해석될 수 있다. 화원면적에 따른 내부유동의 차이를 살펴보면, A_f = 1.44 m², 3.24 m²에서 외기는 구획의 중심부 혹은 심부까지 유입되어 고온기체와 함께 배출된다. 그러나 A_f = 0.36 m²의 경우 외기는 구획심부로 유입되지 않으며, 화원 후방의 순환유동 생성을 확인할 수 있다. Hwang[19]에 따르면 화원의 후방에 생성되는 순환유동은 연소생성물의 구획 내 체류시간을 변화시키게 되며, 이로 인해 O₂, CO₂, CO 및 Soot 등 화학종의 농도가 변화할 수 있다. 이러한 물리량은 ASET의 평가 결과와 직결되므로 화원면적이 화학종 농도에 미치는 영향에 대한 검토가 요구된다.

Figure 7은 개구부 위치(x=1.2m, y=3.75m, z=1.8m)에서N_grid 및 화원면적에 따른 O₂, CO, 온도 및 가시도의 비교를 나타낸다. Figure 7(a)의 O₂ 체적분율(Volume fraction)의 비교를 살펴보면, N₄가 N₆∼N₁₆ 조건의 평균 예측값에 대해 약 9∼17%의 상대오차를 나타내는 것을 통해 N₄의 환기량에 대한 과소예측을 재확인할 수 있다. Figure 7(b)의 CO 농도(ppm) 비교를 통해 확인되는 N₄ 조건의 CO 농도 증가 역시 환기량의 과소예측에 의한 것으로 해석될 수 있다. 이를 통해 N₄를 제외하면 동일 화원면적에서 N₆∼N₁₆ 범위에 대한 N_grid의 변화는 예측값에 큰 영향을 미치지 않음을 알 수 있다. 그러나 동일한 N_grid의 조건에서 화원면적의 변화가 O₂ 및 CO에 미치는 영향은 상당하며, 이는 화원면적에 따른 내부유동 및 유입기의 거동을 통해 설명될 수 있다. 즉, 순환유동에 의해 연소생성물이 구획 내에 체류하거나, 외기가 구획 내 깊은 곳까지 유입되지 않는 경우 O₂ 농도의 감소 및 CO 농도의 증가가 나타날 수 있다. Figure 7(c)에서 화원면적의 감소에 따른 온도의 증가는 화염면 위치의 차이에 의해 설명될 수 있다. Figure 7(d)에는 가시도의 비교가 제시되었으며, N₄∼N₁₆의 조건에 대한 화원면적에 따 른 차이는 약 0.1 m로 그 영향이 크지 않은 것으로 확인된다. 이는 단위시간 당 연료공급량이 동일하다는 점, 연료공급량과 비례하여 사용자가 입력한 Soot yield를 구현하는 FDS의 특성을 통해 설명될 수 있다[1]. 그럼에도 불구하고 화원면적에 따른 O₂, CO 및 온도의 변화를 통해 환기량이 적절히 예측될 때, N_grid의 변화는 큰 영향을 미치지 않지만, 화원면적의 변화는 구획 내부의 열 및 화학적 물리량에 상당한 차이를 가져온다는 결론이 도출될 수 있다.

구획 내부에서 N_grid 및 화원면적의 변화가 미치는 영향이 검토되었으나, 건축물의 화재안전성평가는 한 개의 층 또는 건축물 전체에 대해 수행되는 경우가 일반적이다. 따라서 N_grid 및 화원면적의 변화가 화재실 외부에서의 물리량 예측에 미치는 영향이 검토되어야 한다. 이에 구획외부에서 개구부로부터 거리에 따른 ASET 관련 물리량의 정상상태 예측값이 비교되었다. 구체적으로 거리에 따른 화원면적 별 O₂, CO, 온도 및 가시도가 비교되었으며, Figure 8에 비교결과가 제시되었다. 화원면적의 변화가 미치는 영향을 평가하기 위해 화원면적에 따른 평균값 간의 표준편차가 각 거리에 대해 막대 그래프의 형태로 함께 제시되었다. Figure 8(a)는 O₂의 비교를 나타내며, 구획으로부터 멀어질수록 외기의 희석효과에 의해 O₂ 농도가 증가하는 것을 확인할 수 있다. 거리에 따른 화원면적의 영향을 살펴보면, 상대적으로 가까운 0.6 m∼1.2 m의 거리에서 높은 표준편차를 나타내고 있으며, 이는 구획 내부에서 확인된 화원면적의 영향이 외부에서도 나타난 것으로 해석될 수 있다. 그러나 개구부로부터의 거리가 증가해 1.8 m가 되는 위치에서는 표준편차가 큰 폭으로 감소하여, 화원면적의 변화가 큰 영향을 미치지 않음을 확인할 수 있다. Figure 8(b)에는 화원면적에 따른 CO 농도(ppm)의 비교가 제시되었다. 비교 결과를 살펴보면 CO 농도 역시 1.8 m의 거리에서 화원면적에 의한 표준편차가 크게 감소함을 확인할 수 있다. Figure 8(c)에는 온도에 대한 비교결과가 제시되었다. 온도 역시 O₂ 및 CO와 동일하게 1.8 m 거리에서 표준편차가 급격히 감소하여 화원면적의 영향이 감소하는 경향을 확인할 수 있다. 그러나 Figure 8(d)에 제시된 가시도는 구획에서 가장 가까운 0.6 m 거리에서부터 거리가 멀어질수록 증가하여 2.4 m 지점에서 다시 감소하는 다소 상이한 경향을 나타낸다. Han 등[20]의 연구에서 수행된 실제 성능위주설계 사례에 대한 분석결과를 살펴보면, 최근 5년간의 성능위주설계에서 가시도에 의해 ASET이 결정된 사례는 94%의 비 율을 차지한다. 따라서 화원면적의 변화가 ASET과 관련된 물리량의 예측결과에 미치는 영향이 상세히 분석되기 위해서는 가시도가 상이한 경향을 나타내는 원인에 관한 검토가 이루어져야 할 것이다.

이를 분석하기 위해 300∼600 s의 정상상태에 대한 평균 Soot 질량분율(Mass fraction, Y_soot)의 단면 분포도가 Figure 9(a)에 제시되었으며, 가시도의 측정위치가 거리에 따라 1∼5의 순서로 함께 제시되었다. 그림을 살펴보면 구획 내 상층부에 위치한 고농도의 Soot이 분포하며, 개구부 상단을 통해 배출되는 것을 확인할 수 있다. 이 때 고농도의 Soot을 포함한 고온기체는 1번 측정위치를 거쳐 외부의 천장면을 따라 흐르는 천장제트유동(Ceiling jet flow)의 형태를 갖는다. 2 및 3번 측정위치에서는 천장제트유동에 대한 외기의 유입으로 희석에 의해 Soot 농도가 감소하며, 이를 통해 해당 위치에서 가시도의 증가가 설명될 수 있다. 일반적으로 천장제트유동은 화원으로부터의 거리가 증가함에 따라 열 및 운동에너지의 감소를 겪으며 낙하하게 된다. 본 연구에서는 이러한 현상이 4 및 5번 측정위치에서 나타나며, 이를 통해 Figure 8(d)의 2.4 m∼3.0 m 위치에서 가시도의 감소 원인이 설명될 수 있다. 천장제트유동의 낙하는 Figure 9(b)에 제시된 온도분포에서도 동일하게 확인된다. 1 및 2의 측정위치는 상대적으로 고온영역 내에 위치하지만 3번 위치는 150 ℃ 이하의 온도분포를 갖는 영역에 위치한다. 그러나 4∼5번 측정위치에서는 천장제트유동의 낙하로 150 ℃의 온도분포 영역에 포함되는 것이 확인된다. 이는 Figure 8(c)의 1.8∼2.4 m 구간에서 나타나는 온도 증가에 대한 근거가 될 수 있다. 또한 천장제트유동이 연소생성물로 이루어짐을 고려할 때, Figure 8(a)∼(b)의 동일구간에서 나타나는 O₂의 감소, CO의 증가 역시 설명될 수 있다.

Figure 9를 통해 가시도가 상이한 경향을 나타내는 이유가 검토되었다. 이를 통해 복도가 주거공간보다 높은 일반적 건축물에서는 가시도 역시 화재실의 개구부로부터 1.8 m 이상의 거리가 확보될 때 화원면적의 영향을 크게 받지 않을 것이라 추측할 수 있다. 그러나 이러한 기준이 건축물의 화재안전성평가에 일반적으로 활용되기 위해서는 그 타당성이 검토되어야 한다. 이에 본 연구의 해석영역과 동일한 바닥면적을 갖는 해석영역에서 복도의 층고만이 3.0 m 로 증가할 때, 화원면적의 변화가 가시도에 미치는 영향이 추가로 검토되었다. 화원면적 외 변수의 영향을 최소화하기 위해 가장 조밀한 N₁₆ 조건에서 화원면적이 변화되었다. 한 화원면적에 따른 편차는 앞선 Figure 8에서 확인된 것보다 큰 값을 나타내고 있다. 이러한 현상은 화원면적에 따른 배출유속의 차이에 의한 것으로 해석될 수 있다. 앞서 언급된 것과 같이 동일한 최대 열발생률을 갖는 화원면적의 변화에 따른 연료공급속도의 변화는 개구부를 통한 배출유동의 속도에 차이를 가져온다. 이 때 천장의 높이가 낮다면 천장제트유동은 화원면적과 관계없이 물리량의 측정위치(z=1.8m)와 매우 가까운 거리에서 흐르게 된다. 그러나 복도의 층고가 증가하는 경우 넓은 화원면적 조건에서 배출유동은 느린 유속으로 인해 개구부로부터 가까운 위치에서 천장으로 상승한 뒤 천장면을 따라 흐르게 된다. 반면 좁은 화원면적으로 인한 빠른 유속의 배출유동은 상대적으로 먼 거리에 도달하며, 이로 인해 1.2 m 거리에서 A_f = 0.36 m² 조건의 가시도가 가장 낮게 나타난 것으로 해석될 수 있다.

이를 통해 개구부 근처(0.6∼1.2 m)에서는 화원면적의 변화가 가시도에 미치는 영향이 큰 원인이 설명될 수 있다. 그러나 화재실의 개구부로부터 1.8 m의 거리가 확보될 때, 화원면적에 따른 가시도의 편차는 크게 감소하며, 이는 개구부를 통해 배출된 유동이 부력에 의해 천장면에 도달하여 천장제트유동의 형태로 흐르기 때문이다. 화원면적의 영향은 개구부로부터의 거리가 증가함에 따라 감소하여 3.0 m의 거리에서는 화원면적에 따른 가시도의 편차가 0의 값을 갖는 것으로 확인되었다. 이를 통해 화재실로부터 일정 거리가 확보될 때 화원면적의 변화가 가시도에 미치는 영향은 감소함을 알 수 있다. 또한 FDS가 고온기체의 부력을 온도에 따른 주위류와의 밀도 차에 의해 도출하며[1], Figure 8(c)에서 화원면적이 온도 예측값에 큰 영향을 미치지 않았음을 결과를 고려할 때, 천장제트유동의 하강이 발생하는 거리에서도 화원면적의 영향은 매우 작을 것으로 판단된다. 결론적으로 화재실의 개구부로부터 1.8 m 이상의 거리가 확보될 때, 화원면적의 변화는 ASET과 관련된 화학종농도, 온도 및 가시도의 예측값에 큰 영향을 미치지 않음을 알 수 있다. 또한, 성능위주설계에서 ASET이 화재실 외부의 비상피난계단에서 평가됨을 고려할 때, 이러한 기준은 건축물의 화재안전성평가에 유용하게 활용될 수 있을 것으로 기대된다.

본 연구에서는 화재시뮬레이션 수행 시 화원면적 설정에 따른 불확실성이 수치적 예측결과에 미치는 영향이 평가되었다. 이를 위해 과환기 조건의 ISO 9705 표준 화재실을 대상으로 모델화원의 면적변화에 따른 화재특성 변화가 검토되었다. 또한 화원면적이 변화될 때 화재특성직격 기반의 격자 크기에 대한 의존도가 검토되었다. 주요 결과는 다음과 같다.

1) 과환기 구획화재에 대해 적절한 예측결과를 확보하기 위해서는 특성화재직경(D*) 내에 최소 6개의 격자가 포함되어야 하며, 이는 개구부 유동을 적절히 예측하기 위한 최소 공간해상도의 격자크기와 연관된다.

2) 동일한 최대 열발생률과 화재성장률이 적용될 때, 화원면적의 변화는 화염길이 및 위치 그리고 구획 내부의 유동구조의 변화로 인하여 고온 상층부의 열 및 화학적 특성에 상당한 변화를 가져온다.

3) 그러나 개구부로부터 충분한 거리(1.8 m 이상)가 확보된 구획 외부에서는 화원 면적이 변화되더라도 허용피난시간(ASET)과 연관된 주요 물리량들의 예측결과에 큰 차이가 없음을 확인하였다. 이러한 결과는 국내 성능위주설계에 적용되는 화재시뮬레이션의 화원면적 설정에 유용한 정보를 제공할 것으로 기대된다.