Next Article in Journal
Molecular Detection of Cordycepin-Induced HeLa Cell Apoptosis with Surface-Enhanced Raman Spectroscopy
Next Article in Special Issue
A Steam Ejector Refrigeration System Powered by Engine Combustion Waste Heat: Part 1. Characterization of the Internal Flow Structure
Previous Article in Journal
An Improved Autoencoder and Partial Least Squares Regression-Based Extreme Learning Machine Model for Pump Turbine Characteristics
Previous Article in Special Issue
Numerical Investigations of Micro-Scale Diffusion Combustion: A Brief Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Eddy-Generation Mechanism on the Characteristic of On-Source Fire Whirl

1
School of Mechanical and Manufacturing Engineering, University of New South Wales, Sydney, NSW 2052, Australia
2
Australian Nuclear Science and Technology Organisation (ANSTO), Locked Bag 2001, Kirrawee DC, NSW 2232, Australia
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(19), 3989; https://doi.org/10.3390/app9193989
Submission received: 9 August 2019 / Revised: 6 September 2019 / Accepted: 13 September 2019 / Published: 24 September 2019
(This article belongs to the Special Issue Progress in Combustion Diagnostics, Science and Technology)

Abstract

:
This paper numerically examines the characterisation of fire whirl formulated under various entrainment conditions in an enclosed configuration. The numerical framework, integrating large eddy simulation and detailed chemistry, is constructed to assess the whirling flame behaviours. The proposed model constraints the convoluted coupling effects, e.g., the interrelation between combustion, flow dynamics and radiative feedback, thus focuses on assessing the impact on flame structure and flow behaviour solely attribute to the eddy-generation mechanisms. The baseline model is validated well against the experimental data. The data of the comparison case, with the introduction of additional flow channelling slit, is subsequently generated for comparison. The result suggests that, with the intensified circulation, the generated fire whirl increased by 9.42 % in peak flame temperature, 84.38 % in visible flame height, 6.81 % in axial velocity, and 46.14 % in velocity dominant region. The fire whirl core radius of the comparison case was well constrained within all monitored heights, whereas that of the baseline tended to disperse at 0.5   m height-above-burner. This study demonstrates that amplified eddy generation via the additional flow channelling slit enhances the mixing of all reactant species and intensifies the combustion process, resulting in an elongated and converging whirling core of the reacting flow.

1. Introduction

Fire whirl is a unique swirling diffusion flame that often occurs in urban and wildland fires with disastrous effects [1]. The whirling combustion flow limits the radial flame dispersion and stretches the hot plume to rise upon an elevated vertical angular path, and as a result significantly increases the burning rate, flame height, and flame temperature and radiation flux to the surroundings [2]. Owing to the unique features associated with this particular combustion dynamic behaviour, the formation of the fire whirl in a fire scenario can be devastated. For instance, the intense heat energy contained within the elongated fire whirl core combined with the radiative heat transfer to the surrounding, promote the fire propagation to the neighbouring points. The swirling motion associated with the fire whirl could also intensify the combustion process by enhancing the mixing of the combustible gas mixture with the oxidiser, while at the same time making the unsteady movement of the hot plume region even more spatially unpredictable [3]. The presence of such phenomena associated with fire whirl has been observed and reported in many notorious fire cases [4,5], and it has been identified as the core factor in intensifying the combustion process and making the scenario unmanageable and untenable. For this very reason, fire whirls, in terms of their formation mechanisms, characteristic features, and physical behaviours, are of great interest to both the research and industrial sectors, and hence have been identified as the topic of interest for this study.
Previous experimental and theoretical studies have characterised and identified the three most critical criteria for the formation of a fire whirl, namely, a thermally driven fluid sink, a surface drag force to create a radial boundary layer, which facilitates air entrainment to the generated vortex column, and an eddy-generation mechanism [6,7,8]. In a typical buoyancy-driven diffusion flame, the flame can be considered as the fluid sink and the generated hot plume naturally drives horizontal flow radically towards the vortex column [9,10], which fulfils the first two requirements. Therefore, the most critical element in fire whirl formation is the presence of a vortex generation mechanism. The generation of vortex could occur in a wildland fire situation, induced, for example, by topological obstacles, leeward slope, unpredicted weather conditions, etc. [11,12,13]. The formation of such circulation can also be observed in a compartment fire scenario [14,15,16], particular with the current trend of needing to construct high-density, high-rise buildings with complex geometric configurations. A typical example of a high-risk region for fire whirl formation is an enclosed space with entrainment from the flow channelling slit to introduce the circulation, for example, atrium, lift pit, spiral staircase etc., which are commonly featured in modern constructions. Despite the recently developed bio-based flame-retardant materials that effectively restrain the fire from spreading [17,18,19,20], the occurrence of the spinning flaming flow in skyscraper fires, such as the recent incidents in Beijing Television Cultural Center, Plasco Building, Grenfell Tower, Marina Torch, FR Tower, etc., could aggressively promote an alarming fire propagation within the confined space and cause catastrophic damages [21]. This poses a severe risk to the occupants and makes the fire suppression process extremely challenging [22,23,24]. Therefore, the characterisation of fire whirl formulated for enclosed configurations was targeted as the investigation case.
Fire whirl formulated in enclosed configurations has been intensively studied, both experimentally and numerically [7,25,26,27,28,29,30,31]. The literature has revealed an entangled coupling between the combustion process and the dynamic flow, i.e., the temperature, buoyancy, vortex, the combustion reaction, etc., are interrelated and interact during the formation of the fire whirl. Such interrelation between the involved physical and chemical aspects makes the understanding of the not-well-understood fundamentals of the swirling reacting flow more unfathomable. For instance, the existence of the whirling flame within the enclosure alters the heat feedback from the flame sheet to the fuel surface, and results in an unevenly distributed heating profile over the fuel surface. In a pool fire configuration, such ununiformed heat feedback could affect the burning rate as it influences the evaporation of combustible gas mixture converted from liquid-based parent fuel, and hence varies the ensuing combustion process and heat release rate [32,33,34,35]. Similarly, the residence of the whirling flame escalates the circulation within the enclosure, enhances the mixing of the oxidiser and reactant, and promotes a more completed combustion process. Theoretically speaking, the formation of soot species, the fine particles resulting from incomplete combustion, and the coupled radiative heat transfer to and from the surrounding, should be suppressed, when compared with non-swirling freestanding flame [36,37,38,39]. Nevertheless, the aggravation of the combustion process increases the flame temperature, which often excesses the threshold of soot nucleation, therefore aid the inception of in-flame generated soot species [40,41]. The excess generated soot species could arguably intensify the radiative heat transfer to the fuel surface and further boost the heat release rate, which contradicts with the previous speculation. In addition to the intricate physical and chemical coupling process, the uncertainty caused by the numerical modelling configuration could make the characterisation of fire whirl more cumbersome. For example, the combustion region of the fire whirl is likely to swirl around the domain which may involve tilting, stretching, converging etc. in a highly unpredicted manner. If the numerical domain is discretised non-uniformly, the variation in spatial resolution of the discretised region may significantly affect the evaluated result, particularly for the prediction of the sensitive flame temperature [42].
The ambiguous correlation, involving various aspects, has nevertheless not been decoupled or restrained in many of the previous studies. Strategies, e.g., heavy sooty flame, pool-based flame generation mechanism where the burning rate depends on the evaporation rate, locally refined meshing algorithm with a denser grid at the domain centre, where a freestanding fire is expected to occur and coarse mesh at the fringes, etc., have commonly been adopted in previous studies, both experimental and numerical. Those implementations make quantitative analysis of the fire whirl behaviour unattainable. It is, therefore, critical to construct a numerical framework with the isolation and consideration of the preceding coupling parties, and focus on establishing the correlation between the eddy-generation mechanism and the characteristic of the formulated fire whirl.
In light of the abovementioned gap, the current study aims to characterise fire whirl formulated under various entrainment conditions via a numerical approach. The proposed numerical model, elaborated in more detail in the next section, is deliberately designed to isolate and restrain the coupling effect between some primary physical, chemical and numerical factors. The proposed model adopts soot-free and clean-burning alcohol-based fuel with a set injection speed, and aims to eliminate the soot radiative heat transfer to and from the surrounding, as well as the potential variation in the burning rate. In addition, the numerical domain of the entire enclosure, where both combustion and dynamic flow motion are expected to occur, is discretised with uniformed grid size to minimise the impact of the spatial resolution of computational control volume on the numerical result. The baseline case with one flow channelling slit is built and validated against experimental data. A modified case with the introduction of an additional entrainment slit is also constructed. The results generated by the two are compared and assessed.

2. Numerical Details

A numerical domain of the baseline model replicating the test rig of the previous experimental study [43], as well as that with an additional flow channelling slit, was constructed accordingly. The geometric features of the two models are shown in Figure 1. The fuel pan is converted from circular to square configuration with the same cross-section area as in the experimental setup to achieve a fully structured mesh, as this significantly improves the numerical accuracy and computational efficiency.
Despite the inevitable “flame wander” phenomenon, the domain is discretised to about 800 , 000 elements based on a uniform division algorithm, to limit the numerical uncertainty associated with the spatial resolution of the discretised control volume. A doubling of the number of elements to 1 , 600 , 000 resulted in only a 5 % difference in the centreline temperature, thus achieving a mesh-independent solution.
The fluid flow and heat transfer within the compartment is described through the conservation equations of continuity, Navier-Stokes, and scalar quantities:
ρ ¯ t + x i ρ ¯ U ˜ i = 0
ρ ¯ U i ˜ + U ˜ i U ˜ j + ρ ¯ u i   u j ˜ = p ¯ x i + ( ν + ν t ) 2 ρ ¯ U ˜ i x j x j
ρ ¯ Φ ˜ α + ρ ¯ Φ ˜ α U ˜ i + ρ ¯ ϕ α u i ˜ = ( Γ + Γ t ) 2 ρ ¯ Φ ˜ α x i x i + S ( Φ ) ˜
where ρ ¯ and p ¯ are the mean density and pressure, ν and ν t are the kinematic and turbulent viscosity, U ˜ i and u i are the Favre-averaged and fluctuation value of the mean fluid velocity, Φ ˜ α and ϕ α are the mean and fluctuation value of the α th scalar, Γ and Γ t are molecular and turbulent diffusivity, and S ( Φ ) ˜ is the Favre-averaged chemical reaction source term.
Large eddy simulation (LES) with Wall-Adapting Local Eddy-viscosity (WALE) function, which accounts for both the strain and rotation rate of the turbulent structure is employed as the turbulence model to describe the turbulent reaction flow behaviour. This turbulence model has been extensively validated and proved to be a valid approach for resolving various wall-bounded turbulence flow applications [8,39,44,45,46,47,48].
The chemical reaction source term in the transport equations of the involved reacting scalars is evaluated via the strained laminar flamelet approach. The involved chemical phenomenon is described as the pre-assumed probability density function (pre-PDF) of two quantities, namely the mixture fraction ( f ) and the scalar dissipation ( χ ). Essentially, the mixture fraction governs the proportion of the fuel mixture in each control volume element in the discretised domain, and the first two moments of the mixture fraction, namely the mean mixture fraction, f ¯ , and the mixture fraction variance, g ¯ , are evaluated by resolving their conservation equations, respectively. Scalar dissipation is introduced to describe the strain and extinction of the flame, in which larger values of this quantity depict its departure from chemical equilibrium [49]. The scalar dissipation across the flamelet is determined based on the mean scalar dissipation of the flamelets embedded in turbulent flames at the position where f = f s t , as well as the ratio of the function of the mixture fraction and density at the point of interest and at the position where f is stoichiometric. A detailed GRI-MECH 3.0 chemical reaction mechanism, with the inclusion of 325 reaction steps and 53 chemical species [50], is implemented to establish the flamelet library of the strained laminar flamelet model considering parental fuel. This approach has been validated in previous numerical studies for modelling turbulence chemistry interaction, and has been proved to provide a reasonable result with moderate computational consumption [44,51]. It should be highlighted that methanol ( CH 3 OH ), an alcohol-based fuel, was deliberately selected as the parent fuel to minimise the coupling between radiative energy feedback and the combustion process, due to its feature of clean and soot-free burning behaviour [41]. For the same reason, the concentration of the key intermediate species of soot formation, acetylene ( C 2 H 2 ), in the gas mixture produced by this very parent fuel, is proved to be negligible. The most commonly adopted soot models formulated on C 2 H 2 precursor-based inception and surface growth mechanism are not applicable in this study [52,53,54,55]. As a result, a primitive and computational lightweight two-step soot model is integrated, for concept verification purpose only.
With respect to the boundary conditions, a set flow rate of 0.0216   ms 1 in the direction normal to the fuel surface is applied on the fuel pan. The applied injection velocity is determined on the basis of the cross-sectional area of the fire pan, the density of the parent fuel at reference temperature, and the heat of combustion of the parent fuel in order to match the targeted heat release rate reported in the experimental study. The constant injection rate ensures an evenly distributed burning profile, regardless of the intensity of flow circulation and energy feedback. The fuel surface is set at an elevated height according to the experimental setup. The top and the periphery of the domain is set as the opening to allow naturally convected air entrainment in and out of the system. The base of the domain is set as a non-slip adiabatic wall, across which no heat or matter is allowed to pass. The simulation is initiated with a standard temperature and pressure replicating the ambient conditions for the combustion process to proceed. For convenience, the case with one side entrainment slit is referred to as the 01 slit case, case 1 , or the baseline case/model, and that with two side entrainment slits is referred as the 02 slit case, case 2 , or the comparison case/model, in the following sections.

3. Results

The numerical result generated via the baseline case is firstly validated against the experimental data [43], and the characterisation of the fire whirl evaluated via the two cases is subsequently assessed and discussed in the following section.

3.1. Validation

An accurate description of the fire whirl’s inner structure is intimately tied with a proper description of the temperature and the velocity field of the swirling reacting flow. For this reason, the temperature and velocity profile generated by the baseline model is used to compare with the that reported in the literature for validation purpose.
Due to the constraints related to the physical set up and experimental design, limited information is provided in prior studies. For clarity, herein, temperature measurement at domain centreline at various HABs is selected for comparison. Furthermore, the incoming y-velocity profile through the air gap is also validated against a prior numerical investigation of the similar experimental setup [8]. The result of the validation procedure is presented in Figure 2 and Figure 3.
It can be seen from Figure 2 that the centreline temperature peaks at about HAB 0.15   m and sharply decreases to ambient as the combustion process develops towards downstream. Similarly, high incoming velocity at the air gap is observed in the low-intermediate region, and its quantity gradually drops with the increase of flame height. Despite the fact that the flame temperature in the intermittent region is slightly overestimated, potentially attributed to the absence of radiation correction, as well as the inevitable flame flickering, the comparison shows good agreement between the numerical simulation and experimental data. As presented in Figure 3, a similar agreement could be reached in the comparison of incoming y-velocity through the air gap with the prior study.
In addition, the soot volume fraction evaluated from both cases is proved to be negligible, i.e., on an order of magnitude of 21 and 18   in the baseline and comparison cases, respectively. It was further ascertained that with the adoption of the alcohol-based parent fuel, the in-flame soot formation and the ensuring radiation feedback can be safely neglected in this study.
With the numerical model of the baseline case validated, the model of the two slits cases is constructed accordingly, by implementing the identical numerical setup on the geometry with the introduction of an additional flow channelling slit. The result generated by the two slits model is therefore assumed to be justified. Temperature and velocity profiles, as mentioned, are the primary parameters for characterising the swirling reacting flow involved in the fire whirl. The assorted temperature and velocity data generated numerically by the baseline model are presented alongside those of the comparison model, to assess the effect of the geometrical configuration on the formation of fire whirl.

3.2. Temperature

The thermal component of the flame is undeniably the most critical aspect to be investigated in fire modelling. It is more critical in this study of fire whirl, as the temperature and velocity are interrelated, and such coupling is essential for the understanding of the ensuing formation and sustenance of the swirling combustion phenomenon. The qualitative analysis of the fire whirl therefore starts by assessing the temperature profile of the formulated swirling reacting flow field.
The temperature iso-surface of the flame generated by the two cases, presented in Figure 4, provides a glance of the geometric characteristics of the fire whirl formulated under different entrainment boundary conditions. The figure shows that, with the side flow channelling slit(s) acting as the eddy-generation mechanism, a highly twisted flow field with relatively high temperature can clearly be observed in both geometric configurations. This temperature profile exhibits a helical form, implying the persistence of strong swirl motion within the core structure of both flames. To be more specific, the fire whirl generated from the baseline model is exhibited in the format of a tortuous and highly twisted string. The hot plume accumulates and exits up near exit region in the vertical orientation. Whereas the fire whirl of that generated via the two-slit configuration resembles a convoluted cylinder centred on the lower-mediate region of the computational domain.
The temperature contour of a cross-sectional plane of the computational domains of the two cases is illustrated in Figure 5a,b. The contour plot shows that the maximum flame temperature of the case with the two-slit configuration is higher, i.e., increased by about 9.42 % , when compared with the baseline. Similar to the temperature iso-surface presented above, the cloud of the hot plume, accumulated at the downstream region observed in the baseline case, is cleared out when the additional slit is introduced.
The visible flame height of the two cases determined based on applying the cut the cut-off visible flame tip temperature of 600   ° K is presented in Figure 5c,d. From the comparison, it is evident that the visible flame height for two-slit configuration grows significantly compared with that of the baseline model, i.e., it is increased by 84.38   % from 0.32   m to 0.59   m .
The extension of the flame height of the comparison model is further verified by comparing the spatial temperature distribution, along with the lateral direction, as presented in Figure 6. Theoretically, the peak flame temperature of a typical diffusion flame should occur at the reaction layers, i.e., the intersection where the fuel and the oxidiser encounter each other and react due to diffusion. The reaction layer in the continuous and intermittent region of a diffusion flame is typically away from the flame centreline, and as a result, the spatial distribution of the flame temperature should resemble a hump-like profile. The hump-like profile tends to shift to a Gaussian distribution, with its value maximised at the centreline, as the flame flows further to the downstream, i.e., enters the plume region.
The transition from a hump-like profile to a Gaussian distribution can be observed in the baseline model, as the HAB increased from 0.3   m to 0.5   m . This implies that the domain from HAB 0.5   m and beyond of this very flame can be identified as the plume region. On the other hand, the temperature profile generated by the 02 slit case features a hump-like profile for all monitoring HABs, means that the flaming structure up to HAB 0.5   m of this fire whirl can still be considered to be in the continuous and intermittent region. This therefore indicates that the flame height of the swirling reacting flow triggered by the intensified circulation is significantly higher compared with the baseline model.
It is also worthwhile assessing the temperature core radius b T , defined as the radial location where the excess temperature declines to half of the maximum recorded value, at the monitoring HABs of the two cases. Pioneer studies of the fire whirl suggested that the flame pulsation behaviour is suppressed in the continuous flame region due to the cyclostrophic balance effect, resulting in a near-constant temperature core radius. As the flame enters the plume region, the temperature core radius increases proportionally with the flame height [33].
As reported in Table 1, the temperature core radius of the baseline model expands significantly; for example, it is increased by 153.33 % from 0.0660   m at HAB 0.3   m to 0.1672   m   at HAB 0.5   m , while that of the 02 slit case is reported to be nearly identical at all monitoring HABs. The results suggest an elongated continuous flame region of the 02 slit model when compared with that of the baseline.
The preceding comparison of the spatial temperature distribution demonstrated that, as a buoyant driven diffusion flame, the intensification of the eddy generation promotes the air entrainment from the surroundings, enhances the mixing of fuel and oxidizer. This facilitates the combustion process, hence resulting in the formation of the flame field with higher peak temperature. The additional entrainment induced via the extra slit also amplifies the generation of vorticity, aggravates the swirl motion, and escalates the formulation of fire whirl, which leads to intensified combustion reactions and elongated combustion region. In addition, for this very reason, the 02 slit model promotes natural ventilation at the near exit region, and as a result, the hot plume is entrained out of the computational domain more effectively compared with 01 slit configuration.

3.3. Velocity

To characterise the fire whirl formulated under various entrainment conditions, it is imperative to evaluate and assess the related velocity field. In essence, the inner structure of an on-source fire whirl can be approximated as a 3-dimensional flow. The following vector plot, as well as the peak value of the velocity of U, V and W at all monitoring HABs, presented in Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12, illustrates the key features of the flow behaviour of the swirling flame structure of both cases. This provides compelling evidence that, with the introduction of the slit to trigger eddy generation, the presence of the highly swirled reacting flow field, i.e., the formation of the fire whirl, can be confirmed. Such whirling motion within the flame structure can be observed at all HABs of interest for both geometric configurations.
The peak value of the associated velocity at three monitoring HABs is plotted in Figure 8, Figure 10 and Figure 12. The comparison indicates that the peak values of U, V, and W velocity of the baseline model tend to firstly increase in the region between HAB 0.1   m and HAB 0.3   m , and then decrease sharply as the combustion spreads to downstream HABs. Meanwhile, for the 02 slit configuration, the magnitude of the drop of the assorted velocities from HAB 0.3   m to HAB 0.5   m is observed to be moderated compared with the 01 slit configuration. For example, the results report a 1.06 % , + 36.65 % and 3.83 % variation in U, V, and W velocity between HAB 0.3   m and HAB 0.5   m of the 02 slit case, compared with that of 9.25 % , 26.34 % and −28.34% of the baseline.
The dampened reduction in U and W velocity and the significant increase in V velocity from HAB 0.3   m to HAB 0.5   m demonstrate that, with the introduction of additional air entrainment slit, the formulated flame is still in the developing phase, i.e., in the continuous and intermittent regions at HAB of 0.5   m . Meanwhile, the flame yield for the baseline model is well-developed, and reached the plume region at the same elevated height; hence, it possesses a shorter flame length, which is in good agreement with the speculation deduced from the preceding analysis of the temperature profile.

3.3.1. Axial Velocity

The axial velocity contour of a cross-sectional plane generated by the two cases is presented in Figure 13 for comparison. The contour plot depicts that the maximum axial velocity of the case with 02 slit configuration is higher, i.e., increased by 6.81 % , when compared with that of the baseline model. The region in which axial velocity is dominant was also extended with the introduction of additional air entrainment slit. The axial velocity dominant region, determined as the vertical location at which the excess axial velocity decreases to half of the maximum recorded value, increased by 46.14 % to 0.5739   m , compared with 0.3927   m for the baseline model.
The promotion of axial velocity in the 02 slit case is further justified by assessing the spatial distribution of the axial velocity plotted against the lateral direction at the monitoring HABs, presented in Figure 14. The literature reveals that the axial centreline velocity accelerates upwards to maximum in the continuous flame region. Due to the velocity fluctuation, the axial centreline velocity decreases sharply as the HAB increases to intermittent regions, and further decelerates to ambient conditions as it enters the plume region [56,57]. Such a transition can be clearly seen in the velocity profile yield from the baseline model, indicating that the continuous and intermittent region of the fire whirl generated in the baseline case should be well below HAB 0.5   m . Meanwhile, the peak axial centreline velocity reported in the comparison case is presented in an incremental pattern as the flame developing in the downstream direction, suggesting an extended continuous region when compared with the baseline case.
In addition, experimental data of PIV measurements confirmed that within the continuous flame region, the axial velocity reaches its maximum value in the vicinity of the reaction layer, i.e., away from the flame centreline, and its profile resembles a hump-like distribution along the radial direction. With increasing height from the intermittent flame region and beyond, the location of maximum excess temperature shifts toward the whirl’s centreline and forms a Gaussian profile [27]. The transformation of the axial velocity profile from hump-like to Gaussian form, evidently observed in the baseline model as the HAB increased from 0.3   m to 0.5   m , however, had not yet occurred in the 02 slit model. The hump-like profile at all monitoring HABs up to 0.5   m in the comparison case again confirms the extension of the continuous flame region, i.e., to and beyond HAB 0.5   m , when compared with single slit configuration.

3.3.2. Radial Velocity and Tangential Velocity

The profile of U, V and radial velocity, plotted as a function of the radial position at three monitoring HABs, is presented in Figure 15. The transformation of the profile of U, V and radial velocity from hump-like to a plateau with the increase of flame height was observed in the baseline model, whereas that hump-like velocity profile is inherently persistent in all HABs in the 02 slit case. Furthermore, the radial location at which the velocity of interest peaks shifts away from the centreline axis with an increase in the flame height in the baseline case, whereas this is nearly identical at all HABs in the comparison case. It has been described in prior studies that the maximum radial velocity occurs adjacent to the surface and reverses direction close to the upper parts of the boundary layer, indicating the presence of the circulation zone [26]. The departure of the boundary layer of the circulation zone away from the centreline, observed in the baseline case, implies the dissipation of the fire whirl core region. In other words, the flame departs from the continuous and intermittent region and enters the plume zone at the elevated height of 0.5   m . The near-constant radial position of the peak radial velocity observed in the 02 slit case depicts a well-maintained fire whirl core structure, i.e., the flame at high elevations can still be considered to be within the continuous region. The combustion region of the comparison case, therefore, is demonstrated to be well-extended compared with that of the baseline.
Circulation is indeed the unique feature distinguishing fire whirls from non-swirling fires. The presence of the circulation disturbs the flow field and alters the ensuing combustion process, as a result significantly affecting assorted properties, such as flame temperature, combustion dominant region, as well as the velocity profile. The circulation is characterized by assessing the azimuthal velocity, or tangential velocity, in this study. Herein, the tangential velocity against the radial position at the selected HABs of the cases of interest is plotted in Figure 16.
Similar to the profile of radial velocity, The humped-like profile turns to flat and the radial location where the tangential velocity peaks departs away from domain centreline with the increase of flame height in the baseline model, while those of the comparison case remains relatively unchanged at all monitored heights.
To further elaborate the transition, the mean whirl vortex core radius at the elevated heights of the two cases are evaluated and listed in Table 2. The mean whirl core radius is determined based on the azimuthal velocity profile, i.e., the radial location at which the tangential velocity is at its maximum, and beyond which the circulation becomes nearly constant. The results reveal that the mean whirl core radius of the baseline case departs significantly from the centreline axis, i.e., it increases by 136.91 % from 2.33   cm at HAB 0.3   m to 5.52   cm at HAB 0.5   m , while that reported for the 02 slit case is nearly identical at all monitoring HABs. Previous experimental study examining the correlation between the fire whirl radius and the stage of combustion suggests a reasonably constant whirl core radius in the continuous flame region, which increases with flame height in the intermittent and plume regions [33]. The results presented in Table 2, therefore, again verify the abovementioned speculation that the combustion region is well maintained and is extended by means of the introduction of the additional flow channelling slit.
In addition, by comparing the results presented in Table 2 with those in Table 1, it can be concluded that the radius determined based on the temperature is generally larger than that based on circulation, which again is in good agreement with the experimental observations [1].
Collectively, the preceding assessment of the axial, radial and tangential velocity demonstrates that, with the additional vorticity induced by the entrainment via the additional slit, the swirling motion within the domain is intensified, and the circulation of the dominant region is elongated. The intensified circulation enhances the mixing of the oxidiser with combustible gas mixtures and facilitates the ensuing combustion process. As a result, a quasi-steady, on-source fire whirl with a unique flow core structure that extended to a higher elevated HAB was formulated and well maintained

4. Conclusions

A numerical simulation of a unique combustion process occurring in an enclosed configuration with one or two side entrainment slit(s) was performed to assess the character of a fire whirl formulated by means of various eddy-generation mechanisms. The numerical framework of the proposed model consists of the WALE large eddy simulation turbulence model coupled with reduced detailed chemistry to describe the intricate whirling flaming flow behaviour. The numerical model adopts a uniformed structured discretisation scheme, a fixed fuel injection rate, and a soot-free alcohol-based parent fuel, to isolate and constrain the complicated coupling effect between the combustion process and the flow dynamics.
The numerical simulation of the baseline model consists of one flow channelling slit, replicating the test rig layout in the literature. This was firstly validated against experimental data and achieved satisfactory agreement. A modified model built with the introduction of an additional slit was subsequently constructed, and the results generated were compared with the baseline for assessment. The analysis concluded that, with the intensified eddy-generation mechanism,
  • The peak flame temperature of the formulated fire whirl increased by 9.42 % , from 1380   K to 1510   K ;
  • The visible flame height increased by 84.38 % , from 0.32   m to 0.59   m ;
  • The axial velocity and the axial velocity dominant region increased by 6.81 % and by 46.14 % , respectively;
  • The flame structure was well-maintained as the temperature core radius remained relatively unchanged at all monitoring HABs, compared with the increase by 153.33 % in the baseline case;
  • The circulation core structure was well-sustained, and the mean vortex core radius was observed to be nearly identical at the elevated height of up to 0.5   m , whereas there was an increase of 136.91 % in the baseline case.
In summary, the results generated from this numerical study demonstrate that the addition of an introduced slit intensifies the circulation, enhances the mixing of the oxidiser with the reactants, and hence facilitates the ensuing combustion process. As a result, the numerical model with stronger eddy-generation mechanisms, i.e., featuring two flow channelling slits, formulates and sustains a fire whirl possessing an elongated, constrained and swirling combustion region, compared with those generated from the single-slit configuration. In addition, this numerical study showcases the capability of the proposed numerical framework to capture various features of the whirling reacting flow, which are often difficult to acquire and assess by experimental means. The model also provides visualised information on the combustion process and flow field, thus delivering theoretical insight. The numerical model presented in this study, is therefore proved to be a feasible tool for assisting the understanding of the fundamentals of the formation and evolution of this intriguing combustion phenomenon.

Author Contributions

Conceptualization, A.C.Y.Y., Q.N.C., S.K. and G.H.Y.; Data curation, C.W. and H.L.Y.; Formal analysis, C.W., T.B.Y.C., R.C. and H.L.Y.; Funding acquisition, A.C.Y.Y. and G.H.Y.; Methodology, A.C.Y.Y., Q.N.C., S.K. and G.H.Y.; Resources, S.K.; Supervision, A.C.Y.Y., Q.N.C., S.K. and G.H.Y.; Validation, Q.C.; Visualization, Q.C.; Writing—original draft, C.W.; Writing—review & editing, T.B.Y.C. and R.C.

Funding

This research was funded by Australian Research Council (ARC Industrial Transformation Training Centre IC170100032), and Wuhan Shuanglian-Xingxin Machinery & Equipment Co. Ltd., China.

Acknowledgments

The authors wish to acknowledge the financial support of Australian Research Council (ARC Industrial Transformation Training Centre IC170100032) and the Australian Government Research Training Program Scholarship. It is also sponsored by Wuhan Shuanglian-Xingxin Machinery & Equipment Co. Ltd., China. The authors deeply appreciate all financial and technical supports.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tohidi, A.; Gollner, M.J.; Xiao, H. Fire Whirls. Annu. Rev. Fluid Mech. 2018, 50, 187–213. [Google Scholar] [CrossRef] [Green Version]
  2. Wang, P.; Liu, N.; Bai, Y.; Zhang, L.; Satoh, K.; Liu, X. An experimental study on thermal radiation of fire whirl. Int. J. Wildl. Fire 2017, 26, 693. [Google Scholar] [CrossRef]
  3. Medwell, P.R.; Chan, Q.N.; Kalt, P.A.M.; Alwahabi, Z.T.; Dally, B.B.; Nathan, G.J. Instantaneous Temperature Imaging of Diffusion Flames Using Two-Line Atomic Fluorescence. Appl. Spectrosc. 2010, 64, 173–176. [Google Scholar] [CrossRef] [PubMed]
  4. Soma, S.; Saito, K. Reconstruction of fire whirls using scale models. Combust. Flame 1991, 86, 269–284. [Google Scholar] [CrossRef]
  5. Zhou, K.; Liu, N.; Yin, P.; Yuan, X.; Jiang, J. Fire Whirl due to Interaction between Line Fire and Cross Wind. Fire Saf. Sci. 2014, 11, 1420–1429. [Google Scholar] [CrossRef]
  6. Byram, G.; Martin, R. The Modeling of Fire Whirlwinds. For. Sci. 1970, 16, 386–399. [Google Scholar]
  7. Lei, J.; Liu, N.; Tu, R. Flame height of turbulent fire whirls: A model study by concept of turbulence suppression. Proc. Combust. Inst. 2017, 36, 3131–3138. [Google Scholar] [CrossRef]
  8. Yuen, A.C.Y.; Yeoh, G.H.; Cheung, S.C.P.; Chan, Q.N.; Chen, T.B.Y.; Yang, W.; Lu, H. Numerical study of the development and angular speed of a small-scale fire whirl. J. Comput. Sci. 2018, 27, 21–34. [Google Scholar] [CrossRef]
  9. Church, C.R.; Snow, J.T.; Dessens, J. Intense Atmospheric Vortices Associated with a 1000 MW Fire. Bull. Am. Meteorol. Soc. 1980, 61, 682–694. [Google Scholar] [CrossRef] [Green Version]
  10. Forthofer, J.M.; Goodrick, S.L. Review of Vortices in Wildland Fire. J. Combust. 2011, 2011, 1–14. [Google Scholar] [CrossRef]
  11. Liu, N.; Liu, Q.; Deng, Z.; Kohyu, S.; Zhu, J. Burn-out time data analysis on interaction effects among multiple fires in fire arrays. Proc. Combust. Inst. 2007, 31, 2589–2597. [Google Scholar] [CrossRef]
  12. Emori, R.I.; Saito, K. Model experiment of hazardous forest fire whirl. Fire Technol. 1982, 18, 319–327. [Google Scholar] [CrossRef]
  13. Kuwana, K.; Morishita, S.; Dobashi, R.; Chuah, K.H.; Saito, K. The burning rate’s effect on the flame length of weak fire whirls. Proc. Combust. Inst. 2011, 33, 2425–2432. [Google Scholar] [CrossRef]
  14. Chuah, K.H.; Kuwana, K.; Saito, K.; Williams, F.A. Inclined fire whirls. Proc. Combust. Inst. 2011, 33, 2417–2424. [Google Scholar] [CrossRef]
  15. Lei, J.; Liu, N.; Zhang, L.; Chen, H.; Shu, L.; Chen, P.; Deng, Z.; Zhu, J.; Satoh, K.; De Ris, J.L. Experimental research on combustion dynamics of medium-scale fire whirl. Proc. Combust. Inst. 2011, 33, 2407–2415. [Google Scholar] [CrossRef]
  16. Zhou, K.; Liu, N.; Lozano, J.S.; Shan, Y.; Yao, B.; Satoh, K. Effect of flow circulation on combustion dynamics of fire whirl. Proc. Combust. Inst. 2013, 34, 2617–2624. [Google Scholar] [CrossRef]
  17. Lin, B.; Yuen, A.C.Y.; Li, A.; Zhang, Y.; Chen, T.B.Y.; Yu, B.; Lee, E.W.M.; Peng, S.; Yang, W.; Lu, H.-D.; et al. MXene/chitosan nanocoating for flexible polyurethane foam towards remarkable fire hazards reductions. J. Hazard. Mater. 2020, 381, 120952. [Google Scholar] [CrossRef] [PubMed]
  18. Yang, W.-J.; Yuen, A.C.Y.; Li, A.; Lin, B.; Chen, T.B.Y.; Yang, W.; Lu, H.-D.; Yeoh, G.H. Recent progress in bio-based aerogel absorbents for oil/water separation. Cellulose 2019, 26, 6449–6476. [Google Scholar] [CrossRef]
  19. Si, J.-Y.; Tawiah, B.; Sun, W.-L.; Lin, B.; Wang, C.; Yuen, A.C.Y.; Yu, B.; Li, A.; Yang, W.; Lu, H.-D.; et al. Functionalization of MXene Nanosheets for Polystyrene towards High Thermal Stability and Flame Retardant Properties. Polymers 2019, 11, 976. [Google Scholar] [CrossRef]
  20. Yu, B.; Tawiah, B.; Wang, L.-Q.; Yin Yuen, A.C.; Zhang, Z.-C.; Shen, L.-L.; Lin, B.; Fei, B.; Yang, W.; Li, A.; et al. Interface decoration of exfoliated MXene ultra-thin nanosheets for fire and smoke suppressions of thermoplastic polyurethane elastomer. J. Hazard. Mater. 2019, 374, 110–119. [Google Scholar] [CrossRef]
  21. Chen, T.; Yuen, A.; Yeoh, G.; Yang, W.; Chan, Q. Fire Risk Assessment of Combustible Exterior Cladding Using a Collective Numerical Database. Fire 2019, 2, 11. [Google Scholar] [CrossRef]
  22. Li, A.; Yuen, A.C.Y.; Chen, T.B.Y.; Wang, C.; Liu, H.; Cao, R.; Yang, W.; Yeoh, G.H.; Timchenko, V. Computational Study of Wet Steam Flow to Optimize Steam Ejector Efficiency for Potential Fire Suppression Application. Appl. Sci. 2019, 9, 1486. [Google Scholar] [CrossRef]
  23. Yuen, A.C.Y.; Yeoh, G.H.; Alexander, B.; Cook, M. Fire scene investigation of an arson fire incident using computational fluid dynamics based fire simulation. Build. Simul. 2014, 7, 477–487. [Google Scholar] [CrossRef]
  24. Yuen, A.C.Y.; Yeoh, G.H.; Alexander, R.; Cook, M. Fire scene reconstruction of a furnished compartment room in a house fire. Case Stud. Fire Saf. 2014, 1, 29–35. [Google Scholar] [CrossRef] [Green Version]
  25. Dobashi, R.; Okura, T.; Nagaoka, R.; Hayashi, Y.; Mogi, T. Experimental Study on Flame Height and Radiant Heat of Fire Whirls. Fire Technol. 2016, 52, 1069–1080. [Google Scholar] [CrossRef]
  26. Hartl, K.A.; Smits, A.J. Scaling of a small scale burner fire whirl. Combust. Flame 2016, 163, 202–208. [Google Scholar] [CrossRef]
  27. Wang, P.; Liu, N.; Hartl, K.; Smits, A. Measurement of the Flow Field of Fire Whirl. Fire Technol. 2016, 52, 263–272. [Google Scholar] [CrossRef]
  28. Xiao, H.; Gollner, M.J.; Oran, E.S. From fire whirls to blue whirls and combustion with reduced pollution. Proc. Natl. Acad. Sci. USA 2016, 113, 9457–9462. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Muraszew, A.; Fedele, J.B.; Kuby, W.C. The fire whirl phenomenon. Combust. Flame 1979, 34, 29–45. [Google Scholar] [CrossRef]
  30. Chuah, K.H.; Kushida, G. The prediction of flame heights and flame shapes of small fire whirls. Proc. Combust. Inst. 2007, 31, 2599–2606. [Google Scholar] [CrossRef]
  31. Parente, R.M.; Pereira, J.M.C.; Pereira, J.C.F. On the influence of circulation on fire whirl height. Fire Saf. J. 2019, 106, 146–154. [Google Scholar] [CrossRef]
  32. Lei, J.; Liu, N.; Satoh, K. Buoyant pool fires under imposed circulations before the formation of fire whirls. Proc. Combust. Inst. 2015, 35, 2503–2510. [Google Scholar] [CrossRef]
  33. Lei, J.; Liu, N.; Zhang, L.; Satoh, K. Temperature, velocity and air entrainment of fire whirl plume: A comprehensive experimental investigation. Combust. Flame 2015, 162, 745–758. [Google Scholar] [CrossRef]
  34. Sikanen, T.; Hostikka, S. Modeling and simulation of liquid pool fires with in-depth radiation absorption and heat transfer. Fire Saf. J. 2016, 80, 95–109. [Google Scholar] [CrossRef]
  35. Yao, W.; Yin, J.; Hu, X.; Wang, J.; Zhang, H. Numerical modeling of liquid n-heptane pool fires based on heat feedback equilibrium. Procedia Eng. 2013, 62, 377–388. [Google Scholar] [CrossRef]
  36. Wang, C.; Chan, Q.N.; Zhang, R.; Kook, S.; Hawkes, E.R.; Yeoh, G.H.; Medwell, P.R. Automated determination of size and morphology information from soot transmission electron microscope (TEM)-generated images. J. Nanopart. Res. 2016, 18, 127. [Google Scholar] [CrossRef]
  37. Wang, C.; Chan, Q.N.; Kook, S.; Hawkes, E.R.; Lee, J.; Medwell, P.R. External irradiation effect on the growth and evolution of in-flame soot species. Carbon 2016, 102, 161–171. [Google Scholar] [CrossRef] [Green Version]
  38. Yuen, A.C.Y.; Yeoh, G.H.; Timchenko, V.; Cheung, S.C.P.; Chan, Q.N.; Chen, T. On the influences of key modelling constants of large eddy simulations for large-scale compartment fires predictions. Int. J. Comut. Fluid Dyn. 2017, 31, 324–337. [Google Scholar] [CrossRef]
  39. Yuen, A.C.Y.; Yeoh, G.H.; Timchenko, V.; Chen, T.B.Y.; Chan, Q.N.; Wang, C.; Li, D.D. Comparison of detailed soot formation models for sooty and non-sooty flames in an under-ventilated ISO room. Int. J. Heat Mass Transf. 2017, 115, 717–729. [Google Scholar] [CrossRef]
  40. Tree, D.R.; Svensson, K.I. Soot processes in compression ignition engines. Prog. Energy Combust. Sci. 2007, 33, 272–309. [Google Scholar] [CrossRef]
  41. Wang, H. Formation of nascent soot and other condensed-phase materials in flames. Proc. Combust. Inst. 2011, 33, 41–67. [Google Scholar] [CrossRef]
  42. Jones, W.P.; Whitelaw, J.H. Calculation methods for reacting turbulent flows: A review. Combust. Flame 1982, 48, 1–26. [Google Scholar] [CrossRef]
  43. Chow, W.; Han, S.S. Experimental Data on Scale Modeling Studies on Internal Fire Whirls. Int. J. Eng. Perform. Based Fire Codes 2011, 10, 63–74. [Google Scholar]
  44. Chen, T.B.Y.; Yuen, A.C.Y.; Wang, C.; Yeoh, G.H.; Timchenko, V.; Cheung, S.C.P.; Chan, Q.N.; Yang, W. Predicting the fire spread rate of a sloped pine needle board utilizing pyrolysis modelling with detailed gas-phase combustion. Int. J. Heat Mass Transf. 2018, 125, 310–322. [Google Scholar] [CrossRef]
  45. Yuen, A.; Chen, T.; Yang, W.; Wang, C.; Li, A.; Yeoh, G.; Chan, Q.; Chan, M. Natural Ventilated Smoke Control Simulation Case Study Using Different Settings of Smoke Vents and Curtains in a Large Atrium. Fire 2019, 2, 7. [Google Scholar] [CrossRef]
  46. Yuen, A.C.Y.; Yeoh, G.H.; Timchenko, V.; Cheung, S.C.P.; Chen, T. Study of three LES subgrid-scale turbulence models for predictions of heat and mass transfer in large-scale compartment fires. Numer. Heat Transf. Part A Appl. 2016, 69, 1223–1241. [Google Scholar] [CrossRef]
  47. Yuen, A.C.Y. On the Prediction of Combustion Products and Soot Particles in Compartment Fires. Ph.D. Thesis, UNSW Sydney, Kensington, Australia, 2014. [Google Scholar]
  48. Wang, C. External Irradiation Effect on the Evolution of In-flame Soot Species. Ph.D. Thesis, University of New South Wales, Sydney, Australia, 2015. [Google Scholar]
  49. Yuen, A.C.Y.; Yeoh, G.H.; Timchenko, V.; Barber, T. LES and multi-step chemical reaction in compartment fires. Numer. Heat Transf. Part A Appl. 2015, 68, 711–736. [Google Scholar] [CrossRef]
  50. Kee, R.J.; Rupley, F.M.; Miller, J.A.; Coltrin, M.E.; Grcar, J.F.; Meeks, E.; Moffat, H.K.; Lutz, A.E.; Dixon-Lewis, G.; Smooke, M.D.; et al. Chemkin Collection Release 3.6; Reaction Design: San Diego, CA, USA, 2000. [Google Scholar]
  51. Yuen, A.C.Y.; Yeoh, G.H.; Timchenko, V.; Cheung, S.C.P.; Barber, T.J. Importance of detailed chemical kinetics on combustion and soot modelling of ventilated and under-ventilated fires in compartment. Int. J. Heat Mass Transf. 2016, 96, 171–188. [Google Scholar] [CrossRef]
  52. Chen, T.B.Y.; Yuen, A.C.Y.; Yeoh, G.H.; Timchenko, V.; Cheung, S.C.P.; Chan, Q.N.; Yang, W.; Lu, H. Numerical study of fire spread using the level-set method with large eddy simulation incorporating detailed chemical kinetics gas-phase combustion model. J. Comput. Sci. 2018, 24, 8–23. [Google Scholar] [CrossRef]
  53. Wang, C.; Yuen, A.C.Y.; Chan, Q.N.; Chen, T.B.Y.; Yang, W.; Cheung, S.C.-P.; Yeoh, G.H. Sensitivity Analysis of Key Parameters for Population Balance Based Soot Model for Low-Speed Diffusion Flames. Energies 2019, 12, 910. [Google Scholar] [CrossRef]
  54. Mueller, M.E.; Chan, Q.N.; Qamar, N.H.; Dally, B.B.; Pitsch, H.; Alwahabi, Z.T.; Nathan, G.J. Experimental and computational study of soot evolution in a turbulent nonpremixed bluff body ethylene flame. Combust. Flame 2013, 160, 1298–1309. [Google Scholar] [CrossRef]
  55. Wang, C.; Yuen, A.C.Y.; Chan, Q.N.; Chen, T.B.Y.; Yang, W.; Cheung, S.C.P.; Yeoh, G.H. Characterisation of soot particle size distribution through population balance approach and soot diagnostic techniques for a buoyant non-premixed flame. J. Energy Inst. 2019, in press. [Google Scholar] [CrossRef]
  56. Zukoski, E.E.; Kubota, T.; Cetegen, B. Entrainment in fire plumes. Fire Saf. J. 1981, 3, 107–121. [Google Scholar] [CrossRef]
  57. Lei, J.; Liu, N. Flame precession of fire whirls: A further experimental study. Fire Saf. J. 2016, 79, 1–9. [Google Scholar] [CrossRef]
Figure 1. Geometric features of the computational domain.
Figure 1. Geometric features of the computational domain.
Applsci 09 03989 g001
Figure 2. Numerical prediction of temperature at domain centreline, compared with experimental measurements.
Figure 2. Numerical prediction of temperature at domain centreline, compared with experimental measurements.
Applsci 09 03989 g002
Figure 3. Numerical prediction of incoming y-velocity through the air gap, compared with that generated from prior study.
Figure 3. Numerical prediction of incoming y-velocity through the air gap, compared with that generated from prior study.
Applsci 09 03989 g003
Figure 4. Temperature iso-surface of the fire whirl formulated in 01 slit case (a) and the 02 slit case (b).
Figure 4. Temperature iso-surface of the fire whirl formulated in 01 slit case (a) and the 02 slit case (b).
Applsci 09 03989 g004
Figure 5. Temperature contour of a cross-sectional plane of the 01 slit case (a) and the 02 slit case (b), and the visible flame height of the fire whirl formulated by the 01 slit case (c) and the 02 slit case (d).
Figure 5. Temperature contour of a cross-sectional plane of the 01 slit case (a) and the 02 slit case (b), and the visible flame height of the fire whirl formulated by the 01 slit case (c) and the 02 slit case (d).
Applsci 09 03989 g005
Figure 6. Plot of temperature along with the radial direction at three monitoring HABs: 01 slit case (a) and 02 slit case (b).
Figure 6. Plot of temperature along with the radial direction at three monitoring HABs: 01 slit case (a) and 02 slit case (b).
Applsci 09 03989 g006
Figure 7. Plot of vector field of U-velocity of the two cases, at all monitoring HABs.
Figure 7. Plot of vector field of U-velocity of the two cases, at all monitoring HABs.
Applsci 09 03989 g007
Figure 8. Peak value of U-velocity of the two cases, at all monitoring HABs.
Figure 8. Peak value of U-velocity of the two cases, at all monitoring HABs.
Applsci 09 03989 g008
Figure 9. Plot of the vector field of V-velocity of the two cases, at all monitoring HABs.
Figure 9. Plot of the vector field of V-velocity of the two cases, at all monitoring HABs.
Applsci 09 03989 g009
Figure 10. Peak value of V-velocity of the two cases, at all monitoring HABs.
Figure 10. Peak value of V-velocity of the two cases, at all monitoring HABs.
Applsci 09 03989 g010
Figure 11. Plot of vector field of W-velocity of the two cases, at all monitoring HABs.
Figure 11. Plot of vector field of W-velocity of the two cases, at all monitoring HABs.
Applsci 09 03989 g011
Figure 12. Peak value of W-velocity of the two cases, at all monitoring HABs.
Figure 12. Peak value of W-velocity of the two cases, at all monitoring HABs.
Applsci 09 03989 g012
Figure 13. Axial velocity contour of a cross-sectional plane of the 01 slit case (a), and the 02 slit case (b).
Figure 13. Axial velocity contour of a cross-sectional plane of the 01 slit case (a), and the 02 slit case (b).
Applsci 09 03989 g013
Figure 14. Plot of axial velocity, along with the radial direction at three monitoring HABs: 01 slit case (a), and 02 slit case (b).
Figure 14. Plot of axial velocity, along with the radial direction at three monitoring HABs: 01 slit case (a), and 02 slit case (b).
Applsci 09 03989 g014
Figure 15. Plot of U, V and radial velocity along with the radial direction at three monitoring HABs: U-velocity of the 01 slit case (a), U-velocity of the 02 slit case (b), V-velocity of the 01 slit case (c), V-velocity of the 02 slit case (d), radial velocity of the 01 slit case (e), radial velocity of the 02 slit case (f).
Figure 15. Plot of U, V and radial velocity along with the radial direction at three monitoring HABs: U-velocity of the 01 slit case (a), U-velocity of the 02 slit case (b), V-velocity of the 01 slit case (c), V-velocity of the 02 slit case (d), radial velocity of the 01 slit case (e), radial velocity of the 02 slit case (f).
Applsci 09 03989 g015
Figure 16. Plot of tangential velocity along with the radial direction at three monitoring HABs: 01 slit case (a), and 02 slit case (b).
Figure 16. Plot of tangential velocity along with the radial direction at three monitoring HABs: 01 slit case (a), and 02 slit case (b).
Applsci 09 03989 g016
Table 1. Temperature core radius at three monitoring HABs: 01 slit case and 02 slit case.
Table 1. Temperature core radius at three monitoring HABs: 01 slit case and 02 slit case.
HAB (m)Temperature Core Radius (cm)
01 Slit02 Slit
0.11.761.36
0.36.603.64
0.516.723.64
Table 2. Mean vortex core radius at three monitoring HABs: 01 slit case and 02 slit case.
Table 2. Mean vortex core radius at three monitoring HABs: 01 slit case and 02 slit case.
HAB (m)Mean Vortex Core Radius (cm)
01 Slit02 Slit
0.13.133.01
0.32.333.13
0.55.522.33

Share and Cite

MDPI and ACS Style

Wang, C.; Yuen, A.C.Y.; Chan, Q.N.; Chen, T.B.Y.; Chen, Q.; Cao, R.; Yip, H.L.; Kook, S.; Yeoh, G.H. Influence of Eddy-Generation Mechanism on the Characteristic of On-Source Fire Whirl. Appl. Sci. 2019, 9, 3989. https://doi.org/10.3390/app9193989

AMA Style

Wang C, Yuen ACY, Chan QN, Chen TBY, Chen Q, Cao R, Yip HL, Kook S, Yeoh GH. Influence of Eddy-Generation Mechanism on the Characteristic of On-Source Fire Whirl. Applied Sciences. 2019; 9(19):3989. https://doi.org/10.3390/app9193989

Chicago/Turabian Style

Wang, Cheng, Anthony Chun Yin Yuen, Qing Nian Chan, Timothy Bo Yuan Chen, Qian Chen, Ruifeng Cao, Ho Lung Yip, Sanghoon Kook, and Guan Heng Yeoh. 2019. "Influence of Eddy-Generation Mechanism on the Characteristic of On-Source Fire Whirl" Applied Sciences 9, no. 19: 3989. https://doi.org/10.3390/app9193989

APA Style

Wang, C., Yuen, A. C. Y., Chan, Q. N., Chen, T. B. Y., Chen, Q., Cao, R., Yip, H. L., Kook, S., & Yeoh, G. H. (2019). Influence of Eddy-Generation Mechanism on the Characteristic of On-Source Fire Whirl. Applied Sciences, 9(19), 3989. https://doi.org/10.3390/app9193989

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop