Next Article in Journal
Tissue-Level Flammability Testing: A Review of Existing Methods and a Comparison of a Novel Hot Plate Design to an Epiradiator Design
Previous Article in Journal
A Probabilistic Model for Fire Temperature Rise in High-Rise Residential Buildings under the Action of Uncertain Factors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation of the Impact of Firebrands on the Process of the Wood Layer Ignition

1
Department of Automobile Roads, Tomsk State University of Architecture and Building, 634003 Tomsk, Russia
2
Department of Physical and Computational Mechanics, Tomsk State University, 634050 Tomsk, Russia
*
Author to whom correspondence should be addressed.
Fire 2023, 6(4), 148; https://doi.org/10.3390/fire6040148
Submission received: 11 March 2023 / Revised: 30 March 2023 / Accepted: 4 April 2023 / Published: 6 April 2023

Abstract

:
In this study, a theoretical formulation of the ignition and combustion of the wood layer by burning and smoldering firebrands has been considered. The effect of the firebrands’ length, distances between firebrands and their geometrical parameters on the heat exchange with the wood layer and the ignition process were analyzed. With a decrease in firebrand size, ignition of wood is possible with a decrease in the distance between the firebrands. With an increase in firebrand size at the same distance between them, the ignition regime becomes possible albeit with a longer delay time Δt. With a decrease in the distance between the firebrands, the ignition of wood is possible with an increase in Δt. As a result of mathematical modeling of the process, the following processes are noted: the heat stored in firebrands of small sizes is insufficient to initiate the ignition process; the temperature in the wood layer, due to conductive heat exchange, slightly increases at first, before beginning to decrease as a result of heat exchange with the surrounding air and the wood layer; intensive heat exchange with the environment of small size firebrands leads to the end of firebrand smoldering and its cooling; and, if the firebrand size reaches a critical value, then the pyrolysis process begins in the area adjacent to it.

1. Introduction

Forest ecosystems are subject to numerous destructive impacts, such as wildland fires, mass reproduction of insect pests, deforestation, and other natural and anthropogenic factors. Wildland fires have a dominant negative impact on the state and dynamics of forest ecosystems [1,2], causing significant material and environmental damage.
Burning and smoldering firebrands of natural origin are the particles that have been generated as a result of heating and deformation of flammable natural materials, such as shrubs, trees (or any other combustible natural material), or building materials [3]. One of the effects observed in large-scale wildfires is burning and smoldering firebrands formed in the fire front. They are capable of covering a distance of several kilometers and initiating a new combustion center [4].
Recently, studies have been actively conducted that aim to systematically investigate the formation of firebrands, including both full-scale studies of the ignition of building structures [5,6,7,8] and modeling of the interaction of firebrands with structures based on wood [9,10,11,12,13,14,15,16]. The mechanism of ignition and combustion of the single firebrand model of various shapes, as well as the interaction of such a point source of thermal exposure with forest fuels, has been studied detailly in past studies [17,18,19,20,21,22].
For some time, attempts have been made to quantify the critical conditions for ignition caused by the accumulation of burning and smoldering firebrands. The thermal characteristics of a pile of firebrands on an inert plate were studied in a number of works [16,23,24]. These experiments showed that heating from a pile of firebrands is affected not only by the speed of the surrounding wind, but also by the mass and bulk density of the accumulated firebrands. It is, thus, of interest to study the process of firebrand accumulation near various types of building and structural materials to assess its tendency to ignite. In particular, a series of experiments was carried out in a past study [25] with 9 firebrands (3 by 3 square) applied to birch plywood 6.35 mm thick. The mass of the firebrands was 13–15 g. The distance between firebrands varied from 10 to 30 mm. The firebrands were ignited and placed on plywood. Various wind speeds were taken: 0, 0.05, and 0.75 m/s. The results show that firebrands with a spacing of 20 mm or more were able to burn through only the plywood surface until completely burning out. When the distance between the firebrands is less than 20 mm, plywood ignites and continues to burn even after the firebrands have completely burned out.
It was concluded that a single firebrand cannot cause plywood to ignite under the chosen parameters of the experiment. Firebrands located at a distance of 30 mm or more from each other cannot cause a sustained fire either. When the distance between the firebrands decreased, the plywood ignited and the fire continued to grow even when the firebrands stopped burning. This effect was also observed by Kwon and Liao [26]. In these experiments, firebrands were simulated using nine wooden cubes, with 19 mm on each side. The results showed that the flame height and the sample mass loss rate had non-monotonic dependencies based on the gap spacing. For smoldering combustion, compared to a single burning sample, the smoldering temperature and duration significantly increased due to the thermal interactions between adjacent burning samples [27].
Another study [28] presented an experimental and numerical study of how the thermal interactions of the fuel bed affect the subsequent tendency to ignite wood in contact with two electric heaters acting as idealized firebrands. The authors tried to find out the critical conditions (distance between firebrands, size, etc.) at which two idealized firebrands together can ignite solid wood fuel. Two sizes of heaters were used in the study: width 7.5 and 15 mm; length 25 and 35 mm, respectively.
Experimental studies show that for two heaters of both sizes, reducing the distance between them leads to the appearance of igniting ability. It is shown that the proximity of two heaters not only leads to an increase in the range of conditions under which ignition occurs, but also reduces the process of fuel ignition (dry white pine wood was used as fuel). In addition, the authors experimentally confirm that when several firebrands are placed in close proximity or collected in compact piles, an increase in the number of firebrands creates the opportunity for more heat to be released, which is subsequently transferred to the wood. Moreover, the accumulated firebrands can isolate each other, causing most of their energy to be transferred to the fuel rather than the environment. Thus, the accumulation of combustible mixture is considered to act as a more favorable source of ignition, co-igniting the target fuel and resulting in faster or more probable ignition compared to a similar situation with only one combustible mixture.
Apart from laboratory experiments, large-scale experimental studies were carried out to study the ability to ignite structural fuel by piles of burning and smoldering firebrands [29,30,31,32]. Some studies have used versions of the firebrand generator [33,34,35] to spray firebrands onto the fuel and observe ignition during their accumulation. Under certain conditions, individual or small piles of firebrands can ignite receptive fuels [36]; in other cases, they can accumulate and ignite larger fuels, such as wood covering [29]. The ability of single firebrands [20,37,38,39] or idealized firebrands [40,41,42] to ignite fuel was also studied in laboratory studies. The igniting possibilities of piles of firebrands to ignite a structural fuel have been investigated in laboratory experiments [16,42,43,44,45,46] that characterized the heat transfer from a pile of firebrands to the fuel and the subsequent ignition process. Heat exchange studies showed peak heat fluxes up to 70 kW/m2 transferred from smoldering firebrands [16,23,24,47,48]. Some studies have shown that piles of firebrands not only form a large heated area, but also cause a higher heat flux to the fuel [16,23,24]. Numerical or analytical studies of fuel heating and/or ignition by firebrands have been carried out, but they have mainly focused on single firebrands [21,22,38,49]; one study also considered small compact piles [39]. Several numerical studies assessed firebrands involvement in the co-ignition of various fuel types [50,51,52].
Numerous experiments [27,53,54] have been carried out on full-size structural elements, such as decks, fences, and building corners, subjected to firebrands’ shower using the Fire Dragon. The authors found that these constructions ignited more easily from a pile of firebrands than a single firebrand. It was shown in past studies [23,48,55,56] that the condition for ignition by a pile and a single firebrand differ.
It is necessary to study the ignition process for firebrands of specific sizes; a single hit on the wood layer will not lead to ignition. Indeed, if a single firebrand is able to initiate the process of ignition and combustion, the hit of a combination of such firebrands on the wood layer will also cause its ignition. Thus, in mathematical modeling it was assumed that the ignition of the wood layer by a single firebrand does not occur. However, its ignition becomes possible when two firebrands hit the surface. The fulfillment of this assumption was ensured by the choice of the appropriate firebrand sizes.
The present study focuses on studying wood layer ignition as a result of the accumulation of model firebrands simulating pine twigs. The novelty concludes in the investigation of the ignition process by firebrands falling on the surface non-simultaneously.

2. Materials and Methods

Mathematical Model

The ignition and combustion of the wood layer by burning and smoldering firebrands was considered. It was assumed that the combustible components, which made up the volatile products of pyrolysis, could be modeled by one effective combustible gas with the reactive properties of carbon monoxide. The rectangular coordinate system was used for a mathematical description: the OX and OY axes were directed horizontally, while the OZ axis was directed vertically upwards (Figure 1).
The balance of thermal energy in a smoldering firebrand and a layer of wood is described by the thermal conductivity equation. This equation has the form for a smoldering firebrand [57]:
C f b ρ f b T f b t = x λ f T f b x + y λ f T f b y + z λ f T f b z + Q t d Φ t d + Q p y r Φ p y r + ε * U σ T f b 4 δ .
The thermal conductivity equation for the wood layer can be written as:
C w ρ w T w t = x λ w T w x + y λ w T w y + z λ w T w z + Q t d Φ t d Q d r Φ d r + Q p y r Φ p y r + ε * U σ T w 4 δ ,
where C , ρ , λ , ε * are the heat capacity, density, thermal conductivity, and emissivity of the condensed phase. Index f b is related to the firebrand, w —to the wood layer.
Radiation heat transfer is determined by the P 1 approximation. The radiation intensity is determined from the equation solution [57,58]:
2 U x 2 + 2 U y 2 + 2 U z 2 = 3 β χ U σ T g 4 ,
where β is the integral radiation absorption coefficient, χ is the integral radiation attenuation coefficient, and σ is the black body radiation constant.
Evaporation of moisture associated with the wood layer can be described by the following dependence [59]:
Φ d r = k d r ρ f M h u m T 2.25 e x p T d r T f ,
where k d r = 6 · 10 5 K 2.25 / s 1 , T d r = 6 · 10 4 K is a pre-exponential factor and an evaporation activation temperature. The heat of evaporation reaction was assumed to be Qdr = 3 · 106 J/kg. Mhum is the relative wood humidity.
The pyrolysis mass rate of the condensed phase (firebrand and layer of wood) can be calculated using the dependence [57,59]:
Φ p y r = k p y r ρ f M 1 e x p T p y r T f
where thermokinetic parameters have the following values: k p y r = 3.63 · 10 4 s 1 ,   T p y r = 9.4 · 10 3 K ,   Q p y r = 10 7 J/kg. M 1 is the mass fraction of dry organic material, f = f b for igniting firebrands, and f = w for the wood layer.
The combustion rate of condensed pyrolysis products are described by the following dependence [57,59]:
Φ t d = k t d ρ f M O 2 M p y r e x p T t d T f
where k t d = 10 6 s 1 ,   T t d = 10 4 K —pre-exponential index and temperature of the heterophase combustion reaction activation. The heat of combustion reaction of condensed pyrolysis products is taken equal to Q d r = 1.2 · 10 7 J/kg. M p y r is the mass fraction of condensed pyrolysis products.
The mass change in the wood layer, water change in a liquid-drop state, and change in the condensed pyrolysis products are described by the balance equations:
ρ f M 1 t = Ф p y r , ρ f M 2 t = Ф d r , ρ f M 3 t = Ф p y r Ф t d
The thermal conductivity equation of the gas phase can be written as [60,61,62,63]:
C p , g ρ g T g t + ρ g v x T g x + ρ g v y T g y + ρ g v z T g z = x λ g T g x + y λ g T g y + z λ g T g z + Q g Φ g + h ˙ + β U σ T g 4
The enthalpy change in the gas phase due to the influx of gaseous pyrolysis products and water vapor into it, which are formed during the drying of wood, can be calculated as follows:
h ˙ = h ˙ p y r + h ˙ d r ,   h ˙ p y r = C p , C O T f T g Φ p y r Φ t d ,   h ˙ d r = C p , H 2 O T f T g Φ d r .
The chemical reaction rate in the gas phase is described by the dependence [64,65]:
Ф g = k g ρ g 2 M O 2 M C O e x p E g R g T g
The chemical reaction parameters were taken in accordance with the data of two past studies [66,67]: k g = 2.2 · 10 8 · m 3 · k g 1 · s 1 , E g = 104   k J · m o l 1 .
The mass balance equations for the components O 2 , C O , C O 2 , N 2 , H 2 O were used [68,69,70] to describe the processes of diffusion, mixing, chemical reaction, and combustion in the gas phase, in addition to the energy equation:
ρ g M C O t + ρ g v x M C O x + ρ g v y M C O y + ρ g v z M C O z = x ρ g D M C O x + y ρ g D M C O y + z ρ g D M C O z 2 W C O W O 2 Φ g + Φ p y r ,
ρ g M O 2 t + ρ g v x M O 2 x + ρ g v y M O 2 y + ρ g v z M O 2 z = x ρ g D M O 2 x + y ρ g D M O 2 y + z ρ g D M O 2 z Φ g ,
ρ g M C O 2 t + ρ g v x M C O 2 x + ρ g v y M C O 2 y + ρ g v z M C O 2 z = x ρ g D M C O 2 x + y ρ g D M C O 2 y + z ρ g D M C O 2 z + 2 W C O 2 W O 2 Φ g ,
ρ g M N 2 t + ρ g v x M N 2 x + ρ g v y M N 2 y + ρ g v z M N 2 z = x ρ g D M N 2 x + y ρ g D M N 2 y + z ρ g D M N 2 z ,
ρ g M H 2 O t + ρ g v x M H 2 O x + ρ g v y M H 2 O y + ρ g v z M H 2 O z = x ρ g D M H 2 O x + y ρ g D M H 2 O y + z ρ g D M H 2 O z + Φ H 2 O .
The state equation of the gas phase has the form:
ρ g = p g R T g M C O W C O + M O 2 W O 2 + M C O 2 W C O 2 + M N 2 W N 2 + M H 2 O W H 2 O 1 .
The Navier–Stokes equations were used in the Oberbeck–Boussinesq approximation to describe the flow field, which in the rectangular coordinate system has the form [71,72]:
ρ v x x + ρ v y y + ρ v z z = 0 ,
ρ v x t + ρ ν x 2 x + ρ v y v x y + ρ v z v x z = p x + x μ v x x + y μ v x y + z μ v x z ,
ρ v y t + ρ v x v y x + ρ ν y 2 y + ρ v z v x z = p y + x μ v y x + y μ v y y + z μ v y z ,
ρ v z t + ρ v x v z x + ρ v y v z y + ρ ν z 2 z = p z + x μ v z x + y μ v z y + z μ v z z + ρ ρ 0 g .
It was assumed in the flow structure model that the air movement near the surface occurred in a laminar regime, while the effect of the influx of pyrolysis products and evaporated water into the gas phase on the flow structure was negligibly small.
Boundary conditions for a system of differential equations are determined by the boundary type. Uniform distributions of temperature and mass fraction of the gas phase components are stated for the windward side. Soft boundary conditions were set for the leeward, lateral side, and upper boundaries. The wind speed was assumed to be constant and independent of height. The temperatures of the wood layer and air were equal at the initial moment of time. The firebrand temperature corresponded to the smoldering temperature of tree bark and twigs of conifers.
The equations in this work presented a completely closed system of equations which, under appropriate boundary and initial conditions and known properties of the medium, determined the main characteristics of heat exchange. The equations were solved numerically using the finite volume method. In accordance with this method, finite difference equations were obtained by integrating differential equations over control volumes containing points of a finite difference grid. The calculations were carried out on a grid with 2000 points in the OX direction, 2500 nodes in the OY direction, and 1500 nodes in the OZ direction. The grid was thickened near the firebrands. The continuity equation was satisfied using the SIMPLEC algorithm. It was considered that the iteration convergence was achieved if the root–mean–square discrepancy for all variables did not exceed 1%. A series of calculations were performed on sequences of refining grids to assess the accuracy of the calculations. The test results showed that a 2-fold decrease in the base grid step along the axial and radial coordinates led to a change in the values of the main variables by no more than 1%.
The ignition of the wood layer by two firebrands of forest combustible material was simulated during the calculations. Pine twigs were chosen as the firebrands’ model.
To study the considered process, the evolution of isotherms on the surface of the wood layer was analyzed, along with the change in temperature over time at three characteristic points. Point (A) was chosen at a distance of 3 mm from the windward edge of the firebrand lying upstream. Point (B) was located between the firebrands at an equal distance from their edges. Point (C) was chosen at a distance of 3 mm from the leeward edge of the firebrand lying downstream. All control points were located on the OX axis. The layout of firebrands and control points is shown in Figure 2 (axes were chosen in the wind direction).

3. Results and Discussion

The following analysis of the results of mathematical modeling was carried out. Figure 3 shows the temperature change over time at the control point when two thin firebrands hit the surface simultaneously. One can see from Figure 3 that the heat exchange in thin firebrands is carried out in a low-temperature regime. The temperature of the wood layer increases at first, which is associated with the wood being heated by the firebrand. As a result of heat exchange with the environment, cooling then occurs after the thermal energy reserves are exhausted and the temperature decreases. The results of calculations show that the firebrands’ length has practically no effect on the process of heating the wood layer during the first three seconds. However, since there is an increase in the length of the firebrands, an increase in the thermal energy reserves occurs and the transition time for heating the wood layer increases. As can be seen from Figure 3, the coordinate of the maximum temperature curve shifts to the right with a rise in L. At the same time, an increase in the maximum temperature at the control point is observed. To initiate the ignition process, a glowing firebrand must have significant thermal energy, which depends on the size of the firebrands. The calculation results show that the ignition of the wood layer by two thin firebrands with a diameter of d < 5 mm located at a distance of more than 3 mm is impossible. However, in the case of firebrands with a diameter of more than 5 mm, a high-temperature heat exchange regime becomes possible, leading to ignition of the wood surface.
The effect of the firebrands’ length on the heat exchange with the wood layer and the ignition process is shown in Figure 4. A low-temperature regime is realized for firebrands of small length. As the firebrand length increases, the temperature of the wood layer in its vicinity increases. The wood surface ignites if the firebrand length exceeds a critical value, which depends on other firebrands’ sizes and heat transfer conditions.
A further increase in the firebrand length leads to a decrease in the induction period and earlier ignition. Curves 1–6 characterize the heat exchange in thick firebrands and refer to the high temperature regime. The temperature change over time is close to linear in the initial time period, which corresponds to the heating and drying of the wood layer. An exponential rise in temperature occurs after the wood layer is ignited. The temperature rate increase slows down as the burning increases. The temperature approaches the adiabatic combustion temperature. However, it does not reach its value due to heat exchange. The temperature then decreases at the considered points as the burning proceeds. The cooling process that follows the wood layer burning out occurs at sufficiently long times and is not shown in the Figure.
Figure 5 shows the effect of distances between firebrands on heat exchange and ignition processes. As a single firebrand of the considered sizes does not initiate the ignition of the wood layer, the high-temperature heat exchange regime in the vicinity of the considered point is realized as a result of the energy supply from both firebrands. The calculation results show that two thin twigs located in close proximity to each other ( h = 0 ) cause ignition of the wood layer. At the same time, the processes of heating and drying proceed with high intensity, which leads to rapid wood ignition.
If there is some distance between the firebrands, part of their thermal energy will be spent on heating the gap between them. As a result, the heating rate of the outer domain slows down. However, the energy supplied from them may be sufficient for ignition if the firebrands are located close to each other. Thus, for closely spaced firebrands, their thermal interaction with the wood layer occurs in the ignition regime. The ignition time increases nonlinearly with increasing distance between them and tends to infinity at h h * . The distance between the firebrands h * determines the boundary of the transition from the high-temperature regime to the low-temperature regime. In the low-temperature regime (at h h * ), the wood layer cools down after a slight initial heating.
Figure 6 shows the ignition process of the wood layer with two thick twigs. At the initial moment of time, the wood layers adjacent to it are heated as a result of heat exchange between the wood layer and a twig. The heating area around each twig increases. The interaction processes of smoldering twigs with the heating layer of the wood layer are intensified when the closing of thermal fronts occurs. As a result, the drying process is intensified. The pyrolysis reaction is initiated after the moisture is removed from the wood. The gaseous pyrolysis products mix with atmospheric oxygen and an inflammation occurs. A high-temperature combustion zone is formed near the surface of the firebrand, which increases in size with time. After the pyrolysis products burn out, the afterburning of the condensed products occurs in the smoldering regime.
The non-simultaneous reaching of the wood surface by igniting firebrands is considered (Figure 7). While carrying out the main series of studies, it was assumed that the wood surface is first reached by firebrand 1, and later by firebrand 2, which is located on the leeward side with respect to firebrand 1.
The thermal energy stored in a small firebrand is not enough to initiate ignition. Therefore, the interaction of such firebrand with a wood surface occurs in a low-temperature regime. The adjacent wood layers are heated as a result of heat exchange with the firebrand (Figure 7b–d). The heating area increases in the initial period of time. An increase in the area of the heating zone leads to the cooling of the firebrand. As a result, the intensity of smoldering in the firebrand weakens and eventually stops. Subsequently, heat exchange with the environment leads to the cooling of the firebrand and the wood layer adjacent to it. The heat transfer regime changes significantly after the hitting of the surface by the second firebrand. The wood layer is heated and dried by the time the second firebrand hits the surface. The pyrolysis reaction is intensified in the dried and heated wood layer. The gaseous pyrolysis products mix with atmospheric oxygen, the temperature in the channel wall reaches the adiabatic combustion temperature T a , and an inflammation occurs. A high-temperature combustion zone is formed near the firebrand surface, which increases in size over time.
The effect of geometrical parameters of firebrands on heat exchange and ignition is also considered. Figure 8 shows the temperature change over time at three control points: A, B, C; these points are calculated for various firebrands’ lengths. One can see from the Figure 8 below that the heat exchange in firebrands of small length is carried out in a low-temperature regime. The temperature of the firebrands approaches the ambient temperature in this case. The heat exchange in long firebrands proceeds in a high-temperature regime, leading to the ignition of the wood layer.
The temperature change at the control points in the initial period of time (until the second firebrand lands on the surface) is typical for the low-temperature regime. The wood layers adjacent to the firebrands are cooled after a slight heating. In this case, the firebrand temperature approaches the ambient temperature. The cooling process after the wood layer burns out occurs over a sufficiently long time period and is not shown in the Figure.
After the firebrand hits the wood layer, a temperature boundary layer is formed, which spreads over the wood surface over time. After the temperature front reaches the control point, its vicinity begins to warm up. The lowest heating temperature is reached at point C, which is the furthest point from the firebrand edge; the highest one is at point A, which is closest to the firebrand edge. The temperature rises sharply after the second firebrand hits the wood surface at points B and C. At point A, which is the most distant from the second firebrand, the temperature still decreases for some time. The heating of the wood layer begins only after a period of time equal to δ t = λ W C w ρ w X A 2 ( X A 2 is the distance from point A to the windward edge of the second firebrand). After hitting the surface by the second firebrand, intense heating of the area bounded by these firebrands occurs. Thus, the region of the highest wood temperature is localized. In this case, pyrolysis and ignition occur if the intensity of heating of the wood layer by smoldering firebrands exceeds the intensity of heat losses. Firstly, the wood layers adjacent to the left edge of the second firebrand are ignited, followed by the right edge. The mixing of pyrolysis products with atmospheric oxygen leads to an inflammation, accompanied by a sharp, pike-shaped increase in temperature. After the gaseous pyrolysis products in a current domain burn out, the combustion regime is replaced by a smoldering one. In this case, the pyrolysis area extends into the region external to the firebrand. If the firebrand size is sufficiently large, a flame combustion regime occurs in the gas phase after ignition of the wood layer. After the pyrolysis products burn out, the burning of the condensed products occurs in the smoldering regime.
The analysis of the firebrand diameter effect on the ignition conditions of the wood surface will be considered (Figure 9). The dependences presented in the Figure below are similar to those presented above. A low-temperature heat exchange regime is realized for small-diameter firebrands. Comparison of curves 1, 2, illustrating the low-temperature regime of heat exchange, shows that the temperature of the wood layer slightly increases at all control points with an increase in the firebrand diameter. We also note that for point C, the increase in temperature with increasing diameter becomes noticeable only after the second firebrand hits the surface as the effect of firebrand 1 is small due to its distance from point C. A significant supply of thermal energy is required to ignite the wood layer, which can be possessed by large firebrands (curves 3–5). The change in temperature at point B for firebrands of average diameter (d = 5.5 mm) will be considered in more detail. In this case, before the second firebrand hits the surface, the surface is primarily heated by the first firebrand. The surface is later cooled as a result of heat transfer. After the second firebrand hit the surface, intensive heating and drying of wood occurs in the vicinity of point B. The beginning of the pyrolysis process characterizes the inflection point on the temperature curve and outlet on the horizontal asymptote indicates the transition to ignition. The difference of curve 3 at point C is the absence of the initial cooling stage, since the thermal wave from the first firebrand still continues to propagate. Moreover, the duration of the drying period, which precedes the start of pyrolysis, increases. This is explained by the fact that the temperature at point C at the moment the second firebrand hits the surface is lower than the temperature at point B. The cooling time of the wood layer at point A is much longer than at point B. At the same time, the start of secondary heating at point A occurs when the transition to ignition is carried out at point B. The time preceding the ignition of the firebrands decreases with a further increase in the firebrand diameter. In the case of oversized firebrands (d = 7.5 mm), when the surroundings of the first firebrand are intensively heated, pyrolysis and ignition of the wood layer is carried out immediately after the second firebrand hits the surface: first at point B, then at A, and lastly at C.
Figure 10 shows the change in temperature at the control points over time at various distances between the ignition firebrands. One can see from the Figure that the main regularities of the temperature distribution evolution are preserved. At small distances between the firebrands, the process of the second firebrand’s interaction with the wood layer after hitting the surface passes successively through all the main stages: heating, drying, pyrolysis, ignition, and combustion. A high-temperature combustion zone is formed near the firebrand surface, which increases in size over time. The heating and drying time increases significantly with an increase in the distance between the firebrands. The period duration of the transition to pyrolysis and ignition at points B and, especially, at C increases significantly. At point A, the cooling time increases significantly. In this case, ignition in the vicinity of this point occurs after the formation of the combustion zone at points B and C. A low-temperature heat exchange regime is realized at large distances between firebrands. The thermal energy stored in the firebrand is not sufficient for further heating, drying, and initiation of a chemical reaction. Thermal energy is reduced as a result of heat exchange with air and the wood layer, which leads to cooling and a decrease in temperature. The effect of the distance between the firebrands on the temperature change at the control points becomes insignificant due to their considerable long distance.
Figure 11 shows the change in temperature at the control points for various times of reaching the surface by the second firebrand Δ t.
If the firebrands reach the surface at the same time Δ t = 0, the high-temperature heat exchange regime is realized. Curve 1 (points A and C) clearly shows that there is an increase in temperature in the initial period of time, which indicates the heating of the wood. The drying process then begins and the temperature curves reach saturation, after which pyrolysis and ignition of the wood take place. At point B, which is under the effect of two firebrands, these processes proceed much faster, which leads to a decrease in the ignition time. One can see an interesting effect at small delay times Δ t < 1 s. If an increase in Δ t at point A leads only to a slight increase in the drying time, the heating of the point C vicinity occurs more intensively as the second firebrand falls on a heated and partially dried surface.
As a consequence of this effect, the ignition time at point C decreases compared to the case Δ t = 0. The same effect, however, shows itself in a noticeably weaker form at point B. A further increase in Δ t leads to an increase in the ignition time and a displacement of the temperature curves to the right. This statement is true as long as the high-temperature regime of heat exchange is realized. At long delay times, Δ t > 3c firebrand falls on an already cooled surface and cannot cause heating sufficient for ignition. The low-temperature regime is implemented.
The results of the carried out studies are summarized in Figure 12, which shows the dependence of the critical delay time Δ t on the size of the firebrands and the distance between them. The domain, limited by the coordinate axes and curves in the Figure, corresponds to the implementation of the ignition regime. The outer domain corresponds to the low-temperature regime. One can see from the Figure that the ignition regime becomes possible with an increase in the firebrand size and a longer delay time Δ t. The wood ignition is possible with an increase in the distance between the firebrands and a decrease in Δ t.

4. Conclusions

The model presented considers the heat exchange between the firebrands, wood layer, and the gas phase, as well as moisture evaporation in the firebrands and the diffusion gases of water vapor in the pyrolysis zone. The research results show how the non-simultaneous falling of firebrands onto the wood layer, as well as the distance between them, impacts the ignition and combustion process. The following conclusions can be made as a result of mathematical modeling:
1. The heat stored in firebrands of small sizes is insufficient to initiate the ignition process. The temperature in the wood layer, due to conductive heat exchange, slightly increases at first. Afterwards, it begins to decrease as a result of heat exchange with the surrounding air and the wood layer. Intensive heat exchange with the environment of small size firebrands leads to the end of firebrand smoldering and its cooling. If the firebrand size reaches a critical value, then the pyrolysis process begins in the area adjacent to it.
2. The effect of a set of firebrands (thin twigs) on the process of the wood layer ignition allows one to conclude the following. With a decrease in firebrand size, ignition of wood is possible with a decrease in the distance between the firebrands. With an increase in firebrand size at the same distance between them, the ignition regime becomes possible with a longer delay time Δ t. With a decrease in the distance between the firebrands, the ignition of wood is possible with an increase in Δ t. Results from this study demonstrate the significant influence of spacing between the firebrands on the ignition and the burning behavior of the wood layer. It is note that the narrow spacing between the firebrands significantly increases the likelihood of fire on a building material.

5. Limitations

  • The blown air flow was assumed to be laminar. Turbulence was not considered.
  • The wind speed was assumed to be constant over time.
  • The wind direction was chosen to be perpendicular to the long side of the firebrand.
  • The firebrands were assumed to have a regular shape (rectangular parallelepiped) with uniform thermophysical properties. The temperature distribution in the firebrands at the initial moment of time was not taken into account (the temperature at all points of the firebrand had the same value).
  • Ideal contact of the firebrand surface with the wood layer was assumed.

Author Contributions

Conceptualization, D.K., O.M. and E.L.; methodology, D.K. and O.M.; validation, D.K., E.L. and O.D.; formal analysis, D.K. and O.M.; investigation, D.K., O.M. and A.L.; writing—original draft preparation, D.K., O.D. and O.M.; writing—review and editing, D.K., E.L. and O.M.; visualization, D.K., O.M., O.D. and A.L.; supervision, O.M.; project administration, D.K.; funding acquisition, D.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by RUSSIAN SCIENCE FOUNDATION, grant number 20-71-10068. Link to information about the project: https://rscf.ru/en/project/20-71-10068/ (accessed on 11 January 2022).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

TermMeaning (Units)
Nomenclature
C heat capacity ( J k g 1 K 1 )
ρ density ( k g   m 3 )
T temperature ( K )
λ thermal conductivity ( W m 1 k 1 )
U radiation intensity (W m 2 )
δ is the characteristic particle size ( m )
Q reaction heat of pyrolysis, drying and combustion of coke ( J k g 1 )
Ф the rate of pyrolysis, drying, and combustion (coke)
ε * is the degree of blackness of the wood layer
β integral radiation absorption coefficient, ( m 1 )
χ integral radiation attenuation coefficient, ( m 1 )
σ is the black body radiation constant
k pre-exponential factor
M h u m relative moisture content of wood sample
M 1 m a s s   f r a c t i o n   o f   d r y   o r g a n i c   m a t t e r
M p y r m a s s   f r a c t i o n   o f   c o n d e n s e d   p y r o l y s i s   p r o d u c t s
M m a s s   f r a c t i o n
h ˙ rate of the enthalpy change, ( J k g 1 s 1 )
E activation energy, (J m o l e 1 )
v wind speed, ( m s 1 )
g gravity factor, (m s 1 )
Г coefficient of diffusion, ( m 2 s 1 )
W molecular mass, (kg m o l e 1 )
p pressure, P a
μ viscosity, (Pa s)
R g universal gas constant, ( J m o l e 1 K 1 )
C p , g heat capacity of gas at constant pressure, (Pa)
ttime, (s)
x, y, zcoordinates, (m)
Subscripts
f b f i r e b r a n d
f solid phase (wood, firebrands)
t d t h e r m a l   d e g r a d a t i o n
p y r pyrolysis
d r drying
w wood layer
g gas
1dry organic substance
2water in the liquid-drop condition
3condensed pyrolysis products
O 2 oxygen
C O carbon monoxide
C O 2 carbon dioxide
N 2 nitrogen
H 2 O water
Oambient

References

  1. Katurji, M.; Noonan, B.; Zhang, J.; Valencia, A.; Shumacher, B.; Kerr, J.; Strand, T.; Pearce, G.; Zawar-Reza, P. Atmospheric turbulent structures and fire sweeps during shrub fires and implications for flaming zone behaviour. Int. J. Wildland Fire 2023, 32, 43–55. [Google Scholar] [CrossRef]
  2. Voulgarakis, A.; Field, R.D. Fire Influences on Atmospheric Composition, Air Quality and Climate. Curr. Pollut. Rep. 2015, 1, 70–81. [Google Scholar] [CrossRef] [Green Version]
  3. Babrauskas, V. Firebrands and Embers. In Encyclopedia of Wildfires and Wildland-Urban Interface (WUI) Fires; Manzello, S., Ed.; Springer: Cham, The Switzerland, 2018. [Google Scholar] [CrossRef]
  4. Koo, E.; Pagni, P.J.; Weise, D.R.; Woycheese, J.P. Firebrands and spotting ignition in large-scale fires. Int. J. Wildland Fire 2010, 19, 818–843. [Google Scholar] [CrossRef] [Green Version]
  5. Suzuki, S.; Manzello, S.L. Characteristics of Firebrands Collected from Actual Urban Fires. Fire Technol. 2018, 54, 1533–1546. [Google Scholar] [CrossRef]
  6. Suzuki, S.; Manzello, S.L.; Lage, M.; Laing, G. Firebrand generation data obtained from a full-scale structure burn. Int. J. Wildland Fire 2012, 21, 961–968. [Google Scholar] [CrossRef] [Green Version]
  7. Zhou, A.; Quarles, S.L.; Weise, D.R. Final Report: Fire Ember Production from Wildland and Structural Fuels; JFSP Project ID: 15-1-04-4; North Carolina Agricultural and Technical State University: Greensboro, NC, USA, 2019; pp. 1–66. [Google Scholar] [CrossRef]
  8. Quarles, S.L. Vulnerability of Vents to Wind-Blown Embers; IBHS: Richburg, SC, USA, 2017. [Google Scholar]
  9. Suzuki, S.; Manzello, S.L. Initial study on thatched roofing assembly ignition vulnerabilities to firebrand showers. Fire Saf. J. 2019, 103, 34–37. [Google Scholar] [CrossRef]
  10. Kasymov, D.P.; Paletsky, A.A. Effect of a fire retardant on the ignition of pine wood exposed to smoldering particles of pine bark. In Proceedings of the EPJ Web of Conferences, Crete, Greece, 17–29 August 2017; Volume 159, pp. 1–26. [Google Scholar] [CrossRef] [Green Version]
  11. Suzuki, S.; Manzello, S.L. Firebrand production from building components fitted with siding treatments. Fire Saf. J. 2016, 80, 64–70. [Google Scholar] [CrossRef] [Green Version]
  12. Suzuki, S.; Manzello, S.L. Garnering understanding into complex firebrand generation processes from large outdoor fires using simplistic laboratory-scale experimental methodologies. Fuel 2020, 267, 117154. [Google Scholar] [CrossRef]
  13. Kasymov, D.P.; Agafontsev, M.V.; Perminov, V.V.; Martynov, P.S. Estimation of pine wood fire-resistance exposed to the heat emission using thermography. In Proceedings of the 25th International Symposium on Atmospheric and Ocean Optics: Atmospheric Physics, Tomsk, Russian, 30 June–5 July 2019; Volume 11208. [Google Scholar] [CrossRef]
  14. Suzuki, S.; Manzello, S.L.; Hayashi, Y. The size and mass distribution of firebrands collected from ignited building components exposed to wind. Proc. Combust. Inst. 2013, 34, 2479–2485. [Google Scholar] [CrossRef]
  15. Kasymov, D.P.; Tarakanova, V.A.; Martynov, P.S.; Agafontsev, M.V. Studying firebrands interaction with flat surface of various wood construction materials in laboratory conditions. J. Phys. Conf. Ser. 2019, 1359, 012092. [Google Scholar] [CrossRef]
  16. Salehizadeh, H.; Hakes, R.S.P.; Gollner, M.J. Critical Ignition Conditions of Wood by Cylindrical Firebrands. Front. Mech. Eng. 2021, 7, 630324. [Google Scholar] [CrossRef]
  17. Ellis, P.F. The Aerodynamic and Combustion Characteristics of Eucalypt Bark—A Firebrand Study. Ph.D. Thesis, Australian National University, Canberra, ACT, Australia, 2000. [Google Scholar]
  18. Ganteaume, A.; Lampin-Maillet, C.; Guijarro, M.; Hernando, C.; Jappiot, M.; Fonturbel, T.; Pérez-Gorostiaga, P.; Vega, J. Spot fires: Fuel bed flammability and capability of firebrands to ignite fuel beds. Int. J. Wildland Fire 2009, 18, 951–969. [Google Scholar] [CrossRef] [Green Version]
  19. Krupiczka, R. Analysis of thermal conductivity in granular materials. Int. J. Chem. Eng. 1967, 7, 122–144. [Google Scholar] [CrossRef]
  20. Manzello, S.L.; Cleary, T.G.; Shields, J.R.; Yang, J.C. On the ignition of fuel beds by firebrands. Fire Mater. 2006, 30, 77–87. [Google Scholar] [CrossRef]
  21. Lautenberger, C.; Fernandez-Pello, A.C. Generalized pyrolysis model for combustible solids. Fire Saf. J. 2009, 44, 819–839. [Google Scholar] [CrossRef] [Green Version]
  22. Matvienko, O.V.; Kasymov, D.P.; Filkov, A.I.; Daneyko, O.I.; Gorbatov, D.A. Simulation of fuel bed ignition by wildland firebrands. Int. J. Wildland Fire 2018, 27, 550–561. [Google Scholar] [CrossRef]
  23. Hakes, R.S.P.; Salehizadeh, H.; Weston-Dawkes, M.J.; Gollner, M.J. Thermal characterization of firebrand piles. Fire Saf. J. 2019, 104, 34–42. [Google Scholar] [CrossRef] [Green Version]
  24. Tao, Z.; Bathras, B.; Kwon, B.; Biallas, B.; Gollner, M.J.; Yang, R. Effect of firebrand size and geometry on heating from a smoldering pile under wind. Fire Saf. J. 2021, 120, 103031. [Google Scholar] [CrossRef]
  25. Kwon, B. In Proceedings of the Ignition Propensity of Structural Materials Exposed to Multiple Firebrands in Wildland-Urban Interface (WUI) Fires: Effects of Firebrand Distribution and Ambient Wind Spring Technical Meeting of the Central States Section of the Combustion Institute, Detroit, Michigan, 1 May 2022; pp. 1–10.
  26. Kwon, B.; Liao, Y.-T.T. Effects of spacing on flaming and smoldering firebrands in wildland–urban interface fires. J. Fire Sci. 2022, 40, 155–174. [Google Scholar] [CrossRef]
  27. Manzello, S.L.; Suzuki, S. The world is burning: What exactly are firebrands and why should anyone care? Front. Mech. Eng. 2023, 8, 1072214. [Google Scholar] [CrossRef]
  28. Rivera, J.; Hernández, N.; Consalvi, J.; Reszka, P.; Contreras, J.; Fuentes, A. Ignition of wildland fuels by idealized firebrands. Fire Saf. J. 2021, 120, 103036. [Google Scholar] [CrossRef]
  29. Manzello, S.L.; Suzuki, S. Experimental investigation of wood decking assemblies exposed to firebrand showers. Fire Saf. J. 2017, 92, 122–131. [Google Scholar] [CrossRef] [PubMed]
  30. Manzello, S.L.; Suzuki, S.; Hayashi, Y. Enabling the study of structure vulnerabilities to ignition from wind driven firebrand showers: A summary of experimental results. Fire Saf. J. 2012, 54, 181–196. [Google Scholar] [CrossRef]
  31. Manzello, S.L.; Park, S.-H.; Suzuki, S.; Shields, J.R.; Hayashi, Y. Experimental investigation of structure vulnerabilities to firebrand showers. Fire Saf. J. 2011, 46, 568–578. [Google Scholar] [CrossRef]
  32. Suzuki, S.; Manzello, S.L. Role of accumulation for ignition of fuel beds by firebrands. Appl. Energy Combust. Sci. 2020, 1–4, 100002. [Google Scholar] [CrossRef]
  33. Manzello, S.L.; Shields, J.R.; Cleary, T.G.; Maranghides, A.; Mell, W.E.; Yang, J.C.; Hayashi, Y.; Nii, D.; Kurita, T. On the development and characterization of a firebrand generator. Fire Saf. J. 2008, 43, 258–268. [Google Scholar] [CrossRef]
  34. Manzello, S.L.; Suzuki, S. Experimentally Simulating Wind Driven Firebrand Showers in Wildland-urban Interface (WUI) Fires: Overview of the NIST Firebrand Generator (NIST Dragon) Technology. Procedia Eng. 2013, 62, 91–102. [Google Scholar] [CrossRef] [Green Version]
  35. Kasymov, D.P.; Agafontsev, M.V.; Perminov, V.V.; Loboda, E.L.; Loboda, Y.A.; Reyno, V.V.; Orlov, K.E. Study of resistance to ignition of wood structures under the thermal impact of the flow of burning particles. Combust. Explos. Shock. Waves 2023, 2. in press. [Google Scholar]
  36. Manzello, S.L.; Park, S.-H.; Cleary, T.G. Investigation on the ability of glowing firebrands deposited within crevices to ignite common building materials. Fire Saf. J. 2009, 44, 894–900. [Google Scholar] [CrossRef]
  37. Viegas, D.; Almeida, M.; Raposo, J.; Oliveira, R.; Viegas, C. Ignition of Mediterranean Fuel Beds by Several Types of Firebrands. Fire Technol. 2014, 50, 61–77. [Google Scholar] [CrossRef]
  38. Urban, J.L.; Song, J.; Santamaria, S.; Fernandez-Pello, C. Ignition of a spot smolder in a moist fuel bed by a firebrand. Fire Saf. J. 2019, 108, 102833. [Google Scholar] [CrossRef]
  39. Wessies, S.S.; Chang, M.K.; Marr, K.C.; Ezekoye, O.A. Experimental and Analytical Characterization of Firebrand Ignition of Home Insulation Materials. Fire Technol. 2019, 55, 1027–1056. [Google Scholar] [CrossRef]
  40. Hernández, N.; Fuentes, A.; Consalvi, J.; Elicer-Cortés, J. Spontaneous ignition of wildland fuel by idealized firebrands. Exp. Therm. Fluid Sci. 2018, 95, 88–95. [Google Scholar] [CrossRef]
  41. Koo, E.; Rodman, R.; Linn, P.J.; Pagni, C.B. Edminster Modelling firebrand transport in wildfires using HIGRAD/FIRETEC. Int. J. Wildland Fire 2012, 21, 396–417. [Google Scholar] [CrossRef] [Green Version]
  42. Reszka, P.; Cruz, J.; Valdivia, J.; González, F.; Rivera, J.; Carvajal, C.; Fuentes, A. Ignition delay times of live and dead pinus radiata needles. Fire Saf. J. 2020, 112, 102948. [Google Scholar] [CrossRef]
  43. Manzello, S.L.; Cleary, T.G.; Shields, J.R.; Maranghides, A.; Mell, W.; Yang, J.C. Experimental investigation of firebrands: Generation and ignition of fuel beds. Fire Saf. J. 2008, 43, 226–233. [Google Scholar] [CrossRef]
  44. Manzello, S.L.; Cleary, T.G.; Shields, J.R.; Yang, J.C. Ignition of mulch and grasses by firebrands in wildland—Urban interface fires. Int. J. Wildl. Fire 2006, 15, 427–431. [Google Scholar] [CrossRef] [Green Version]
  45. Meerpoel-Pietri, K.; Tihay-Felicelli, V.; Santoni, P.-A. Determination of the critical conditions leading to the ignition of decking slabs by flaming firebrands. Fire Saf. J. 2021, 120, 103017. [Google Scholar] [CrossRef]
  46. Filkov, A.; Kasymov, D.; Zima, V.; Matvienko, O. Experimental investigation of surface litter ignition by bark firebrands. AIP Conf. Proc. 2016, 1698, 060004. [Google Scholar] [CrossRef] [Green Version]
  47. Bearinger, E.D.; Hodges, J.L.; Yang, F.; Rippe, C.M.; Lattimer, B.Y. Localized heat transfer from firebrands to surfaces. Fire Saf. J. 2020, 120, 103037. [Google Scholar] [CrossRef]
  48. Abul-Huda, Y.M.; Bouvet, N. Thermal dynamics of deposited firebrands using phosphor thermometry. Proc. Combust. Inst. 2021, 38, 4757–4765. [Google Scholar] [CrossRef]
  49. Warey, A. Influence of thermal contact on heat transfer from glowing firebrands. Case Stud. Therm. Eng. 2018, 12, 301–311. [Google Scholar] [CrossRef]
  50. Glushkov, D.O.; Kuznetsov, G.V.; Strizhak, P.A. Heat and mass transfer at gas-phase ignition of grinded coal layer by several metal particles heated to a high temperature. Thermophys. Aeromechanics 2017, 24, 593–604. [Google Scholar] [CrossRef]
  51. Glushkov, D.; Kuznetsov, G.; Strizhak, P. Ignition of a metallized composite solid propellant by a group of hot particles. Combust. Explos. Shock. Waves 2016, 52, 694–702. [Google Scholar] [CrossRef]
  52. Glushkov, D.O.; Strizhak, P.A. Mathematical Modeling of Heat and Mass Transfer Processes with Chemical Reaction at Polymeric Material Ignition by Several Small-Size Hot Particles. Math. Probl. Eng. 2015, 2015, 613143. [Google Scholar] [CrossRef] [Green Version]
  53. Manzello, S.L.; Suzuki, S.; Gollner, M.; Carlos Fernandez-Pello, A. Role of firebrand combustion in large outdoor fire spread. Prog. Energy Combust. Sci. 2020, 76, 100801. [Google Scholar] [CrossRef] [PubMed]
  54. Manzello, S.L.; Suzuki, S. Influence of board spacing on mitigating wood decking assembly ignition. Fire Saf. J. 2019, 110, 102913. [Google Scholar] [CrossRef]
  55. Suzuki, S.; Manzello, S.L. Investigating Coupled Effect of Radiative Heat Flux and Firebrand Showers on Ignition of Fuel Beds. Fire Technol. 2020, 57, 683–697. [Google Scholar] [CrossRef]
  56. Manzello, S.; Suzuki, S. Exposing Decking Assemblies to Continuous Wind-Driven Firebrand Showers. Fire Saf. Sci. 2014, 11, 1339–1352. [Google Scholar] [CrossRef] [Green Version]
  57. Ozicik, M.N. Radiative Transfer and Interactions with Conduction and Convection; J. Wiley: Hoboken, NJ, USA, 1973. [Google Scholar]
  58. Grishin, A.M. Mathematical Modeling of Forest Fires and New Methods of Fighting Them; Albini, F.A., Ed.; Tomsk Publishing House: Tomsk, Russia, 1997. [Google Scholar]
  59. Grishin, A.M.; Matvienko, O.V. Numerical Investigation of the Formation of a Convective Column and a Fire Tornado by Forest Fires. J. Eng. Phys. Thermophys. 2014, 87, 1080–1093. [Google Scholar] [CrossRef]
  60. Patankar, S. Numerical Methods of Solving Problems on Heat and Mass Exchange and Fluid Dynamics; Énergoatomizdat: Moscow, Russia, 1983. [Google Scholar]
  61. Arkhipov, V.A.; Matvienko, O.V.; Trofimov, V.F. Combustion of sprayed liquid fuel in a swirling flow. Combust. Explos. Shock. Waves 2005, 41, 140–150. [Google Scholar] [CrossRef]
  62. Ushakov, V.M.; Matvienko, O.V. Numerical Investigation of the Heat Exchange and Firing of Reactive Channel Walls by a High-Temperature Swirling-Gas glow. J. Eng. Phys. Thermophys. 2005, 78, 541–547. [Google Scholar] [CrossRef]
  63. Arkhipov, V.A.; Egorov, A.G.; Ivanin, S.V.; Maslov, E.A.; Matvienko, O.V. Numerical simulation of aerodynamics and combustion of a gas mixture in a channel with sudden expansion. Combust. Explos. Shock Waves 2010, 46, 647–655. [Google Scholar] [CrossRef]
  64. Matvienko, O.V. Mathematical Modeling of the Heat Transfer and Conditions of Ignition of a Turbulent Flow of a Reactive Gas. J. Eng. Phys. Thermophys. 2016, 89, 212–220. [Google Scholar] [CrossRef]
  65. Egorov, A.G.; Tizilov, A.S.; Niyazov, V.Y.; Arkhipov, V.A.; Matvienko, O.V. Effect of the swirl of cocurrent high-velocity wind speed on the geometry of an aluminum-air flame. Russ. J. Phys. Chem. 2014, 8, 712–715. [Google Scholar] [CrossRef]
  66. Westbrook, C.K.; Dryer, F.L. Chemical kinetic modeling of hydrocarbon combustion. Prog. Energy Combust. Sci. 1984, 10, 1–57. [Google Scholar] [CrossRef]
  67. Loboda, E.; Matvienko, O.; Vavilov, V.; Reyno, V. Infrared thermographic evaluation of flame turbulence scale. Infrared Phys. Technol. 2015, 72, 1–7. [Google Scholar] [CrossRef]
  68. Matvienko, O.V.; Baigulova, A.I.; Bubenchikov, A.M. Mathematical Modeling of Catalytic Oxidation of Methane in a Channel with a Porous Insert. J. Eng. Phys. Thermophys. 2014, 87, 1298–1312. [Google Scholar] [CrossRef]
  69. Grishin, A.M.; Matvienko, O.V.; Rudi, Y.A. Mathematical modeling of gas combustion in a twisted jet and of the formation of a fiery whirlwind. J. Eng. Phys. Thermophys. 2009, 82, 906–913. [Google Scholar] [CrossRef]
  70. Loitsyanskii, L.G. Mekhanika Zhidkosti i Gaza [Fluid and Gas Mechanics]; Nauka Publ.: Moskow, Russia, 1987; 840p. [Google Scholar]
  71. Gebhart, B.; Jaluria, Y.; Mahajan, R.L.; Sammakia, B. Free-Convective Flows, Heat and Mass Transfer; MIR: Moscow, Russia, 1991; Volume 1, 678p. [Google Scholar]
  72. Kasymov, D.; Agafontsev, M.; Perminov, V.; Martynov, P.; Reyno, V.; Loboda, E. Experimental Investigation of the Effect of Heat Flux on the Fire Behavior of Engineered Wood Samples. Fire 2020, 3, 61. [Google Scholar] [CrossRef]
Figure 1. Process diagram.
Figure 1. Process diagram.
Fire 06 00148 g001
Figure 2. Scheme of control points.
Figure 2. Scheme of control points.
Fire 06 00148 g002
Figure 3. Temperature change over time: Δ t = 0, d = 4.5 mm, h = 40 mm, 1—L = 140 mm, 2—L = 100 mm, 3—L = 84 mm, 4—L = 76 mm, 5—L = 68 mm, 6—L = 62 mm, 7—L = 58 mm.
Figure 3. Temperature change over time: Δ t = 0, d = 4.5 mm, h = 40 mm, 1—L = 140 mm, 2—L = 100 mm, 3—L = 84 mm, 4—L = 76 mm, 5—L = 68 mm, 6—L = 62 mm, 7—L = 58 mm.
Fire 06 00148 g003
Figure 4. Temperature change over time: Δ t = 0, d = 15 mm, h = 40 mm, 1—L = 59 mm, 2—L = 49 mm, 3—L = 39 mm, 4—L = 29 mm, 5—L = 26 mm, 6—L = 22 mm, 7—L = 18 mm, 8—L = 14 mm.
Figure 4. Temperature change over time: Δ t = 0, d = 15 mm, h = 40 mm, 1—L = 59 mm, 2—L = 49 mm, 3—L = 39 mm, 4—L = 29 mm, 5—L = 26 mm, 6—L = 22 mm, 7—L = 18 mm, 8—L = 14 mm.
Fire 06 00148 g004
Figure 5. Temperature change over time: Δ t = 0, d = 15 mm, L = 22 mm, h = 40 mm, 1—h = 32 mm, 2—h = 45 mm, 3—h = 59 mm, 4—h = 12 mm, 5—h = 26 mm, 6—h = 55 mm.
Figure 5. Temperature change over time: Δ t = 0, d = 15 mm, L = 22 mm, h = 40 mm, 1—h = 32 mm, 2—h = 45 mm, 3—h = 59 mm, 4—h = 12 mm, 5—h = 26 mm, 6—h = 55 mm.
Fire 06 00148 g005
Figure 6. Isotherms on the wood layer surface: at different moment L = 69 mm, d = 7.5 mm, h = 16.6 mm, Δ t = 0.
Figure 6. Isotherms on the wood layer surface: at different moment L = 69 mm, d = 7.5 mm, h = 16.6 mm, Δ t = 0.
Fire 06 00148 g006aFire 06 00148 g006b
Figure 7. Isotherms on the wood layer surface at difference moment: L = 69 mm, d = 6.5 mm, h = 8.5 mm.
Figure 7. Isotherms on the wood layer surface at difference moment: L = 69 mm, d = 6.5 mm, h = 8.5 mm.
Fire 06 00148 g007aFire 06 00148 g007bFire 06 00148 g007c
Figure 8. Temperature change over time: Δ t = 2 s, d = 5.5 mm, h = 6.5 mm, 1 —L = 19 mm, 2 – L = 29, 3– L = 39, 4 – L = 49, 5 – L = 59, 6 – L = 69.
Figure 8. Temperature change over time: Δ t = 2 s, d = 5.5 mm, h = 6.5 mm, 1 —L = 19 mm, 2 – L = 29, 3– L = 39, 4 – L = 49, 5 – L = 59, 6 – L = 69.
Fire 06 00148 g008aFire 06 00148 g008b
Figure 9. Temperature change over time: Δ t = 3 s, h = 6.5 mm, 1—L = 69 mm, d = 3.5 mm, 2—d = 4.5 mm, 3—d = 5.5 mm, 4—d = 6.5 mm, 5– d = 7.5 mm.
Figure 9. Temperature change over time: Δ t = 3 s, h = 6.5 mm, 1—L = 69 mm, d = 3.5 mm, 2—d = 4.5 mm, 3—d = 5.5 mm, 4—d = 6.5 mm, 5– d = 7.5 mm.
Fire 06 00148 g009aFire 06 00148 g009b
Figure 10. Temperature change over time: Δ t = 4 s, L= 69 mm, d = 6.5 mm: 1—h = 4.5 mm, 2—h = 6.5 mm, 3—h = 8.5 mm, 4—h = d = 10.5 mm, 5—h = 12.5 mm, 6—h = 14.5 mm, 7—16.5 mm, 8—18.5 mm, 9—h = 20.5 mm, 10—h = 22.5 mm, 11—h = 24.5 mm –, 12—h = 26.5 mm, 13—h = 28.5 mm, 14—h = 30.5 mm, 15—h = 32.5 mm, 16—h = 34.5 mm, 17—36.5 mm, 18—38.5 mm, 19—h=40.5 mm, 20–42.5 mm.
Figure 10. Temperature change over time: Δ t = 4 s, L= 69 mm, d = 6.5 mm: 1—h = 4.5 mm, 2—h = 6.5 mm, 3—h = 8.5 mm, 4—h = d = 10.5 mm, 5—h = 12.5 mm, 6—h = 14.5 mm, 7—16.5 mm, 8—18.5 mm, 9—h = 20.5 mm, 10—h = 22.5 mm, 11—h = 24.5 mm –, 12—h = 26.5 mm, 13—h = 28.5 mm, 14—h = 30.5 mm, 15—h = 32.5 mm, 16—h = 34.5 mm, 17—36.5 mm, 18—38.5 mm, 19—h=40.5 mm, 20–42.5 mm.
Fire 06 00148 g010aFire 06 00148 g010b
Figure 11. Temperature change over time: L = 49 mm, d = 5.5 mm, h = 6.5 mm, 1— Δ t = 0, 2—0.5 s, 3—1 s, 4—2 s, 5—3 s, 6—4 s, 7—5 s, 8—6 s, 9—one firebrand.
Figure 11. Temperature change over time: L = 49 mm, d = 5.5 mm, h = 6.5 mm, 1— Δ t = 0, 2—0.5 s, 3—1 s, 4—2 s, 5—3 s, 6—4 s, 7—5 s, 8—6 s, 9—one firebrand.
Fire 06 00148 g011aFire 06 00148 g011b
Figure 12. Dependence of the critical delay time on the size of the firebrands and the distance between them.
Figure 12. Dependence of the critical delay time on the size of the firebrands and the distance between them.
Fire 06 00148 g012
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Matvienko, O.; Kasymov, D.; Loboda, E.; Lutsenko, A.; Daneyko, O. Simulation of the Impact of Firebrands on the Process of the Wood Layer Ignition. Fire 2023, 6, 148. https://doi.org/10.3390/fire6040148

AMA Style

Matvienko O, Kasymov D, Loboda E, Lutsenko A, Daneyko O. Simulation of the Impact of Firebrands on the Process of the Wood Layer Ignition. Fire. 2023; 6(4):148. https://doi.org/10.3390/fire6040148

Chicago/Turabian Style

Matvienko, Oleg, Denis Kasymov, Egor Loboda, Anastasia Lutsenko, and Olga Daneyko. 2023. "Simulation of the Impact of Firebrands on the Process of the Wood Layer Ignition" Fire 6, no. 4: 148. https://doi.org/10.3390/fire6040148

APA Style

Matvienko, O., Kasymov, D., Loboda, E., Lutsenko, A., & Daneyko, O. (2023). Simulation of the Impact of Firebrands on the Process of the Wood Layer Ignition. Fire, 6(4), 148. https://doi.org/10.3390/fire6040148

Article Metrics

Back to TopTop