Next Article in Journal
Tailoring Microstructure and Mechanical Properties of Additively-Manufactured Ti6Al4V Using Post Processing
Next Article in Special Issue
Ceramic Biomaterial Pores Stereology Analysis by the Use of Microtomography
Previous Article in Journal
Modeling Bainitic Transformations during Press Hardening
Previous Article in Special Issue
Modelling of the Guillotine Cutting Process by Means of a Symmetrical Blade with the Defined Geometry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Analysis of the Activated Combustion High-Velocity Air-Fuel Spraying Process: A Three-Dimensional Simulation with Improved Gas Mixing and Combustion Mode

School of Materials Science and Engineering, North University of China, Taiyuan 030051, China
*
Author to whom correspondence should be addressed.
Materials 2021, 14(3), 657; https://doi.org/10.3390/ma14030657
Submission received: 2 January 2021 / Revised: 22 January 2021 / Accepted: 27 January 2021 / Published: 31 January 2021

Abstract

:
Owing to its low flame temperature and high airflow velocity, the activated combustion high-velocity air-fuel (AC-HVAF) spraying process has garnered considerable attention in recent years. Analyzing the velocity field, temperature field, and composition of AC-HVAF spray coatings plays a vital role in improving the quality of coatings. In this study, an actual spray gun is adopted as a prototype, and the radial air inlets are introduced to improve the reaction efficiency so that the chemical reaction can be completed in the combustion chamber. Furthermore, a complete three-dimensional (3D) model is established to examine the effects of radial inlets and porous ceramic sheet on the combustion and flow fields. The hexahedral cells are used to discretize the entire model for reducing the influence of false-diffusion on the calculation results. The gas flow field is simulated by the commercial Fluent software, and the results indicate that the porous ceramic sheet effectively reduces the turbulent dissipation of the airflow with a good rectification effect (the ceramic sheet ensures a consistent airflow direction). The radial inlets and the porous ceramic sheet promote the formation of vortex in the combustion chamber, increase the residence time and stroke of the gas in the combustion chamber, and improve the probability of chemical reactions. In addition, it is observed that the stability of velocity for the airflow is strongly related to the airflow density.

1. Introduction

The activated combustion high-velocity air-fuel (AC-HVAF) spraying is a novel surface spraying process developed on the basis of high velocity oxygen-fuel (HVOF) spray coating process [1]. Compared with HVOF spraying, AC-HVAF spraying has the advantages of low flame temperature, low cost, and efficient processing [2,3,4,5]. Besides that, AC-HVAF coatings have a lower oxygen content, lower porosity, denser layer, and beneficial compressive stress on the surface as compared with HVOF coatings [6,7,8]. Therefore, AC-HVAF spraying has received considerable attention and can be used to prepare a wide range of coatings, including pure metals, alloys, intermetallic compounds, nano-coatings, and amorphous coatings [2,9,10,11,12,13].
Although several studies have focused on understanding the formation mechanism of AC-HVAF coatings, the high airflow velocity and the difficulty in measuring physical quantities pose a significant challenge for further research. Computer simulation technology has become an effective means to solve this problem. Over the past decades, there have been relatively few simulation studies on HVAF spraying. Gorlach simulated the thermal stress and temperature changes of different materials of sandblasting nozzles on spray gun over time using a two-dimensional (2D) swirl axisymmetric model and provided useful suggestions for the selection of HVAF nozzle materials [14]. To realize cold spraying through HVAF coating, Yuan et al. injected liquid water into the high-velocity airflow through a small radial hole downstream of the throat of the spray gun, which reduced the temperature of the airflow. Further, they established a three-dimensional (3D) model to obtain the optimum mass flow rate of liquid water and the flame characteristics at such a flow rate [15]. Jiang et al. assumed that fuel and air were fully mixed and regarded the small holes on the ceramic sheet in the spray gun as the gas inlets to establish a 3D model. Based on this model, the airflow velocity, temperature, pressure, and other parameters under a certain mass flow rate were obtained, which were used to simulate the flight process of particles [16]. Based on a 2D swirl axisymmetric model, Gao et al. examined the influence of different mass flow rates and nozzle geometries on the flow field of HVAF spraying [17]. In addition, although the structure of spray gun and the spraying principle of other spraying such as HVOF spraying [18,19,20,21,22,23,24,25,26,27,28,29,30,31,32], plasma spraying [33,34,35], and cold spraying [36] are extremely different from HVAF spraying, the related simulation studies also provided a useful reference for this study on AC-HVAF spraying.
The above studies have provided interesting results on AC-HVAF spraying. However, because of their diverse research objectives, different studies used different models, and the actual spray gun was simplified to various degrees, which may lead to errors in the simulated results with respect to the actual situation. In particular, the chemical reaction in the combustion chamber plays a decisive role in the entire flow field. However, the earlier studies ignored the following question: Does the chemical reaction continue within the entire spray gun? Apparently, the degree of chemical reaction determines the heat released in the reaction, which profoundly affects the composition, pressure, temperature, and velocity of the gas. These physical quantities have a significant impact on the spraying process, and they should be different if the end position of chemical reaction is different. Owing to the lack of reliable testing methods and experimental data, it is difficult to determine the position at which the chemical reaction in the actual spray gun is completed. The chemical reaction in the combustion chamber is governed by several factors. Among them, as a porous medium, the ceramic sheet has a considerable effect on the combustion reaction [37,38,39], but its influence mechanism needs to be further examined.
According to the flame characteristics of actual spray guns, the chemical reaction may well be completed inside the spray gun. Therefore, it is assumed in this study that the combustion reaction is completed inside the combustion chamber. Through a large number of simulation experiments, it was found that the combustion reaction may not be completed quickly in an extremely small space in the combustion chamber. This may be attributed to the fact that the existing chemical reaction models are governed by the generation, movement, and dissipation of turbulence [40,41,42]. The actual turbulence might be much more complicated than the simulated turbulence because the smooth wall surface used in the simulation cannot produce an identical vortex as the microscopic uneven surface of the actual spray gun. To resolve the above issues, the porous ceramic sheet is used in the simulation to examine its influence on the combustion, and radial compressed air inlets are added to the wall of the combustion chamber to enhance the vortexes, thereby promoting the chemical reaction. Further, a full 3D model is established to characterize the ceramic sheet and the radial inlets. Four kinds of meshes are used to analyze the influence of the type of mesh cells, ceramic sheet, and radial inlets on the flow field. The commercial software Fluent is used to simulate the combustion and species transportation process, as well as the distribution of physical quantities such as velocity, temperature, and composition. The effect of chemical reactions on the flow field is examined through two chemical reactions (one-step method and two-step method). The present research work tried to provide theoretical reference for understanding the mechanism of the AC-HVAF spraying process by a simulation study.

2. Configuration of the Model

2.1. Model Structure

The model structure and spraying mechanism are shown in Figure 1. When the spray gun is working, the compressed air and propane are separately injected into the combustion chamber in a discrete form. After ignition of the AC-HVAF gun, the ceramic sheet is preheated up to the auto-ignition temperature of the mixture gas; then, the high temperature ceramic sheet is used as catalytic medium to continuously ignite the mixture, which is called activate combustion. The combustion products enter the rectifying chamber through hundreds of orifices of ceramic sheets and are axially injected with the nitrogen pre-mixed powder particles, which are accelerated in the nozzle to hit the substrate and form a coating.
The model structure mainly includes gas inlets, combustion chamber, rectifying chamber, porous ceramic sheet, and Laval nozzle. The most significant difference between the model and the actual spray gun structure is the additional radial compressed air inlets. All gases are discretely injected into the combustion chamber through 74 small holes (compressed air inlets: 16 in the axial direction, 40 in the radial direction; propane inlets: 18 in the axial direction).
The square holes are evenly distributed on the ceramic sheet, and there is a hollow circular hole with a diameter of 12 mm at the axial center of the combustion chamber and the ceramic sheet, which is the installation position of the powder injector for conveying powder particles and carrier gas (N2). The detailed dimensions of the model are shown in Figure 2.

2.2. Computational Domain and Boundary Conditions

As shown in Figure 2, the flow field includes internal flow field and external flow field. In the internal flow field, the mass flow inlets are adopted as the gas inlets. Here, the temperatures of the compressed air and propane are assumed to be 320 K and 340 K, respectively (propane is preheated). The wall temperature of the spray gun is assumed to be 600 K due to the cooling effect of the compressed air. As the ceramic sheet has high heat capacity, high wall temperature is employed, which has a tendency of linear increase along the axis from 1500 K to 1800 K.
The external flow field is assumed to be a cylinder with a diameter of 400 mm to ensure the reliability of the simulation results. The pressure outlet is located at the cylindrical surface, and is set to 1 atm. The left boundary is defined as the velocity inlet of the external flow field, and the velocity decreases gradually from the muzzle to the pressure outlet, i.e., the velocity increases with decrease in the distance from the jet center. The maximum and minimum velocities are 8 m/s and 0.1 m/s, respectively. The temperature of the ceramic sheet and the velocity at the velocity inlet of the external flow field change with the position, which is realized by the following custom functions:
T ( x ) = ( T m a x T m i n ) · x x m i n x m a x x m i n + T m i n
u ( r ) = u m a x · ( 1 r r m i n r m a x r m i n )
T m a x = 1800   K ,   T m i n = 1500   K ,   x m a x = 24   mm ,   x m i n = 20   mm u m a x = 8   m / s ,   u m i n = 0.1   m / s ,   r m a x = 200   mm ,   r m i n = 4   mm .
where T and x represent the temperature and the absolute value of the x-axis coordinate of each point on the ceramic sheet, respectively; and u and r represent the velocity and the absolute value of the radial coordinate in the left boundary of the outflow field, respectively. The corresponding angular symbols represent the maximum and minimum value, respectively.

2.3. Meshing

To examine the effect of the type of mesh cells, ceramic sheet, the size of ceramic holes, and radial inlets on the simulation results, four types of meshes are used, and the detailed parameters are listed in Table 1. Figure 3 shows the layout of the third meshes, which consist of 1,082,096 hexahedral cells and 1,149,140 nodes, and the mesh independence is tested. As mentioned earlier, the small holes in the model such as gas inlets and the ceramic holes are considered to be square holes (different from the round holes in the actual situation) for improving the quality of the meshes and reducing the number of meshes. Although the geometry is not completely consistent with the actual spray gun, it is verified that there is almost no difference with the simulation results based on round holes.

2.4. Flow Dynamic Model

The AC-HVAF spraying process is accompanied with complex physical, chemical, metallurgical, and species transport processes, where the gas composition, momentum, and energy transport are the primary difficult points of simulation. From the perspective of fluid mechanics, as a continuous phase, gas follows the conservation of mass, momentum, and energy. Therefore, the governing equations of gas flow include mass conservation equation, momentum conservation equation, energy transport equation, and species transport equation [43,44].
(1)
Mass conservation equation:
ρ t + · ( ρ v ) = 0
(2)
Momentum conservation equation:
t ( ρ v ) + · ( ρ v v ) = p + · ( τ ̿ e f f ) + · ( ρ v v ¯ )
The stress tensor τ ̿ e f f is given by
τ ̿ e f f = μ e f f [ ( v + v T ) 2 3 v I ]
(3)
Energy transport equation:
t ( ρ E ) + · [ v ( ρ E + p ) ] = · [ k e f f T j h i J i + ( τ ̿ e f f · v ) ] + S h
(4)
Species transport equations:
t ( ρ Y i ) + · ( ρ v Y i ) = · J i + R i
where ρ is the density of airflow,   p is the static pressure, v is the velocity, ρ v v ¯ is the Reynolds stress, μ e f f is the effective viscosity, I is the unit tensor, E is the enthalpy, k e f f is the effective conductivity, T is the temperature, J i is the diffusion flux of species i , h i is the sensible enthalpy of species i , S h is the heat generation rate from chemical reactions, Y i is the mass fraction of species i , and R i is the net rate of production of species i by chemical reaction.
The realizable k-ε model based on the standard k-ε model is selected to characterize the turbulent motion of airflow, where k is the turbulent kinetic energy and ε is the rate of dissipation of turbulent energy. The realizable k-ε model modifies the standard k-ε model in two important ways: it contains an alternative formulation for the turbulent viscosity; and a modified transport equation for the dissipation rate, ε, has been derived from an exact equation for the transport of the mean-square vorticity fluctuation [45]. Therefore, compared with the standard k-ε model, the realizable k-ε model has advantages in the following simulation: free flows including jets and mixing layers, channel and boundary layer flows, and separated flows. Especially noteworthy is the fact that the realizable k-ε model resolves the round-jet anomaly [46,47,48,49]. The local values of k and ε are obtained by solving the following equations [45]:
t ( ρ k ) + X j ( ρ k u j ) = X j [ ( μ + μ t σ k ) k X j ] + G k + G b ρ ε Y M + S k  
t ( ρ ε ) + X j ( ρ ε u j ) = X j [ ( μ + μ t σ ε ) k X j ] + ρ C 1 S ε ρ C 2 ε 2 k + v ε + C 1 ε ε k C 3 ε G b + S ε
G b = β g i μ t P r t T x i
β = 1 ρ ( ρ T ) p
For ideal gases, Equation (9) reduces to
G b = g i μ t ρ P r t T x i
C 1 = max [ 0.43 , η η + 5 ] ,   η = S k ε ,   S = 2 S i j S i j ,   P r t = 0.85
Here, G k represents the generation of turbulence kinetic energy due to the mean velocity gradients. G b is the generation of turbulence kinetic energy due to buoyancy. Y M represents the contribution of the fluctuating dilatation in compressible turbulence to the overall dissipation rate. P r t is the turbulent Prandtl number for energy and g i is the component of the gravitational vector in the ith direction. C 2 and C 1 ε are constants. σ k and σ ε are the turbulent Prandtl numbers for k and ε , respectively. S k and S ε are user-defined source terms.

2.5. Combustion Model

The Eddy dissipation model (EDM) is used for the combustion reaction between air and propane. In this model, it is assumed that the rate of chemical reaction is higher than that of the turbulent mixing, which implies that the reaction progress is determined by the mixing rate between fuel and oxidant. In fact, the reaction process of propane and oxygen is rather complicated, which involves much intermediate products such as CO and NOx. However, according to the earlier reported simulation results [16], the concentration of other intermediate products in the full flow field is extremely low, except for CO. To simplify the calculation, the reaction between propane and oxygen is simulated by the one-step method and two-step method [50].
One-step method:
  C 3 H 8 + 5 O 2 3 CO 2 + 4 H 2 O
Two-step method:
2 C 3 H 8 + 7 O 2 6 CO + 8 H 2 O
2 CO + O 2 2 CO 2
Theoretically, three situations may occur in the chemical reaction: excess air, excess fuel, and stoichiometric numbers that meet the reaction formula. Excess air can cause oxidation of the metal parts of the spray gun and particles; hence, this situation need not be discussed as it has little practical significance. The slight excess fuel has certain advantages because the intermediate product CO can prevent the oxidation of particles.
To analyze the influences of two chemical reactions on the flow field, the parameters listed in Table 2 are examined, and the other parameters remain constant. In the experiment name, the Roman numerals (e.g., I in I-9.5) denote the chemical reaction mode, and the following numbers (e.g., 9.5 in I-9.5) denote the mass flow ratio of air to fuel.

3. Results and Discussion

3.1. Influence of Mesh Cells Pattern on Simulations

Figure 4 shows the temperature fields of the combustion chamber and the substrate formed using different mesh cells. Comparing Figure 4a,c it can be seen that the temperature field obtained using the fourth meshes has better symmetry than that using the first meshes. After comprehensive comparison of Figure 4a-d, it becomes clear that the asymmetry of temperature field in the combustion chamber affects the subsequent flow field when the first meshes are adopted, and this effect decreases with the large-gradient change in the subsequent flow field. In fact, the other physical quantities such as velocity and composition in the combustion chamber are also distributed asymmetrically when the first meshes are adopted.
Obviously, the model is structurally symmetrical, and the flow field should also be at least somewhat symmetrical. Hence, compared with the hexahedral cells, the simulation results exhibit more errors when tetrahedral cells are adopted. The fourth meshes are structured meshes, whose direction is highly consistent with the fluid flow direction, while tetrahedral cells of the first meshes cannot be aligned with the flow direction. The inconsistency between the meshes’ direction and the flow direction leads to the false-diffusion, which affects the accuracy of the simulation results. The false-diffusion is small when low-velocity flow is simulated, which has a minor impact on the simulation results, but its effect cannot be ignored if the flow rate increases [51]. The airflow velocity of AC-HVAF spraying is rather high, especially the supersonic airflow behind the nozzle throat. Therefore, the false-diffusion caused by hexahedral mesh cells is smaller than that caused by tetrahedral mesh cells in the supersonic flow field, and the former provides higher simulation accuracy.

3.2. Effect of Ceramic Sheet on Gas Flow Field

Figure 5 and Figure 6 show the comparisons of velocity and temperature with and without a ceramic sheet in the case wherein the other parameters are the same. It is clear that the ceramic sheet has a rectification effect on the airflow. The velocity and temperature of the airflow through the ceramic sheet are more uniform. Further, it can be seen in Figure 5 that the ceramic sheet increases the vortexes in the combustion chamber and prevents the fluid element with high radial velocity from passing, which increases the residence time and stroke of these fluid elements in the combustion chamber as well as the reaction probability between the fuel and the oxidant. This reduces the production of harmful gases such as CO and NOx that are discharged from the gun, which leads to efficient combustion and reduces the energy consumption.

3.3. Effect of Hole Sizes of Ceramic Sheet on Gas Flow Field

The turbulent kinetic energy k is given by the following [52]:
k = 1 2 ( u 2 + v 2 + w 2 )
In Equation (14), the right-hand side of the equation represents one half of the sum of velocity variances. The greater the turbulent kinetic energy, the greater the velocity inhomogeneity, resulting in more energy consumption between the airflows.
Figure 7 shows the difference of the turbulent kinetic energy inside the spray gun of three different configurations. As shown in this figure, the ceramic sheet effectively reduces the turbulent kinetic energy inside the spray gun, and the smaller the holes in the ceramic sheet, the smaller the turbulent kinetic energy after the region of the ceramic sheet. The lower turbulent kinetic energy of the flow is conducive to reducing the internal kinetic energy consumption of the airflow, so the higher airflow velocity appears at the muzzle, which is demonstrated by the curves in Figure 8b.

3.4. Effect of the Radial Air Inlets on Gas Flow Field

The flow field of the model is completely different if there are no radial inlets. The experimental results in Figure 8 are obtained according to I-9.5. In Figure 8a, the radial inlets are closed and their corresponding mass flow is added to the axial air inlets. As shown in the figure, the combustion reaction occurs over larger area when there are no radial inlets, which results in a more uneven temperature field inside the gun. Figure 8b shows that the airflow velocity is related to both the ceramic sheet and the radial inlets. When there are radial inlets, the smaller the hole in the ceramic sheet, the higher the jet velocity of the gun; however, when there are no radial inlets, the jet velocity of the gun is lower even if there is a ceramic sheet.
It is noteworthy that the radial air inlets are close to the ceramic sheet in the all models; the low-temperature compressed air that enters the combustion chamber is heated by the high-temperature ceramic sheet, and its temperature rises rapidly, which is extremely beneficial for promoting chemical reactions. Combining Figure 5, Figure 6, Figure 7 and Figure 8, it is clear that the gas from radial air inlets increases the radial velocity of the fluid element in front of the ceramic sheet, so they have a low probability of passing through the ceramic sheet. In other words, the radial inlets not only promote combustion, but also enhance the rectification effect of the ceramic sheet. As mentioned above, this can promote the consistency between airflow direction and meshes’ direction, reduce the false-diffusion, and improve the accuracy of the simulation results.

3.5. Influence of Different Chemical Reactions on the Simulation Results

3.5.1. Influence of the Different Chemical Reactions on Velocity Field

The velocity of the airflow is the most important indicator of AC-HVAF spraying, and it determines the velocity and heated history of the particles. For the convenience of discussion, the flow field is divided into three zones along the axial direction: zone I indicates the internal flow field, zone II represents the shock zone, and zone III is the region from the end position of shock wave to the substrate. Figure 9 shows an example of these three zones, and the partial contours of the airflow velocity under eight experimental schemes. Table 3 lists the maximum velocity for these experimental schemes. Figure 10 clearly shows the variation in velocity along the axis. Here, it is clear that the velocity change process of each experimental scheme is basically the same.
In zone I, the airflow in the nozzle continues to accelerate and reaches the maximum velocity until the outlet. As shown in Figure 10 and Table 2, the velocities of different experimental schemes have a negligible difference in zone I, which indicates that the different chemical reactions have a minor effect on the velocity in the internal flow field.
In zone II, the airflow is in an over-expanded state, and some diamond shock waves are formed because the pressure at the outlet is lower than the ambient pressure. Because of the impact of the shock waves, the velocity increases and decreases periodically. The momentum exchange between the jet and the surrounding air leads to a gradual decrease in the average velocity and the amplitude of the airflow. The difference is observed at the end positions of the shock zone only between I-9.5 and II-9.5. The end position of the shock zone for I-9.5 is later than that of II-9.5. This indicates that the airflow of I-9.5 has more stability of the velocity, and its high velocity can be maintained for a long distance.
In zone III, the velocity of each experiment drops rapidly, and the difference between them lies in decline rate of the velocity. However, there is no obvious difference between the two chemical reactions.

3.5.2. Influence of the Different Chemical Reactions on Temperature Field

Temperature field is another important factor affecting the coating quality. The temperature distribution near the axis of the spray gun governs the heated history of the particles. However, the temperature field near the substrate influences the final quality of the coating deeply.
Figure 11 shows the change of airflow temperature along the axis. The negligible difference in temperature between different chemical reactions can be seen in all zones except for I-9.5 and II-9.5. However, the differences between them are limited to zone III. As shown in Figure 11, the temperature fields near substrates in zone III show an obvious difference in different experiments, which is mainly attributed to two factors: different degrees of compression effect and released heat by the combustion between the remaining fuel and the oxygen from the external. Figure 12 provides a good proof of this phenomenon.

3.5.3. Influence of the Different Chemical Reactions on Gas Species Transport

Figure 12a,b show that the degree of fuel excess of each experiment, and Figure 12c implies that the chemical reaction basically ends in the combustion chamber, which satisfies the previous assumption. Obviously, the reaction products of different chemical reactions are different because of the excessive fuel. If the two-step method is used, the excess C3H8 generates CO, so the average density of the gas is lower than the one-step method. The higher the airflow density, the higher the momentum at the same velocity, and the state of the airflow motion is less susceptible to interference by the air from the external flow field. Therefore, the difference in the density of the airflow causes a significant velocity variation after the shock wave, which is the reason that the airflow of I-9.5 has the more stability of the velocity than the one of II-9.5.
Figure 12c shows that the oxygen from the external field rapidly enters the center of the airflow through diffusion and results in a rapid increase of the jet oxygen content. This may lead to the oxidation of the metal powder and increase of the oxygen content in the coating, and decrease the quality of the coating. However, this can only occur when the gas inlets flow rate is insufficient, such as I-15.3-S and II-15.3-S. It should be noted that the combustion reaction of excess fuel is very violent near the substrate because the oxygen from the external field is carried by the jet and concentrates in a small area. Therefore, the temperature of the gas can rise rapidly near the substrate, which will cause the substrate to form a very large temperature gradient (e.g., up to 200 K/mm for I-11.2), as shown in Figure 11.
According to the stoichiometry in the chemical reaction equation, in order to fulfill a complete combustion process, the mass flow ratio of air to fuel should be about 15.6. The simulation results of the two chemical reactions are relatively close when the degree of fuel excess is less than 24% (corresponding to I-11.2 and II-11.2). However, the simulation results were significantly different when the degree of fuel excess was 68% (corresponding to I-9.8 and II-9.8).
In a word, on the premise that the combustion can end in the combustion chamber, the influence of chemical reaction method on the simulation results is significant only when the fuel excess degree is large.

4. Conclusions

In this study, a complete 3D model was established using an actual spray gun as a prototype. It is assumed that the combustion reaction is completed inside the combustion chamber. The influence of the mesh cells’ type on the simulation results was analyzed, and the effects of porous ceramic sheet, radial inlets, and chemical reaction process were thoroughly investigated. The main results of the study are summarized as follows:
  • The false-diffusion caused by hexahedral mesh cells is smaller than that caused by tetrahedral mesh cells in the supersonic flow field, and the former provides higher simulation accuracy.
  • The porous ceramic sheet has a rectification effect. The radial inlet and porous ceramic sheets increase the residence time and stroke of the gas in the combustion chamber, and increase the probability of chemical reactions.
  • Under the assumption that the combustion reaction ends in the spray combustion chamber, the influence of chemical reaction method on the simulation results is significant only when the fuel is excessive.

Author Contributions

Conceptualization, F.L. and Z.L.; methodology, F.L.; software, F.L.; validation, F.L., Z.L., and M.F.; formal analysis, F.L.; investigation, H.H.; resources, M.F.; data curation, F.L.; writing—original draft preparation, H.H.; writing—review and editing, F.L.; visualization, H.H.; supervision, F.L.; project administration, Z.L.; funding acquisition, Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (51805492).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Berger, L.; Puschmann, R.; Spatzier, J.; Matthews, S. Potenzial von HVAF-Spritzprozessen; Potential of HVAF Spray Processes. Therm. Spray Bull. 2013, 6, 16–20. Available online: https://www.researchgate.net/publication/257656151_Potenzial_von_HVAF-Spritzprozessen_Potential_of_HVAF_Spray_Processes (accessed on 30 January 2021).
  2. Sadeghimeresht, E.; Markocsan, N.; Nylén, P. Microstructural characteristics and corrosion behavior of HVAF- and HVOF-sprayed Fe-based coatings. Surf. Coat. Technol. 2017, 318, 365–373. [Google Scholar] [CrossRef]
  3. Hulka, I.; Serban, V.; Secoşan, I.; Vuoristo, P.; Niemi, K. Wear properties of CrC–37WC–18M coatings deposited by HVOF and HVAF spraying processes. Surf. Coat. Technol. 2012, 210, 15–20. [Google Scholar] [CrossRef]
  4. Milanti, A.; Koivuluoto, H.; Vuoristo, P. Influence of the Spray Gun Type on Microstructure and Properties of HVAF Sprayed Fe-Based Corrosion Resistant Coatings. J. Therm. Spray Technol. 2015, 24, 1312–1322. [Google Scholar] [CrossRef]
  5. Sadeghi, E.; Reddy, L.; Hussain, T.; Huhtakangas, M.; Markocsan, N.; Joshi, S. Influence of KCl and HCl on high temperature corrosion of HVAF-sprayed NiCrAlY and NiCrMo coatings. Mater. Des. 2018, 148, 17–29. [Google Scholar] [CrossRef]
  6. Jafari, R.; Sadeghi, E.; Farahani, T.S.; Huhtakangas, M.; Markocsan, N.; Joshi, S. KCl-Induced High-Temperature Corrosion Behavior of HVAF-Sprayed Ni-Based Coatings in Ambient Air. J. Therm. Spray Technol. 2018, 27, 500–511. [Google Scholar] [CrossRef] [Green Version]
  7. Liu, Y.; Liu, W.; Ma, Y.; Meng, S.; Liu, C.; Long, L.; Tang, S. A comparative study on wear and corrosion behaviour of HVOF- and HVAF-sprayed WC–10Co–4Cr coatings. Surf. Eng. 2016, 33, 1–9. [Google Scholar] [CrossRef]
  8. Janka, L.; Norpoth, J.; Trache, R.; Thiele, S.; Berger, L.-M. HVOF- and HVAF-Sprayed Cr3C2-NiCr Coatings Deposited from Feedstock Powders of Spherical Morphology: Microstructure Formation and High-Stress Abrasive Wear Resistance Up to 800 °C. J. Therm. Spray Technol. 2017, 26, 1720–1731. [Google Scholar] [CrossRef]
  9. Zeng, Z.; Sakoda, N.; Tajiri, T.; Kuroda, S. Structure and corrosion behavior of 316L stainless steel coatings formed by HVAF spraying with and without sealing. Surf. Coat. Technol. 2008, 203, 284–290. [Google Scholar] [CrossRef]
  10. Gorlach, I. A new method for thermal spraying of Zn–Al coatings. Thin Solid Films 2009, 517, 5270–5273. [Google Scholar] [CrossRef]
  11. Sadeghi, E.; Markocsan, N.; Nylén, P. Microstructural and electrochemical characterization of Ni-based bi-layer coatings produced by the HVAF process. Surf. Coat. Technol. 2016, 304, 606–619. [Google Scholar] [CrossRef]
  12. Liang, D.; Ma, J.; Cai, Y.; Liu, X.; Xie, S.; Wei, X.; Xu, G.; Shen, J. Characterization and elevated-temperature tribological performance of AC–HVAF-sprayed Fe-based amorphous coating. Surf. Coat. Technol. 2020, 387, 125535. [Google Scholar] [CrossRef]
  13. Fauchais, P.; Vardelle, A. Thermal Sprayed Coatings Used Against Corrosion and Corrosive Wear. Adv. Plasma Spray Appl. 2012, 10, 34448. [Google Scholar] [CrossRef] [Green Version]
  14. Gorlach, I.A. Thermal Stress Evaluation of Thermo-Blast Jet Nozzle Materials. Ph.D. Thesis, North-West University, Potchefstroom, South Africa, July 2004. Available online: http://repository.nwu.ac.za/handle/10394/352 (accessed on 30 January 2021).
  15. Yuan, X.; Wang, H.; Hou, G.; Zha, B. Numerical Modeling of a Low Temperature High Velocity Air Fuel Spraying Process with Injection of Liquid and Metal Particles. J. Therm. Spray Technol. 2006, 15, 413–421. [Google Scholar] [CrossRef]
  16. Jiang, H.-R.; Li, M.-L.; Wei, X.-S.; Ma, T.-C.; Dong, Y.; Ying, C.-X.; Liao, Z.-Y.; Shen, J. Numerical Investigation of In-Flight Behavior of Fe-Based Amorphous Alloy Particles in AC-HVAF Thermal Spray Process. J. Therm. Spray Technol. 2019, 28, 1146–1159. [Google Scholar] [CrossRef]
  17. Gao, X.; Li, C.; Zhang, D.; Gao, H.; Han, X. Numerical analysis of the activated combustion high-velocity air-fuel (AC-HVAF) thermal spray process: A survey on the parameters of operation and nozzle geometry. Surf. Coat. Technol. 2021, 405, 126588. [Google Scholar] [CrossRef]
  18. Zhou, Z.; Chen, Y.; Hu, Z.; Shen, B.; Liang, X. Numerical Investigation of Particles in Warm-Particle Peening-Assisted High-Velocity Oxygen Fuel (WPPA-HVOF) Spraying. J. Therm. Spray Technol. 2020, 1–13. [Google Scholar] [CrossRef]
  19. Jafari, H.; Emami, S.; Mahmoudi, Y. Numerical investigation of dual-stage high velocity oxy-fuel (HVOF) thermal spray process: A study on nozzle geometrical parameters. Appl. Therm. Eng. 2017, 111, 745–758. [Google Scholar] [CrossRef] [Green Version]
  20. Pan, J.; Hu, S.; Lijun, Y.; Ding, K.; Ma, B. Numerical analysis of flame and particle behavior in an HVOF thermal spray process. Mater. Des. 2016, 96, 370–376. [Google Scholar] [CrossRef]
  21. Tabbara, H.; Gu, S. Computational simulation of liquid-fuelled HVOF thermal spraying. Surf. Coat. Technol. 2009, 204, 676–684. [Google Scholar] [CrossRef] [Green Version]
  22. Kamnis, S.; Gu, S. Study of In-Flight and Impact Dynamics of Nonspherical Particles from HVOF Guns. J. Therm. Spray Technol. 2009, 19, 31–41. [Google Scholar] [CrossRef] [Green Version]
  23. Vengeepuram Rajagopalan, S. Numerical analysis of a High Velocity Oxygen Fuel (HVOF) Process. Master’s Thesis, Concordia University, Montreal, QC, Canada, August 2007. Available online: https://spectrum.library.concordia.ca/975504/ (accessed on 13 July 2020).
  24. Kamnis, S.; Gu, S. Numerical modelling of propane combustion in a high velocity oxygen–fuel thermal spray gun. Chem. Eng. Process. Process. Intensif. 2006, 45, 246–253. [Google Scholar] [CrossRef]
  25. Kamnis, S.; Gu, S. 3-D modelling of kerosene-fuelled HVOF thermal spray gun. Chem. Eng. Sci. 2006, 61, 5427–5439. [Google Scholar] [CrossRef]
  26. Cheng, D.; Trapaga, G.; McKelliget, J.W.; Lavernia, E.J. Mathematical modelling of high velocity oxygen fuel thermal spraying of nanocrystalline materials: An overview. Modelling Simul. Mater. Sci. Eng. 2003, 11, R1–R31. Available online: https://iopscience.iop.org/article/10.1088/0965-0393/11/1/201/meta (accessed on 30 January 2021). [CrossRef]
  27. Dolatabadi, A.; Mostaghimi, J.; Pershin, V. Effect of a cylindrical shroud on particle conditions in high velocity oxy-fuel spray process. Sci. Technol. Adv. Mat. 2002, 3, 245–255. [Google Scholar] [CrossRef]
  28. Eastwick, C.N.; Simmons, K.; McCartney, D.G. Computational Fluid Dynamic Modeling of Gas Flow Characteristics in a High-Velocity Oxy-Fuel Thermal Spray System. J. Therm. Spray Technol. 2001, 10, 461–469. [Google Scholar] [CrossRef]
  29. Hassan, B.; Lopez, A.R.; Oberkampf, W.L. Computational analysis of a three-dimensional High-Velocity Oxygen-Fuel (HVOF) Thermal Spray torch. J. Therm. Spray Technol. 1998, 7, 71–77. Available online: https://link.springer.com/article/10.1361/105996398770351061 (accessed on 30 January 2021). [CrossRef]
  30. Eidelman, S.; Yang, X. Three dimensional simulation of HVOF spray deposition of nanoscale materials. Nanostructured Mater. 1997, 9, 79–84. [Google Scholar] [CrossRef]
  31. Yang, X.; Eidelman, S. Numerical analysis of a high-velocity oxygen-fuel thermal spray system. J. Therm. Spray Technol. 1996, 5, 175–184. [Google Scholar] [CrossRef]
  32. Oberkampf, W.L.; Talpallikar, M. Analysis of a high-velocity oxygen-fuel (HVOF) thermal spray torch part 1: Numerical formulation. J. Therm. Spray Technol. 1996, 5, 53–61. [Google Scholar] [CrossRef]
  33. Li, H.-P.; Pfender, E. Three Dimensional Modeling of the Plasma Spray Process. J. Therm. Spray Techn. 2007, 16, 245–260. Available online: https://link.springer.com/article/10.1007/s11666-007-9023-x (accessed on 30 January 2021). [CrossRef] [Green Version]
  34. Xiong, H.-B.; Zheng, L.-L.; Sampath, S.; Williamson, R.; Fincke, J.R. Three-dimensional simulation of plasma spray: Effects of carrier gas flow and particle injection on plasma jet and entrained particle behavior. Int. J. Heat Mass Transf. 2004, 47, 5189–5200. [Google Scholar] [CrossRef]
  35. Ahmed, I.; Bergman, T.L. Three-Dimensional Simulation of Thermal Plasma Spraying of Partially Molten Ceramic Agglomerates. J. Therm. Spray Technol. 2000, 9, 215–224. [Google Scholar] [CrossRef]
  36. Samareh, B.; Dolatabadi, A. A Three-Dimensional Analysis of the Cold Spray Process: The Effects of Substrate Location and Shape. J. Therm. Spray Technol. 2007, 16, 634–642. [Google Scholar] [CrossRef]
  37. Trimis, D.; Durst, F. Combustion in a Porous Medium-Advances and Applications. Combust. Sci. Technol. 1996, 121, 153–168. [Google Scholar] [CrossRef]
  38. Kamal, M.M.; Mohamad, A.A. Combustion in Porous Media. Proc. Inst. Mech. Eng. Part A: J. Power Energy 2006, 220, 487–508. [Google Scholar] [CrossRef]
  39. Mujeebu, M.A.; Abdullah, M.; Abu Bakar, M.; Mohamad, A. Applications of porous media combustion technology—A review. Appl. Energy 2009, 86, 1365–1375. [Google Scholar] [CrossRef]
  40. Kärrholm, F.P. Numerical Modelling of Diesel Spray Injection, Turbulence Interaction and Combustion; Chalmers University of Technology: Gothenburg, Sweden, 2008; Available online: https://www.researchgate.net/profile/Fabian_Peng_Kaerrholm/publication/230605808_Numerical_Modelling_of_Diesel_Spray_Injection_Turbulence_Interaction_and_Combustion/links/5b8f735492851c6b7ec054d8/Numerical-Modelling-of-Diesel-Spray-Injection-Turbulence-Interaction-and-Combustion.pdf (accessed on 30 January 2021).
  41. Giordano, P.; Lentini, D. Combustion-Radiation-Turbulence Interaction Modeling in Absorbing/Emitting Nonpremixed Flames. Combust. Sci. Technol. 2001, 172, 1–22. [Google Scholar] [CrossRef]
  42. Bilger, R.; Pope, S.; Bray, K.; Driscoll, J. Paradigms in turbulent combustion research. In Proceedings of the Combustion Institute; Elsevier: Amsterdam, The Netherlands, 2005; Volume 30, pp. 21–42. [Google Scholar]
  43. Anderson, J.D. Governing Equations of Fluid Dynamics. In Computational Fluid Dynamics; Springer Nature: London, UK, 2009; pp. 15–51. [Google Scholar]
  44. ANSYS, Fluent (2020 R1), Theory Guide. Available online: https://ansyshelp.ansys.com/account/secured?returnurl=/Views/Secured/corp/v201/en/flu_th/flu_th_sec_turb_realizke.html (accessed on 31 January 2021).
  45. Shih, T.-H.; Liou, W.W.; Shabbir, A.; Yang, Z.; Zhu, J. A new k-ϵ eddy viscosity model for high reynolds number turbulent flows. Comput. Fluids 1995, 24, 227–238. [Google Scholar] [CrossRef]
  46. Imanian, H.; Mohammadian, A. Numerical simulation of flow over ogee crested spillways under high hydraulic head ratio. Eng. Appl. Comput. Fluid Mech. 2019, 13, 983–1000. [Google Scholar] [CrossRef] [Green Version]
  47. Shaheed, R.; Mohammadian, A.; Gildeh, H.K. A comparison of standard k–ε and realizable k–ε turbulence models in curvedand confluent channels. Environ. Fluid Mech. 2019, 19, 543–568. [Google Scholar] [CrossRef]
  48. Kucukgokoglan, S.; Aroussi, A.; Pickering, S.; Menacer, M. Prediction of interactions between burners in multi-burner systems. In Proceedings of the International Symposium on Multiphase Flow and Transport Phenomena, Antalya, Turkey, 3 November 1997; Begel House Inc.: Antalya, Turkey, 1997. Available online: http://old.ichmt.org/abstracts/MFTP-2000/Abstracts/71.pdf (accessed on 31 January 2021).
  49. Reynolds, W. Fundamentals of Turbulence for Turbulence Modeling and Simulation; Stanford Univ CA Dept of Mechanical Engineering: Palo Alto, CA, USA, 1987; Available online: https://apps.dtic.mil/sti/citations/ADP005793 (accessed on 31 January 2021).
  50. Kundu, K.; Penko, P.; Yang, S. Reduced reaction mechanisms for numerical calculations in combustion of hydrocarbon fuels. In Proceedings of the 36th AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 12–15 January 1998; p. 803. [Google Scholar] [CrossRef]
  51. Patankar, S.V. Numerical Heat Transfer and Fluid Flow; Hemisphere Publishing Corp: Washington, DC, USA, 1980; Volume 210, p. 108. Available online: https://books.google.de/books?hl=zh-CN&lr=&id=Y2G1DwAAQBAJ&oi=fnd&pg=PP1&dq=Numerical+Heat+Transfer+and+Fluid+Flow&ots=7TF54WA5UJ&sig=6e6njvYn8gb5GkAY4a61JmJo80s&redir_esc=y#v=onepage&q=Numerical%20Heat%20Transfer%20and%20Fluid%20Flow&f=false (accessed on 31 January 2021).
  52. Harsha, P.T.; Lee, S.C. Correlation between turbulent shear stress and turbulent kinetic energy. AIAA J. 1970, 8, 1508–1510. [Google Scholar] [CrossRef]
Figure 1. Structural schematic of the activated combustion high-velocity air-fuel (AC-HVAF) spray gun model and the spraying principle: (a) left view of the spray gun; (b) ceramic sheet and powder injector mounting hole; (c) powder injector of actual spray gun; (d) propane distributor of actual spray gun; (e) porous ceramic sheet of actual spray gun.
Figure 1. Structural schematic of the activated combustion high-velocity air-fuel (AC-HVAF) spray gun model and the spraying principle: (a) left view of the spray gun; (b) ceramic sheet and powder injector mounting hole; (c) powder injector of actual spray gun; (d) propane distributor of actual spray gun; (e) porous ceramic sheet of actual spray gun.
Materials 14 00657 g001
Figure 2. Model size and boundary conditions (unit: mm).
Figure 2. Model size and boundary conditions (unit: mm).
Materials 14 00657 g002
Figure 3. Configuration of the third meshes: (a) internal flow field and (b) external flow field.
Figure 3. Configuration of the third meshes: (a) internal flow field and (b) external flow field.
Materials 14 00657 g003
Figure 4. Temperature fields formed using different mesh cells. (a,c) Temperature field at the middle section of the combustion chamber; (b,d) temperature field on the surface of the substrate. (a,b) The fourth meshes are adopted; (c,d) the first meshes are adopted (the pictures are scaled to different degrees).
Figure 4. Temperature fields formed using different mesh cells. (a,c) Temperature field at the middle section of the combustion chamber; (b,d) temperature field on the surface of the substrate. (a,b) The fourth meshes are adopted; (c,d) the first meshes are adopted (the pictures are scaled to different degrees).
Materials 14 00657 g004
Figure 5. Internal velocity vector diagram of the spray gun: (a) second meshes and (b) fourth meshes.
Figure 5. Internal velocity vector diagram of the spray gun: (a) second meshes and (b) fourth meshes.
Materials 14 00657 g005
Figure 6. Temperature distribution inside the spray gun: (a) second meshes and (b) fourth meshes.
Figure 6. Temperature distribution inside the spray gun: (a) second meshes and (b) fourth meshes.
Materials 14 00657 g006
Figure 7. Contours of turbulent kinetic energy inside the spray gun: (a) second meshes, (b) third meshes, and (c) fourth meshes.
Figure 7. Contours of turbulent kinetic energy inside the spray gun: (a) second meshes, (b) third meshes, and (c) fourth meshes.
Materials 14 00657 g007
Figure 8. Contours of temperature without radial air inlets (a). The velocity at the central axis of the model with no radial inlets, big holes, and small holes in the ceramic sheet (b).
Figure 8. Contours of temperature without radial air inlets (a). The velocity at the central axis of the model with no radial inlets, big holes, and small holes in the ceramic sheet (b).
Materials 14 00657 g008
Figure 9. The flame of the actual spray gun and the partial contours of the airflow velocity: (a) actual spray gun, (bi) experimental scheme corresponding to Table 1.
Figure 9. The flame of the actual spray gun and the partial contours of the airflow velocity: (a) actual spray gun, (bi) experimental scheme corresponding to Table 1.
Materials 14 00657 g009
Figure 10. Velocity at the axis of the model (a) and a partially enlarged view of the red areas (b).
Figure 10. Velocity at the axis of the model (a) and a partially enlarged view of the red areas (b).
Materials 14 00657 g010
Figure 11. Temperature at the axis of the model.
Figure 11. Temperature at the axis of the model.
Materials 14 00657 g011
Figure 12. Mass fraction of gas species at the axis of the model: (a) C3H8, (b) CO, and (c) O2.
Figure 12. Mass fraction of gas species at the axis of the model: (a) C3H8, (b) CO, and (c) O2.
Materials 14 00657 g012
Table 1. Four types of the meshes used in the model.
Table 1. Four types of the meshes used in the model.
Name of MeshWith Ceramic Sheet (Yes/No)Size of Ceramic HolesNumber of Ceramic HolesType of Cells
Combustion ChamberRest of the Model
#1Yes2 × 2 mm2452TetrahedronHexahedron
#2No--HexahedronHexahedron
#3Yes2 × 2 mm2452HexahedronHexahedron
#4Yes1 × 1 mm21472HexahedronHexahedron
Table 2. Experimental program and parameters.
Table 2. Experimental program and parameters.
Experiment NameChemical ReactionMass Flow Rate at Axial Air Inlets (g/s)Mass Flow Rate at Radial Air Inlets (g/s)Mass Flow Rate at Fuel Inlets (g/s)Air-Fuel Mass Flow RatioTotal Mass Flow Rate (g/s)
I-9.5One-step11.0024.203.709.5038.90
II-9.5Two-step11.0024.203.709.5038.90
I-11.2One-step8.8022.602.8011.2034.20
II-11.2Two-step8.8022.602.8011.2034.20
I-15.3One-step11.0024.202.3015.3037.50
II-15.3Two-step11.0024.202.3015.3037.50
I-15.3-SOne-step8.8019.361.8415.3030.00
II-15.3-STwo-step8.0020.001.8315.3029.83
Table 3. Maximum velocity corresponding to different experimental schemes.
Table 3. Maximum velocity corresponding to different experimental schemes.
ExperimentI-9.5II-9.5I-11.2II-11.2I-15.3II-15.3I-15.3SII-15.3S
Maximum Velocity (m/s)15961596160316081668167116991699
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, F.; Li, Z.; Fang, M.; Hou, H. Numerical Analysis of the Activated Combustion High-Velocity Air-Fuel Spraying Process: A Three-Dimensional Simulation with Improved Gas Mixing and Combustion Mode. Materials 2021, 14, 657. https://doi.org/10.3390/ma14030657

AMA Style

Liu F, Li Z, Fang M, Hou H. Numerical Analysis of the Activated Combustion High-Velocity Air-Fuel Spraying Process: A Three-Dimensional Simulation with Improved Gas Mixing and Combustion Mode. Materials. 2021; 14(3):657. https://doi.org/10.3390/ma14030657

Chicago/Turabian Style

Liu, Fuqiang, Zhiyong Li, Min Fang, and Hua Hou. 2021. "Numerical Analysis of the Activated Combustion High-Velocity Air-Fuel Spraying Process: A Three-Dimensional Simulation with Improved Gas Mixing and Combustion Mode" Materials 14, no. 3: 657. https://doi.org/10.3390/ma14030657

APA Style

Liu, F., Li, Z., Fang, M., & Hou, H. (2021). Numerical Analysis of the Activated Combustion High-Velocity Air-Fuel Spraying Process: A Three-Dimensional Simulation with Improved Gas Mixing and Combustion Mode. Materials, 14(3), 657. https://doi.org/10.3390/ma14030657

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

Article Metrics

Back to TopTop