Simulating Smoke Filling in Big Halls by Computational Fluid Dynamics

Many tall halls of big space volume were built and, to be built in many construction projects in the Far East, particularly Mainland China, Hong Kong, and Taiwan. Smoke is identified to be the key hazard to handle. Consequently, smoke exhaust systems are specified in the fire code in those areas. An update on applying Computational Fluid Dynamics (CFD) in smoke exhaust design will be presented in this paper. Key points to note in CFD simulations on smoke filling due to a fire in a big hall will be discussed. Mathematical aspects concerning of discretization of partial differential equations and algorithms for solving the velocity-pressure linked equations are briefly outlined. Results predicted by CFD with different free boundary conditions are compared with those on room fire tests. Standards on grid size, relaxation factors, convergence criteria, and false diffusion should be set up for numerical experiments with CFD.


Introduction
Many fire models have been developed for studying preflashover fire and reported in the literature [1][2][3].Computational Fluid Dynamics (CFD) or known as field models [4,5], originally developed for fire simulations, are widely applied in smoke management system design.Development of fire field models started 30 years ago [6][7][8][9][10] at leading institutes all over the world.These are the Building Research Institute [11] and Tokyo University of Science in Japan; Notre Dame University [12] and the Fire Research Centre, National Bureau of Standards [13] in the USA; and the Fire Research Station [14] in the UK.Later, it was known as application of CFD [15] in mid 1990s.This was much widely applied when some of the research CFD softwares became commercial packages.Progress in CFD research for building fire applications was very slow in the past 10 years as explained in above.Most of the efforts were on applying CFD in performance-based design (PBD) on fire safety provisions.Fire Dynamics Simulator (FDS) [9,10] is now commonly used.
Applications are now restricted to only smoke exhaust systems including both static smoke extract system or natural venting, and dynamic smoke extract system with mechanical fans.Building fire is a low-speed chemically reacting flow [16] dominated by buoyancy [17].Three key fluid flow phenomena [4,5,16,17] on scalar transport are dominated by buoyancy-induced turbulence; heat transfer with thermal radiation at the burning object and the associated combustion chemistry involving hundreds of intermediate chemical reactions [18,19] should be considered at the same time.It is difficult to develop a simple CFD model to include all these three parts realistically without using any simpler models with empirical parameters.
However, there are always difficulties in finding resources for in-depth research as pointed out [20,21].On the other hand, many engineers with some CFD training were recruited as fire engineers for handling PBD [22,23] on building fire safety.Most of research efforts were therefore spent on applying CFD in building applications and generating beautiful graphics presentation, rather than working or getting more advanced CFD fires.
Therefore, CFD model was only used for designing smoke movement and design of air handling system such as smoke exhaust systems [24] or tunnel ventilation system.
Although the computing time required is still long, CFD has been improved with efficient numerical schemes.Faster personal computers with larger memory capacity are now available.Fluid flow problems up to 10 millions of computing cells can be accommodated with true 64-bit operating system.However, the computing time still takes several days.Efforts were even made on developing standards on using fire models [25,26], but how?
It was believed [3,15,27] in 1990s that CFD field model has the greatest potential to simulate building fires as discussed in several classical publications [6][7][8].A fire was taken as a heat source generating convective flow by buoyancy.A system of coupled, nonlinear, partial differential equations was set up on conservation of momentum, heat, and mass of smoke induced by fire.Temperature and velocity fields induced by the fire source were predicted.This approach is physical and incorporates basic principles supposed to predict more detailed information with minimum experimental data input.CFD model is good for predicting the thermal environment in big enclosures such as an atrium or a tunnel, though it is less suitable for regions near to a fire.The model should be developed by studying further on thermal radiation, combustion chemistry for common building materials, and turbulent mixing with air.Further, very few experimental verification and validation [28] were reported [28][29][30][31].Most of the experimental data used for justifying CFD results are originally for studying fire zone models [32].In fact, very limited systematic full-scale experimental works, except some by Ingasson and Olsson [33] on sprinkler fires and natural vents, were specially designed for justifying CFD models.The data reported by Steckler et al. [34] was performed originally for studying the doorway flow, but always used for verifying field models because of having many data at different positions.However, such limited studies did not give any promising results on CFD predictions as indicated by the VVT-FDS [35] study.Validation works were started as a longer-term project by the US Nuclear Regulatory Commission [36].Therefore, CFD are limited on studying smoke movement in regions far away from the fire, to avoid simulating combustion and thermal radiation effects.Even so, there are always challenges on accepting design based on CFD predictions.Consequently, field tests are required to evaluate whether the system will perform as the design.Atrium hot smoke test [37][38][39][40][41] is then required in many places of the Far East.There are also problems in training sufficient number of government officers, engineers and facilities managers to understand CFD.An updated review on applying CFD in designing smoke exhaust system is presented in the paper.

Smoke Filling
Fire is taken [6][7][8]10] as a physiochemical system of momentum, heat and mass transfer.A set of nonlinear and coupled partial-differential equations governing the conservation of momentum, heat, and mass are set up to describe the fire system.Whenever a fire occurs, natural convective flow fields would be generated, with the general form of the conservation law for a time-averaged fluid property.The thermal power from a fire in an enclosure is of the order of megawatt, such as that produced by burning a polyurethane foam sofa.The resultant flow is turbulent and appropriate for the use of turbulence models.
There are three main parts [4,5] in a CFD field model: turbulence in simulating the buoyancy induced motion, discretization of the set of partial differential equations, and algorithms in solving the velocity-pressure linked equations.
CFD predictions are mainly applied for designing smoke exhaust system in big space, particularly in atrium and tunnel.Static smoke extraction system is not so welcome by the Authority.Any change in the surrounding air pressure distribution due to wind or any other thermal perturbations will give very different system performance.Therefore, it is only suitable for small buildings with very low occupant loading, such as air cargo terminals.Performance of static smoke extraction system in crowded spaces such as shopping malls and public transport interchanges, particularly those located in deep underground, should be justified seriously by fire science and engineering.The scenario of pushing smoke down from air above must be demonstrated by full-scale burning tests that it would not happen.The design should not be accepted without such justification.
Other parameters including numerical parameters for stability and convergence such as relaxation factors, convergence criteria, residuals, and time steps, and parameters concerned in free boundaries, orifice coefficient of the openings, and how long the computing domain has to be extended to outside from an opening must be justified.Values of all these parameters should be selected very carefully with experiment.Different problems should be simulated by "trial and error" numerical experiments first.

Turbulence Model
Two common approaches [6,10] were used to simulate turbulent flow.These are the Reynolds Averaging of Navier-Stokes Equation (RANS) method and Large Eddy Simulation (LES).
The time-averaged values of the air flow variables including the momentum, density, pressure, enthalpy, and smoke concentration are solved using RANS.Any instantaneous value for these variables φ t is expressed as its time-averaged value φ plus the fluctuation φ .The set of equations describing conservation laws on φ t can be transformed in terms of φ with the product of the fluctuation terms φ for different flow variables separated out.Different turbulent models (e.g., [6]) have been proposed to close the set of equations to make it look similar to those for φ t , but the diffusion coefficient concerned becomes an effective value which is much greater than the laminar one.The k-ε model family of turbulent model is commonly used in simulating fire-induced flow fields.Time averages of the velocity vector components u, v, w (in a Cartesian coordinate system), the enthalpy h and the turbulence parameters k and ε are then taken as φ.There are many CFD packages based on RANS available in the literature in the 1990s.Commercial software was used to have a preprocessor menu for putting in building geometry.Because of the hardware limitation, CFD model based on RANS had been applied for designing smoke management systems in big halls and tunnels.However, there are criticisms on the accuracy due to grid size and boundary condition.The model cannot simulate combustion realistically and hence to predict carbon monoxide and soot concentrations in a fire.
In LES, turbulent motion is decomposed as large-and small-scale motions by filtering.The large-scale flow structures are calculated numerically by solving the differential equations concerned.Small-scale motions are modeled by subgrid-scale Reynolds stresses (e.g., [10]).The first step of LES is filtering.A flow variable φ( x, t) at position x and time t is decomposed into a large-scale component φ( x, t) and a small-scale (subgrid-scale) component φ ( x, t).The large-scale component φ( x, t) is a spatial average through a function G( x − x , Δ) in the domain of interest Ω through a filter width Δ as The choice of the filter function is important.LES approach would require much more computing cells than the RANS.There were discussions before that RANS with similar number of computing cells would give similar predictions by LES.LES is now commonly used in fire simulations.However, treatment of fire boundary should be watched.There are many physical empirical parameters in both approaches.Values of empirical parameters C 1 , C 2 , C 3 , C μ and C D used in RANS k-ε model [6] had already led to argument.These parameters are suggested to be tuned for different flow problems.But doing this would take CFD as a curve-fitting exercise [42].The parameter C s in the subgridscale model in using LES [43] is important.

Finite Volume Method
The set of partial differential equations describing the flow variable φ is commonly solved by the Spalding-Patankar control volume method [6,7].This is a finite difference method with the equations discretized by integration over a control volume.Special attention is paid to the convective part which depends on the flow direction.All these have been discussed clearly in the literature and will not be repeated in detail here.
Different schemes are available in discretizing the equations with the finite volume (or control volume) method.Common differencing schemes in fire simulations [6,20,44] are the upwind scheme, the hybrid scheme, the power law scheme, and the second-order Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme [44].The central differencing scheme was tested to have difficulties in getting converged results in fire simulations [15,45].
The upwind differencing scheme is equivalent to the zeroth-order interpolation with the choice of interfacial values depending on the sign of velocity.It is first-order accurate in terms of the Taylor series truncation error.Solution of the discretized equations based on this scheme is comparatively stable, seldom diverge or oscillate.Therefore, numerical schemes based on stable first-order upwinding are popular in studying convective transport.However, the firstorder accuracy makes it prone to numerical diffusion errors.False diffusion is a multidimensional phenomenon explained by taking the three-dimensional flow on each cell boundary as a one-dimensional flow.It occurs when the flow is oblique to the grid lines and when there is a nonzero gradient of the dependent variable in the direction normal to the flow.Such errors appear when the difference scheme fails to account for the true direction of the flow.Both the upwind and central difference schemes are not able to satisfy this requirement.
The false diffusion errors can be minimized by employing higher-order discretization.Scheme with second-order accuracy would involve a third-derivative truncation error.This becomes a source of unphysical oscillations under high convection conditions.Many schemes having higher order accuracy have been proposed to handle false diffusion.Those schemes involve discretization improvements in the simple upwind representation of the convection term.The QUICK scheme (e.g., [44]) is a popular one which uses a threepoint upstream-weighted quadratic interpolation for each interface.The resulting finite difference equation is thirdorder accurate based on Taylor series analysis.This is a thirdorder accurate upwind difference scheme with additional control volume cells involved in the interpolation which will reduce diffusion errors to a certain extent by bringing in a wider influence.The scheme possesses the stability of firstorder upwinding but is free of its second-order numerical diffusion.

Velocity-Pressure-Linked Equations
There is no explicit equation to calculate air pressure distribution in a CFD model.However, pressure gradients appear in the momentum conservation equation.Staggered grid system with the pressure (and other scalar quantities) nodes coincided with the cell faces of the control volume for velocity component is used for solving the equations.In this way, the pressure gradient can be represented properly for flows with oscillating pressure.Further, the velocity components are calculated at the exact locations for computing other scalar quantities.
The velocity components u, v, and w can be computed once the pressure field is known.But these computed velocity components might not satisfy the continuity equation.The Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) algorithm was proposed for solving this problem [6,44].This is an important step in developing CFD.Since there is no explicit equation for calculating the pressure field, pressure distribution pattern has to be guessed at the beginning.Therefore, SIMPLE required long computing time if the guessed pressure distribution is different from the reality.Modifications of the scheme were proposed for fire simulations.Three common algorithms [46,47] including the SIMPLE-Revised (SIMPLER), the noniterative Pressure Implicit with Splitting of Operators (PISOs); and the iterative PISOs were tested before on fire simulation for solving the velocity-pressure linked equations of the above problem [45].

Modelling and Simulation in Engineering
The standard k-ε turbulent model was employed for describing the turbulent effect.Numerical experiments were performed in simulating the fire-induced flow and temperature in the large test room reported by Ingason and Olsson [33].The predicted results using the three algorithms were compared with the experimental data.The results show that iterative-type solution method is good for fire field modeling.Solutions predicted from the original noniterative PISO algorithm might not satisfy the continuity equation; but give preliminary estimation for fire engineering application as computation time required is very much reduced.The treatment of free boundary conditions for the cases studied is important in achieving a reasonably steady-state solution.
For reducing computing time and memory demand, the solution domains were usually confined in the fire compartment.The boundary condition for pressure P along x-direction was set to However, reasonably steady-state results could not be obtained in this way.The flow field predicted used to be strong.In fact, the predicted mass drawn in and flow out from the compartment seems to be increasing indefinitely.One explanation for this phenomenon is that the fire source is located very close to the door opening where the free boundary conditions cannot be assumed.In order to overcome this problem, the solution domain must be extended to include the free space near the left door to a sufficient distance that the error due to the inaccuracy of boundary conditions was kept to minimum.In that study, the free boundary was extended to the outside air of the compartment.In this way, the neutral plane height can be predicted.But the computing time and memory requirement would be much longer and larger.The noniterative PISO scheme, though much less computation time is required for the simulation, is not reliable enough where accurate results are desirable.The iterative PISO method is commonly applied for fire simulations though SIMPLER [7] also gives good results.

The Fire Dynamics Simulator
After releasing version 3.01 [9,10] by the Building and Fire Research Laboratory, National Institute of Standards and Technology (NIST), USA in late 2002 [48,49] after almost 30 years of research, FDS is then actively upgraded.It is now a popular fire engineering tool.Compressible fire-induced air flow can be simulated with hot smoke taken as a thermally expandable gas.A set of governing equations suitable for simulating fluid flow induced by buoyancy with a low Mach number have been derived.The Boussinesq approximation is no longer necessary and constraints on inviscid fluid have been removed.The FDS model allows larger variations of both density and temperature.An important assumption is that the pressure variation is small, so that the ideal gas law can be approximated by decomposing the pressure into a "background" component, a hydrostatic component, and a flow-induced perturbation.
This assumption is reasonable for flows with low Mach number because the temperature is inversely proportional to density.The pressure in the state and energy equation is replaced by the background pressure to filter out sound waves of speed much faster than the typical flow speeds expected in a compartment fire.The momentum equations were rearranged based on the assumption that vorticity generation due to buoyancy prevails over its generation due to the baroclinic effect.Short turbulence length scale will be simulated by LES with the Smagorinsky subgrid models.
Turbulence, intermediate combustion chemistry and thermal radiation have to be included instantaneously to study the heat released from combustion.Combustion chemistry should be incorporated properly in FDS.There have been studies on intermediate chemistry for polymethyl methacrylate (PMMA), polyvinyl chloride (PVC), polyurethane foam (PUF), and wood.Possible reaction kinetics from the predicted mixing with oxygen and temperature were studied.The chemical reactions of combustion for PMMA are very complex with intermediate reactions classified under the seven groups on thermal decomposition; thermal oxidative decomposition; decomposition of monomer MMA; methane combustion; methanol combustion; formaldehyde oxidation and acetylene combustion.For the group of methane combustion [19], 77 intermediate reactions were identified.This reaction identification process had been reported for the other six groups.Such reactions will be reviewed carefully for identifying the important ones such as those generating heat.Only those key reactions are considered to be put into FDS.
The Smagorinsky constant C s in the subgrid-scale model used in FDS has to be evaluated [43].The boundary layer is not so important for this scale.The coefficient C s is varying from 0.1 to 0.25.Values of the Smagorinsky constant C s would affect the results related to the grid refinement predicting indoor airflow.For coarser grids, a larger value of C s would give better results.For fire simulations, C s was selected as 0.2.Another important coefficient, the subgridscale turbulent Prandtl number P rt , is varying from 0.2 to 0.9 as reported [50,51].Results for fire simulations are not very sensitive to changes in P rt , taken as 0.5 [51] for simulating the indoor airflow.
The mixture fraction model [48,49] is used in FDS to describe the fire burning process, and the model is based on the assumption that combustion is mixing-controlled.Each chemical specie of interest is described by a mixture fraction, which is a conserved quantity representing the fraction of its specie at a given position originating from the fuel.The mass fraction of each species is related to the mixture fraction by "state relations" from where necessary information can be extracted.For example, the state relation for the oxygen mass fraction would be used to calculate the local oxygen mass consumption rate.The local heat release rate can then be computed from the local oxygen consumption rate by assuming that it is directly proportional to the oxygen consumption rate and independent of the fuel involved.This procedure will be reviewed for its suitability in dealing with combustion chemistry.

Verification and Validation
CFD predicted results are always queried.A detailed analysis on quantifying the modeling and numerical uncertainties in those fire simulations is therefore required [28].Within the last decade, the quality or credibility of CFD fire simulations has drawn increasing attention from the society of fire engineering.The procedures for assessing the credibility of fire simulation codes are not agreed [52].
The term verification and validation (V&V) was used in operations research with definitions of V&V reviewed [28].Implementation of V&V for credibility quantification in CFD field model is studied in the preliminary stage.Some available CFD fire codes do not go through cautious tests of V&V, their work can only be regarded as just a demonstration that some special fire scenarios are probably simulated.This situation persists in the arena of CFD fire modeling, especially for the sector of CFD commercial codes [53].One of the possible reasons is that there is no commonly accepted V&V methodology available, though at least three different standards on the credibility of CFD have been published [28,51].
In contrast to zone models which have been well validated by experiments [15], to the best knowledge of the author, experimental validations of field models have not been carried out to the same extent as for zone models.In fact, very limited experimental work, except some by Ingasson and Olsson [33] on sprinkler fires and natural vents, have been performed for verifying field models.Experimental data reported by Steckler et al. [34] originally for studying doorway flow have been used many times for verifying field models [29].
There are validation documents on FDS and on other models released by the US Nuclear Regulatory Commission.Seven volumes of the draft report have just been made available [36].A series of 25 fire tests were conducted with primary purpose to provide data to validate computer models for nuclear power plant.The experiments were conducted in an enclosure of size 18 m by 12 m with height 6 m.The fires consisted of a simple gas burner, heptane or methanol liquid pools, or a PMMA solid fire.Four of the tests were conducted with a full-scale control room mockup in place.Parameters including the fire intensity, enclosure ventilation rate, and fire location were varied during the experiments.
Full-scale burning facilities should be developed for verifying and validating CFD results.Some on FDS had been reported [54].An exhaust hood fitted with a fanduct system is necessary to measure the heat release rate of burning combustibles by the oxygen consumption method.Different scenarios should be set up by burning combustibles of different sizes at different locations and ventilation conditions.Air temperature and velocity of the induced fire flow should be measured together with the heat release rate of the fire.

Functional Analysis
Functional analysis [55] was proposed to justify the CFD predicted results better.This is better than just saying the agreement between predictions and measured values is "good" or "bad."This method was proposed as a quantitative, statistical and analytical approach to evaluating fire zone models [32,56] first, and applied to verify CFD models [57].The transient measured (t i , m i ) and predicted (t i , p i ) curves are denoted by two vectors m and p by summing up all the n data points m i and p i at different time intervals t i , i = 1, . . ., n vectorally.The vector m − p will give the deviation.Two parameters were defined to give a quantitative measure on the similarity of the curve shape.
(i) The "norm" is the normalized relative distance between the vectors m and p.As recommended by Peacock et al. [55], the Euclidean "norm" for n data is calculated by m i and the corresponding p i at t i with s data points used to smoothen the fluctuating curves to provide better estimates inside as (ii) The parameter "cosine" (or "inner product") is defined as the angular difference between the resultant vectors as where A denotes (t i − t i−1 ).
The "norm" approaches zero and the "cosine" approaches 1 when the two curves are very close to each other.

An Example Case Study
Quality experiments on smoke filling were reported by Hägglund et al. [58], Li and Chow [59].A single-storey building of floor area 11 m by 39 m and height 8 m as shown in Figure 1 was considered.The wall and ceiling were constructed by mineral wool, and the floor was concrete.There were three vents labeled D1, D2, and D3, each of area 2.2 m by 0.6 m on the top of the front wall.The upper edge of the openings was positioned 0.2 m below the ceiling level.Two air inlets A1 and A2 of area 1 m by 2 m were at the lower edge of the left wall and front wall as shown in Figure 1.
The fires were positioned in the center of one half of the building.Methanol was used as the fuel and was burnt in square steel pans of sizes 0.56 m 2 and 2.0 m 2 .The heat release rate per unit horizontal area of the fuel was 440 kWm −2 .Hot air movement was observed by discharging smoke.The gas temperature was measured by thermocouples distributed over four vertical lines with 50 cm spacing as shown in   Two tests on smoke filling were selected to be simulated by CFD model FDS4.Smoke filling of two large room fire tests with and without wind outside the test room were studied.The numerical experiments are labeled as Cases 1 and 2 as follows.
Case 1. Fire size: 2.0 m 2  Internal air temperature: 15 The heat release rate per unit area of the fuel was 440 kWm −2 .The growth period of the mass burning rate was set to 60 s.The input heat release rate curves are shown in Figure 2. The maximum heat release rates were 880 kW and 246.4 kW for Cases 1 and 2, respectively.
Four thermocouple trees labeled T1 to T4 were set as in Figure 1(c).Predicted air temperatures for Case 1 are compared with experimental results in Figures 3(a  It can be seen that the two curves matched better with each other at height 4 m above floor, but not so good for 6 m above floor in the four thermocouple trees.At 6 m above floor, the predicted temperature increased more rapidly than experiment. Hot smoke layer height results at M1, M2, and M3 are shown in Figure 4(a).Three predicted curves were all kept steady at about 3 m above floor and fluctuated in a small scale.The smoke layer height measured in experiment was  As the heat release rate of Case 2 was smaller than that of Case 1, two positions at 6 m and 7 m above floor on each thermocouple tree were compared.The predicted air temperature curve agreed well with experiment.The CFD tool FDS appears to be good for simulating the fire scenario with wind effect.
Predicted smoke layer heights at M1 to M3 are shown in Figure 4(b).The smoke layer of Case 2 fluctuated more vigorously than that of Case 1.A possible explanation is due to wind outside.Smoke layer height at M1 fluctuated, particularly near to the left inlet.
Observed descending of the smoke layer in the experiment was reported [58] to be about 1 m.Smoke fell down to the top of the vertical inlet vent at height 1 m.There were three vents near to the ceiling.Hot smoke exhausted through two vents at the ceiling, and air was sucked into the room through the third one close to the fire as shown in Figure 6.Predicted velocity vector profile across vents is also shown in Figure 6.As shown in Figure 6(b), the air flow velocity vector at the two vents on the right was towards the inside of the room, while at the left vent it was opposite.This was the same as what reported in the tests [58].Predicted air temperature and smoke layer interface height were compared with the test results by applying functional analysis.Air temperature predictions in Figures 3 and 5     justified by taking n to be 8 and s as 2. Results on norm and cosine are shown also in the figures.For Case 1, FDS predictions gave norms from 0.08 to 0.37, and cosines from 0.71 to 0.99.In Case 2, the norms were from 0.07 to 0.14, and cosines from 0.82 to 0.93.CFD results predicted at some positions were satisfactory and very close to those measured values.For example, norm of T4 at 4 m above floor is 0.01, approaching to 0. Cosine is 0.95 with value close to 1.However, the predicted results were fairly acceptable at some positions with large deviations.For predictions at T1 with 6 m above the floor, norm is 0.37 and cosine is 0.71.Further, predicted air temperatures increased faster than the observed rate in the test.More vigorous justification of fire models is necessary before using them.Therefore, CFD fire models should be improved for application on fire protection design and evaluation [60].
A scaling system can be set up based on the values of norm and cosine on key flow variables to denote how good the CFD predicted results agree with experiment.However, this cannot be only set up by research workers.Authorities with responsibility to accept fire safety provisions for PBD [22,23] should play the key role in deciding whether to use a 3-point scale (bad, neutral, and good) or a 5-point scale (very bad, bad, fairly good, good and very good).

Free Boundary Conditions
In applying CFD to building fire hazard assessment, there are always windows and doors open to outside.Bi-directional flow would be observed experimentally, with hot gas flowing out and cool air coming into the room.There are many debates on specifying free boundary conditions for such openings.Flow parameters, particularly the pressure, have to be specified carefully.There are empirical correlations relating the pressure profiles across the vertical openings such as windows or doors, under different room geometries, heat release rates of the fires, and opening sizes.However, such boundary conditions might not give proper specification.A better approach is to extend the computing domains outside, as pointed out years ago by Schaelin et al. [60] in 1992 for simulating plume flow, and Mawhinney et al. [29] on fire studies in 1994.This part should be watched carefully as pointed out again before [21].The fire tests reported by Steckler et al. [34], Li and Chow [61] in a room of area 2.8 m by 2.8 m and height 2.18 m as shown in Figure 7 were taken as an example.A 0.3 m square methane burner of heat release rate 62.9 kW was placed at the centre of the test room.A door of width 0.9 m and length 1.85 m was located at the centre of the side wall.Two thermocouple trees, C and D, each with 22 probes, were placed at the corner and the doorway, respectively.Vertical air temperature profiles from floor to ceiling were measured at these two thermocouple trees.
Taking into account the accuracy and the computing time, cubic grids of size 0.05 m by 0.05 m by 0.05 m were used in CFD simulations with FDS.The room was divided into a grid system of 246,400 cells with 100 by 56 by 44 parts along the x-, y-, and z-directions, respectively, in the FDS input file.Three free boundary conditions labeled as FB1, FB2, and FB3 were tested in FDS simulations.
For free boundary condition FB1, FDS menu for free openings was used without extending the computing domain to outside.Air flow parameters at the open boundary depend on whether the flow is outgoing or incoming [34].Taking air pressure P as an example, P is related to horizontal velocity u by , outgoing, 0, incoming.
( The outgoing boundary condition assumes that the pressure perturbation is zero at an outgoing boundary and that P is constant along streamlines.The incoming boundary condition assumes that P is zero infinitely far away.At the boundary between two meshes, the pressure boundary condition is similar to that at an external open boundary.But P is taken from the adjacent mesh where the flow is incoming.
The lengths and heights of the extended space were once or twice of the room length and height, respectively, for FB2 and FB3.The input geometries of FDS simulations for these three boundary conditions are shown in Figure 8.
Relative air pressures with respect to ambient values at the doorway and horizontal components of velocity at 100 s, 150 s, and 200 s are shown in Figures 9 and 10.Air temperatures at different height at the corner in the room are shown in Figure 11.In applying functional analysis to compare with measured air temperature, values of norm for simulations FB1, FB2, and FB3 are 0.015, 0.013, and 0.013 respectively, and values of cosine are 0.59, 0.58, and 0.57 respectively.Predicted air temperatures are very close to experiments using the three boundary conditions.All have low values of cosine.However, the pressure profiles across the door were entirely different in view of Figure 9.It is difficult to accept the predicted pressure profiles of FB1, though experimental pressure data were not available on that set of tests.Predicted values are not so satisfactory as in the first example of simulating smoke filling in a bigger hall.A possible explanation is because combustion plays a key role in simulating fire in a small room.Further investigation would be worked on simulating the combustion process in CFD fire models.

Conclusions
Although CFD has been used in designing smoke exhaust systems in those atria [37][38][39][40][41] for over 30 years, not many studies were devoted to simulate turbulent mixing, intermediate combustion chemistry, and thermal radiation realistically without using empirical formula.Further, there were very few verification and validation works on the predicted results with limited full-scale burning tests, especially designed for CFD fire models.Fire engineers are then required to justify the CFD results while testing and commissioning the smoke exhaust systems [62].Hot smoke tests in the atrium site [37][38][39][40][41] are therefore required to evaluate the performance of atrium smoke exhaust systems in many places including Hong Kong.Fire models had been used in many projects but led to many arguments and debates on some submodels of those available tools.If the results are not demonstrated clearly, say by full-scale burning test, the process would appear as a "curve-fitting exercise.""Verification" exercises might be necessary, in addition to "engineering judgment."Although the associated costs for carrying out the burning tests might be high, it is necessary for some projects to demonstrate that adequate fire safety is provided.
Engineering tools with CFD models are applied in fire safety engineering.There are many reasons why the CFD field model is so attractive.The predicted "microscopic" picture of the thermal environment described by the velocity vector diagram, the temperature, pressure and smoke concentration (taken as mass of chemical specie) contours are useful for deriving the relevant macroscopic parameters for engineering purposes.The calculated vertical distribution of air temperature is useful in providing the appropriate thermal sensitivity for sprinkler heads.Also, the air entrainment rates for a fire plume were calculated from the predicted horizontal velocity components.There is still room for improving the accuracy of those engineering tools.Verification and validation of those models is particularly important.As pointed out in Hong Kong [62], both building designers and relevant government officials should be properly trained in using fire models.Note that the development of fire model itself is very rapid.For example, three-dimensional simulations with CFD models can be handled readily with a personal computer now.But 20 years ago, only two-dimensional simulations could be carried out in a mainframe computer.

Modelling and Simulation in Engineering
Relevant education and training must be provided [63,64].Offering degree programmes up to MSc level is necessary for upgrading the quality of fire engineering personnel.At least, continued professional development programmes are recommended to be offered regularly.
FDS [48,49] is very suitable in fire safety engineering.However, there is still much development to do and careless application of the technique would lead to wrong results.Verification and validation studies such as those works started by the US Nuclear Regulatory Commission [36] are important.Full-scale burning tests [61,62] for validating CFD fire field models are necessary.Filtering function h: E n t h a l p y k:

Nomenclature
Turbulence parameter m: A vector representing measured data m i : M e a s u r e d d a t a p o i n t s n: Index in functional analysis

P:
Air pressure P rt : Subgrid-scale turbulent Prandtl number p i : Predicted data point p: A vector representing predicted data t: T i m e u, v, w: Velocity vector components x: Position vector ε: Turbulence parameter Δ: Filter width φ: Time averaged value of flow variables φ t or φ( x, t): Transient flow variables φ or φ ( x, t): Fluctuation of flow variables Ω: Domain of interest.
Locations of thermocouples and measuring scales

Figure 1 (
Figure 1(c).Two video cameras were used to record the descending of smoke layer.Two tests on smoke filling were selected to be simulated by CFD model FDS4.Smoke filling of two large room fire tests with and without wind outside the test room were studied.The numerical experiments are labeled as Cases 1 and 2 as follows.

Figure 6 :
Figure 6: Predicted velocity vector across a plane with vents for Case 2.

Figure 11 :
Figure 11: Air temperatures at 200 s by thermocouple tree C.
g o r i n s k y c o e fficient in LES model C 1 , C 2 , C 3 , C μ , C D : Empirical parameters in RANS model G( x − x , Δ):