Evaluation of Aerosol Fire Extinguishing Agent Using a Simple Diffusion Model

Abstract

Aerosol fire extinguishing agents have been recognized as an effective replacement of Halon. The uniform diffusion time and the effective concentration of the aerosol fire extinguishing agent are important parameters for putting out a fire. In this work, the effective concentration is derived based on the diffusion equation, and its variation with changing diffusion coefficient and diffusion time is analyzed. The uniform diffusion time could then be conveniently estimated using this equation. Based on experimental data, the concentration kinetics of the aerosol is drawn and the relation between the diffusion coefficient and the uniform diffusion time is analyzed. It was found that the uniform diffusion time is not dependent on the shape of the closed room but dependent on the total room volume and the position of its diffusion source. This model is demonstrated as a facile tool for the convenient evaluation and reasonable application of the aerosol fire extinguishing agent by predicting the uniform diffusion times of extinguishing aerosol in closed rooms.

1. Introduction

Halon has been banned by most countries worldwide because of its adverse effect on the ozone layer. This has brought the development of alternative fire-fighting materials into attention, among which the aerosol fire extinguisher has gained popularity due to its advantages such as high efficiency, low cost, and easy maintenance [1, 2]. The pyrotechnic fire extinguishing aerosol consists of 0.001–1 𝜇m microparticle agglomerates generated from pyrotechnic agent burning [3, 4]. The microparticles have very large surface area, which increases the efficiency of flame quenching. Besides, the microparticles in violent turbulent motion can rapidly fill up available space in a “total flooding” manner, which enables special protection against particular risks [5]. The aerosol fire extinguishing agents have been widely used in various settings such as aircraft, ship, land combat vehicle, computer lab, and control room.

The turbulent diffusion of extinguishing aerosol is complicated because the diffusion process is a continuous flux during which slow physical and chemical changes take place, and the distribution of temperature and humidity is not uniform [6]. The uniform diffusion time is the time span from the initial release of the aerosol until the space is completely filled by the extinguishing aerosol, and the effective concentration is the critical concentration of the extinguishing aerosol that can put out a fire. Therefore, the uniform diffusion time and the effective concentration of the extinguishing aerosol are important parameters to ensure timely fire extinguishing. The diffusion course is influenced by not only the intrinsic aerosol properties but also environmental factors such as temperature, humidity, and air current. The diffusion involves mass transfer among multiple gaseous, liquid, and solid components with complicated chemical equilibrium, which makes it difficult to precisely calculate when a fire will die out, especially in a room of complicated structure. Therefore, the aim of this work is to establish a technical guideline for the evaluation and application of aerosol fire extinguishing agent. Based on performance tests, the diffusion model of the extinguishing aerosol is first built to explore the uniform diffusion time of the fire extinguishing aerosol and study the time of fire. Note that this paper relies on a statistic analysis and circumvent complicated issues such as the fluid mechanics and the reaction dynamics between flame and extinguishing agent. Instead, the objective is to find a simple method to estimate the time of fire to die out, which can then suggest the optimal allocation of the fire extinguishing aerosol, that is, the best structure, overall arrangement, and manner.

2. Diffusion Model of Extinguishing Aerosol

After the release of the aerosol, the extinguishing aerosol rapidly (within 0.3 s) generates an approximately spherical cloud of radius 𝑟₀ before its spontaneous diffusion, and then for a considerably long period the diffusion coefficient attenuates only very slowly. There are many models for aerosol diffusion [6, 7]. In this paper, the movement law of the extinguishing aerosol cloud can be described by the second order partial differential equation as 𝜕𝐶𝜕𝜕𝑡=𝐷²𝐶𝜕𝑥²+𝜕²𝐶𝜕𝑦²+𝜕²𝐶𝜕𝑧², (2.1) where (𝑥,𝑦,𝑧) is an arbitrary point in the infinite space at time 𝑡, 𝐶(𝑥,𝑦,𝑧,𝑡) is the concentration of the extinguishing aerosol at the point (𝑥,𝑦,𝑧), and 𝐷 is the diffusion coefficient. Idealization of 𝐶(𝑥,𝑦,𝑧,𝑡), 𝑅²=𝑥²+𝑦²+𝑧² as equivalent to a spherical surface then gives [8, 9] 𝐶()=𝑄𝑅,𝑡(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2/4𝐷𝑡}⇒𝑅²𝑄=4𝐷𝑡ln𝐶⋅(4𝜋𝐷)^3/2𝑡^3/2, (2.2) where 𝑄 is the amount of released extinguishing aerosol.

According to (2.2), a series of curves of aerosol concentration, time and aerosol cloud radius can be plotted. The aerosol concentration is normally distributed with regard to radius and rapidly diminishes as time elapses.

The concentration profile of 5000 g extinguishing aerosol released in an infinite space is shown in Figure 1. As clearly shown in Figure 1, at 1 m from the diffusion source, the aerosol concentration fell notably from 238 g·m⁻³ at 0.4 s to 35 g·m⁻³ at 2.0 s. In contrast, at 2 m from the diffusion source, the aerosol concentration initially increased to about 50 g·m⁻³ but then started to decline.

Figure 1 

Distribution of extinguishing aerosol concentration at different positions and diffusion times.

3. Practical Application of the Extinguishing Aerosol Model

3.1. Measurement of Diffusion Coefficient

Equation (2.2) can be transformed into 𝑄𝐶⋅(4𝜋𝐷)^3/2𝑡^3/2=𝑒^{𝑅2/4𝐷𝑡}. (3.1)

Let two fire sources be placed at different positions along the diffusion direction, and two groups of data are then measured. However, the diffusion coefficient does not remain constant during the turbulence diffusion of aerosol; rather, it is influenced by temperature, humidity, and air current and weakens gradually as time elapses. To reduce error in the evaluation, the best solution is to take the average of the diffusion coefficient.

According to (3.1) 𝑄𝐶₁⋅(4𝜋𝐷)^3/2𝑡₁^3/2=𝑒^{𝑅21/4𝐷𝑡1},𝑄𝐶₂⋅(4𝜋𝐷)^3/2𝑡₂^3/2=𝑒^{𝑅22/4𝐷𝑡2}, (3.2)(3.3) where 𝑅₁ and 𝑅₂ are distances to the diffusion source and 𝑡₁ and 𝑡₂ are the time when the two same fire sources die out, so 𝐶₁=𝐶₂.

Divide (3.3) by (3.2), then 𝑅𝐷=²¹𝑡₁−𝑅²₂𝑡₂𝑡4ln^23/2𝑡₁^3/2−1. (3.4)

The diffusion coefficient 𝐷 can thus be calculated from the measured values of 𝑅₁, 𝑅₂, 𝑡₁, 𝑡₂.

Due to the turbulence diffusion of the aerosol fire extinguishing agent, there are large errors in measuring 𝑅₁ and 𝑅₂. To avoid this problem, (3.4) can be transformed into 1𝐷=𝜋𝜋𝑅²₁𝑡₁−𝜋𝑅²₂𝑡₂𝑡4ln^23/2𝑡₁^3/2−1⇒𝑆₁𝑡₁−𝑆₂𝑡₂𝑡4𝜋ln^23/2𝑡₁^3/2−1, (3.5) where 𝑆₁ is the area of the circle with radius 𝑅₁, and 𝑆₂ is the area of the circle with radius 𝑅₂. Because the change in total area always has a well-defined pattern regardless of whether the aerosol cloud has an ideal spherical shape, the measurement errors of 𝑅₁ and 𝑅₂ are thus not introduced into the equation. The area can be determined by standardized image processing methods, such as edge detection. Figure 2 illustrates the formation of aerosol cloud.

Figure 2 

The expansion and formation of aerosol cloud.

The parameters 𝑆₁, 𝑆₂, 𝑆₃, 𝑆₄, 𝑆₅, 𝑆₆ and 𝑡₁, 𝑡₂, 𝑡₃, 𝑡₄, 𝑡₅, 𝑡₆ can be determined from the individual figures in Figure 2, thus allowing to calculate a series of diffusion coefficient as well as the average diffusion coefficient, shown in Table 1 by (3.5).

Table 1 

Calculation of diffusion coefficient.

The magnesium powder and strontium nitrate-based pyrotechnic aerosol fire extinguishing agents have diffusion coefficient in the range of 0.1∼1.0m²·s⁻¹, and the value of diffusion coefficient depends on the type of formulation. To facilitate calculation, 𝐷=1.0m²·s⁻¹ is assumed in the subsequent sections.

The calculated results are in good agreement with the measured ones [8, 9]. Note that the fast formation of the extinguishing aerosol causes vigorous turbulent motion of the aerosol cloud, which results in the much greater diffusion coefficient compared with the typical molecular diffusion of gases (𝐷≈1.0×10⁻⁵m²·s⁻¹).

3.2. Effective Concentration of the Extinguishing Aerosol

The effective concentration of the extinguishing aerosol that can put out a fire is 𝐶_𝑒. Assuming the effective concentration 𝐶_𝑒=100g·m⁻³ and having 5000 g extinguishing aerosol released in an infinite space, only an area of less than 1.5 m in radius can reach the effective concentration based on (2.2).

𝑅_𝑒 is defined as the spherical radius where 𝐶=𝐶_𝑒. By (2.2), the radius 𝑅_𝑒 of the spherical surface will change during diffusion. The curve describing the kinetic variation of radius is drawn in Figure 3. It can be seen that the radius 𝑅_𝑒 grows from 0.1 s to 0.4 s to reach a maximum of 1.5 m but then decreases to zero from 0.4 s to 1.09 s.

Figure 3 

The variation curve of 𝑅𝑒 with diffusion time based on (2.2).

Essential differences exist in the extinguishing behavior between the area inside 𝑅_𝑒 and the area outside 𝑅_𝑒. The aerosol concentration in the area inside 𝑅_𝑒 is always higher than the effective concentration, as shown in Figure 4, and the maximum radius 𝑅_𝑒 is less than 1.5 m. In this zone, the extinguishing aerosol can directly swallow the flame (Figure 4).

Figure 4 

Photographs of flame quenching inside the immediate effective concentration zone. Because of the effects of high concentration, the flame is suppressed directly. The photographs are from high speed camera, the time interval is known, and the height of flame can be determined, so the local diffusion coefficient can be calculated based on Figure 4.

In contrast, the aerosol concentration in the area outside 𝑅_𝑒 cannot reach the effective concentration from diffusion alone but gradually grows to the effective concentration in the closed room. The flame in this area becomes increasingly thin and faint and finally dies out, as shown in Figure 5. The uniform diffusion time is the time span between the release of the aerosol of the extinguishing agent and the final flame quenching.

Figure 5 

Photographs of flame quenching within the buffer area. It is different from the immediate effective concentration zone in Figure 4 as the flame cannot be suppressed directly, and the flame is smothered by the gradual increase of the aerosol concentration in the buffer area.

3.3. The Buffer Zone and Its Interface of Extinguishing Aerosol Diffusing in a Closed Room

A relatively closed room with 50 m³ spherical space has a radius of about 2.3 m. Because the effective concentration of extinguishing aerosol 𝐶_𝑒=100g⋅m⁻³, the space will require 5000 g extinguishing aerosol. Figure 1 is based on (2.2), and the application of (2.2) in a relatively closed room is now further explored.

Aerosol diffusion in a closed room is subject to the wall effect. Because the aerosol particles can hit and bounce from the wall, the aerosol concentration near the wall can rapidly increase. Figure 1 and (2.2) give the concentration gradients of the aerosol diffusion. Note that the aerosol concentrations near the wall cannot exceed the lowest concentration gradient from the diffusion source in the closed room. As a result, a buffer area forms near the wall, and the aerosol concentration increases uniformly in this area. As the aerosol diffuses, the buffer area will gradually widen and eventually reach the center of the diffusion source, as shown in Figure 6. At this point, the room becomes uniformly filled with the extinguishing aerosol.

Figure 6 

Schematic of extinguishing aerosol diffusion in a spherical closed room. 𝑅𝑡 is the radius of the interface between the buffer area and the concentration gradients, and 𝑅𝑐 is the radius to the wall of the closed room.

Equation (2.2) shows how the aerosol concentration varies from the diffusion source to the surrounding. The mechanism of the diffusion process is the mass transfer from high concentration area to low concentration area. Therefore, based on Figure 6, in which 𝑅_𝑡 is the radius of interface between the buffer area and concentration gradients, a range that is less than 𝑅_𝑡 in radius satisfies (2.2). During the diffusion, 𝑅_𝑡 becomes increasingly smaller and eventually reaches zero. Meanwhile, the buffer area expands gradually until a uniform aerosol concentration is reached throughout the room, at which time this uniform concentration is higher than or equal to the effective concentration.

The distribution of extinguishing aerosol can be calculated at any time and any point by the following, which is derived from (2.2)𝑄_𝜏=^𝑅𝜏₀4𝜋𝑅²⋅𝐶(𝑅,𝜏)𝑑𝑅=^𝑅𝜏₀4𝜋𝑅²⋅𝑄(4𝜋𝐷𝜏)^3/2𝑒^{−𝑅2/4𝐷𝜏}𝑑𝑅, (3.6) where 𝑅_𝜏 is the radius of a named point at time 𝜏, and 𝑄_𝜏 is the total quantity of aerosol within a radius 𝑅_𝜏. Therefore, the kinetic variation of aerosol quantity can be calculated and the variation of 𝑅_𝑡 can be analyzed.

Thus (3.7) can be derived based on (2.2) and (3.6) and Figure 3 to describe the variation of the 𝑅_𝑡 interface, which is the interface where forward diffusion concentration equals to the counter flow concentration: 𝑄(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2𝑡/4𝐷𝑡}=𝑄−∫^𝑅𝑡₀4𝜋𝑅²⋅𝑄/(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2/4𝐷𝑡}𝑑𝑅(4/3)𝜋𝑅³_𝑐−(4/3)𝜋𝑅³_𝑡. (3.7) Simplifying (3.7), the following is derived 43𝜋𝑅³_𝑐−43𝜋𝑅³_𝑡1(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2𝑡/4𝐷𝑡}+^𝑅𝑡₀4𝜋𝑅²⋅1(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2/4𝐷𝑡}𝑑𝑅=1, (3.8) where 𝑅_𝑐 is the radius of the spherical room shown in Figure 6. Equation (3.8) is the relationship function between the interface 𝑅_𝑡 and time. Figure 7 can be drawn from Figure 6 and (2.2) and (3.7).

Figure 7 

The variation curve of the interface between the diffusion area and the buffer area.

According to Figure 7, at 2 m from the diffusion source and 𝜏=0.3s, the concentration of counter flow reaches a balance with the concentration of forward diffusion because of the wall effect. The same balance is observed at 1.8 m from the diffusion source and 𝜏=0.4s. According to the total aerosol quantity and the room size, the average concentration is 100g⋅m⁻³, and the extinguishing aerosol will reach uniform diffusion when the concentration of counter flow reaches the average concentration. When 𝑅_𝑡=0, the aerosol reaches uniform diffusion in the spherical room, and (3.8) becomes 1𝑡=(36𝜋)^1/3𝑅²_𝑐𝐷=³𝑉².4𝜋𝐷 (3.9) Taking 𝑅_𝑐=2.3m and 𝐷=1m²·s⁻¹, the uniform diffusion time is 𝑡=1.09s, which is consistent with the end point in Figures 3 and 7.

3.4. Comparison of Diffusion Time in Different Room

Note that the above calculations consider a spherical room. If the room is rectangular with edges and corners, the wall effect would become very complicated. Nevertheless, for the rectangular room, obviously the extinguishing aerosol diffuses to vertically opposite corners (angles) regardless of the position of the diffusion source. As the mass transfer from the diffusion source proceeds, the concentration of extinguishing aerosol of the farthest corners continues to increase and the interface of the concentration balance will return to the diffusion source.

For example, assume a closed rectangular room with the dimensions of 𝐿×𝑊×𝐻 and with the source of extinguishing aerosol placed at the center. To simplify, the room is separated into eight symmetric parts as shown in Figure 8 and only one octant needs to be considered and analyzed. Based on Figure 8, no matter how the diffusion proceeds, it has an added wall effect from the diffusion source 𝑂 to point 𝐴. In fact, the process is very similar to the diffusion from the top angle of a tapered container to the top angle of another tapered container, both sitting on the bottom of each other. Assuming walls are rigid with smooth surfaces, two identical tapered containers are used to simplify the calculation process (Figure 9).

Figure 8 

Schematic of diffusion in the rectangular room.

Figure 9 

Schematic of diffusion in two tapered containers.

According to (3.8) and Figure 8, the height of the tapered container equals 𝑂𝐴/2 and the volume of the tapered container is one-sixteenth of the room. Using this setting, the total volume of the two tapered containers is equal to one-eighth of the rectangular space. Besides, assuming the diffusion source is at the center of the rectangular space, the diffusion modes of the two cases are similar, as the diffusion proceeds from the apex of one tapered container to that of the other. The aerosol diffusion in the two tapered containers is then calculated, as shown in Figure 9. Because the room size is 𝐿×𝑊×𝐻, ℎ=𝐿²+𝑊²+𝐻²4,𝑟²=3×𝐿⋅𝑊⋅𝐻4𝜋𝐿²+𝑊²+𝐻², (3.10) where ℎ is the height of the tapered container, and 𝑟 is the radius of the bottom. Let 𝑙 denote the side length of the tapered container, so 𝑙=𝑟²+ℎ². When 2ℎ>𝑅>𝑙, the concentration balance can be expressed as 𝑄(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2/4𝐷𝑡}=𝑄−∫^𝑅₀4𝜋𝑅²⋅𝑄/(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2/4𝐷𝑡}𝑑𝑅8×(1/3)𝜋(((2ℎ−𝑅+𝑥)/ℎ)𝑟)²(2ℎ−𝑅+𝑥)−(𝜋𝑥/6)3(((2ℎ−𝑅+𝑥)/ℎ)𝑟)²+𝑥², (3.11) where 𝑥 is a chord height changing with 𝑅 following 𝑅²=2ℎ−𝑅+𝑥ℎ𝑟²+(𝑅−𝑥)²⇒2ℎ−𝑅+𝑥ℎ𝑟²−2𝑅𝑥+𝑥²=0. (3.12) In contrast, when 𝑅<𝑙, the concentration balance gives 𝑄(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2/4𝐷𝑡}=𝑄−∫^𝑅₀4𝜋𝑅²⋅𝑄/(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2/4𝐷𝑡}𝑑𝑅8×(2/3)𝜋𝑟²ℎ−𝜋𝑥63(𝑅−𝑥)²𝑟²/ℎ+(𝑅−𝑥)²−(𝜋/3)(𝑅−𝑥)³/ℎ²𝑟², (3.13) where 𝑥 is a chord height changing with 𝑅 following (𝑅−𝑥)²+𝑅−𝑥ℎ𝑟²=𝑅²⇒𝑥=𝑅−𝑅ℎ𝑟²+ℎ². (3.14) Because when 𝑅<𝑙 a value may be assumed, namely, 𝑅=0, so 𝑥=0, and then (3.13) becomes 𝑡³=4𝑟⁴ℎ²9𝜋𝐷³. (3.15)

In practice, assuming the volume of the room is still 50m³ and the dimensions are 𝐿=3m, 𝑊=4m, and 𝐻=4.167m, then ℎ=1.63m and 𝑟=1.164m and it can be calculated that 𝑡=1.09s.

This result is not coincidental. In fact, regardless of the container shape, the Fick diffusion theory states that the diffusion proceeds to the farthest position. For either the spherical room or the two tapered containers, ignoring the intermediate processes of complicated calculation, the final results are always the same with (3.9), and the end points of Figures 3 and 7. Again notice the above rectangular room in Figure 8 with dimensions of 𝐿×𝑊×𝐻. Because the result relates only to the final equation based on the models of the spherical and the two tapered containers, the equation only needs to consider the final stage near the diffusion source. Diffusing from point 𝑂 to point 𝐴 is exactly equivalent to diffusing through one-eighth of the rectangular room. So 𝑄(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2/4𝐷𝑡}=𝑄−∫^𝑅₀4𝜋𝑅²⋅𝑄/(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2/4𝐷𝑡}𝑑𝑅8×𝑉/8−(1/8)×(4/3)𝜋𝑅³,4𝑉−3𝜋𝑅³1(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2/4𝐷𝑡}+^𝑅₀4𝜋𝑅²⋅1(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2/4𝐷𝑡}𝑑𝑅=1. (3.16) When 𝑅=0, 𝑉=(4𝜋𝐷𝑡)^3/2⇒𝑡=³𝑉².4𝜋𝐷 (3.17)

The result is identical with (3.9). Despite the different room shapes, when 𝑅=0, the uniform diffusion time is the same. The results of the above models are analyzed again.

The result of the spherical room is (3.9) 1𝑡=(36𝜋)^1/3𝑅²_𝑐𝐷⇒𝑡³=1𝑅36𝜋𝐷^6𝑐3⇒𝑡³=164𝜋²(4/3)𝜋𝑅³_𝑐²𝐷³⇒𝑡³=164𝜋²𝑉²𝐷³⇒𝑡=³𝑉².4𝜋𝐷 (3.18)

The result of the two tapered containers room is (3.9) 𝑡³=4𝑟⁴ℎ²9𝜋𝐷³=(2/3)𝜋𝑟²ℎ²(𝜋𝐷)³. (3.19)

Because the room size 𝑉=𝐿×𝑊×𝐻, and then (1/3)𝜋𝑟²ℎ=𝑉/16, so 𝑡³=4(𝑉/16)²(𝜋𝐷)³⇒𝑡=³4𝑉²/16²𝜋𝐷⇒𝑡=³𝑉²/64=𝜋𝐷³𝑉².4𝜋𝐷 (3.20)

Therefore the obtained diffusion time in the rectangular room is the same in (3.17), and all results are the same as in (3.9).

In all the three models above, the diffusion source is placed at the center of the closed room. If the diffusion source is at the corner of the closed room, the result will be different. For example, if the diffusion source is placed at the corner, such as at point 𝐴 in Figure 8, letting the aerosol diffuse through one octant, this will be equivalent to having eight times the initial aerosol mass when the diffusion source is at the room center. The final equation then becomes 8𝑄(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2/4𝐷𝑡}=𝑄−∫^𝑅₀4𝜋𝑅²⋅8𝑄/(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2/4𝐷𝑡}𝑑𝑅𝑉−(1/8)×(4/3)𝜋𝑅³ (3.21) and can be simplified as 1𝑉−6𝜋𝑅³8(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2/4𝐷𝑡}+^𝑅₀4𝜋𝑅²⋅8(4𝜋𝐷𝑡)^3/2𝑒^{−𝑅2/4𝐷𝑡}𝑑𝑅=1. (3.22) When 𝑅=0, 𝑡=³𝑉².𝜋𝐷 (3.23)

Compared with (3.9), the diffusion time is four times than when the diffusion source is at room center. Clearly, if the diffusion source is placed at the middle of the line 𝐴𝐵 in Figure 7, the result becomes 𝑡=³²𝑉₃√.4𝜋𝐷 (3.24) If the diffusion source is placed at the center of the 𝐴𝐵𝐶𝐷 plane in Figure 8, the result becomes 𝑡=³²𝑉₃√.16𝜋𝐷 (3.25)

The above results indicate that the differences in uniform diffusion time are due to the position of the diffusion source, and the shape of the closed room does not affect the results directly.

4. Conclusion

It is very difficult to estimate the time from the release of the aerosol to extinguishment due to the turbulent diffusion of aerosol fire extinguishing agent, the complicated room structure, and the environmental conditions. A simple method is developed to determine the uniform diffusion time of pyrotechnic aerosol fire extinguishing agent. First, the diffusion coefficient is experimentally determined, then a series of models are established using the Fick diffusion theory, and finally the uniform diffusion time is analyzed in detail for different shapes of closed room and different diffusion source positions. The conclusion is shown that(i)This method can be used to easily estimate the extinguishing time. (ii)It is found that the uniform diffusion time depends on the total room volume, the diffusion coefficient, and the diffusion source position, but not on the shape of the room. (iii)The equations derived from the above models have illustrated the potential of the proposed method in predicting the diffusion time of extinguishing aerosols. (iv)Because the diffusion time involves the effective concentration, the diffusion coefficient, and the diffusion source position, this method can be applied to determine the optimal allocation and quantity of aerosol fire extinguishing agent and evaluate the effect of fire extinguishment.

Acknowledgment

The authors thank the National Natural Science Foundation of the People’s Republic of China for financial support (Grant no. 51076066).