The Simulation of Different Combustion Stages of Micron-Sized Aluminum Particles

Featured Application: As a kind of high-energy green energy source with great potential, aluminum particles are still not well understood in their combustion mechanism. The research results of this paper propose a new type of combustion model, which can more accurately simulate the combustion characteristics of aluminum particles, and provide theoretical guidance for the large-scale application of aluminum particles. Abstract: In this study, a quasi-steady combustion model of an aluminum particle is established, which is more accurate to simulate the physical combustion process. Detailed gas-phase reaction mechanism and surface reaction mechanism are considered. Moreover, the particle temperature is not constant in this work, which is calculated in different combustion stages. The judgement standard of each combustion stage is from observational data in the experiment and the simulation results of combustion durations of each stage, and distribution of ambient temperature and gas-phase species proﬁles are in good agreement with experimental results. The calculation results show that in the ﬁrst stage, burning rate of the particle is the fastest, and in the second stage, particle temperature can drop to more than 100 K below the boiling point for the large particles, which is slightly below the boiling point for small ones. As the combustion stage changes, Da number is going to keep going down, which will lead to the transition of combustion method from diffusion-limited control to kinetic-limited control for an aluminum particle.


Introduction
Aluminum is widely used as a metal additive in solid propellants due to its good energy characteristics [1]. Due to its advantages of wide source, low cost, combustion products containing only alumina and easy for recycling, aluminum has the potential to be a new type of renewable green energy [2][3][4][5][6]. Nano-sized aluminum particles are a new kind of high-energy material, which is very suitable for microscale energy-demanding systems [7]. In addition, aluminum is easy to store and transport over long distances, which is the basis for the large-scale use of aluminum particles as a new type of green energy. However, aluminum is a kind of metal fuel; its melting point and boiling point are much higher than fossil fuels, and during the combustion process, aluminum oxide will deposit on the surface to form an oxide cap [8], which hinders the evaporation process of high-temperature aluminum droplets. The combustion process of aluminum particles is complicated and requires further in-depth study of its combustion mechanism.
The experimental and theoretical studies on the combustion characteristics of the aluminum particles have been extensively carried out since the 1950s. Beckstead [9] t = D 1.8 6.27 × 10 −7 (X O 2 + 0.6X H 2 O + 0.22X CO 2 )p 0.1 T 0.2 (1) Though this empirical correlation is widely used, Lynch, et al. [10] has reported it to be unsuitable for small particle sizes and has proposed another correlation for burning time for small particle sizes. Feng, et al. [11,12] accurately observed a single aluminum particle and obtained the burning time of a large number of micron aluminum particles in a complex high-temperature gas. Bucher, et al. [13] measured the distribution of ambient temperature and gas phase composition during the combustion process and suggested that the position of the peak of the AlO concentration represented the position of the flame. Dreizin [14] measured the changes in the aluminum particle temperature and divided the combustion process of the aluminum particles in air into three stages just as Figure 1 shows. Recently, Marsh, et al. [15] observed a similar phenomenon by the holography technique. Though this empirical correlation is widely used, Lynch, et al. be unsuitable for small particle sizes and has proposed another co time for small particle sizes. Feng, et al. [11,12] accurately observe particle and obtained the burning time of a large number of micro in a complex high-temperature gas. Bucher, et al. [13] measured the ent temperature and gas phase composition during the combustion p that the position of the peak of the AlO concentration represente flame. Dreizin [14] measured the changes in the aluminum particle vided the combustion process of the aluminum particles in air int Figure 1 shows. Recently, Marsh, et al. [15] Figure 1. The three stages during the combustion process of an aluminum p Brzustowski and Glassman [16] proposed that the combustion o was similar to the hydrocarbon fuels. Liang and Beckstead [17,18] pr sional unsteady state/quasi-steady state aluminum combustion mod rate chemical reactions and oxide cap. Further, Widener, et al. [19] by changing the ambient atmosphere to the rocket environment and lated the burning time based on the assumption that the mass consum linearly. Washburn, et al. [20,21] incorporated the detailed mechan chemical reactions and studied the transition from diffusion-limi combustion control. Glorian, et al. [22] further supplemented the d the surface reactions and regarded the evaporation process as a surf al. [23,24] suggested that the conversion rate of AlO and Al2O at th termined the combustion control method.
Most of the previous aluminum combustion models are quasirect to use a quasi-steady state result to represent the entire combu actual combustion process, the temperature of the particles is constan panied by the growth of oxidation caps, and sometimes the particles explode. According to the experimental results of Dreizin [14], this w particle temperature in Stage I is kept at the boiling point, which is r pressure [19] as:  Brzustowski and Glassman [16] proposed that the combustion of aluminum particles was similar to the hydrocarbon fuels. Liang and Beckstead [17,18] proposed a twodimensional unsteady state/quasi-steady state aluminum combustion model, including the finite rate chemical reactions and oxide cap. Further, Widener, et al. [19] improved the model by changing the ambient atmosphere to the rocket environment and subsequently calculated the burning time based on the assumption that the mass consumption rate decreased linearly. Washburn, et al. [20,21] incorporated the detailed mechanism of the gas phase chemical reactions and studied the transition from diffusion-limited to kinetic-limited combustion control. Glorian, et al. [22] further supplemented the detailed mechanism of the surface reactions and regarded the evaporation process as a surface reaction. Bojko, et al. [23,24] suggested that the conversion rate of AlO and Al 2 O at the particle surface determined the combustion control method.
Most of the previous aluminum combustion models are quasi-steady, but it is incorrect to use a quasi-steady state result to represent the entire combustion process. In the actual combustion process, the temperature of the particles is constantly changing, accompanied by the growth of oxidation caps, and sometimes the particles will rotate and micro-explode. According to the experimental results of Dreizin [14], this work assumes that the particle temperature in Stage I is kept at the boiling point, which is related to the ambient pressure [19] as: (2) p is united in atm in the above equation and the particle temperature in Stage II is calculated by the energy balance equation:

of 16
The heat absorbed by the evaporation of aluminum is regarded as part of the heat of the surface reaction, and the radiative heat absorption of the particle is not considered in the energy equation, because the calculation results in [8] suggested that the influence of radiant heat was small and could be ignored. The goal of this article is to establish a more accurate quasi-steady combustion model, considering the different stages of combustion, as well as detailed gas phase and surface chemical reaction mechanisms.

Model Formation
This article focuses on the combustion process of a burning aluminum particle and tries to use different quasi-steady states to model the different combustion stages. The particle is assumed to be spherically symmetric as Figure 2 shows and oxide cap is ignored. A lumped capacitance approach is used because the Biot number (h p D/λ p ) is estimated to be no larger than 0.01 for all cases considered. Most of the oxidizers in the environment will react with the vaporized aluminum vapor and form a high temperature diffusion flame zone, while part of them will cross the flame zone and react with the liquid aluminum at particle surface. Both gas phase reactions and surface reactions considered in this work are regarded as finite rate reactions and detailed reaction mechanisms are used.
The heat absorbed by the evaporation of aluminum is regarded the surface reaction, and the radiative heat absorption of the particl the energy equation, because the calculation results in [8] suggested radiant heat was small and could be ignored. The goal of this article accurate quasi-steady combustion model, considering the different as well as detailed gas phase and surface chemical reaction mechan

Model Formation
This article focuses on the combustion process of a burning al tries to use different quasi-steady states to model the different co particle is assumed to be spherically symmetric as Figure 2 shows nored. A lumped capacitance approach is used because the Biot nu mated to be no larger than 0.01 for all cases considered. Most of the ronment will react with the vaporized aluminum vapor and form a fusion flame zone, while part of them will cross the flame zone and aluminum at particle surface. Both gas phase reactions and surface in this work are regarded as finite rate reactions and detailed reac used.

Gas-Phase Analysis
The diffusion and reaction of gas phase substances follow the of conservation of mass, species, momentum and energy, as shown respectively: the right side of Equation (5) successively represents the mass chang and chemical reaction. Schematic diagram of gas phase reaction and surface reaction process between oxidizer with aluminum vapor and aluminum particle.

Gas-Phase Analysis
The diffusion and reaction of gas phase substances follow the governing equations of conservation of mass, species, momentum and energy, as shown in Equations (4)-(7), respectively: ∂ρ ∂t the right side of Equation (5) successively represents the mass change caused by diffusion and chemical reaction. ∂ρ the right side of Equation (6) successively represents the momentum change caused by pressure and viscous force.
the right side of Equation (7) successively represents the total energy change caused by the heat conduction, species diffusion and surface force work. Since the model is quasi-steady state, the change of each physical quantity in the equation with time is not considered. The detailed gas-phase mechanism employed to model the aluminum combustion is taken from [22,25] and is shown in Table 1. As the aluminum particles burned in the air react weakly with nitrogen, reaction and species involving nitrogen are discarded and argon is used as a dilution gas. The gas-phase mechanism includes K g = 10 gas-phase species and N g = 16 reactions. Any reaction in the gas-phase mechanism can be represented as: where A k is the chemical symbol for the k th species and the rate-of-progress variable q i for the i th gas-phase reaction is calculated by Equations (9)-(11): where the preexponential factor A i , temperature exponent b i and activation energy E i are given in Table 1. Thus, source term in Equation (5) can be obtained by following equation: . Condensation reaction is vital for the gas-phase reactions due to the condensation heat and it has been regarded as a reversible reaction with high kinetic parameters, as the liquefaction process of the high-temperature gaseous alumina is extremely rapid. However, the combustion flame temperature of aluminum is close to 4000 K; thus, the reverse reaction of the condensation process should not be ignored. For gas-phase species, the transport properties are needed to model the diffusion flame. In this work, the coefficients of molecular diffusion, thermal conductivity and viscosity were estimated with the usage of known formulas [26]. The dipole moment µ, polarizability α, Van-der-Waals collision diameters ε and well depth of Lennard-Jones potential, needed for these estimations were taken from [27,28] and are presented in Table 2. As the particle size of the generated alumina smoke is small (~1 µm), the transport properties of condensed phase Al 2 O 3 are identical to its gas phase [29]. Table 2. Parameters taken for the calculation of species transport properties [27,28].

Species Molecular Geometry
Lennard-Jones Potential Well Depth,

Surface Formalism
The mathematical formalism to describe surface kinetics for events such as adsorption, desorption, surface reactions and deposition is developed in this section. There are three types of species in surface formalism: gas-phase, surface and bulk. The first is a species in the gas phase above the particle surface. A surface species denoted by (S), is defined to be the chemical species on the top-most layer of the particle which means the liquid-gas interface. For example, the evaporation and absorption process (Al + Al(L) = Al(B) + Al(L)) of aluminum are shown in Figure 3. The symbol Al represents the gas atom of Al; Al(L) represents the open site at the surface, which will be occupied by gas species to form surface species; Al(B) represents the liquid aluminum particle. During the combustion of the aluminum, the two processes are going on at the same time, and when equilibrium is reached, the evaporation rate will be much greater than the absorption rate.  The detailed surface reaction mechanism employed to model t tion is taken from [22] and is shown in Table 3. The surface reaction in the general form: The detailed surface reaction mechanism employed to model the aluminum oxidation is taken from [22] and is shown in Table 3. The surface reactions can be represented in the general form: In any surface reaction, gas-phase species can participate in surface chemistry as well. The rate-of-progress variable for the i th surface reaction can also be written in the similar form as gas-phase reaction: The only difference is that [X k ] can represent either a volume (mol/m 3 ) or a surface molar concentration (mol/m 2 ). Surface molar concentration is calculated as: Γ and σ k is obtained from quantum chemistry [30,31]. For bulk species, [X k ] is taken as 1. For absorption reactions, their forward rate k f ,i is calculated by stick coefficient r in following equation:

Numerical Method
An axisymmetric two-dimensional grid was established to model the particle and the surroundings. The size of the computation domain is 50 times of the particle radius. The number of radial grids is N1 and the grid is refined close the center, and the number of circumferential grids is N2, which are evenly distributed. Half of a typical grid with N1 and N2 of 100 is shown in Figure 4. In this study, the grids with N1 × N2 being 50 × 50, 100 × 50, 100 × 100, and 200 × 100 were calculated to study the influence of the number of grids. The calculation model is a laminar component transport model with finite rate of chemical reaction, and a steady pressure-based solver is used for the numerical calculations. The coupled method is used for the pressure-velocity coupling, and the second-order upwind difference scheme is used for the continuity, component, momentum and energy equations. The simulation results of the ambient temperature in Y direction are shown in Figure 5. It can be seen that grid numbers have resulted in minor differences, thus, N1 × N2 of 100 × 50 was chosen to study the combustion of aluminum particles.
Appl. Sci. 2021, 11, x FOR PEER REVIEW shown in Figure 5. It can be seen that grid numbers have resulted in minor d thus, 1 × 2 of 100 × 50 was chosen to study the combustion of aluminum    [22,31].

Temperature and Species Profile
Bucher, et al. [13] studied the combustion of a free falling aluminu μm) in air. The authors measured the temperature and composition d the particle. The data was obtained at the stable radiation signal durin

Temperature and Species Profile
Bucher, et al. [13] studied the combustion of a free falling aluminum particle (D = 210 µm) in air. The authors measured the temperature and composition distribution around the particle. The data was obtained at the stable radiation signal during the early stage of combustion. The combustion at this moment could be considered to be in the boiling combustion stage. As the data was acquired after the particles have been burned for a while, thus, in order to achieve an effective comparison with the experimental results, the combustion of an aluminum particle (D = 200 µm) in an O2/Ar = 21/79 environment was measured. The particle temperature was set at 2713 K (boiling point at 1 atm), with the ambient temperature of 300 K and flow velocity of 0.3 m/s. The flow field and temperature distribution around the particle are shown in Figure 6a, whereas the experimental and computed profiles are presented in Figure 6b as a function of the non-dimensional distance to the surface Y/R (R is the particle radius).
Appl. Sci. 2021, 11, x FOR PEER REVIEW 9 of 19 bution around the particle are shown in Figure 6a, whereas the experimental and computed profiles are presented in Figure 6b as a function of the non-dimensional distance to the surface Y/R (R is the particle radius).
AlO radical is often interpreted as a marker for the flame position, and the simulated AlO distribution is noted to be in good agreement with the experimentally measured one. The simulated peak is at Y~2.2R, whereas the experiments suggest it to lie further from the surface (at Y~2.5R). The calculated maximum flame temperature (around 4000 K) is also in good agreement with the experimental result. In addition, it can be seen that the simulation results obtained in this study are more accurate than the previous literature reports. Figure 7 shows the fraction of sites occupied by the surface species, Al(L), Al2O(S) and Al2O2(S), are the most important surface intermediates, whereas the fractions of the other three species are close to 0. The flow field and temperature distribution around the particle (a) and comparison of the simulated AlO normalized mole fraction and temperature (K) profiles with the experimental measurements reported by Bucher [13] and simulated results reported by Glorian [22] (b).  AlO radical is often interpreted as a marker for the flame position, and the simulated AlO distribution is noted to be in good agreement with the experimentally measured one. The simulated peak is at Y~2.2R, whereas the experiments suggest it to lie further from the surface (at Y~2.5R). The calculated maximum flame temperature (around 4000 K) is also in good agreement with the experimental result. In addition, it can be seen that the simulation results obtained in this study are more accurate than the previous literature reports. Figure 7 shows the fraction of sites occupied by the surface species, Al(L), Al 2 O(S) and Al 2 O 2 (S), are the most important surface intermediates, whereas the fractions of the other three species are close to 0.

Burning Duration of Different Combustion Stages
The experimental measurements involving the burning duration of the different combustion stages are scarce. Only Dreizin [14] analyzed the burning of an aluminum particle with a diameter D = 165 µm in air. Figure 8 shows the results from [14]. In order to estimate the burning time from mass consumption . m Al , some ancillary assumption must be applied: 1.
. m Al decreased linearly with the residual aluminum mass; 2.
Combustion ends as the residual aluminum mass reaches 0.01 of the initial mass. (a) (b) Figure 6 The flow field and temperature distribution around the particle (a) and simulated AlO normalized mole fraction and temperature (K) profiles with the e measurements reported by Bucher [13] and simulated results reported by Glori  The experimental measurements involving the burning duration combustion stages are scarce. Only Dreizin [14] analyzed the burning particle with a diameter D = 165 μm in air. Figure 8 shows the results fr to estimate the burning time from mass consumption , some anci must be applied: 1.

Burning Duration of Different Combustion Stages
decreased linearly with the residual aluminum mass; 2. Combustion ends as the residual aluminum mass reaches 0.01 of th At the beginning of the combustion process, the particle diameter experiment was the diameter of the pure aluminum particle. However, i due to the growth of the oxidation cap, the initial diameter in Stage III w from 110 μm to 100 μm. The duration and initial diameter ranges for th are shown in Table 4. In this study, the burning time of each stage is ca Equation (4): the density of aluminum was obtained by using Marion's correlatio 3236 − 0.604 , and the density at 2600 K is taken as the average density The mass consumption rate with the aluminum mass in the whole cess for an aluminum particle (D = 165 μm) is shown in Figure 9, and th experimental durations for each stage are shown in Figure 10. It can be Figure 8. The temporal variation of particle diameter in different combustion stages in [14].
At the beginning of the combustion process, the particle diameter measured in the experiment was the diameter of the pure aluminum particle. However, in the third stage, due to the growth of the oxidation cap, the initial diameter in Stage III was corrected here from 110 µm to 100 µm. The duration and initial diameter ranges for the different stages are shown in Table 4. In this study, the burning time of each stage is calculated by using Equation (4): the density of aluminum was obtained by using Marion's correlation [32]: ρ Al (T) = 3236 − 0.604T, and the density at 2600 K is taken as the average density. The mass consumption rate with the aluminum mass in the whole combustion process for an aluminum particle (D = 165 µm) is shown in Figure 9, and the computed and experimental durations for each stage are shown in Figure 10. It can be observed that the calculated durations of Stage II and Stage III are highly consistent with those measured experimentally; however, the time of Stage I exhibits a specific error. This may be due to the fact that the linear reduction of . m Al considered in this study is different from the actual combustion process in Stage I. However, the consistency for the other two stages shows that the combustion model still has a high accuracy and can be used to analyze the characteristics of the different combustion stages.
Appl. Sci. 2021, 11, x FOR PEER REVIEW the fact that the linear reduction of considered in this study is diffe tual combustion process in Stage I. However, the consistency for the o shows that the combustion model still has a high accuracy and can be use characteristics of the different combustion stages.

Combustion Time
for each particle size in the two combustion states was calculate and subsequently the fitting relationship between and mass of the obtained (Figure 11). It can be seen that decreases as the diam and in the boiling combustion state is always greater than that in the bustion state. In addition, for particles with same diameter, the temperatu in the evaporative combustion state is lower than the boiling point and diameter decreases. As the diameter is decreased, the surface reaction b tense and the particles absorb more heat; thus, the temperature is enha been experimentally observed that the particle temperature exhibits a per in the later stages of combustion [7]. As there is no basis for dividing the c Figure 9. Mass consumption rate with the aluminum mass in the whole combustion process for an aluminum particle (D = 165 µm).
Appl. Sci. 2021, 11, x FOR PEER REVIEW the fact that the linear reduction of considered in this study is differ tual combustion process in Stage I. However, the consistency for the o shows that the combustion model still has a high accuracy and can be use characteristics of the different combustion stages.

Combustion Time
for each particle size in the two combustion states was calculate and subsequently the fitting relationship between and mass of the obtained (Figure 11). It can be seen that decreases as the diame and in the boiling combustion state is always greater than that in the bustion state. In addition, for particles with same diameter, the temperatu in the evaporative combustion state is lower than the boiling point and diameter decreases. As the diameter is decreased, the surface reaction be tense and the particles absorb more heat; thus, the temperature is enhan been experimentally observed that the particle temperature exhibits a per in the later stages of combustion [7]. As there is no basis for dividing the co

Combustion Time
. m Al for each particle size in the two combustion states was calculated first (Table 5), and subsequently the fitting relationship between . m Al and mass of the two states was obtained ( Figure 11). It can be seen that . m Al decreases as the diameter is reduced, and . m Al in the boiling combustion state is always greater than that in the vaporized combustion state. In addition, for particles with same diameter, the temperature of the particle in the evaporative combustion state is lower than the boiling point and increases as the diameter decreases. As the diameter is decreased, the surface reaction becomes more intense and the particles absorb more heat; thus, the temperature is enhanced. It has also been experimentally observed that the particle temperature exhibits a period of increment in the later stages of combustion [7]. As there is no basis for dividing the combustion stage of particles with different diameters, all calculations are based on the classification basis presented in Table 4. Figure 12 shows the calculated duration of each combustion stage for the particles with different diameters. The total combustion time is proportional to the D 1.94 , which is consistent with the trend observed for the D 2 law.  Table 4. Figure 12 shows the calculated duration of each c for the particles with different diameters. The total combustion time is pr D 1.94 , which is consistent with the trend observed for the D 2 law.

Large Particles
As mentioned earlier, for an aluminum particle with D = 165 µm, Stage I lasts until 140 µm, whereas Stage II lasts until 100 µm. The calculation results for the particles with 160 µm and 120 µm diameters are used to represent Stage I and Stage II, respectively. Figure 13a shows the temperature and gas composition distributions of the two combustion stages. The highest temperature of the flame in Stage I is noted to be about 100 K higher than that in Stage II, and the temperature in Stage II drops faster after the flame. The peak position of AlO in Stage II is observed to be closer to the particle surface than that in Stage I, and the mole fraction of Al 2 O near the surface is much larger than that of AlO, where the content of AlO is close to 0.
As mentioned earlier, for an aluminum particle with D = 165 μm, Stage I lasts until 140 μm, whereas Stage II lasts until 100 μm. The calculation results for the particles with 160 μm and 120 μm diameters are used to represent Stage I and Stage II, respectively. Figure 13a shows the temperature and gas composition distributions of the two combustion stages. The highest temperature of the flame in Stage I is noted to be about 100 K higher than that in Stage II, and the temperature in Stage II drops faster after the flame. The peak position of AlO in Stage II is observed to be closer to the particle surface than that in Stage I, and the mole fraction of Al2O near the surface is much larger than that of AlO, where the content of AlO is close to 0.
To illustrate the trade-off between the near-surface reaction and diffusion, Figure 13b shows the trend of number vs. normalized radius: = ⁄ , in which = / and = / , is the reaction rate of Al2O = AlO + Al. The number values near the wall in Stage I and Stage II are observed to be much higher than 1, indicating that the two stages are in a diffusion-controlled combustion mode. Further, from Stage I to Stage II, the number decreases significantly, indicating that the diffusioncontrolled combustion eventually transforms into the chemical kinetics-controlled combustion. To represent the importance of the surface reaction, the ratio of the Al2O3 deposition rate and Al consumption rate was used. The calculated results for the surface reaction are shown in Figure 14. The proportion of the surface reaction is observed to be very small for large particle sizes; however, the importance of surface reaction increases with the combustion stage. The fraction of the sites occupied by the three major surface species is almost similar. On changing the combustion stage, the proportion of Al (L) increases slightly, while the proportion of the other two species decreases to some extent. To illustrate the trade-off between the near-surface reaction and diffusion, Figure 13b shows the trend of Da number vs. normalized radius: Da = τ di f f /τ rea , in which τ di f f = D 2 /D AlO and τ rea = [X AlO ]/q, q is the reaction rate of Al 2 O = AlO + Al. The Da number values near the wall in Stage I and Stage II are observed to be much higher than 1, indicating that the two stages are in a diffusion-controlled combustion mode. Further, from Stage I to Stage II, the Da number decreases significantly, indicating that the diffusion-controlled combustion eventually transforms into the chemical kinetics-controlled combustion.
To represent the importance of the surface reaction, the ratio of the Al 2 O 3 deposition rate . m Al2O3 and Al consumption rate . m Al was used. The calculated results for the surface reaction are shown in Figure 14. The proportion of the surface reaction is observed to be very small for large particle sizes; however, the importance of surface reaction increases with the combustion stage. The fraction of the sites occupied by the three major surface species is almost similar. On changing the combustion stage, the proportion of Al (L) increases slightly, while the proportion of the other two species decreases to some extent.

Small Particles
A particle with D = 30 µm was chosen to study the combustion of the small particles. The calculation results for the particle with 30 µm and 20 µm diameters are used to represent Stage I and Stage II, respectively. Figure 15a shows the temperature and gas composition distributions of the two combustion stages. Compared with the particle with D = 165 µm, the temperature around the particle is lower. The content of AlO near the surface cannot be ignored, and the peak position of AlO hardly changes from Stage I to Stage II. From Figure 15b, it can be seen that the Da number is much smaller than that for the large particles, indicating that the rate of the chemical reaction is similar to the diffusion rate, and the combustion is undergoing transition stage from the diffusion-controlled combustion to the kinetics-controlled combustion. The calculated results for the surface reaction of the small particles are shown in Figure 16.

Small Particles
A particle with D = 30 μm was chosen to study the combustion of the small partic The calculation results for the particle with 30 μm and 20 μm diameters are used to r resent Stage I and Stage II, respectively. Figure 15a shows the temperature and gas co position distributions of the two combustion stages. Compared with the particle with 165 μm, the temperature around the particle is lower. The content of AlO near the surf cannot be ignored, and the peak position of AlO hardly changes from Stage I to Stag From Figure 15b, it can be seen that the number is much smaller than that for large particles, indicating that the rate of the chemical reaction is similar to the diffus rate, and the combustion is undergoing transition stage from the diffusion-contro combustion to the kinetics-controlled combustion. The calculated results for the surf reaction of the small particles are shown in Figure 16.
/ is more than a quar indicating the importance of the surface reaction in small diameter particles. The fract of sites of Al(L) is higher than the sum of other two, indicating intense evaporation.

Small Particles
A particle with D = 30 μm was chosen to study the combustion of the small particles. The calculation results for the particle with 30 μm and 20 μm diameters are used to represent Stage I and Stage II, respectively. Figure 15a shows the temperature and gas composition distributions of the two combustion stages. Compared with the particle with D = 165 μm, the temperature around the particle is lower. The content of AlO near the surface cannot be ignored, and the peak position of AlO hardly changes from Stage I to Stage II. From Figure 15b, it can be seen that the number is much smaller than that for the large particles, indicating that the rate of the chemical reaction is similar to the diffusion rate, and the combustion is undergoing transition stage from the diffusion-controlled combustion to the kinetics-controlled combustion. The calculated results for the surface reaction of the small particles are shown in Figure 16.
/ is more than a quarter, indicating the importance of the surface reaction in small diameter particles. The fraction of sites of Al(L) is higher than the sum of other two, indicating intense evaporation.

Conclusions
In this study, the different combustion stages of an aluminum particle were sim lated by using a quasi-steady combustion model. The model considered the detailed g

Conclusions
In this study, the different combustion stages of an aluminum particle were simulated by using a quasi-steady combustion model. The model considered the detailed gas-phase and surface reaction mechanisms, which is more detailed and accurate than previous quasi-steady combustion models. This is the first report employing the quasi-steady-state model to simulate the different combustion stages and the results are observed to be in good agreement with the experimental data. During the aluminum combustion, the particle temperature keeps going down, which slows down the combustion rate. Particle temperature is a key factor affecting combustion intensity and cannot been seen as constant like in the previous model.
As the combustion process goes on, the flame temperature keeps decreasing, and the position of the flame is close to the surface and the surface reaction plays an increasingly importance role. The particle temperature in Stage II is lower than that in Stage I. For large particles, the temperature in stage II would be more than 100 K lower than the boiling point, while the value of small particles is slightly lower than the boiling point, indicating that the surface reaction of small particles of aluminum is more intense. The change of Da number reveals that the combustion model gradually transforms from diffusion control to kinetic control.
During the combustion process of aluminum particles, molten aluminum oxide will continuously deposit to form an oxidation cap, which may cause the particles to rotate and break, thereby affecting the combustion process. However, the model in this work ignores the influence of the oxidation cap. In the future, research on the growth process of the oxidation cap will be the focus of the aluminum particle combustion model.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Data Availability Statement:
The study did not report any data.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations h p
Convection coefficient of particle (W·m −2 ·K −1 ) D Diameter (m) D k Mass diffusion coefficient of the k th species in the mixture (m 2 ·s −1 ) λ p Thermal conductivity of particle (W·m −1 ·K −1 ) λ Thermal conductivity of the gas mixture (W·m −1 ·K −1 ) X Molar fraction (X k ) Molar concentration of the species (gas or surface) (mol m −3 or mol·m −2 ) p Pressure (Pa) Entropy of reaction (J·mol −1 ·K −1 ) K p,i Equilibrium constant in pressure units for the ith reaction K g Total number of gas phase species K s Total number of surface species N g Total number of gas phase reactions N s Total number of surface reactions v ki Forward stoichiometric coefficients for the k th species in reaction i v ki Backward stoichiometric coefficients for the k th species in reaction i k f ,i Forward rate constant for reaction i k r,i Backward rate constant for reaction i A i Preexponential factor of reaction i (cm-mol-s-K) b i Temperature exponent of reaction i q i Rate of progress of the i th reaction (mol·m −3 ·s −1 or mol·m −2 ·s −1 ) . ω g k Mass production rate for the k th species due to gas phase reactions (kg·m −3 ·s −1 ) . ω s k Mass production rate for the k th species due to surface reactions (kg·m −2 ·s −1 ) γ i Sticking coefficient for the i th reaction Γ Total density of open sites (mol·m −2 ) Z k Surface site fraction for the k th species Dimensionless number representing the speed of diffusion and reaction Y/R Ratio of Y coordinate value to particle radius