A Study on Heat Release Rates of Furniture Under Well-Developed Fire

With the increasing interest in well-developed fires, the burning behavior of furniture under flashover should be studied. Full-scale burning tests on the heat release rates of furniture under flashover with oxygen calorimetry will be reported in this article.

Nine tests on sofas with and without the treatment of fire retardant, wood desks, and the fire source itself to start the burning were carried out. The heat release rate, oxygen concentration, thermal radiative heat fluxes at floor level, and air temperatures at some selected positions were measured.

Response of furniture under small accidental fires and bigger pool fires to onset flashover will be discussed. Obviously, burning furniture under a flashover fire will give a much higher heat release rate. Sofas treated with fire retardant might burn even more vigorously. This point should be considered carefully in designing fire-safe furniture.     

Towards large eddy simulations of flame extinction and carbon monoxide emission in compartment fires

The general objective of this research is to adapt current combustion modeling capabilities used in computational fluid dynamics solvers to the treatment of under-ventilated compartment fires. More specifically, we consider in the present study two models proposed to describe: diffusion flame extinction due to air vitiation; and the emission of carbon monoxide (CO) and unburnt hydrocarbon (HC) mass in a compartment fire. The flame extinction model is based on a flammability diagram parametrized in terms of vitiated air properties. The CO/HC mass model is based on: a transport equation for fuel mass; a comparison of this fuel mass to a Burke–Schumann chemical-equilibrium expression; and an interpretation of the difference as a measure of incomplete combustion. Both models are implemented into a large eddy simulation solver developed by the National Institute of Standards and Technology, USA. The models performance is tested via detailed comparisons with an experimental database corresponding to reduced-scale compartment fires. The study considers two cases that correspond to different values of the fire room global equivalence ratio and are representative of strikingly different flame behaviors. The comparative tests serve to evaluate the general ability of the models to describe the transition from extinction-free conditions to conditions in which the flame experiences partial or total quenching, as well as the transition from fire regimes with no or little CO emission to regimes that emit hazardous levels.     

Behaviour of a confined fire located in an unventilated zone

The behaviour of a fire in an enclosure is studied for a configuration where the fuel source is located in the upper hot unventilated zone trapped by a soffit. The experimental study, undertaken in a laboratory-scale compartment with a fuel source above the level of a soffit, included the determination of the parameters (ventilation factor, rate of fuel supply) controlling the combustion or leading to extinction. Measurements (PIV, thermocouples, gas sampling and analysis) were performed to propose a hypothesis on the structure of the flame (flame stabilisation mechanisms, premixed or diffusion types). Video photography is used to determine the area covered by the flames. This information is used as a criterion to identify the combustion regimes. The results show that the gaseous fuel is diluted in the combustion products (CP) in the upper layer and that a recirculatory motion is formed, driven by buoyancy forces, which enhances the mixing of fuel and CP. These then travel horizontally towards the vent along the interface between the lower fresh air and upper zones, and are premixed with the convected air in the enclosure, before entering the reaction zone and being burnt. The flame stabilises at the interface between the upper hot and lower ventilated layers in the compartment. The observed “ghosting flame” is stabilised by a triple flame if the flame speed of the premixed flame is higher than the natural convection velocity induced in the compartment. The flame stability is quantified by a criterion based on the area of the horizontal flame. It has been observed that the combustion is controlled by the available mass fuel flux at the reaction zone if the ventilation is sufficient. This information is essential for the modelling of the phenomena involved in fires with such an underventilated fuel source.     

Nonlinear simulations of combustion instabilities with a quasi-1D Navier-Stokes code

As lean premixed combustion systems are more susceptible to combustion instabilities than non-premixed systems, there is an increasing demand for improved numerical design tools that can predict the occurrence of combustion instabilities with high accuracy. The inherent nonlinearities in combustion instabilities can be of crucial importance, and we here propose an approach in which the one-dimensional (1D) Navier-Stokes and scalar transport equations are solved for geometries of variable cross-section. The focus is on attached flames, and for this purpose a new phenomenological model for the unsteady heat release from a flame front is introduced. In the attached flame method (AFM) the heat release occurs over the full length of the flame. The nonlinear code with the use of the AFM approach is validated against analytical results and against an experimental study of thermoacoustic instabilities in oxy-fuel flames by Ditaranto and Hals [Combustion and Flame 146 (2006) 493-512]. The numerical simulations are in accordance with the experimental measurements and the analytical results and both the frequencies and the amplitudes of the resonant acoustic pressure modes are reproduced with good accuracy.     

Temperature evolution,flame behavior and extinction inside compartment with a ceiling vent under the ambient wind

This article presents an experimental and CFD simulation investigation and analysis on temperature evolution and flame behavior inside compartment with a ceiling vent under ambient wind. Experiments were conducted employing a reduced-scale model containing a cubic fire compartment with a ceiling vent under the external wind generated by a wind tunnel. The temperature of windward and leeward inside compartment as well as the flame behavior were recorded for various vent dimensions, heat release rates and wind speeds for a total of 720 test conditions. Results show that there are two types of fire behavior regimes inside the compartment for relatively small- (Bernoulli flow regime) or large vents (oscillatory exchange flow regime). With a relatively small vent, the flame is located at the center of compartment and the temperatures of windward and leeward are almost the same with or without the ambient wind. Then the flame turns to extinct due to lack of oxygen with increasing of fuel supply. However, with a relatively large vent, there is a flame transition from windward to leeward accompanied by the transition of temperature distribution with increasing of fuel supply when subject to ambient wind. Such extinction and transition mechanisms are interpreted by the aid of CFD simulation of global equivalence ratio and flow/oxygen field inside compartment. The critical heat release rate for the occurrence of flame transition decreases with raising wind speed while increases with vent size. Their complex dependence is found to be well represented in terms of a non-dimensional heat release rate as a function of the wind Froude number, employing the vent area-equivalent characteristic diameter as length scale. These new findings facilitate the understanding of the compartment fire evolution with a ceiling vent subject to ambient wind.     

Attenuation of combustion instability in a fuel-staged dual-nozzle gas turbine combustor with asymmetric hydrogen composition

The instability attenuation mechanism of fuel staging was investigated in a CH₄/H₂ fueled dual-nozzle gas turbine combustor. Fuel staging was implemented using an asymmetry in fuel composition between the two nozzles. The fuel composition of the upper nozzle was varied while keeping that of the lower nozzle constant. Under these conditions, the self-excited and forced responses of fuel-staged flames were analyzed using OH* chemiluminescence imaging, OH planar laser-induced fluorescence, and particle image velocimetry. In the self-excited measurements, although strong combustion instability was exhibited in the symmetric condition, it weakened gradually with increasing asymmetry in fuel composition. The symmetric flame exhibited significant fluctuations in the heat release rate around the flame tip, which acted as the primary cause of driving combustion instability. However, in asymmetric flames, the H₂ addition induced phase leads in heat release rate fluctuations at the upper region, which damped combustion instability. Thus, our observations revealed a high correlation between the phase leads and the attenuation of combustion instability. Analyses of the forced responses showed that the heat release rate fluctuations were induced by interactions between the flame and the shedding vortex released from the nozzle tip into the downstream. Although these characteristics of shedding vortices did not depend on the H₂ addition, the change in the axial position of the flame caused by the H₂ addition induced the relocation of the site, at which the flame interacted with the vortex. Subsequently, it induced phase leads in the heat release rate fluctuations. The phase difference of heat release rate fluctuations between the two flames due to this phase leads enlarged progressively with increasing asymmetry in fuel composition, leading to the attenuation of combustion instability in asymmetric conditions.     

An experimental study on flame geometry and radiation flux of line-source fire over inclined surface

An experimental study was performed on line-source fire over an inclined surface (ground) to simulate downhill fire spread behavior. The flame geometry and the thermal radiation to both far-field surroundings and near-field inclined surface were investigated. As a basic configuration for wildland fire over a slope, the buoyancy induced natural convection flow along the inclined surface and the constraint of air entrainment by the inclined surface change the flame geometry as well as its radiation emission. Various surface (ground) inclination angles (from 0°-80°), fire source heat release rates and fuels were considered comprehensively with a total of 126 test conditions. Results showed that the flame perpendicular height decreased, while both the flame parallel length and base drag length along the inclined surface increased, with the increased inclination angle. A dimensional analysis was then performed based on the controlling mechanisms, with the dimensionless heat release rate, the density ratio of fuel vapor to air, along with sinα and cosα involved to represent the components in the parallel and perpendicular directions. The flame geometry parameters were well represented by the proposed dimensional analysis. Both the radiation fluxes to far-field surroundings and to near-field inclined surface decreased with the increased inclination angle. The far-field radiation was found to be well characterized by a model based on the soot volume fraction analysis according to single point source model. Concerning the near-field radiation to inclined surface, an inclined cuboid radiative modeling was developed. The predicting results by the proposed model and the experimental values showed good agreement. The present study has explained the controlling physics and proposed non-dimensional functions for flame geometry and modeling the downslope radiation of the line-source fire over the inclined surface, which facilitates the understanding of the wildland fire spread behavior over a slopping ground in the downhill direction.     

An investigation of oscillating instabilities in deflagrations with edges

We examine the Lewis-number-greater-than-1 stability of a deflagration sitting on a porous-plug burner with an inert coflow. The flame edges generated by the coflow influence the stability, and this influence is examined. Very wide flames display the same stability characteristics as unbounded flames (flames sans edges), but for moderately wide flames the instability is suppressed. A new two-dimensional instability can occur for narrow flames. There is a range of mass fluxes for which a monotonic decrease in burner (flame) width generates a transition from unstable flames to stable flames, to unstable flames, to quenching. The insertion of a cold probe into the combustion field can stabilize an unstable flame or destabilize a stable flame, depending on the point of insertion.     

Testing the acoustical corrections for reflections on a façade

The assessment of noise levels, in the proximity of a building or on its façade, is a requirement of the European Environmental Noise Directive 2002/49/EC concerning environmental noise produced by road and railway traffic, airports and industries. The corrected values for the noise levels and the problems of measuring noise near a façade are discussed here for the case of road traffic. A complete set of measurements in different situations was performed along several roads in an urban environment as well as in a controlled situation using a loudspeaker. The experimental results are then compared against theoretical models and international standards, in particular those that suggest to use a +3 or +6 dB correction as a function of the microphone position and the NordTest method. Some suggestions are given for the different corrective factors to apply when measuring environmental noise between 0 and 2 m away from a building façade, and practical solutions identified.     

On the importance of the heat release rate in numerical simulations of fires in mechanically ventilated air‐tight enclosures

A fire in a compartment with limited ventilation can cause a significant pressure rise, up to hundreds of Pascal. This is important in practice, as the pressure rise can cause damage or hinder evacuation, but also from the perspective of fire safety science. From the energy balance, taking into account the interaction between compartment pressure, fire dynamics and mechanical ventilation, the importance of the net heat gained per unit time in the gas phase is well recognized. This leads to the need to accurately quantify the heat release rate inside the compartment as a function of time. It is explained that scaling of the transient phenomena is not straightforward. The paper then focuses on numerical simulations, in particular on CFD in the gas phase. An overview is presented of different existing approaches for turbulent combustion modelling in turbulent buoyancy-driven flames with low values of scalar dissipation rate, typical for fire flames. A dynamic approach for modelling turbulent combustion, and the coupling with radiation modelling, is briefly discussed. Extinction and re-ignition are discussed extensively, in the context of reduced ventilation conditions. Finally, low-frequency oscillatory behaviour in mechanically ventilated air-tight compartments is addressed. It is argued that CFD simulations are a very valuable tool to gain further insight in this phenomenon. Suggestions for future research are formulated.     

Modelling thermal radiation in buoyant turbulent diffusion flames

This work focuses on the numerical modelling of radiative heat transfer in laboratory-scale buoyant turbulent diffusion flames. Spectral gas and soot radiation is modelled by using the Full-Spectrum Correlated-k (FSCK) method. Turbulence-Radiation Interactions (TRI) are taken into account by considering the Optically-Thin Fluctuation Approximation (OTFA), the resulting time-averaged Radiative Transfer Equation (RTE) being solved by the Finite Volume Method (FVM). Emission TRIs and the mean absorption coefficient are then closed by using a presumed probability density function (pdf) of the mixture fraction. The mean gas flow field is modelled by the Favre-averaged Navier–Stokes (FANS) equation set closed by a buoyancy-modified k-? model with algebraic stress/flux models (ASM/AFM), the Steady Laminar Flamelet (SLF) model coupled with a presumed pdf approach to account for Turbulence-Chemistry Interactions, and an acetylene-based semi-empirical two-equation soot model. Two sets of experimental pool fire data are used for validation: propane pool fires 0.3 m in diameter with Heat Release Rates (HRR) of 15, 22 and 37 kW and methane pool fires 0.38 m in diameter with HRRs of 34 and 176 kW. Predicted flame structures, radiant fractions, and radiative heat fluxes on surrounding surfaces are found in satisfactory agreement with available experimental data across all the flames. In addition further computations indicate that, for the present flames, the gray approximation can be applied for soot with a minor influence on the results, resulting in a substantial gain in Computer Processing Unit (CPU) time when the FSCK is used to treat gas radiation.     

Flame spread and interactions in an array of thin solids in low-speed concurrent flows

Flame spread in an array of thin solids in low-speed concurrent flows was investigated and numerical solved. A previous steady, two-dimensional flame-spread model with flame radiation was employed and adapted in this work. The flame structures of spreading flames between parallel solids were demonstrated and some of the features were presented, including flow channelling effect and flame radiation interactions. The channelling effect is caused by flow confinement by the presence of the other solids; the flows through the hot combustion gases are accelerated downstream drastically. Radiation interactions between flames and solids contributed to a less heat-loss system, and radiation re-absorption by flames resulted in a larger flame with higher temperature, which increased the conductive heat fluxes to the solids and flame spread rate. Consequently, the extinction limit for the interacting flames is extended beyond the low-speed quenching limit for a single flame. The influence of the separation distance on the flame spread rate was also studied, which exhibits a non-monotonic behaviour. At larger separation distance, the flame spread rate increases with decreasing the separation distance owing to the channelling effect and radiation interactions. However, at very small separation distance, the flame spreading rate decreases with decreasing the distance because of the limited space for thermal expansion and flow résistance between solids.     

Thermal Radiation Modelling in Tunnel Fires

Modelling based on Computational Fluid Dynamics (CFD) is by now effectively used in fire research and hazard analysis. Depending on the scenario, radiative heat transfer can play a very important role in enclosure combustion events such as tunnel fires. In this work, the importance of radiation and the effect of the use of different approaches to account for it were assessed. Firstly, small-scale tunnel fire simulations were performed and the results compared with experimental data, then realistic full-scale scenarios were simulated. The results show up the capability of CFD modelling to reproduce with good approximation tunnel fires. Radiation proved to be noteworthy mainly when the scale of the fire is relatively large. Among the various approaches employed to simulate radiation, the use of the Discrete Transfer model gave the most accurate results, mainly when the absorption-emission characteristics of the combustion products were taken into account. Finally, the suitability of the use of CFD in quantitative Fire Hazard Analysis is discussed.     

The two-dimensional structure of low strain rate counterflow nonpremixed-methane flames in normal and microgravity

The structure and extinction of low strain rate nonpremixed methane–air flames was studied numerically and experimentally. A time-dependent axisymmetric two-dimensional (2D) model considering buoyancy effects and radiative heat transfer was developed to capture the structure and extinction limits of normal gravity (1-g) and zero gravity (0-g) flames. For comparison with the 2D modelling results, a one-dimensional (1D) flamelet computation using a previously developed numerical code was exercised to provide information on the 0-g flames. A 3-step global reaction mechanism was used in both the 1D and 2D computations to predict the measured extinction limit and flame temperature. Photographic images of flames undergoing the process of extinction were compared with model calculations. The axisymmetric numerical model was validated by comparing flame shapes, temperature profiles, and extinction limits with experiments and with the 1D computational results. The 2D computations yielded insight into the extinction mode and flame structure. A specific maximum heat release rate was introduced to quantify the local flame strength and to elucidate the extinction mechanism. The contribution by each term in the energy equation to the heat release rate was evaluated to investigate the multi-dimensional structure and radiative extinction of the 1-g flames. Two combustion regimes depending on the extinction mode were identified. Lateral heat loss effects and multi-dimensional flame and flow structure were also found. At low strain rates in 1-g flames (‘regime A’), the flame is extinguished from the weak outer edge of the flame, which is attributed to a multi-dimensional flame structure and flow field. At high strain rates, (‘regime B’), the flame extinction initiates near the flame centreline owing to an increased diluent concentration in the reaction zone, similar to the extinction mode of 1D flames. These two extinction modes can be clearly explained by consideration of the specific maximum heat release rate.     

Simulating fire whirls

A numerical investigation of swirling fire plumes is pursued to understand how swirl alters the plume dynamics and combustion. One example is the ‘fire whirl’ which is known to arise naturally during forest fires. This buoyancy-driven fire plume entrains ambient fluid as heated gases rise. Vorticity associated with a mechanism such as wind shear can be concentrated by the fire, creating a vortex core along the axis of the plume. The result is a whirling fire. The current approach considers the relationship between buoyancy and swirl using a configuration based on fixing the heat release rate of the fire and imposing circulation. Large-eddy methodologies are used in the numerical analyses. Results indicate that the structure of the fire plume is significantly altered when angular momentum is imparted to the ambient fluid. The vertical acceleration induced by buoyancy generates strain fields which stretch out the flames as they wrap around the nominal plume centreline. The whirling fire constricts radially and stretches the plume vertically.(Some figures in this article are in colour only in the electronic version; see www.iop.org)     

Detailed measurements of local scalar front structures in stagnation-type turbulent premixed flames

Local scalar front structures of OH mole fraction, reaction progress variable, and its three-dimensional gradient have been measured in stagnation-type turbulent premixed flames. The reaction progress variable front is observed to change with increasing turbulence from parallel iso-scalar contours but reduced progress variable gradients, called the lamella-like front, to disrupted non-parallel iso-contours that deviate substantially from those of wrinkled laminar flamelets, called the non-flamelet front. This transition is attributed to the different scales of interaction between the flame internal structure and a spectrum of turbulence extending from the integral scale to the Kolmogorov scale. The lamella-like front pattern occurs when the length scales of interaction are smaller than the laminar flame thickness but the time scales are greater than the flame residence time. The non-flamelet front pattern occurs when the length scales of interaction are greater than the laminar flame thickness but the time scales are smaller than the flame residence time. This difference corresponds to the change of combustion regime from complex-strain flame front to turbulent flame front on a revised regime diagram. A correlation is also proposed for the turbulent flame brush thickness as a function of turbulent Reynolds number and heat release parameter. The heat release parameter is considered to arise from the non-passive effects of flame-surface wrinkling.     

Concept and combustion characteristics of ultra-micro combustors with premixed flame

Under micro-scale combustion influenced by quenching distance, high heat loss, shortened diffusion characteristic time, and flow laminarization, we clarified the most important issues for the combustor of ultra-micro gas turbines (UMGT), such as high space heating rate, low pressure loss, and premixed combustion. The stability behavior of single flames stabilized on top of micro tubes was examined using premixtures of air with hydrogen, methane, and propane to understand the basic combustion behavior of micro premixed flames. When micro tube inner diameters were smaller than 0.4 mm, all of the fuels exhibited critical equivalence ratios in fuel-rich regions, below which no flame formed, and above which the two stability limits of blow-off and extinction appeared at a certain equivalence ratio. The extinction limit for very fuel-rich premixtures was due to heat loss to the surrounding air and the tube. The extinction limit for more diluted fuel-rich premixtures was due to leakage of unburned fuel under the flame base. This clarification and the results of micro flame analysis led to a flat-flame burning method. For hydrogen, a prototype of a flat-flame ultra-micro combustor with a volume of 0.067 cm³ was made and tested. The flame stability region satisfied the optimum operation region of the UMGT with a 16 W output. The temperatures in the combustion chamber were sufficiently high, and the combustion efficiency achieved was more than 99.2%. For methane, the effects on flame stability of an upper wall in the combustion chamber were examined. The results can be explained by the heat loss and flame stretch.     

Problems of predicting turbulent burning rates

Different approaches to the modelling of turbulent combustion first are reviewed briefly. A unified, stretched flamelet approach then is presented. With Reynolds stress modelling and a generalized probability density function (PDF) of strain rate, it enables a source term, in the form of a probability of burning function, P_b, to be expressed as a function of Markstein numbers and the Karlovitz stretch factor. When P_b is combined with some turbulent flame fractal considerations, an expression is obtained for the turbulent burning velocity. When it is combined with the profile of the unstretched laminar flame volumetric heat release rate plotted against the reaction progress variable and the PDF of the latter, an expression is obtained for the mean volumetric turbulent heat release rate. Through these relationships experimental values of turbulent burning velocity might be used to evaluate P_b and hence the CFD source term, the mean volumetric heat release rate.

Different theoretical expressions for the turbulent burning velocity, including the present one, are compared with experimental measurements. The differences between these are discussed and this is followed by a review of CFD applications of these flamelet concepts to premixed and non-premixed combustion. The various assumptions made in the course of the analyses are scrutinized in the light of recent direct numerical simulations of turbulent flames and the applications to the flames of laser diagnostics. Remaining problem areas include a sufficiently general combination of strain rate and flame curvature PDFs to give a single PDF of flame stretch rate, the nature of flame quenching under positive and negative stretch rates, flame responses to changing stretch rates and the effects of flame instabilities.     

Turbulent scalar fluxes in H2-air premixed flames at low and high Karlovitz numbers

The statistical behaviour and closure of sub-grid scalar fluxes in the context of turbulent premixed combustion have been assessed based on an a priori analysis of a detailed chemistry Direct Numerical Simulation (DNS) database consisting of three hydrogen-air flames spanning the corrugated flamelets (CF), thin reaction zones (TRZ) and broken reaction zones (BRZ) regimes of premixed turbulent combustion. The sub-grid scalar fluxes have been extracted by explicit filtering of DNS data. It has been found that the conventional gradient hypothesis model is not appropriate for the closure of sub-grid scalar flux for any scalar in the context of a multispecies system. However, the predictions of the conventional gradient hypothesis exhibit a greater level of qualitative agreement with DNS data for the flame representing the BRZ regime. The aforementioned behaviour has been analysed in terms of the properties of the invariants of the anisotropy tensor in the Lumley triangle. The flames in the CF and TRZ regimes are characterised by a pronounced two-dimensional anisotropy due to strong heat release whereas a three-dimensional and more isotropic behaviour is observed for the flame located in the BRZ regime. Two sub-grid scalar flux models which are capable of predicting counter-gradient transport have been considered for a priori DNS assessment of multispecies systems and have been analysed in terms of both qualitative and quantitative agreements. By combining the latter two sub-grid scalar flux closures, a new modelling strategy is suggested where one model is responsible for properly predicting the conditional mean accurately and the other model is responsible for the correlations between model and unclosed term. Detailed physical explanations for the observed behaviour and an assessment of existing modelling assumptions have been provided. Finally, the classical Bray–Moss–Libby theory for the scalar flux closure has been extended to address multispecies transport in the context of large eddy simulations.