Mathieu Duny

Dhionis Dhima

Jean-Pierre Garo

Hui-Ying Wang
email: wang@ensma.fr

Numerical and theoretical evaluations of a full-scale compartment fire with an externally venting flame

Keywords: compartment fire, externally venting flame, heat flux, façade, parallel sidewalls, opening 1. INTRODUCTION

Understanding of the physics and mechanisms of externally venting flames from enclosure fires with parallel sidewalls at the opening is fundamental to investigating fire spread of U-shape building façade. This work aims to identify the impact of separation distance of sidewalls on façade flame height, heat flux and sustained internal combustion. A series of numerical simulations is conducted in a full-scale compartment with 3.6 m in length/width and 3.1 m in height. An external façade wall is measured 6 m in height and 3.6 m in width. A pyrolysis surface with an area of 3.2x3.2 m 2 is placed in the middle of the enclosure with a theoretical heat lease rate varying from 1 to 8 MW. The semi-empirical correlations, derived from a cubic enclosure of 0.4 m with a HRR lower than 250 kW, are evaluated for a full scale façade fire for identifying the flame height dependence on the separation distance of sidewalls. It is found that an enhanced buoyancy-induced flow between sidewalls leads to an increase of the flame height by a factor of 50% and a rise of the peak of radiation heat flux by a factor of 20% as compared to a flat façade fire. The height of externally venting flame is sensitively reduced by a factor of 40-50% with an increase of the normalized sidewalls distance by opening width up to 3. The predicted flame height is more pronounced with an increase of flame height by a factor of 40% for a full scale façade than for a small scale one due to turbulence development. The predicted internal HRR is in quantitative agreement with the one from the correlation for a Global Equivalence Ratio (GER) below 3. The predicted temperatures are in good agreement with the experimental data except in the intermittent-like regions where an over-prediction of the temperature by a factor of 50% is found due to a reduction in accuracy of LES combustion model.

As an important fundamental topic in fire and combustion research, compartment fire phenomena have been studied for decades. The ventilated compartment fires have always occupied a large proportion of the research effort in fire science. The works [START_REF] Lu | Facade flame heights from enclosure fires with side wall at the opening[END_REF][START_REF] Lu | Experimental investigation on window ejected facade flame heights with different constraint side wall lengths and global correlation[END_REF][START_REF] Tang | Experimental study on flame height and temperature profile of buoyant window spill plume from an under-ventilated compartment fire[END_REF] verified theoretically and experimentally that the heat release rate depends on the amount of air consumed inside the compartment, which is proportional to the ventilation factor. Recent publications [START_REF] Merci | Numerical simulations of full-scale enclosure fires in a small compartment with natural roof ventilation[END_REF][START_REF] Wegrzynski | Experimental and numerical evaluation of the influence of the soot yield on the visibility in smoke in CFD analysis[END_REF][START_REF] Hostikka | Fire-induced pressure and smoke spreading in mechanically ventilated buildings with air-tight envelopes[END_REF][START_REF] Sikanen | Predicting the heat release rates of liquid pool fires in mechanically ventilated compartments[END_REF] using numerical simulations for the well ventilated enclosure fire are specially focused on studying qualitatively the thermal behaviour [START_REF] Merci | Numerical simulations of full-scale enclosure fires in a small compartment with natural roof ventilation[END_REF], visibility [START_REF] Wegrzynski | Experimental and numerical evaluation of the influence of the soot yield on the visibility in smoke in CFD analysis[END_REF] and pressure [START_REF] Hostikka | Fire-induced pressure and smoke spreading in mechanically ventilated buildings with air-tight envelopes[END_REF]. Sikanen et al. [START_REF] Sikanen | Predicting the heat release rates of liquid pool fires in mechanically ventilated compartments[END_REF] performed fully predictive simulations of liquid pool fires in mechanically ventilated compartments in a full-scale room. The quantitative agreement is fair with an overprediction of the burning rate of liquid fuel [START_REF] Sikanen | Predicting the heat release rates of liquid pool fires in mechanically ventilated compartments[END_REF].

When a fire inside a compartment develops to be at its most vigorous stage, flames may spill out of openings, forming externally venting flames, also known as façade fires. The structural barrier (for example, balconies) affects significantly the vertical spread of opening spill plume [START_REF] Xing | A computational study on structural barrier to vertical spread of window spill plume along building exterior facade[END_REF][START_REF] Cao | Numerical simulation of external fire spread through openings[END_REF][START_REF] Xin | Large eddy simulation on effects of wall pier and overhanging eave in preventing vertical fire spread[END_REF]. Numerical investigations show that for a compartment fire, opening geometry [START_REF] Duny | Numerical investigation on window ejected façade flames[END_REF] and prevailing ventilation conditions [START_REF] Pope | Quantitative comparison of FDS and parametric fire curves with post-flashover compartment fire test data[END_REF][START_REF] Martinez | Contribution to the study by numerical simulation of a fire plume coming out through a window on façade -Application to the comprehension of the wind influence on an eventually stack effect[END_REF] severely impact the interior fire behaviour via oxygen availability and the development of externally venting flame. In the work of Asimakopoulou et al. [START_REF] Asimakopoulou | Thermal characteristics of externally venting flames and their effect on the exposed façade surface[END_REF], an overall qualitative assessment of the correlations [START_REF] Beyler | Fire plumes and ceiling jets[END_REF][START_REF] Himoto | Modeling thermal behaviors of window flame ejected from a fire compartment[END_REF][START_REF] Lee | Heat flux distribution and flame shapes on the inert façade[END_REF] is presented as compared to measurements obtained in large-scale façade fire tests. It is found that the correlations [START_REF] Beyler | Fire plumes and ceiling jets[END_REF] originating from the experimental investigation of fire plumes significantly under-predict the external venting flame temperature and heat flux near the opening. The correlations [START_REF] Himoto | Modeling thermal behaviors of window flame ejected from a fire compartment[END_REF] exhibit a more conservative behaviour with a significant over-prediction of the temperature along a large-scale façade fire. Predictions of the heat flux to the façade, using various correlations [START_REF] Lee | Heat flux distribution and flame shapes on the inert façade[END_REF], highlighted the sensibility to the empirical extinction coefficient. Effects of façade inclination [START_REF] Tang | Window ejected flame height and heat flux along façade with air entrainment constraint by a sloping facing wall[END_REF][START_REF] Quinn | Experimental and numerical investigation of localised fire on glazing façades having different orientations[END_REF] on ventilated façade and their thermal behavior have been also a topic of research over the last years. Tang et al. [START_REF] Tang | Window ejected flame height and heat flux along façade with air entrainment constraint by a sloping facing wall[END_REF] studied the opening ejected façade flame characteristics with air entrainment constraint by a sloping wall. It is found that the façade flame height and the heat flux increased with an increase in sloping wall angle. Quinn et al. [START_REF] Quinn | Experimental and numerical investigation of localised fire on glazing façades having different orientations[END_REF] performed an experimental investigation of localized fire on façade with different inclination. It was observed that the inclined façade had a faster initial growth rate in terms of both heat flux and temperature for a similar heat release rate. The effects of external wind on the ejected flame along the façade wall from a compartment with multiple openings was also investigated [START_REF] Lu | Merging behavior of façade flames ejected from two windows of an under-ventilated compartment fire[END_REF][START_REF] Li | Wind effects on flame projection probability from a compartment with opposing openings[END_REF]. Lu et al. [START_REF] Lu | Merging behavior of façade flames ejected from two windows of an under-ventilated compartment fire[END_REF] described experimentally the merging behavior of flames ejected from two parallel openings of an under-ventilated compartment fire. Their results show that the flame height decreases as the separation distance between openings increases, and reaches the value corresponding to a single opening as the separation distance is beyond a critical value. The propagation of ejected externally venting flame is strongly accelerated due to the external wind from an opposing opening [START_REF] Li | Wind effects on flame projection probability from a compartment with opposing openings[END_REF].

Limit work [START_REF] Lu | Facade flame heights from enclosure fires with side wall at the opening[END_REF][START_REF] Lu | Experimental investigation on window ejected facade flame heights with different constraint side wall lengths and global correlation[END_REF] has been conducted to describe the effects of the presence of sidewalls adjacent to an opening on the ventilation factor, and consequently, the propagation of ejected externally venting flame. A real ejected fire propagation [START_REF] Asimakopoulou | Thermal characteristics of externally venting flames and their effect on the exposed façade surface[END_REF] from the openings of the vertical U-shape building façade to the upper floors is shown in Fig. 1. A buoyancy-induced flow between sidewalls significantly increase radiation and convection of the flames associated directly to the ejected spill flame height. The semi-empirical correlations [START_REF] Lu | Facade flame heights from enclosure fires with side wall at the opening[END_REF][START_REF] Lu | Experimental investigation on window ejected facade flame heights with different constraint side wall lengths and global correlation[END_REF] are commonly derived in conjunction with experimental data from a reduce-scale enclosure with sidewalls. All these experiments were performed using a heat release rate (HRR) lower than 250 kW. The correlations derived from large scale test are not available for the effects of sidewalls on the flame height and heat flux. In the current study, a real case of U-shape building façade fire (cf. Fig. 1) is simplified as a full-scale compartment fire with an externally venting flame between two vertical parallel walls, as schematized in Figure 2. Quantitative comparisons were first made between several computational results and well-instrumented experimental databases from a full-scale enclosure fire scenario, allowing to identify the strengths and weaknesses of the version 6.2 of FDS [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF]. A wide variety of possible configurations, including opening geometry and distance between opposite walls adjacent to an opening with various heat release rate has been numerically studied. This paper aims to broaden the scope of the validation of the available correlations [START_REF] Lu | Facade flame heights from enclosure fires with side wall at the opening[END_REF][START_REF] Lu | Experimental investigation on window ejected facade flame heights with different constraint side wall lengths and global correlation[END_REF] for the fire propagation behaviours of full-scale U-shape geometry through flame height and heat flux analysis with various separation distance of sidewalls.

NUMERICAL MODELLING

A picture of 3D view of the FDS computational domain and the coordinate system in the numerical simulation are shown in Figure 2. Advanced fire physics model, which is the subject of the current work, requires state-of-the-art submodels (combustion, multidimensional participating radiation, soot formation, heat transfer, condensed fuel vaporization, etc.) coupled with the flowfield governing momentum solution. The basis of the analysis is the conservation equations of mass, momentum, energy and species, a set of three-dimensional elliptic, timedependent Navier-Stokes equations. The finite-difference technique is used to discretize the mathematical representation of the reacting flow phenomena of interest here, and can be found elsewhere [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF].

In a Large Eddy Simulation, turbulence is modelled using a Deardorff's model via eddy viscosity as,

sgs v t k C     (1)
The model constant Cv is set to the literature value of 0.1 [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF]. The algebraic form of subgrid kinetic energy, ksgs, is obtained from a scale-similarity model as proposed in FDS [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF]. The filter width,  , is attached to the grid size. The turbulent Prandtl number is assumed to be constant with the value of 0.5 [START_REF] Zhang | Modeling of the combustion in compartment fires using Large Eddy Simulation Approach[END_REF].

The combustion processes are governed by the conservation equations for the mass fraction of the six major chemical species, such as fuel, O2, CO2, H2O, N2 and soot via a single global reaction step for wood. 𝐶 3.4 𝐻 6.2 𝑂 2.5 + 3.62(𝑂 2 + 3.76𝑁 2 ) → 3.32𝐶𝑂 2 + 3.10𝐻 2 𝑂 + 0.08𝐶 𝑠𝑜𝑜𝑡 + 13.61𝑁 2

(2)

The pine type of wood is chosen for the chemical reaction based on the available experimental database [START_REF] Kruppa | Comportement au feu des poteaux extérieurs en acier[END_REF]. The heat release rate is calculated from the local reaction rate of fuel via an Eddy Dissipation Concept (EDC) [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF].

c lim , 2 O 2 O eav 2 O i mix H ) Y Y ( H ) s Y , Y min( Q          (3)
Where c H  is the amount of energy released per unit mass of fuel consumed, and taken as 17.5 MJ/kg for pine type of wood [START_REF] Kruppa | Comportement au feu des poteaux extérieurs en acier[END_REF][START_REF] Desanghere | Détermination des conditions d'échauffement de structure extérieure à un bâtiment en situation d'incendie[END_REF]. The local key timescale,  mix , is supposed to relate approximately to the three processes such as diffusion, subgrid-scale advection and buoyant acceleration [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF]. The source term is multiplied by Heav(YO2-YO2,lim), where Heav is the Heaviside unit step function, which is zero when its argument is negative (YO2<YO2,lim) and 1 when it is positive (YO2>YO2,lim). This implies that flame extinction occurs as the local oxygen concentration is below a critical value, YO2,lim, evaluated as,

  o m f p lim , 2 O H T T C Y    (4)
Where Ho  (=13100 kJ/kg) is amount of energy released per unit mass of oxygen consumed, Tf adiabatic flame temperature, m T bulk temperature of a control volume and Cp specific heat. Note that all these parameters affect the history of externally venting flame propagation. In the current work, it is found that combustion at a stoichiometric fuel-air mixture allows to reproduce the experimental trend, signifying a positive argument of Heaviside function.

In FDS, soot production is derived from the fraction of the fuel mass that is converted into soot without taking into account the essential physical processes of soot formation and oxidation. With such soot yield approach, the predicted results are not improved using six radiation bands of CO2 and H2O due to an important production of soot even for a wellventilated fire. In the current simulations, it is assumed that soot is the most important combustion product controlling the thermal radiation from the flame and hot smoke. The gas behaves as a grey medium with one absorption coefficient without the spectral dependence. The grey gas assumption can lead to an over-prediction of the emitted radiation when soot production is small compared to the yields CO2 and water vapour. A radiative transfer equation is solved by using a discrete expression adapted to a finite volume method [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF].

A wall-damping is treated by using a wall function for convective motion near the wall due to the coarse meshes used here [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF]. The enclosure/façade walls with a thickness of 0.2 m are assumed to be thermally-thick. A one-dimensional heat conduction equation for the material temperature is solved using the thermo-physical properties in Table 1. The heat transfer equation is discretized by finite difference. The size of the cell size at the interior of the wall is selected automatically and is smaller than p c / k  where k denotes thermal conductivity,  density and cp specific heat. By default, the solid mesh cells increase towards the middle of the material layer and are smallest on the layer boundaries. The cells number is chosen automatically according to the thickness of the material layer [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF].

RESULTS AND DISCUSSIONS

In LES as used in FDS6 [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF], the filter width of sub-grid viscosity is attached to the grid size, and thus, the numerical model uses the quantities that are mesh dependent. In order to have a reliable prediction, it is important to use an adequate mesh for capturing the primary momentum transport of turbulent pool fire. For an enclosure fire, the viscous sub-layer near the walls is thin, this constraint translates into large computational grid requirement with a mesh size of mm that is generally incompatible with full-scale fire simulations [START_REF] Wegrzynski | Experimental and numerical evaluation of the influence of the soot yield on the visibility in smoke in CFD analysis[END_REF][START_REF] Hostikka | Fire-induced pressure and smoke spreading in mechanically ventilated buildings with air-tight envelopes[END_REF][START_REF] Sikanen | Predicting the heat release rates of liquid pool fires in mechanically ventilated compartments[END_REF]. Nevertheless, a priori estimates of the grid size using a characteristic length of pool fire D* are not trivial for starting the simulation [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF]. By treating the expendable fuel (wood) source as a pool fire, a characteristic length of pool fire D* is related to the heat release rate (HRR) as below [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF]:

𝐷 * = ( 𝑄 ρ∞ 𝑐 𝑝,∞ 𝑇 ∞ √𝑔 ) 2 5 ⁄ (5)
The HRR can be derived from the measured mass loss rate and the energy released per kilogram of fuel (wood). For this full-scale compartment fire, the peak in HRR is approximately up to 8 MW, giving a characteristic length, D* of 2.21 m. Generally, the large-scale energy-containing eddies that are controlled by the inviscid terms can be completely described when the characteristic length D* is spanned by roughly sixteen computational cells. The grid size depends considerably on what we are trying to accomplish. In fact, the compartment quickly fills with smoke, and the layered structure in a compartment fire can be resolved using significantly fewer grid points compared to the structure of the flame. A relatively coarse mesh size of 10 cm was estimated, and then gradually refine the mesh until we do not see appreciable difference in the predicted results. A cell size of 10 cm is adequate solely for evaluation of the smoke temperature field of a sizable enclosure fire.

The work of Zhang [START_REF] Zhang | Effects of computational domain on numerical simulation of building fires[END_REF] indicated that the domain extension is more efficacious in the wall fire normal x-direction, while is limited in the vertical z-direction and in the spanwise ydirection. It is found that extension of the computational domain beyond the vertical opening by one hydraulic diameter of the vent opening is sufficient for taking into account the interactions of internal and external flames. In the current compartment fire simulations, the overall computational domain is 6.5 m in the wall fire normal x-direction, 12 m in the spanwise y-direction and 12.5 m in the vertical z-direction, which is practically 3 -6 times of the hydraulic diameter of the vent opening of 2 m, large enough to assume a zero gradient condition at the free boundaries. An excessive domain extension needs the use of a highly compressed grid system, and build-up of numerical error could produce spurious results over the course of a LES calculation due to commutation of the filtering operation [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF]. The uniform grid system with a mesh size of 10 cm and a moderate computational domain are extensively used [START_REF] Merci | Numerical simulations of full-scale enclosure fires in a small compartment with natural roof ventilation[END_REF][START_REF] Hostikka | Fire-induced pressure and smoke spreading in mechanically ventilated buildings with air-tight envelopes[END_REF][START_REF] Zhang | Effects of computational domain on numerical simulation of building fires[END_REF], which offered the best tradeoff between accuracy and cost for large scale fire simulations.

Validation from a full-scale flat façade fire

A full-scale enclosure fire test conducted at CTICM [START_REF] Kruppa | Comportement au feu des poteaux extérieurs en acier[END_REF] in France, as schematized in Figure 3, has been largely served for numerical model validation [START_REF] Desanghere | Détermination des conditions d'échauffement de structure extérieure à un bâtiment en situation d'incendie[END_REF]. The internal compartment dimensions were 3.6 m in length (x), 3.6 m in width (y) and 3.1 m in height (z) equipped with an external flat façade wall measured 6 m in height (z) and 3.6 m in width (y). The fire compartment opening measured 1.4 m in height (h) and 2.6 m in width (w), corresponding to a window, is located at the façade wall.

The compartment and façade lining walls are made of the resistant gypsum plasterboards with a thickness of 0.2 m. The physical properties, as presented in Table 1, of the experimental setup, such as density, thermal conductivity and specific heat in ambient conditions are used in the numerical simulation. Aiming to simulate realistic building conditions, an expendable solid fuel (wood) is employed. Wooden cribs are typically used as the fuel load on compartment fires, as they have been shown to be representative of typical room furnishings and have good repeatability [START_REF] Kruppa | Comportement au feu des poteaux extérieurs en acier[END_REF][START_REF] Desanghere | Détermination des conditions d'échauffement de structure extérieure à un bâtiment en situation d'incendie[END_REF]. There are six vertical profiles (A-F), and each of the profiles comprises 6 thermocouples positioned at regular intervals over the height of the façade wall (cf. Fig. 3a,b). Predicting mass loss rate of wood cribs with the potential of burning away via heat and mass transfers is complex, and generally incompatible with full-scale fire simulations due to large computational grid requirement [START_REF] Desanghere | Détermination des conditions d'échauffement de structure extérieure à un bâtiment en situation d'incendie[END_REF]. In the current numerical simulations, the wood cribs are in a simplified configuration in which the mass loss rate collected experimentally (cf. Fig. 4) is prescribed with a pyrolysis surface area of 3.2x3.2 m 2 which is placed in the middle of the enclosure, slightly elevated at a height of 0.3 m. The pine type of wood is chosen as in the experiments. FDS [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF] manages to correctly predict the evolutions of the total HRR and the internal one inside the compartment, as presented in Figure 5. The heat release rate is calculated from the combustion model based on an Eddy Dissipation Concept (Eq.3). Integration of the local heat release rate over all cells with a volume V  inside the enclosure allows to monitor the internal heat release rate.

Evolution of the HRR (cf. Fig. 5) can be said to follow the four characteristic stages : 1) The initial growth, pre-flashover phase up to 300 s, during which combustion is limited at the interior of the compartment with an identical value between the total and internal HRRs. 2) The fully developed, flashover phase, defined by the rapid increase in HRR to 8 MW, during which the flame front moves gradually from the vicinity of the fuel source due to lack of oxygen, expanding radially and horizontally towards the opening. This signifies the beginning of the external flame jets and quick flashes appearing at the exterior of the fire compartment. 3) The post-flashover phase during a long period of about 20 min, the sustained external combustion of unburnt volatiles consistently covers the region above the opening; a large part (4 MW) of about 50% of the total energy is released outside along the façade. 4) The decay phase of external flaming during which flame jets appear intermittently outside the compartment, characterized by a burning transition back to fuel-control before auto-extinction is achieved to burnout of the fuel load. Histories of the gas temperature along the façade wall are analysed. In fire environment, the temperature experimental data from thermocouple bead are radiation corrected. In order to compare the simulated and measured temperature data, the thermocouple model is incorporated into FDS6 [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF]. The thermocouple model (𝑇 𝑇𝐶 ) is based on a simple balance for thermocouple bead and is informed by LES-based estimates of the rates of convective and radiative heat transfer, as follows :

𝜌 𝑇𝐶 𝑐 𝑇𝐶 𝑑𝑇 𝑇𝐶 𝑑𝑡 = 𝜀 𝑇𝐶 (𝑈 4 ⁄ -𝜎𝑇 𝑇𝐶 4 ) + ℎ(𝑇 𝑔 -𝑇 𝑇𝐶 ) = 0 (6)
Where Tg denotes the gas temperature derived from a state equation, as RTM / P Tg 

, 𝜀 𝑇𝐶 the emissivity, U the integrated radiative intensity and h the heat transfer coefficient. Since the type of thermocouple used in the tests [START_REF] Kruppa | Comportement au feu des poteaux extérieurs en acier[END_REF] is unknown, the default values proposed in FDS [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF] for nickel are taken into account in the simulation.

The two methods based on both the gas temperature and thermocouple model are compared with the experimental data in Figure 6. The radiation heat flux to the thermocouples bead depends mainly on three factors such as flame size, concentrations of gaseous and particulate soot emitting species, and view factor from flame to the exposed thermocouples bead. The gas temperature near the wall at position P1 or P6 (cf. Fig. 3b) is mainly controlled by convection, and contribution by radiation appears less important due to a decrease in view factor. The temperature form thermocouple model near the façade wall (position P) is slightly higher than the gas temperature. The positions Q1 and R1 are close to the opening with a normal distance beyond 0.5 m from the façade wall. Here the gas temperatures, Tg, are much lower as compared to thermocouple temperatures, TTC. This is mainly attributed to the important contribution of radiation flux received by the thermocouples bead with an increase of view factor. However, any attempt to draw a general conclusion for some metrics is discouraged due to the uncontrollable factors such as fire size, chemical species, etc.

Generally, the thermocouple model (𝑇 𝑇𝐶 ) shows a remarkable agreement with measured values for a consistent and continuous behavior of externally venting flame. In fact, an intermittent zone of externally venting fire which coincides with forward pulsations of flame with a fluctuating, turbulent character is experimentally identified [START_REF] Kolb | Cross flow effects on the flame height of an intermediate scale diffusion flame[END_REF][START_REF] Tang | An experimental study on the intermittent extension of flames in wind-driven fires[END_REF]. The EDC type [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF] is temperature independent combustion model, which is incapable of reproducing the frequencies of flame presence/absence in the intermittent zone at the top of fire plume. An overestimate of the temperature close to the opening ejected flame position Q6 is attributed to a reduced accuracy of LES combustion model (EDC) under intermittent-like conditions. A rigorous numerical study [START_REF] Ning Ren | Large eddy simulation of vertical turbulent wall fires, The 9th Asia-Oceania Symposium on Fire Science and Technology[END_REF] using a similar combustion model on a vertical wall fire with a mesh size of 3 mm is incapable of identifying the intermittency values of the oscillation behavior. Note that both the local gas temperature [START_REF] Asimakopoulou | Thermal characteristics of externally venting flames and their effect on the exposed façade surface[END_REF][START_REF] Wang | Large-eddy-simulation of buoyancy-driven fire propagation behind a pyrolysis zone along a vertical wall[END_REF] and thermocouple model [START_REF] Ning Ren | Large eddy simulation of vertical turbulent wall fires, The 9th Asia-Oceania Symposium on Fire Science and Technology[END_REF] are largely used for evaluation of limitations of fire models. Quantifying the uncertainty and accuracy of both the two methods in position Q6 ensures that overestimate of the temperature is attributed to the predictive capabilities of the fire models. 6a) at position P1 and P6 with a discrepancy within 5% 6b) at position Q1 and Q6 with a discrepancy within 50% 6c) at position R1 and R6 with a discrepancy up to 300% The deviations of thermocouple model from gas temperature is an alternative to the classical experimental approach in which raw thermocouple data are corrected by radiation heat transfer effects. The thermocouple model is used to check the grid refinement from 20, 10 to 5 cm on the predicted temperatures for the positions P, Q and R (cf. Fig. 3b) in Figures 7, 8 and9. The measured temperatures which are averaged from the points A to F at each given elevation are also provided in these figures. The computation work was performed with a single mesh on a personal computer with a 3.2 GHz Intel Xeon processor. With a mesh size of 10 cm, the CPU time for a single simulation with a physical time of 1800 s and a time-step of about 10 -3 s is approximately 50 h. With halving of the grid size to 5 cm, the numerical simulation was achieved only at 1100 s because the time required for the simulation increases by a factor of 2 4 (a factor of two for each spatial coordinate, plus time) [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF]. The fine mesh size of 5 cm is used just for judging the grid sensibility, and computations of the large-scale fire are limited to relatively coarse meshes of 10 cm everywhere [START_REF] Wegrzynski | Experimental and numerical evaluation of the influence of the soot yield on the visibility in smoke in CFD analysis[END_REF][START_REF] Hostikka | Fire-induced pressure and smoke spreading in mechanically ventilated buildings with air-tight envelopes[END_REF][START_REF] Sikanen | Predicting the heat release rates of liquid pool fires in mechanically ventilated compartments[END_REF][START_REF] Zhang | Effects of computational domain on numerical simulation of building fires[END_REF].

It can be noticed that the numerical results with a large mesh size of 20 cm do not match the measured ones near the wall region (P) particularly far away from the opening. As compared to the experimental data, change of the predicted temperature plots along the façade by using the medium (10 cm) and fine (5 cm) mesh sizes is less than 5%. This confirms that a mesh size of 10 cm gives practically a grid independent resolution, and thus is considered as appropriate for the purposes of this work. The predictions are in remarkable agreement with experimental results for positions P near the façade wall (cf. Fig. 7) and R in the plume region (cf. Fig. 9), implying that the flame thickness is correctly reproduced. A good estimation of near-wall thermal field (position P) allows to calculate properly the contribution of the heat convection near façade wall. Whereas, at each height, the temperatures are overestimated at position Q (cf. Fig. 8), corresponding to the intermittent flame zone. In the other hand, the experimental database in Ref. [START_REF] Kruppa | Comportement au feu des poteaux extérieurs en acier[END_REF] does not provide a complete description of the total uncertainty in the experimental temperature measurements. To the best knowledge of the authors, the response time of the thermocouples used in the experiment is close to a second by including thermal mass inertia and accumulation of soot particles. Thus, the uncertainty in the measurement of gas temperature by using thermocouples is estimated within 5-10% by taking into account a multitude of potential errors, including surface reactions, radiation, stem loss, etc. Besides, the humidity of pine type of wood [START_REF] Nolan | Gamme-radiographic study of wood combustion[END_REF] used in the tests [START_REF] Kruppa | Comportement au feu des poteaux extérieurs en acier[END_REF] was not experimentally identified. Consequently, use of the measured wood mass loss rate (cf. Fig. 4) in the numerical simulation may cause an overestimation of the HRR (cf. Fig. 5), which affects the predicted temperature distribution. It seems most likely that the discrepancies are due to a combination of the experimental uncertainties and the possible error in the numerical simulation by neglecting the water evaporation rate. Figure 10 shows the evolution of the predicted temperature, averaged from the 12 measurement points inside the compartment (cf. Fig. 3b). As expected, the predicted temperature with a mesh size of 20 cm is far away from the experimental data. Globally, the model appears to accurately estimate the gas temperature during the quasi steady state period between 400 and 1200 seconds with a grid size of 10 cm. These results show a similar trend to the findings in the Desanghere's work [START_REF] Desanghere | Détermination des conditions d'échauffement de structure extérieure à un bâtiment en situation d'incendie[END_REF] using the version FDS2 for the same compartment fire and experimental data. The subgrid model in FDS is known as dissipative with a mesh size of cm [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF], and lacks the ability to predict the transient phenomena under laminar-like conditions during the fire decay phase after 1500 s. A more rapid decrease of the predicted temperature inside the enclosure may be attributed to the overestimate of the fresh-air entrainment towards the enclosure during the fire extinction phase. Figure 10. Influence of the grid size on the mean temperature inside compartment with a difference varying from 50% to 5% for a mesh size of 20 cm and 10 cm, respectively.

Influence of the sidewalls on the flame ejecting behavior

In the numerical simulation, a time-dependent heat release rate, similar to the one in Fig. 5, is used with a prescribed peak value varying from 1 to 8 MW. Scaling laws [START_REF] Kolbrecki | Model of fire spread out on outer building surface[END_REF][START_REF] Kawagoe | Fire Behaviour in Rooms, Report of the Building Research Institute[END_REF] via the concept of Global Equivalence Ratio (GER), are proposed for comparing the full-scale simulations results with the correlations [START_REF] Lu | Facade flame heights from enclosure fires with side wall at the opening[END_REF][START_REF] Lu | Experimental investigation on window ejected facade flame heights with different constraint side wall lengths and global correlation[END_REF] derived from small-scale tests. The global equivalent ratio (GER), m / s m air fuel

   
, is calculated from the fuel supply rate (𝑚̇𝑓 𝑢𝑒𝑙 ), air inflow rate (𝑚̇𝑎 𝑖𝑟 ) and stoichiometric coefficient s.

The semi-empirical correlation [START_REF] Kawagoe | Fire Behaviour in Rooms, Report of the Building Research Institute[END_REF] is derived from a mass balance in conjunction with the experimental data which allows to adjust the constant related to the flow coefficient (0.41). At a certain opening geometry, the mass outflow rate (𝑚̇𝑜 𝑢𝑡 ) from the compartment is calculated as follows:

𝑚̇𝑜 𝑢𝑡 = 0.41𝜌 𝑎 𝐴√ℎ [START_REF] Sikanen | Predicting the heat release rates of liquid pool fires in mechanically ventilated compartments[END_REF] This implies that the height of the opening (h) has an influence more significant than the opening width (w). The mass inflow rate can be derived from the mass conservation at the opening:

𝑚̇𝑎 𝑖𝑟 = 𝑚̇𝑜 𝑢𝑡 -𝑚̇𝑓 𝑢𝑒𝑙 (8)
The dimensionless internal heat release rate (𝑄 ̇𝑖𝑛 * ) is defined as follows:

𝑄 ̇𝑖𝑛 * = 𝑄 ̇𝑖𝑛 𝑚̇𝑎 𝑖𝑟 𝑌 𝑂 2 ∆𝐻 𝑂 2 (9)
Here, the internal heat release rate (𝑄 ̇𝑖𝑛 ) is derived from the air entrainment inside compartment: Here, 𝑌 𝑂 2 and ∆𝐻 𝑂 2 are respectively the ambient mass fraction of oxygen and the amount of energy released per unit mass of oxygen consumed. The external heat release rate, ex Q  , is derived from the total heat release rate, tot Q  , and internal one, in Q  , as follows:

𝑄 ̇𝑖𝑛 = 𝑚̇𝑎 𝑖𝑟 𝑌 𝑂 2 ∆𝐻 𝑂 2 = (𝑚̇𝑜 𝑢𝑡 -𝑚̇𝑓 𝑢𝑒𝑙 )𝑌 𝑂 2 ∆𝐻 𝑂 2 (10)
in tot ex Q Q Q      (11)
Influence of the sidewalls distance (cf. Fig. 2) on the spill flame height due to burning the excess fuel gas ) Q ( ex  ejected from an opening was experimentally investigated by Lu [START_REF] Lu | Facade flame heights from enclosure fires with side wall at the opening[END_REF] from a reduced scale enclosure with a HRR below 60 kW. A global dimensionless parameter, K, is assumed to describe the mean flame height as a function of the separation distance D, the characteristic length scale of the opening size (h, w) in conjunction with the competition of momentum and buoyancy flux at the opening under ventilated condition. The parameter, K, represents the flame height (H D ) for a distance D normalized by the one (H 0 ) of a flat façade, as follows:

) L , D , h , w ( f H H K 0 D   ( 12
)
The detailed formulation for the parameter K can be found elsewhere [START_REF] Lu | Facade flame heights from enclosure fires with side wall at the opening[END_REF].

In order to establish various under-ventilated fire conditions, seven scenarios with different fire load and opening dimensions were selected, as summarized in Table 2. For each Case 1 to 7, the separation distance D varies from 2.6, 2.8, 3.6, 4.4, 6, 8 m up to  corresponding to a flat façade. The depth of the opposed sidewalls (cf. Fig. 2) is fixed at L=3 m, and this allows to focus solely on the influence of the separation distance on the flame ejecting behavior. Turbulent flames exhibit a pulsing behavior. A mapping flame luminosity technique [START_REF] Lu | Facade flame heights from enclosure fires with side wall at the opening[END_REF][START_REF] Kolb | Cross flow effects on the flame height of an intermediate scale diffusion flame[END_REF][START_REF] Coutin | Flame heights in wall fires : effects of width, confinement and pyrolysis length[END_REF] was developed to measure the mean flame height through a series of continuous images recorded using a CCD camera. The selected luminosity threshold corresponding to a continuous flame presence probability of 0.95 is used. It was checked that the so-determined persistent flame corresponds to a gas temperature of about 500°C [START_REF] Kolb | Cross flow effects on the flame height of an intermediate scale diffusion flame[END_REF][START_REF] Coutin | Flame heights in wall fires : effects of width, confinement and pyrolysis length[END_REF]. Note that an alternative for the description of the flame structure is based on the CH*, OH* radicals formation zone [START_REF] Legros | Soot volume fraction measurements in a three-dimensional diffusion flame established in microgravity[END_REF], but it is not applicable to a full-scale fire due to both the cost and difficulty in the measurement of such species.

LES approach provides an instantaneous gas temperature output,  i , and the time- averaged quantity,  ˆ over a range of the computational time is obtained as,

         N 1 i i N 1 i i i t t ˆ ( 13
)
For each case, the temperature is averaged over a steady state period of 600 s during a physical time of 1800 s for determining the flame shape. As an illustration, the predicted, time-averaged temperature fields for Case 4 are shown in Figure 11(a-c) for three separation distances. The air supply rate is much less than the stoichiometric requirements inside the enclosure, inducing a more important ejection of flammable fuel-rich vapour cloud from the opening. There is quasi-steady spill burning characterized by the two regions of the flame one inside the box near the fuel source, and the other outside the box. Outside the enclosure, the height, Hf, of the visible flame past the opening from the ground level, corresponds to the zone where the temperature is higher than 500°C which is of particular concern for the radiation heat flux [START_REF] Lee | Heat flux distribution and flame shapes on the inert façade[END_REF]. With an increase of the separation distance, D, the flame is detached from the façade wall, accompanied by a remarkable decrease of the flame height (cf. Fig. 11a). As the sidewalls are adjacent directly to the opening (cf. Fig. 11c), the entrainment of air mostly occurs from front of the opening and the flame is attached to the façade wall. Otherwise, the increase of entrainment from side of the opening makes the flame detachment from the façade wall (cf. Fig. 11b). The predicted ratio, K=HD/H0, is used to make quantitative comparison with the correlation [START_REF] Pope | Quantitative comparison of FDS and parametric fire curves with post-flashover compartment fire test data[END_REF] in Figure 12 as a function of the normalized separation distance D/w varying from 1 to 7. The predicted ratio, K, follows closely the experimental or correlation trend with a nonlinear relation between the ratio K and the separation distance. It is found that an increase of the separation distance between the sidewalls, D w , from 1 to 3 causes sensitively a reduction in the external flame height. A further increase of the distance between the sidewalls, D w , beyond 3 brings about a progressive decrease of the ratio K, consistent with the experimental observation [START_REF] Lu | Facade flame heights from enclosure fires with side wall at the opening[END_REF]. The influence of the sidewalls on the externally venting flame height becomes negligible for the normalized distance, D w , higher than 5. Globally, both the correlation (12) and the prediction match approximately the experimental trend within 10% uncertainty. The correlation [START_REF] Pope | Quantitative comparison of FDS and parametric fire curves with post-flashover compartment fire test data[END_REF] proved to be valid for a wide range of conditions using experimental data and are also validated by CFD simulation results for a full-scale façade fire. Note that these results are valid only for inert façade, and should not be directly extrapolated to fire propagation over a façade assembly including an insulation inflammable materials. In Eq.( 12) for the ratio K, the spill flame height of a flat façade, H0, can be determined using the correlation ( 14) which was developed from a reduced scale compartment fire in conjunction with the experimental data [START_REF] Lu | Experimental investigation on window ejected facade flame heights with different constraint side wall lengths and global correlation[END_REF]. This correlation is typically a function, f, related to the opening factor of a compartment (w, h) and the dimensionless external heat release rate,

Q * ex  . ) Q , h , w ( f H * ex 0   (14)
The detailed formulation for H0 can be found elsewhere [START_REF] Lu | Experimental investigation on window ejected facade flame heights with different constraint side wall lengths and global correlation[END_REF]. The parameter Q * ex  is related also to the opening factor and the ambient air properties (𝜌 ∞ , 𝐶 𝑝 , 𝑇 ∞ ) as follows:

gH A T C Q Q p ex * ex       (15)
The external heat release rate, ex Q  , due to the burning excess fuel gas ejected from opening, is determined from Eq.(11).

In the numerical simulation, a range of 1 to 3 m with a regular interval of 0.2 m is chosen for the opening height (h)/width (w) in addition to a HRR ( Q  ) varying from 1 to 8 MW with a uniform interval of 0.5 MW. Different fire scenarios were simulated by combining the three parameters such as Q  , h and w for providing various external HRR values which influence directly the façade fire development. A few numerical data points compare with available correlation results in Fig. 13 for identifying the façade flame height dependence on the dimensionless external HRR. It is clear that the correlation ( 14) ignores the flame appearance outside the compartment when ex Q  tends to zero for GER below 1. Whereas, it was experimentally observed [START_REF] Lu | Facade flame heights from enclosure fires with side wall at the opening[END_REF][START_REF] Lu | Experimental investigation on window ejected facade flame heights with different constraint side wall lengths and global correlation[END_REF][START_REF] Tang | Experimental study on flame height and temperature profile of buoyant window spill plume from an under-ventilated compartment fire[END_REF] that when the fire is still fuel controlled, i.e. the GER is close but less than unity, the further evolved flames in the ceiling jet may become long enough to eject from the compartment opening. The numerical model is found to be able to predict the presence of the external flame in the situation for GER below 1. Globally, the overall evolution of predicted flame height follows closely the correlation trend as a function of the external HRR. Note that the numerical model exhibits a more conservative behavior with an overprediction of 40% as compared to the correlation [START_REF] Asimakopoulou | Thermal characteristics of externally venting flames and their effect on the exposed façade surface[END_REF]. This may be attributed to the overestimation of the gas temperature (cf. Fig. 8), which affects the externally venting flame behaviour via acceleration of the buoyancy-induced flow along the façade. Moreover, a similitude between a full scale and a reduced scale compartments cannot be totally preserved due to the confinement effect [START_REF] Lassus | Ventilation effects in confined and mechanically ventilated fires[END_REF] influencing the turbulence development of the ejected fire. Besides, the discrepancy may partly be due to different identification of the visible flame shape in the experiment and the numerical simulation. As an illustration in Figure 14(a,b) for a heat release rate of 1.5 MW, the opening geometries affect not only the flame height, but also the fire plume structure. With a reduction of opening area from 4.2 m 2 (w=3.0 m, h=1.4 m ) to 3.6 m 2 (w=2.0 m, h=1.8 m), the flame (cf. Fig. 14b) is detached from the façade wall with a remarkable volume growth of the thermal plume. The excess heat release rate of the fuel burning outside the opening with a reduced area becomes high enough to produce flames controlled by three-dimensional air entrainment. It is found from Fig. 14 that the opening geometries control the supply of air towards the compartment, influencing directly the excess heat release rate inside enclosure. Figure 15 shows the evolution of dimensionless internal HRR, 𝑄 ̇𝑖𝑛 * (Eq.9) as a function of the GER from both the numerical simulation and the correlation. As for the flame height, H0 (cf. Fig. 13), a series of fire simulations were conducted by combining also the three parameters such as Q  , h and w for providing various global equivalent ratios (GER) which influence significantly the sustained internal combustion and unburnt volatiles formation. A few numerical data points compare with available correlation results in Fig. 15 for identifying the dependence of the dimensionless internal HRR on GER. The experimental investigations [START_REF] Lu | Facade flame heights from enclosure fires with side wall at the opening[END_REF][START_REF] Lu | Experimental investigation on window ejected facade flame heights with different constraint side wall lengths and global correlation[END_REF][START_REF] Tang | Experimental study on flame height and temperature profile of buoyant window spill plume from an under-ventilated compartment fire[END_REF][START_REF] Kolbrecki | Model of fire spread out on outer building surface[END_REF][START_REF] Kawagoe | Fire Behaviour in Rooms, Report of the Building Research Institute[END_REF] indicate that with GER below 1, there is enough air supply for burning all the released fuel inside enclosure, corresponding to an over-ventilated fire. The unity of GER implies that the amount of air corresponds exactly to the amount needed to react with all the fuel released from the pyrolysis. The 𝑄 ̇𝑖𝑛 * below 1 implies that a little part of the combustible gases is ejected through the opening and burnt outside. With GER higher than 1, air inside the compartment no longer allows combustion of all pyrolysis gases. The fundamental research [START_REF] Tewarson | Generation of Heat and Chemical Compounds in Fires[END_REF][START_REF] Babrauskas | Ignition Handbook[END_REF] shows that fire exhaust occurs as fuel-air mixture reaches a rich limit of flammability. The correlation (9) ignores the fire dynamics under lowered oxygen vitiation enclosure. By using this type correlation, at a certain opening geometry, the critical internal HRR is limited only by air supply rate, and it remains constant with a unity value even the fuel-air mixture is beyond the rich limit of flammability. In fact, the correlation requires that the oxygen concentration in the outflow is uniform (or well-mixed), which is the case for fully developed (post flashover) fires. For under- 9) may result in an overestimation of the HRR due to the lower uniformity of oxygen concentration in the outflow. In the numerical simulation, when the value of GER is close to 3, the fire source inside the compartment is no longer sufficiently supplied with fresh air. The predicted GER value of 3 corresponds to the rich limit of flammability of wood fuel-air mixture beyond which fire exhaust occurs inside compartment with a significant reduction in HRR (cf. Fig. 15). Indeed, the predicted critical GER value for fire exhaust is somewhat below to the Tewarson's data [START_REF] Tewarson | Generation of Heat and Chemical Compounds in Fires[END_REF] primarily from the fire propagation apparatus. It is believed that the Tewarson's data may be not completely applicable to a compartment with an opening that is elevated above the floor. In fact, fire exhaust inside enclosure remains quite complex since it depends strongly on the temperature-dependent reaction kinetics. The EDC type is temperature independent combustion model, and chemical reaction takes place at a stoichiometric fuel-air mixture. So the predicted critical GER value of 3 is qualitatively promising as compared to the experimentally determined one [START_REF] Tewarson | Generation of Heat and Chemical Compounds in Fires[END_REF].

Figure 15. Evolution of the dimensionless internal HRR by combining the three parameters such as Q  , h and w with a discrepancy below 5% between the prediction and correlation for GER below 3.

As an illustration, Case 4 (cf. Tab.2) is chosen for understanding evolution of the compartment mean temperature and the internal heat release rate as a function of the normalized separation distance, D/w, varying from 1 to 4. According to this numerical study, the fire behavior inside compartment with the seven scenarios (cf. Tab.2) is substantially identical to Case 4. The temperature inside the enclosure is spatially uniform (cf. Fig. 14). As shown in Figure 16, the presence of the sidewalls positioned at the opening does not affect the mean temperature inside the compartment. This trend is similar to the experimentally observed one [START_REF] Lu | Facade flame heights from enclosure fires with side wall at the opening[END_REF][START_REF] Lu | Experimental investigation on window ejected facade flame heights with different constraint side wall lengths and global correlation[END_REF]. The façade fire ejecting behavior without sidewalls is closely related to the air entrainment from both front and side. With presence of the sidewalls adjacent directly to the opening, the fresh-air entrainment from side towards the compartment is fully constrained. Thus, the ventilation controlled internal HRR for the distance of D/w=1 is substantially lower than the one for the separation distance (D/w) beyond 1. This induces an increase in the amount of combustible gases ejected from the opening, and consequently in flame height (cf. Fig. 12). Whereas, the heat released inside the compartment varies slightly with a further increase of the sidewalls separation distance (D/w>1.5). This suggests that the fresh air inflow rate from front and side is not affected once the sidewalls are far away from the opening, and the same ventilation-controlled fire regime inside the compartment is maintained. In fire plume, the majority of the radiation is derived from the visible part of the flame with a temperature higher than 500°C, where soot particles are radiating heat [START_REF] Kolb | Cross flow effects on the flame height of an intermediate scale diffusion flame[END_REF][START_REF] Rasbash | Fundamentals of smoke production[END_REF]. The radiation heat flux is a function of both the gas temperature with the T 4 dependence and chemical composition (soot, CO2, H2O). Neither of which can be reliably calculated in a fullscale fire calculation because the flame sheet is not well-resolved on a relatively coarse numerical grid. The measured radiative heat flux is not available for a full-scale façade fire with sidewalls due to both the cost and difficulty in the measurement of heat flux. The numerical simulations appear useful to evaluate qualitatively the modes of heat transfer.

In FDS [START_REF] Mcgrattan | Fire Dynamics Simulator (Version 6), user's guide[END_REF], soot production is derived solely from the fraction of the fuel mass that is converted into soot. Such soot yield approach does not incorporate the essential physical processes of soot inception, coagulation and surface growth in addition to the influences of turbulent fluctuations, temperature and fuel type [START_REF] Moss | Modelling soot formation and burnout in a high temperature laminar diffusion flame burning under oxygen-enriched conditions[END_REF][START_REF] Lindstedt | Soot formation in Combustion[END_REF]. For the pine type of wood, the soot yield value of 0.015 is available [START_REF] Tewarson | Generation of Heat and Chemical Compounds in Fires[END_REF][START_REF] Rasbash | Fundamentals of smoke production[END_REF] only for well-ventilated fires. Therefore, sensitivity of the soot yield value on the radiative predictions at a given elevation of 5.5 m along the façade wall was first evaluated for a under-ventilated fire as Case 4. Two specific soot yield values of 0.01 and 0.04 were chosen to evaluate its sensitivity on the predicted radiative heat flux in Figure 17. For a given soot yield value, the oscillation behavior is attributed to the buoyancyinduced instability. Use of two soot-yield values induces the different transient peak values for the radiation heat flux via a radiative transfer equation. It is found that with multiplying the soot yield value by a factor of 4, the radiation heat flux increases by a factor of 30% with a rise of the peak from 15 to 20 kW/m². Therefore, use of such soot yield approach may induce an uncertainty within 30% by taking into account the transition from well-to under-ventilated enclosure fire.

As an illustration, the effect of reducing the separation distance between sidewalls on the radiation heat flux to the façade surface at centerline is examined in Figure 18 for Case 4 using a soot yield value of 0.015 for wood. Evolution of the radiation flux along the flat façade ) D (   shows a similar trend to the Chotzoglou's finding [START_REF] Chotzoglou | Experimental investigation of externally venting flames geometric characteristics and impact on the façade of corridor-like enclosures[END_REF] from a small-scale corridor-like enclosure. As expected, radiation energy absorbed by the façade wall depends mainly on the size of the fire and the soot loading. The peak of the heat flux from the external flame to the façade wall rises from 70 to 90 kW/m 2 with a reduction in the separation distance. Presence of the sidewalls adjacent closely to the opening (D=2.8 m) clearly enhances the radiative heat transfer with the strongest heat flux close to the opening due to a large flame thickness there. In all the cases, it is ensured that if the façade assembly contains inflammable materials, the heat flux is high enough for inducing an ignition of façade fuel package exposed to the external flame near the opening. Particularly, a reduction in the distance between the sidewalls will significantly enhance the fire spreading rate and its size due to the increasing contribution by radiation. An experimental database on heat flux, owing to pyrolysis of an inflammable façade, over the full-scale façade with sidewalls constrains at the opening is needed for validation of the numerical tool.

CONCLUSIONS

A series of numerical simulations was carried out for highlighting the flame ejecting behavior of full-scale enclosure fires with sidewalls constrains at the opening. The numerical model exhibits a more conservative behavior (over-prediction) in the intermittence zone of the spill flame which coincides with the periodic forward pulsations of flame due to the buoyant instability. Nevertheless, the numerical tool, is useful after detailed validation for full-scale fire to evaluate the correlations which are commonly derived from small-scale tests. The height of externally venting flame with sidewalls is much higher than that of flat façade due to a severe constrain of the air entrainment from side. The predicted spill flame height in full-scale is generally higher than that in small-scale from the correlation [START_REF] Asimakopoulou | Thermal characteristics of externally venting flames and their effect on the exposed façade surface[END_REF], owing probably to the turbulence development. The correlation (9) for the internal HRR might be used only in the case where the GER is lower than 3 corresponding to the rich limit of flammability of fuel-air mixture. The predicted dimensionless spill flame height as a function of sidewalls distance are shown to be well collapsed by the correlation [START_REF] Pope | Quantitative comparison of FDS and parametric fire curves with post-flashover compartment fire test data[END_REF]. An inverse dependence of the flame height with the separation distance of sidewalls is predicted, following the correlation or experimental trend. Subsequently, the correlation [START_REF] Pope | Quantitative comparison of FDS and parametric fire curves with post-flashover compartment fire test data[END_REF] for dimensionless façade flame heights with sidewalls constrains can be used to enhance fire safety design methodologies for U-shape building façade in spite of scale effect.

Figure 1 .Figure 2 .

 12 Figure 1. Real fire propagation of vertical U-shape building façade in France (left) and in China (right) both occurred in 2010

Figure 3 .

 3 Figure 3. General layout and characteristic configuration of a full-scale façade fire test: a) Front view for the thermocouple locations; b) Side view for the thermocouple locations.

Figure 4 .

 4 Figure 4. Evolution of the measured mass loss rate of the pine type of wood used as input data in numerical simulation

Figure 5 .

 5 Figure 5. Predicted total heat release rate and internal one inside the compartment based on the measured mass loss rate

Figure 6 .

 6 Figure 6. Comparison between the thermocouple model and gas temperature

Figure 7 .Figure 8 .TimeFigure 9 .

 789 Figure 7. Influence of the mesh size varying from 5, 10 to 20 cm on the predicted external temperatures at position P.

Figure 11 .

 11 Figure 11. Distribution of the predicted time-averaged gas temperature on the vertical symmetric plane at y=1.8 m for Case 4.

Figure 12 .

 12 Figure 12. Comparison of dimensionless flame height between correlation from small-scale tests and numerical simulation from a full-scale compartment as a function of the normalized separation distance (D/w) within 10% difference

Figure 13 .

 13 Figure 13. Comparison of the flame height between prediction and correlation as a function of the dimensionless external HRR within 40% difference by combining the three parameters such as Q  , h and w.

Figure 14 .

 14 Figure 14. Influence of the opening geometries as h and w on distribution of the predicted timeaveraged gas temperature for a heat release rate of 1.5 MW.

Figure 16 .

 16 Figure 16. Evolution of the compartment mean temperature and heat release rate inside the enclosure as a function of the normalized separation distance (D/w) for Case 4

Figure 17 .

 17 Figure 17. Comparison of the radiative heat flux over façade surface at a height of 5.5 m for Case 4 with a difference of 30% between the soot yield values of 0.01 and 0.04

Figure 18 .

 18 Figure 18. Evolution of heat flux along the façade centerline surface for three separation distances (D) within 30% uncertainty for Case 4.

Table 1 :

 1 Material thermo-physical properties for the façade and enclosure walls

	Parameters	Thermal	Density	Thermal	Thickness (m)
		conductivity	(kg/m 3 )	specific capacity
		(W/m.K)		(J/kg.K)
	Façade	0.2	450	1000	0.2
	Ceiling	0.2	450	1000	0.2
	Wall	0.45	1050	840	0.2
	Floor	1.6	2300	1000	0.2

Table 2 .

 2 Summary of simulation scenarios with different opening dimension and HRR

	Parameter	Case 1	Case 2	Case 3	Case 4	Case 5	Case 6	Case 7
	w (m)	2.0	2.0	2.0	2.6	2.6	2.6	3
	h (m)	1.8	1.4	1.4	1.4	1.5	1.4	1.4
	HRR (MW) 6.5	6.5	5.5	6.5	6.5	5.5	6.5