# A Review of Fire Interactions and Mass Fires

Authors: Mark A. Finney; Sara S. McAllister
Journal: Journal of Combustion
 (2011)

---

## Abstract

The character of a wildland fire can change dramatically in the presence of another nearby fire. Understanding and predicting the changes in behavior due to fire-fire interactions cannot only be life-saving to those on the ground, but  also be used to better control a prescribed fire to meet objectives. In discontinuous fuel types, such interactions may elicit fire spread where none otherwise existed. Fire-fire interactions occur naturally when spot fires start ahead of the main fire and when separate fire events converge in one location. Interactions can be created intentionally during prescribed fires by using spatial ignition patterns. Mass fires are among the most extreme examples of interactive behavior. This paper presents a review of the detailed effects of fire-fire interaction in terms of merging or coalescence criteria, burning rates, flame dimensions, flame temperature, indraft velocity, pulsation, and convection column dynamics. Though relevant in many situations, these changes in fire behavior have yet to be included in any operational-fire models or decision support systems.


---

## Body

## 1. Introduction

Some interactions of wildland fires are experienced routinely under field conditions. Firefighters and prescribed fire personnel see flames tilting towards adjacent ignition points or fire edges, particularly as the sources advance closer together [1, 2]. In the extreme case, interactions occurring when large areas are ignited and burning simultaneously are described as mass fires, area fires, or “fire storms” [3]. Hundreds or thousands of individual fires may interact over an area and exhibit some “unified” behavior. Such fires are generally described as having such strong indrafts that outward propagation is minimal, extremely tall convection columns or smoke plumes, and burn for long durations until all the fuel within its perimeter is consumed. Good reviews of mass or large area fires can be found in [4–6]. Mass fires were responsible for tremendous burning rates and tornado-strength winds [7] witnessed after the fire bombings of cities in Germany and Japan during World War II [8, 9] and have been studied mainly in relation to consequences of nuclear attacks [3, 10–34]. Many of these studies were through “Project Flambeau,” a joint effort between the U.S. Office of Civil Defense-Defense Atomic Support Agency and the U.S. Department of Agriculture—Forest Service in the mid-1960s. These fires were designed to mimic a suburb fire. Each square fuel bed was constructed with a mixture of pinyon pine and juniper and was approximately the same size and fuel load as a typical suburban house (185.8 m2 and about 18,000 kg of fuel). The spacing between fuel beds was either 7.6 m or 35.1 m and fire sizes were 2, 6, 12, and 20 hectares. Airflow velocities and temperatures were measured inside and just outside the fire area along with thermal radiation just outside the fire area, oxygen and carbon dioxide concentrations inside the fire area, and the mass loss rate of the fuel beds [3, 12, 13].Wildland fire interactions are intentionally manipulated for ignition or firing operations (see Figures1 and 2) to orient spread directions [35], use indrafts for backfire operations [36], increase the development of convection columns on prescribed fires through center-firing techniques [1], and limit spread and intensity with spot fire ignitions [37–39]. Rapid increases in fire growth and energy release—termed “blowup”—are sometimes associated with fire interactions [40]. Yet, despite the common usage and practical familiarity with interactions that fire personnel often acquire, there is very little quantitative physical understanding of these behaviors and no operational models that can predict them. By comparison to other fire behavior characteristics, such as fire spread rates, fire interactions at any scale have been subject to limited study.Figure 1
Prescribed fire ignition patterns designed to restrict or enhance fire front interactions [1].Figure 2
Indrafting and flame response of sequential line fires in prescribed burning [2].In this paper, we endeavored to obtain literature from many sources, including wildland fire, structural fire, as well as combustion engineering and fluid dynamics, in order to cover the range of research on fire-fire interactions and the state of knowledge. Our search revealed that the topic of fire interactions overlaps considerably with other fire behaviors that are distinguished individually, such as vortices and terrain effects. These behaviors will be mentioned when appropriate, but their full discussion is beyond the scope of this paper.


## 2. Background: Time-Dependent Fire Behaviors

For a constant set of environmental conditions, fire behavior is known to change with time. These changes are not expressly considered interactions, but spread and intensity changes within individual fires are also affected during interaction among fires, and may contribute to later development of interactions. Thus, such behaviors provide useful background material for discussion of fire-fire interactions, although studies of fire acceleration have not directly addressed interactions of multiple fires. Many of the time-dependent changes in fire behavior are associated with fire growth or expansion in two dimensions. Changes are observed in spread rates (acceleration), frontal geometry (width, curvature), and heat transfer indicated by the orientation and size of flames. These fire characteristics are interrelated with spread processes, and the literature does not discern the causes of observable features from their probable effects.


### 2.1. Fire Acceleration

Fire acceleration is defined as the time-dependent changes in spread and intensity occurring under constant weather and uniform fuel conditions. The notion of acceleration is implicitly applied to fires that are already capable of spreading as compared to combinations of threshold conditions where spread only occurs above some limit. Various mathematical representations of acceleration (Figure3) have been proposed from a theoretical standpoint that express spread rate from a point-source fire as a negative exponential function of time [41, 42]. Parameters of these equations were fit to empirical data from wind tunnel experiments by [43]. These functions asymptote toward a final equilibrium rate and are, thus, commonly communicated in terms of the time to reach some fixed fraction of equilibrium (e.g., 90%). A similar result was developed by Weber [44] which represented acceleration of fires expanding as a circle from a point ignition and depended on the curvature of the fire front.Figure 3
Theoretical fire spread rate acceleration curves from point-source ignitions show asymptotic increase in spread rate over time toward an equilibrium rate (from [43]).Studies of acceleration typically report time elapsed from ignition to a near-steady spread rate. Values of 20–30 minutes for point-source ignitions in slash fuels for prescribed fire conditions [45] and in pine litter and feather moss [46] have been reported. Wind-driven grass fires in Australia [47] showed large variation in acceleration times (about 6 minutes under slow wind conditions to over 45 minutes with faster winds) and a strong dependency on the width of the fire front. Wind tunnel burns of shallow (8 cm deep) pine needle and excelsior beds suggested time to equilibrium of only a few minutes [43] and largely independent of windspeed. Data from point-ignitions in pine needle litter reported by Curry and Fons [48] suggested wind speed affected acceleration rate (increased time to equilibrium) as well as a final spread rate. Wind speed may also affect acceleration times for conflagrations involving structures at urban densities. Chandler et al. [49] referenced much longer time estimates than wildland fuels, including 1 hour to achieve near-steady spread rates with wind speed up to 6.7 m/s, 2 hrs for winds to 17.9 m/s and possibly much longer times for stronger winds. A long acceleration period, exceeding the 36 minute observation time, was described for line ignitions in heavy fuel loadings associated with felled eucalyptus slash [50]. By contrast, rapid acceleration to near-steady burning after line ignition was reported for experimental crown fires in jack pine forests [51]. Implications of a theoretical analysis by Albini [47] suggests that line ignitions in surface fuels could accelerate very rapidly, initially overshooting the steady rate, but slow and exhibit damped oscillations toward the steady value as the increasing vertical buoyancy of the combustion zone offsets horizontal wind force. From the existing literature, it is not clear what influences the various factors of fuel loading, fuel sizes or burning duration, and final spread rates have on acceleration time, nor more complicated interactions among multiple flame zones or heat sources.It should be mentioned that acceleration of fires can also occur when air inflow is asymmetrically restricted by surface topography, either in canyons [52] or inclined channels [53] and slopes [54, 55]. Detailed treatment of these important fire-topographic interactions, however, is beyond the scope of this paper of fire-fire interactions.


### 2.2. Length of Fire Front

Fire acceleration and final spread rate appear to be dependent on fire size. Fires accelerate slowly from point-ignition sources [43, 45, 56] relative to line-source ignitions [35, 56]. At the small scale of laboratory stick arrays, fuel bed width and proportion of edge on the curvature of the head fire had significant effects on spread rate [57]. In wind driven grass fires, fire spread rates were found to be dependent on the length of the ignition for lines shorter than 50–75 m [56] and required longer acceleration times for higher winds (Figure 4). Experiments and modeling by Wotton et al. [58] for fires in red pine litter, however, showed no increase in radiation from flames for ignition lines longer than about 2 m and no effect of line width on spread rate beyond about 1 m. Dold et al. [59] offered an explanation for fire size effect on forward spread rate. As fires expand in two dimensions, the distance between the fire edges increases, meaning that buoyancy-induced inflow along segments of flaming front comes from a wider area. This allows ambient winds from behind the front to penetrate to the heading portion of the flame zone. Such effects on narrow combustion zones of expanding fires are presumably different than for mass-fires or large-area ignitions which create indrafts from all directions [60, 61] and strong buoyancy-driven convection may deflect ambient air-flow around the column [3].Figure 4
Fire spread rates in grass fuels were found to increase with the width of head fires and depend on the final spread rates determined by wind speed [62].


### 2.3. Flame Tilt

Flame angle orientation relative to the unburned fuel is related to acceleration and is affected by fire size and stage of growth. Flames can tilt due to wind, slope, or the interaction with other fires. Flames tilted away from the direction of spread are referred to as backing fires, and flames tilting toward the direction of spread are referred to as heading fires. Flames tilt toward the interior of the burned area in small fires or point-source fires, producing backing spread [37, 57, 62]. Spread rate of backing fires spreading downslope has been shown to be only weakly diminished as slope increases [63] and little affected by wind [43, 64]. Backing fires have been reported to increase fuel consumption and residence times. As fires grow larger, backing fire remains only at the rear of the perimeter (upwind or downslope) and flames for the heading portion of the fire tilt toward the unburned fuel. The very large differences in spread rate and intensity between backing and heading fires (and flanking fires) can be estimated assuming elliptical fire shapes [65]. Numerous studies of flame tilt angle in wildland fuel bed on flat terrain in wind have consistently found a strong relationship to the Froude number calculated from ratios of wind speed to intensity or flame length [47, 66–68]. Similar experimental results were found using liquid pool fires [69–72] and explained as the counteraction of upward buoyant forces by cross flow, including flame trailing  (lateral deflection of combustion products and flames) with high wind speeds. Recent numerical modeling [73] has also reported Froude number relationships for both line-source and point-source simulated fires. Although slope effects were deemed significant [68], they are not accounted for in such formulations. When fires are in close proximity, the interaction between them can change the flame tilt angle and rates of spread [5, 74, 75]. In these cases, the flame tilt angles can be correlated with a modified Froude number that includes the separation distance of the fires [5, 74, 75]. In the case of no wind, a modified Grashof number is used [76, 77] to describe the flame tilt purely due to flame interaction.


### 2.4. Spread Thresholds

Thresholds describe a point of near-instantaneous acceleration that delineates when fire will and will not spread. Threshold-crossing for fire spread has been documented for many discontinuous fuel types including grasses [78], shrubs [79–82], and trees [83, 84]. Laboratory scale fires reveal similar spread thresholds in arrays of small sticks [85–87] and taller beds of excelsior [88].  These studies reveal threshold dependencies on multiple environmental, fuel, and fire variables, such as wind speed, fuel moisture, slope, horizontal fuel gap dimensions, fuel bed depth, fuel combustion rate, and flame size. Chandler et al. [49] proposed combinations of ranges of wind speed, humidity, and rainfall by fuel type to define spread thresholds for significant growth of large fires. Recent studies of fire spread sustainability provide empirical evidence on the importance of fuel moisture, wind, and fuel loading [89, 90]. As described in later sections of this chapter, fire-interactions exert strong influences over many of these same environmental and fire variables, and thus, may elicit threshold-crossing spread for fires burning in discontinuous fuels.


## 2.1. Fire Acceleration

Fire acceleration is defined as the time-dependent changes in spread and intensity occurring under constant weather and uniform fuel conditions. The notion of acceleration is implicitly applied to fires that are already capable of spreading as compared to combinations of threshold conditions where spread only occurs above some limit. Various mathematical representations of acceleration (Figure3) have been proposed from a theoretical standpoint that express spread rate from a point-source fire as a negative exponential function of time [41, 42]. Parameters of these equations were fit to empirical data from wind tunnel experiments by [43]. These functions asymptote toward a final equilibrium rate and are, thus, commonly communicated in terms of the time to reach some fixed fraction of equilibrium (e.g., 90%). A similar result was developed by Weber [44] which represented acceleration of fires expanding as a circle from a point ignition and depended on the curvature of the fire front.Figure 3
Theoretical fire spread rate acceleration curves from point-source ignitions show asymptotic increase in spread rate over time toward an equilibrium rate (from [43]).Studies of acceleration typically report time elapsed from ignition to a near-steady spread rate. Values of 20–30 minutes for point-source ignitions in slash fuels for prescribed fire conditions [45] and in pine litter and feather moss [46] have been reported. Wind-driven grass fires in Australia [47] showed large variation in acceleration times (about 6 minutes under slow wind conditions to over 45 minutes with faster winds) and a strong dependency on the width of the fire front. Wind tunnel burns of shallow (8 cm deep) pine needle and excelsior beds suggested time to equilibrium of only a few minutes [43] and largely independent of windspeed. Data from point-ignitions in pine needle litter reported by Curry and Fons [48] suggested wind speed affected acceleration rate (increased time to equilibrium) as well as a final spread rate. Wind speed may also affect acceleration times for conflagrations involving structures at urban densities. Chandler et al. [49] referenced much longer time estimates than wildland fuels, including 1 hour to achieve near-steady spread rates with wind speed up to 6.7 m/s, 2 hrs for winds to 17.9 m/s and possibly much longer times for stronger winds. A long acceleration period, exceeding the 36 minute observation time, was described for line ignitions in heavy fuel loadings associated with felled eucalyptus slash [50]. By contrast, rapid acceleration to near-steady burning after line ignition was reported for experimental crown fires in jack pine forests [51]. Implications of a theoretical analysis by Albini [47] suggests that line ignitions in surface fuels could accelerate very rapidly, initially overshooting the steady rate, but slow and exhibit damped oscillations toward the steady value as the increasing vertical buoyancy of the combustion zone offsets horizontal wind force. From the existing literature, it is not clear what influences the various factors of fuel loading, fuel sizes or burning duration, and final spread rates have on acceleration time, nor more complicated interactions among multiple flame zones or heat sources.It should be mentioned that acceleration of fires can also occur when air inflow is asymmetrically restricted by surface topography, either in canyons [52] or inclined channels [53] and slopes [54, 55]. Detailed treatment of these important fire-topographic interactions, however, is beyond the scope of this paper of fire-fire interactions.


## 2.2. Length of Fire Front

Fire acceleration and final spread rate appear to be dependent on fire size. Fires accelerate slowly from point-ignition sources [43, 45, 56] relative to line-source ignitions [35, 56]. At the small scale of laboratory stick arrays, fuel bed width and proportion of edge on the curvature of the head fire had significant effects on spread rate [57]. In wind driven grass fires, fire spread rates were found to be dependent on the length of the ignition for lines shorter than 50–75 m [56] and required longer acceleration times for higher winds (Figure 4). Experiments and modeling by Wotton et al. [58] for fires in red pine litter, however, showed no increase in radiation from flames for ignition lines longer than about 2 m and no effect of line width on spread rate beyond about 1 m. Dold et al. [59] offered an explanation for fire size effect on forward spread rate. As fires expand in two dimensions, the distance between the fire edges increases, meaning that buoyancy-induced inflow along segments of flaming front comes from a wider area. This allows ambient winds from behind the front to penetrate to the heading portion of the flame zone. Such effects on narrow combustion zones of expanding fires are presumably different than for mass-fires or large-area ignitions which create indrafts from all directions [60, 61] and strong buoyancy-driven convection may deflect ambient air-flow around the column [3].Figure 4
Fire spread rates in grass fuels were found to increase with the width of head fires and depend on the final spread rates determined by wind speed [62].


## 2.3. Flame Tilt

Flame angle orientation relative to the unburned fuel is related to acceleration and is affected by fire size and stage of growth. Flames can tilt due to wind, slope, or the interaction with other fires. Flames tilted away from the direction of spread are referred to as backing fires, and flames tilting toward the direction of spread are referred to as heading fires. Flames tilt toward the interior of the burned area in small fires or point-source fires, producing backing spread [37, 57, 62]. Spread rate of backing fires spreading downslope has been shown to be only weakly diminished as slope increases [63] and little affected by wind [43, 64]. Backing fires have been reported to increase fuel consumption and residence times. As fires grow larger, backing fire remains only at the rear of the perimeter (upwind or downslope) and flames for the heading portion of the fire tilt toward the unburned fuel. The very large differences in spread rate and intensity between backing and heading fires (and flanking fires) can be estimated assuming elliptical fire shapes [65]. Numerous studies of flame tilt angle in wildland fuel bed on flat terrain in wind have consistently found a strong relationship to the Froude number calculated from ratios of wind speed to intensity or flame length [47, 66–68]. Similar experimental results were found using liquid pool fires [69–72] and explained as the counteraction of upward buoyant forces by cross flow, including flame trailing  (lateral deflection of combustion products and flames) with high wind speeds. Recent numerical modeling [73] has also reported Froude number relationships for both line-source and point-source simulated fires. Although slope effects were deemed significant [68], they are not accounted for in such formulations. When fires are in close proximity, the interaction between them can change the flame tilt angle and rates of spread [5, 74, 75]. In these cases, the flame tilt angles can be correlated with a modified Froude number that includes the separation distance of the fires [5, 74, 75]. In the case of no wind, a modified Grashof number is used [76, 77] to describe the flame tilt purely due to flame interaction.


## 2.4. Spread Thresholds

Thresholds describe a point of near-instantaneous acceleration that delineates when fire will and will not spread. Threshold-crossing for fire spread has been documented for many discontinuous fuel types including grasses [78], shrubs [79–82], and trees [83, 84]. Laboratory scale fires reveal similar spread thresholds in arrays of small sticks [85–87] and taller beds of excelsior [88].  These studies reveal threshold dependencies on multiple environmental, fuel, and fire variables, such as wind speed, fuel moisture, slope, horizontal fuel gap dimensions, fuel bed depth, fuel combustion rate, and flame size. Chandler et al. [49] proposed combinations of ranges of wind speed, humidity, and rainfall by fuel type to define spread thresholds for significant growth of large fires. Recent studies of fire spread sustainability provide empirical evidence on the importance of fuel moisture, wind, and fuel loading [89, 90]. As described in later sections of this chapter, fire-interactions exert strong influences over many of these same environmental and fire variables, and thus, may elicit threshold-crossing spread for fires burning in discontinuous fuels.


## 3. Conditions where Fire Interactions Occur

Interactions are possible when many separate fires grow together or multiple segments of a single continuous fire are oriented in close proximity. In natural wildland fires, multiple fronts often occur because of spotting from a single main fire. Spot fires are relatively common under dry and windy conditions and even long-distance spotting contributes to fire movement [91]. But, massive deposition of firebrands at relatively short distances from the fire front (a few kilometers) can substantially increase spread rate and create simultaneous area ignition [41]. On wildfires Cheney and Bary observed the highest concentration of fire brands fell within a fan-shaped zone about 9 degrees in angle on either side of the primary wind direction and theorized that mass fire behavior could be achieved for certain unspecified combinations of fire brand density and acceleration time for individual ignitions. Johansen [38] made similar observations for spot ignition patterns on prescribed burns where higher spot densities increased the numbers and frequencies of junction or merger zones. The increase in intensity at such junction zones have been documented empirically [38, 92] and modeled [93] leading to recommendations for wide separation of ignitions [62, 94] unless area ignition is desired [95]. Mass ember deposition and area ignition has been documented by McArthur [96] for Tasmanian fires which resulted in near-simultaneous ignition of hillsides. A similar process was proposed for the Air Force Bomb Range Fire [97] which periodically caused area ignition ahead of the main front and vertical development of a convection column. Modeling by Weihs and Small [98] showed that interactions between large mass fires can even cause these typically nonspreading fires to propagate toward one another.How close together fires must be before flames visibly interact and subsequently merge is not clear. There have been many empirically derived merging criteria in the literature. Correlations exist for the critical parameters for both flame interaction [99–101] and merging [24, 102–104]. These correlations take many forms—some define a critical ratio between the fire spacing and fire diameter [24, 100] or flame height [99, 101, 104], some define a critical ratio between the flame height and fire diameter [24], and some define a critical dimensionless heat release rate [102, 103]. Upon close examination, however, it becomes clear that fire spacing, fire diameter, flame height, and dimensionless heat release rate have interdependencies and, thus, these different correlations are not necessarily contradictory. The discussion here will focus on the relations between spacing, diameter, and flame height because they are the most intuitive.Using both gas diffusion burners and pool fires, Sugawa and Takahashi [100] report that flames begin to interact when the ratio of the spacing distance to the fire diameter is less than four. In other words, flames can interact, here defined as visually tilting, over distances four times their diameter. Baldwin [99] considered the onset of flame interaction in terms of flame height. Flames were considered to be interacting if the flame heights increased more than 10% above the independent flame height. Using square and round gas burners, wood cribs, and large timber yard fires Baldwin [99] (and [105–107]) correlated experimental data over a wide range of scales and configurations found in the literature and determined that the flames would interact if the spacing was less than 0.22 times the flame length. For a characteristic dimension D and height L, this correlation holds for 1<L/D<300. Liu et al. [101] also found the same dependency but with a slightly different constant of proportionality for merging of round pool fires. In their experiments, flame merging was likely to occur when closer than 0.29 to 0.34 times the merged flame length. Delichatsios [104] also found that flames began to merge at spacing less than 0.33 times the actual flame length for gaseous burners. The discrepancy in these constants may be due to different definitions of flame interaction (tilting versus change in flame height) and flame merging (using completely merged flame height versus actual flame height), different fuels, and possibly uncertainty of measuring flame dimensions. In comparing the results of the Project Flambeau fires to those using a sand-filled pan burner, Wood et al. [24] reported that flames merged if the flame height was at least half of the fire diameter. Heskestad [6] clarified that this occurs when the nondimensional group N~Q2/D5 is near 10−5 (Q is the heat release rate and D is the fire diameter). Clearly there is no definitive criterion for when flames begin to interact and merge, and these relations will remain qualitative guidelines until there is some sort of unifying theory.An opposing effect may occur with area fires over large homogenous fuel beds (small flame height compared to fire diameter). For a sufficiently large fuel bed, it may be impossible for a continuous flame to exist over the entire bed. Instead of one continuous flame, the fire may break up into many distributed flamelets [13, 24, 108]. Heskestad [108] showed that the breakup of continuous flames occurs when the nondimensional group N~Q2/D5 is near 10−6. The convection column for these cases has been described as having two modes: Bénard cell convection near the surface which merges and transitions to a more organized convective plume [109].


## 4. Specific Effects of Fire Interaction

Studies of fire interactions involve specific types of behavior of the combustion and observable fire characteristics. Most of the research on these behaviors comes from laboratory experiments with artificial fuel sources and attempt to isolate the particular response of interest.


### 4.1. Burning Rate

When fire fronts are close enough to interact and merge, such as in a mass fire, the mass of fuel burned as a function of time, or burning rate, of the fire can change dramatically. Though much of the research on fire interactions has been done using gas burners with a fixed burning rate, there has been some work on the interaction of flames over liquid pool fires and wood crib fires. Though the geometry and heat transfer mechanisms inside the fuel bed are different, liquid pool fires are much like fires burning over solid fuel in that the heat transfer from the fire back to the fuel controls the burning rate. In contrast, the burning rate of a gas burner is controlled by using a fixed fuel supply rate. Results from pool and crib fire experiments can often be extended to larger fuel beds using appropriate scaling laws [110].The experiments by Huffman et al. [111] clearly reveal the effect of spacing on the burning rate of pool fires. In this work, the burning rate of an array of liquid pools was measured while keeping a constant fuel depth and varying the number of pools, pool diameter, fuel, and pool separation distance. In general, the burning rate of each individual pool burner increases as the burners are brought closer together and the flames began to interact. In particular, the pools in the middle of the array show a very dramatic increase. For example, Figure 5 shows that the burning rate of 4 inch diameter pools of cyclohexane experienced over a 400% increase in burning rate when the separation distance was halved. At the onset of flame merging, the burning rate is at its maximum. As the flames merge, the burning rate decreases as the separation distance continues to decrease. In the limit of zero separation distance, however, the burning rate of the individual fires is still larger than if they were burning independently with no interaction effects. These trends were also seen by Grumer and Strasser [112] with solid fuel beds.Figure 5
Burning rate as a function of separation distance for 10.1 cm diameter cyclohexane burners from [111].Kamikawa et al. [113] studied the effect of flame merging on heat release rates (heat released per time). Heat release rate is calculated by multiplying the burning rate (mass of fuel burned per time) by the heat of reaction (heat released per mass of fuel burned). However, the heat of reaction is dependent on the fuel and the mixture ratio of fuel to air. In large fire arrays, the inner regions of the array typically experience a shortage of air. Without sufficient air, the fuel cannot completely react and release the full potential heat, that is, the combustion efficiency is low and less heat is released per mass of fuel. Not surprisingly, Kamikawa et al. [113] saw the same trend with heat release rates as Huffman et al. [111] with burning rates. When the flames are merged, the heat release rate increases with separation distance. As the burners are moved further apart, more air can penetrate into the inner regions of the array. More air entrainment means greater combustion efficiency and greater heat release. This in turn heats up and evaporates the unburned fuel more quickly, increasing the burning rate.Liu et al. [114] explain the mechanisms behind these trends in burning and heat release rate with separation distance. The nonmonotonic behavior seen in Figure 5 is the result of two competing mechanisms: heat feedback enhancement and air entrainment restriction. As the burners are moved closer, the view factor between neighboring fires increases. In other words, the fires can “see” each other better, increasing the radiative heat transfer in addition to the convective heat transfer [112]. Because the burning rate is dictated by the heat feedback from the flame, this increased radiative heat seen by the fuel will evaporate the fuel more quickly and increase the burning rate. Conversely, as the fires get sufficiently close there is less room to entrain air inside the array and the flames become “choked.” When the flames are merely interacting, the heat feedback mechanism is more important than the air restriction and the burning rate increases. When the flames have merged, the air restriction is the dominant mechanism and the burning rate decreases.Because the experiments by Kamikawa et al. [113] used wood crib fires, they were also able to examine the release rate as a function of time for merged flames. As with most wildland fires, the heat release rate (and burning rate) of wood crib fires increases as the fire builds, reaches a maximum, then begins to decrease as the fuel is depleted. Kamikawa et al. [113] made the observation that as the number of fires increases, the peak heat release rate increases above that expected by multiplying the independent fire heat release rate by the number of fires. This discrepancy grows as the number of fires increases. So the burning and heat release rates of interacting and merging fires not only are dependent on the spacing of the fires, but also on the total number of fires (see also [114]).Fire interactions can increase burning rates by another mode as well. If the fires interact such that vorticity is generated, fire whirls can form. Though not discussed further here, it has been shown that fire whirls have dramatically increased burning rates in comparison to an equivalent, nonrotating fire (see, e.g., [115, 116]).


### 4.2. Flame Dimensions

Flame height trends for a nonpremixed flame, such as those in a wildfire, are usually discussed in terms of two dimensionless parameters: the dimensionless flame height and the dimensionless heat release rate. The dimensionless flame height is usually defined as the flame height divided by the characteristic burning area diameter (D). The characteristic burning area diameter is a dimensioned parameter frequently introduced in fire arrays and is usually some function of the number of fires, fire diameter, and the fire arrangement (separation distance). The dimensionless heat release rate (Q*) is usually defined as the total heat release rate of the group divided by the characteristic burning area diameter to the five-halves power (material property constants are used to make the ratio dimensionless: Q*~Qtot/D5/2). The dimensionless heat release rate for natural fires tends to fall between 0.05 and 5 [117].Much of the research on flame height has been performed using gas burners. However, two regimes of flow from a gas burner can be identified. When the flow velocity is low or the burner diameter is large, the momentum of the gaseous fuel is due primarily to its buoyancy.  When the flow velocity is high or the burner diameter is small, the flow is like a jet. Putnam and Speich [102] have a method for determining whether the flow from a gas burner is a high momentum jet or buoyancy controlled. The discussion here will be limited to turbulent, buoyancy-driven flames, as this situation better describes what occurs during a wildfire.In general, the flame height increases as the fires are moved closer. When the flames begin to merge, the flame height will dramatically increase with further changes in separation distance. However, once the flames are fully merged, further decreases in separation distance will have little effect [102, 103, 119]. The dimensionless flame height has successfully been correlated to the dimensionless heat release rate raised to some power a. Because the dimensionless heat release rate can vary over at least seven orders of magnitude, this power “a” can take on three different values depending on the range of the dimensionless heat release rate. As shown in Figure 6 [118], the dimensionless flame height increases with the dimensionless heat release rate. These correlations were originally developed for the flame height of a single independent burner where the characteristic dimension is the burner diameter, and hold for buoyancy-driven gas burners, liquid pool fires, and wood crib fires. However, there is an indication that these correlations also apply to interacting flames when the characteristic burning area dimension is given as discussed above. For example, for the interaction of relatively tall flames compared to the actual burner diameter (Lf/D>1, or high values of Q*), Putnam and Speich [102] and Sugawa and Takahashi [100] showed that the dimensionless flame height correlates well with the dimensionless heat release rate to the two-fifths power (Lf/D~Q*2/5). Delichatsios [104] successfully correlated the dimensionless flame height to the dimensionless heat release rate to the two thirds power (Lf/D~Q*2/3) for Q* between 0.1 and 1. On the other hand, Weng et al. [120] and Kamikawa et al. [113] showed that the data for merged flame height is better correlated with the exponent “a” varying with the number of burners.Figure 6
Dimensionless flame height(Lf/D) correlations with dimensionless heat release rate (Q*) [118].With all else remaining constant, these correlations suggest that an increase in either the number of fires or the individual fire heat release rate will increase the interacting or merged flame height. Increases in the separation distance or the fire diameter will result in a decrease in the interacting or merged flame height. An interesting caveat to these correlations is that the burning rate for individual pool or crib fires is not constant, but is a function of the separation distance as discussed above. This trend is not necessarily captured in Figure6 or by Putnam and Speich [102] (gas burners), Kamikawa et al. [113], Fukuda et al. [103], or Delichatsios [104] (all fully merged flames). Also, vorticity can greatly increase flame height as well [115].This literature suggests that in a mass fire situation, as the flames grow closer together, the heat release rate and characteristic “burner” diameter should increase. The net effect is most likely an increase in the flame height. If more spot fires were ignited in the burning area, for example, the flame height would increase further. This is consistent with the observations of spot ignitions on prescribed burns [38] and mass spotting in wildfires [41]. However, for a sufficiently large area or mass fire, when the nondimensional group N~Q2/D5 is near 10−6, the fire is not expected to burn as a continuous flame but will break up into many distributed flamelets [13, 24, 108]. In this case, the flame height will be less than that predicted for a fully-merged, continuous flame but larger than that of isolated flames [121].


### 4.3. Flame Temperatures and Pollutants

As discussed in relation to flame height, as fires are moved closer together, air entrainment is blocked and the gaseous fuel must travel higher to find sufficient air for combustion. Experiments by Chigier and Apak [119] indicated that a fuel particle on its journey from the base to the tip of an interacting turbulent flame would experience delayed combustion compared to an independent flame (see Figure 7(a)). The delay means that the maximum temperature of the interacting flames would occur further from the flame base. With limited mixing of fresh air into the flame to provide cooling, the temperatures inside an interacting flame decay more slowly with height so the flame is hot over a greater portion. In addition, limited mixing of air into the flames causes the formation of more carbon monoxide inside the flame zone. This prompted Countryman [13] to speculate that the lack of oxygen in conjunction with elevated carbon monoxide could be fatal to ground personnel trapped inside the burning area.(a) Effect of nearby burners on flame temperature from [119]. DT is throat diameter, DE is exit diameter, a is separation distance, Tm is merged flame temperature, and Ts is single flame temperature. (b) Temperature compared to independent flame for varying axial distance along flame, number of burners, burner arrangement, and burner spacing from [119].
(a)(b)Chigier and Apak [119] also showed that the maximum temperature achieved by interacting turbulent flames is also a function of the separation distance and the number of burners (see Figure 7(b)). When the flames are close enough to interact, they lose less heat from radiation (the surroundings are at the same temperature) and by mixing with cool, fresh air. The maximum temperatures inside interacting flames therefore increase as the number of fires increases and as the burners get closer together. These increased temperatures could produce more of the smog-forming nitrogen oxide emissions [122].


### 4.4. Indraft Velocity

In typical fire situations where the flame height is relatively tall compared to the fire diameter, standard correlations exist to predict the mass of air entrained by the fire and its plume due to the velocity difference between the plume gases and the ambient air. This air entrainment causes an inflow into the fire and is generally responsible for the bending of two flames in relatively close proximity. However, the standard correlations of plume theory are valid only above the flame. Though several plume theories exist in the literature (see paper in [123]), there is general agreement that the total mass of air entrained can be estimated as proportional to the convective heat release rate (heat release rate minus radiative and other losses) raised to the one-third power and to the height above the fire source to the five-thirds power. Fires with greater heat release rate entrain more air and the total amount of air entrained increases with height above the plume. Note, however, that the velocity of the flow inside the plume decreases with height, so at some point near the top of the plume no further air is entrained (no velocity difference). Current research on the indraft caused by entrainment as related to fire interactions is focused mainly providing better quantitative predictions with CFD modeling [93, 124].However, plumes from wildfires can interact with local meteorology [125] such as wind and atmospheric conditions. Additionally, classic plume theory for entrainment rates may not hold for small ratios of the flame height to fire diameter (Lf/D). Though the exact threshold is not known, Heskestad [123] contends that the standard plume theory falls apart for Lf/D somewhere between 0.14 to 0.9. The perimeter of the plume where entrainment occurs becomes too small in relation to the volume of air inside and the slow moving entrained air will not have much effect on the momentum of the entire plume. Mass fires by definition fall into the range of flame height to fire diameter ratios where classic plume theory does not hold. The results of the Project Flambeau burns confirm that there is little entrainment into the plume core [126]. Many authors (e.g., [30, 60, 127]) also argue that the entrainment of plume theory does not account for the reported high velocity winds associated with mass fires. As discussed earlier, mass fires are characterized by such strong indrafts that the fire does little outward propagation. In their paper of the range of possible indraft velocities, Trelles and Pagni [128] showed that indraft velocities of large fires can range from about 2 to 40 m/s. In the Project Flambeau burns, Countryman [3, 12, 13] also reports complicated airflow patterns and strong downdrafts which cannot be accounted for with simple plume theory.There seem to be two main theories in the literature as to what causes the high velocity inflows. One theory, advanced by Baum and McCaffrey [61] and Carrier et al. [7] is that large-scale vorticity in conjunction with heat release is responsible. These models contend that the entire fire plume slowly rotates. Note, however, that Church et al. [129] and McRae and Flannigan [130] characterize this type of motion as one type of fire whirl. In Baum and McCaffrey’s model (also used by [128, 131]), this rotation is caused by density gradients from the high heat release, and not necessarily by any imposed swirling caused by the ambient environment. The slow rotation of such a large mass of air above the ground translates to high velocity, purely horizontal and nonrotating flow at the ground. One unique feature of the Baum and McCaffrey model is that it treats the large area fire as an ensemble of randomly distributed individual fires of varying strengths. Because of the method chosen to represent the fire, the model is only valid for heights above the fuel bed where the plumes of the individual fires have not merged. The model of Carrier et al. [7] was intended to determine how long it would take to spin up the convective column and under what conditions this occurred. Based on the fact that the fire in Hamburg, Germany took two hours to develop, they concluded that the growth of swirl, at least in this case, was most likely due to the intensification of a preexisting vortex from earlier fires and bombings. Though this contradicts the Baum and McCaffrey model, the experiments and discussion by Church et al. [129] support this argument. The spatial orientation of individual fires may cause a swirling flow due to the interaction of the indrafts to each fire [132]. Carrier et al. [7] found that large diameter plumes spin up faster and proposed a set of four criteria that must be met for a “firestorm” to develop: heat release of 106 MW over a localized area for two to three hours, a preexisting weak vortex, low ambient winds, and a nearly dry-adiabatic lapse rate over the first few kilometers of the atmosphere.Because it seems unlikely that all the criteria for spin up of a convective column will be met, another theory, advanced by Smith et al. [60] and Small et al. [30] is proposed. These authors claim that buoyancy-induced pressure gradients are responsible for the large indrafts. Smith et al. [60] used a simple two-dimensional model of a convective column over a hot area to effectively show that near the fire a dynamic pressure gradient can cause high velocity inflow. This dynamic pressure gradient is caused by a balance between hydrostatic pressure and buoyancy. Buoyancy pushes the hot gases up while atmospheric pressure pushes fresh air at the ground in toward the fire horizontally to fill the gap left by the rising gases. Smith et al. [60] also suggest that the traditional “weakly buoyant” plume theories described above may be valid for a small range of plume heights sufficiently far away from the fire and any inversion layer above. Small et al. [30] use a similar model to Smith et al. [60] but include a volume heat addition and large density and temperature gradients. Small et al. [30] also numerically match their model results of the area near the fire to the results of traditional plume theory for the region far from the fire. In both the Smith et al. [60] and Small et al. [30] models, the fire is treated as a single large heat source (Figure 8). Small et al. [30] used their model to demonstrate how the maximum indraft velocity varies with fire radius, burning rate, and fire height (Figure 8). They showed that the maximum indraft velocity at first increases but eventually levels off (to approximately 40 m/s) with increases in both the fire radius and the burning rate. On the other hand, the maximum indraft velocity appears to be linear with fire height.Model results for flow-field streamlines for three fires in close proximity (from [98]).
(a)(b)A third, yet not well-explored explanation is proposed by Carrier et al. [134]. In this work, they use classic plume theory, but assume that the fire does not burn as a single fire, but a collection of individual fires. They hypothesize that the high indraft velocities are then due to the increased fire perimeter from this “multicellular burning zone.” This hypothesis was not further developed and in later works, these authors treated the fire as a subterranean point source. Interestingly, both the Baum and McCaffrey [61] and Small et al. [30] models reasonably replicate what little experimental data is available. However, the theories differ slightly in their predictions of the distance away from the fire that these indrafts extend [5]. The model of Baum and McCaffrey [61] predicts that the high velocity indrafts will extend much further from the fire compared to the model of Small et al. [30]. Without more detailed experimental data, it is impossible to say which model is more accurately portraying the physics.


### 4.5. Pulsation

Though not an effect of flame interactions, flame pulsation (or puffing) is an interesting phenomenon that can occur in stationary fires, such as a mass fire. This pulsation typically occurs in circular or axisymmetric  fires in weak ambient wind and is periodic in nature. Flame pulsation is important to many researchers because it can have a large influence on air entrainment rates and therefore heat release rates and pollution formation [135]. Observations of this phenomenon reveal the expansion of the flame near the base of the fire as a toroidal vortex, about the size of the fire diameter. As this vortex is shed and propagates upward, the flame necks inward giving the appearance of a “mushroom” shape. Figure 9 illustrates the process with time sequence of photos. Not all circular flames pulsate, however. Using dimensional analysis, Byram and Nelson [136] attempt to describe what type of fires will pulsate. They define a dimensionless “buoyancy” number, π2=Q̇c/(gD)0.5ρcpT, where Q̇c is the rate of convective heat release per area, D is the fire diameter, and ρ, cp, and T are the density, specific heat and temperature of the ambient air. Though no quantitative values are given, they argue that a fire will not pulsate if π2 is either too small (low heat release rate relative to large fire diameter) or too large (large heat release rate relative to small fire diameter).Thermal images of flame pulsation from [133].
(a)(b)Because this puffing occurs in nonreacting helium plumes, it is actually not caused by a combustion instability, but instead is produced by a fluid dynamic instability [137]. Though there is disagreement about the actual cause of the instability [138], the vortex is generally thought to be formed because of the interaction between gravity and the density gradient between the flame and ambient [135].Most of what has been learned about the characteristics of pulsation has been learned through experiments. Cetegen and Ahmed [137] showed that the toroidal vortex forms within one fire diameter above the flame base and that the frequency of the puffing is insensitive to the fuel or the heat release rate. By plotting the available data in the literature, Cetegen and Ahmed, and later Malalasekera et al.[133], showed that the pulsation frequency is proportional to the fire diameter raised to the negative one-half power (f~D-1/2) so that large fires pulsate at a much lower frequency than small fires. Though this correlation was developed using data from fires ranging from 0.1 m to 100 m in diameter (four orders of magnitude) using gaseous, liquid and solid fuels, Baum and McCaffrey [61] suggest that it may well hold for much larger fires as well. For a large fire with a diameter on the order of 20 km, Larson et al. [27] estimate that the pulsation will occur every 20 minutes. Though not accounted for in the above correlation, Malalasekera et al. [133] show that increasing fuel flow rates also result in a small increase in puffing frequency, especially for small fire sizes. Because of this, Malalasekera et al. [133] correlated the puffing frequency in a slightly different manner using the dimensionless Strouhal number (ratio of oscillation frequency to one over the characteristic time of convection) and Froude number (ratio of inertia force to gravitational force) which retains the same dependency on fire diameter but allows for a correction due to changes in fuel flow velocity.


### 4.6. Convection Column

Mass fires are also described as having very tall convection columns, or smoke plumes with large cloud structures due to the moisture release from combustion [32]. As discussed in the section on indraft velocities, the entrainment of cold, ambient air slows the rise of the hot gases by cooling them. Additionally, the density of the ambient air itself decreases with elevation. As the hot gases rise and cool, the density difference driving their upward motion disappears. It follows then that the top of the smoke plume corresponds to the where the combustion products stop moving. As the fire diameter grows, however, the entrainment predicted by classic plume theory becomes less effective. Entrainment occurs at the perimeter of the plume and with large fire sources there is such as large core of hot gases that entrainment is less effective at slowing the rise of the combustion products [126] (Palmer 1981). Thus, it takes longer to entrain enough cold air to slow the combustion products and therefore the smoke plume becomes taller. For example, a lack of entrainment to the convection column was noted and discussed by Taylor et al. [95] on a large prescribed burn. In fact, the plume from a sufficiently large mass fire may be almost as wide as it is tall, so Brode and Small [31] and Palmer [126] contend that air entrainment is not likely to be a major influence on plume height and that it is the structure of the atmosphere itself that is the limiting factor. The plume of large mass fires is therefore more sensitive to atmospheric gradients, inversion heights, and upper atmosphere crosswinds (see also [33]). Brode and Small [31] show that the tropopause/stratosphere transition may be what actually caps the smoke plume. Note, these theories contradict the suggestion of Smith et al. [60] that the traditional plume theory holds at some intermediate height above the ground. Perhaps the scale of the fires modeled by Smith et al. [60] was not large enough to see this effect.Palmer [126] describes the interesting structure of the convection columns that formed during the Project Flambeau tests. In the first few minutes of these large-scale burns, the majority of the gaseous combustion products were contained in a “bubble” near the fire. Once the “bubble” got sufficiently hot, the associated buoyancy was enough to overcome the surface drag forces and the bubble rose. As the bubble rose, a vortex ring would form in a similar manner described above with respect to flame pulsations. Regardless of the atmospheric stability, this vortex ring would rise until it encountered a region of vertical wind shear. The vertical wind shear weakens the vortex enough for the plume to then follow the prevailing horizontal winds. Palmer [126] also notes that the “exterior form of the convection column at a particular altitude was determined by the initial vortex bubble as it passed that altitude.” Most of the plumes in these fires began to rotate as a single vertical vortex, as suggested by the Baum and McCaffrey [61] model. This rotation further inhibits entrainment, which would also prevent the use of classic plume models for mass fires [139].


### 4.7. Summary of Interaction Effects

As the individual spot fires grow together, they will begin to interact. This interaction will increase the burning rates, heat release rates, and flame height until the distance between them reaches a critical level. At the critical separation distance, the flames will begin to merge together and burn with the maximum rate and flame height. As these spot fires continue to grow together, the burning and heat release rates will finally start to decrease but remain at a much elevated level compared to the independent spot fire. The flame height is not expected to change significantly. The more spot fires, the bigger the increase in burning rate and flame height.


## 4.1. Burning Rate

When fire fronts are close enough to interact and merge, such as in a mass fire, the mass of fuel burned as a function of time, or burning rate, of the fire can change dramatically. Though much of the research on fire interactions has been done using gas burners with a fixed burning rate, there has been some work on the interaction of flames over liquid pool fires and wood crib fires. Though the geometry and heat transfer mechanisms inside the fuel bed are different, liquid pool fires are much like fires burning over solid fuel in that the heat transfer from the fire back to the fuel controls the burning rate. In contrast, the burning rate of a gas burner is controlled by using a fixed fuel supply rate. Results from pool and crib fire experiments can often be extended to larger fuel beds using appropriate scaling laws [110].The experiments by Huffman et al. [111] clearly reveal the effect of spacing on the burning rate of pool fires. In this work, the burning rate of an array of liquid pools was measured while keeping a constant fuel depth and varying the number of pools, pool diameter, fuel, and pool separation distance. In general, the burning rate of each individual pool burner increases as the burners are brought closer together and the flames began to interact. In particular, the pools in the middle of the array show a very dramatic increase. For example, Figure 5 shows that the burning rate of 4 inch diameter pools of cyclohexane experienced over a 400% increase in burning rate when the separation distance was halved. At the onset of flame merging, the burning rate is at its maximum. As the flames merge, the burning rate decreases as the separation distance continues to decrease. In the limit of zero separation distance, however, the burning rate of the individual fires is still larger than if they were burning independently with no interaction effects. These trends were also seen by Grumer and Strasser [112] with solid fuel beds.Figure 5
Burning rate as a function of separation distance for 10.1 cm diameter cyclohexane burners from [111].Kamikawa et al. [113] studied the effect of flame merging on heat release rates (heat released per time). Heat release rate is calculated by multiplying the burning rate (mass of fuel burned per time) by the heat of reaction (heat released per mass of fuel burned). However, the heat of reaction is dependent on the fuel and the mixture ratio of fuel to air. In large fire arrays, the inner regions of the array typically experience a shortage of air. Without sufficient air, the fuel cannot completely react and release the full potential heat, that is, the combustion efficiency is low and less heat is released per mass of fuel. Not surprisingly, Kamikawa et al. [113] saw the same trend with heat release rates as Huffman et al. [111] with burning rates. When the flames are merged, the heat release rate increases with separation distance. As the burners are moved further apart, more air can penetrate into the inner regions of the array. More air entrainment means greater combustion efficiency and greater heat release. This in turn heats up and evaporates the unburned fuel more quickly, increasing the burning rate.Liu et al. [114] explain the mechanisms behind these trends in burning and heat release rate with separation distance. The nonmonotonic behavior seen in Figure 5 is the result of two competing mechanisms: heat feedback enhancement and air entrainment restriction. As the burners are moved closer, the view factor between neighboring fires increases. In other words, the fires can “see” each other better, increasing the radiative heat transfer in addition to the convective heat transfer [112]. Because the burning rate is dictated by the heat feedback from the flame, this increased radiative heat seen by the fuel will evaporate the fuel more quickly and increase the burning rate. Conversely, as the fires get sufficiently close there is less room to entrain air inside the array and the flames become “choked.” When the flames are merely interacting, the heat feedback mechanism is more important than the air restriction and the burning rate increases. When the flames have merged, the air restriction is the dominant mechanism and the burning rate decreases.Because the experiments by Kamikawa et al. [113] used wood crib fires, they were also able to examine the release rate as a function of time for merged flames. As with most wildland fires, the heat release rate (and burning rate) of wood crib fires increases as the fire builds, reaches a maximum, then begins to decrease as the fuel is depleted. Kamikawa et al. [113] made the observation that as the number of fires increases, the peak heat release rate increases above that expected by multiplying the independent fire heat release rate by the number of fires. This discrepancy grows as the number of fires increases. So the burning and heat release rates of interacting and merging fires not only are dependent on the spacing of the fires, but also on the total number of fires (see also [114]).Fire interactions can increase burning rates by another mode as well. If the fires interact such that vorticity is generated, fire whirls can form. Though not discussed further here, it has been shown that fire whirls have dramatically increased burning rates in comparison to an equivalent, nonrotating fire (see, e.g., [115, 116]).


## 4.2. Flame Dimensions

Flame height trends for a nonpremixed flame, such as those in a wildfire, are usually discussed in terms of two dimensionless parameters: the dimensionless flame height and the dimensionless heat release rate. The dimensionless flame height is usually defined as the flame height divided by the characteristic burning area diameter (D). The characteristic burning area diameter is a dimensioned parameter frequently introduced in fire arrays and is usually some function of the number of fires, fire diameter, and the fire arrangement (separation distance). The dimensionless heat release rate (Q*) is usually defined as the total heat release rate of the group divided by the characteristic burning area diameter to the five-halves power (material property constants are used to make the ratio dimensionless: Q*~Qtot/D5/2). The dimensionless heat release rate for natural fires tends to fall between 0.05 and 5 [117].Much of the research on flame height has been performed using gas burners. However, two regimes of flow from a gas burner can be identified. When the flow velocity is low or the burner diameter is large, the momentum of the gaseous fuel is due primarily to its buoyancy.  When the flow velocity is high or the burner diameter is small, the flow is like a jet. Putnam and Speich [102] have a method for determining whether the flow from a gas burner is a high momentum jet or buoyancy controlled. The discussion here will be limited to turbulent, buoyancy-driven flames, as this situation better describes what occurs during a wildfire.In general, the flame height increases as the fires are moved closer. When the flames begin to merge, the flame height will dramatically increase with further changes in separation distance. However, once the flames are fully merged, further decreases in separation distance will have little effect [102, 103, 119]. The dimensionless flame height has successfully been correlated to the dimensionless heat release rate raised to some power a. Because the dimensionless heat release rate can vary over at least seven orders of magnitude, this power “a” can take on three different values depending on the range of the dimensionless heat release rate. As shown in Figure 6 [118], the dimensionless flame height increases with the dimensionless heat release rate. These correlations were originally developed for the flame height of a single independent burner where the characteristic dimension is the burner diameter, and hold for buoyancy-driven gas burners, liquid pool fires, and wood crib fires. However, there is an indication that these correlations also apply to interacting flames when the characteristic burning area dimension is given as discussed above. For example, for the interaction of relatively tall flames compared to the actual burner diameter (Lf/D>1, or high values of Q*), Putnam and Speich [102] and Sugawa and Takahashi [100] showed that the dimensionless flame height correlates well with the dimensionless heat release rate to the two-fifths power (Lf/D~Q*2/5). Delichatsios [104] successfully correlated the dimensionless flame height to the dimensionless heat release rate to the two thirds power (Lf/D~Q*2/3) for Q* between 0.1 and 1. On the other hand, Weng et al. [120] and Kamikawa et al. [113] showed that the data for merged flame height is better correlated with the exponent “a” varying with the number of burners.Figure 6
Dimensionless flame height(Lf/D) correlations with dimensionless heat release rate (Q*) [118].With all else remaining constant, these correlations suggest that an increase in either the number of fires or the individual fire heat release rate will increase the interacting or merged flame height. Increases in the separation distance or the fire diameter will result in a decrease in the interacting or merged flame height. An interesting caveat to these correlations is that the burning rate for individual pool or crib fires is not constant, but is a function of the separation distance as discussed above. This trend is not necessarily captured in Figure6 or by Putnam and Speich [102] (gas burners), Kamikawa et al. [113], Fukuda et al. [103], or Delichatsios [104] (all fully merged flames). Also, vorticity can greatly increase flame height as well [115].This literature suggests that in a mass fire situation, as the flames grow closer together, the heat release rate and characteristic “burner” diameter should increase. The net effect is most likely an increase in the flame height. If more spot fires were ignited in the burning area, for example, the flame height would increase further. This is consistent with the observations of spot ignitions on prescribed burns [38] and mass spotting in wildfires [41]. However, for a sufficiently large area or mass fire, when the nondimensional group N~Q2/D5 is near 10−6, the fire is not expected to burn as a continuous flame but will break up into many distributed flamelets [13, 24, 108]. In this case, the flame height will be less than that predicted for a fully-merged, continuous flame but larger than that of isolated flames [121].


## 4.3. Flame Temperatures and Pollutants

As discussed in relation to flame height, as fires are moved closer together, air entrainment is blocked and the gaseous fuel must travel higher to find sufficient air for combustion. Experiments by Chigier and Apak [119] indicated that a fuel particle on its journey from the base to the tip of an interacting turbulent flame would experience delayed combustion compared to an independent flame (see Figure 7(a)). The delay means that the maximum temperature of the interacting flames would occur further from the flame base. With limited mixing of fresh air into the flame to provide cooling, the temperatures inside an interacting flame decay more slowly with height so the flame is hot over a greater portion. In addition, limited mixing of air into the flames causes the formation of more carbon monoxide inside the flame zone. This prompted Countryman [13] to speculate that the lack of oxygen in conjunction with elevated carbon monoxide could be fatal to ground personnel trapped inside the burning area.(a) Effect of nearby burners on flame temperature from [119]. DT is throat diameter, DE is exit diameter, a is separation distance, Tm is merged flame temperature, and Ts is single flame temperature. (b) Temperature compared to independent flame for varying axial distance along flame, number of burners, burner arrangement, and burner spacing from [119].
(a)(b)Chigier and Apak [119] also showed that the maximum temperature achieved by interacting turbulent flames is also a function of the separation distance and the number of burners (see Figure 7(b)). When the flames are close enough to interact, they lose less heat from radiation (the surroundings are at the same temperature) and by mixing with cool, fresh air. The maximum temperatures inside interacting flames therefore increase as the number of fires increases and as the burners get closer together. These increased temperatures could produce more of the smog-forming nitrogen oxide emissions [122].


## 4.4. Indraft Velocity

In typical fire situations where the flame height is relatively tall compared to the fire diameter, standard correlations exist to predict the mass of air entrained by the fire and its plume due to the velocity difference between the plume gases and the ambient air. This air entrainment causes an inflow into the fire and is generally responsible for the bending of two flames in relatively close proximity. However, the standard correlations of plume theory are valid only above the flame. Though several plume theories exist in the literature (see paper in [123]), there is general agreement that the total mass of air entrained can be estimated as proportional to the convective heat release rate (heat release rate minus radiative and other losses) raised to the one-third power and to the height above the fire source to the five-thirds power. Fires with greater heat release rate entrain more air and the total amount of air entrained increases with height above the plume. Note, however, that the velocity of the flow inside the plume decreases with height, so at some point near the top of the plume no further air is entrained (no velocity difference). Current research on the indraft caused by entrainment as related to fire interactions is focused mainly providing better quantitative predictions with CFD modeling [93, 124].However, plumes from wildfires can interact with local meteorology [125] such as wind and atmospheric conditions. Additionally, classic plume theory for entrainment rates may not hold for small ratios of the flame height to fire diameter (Lf/D). Though the exact threshold is not known, Heskestad [123] contends that the standard plume theory falls apart for Lf/D somewhere between 0.14 to 0.9. The perimeter of the plume where entrainment occurs becomes too small in relation to the volume of air inside and the slow moving entrained air will not have much effect on the momentum of the entire plume. Mass fires by definition fall into the range of flame height to fire diameter ratios where classic plume theory does not hold. The results of the Project Flambeau burns confirm that there is little entrainment into the plume core [126]. Many authors (e.g., [30, 60, 127]) also argue that the entrainment of plume theory does not account for the reported high velocity winds associated with mass fires. As discussed earlier, mass fires are characterized by such strong indrafts that the fire does little outward propagation. In their paper of the range of possible indraft velocities, Trelles and Pagni [128] showed that indraft velocities of large fires can range from about 2 to 40 m/s. In the Project Flambeau burns, Countryman [3, 12, 13] also reports complicated airflow patterns and strong downdrafts which cannot be accounted for with simple plume theory.There seem to be two main theories in the literature as to what causes the high velocity inflows. One theory, advanced by Baum and McCaffrey [61] and Carrier et al. [7] is that large-scale vorticity in conjunction with heat release is responsible. These models contend that the entire fire plume slowly rotates. Note, however, that Church et al. [129] and McRae and Flannigan [130] characterize this type of motion as one type of fire whirl. In Baum and McCaffrey’s model (also used by [128, 131]), this rotation is caused by density gradients from the high heat release, and not necessarily by any imposed swirling caused by the ambient environment. The slow rotation of such a large mass of air above the ground translates to high velocity, purely horizontal and nonrotating flow at the ground. One unique feature of the Baum and McCaffrey model is that it treats the large area fire as an ensemble of randomly distributed individual fires of varying strengths. Because of the method chosen to represent the fire, the model is only valid for heights above the fuel bed where the plumes of the individual fires have not merged. The model of Carrier et al. [7] was intended to determine how long it would take to spin up the convective column and under what conditions this occurred. Based on the fact that the fire in Hamburg, Germany took two hours to develop, they concluded that the growth of swirl, at least in this case, was most likely due to the intensification of a preexisting vortex from earlier fires and bombings. Though this contradicts the Baum and McCaffrey model, the experiments and discussion by Church et al. [129] support this argument. The spatial orientation of individual fires may cause a swirling flow due to the interaction of the indrafts to each fire [132]. Carrier et al. [7] found that large diameter plumes spin up faster and proposed a set of four criteria that must be met for a “firestorm” to develop: heat release of 106 MW over a localized area for two to three hours, a preexisting weak vortex, low ambient winds, and a nearly dry-adiabatic lapse rate over the first few kilometers of the atmosphere.Because it seems unlikely that all the criteria for spin up of a convective column will be met, another theory, advanced by Smith et al. [60] and Small et al. [30] is proposed. These authors claim that buoyancy-induced pressure gradients are responsible for the large indrafts. Smith et al. [60] used a simple two-dimensional model of a convective column over a hot area to effectively show that near the fire a dynamic pressure gradient can cause high velocity inflow. This dynamic pressure gradient is caused by a balance between hydrostatic pressure and buoyancy. Buoyancy pushes the hot gases up while atmospheric pressure pushes fresh air at the ground in toward the fire horizontally to fill the gap left by the rising gases. Smith et al. [60] also suggest that the traditional “weakly buoyant” plume theories described above may be valid for a small range of plume heights sufficiently far away from the fire and any inversion layer above. Small et al. [30] use a similar model to Smith et al. [60] but include a volume heat addition and large density and temperature gradients. Small et al. [30] also numerically match their model results of the area near the fire to the results of traditional plume theory for the region far from the fire. In both the Smith et al. [60] and Small et al. [30] models, the fire is treated as a single large heat source (Figure 8). Small et al. [30] used their model to demonstrate how the maximum indraft velocity varies with fire radius, burning rate, and fire height (Figure 8). They showed that the maximum indraft velocity at first increases but eventually levels off (to approximately 40 m/s) with increases in both the fire radius and the burning rate. On the other hand, the maximum indraft velocity appears to be linear with fire height.Model results for flow-field streamlines for three fires in close proximity (from [98]).
(a)(b)A third, yet not well-explored explanation is proposed by Carrier et al. [134]. In this work, they use classic plume theory, but assume that the fire does not burn as a single fire, but a collection of individual fires. They hypothesize that the high indraft velocities are then due to the increased fire perimeter from this “multicellular burning zone.” This hypothesis was not further developed and in later works, these authors treated the fire as a subterranean point source. Interestingly, both the Baum and McCaffrey [61] and Small et al. [30] models reasonably replicate what little experimental data is available. However, the theories differ slightly in their predictions of the distance away from the fire that these indrafts extend [5]. The model of Baum and McCaffrey [61] predicts that the high velocity indrafts will extend much further from the fire compared to the model of Small et al. [30]. Without more detailed experimental data, it is impossible to say which model is more accurately portraying the physics.


## 4.5. Pulsation

Though not an effect of flame interactions, flame pulsation (or puffing) is an interesting phenomenon that can occur in stationary fires, such as a mass fire. This pulsation typically occurs in circular or axisymmetric  fires in weak ambient wind and is periodic in nature. Flame pulsation is important to many researchers because it can have a large influence on air entrainment rates and therefore heat release rates and pollution formation [135]. Observations of this phenomenon reveal the expansion of the flame near the base of the fire as a toroidal vortex, about the size of the fire diameter. As this vortex is shed and propagates upward, the flame necks inward giving the appearance of a “mushroom” shape. Figure 9 illustrates the process with time sequence of photos. Not all circular flames pulsate, however. Using dimensional analysis, Byram and Nelson [136] attempt to describe what type of fires will pulsate. They define a dimensionless “buoyancy” number, π2=Q̇c/(gD)0.5ρcpT, where Q̇c is the rate of convective heat release per area, D is the fire diameter, and ρ, cp, and T are the density, specific heat and temperature of the ambient air. Though no quantitative values are given, they argue that a fire will not pulsate if π2 is either too small (low heat release rate relative to large fire diameter) or too large (large heat release rate relative to small fire diameter).Thermal images of flame pulsation from [133].
(a)(b)Because this puffing occurs in nonreacting helium plumes, it is actually not caused by a combustion instability, but instead is produced by a fluid dynamic instability [137]. Though there is disagreement about the actual cause of the instability [138], the vortex is generally thought to be formed because of the interaction between gravity and the density gradient between the flame and ambient [135].Most of what has been learned about the characteristics of pulsation has been learned through experiments. Cetegen and Ahmed [137] showed that the toroidal vortex forms within one fire diameter above the flame base and that the frequency of the puffing is insensitive to the fuel or the heat release rate. By plotting the available data in the literature, Cetegen and Ahmed, and later Malalasekera et al.[133], showed that the pulsation frequency is proportional to the fire diameter raised to the negative one-half power (f~D-1/2) so that large fires pulsate at a much lower frequency than small fires. Though this correlation was developed using data from fires ranging from 0.1 m to 100 m in diameter (four orders of magnitude) using gaseous, liquid and solid fuels, Baum and McCaffrey [61] suggest that it may well hold for much larger fires as well. For a large fire with a diameter on the order of 20 km, Larson et al. [27] estimate that the pulsation will occur every 20 minutes. Though not accounted for in the above correlation, Malalasekera et al. [133] show that increasing fuel flow rates also result in a small increase in puffing frequency, especially for small fire sizes. Because of this, Malalasekera et al. [133] correlated the puffing frequency in a slightly different manner using the dimensionless Strouhal number (ratio of oscillation frequency to one over the characteristic time of convection) and Froude number (ratio of inertia force to gravitational force) which retains the same dependency on fire diameter but allows for a correction due to changes in fuel flow velocity.


## 4.6. Convection Column

Mass fires are also described as having very tall convection columns, or smoke plumes with large cloud structures due to the moisture release from combustion [32]. As discussed in the section on indraft velocities, the entrainment of cold, ambient air slows the rise of the hot gases by cooling them. Additionally, the density of the ambient air itself decreases with elevation. As the hot gases rise and cool, the density difference driving their upward motion disappears. It follows then that the top of the smoke plume corresponds to the where the combustion products stop moving. As the fire diameter grows, however, the entrainment predicted by classic plume theory becomes less effective. Entrainment occurs at the perimeter of the plume and with large fire sources there is such as large core of hot gases that entrainment is less effective at slowing the rise of the combustion products [126] (Palmer 1981). Thus, it takes longer to entrain enough cold air to slow the combustion products and therefore the smoke plume becomes taller. For example, a lack of entrainment to the convection column was noted and discussed by Taylor et al. [95] on a large prescribed burn. In fact, the plume from a sufficiently large mass fire may be almost as wide as it is tall, so Brode and Small [31] and Palmer [126] contend that air entrainment is not likely to be a major influence on plume height and that it is the structure of the atmosphere itself that is the limiting factor. The plume of large mass fires is therefore more sensitive to atmospheric gradients, inversion heights, and upper atmosphere crosswinds (see also [33]). Brode and Small [31] show that the tropopause/stratosphere transition may be what actually caps the smoke plume. Note, these theories contradict the suggestion of Smith et al. [60] that the traditional plume theory holds at some intermediate height above the ground. Perhaps the scale of the fires modeled by Smith et al. [60] was not large enough to see this effect.Palmer [126] describes the interesting structure of the convection columns that formed during the Project Flambeau tests. In the first few minutes of these large-scale burns, the majority of the gaseous combustion products were contained in a “bubble” near the fire. Once the “bubble” got sufficiently hot, the associated buoyancy was enough to overcome the surface drag forces and the bubble rose. As the bubble rose, a vortex ring would form in a similar manner described above with respect to flame pulsations. Regardless of the atmospheric stability, this vortex ring would rise until it encountered a region of vertical wind shear. The vertical wind shear weakens the vortex enough for the plume to then follow the prevailing horizontal winds. Palmer [126] also notes that the “exterior form of the convection column at a particular altitude was determined by the initial vortex bubble as it passed that altitude.” Most of the plumes in these fires began to rotate as a single vertical vortex, as suggested by the Baum and McCaffrey [61] model. This rotation further inhibits entrainment, which would also prevent the use of classic plume models for mass fires [139].


## 4.7. Summary of Interaction Effects

As the individual spot fires grow together, they will begin to interact. This interaction will increase the burning rates, heat release rates, and flame height until the distance between them reaches a critical level. At the critical separation distance, the flames will begin to merge together and burn with the maximum rate and flame height. As these spot fires continue to grow together, the burning and heat release rates will finally start to decrease but remain at a much elevated level compared to the independent spot fire. The flame height is not expected to change significantly. The more spot fires, the bigger the increase in burning rate and flame height.


## 5. Needs for Further Research and Application

The characteristics of many fire interactions have been examined and reported in the research literature, leaving little doubt that local spread and behavior experienced by wildland fire personnel can be greatly influenced by fire configurations at larger scales. The ignition patterns and “suppression fire” tactics used in firefighting [36, 140] depend on understanding these interactions. However, questions remain about how to extend the findings of fundamental research, to the field-scale for wildland fires and mass fires. In particular, there is no clear method to determine the minimum separation distance between two fires for interaction and merging to occur. The influence of ambient winds or topography on interactions is directly relevant to wildfire management activities and tactics but has not been explored. Large area fires were discussed as an extreme case of fire interactions and often behave quite differently than propagating line fires. Just how much area must be ignited to display “mass fire” characteristics is unknown. Even in the Project Flambeau experiments, Countryman [3] argues that even these large fires were not large enough to be considered mass fires. Both Byram [141] and Thomas et al. [142] developed scaling laws in an attempt to answer this question, but many potentially limiting assumptions were made in the development and the laws were not validated. Baldwin and North [15] attempt to quantify the minimum area for urban applications based on city layout and historical fires, but their estimations are admittedly crude. As discussed, there is no consensus in the literature about the convection column dynamics of mass fires and what mechanism is responsible for the reported strong indrafts. These suggestions are merely a starting point, as the subjects of fire interactions and mass fires clearly involve a great deal of physics and require the union of many fields of study.



---

*Source: 548328-2011-04-06.xml*