Next Article in Journal
Xerogel-like Materials from Sustainable Sources: Properties and Electrochemical Performances
Next Article in Special Issue
Explosive Burning of a Mechanically Activated Al and CuO Thermite Mixture
Previous Article in Journal
Intelligent Reconfigurable Photovoltaic System
Previous Article in Special Issue
Experimental Study on Ignition Characteristics of RP-3 Jet Fuel Using Nanosecond Pulsed Plasma Discharge
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Shock Wave Propagation and Flame Kernel Morphology in Laser-Induced Plasma Ignition of CH4/O2/N2 Mixture

State Key Laboratory of Multiphase Flow in Power Engineering, Xi’an Jiaotong University, Xi’an 710049, China
*
Author to whom correspondence should be addressed.
Energies 2021, 14(23), 7976; https://doi.org/10.3390/en14237976
Submission received: 29 October 2021 / Revised: 24 November 2021 / Accepted: 25 November 2021 / Published: 29 November 2021
(This article belongs to the Special Issue Challenges and Research Trends of Combustion Mechanism)

Abstract

:
The application of laser ignition in the aerospace field has promising prospects. Based on the constant volume combustion chamber, the laser ignition of CH4/O2/N2 mixture with different initial pressure, different laser energy, different equivalence ratio and different oxygen content has been carried out. The development characteristics of the flame kernel and shock wave under different conditions are analyzed. In addition, the Taylor model and Jones model are also used to simulate the development process of the shock wave, and a new modified model is proposed based on the Jones model. The experimental results show that under pure oxygen conditions, the chemical reaction rate of the mixture is too fast, which makes it difficult for the flame kernel to form the ring and third-lobe structure. However, the ring structure is easier to form with the pressure and laser energy degraded; the flame kernel morphology is easier to maintain at a rich equivalence ratio, which is caused by the influence of the movement of hot air flow and a clearer boundary between the ring and the third-lobe. The decrease of the initial pressure or the increase of the laser energy leads to the increase in shock wave velocity, while the change of the equivalence ratio and oxygen content has less influence on the shock wave.

1. Introduction

Methane/oxygen is one of the most promising propellants of spacecraft, which has the advantages of excellent performance, low price, and low carbon emissions [1]. System security issues in the energy production and aerospace fields cannot be ignored [2,3]. In order to improve the safety and reliability of ignition, it is necessary to study the propagation of the flame kernel and the shock wave in the combustion process of engines.
In recent years, more and more attention has been paid to laser-induced plasma ignition (LIPI) thanks to the great innovation of laser technology. Compared to conventional spark ignition, this new technique enables the high-precision control of ignition timing. Besides, multi-point ignition is feasible due to the flexible ignition position, which reduces the possibility of misfire and ignition delay [4]. Furthermore, the non-intrusive configuration eliminates the radiation and heat transfer caused by the existence of electrodes [5]. Several studies of LIPI have been conducted. Bradley et al. [6] studied the laser ignition of gaseous fuels and recorded the development of the flame kernel and the propagation of the shock wave for the first time using a high-speed camera, and analyzed the formation of the third-lobe with gas dynamics theory in detail. Xu et al. [7] studied the flame kernel of ethanol/oxygen mixture and compared the difference of velocities in different directions. The rarefaction wave generated by plasma caused the change in gas density gradient and pressure near the plasma, interacts with the high-temperature gas and produces local flow, leading to an annular structure from the left and right sides [6,8,9]. The structure created by the rarefaction wave forms the third-lobe, and the leading edge of the structure decays rapidly. In addition, numerical simulation of the annular structure was also performed [10,11,12]. From a gas dynamics perspective, Spiglanin et al. [13] believed that the third-lobe was created by the change of local flow velocity due to excessive expansion. Raizer et al. [14] found that the laser-induced plasma made the ionization front depart from the focus and move in the direction of laser.
Wohlhüter et al. [15] conducted numerical simulation and experimental study on the peak pressure of methane/oxygen mixture using laser ignition in a combustion chamber. Dumitrache et al. [16] studied the flammability limit and minimum ignition energy of methane/air mixture in a rapid compression machine. Optical diagnostics, including the schlieren technique, shadow technique and spectral analysis were used to study the kernel structure and combustion of laser ignition; thereby, the detailed process of flame development could be obtained [17]. Leela et al. [18] used an Intensified Charge Coupled Device (ICCD) camera to study the development of laser-induced plasma in air and observed the self-radiation of plasma. They believed that when the energy is high enough, the leading edge of the laser beam first reached the breakdown threshold and produced plasma and the beam interacted with the trailing edge of the plasma later, creating two kernels. Shen et al. [19] found that the leading and trailing edges of the laser beam focused in two different positions in the direction of propagation, thus created two breakdown points.
In addition to the propagation of flame kernel, the theory of explosive waves when LIPI occurs has been studied. The strong blast wave formed by solar flares will decay into a medium-strength shock wave after reaching the earth. Besides, the strong shock wave formed by natural activities such as lightning discharge will attenuate to a medium-strength shock wave after propagating a distance of 5 mm. In fact, in about 99.5% of the propagation distance, the shock wave propagates at a medium intensity.
Many investigations have studied the generation and propagation of explosive waves. Taylor [9] proposed the similarity solution theory of strong explosion waves in 1951. The study ignored the impact of ionization and excitation on shock waves, and then obtained the relationship between radius, velocity, temperature and pressure of the shock wave [5,8,9], but this theory was only applicable when the ratio of the shock wave pressure (PS) to the ambient pressure (P0) was greater than 10. Then, Taylor, Sakurai and Lin et al. [20,21] respectively proposed the similarity to solve the cylindrical explosion wave theory. Many researchers have verified the accuracy of Taylor’s explosive wave theory through experiments, which can predict the early development of LIPI hot gas nuclei and shock waves. After that, Gatti et al. [22], in 1988, conducted a study on the evolution of LIPI spherical shock waves with time in different medium gases. Based on the Taylor blast wave model, they proposed a simplified relationship for calculating the shock radius and shock Mach number. After comparing the calculation results of the above relationship with the experimental data, they found that they were in good agreement. After research by Sedov [23], it was found that the self-similarity solution theory was not suitable for medium-intensity explosion waves; in fact, the best method was a series expansion theory. However, Sakurai [21] proposed the series expansion of cylindrical medium-strength explosion waves and found that the series converged very slowly, especially when the pressure of the explosion wave was low. Jones [24] proposed a relational formula suitable for the calculation of medium-intensity spherical shock waves R and PS. Peters et al. [25] also verified the accuracy of Jones’s explosion wave theory through experiments. They found that, compared with Taylor’s explosion wave theory, Jones’s explosion wave theory could still perform well at 10 μs after the formation of hot gas nuclei and shock wave.
It is worth noting that Jones [24] believed that the key problem of solving strong explosion waves was determining the energy actually used to drive the shock wave propagation (the energy coupled into the shock, E0). Only when E0 is determined is it possible to accurately calculate the wave radius R and wave front pressure PS. After certificating the value of R and PS, E0 can be deduced. In the simplified relationship proposed by Gattihe et al., it is believed that all the energy absorbed by the plasma was used to drive the shock wave propagation, and other energy losses were not considered. When Peter et al. used Taylor blast wave theory and Jones blast wave theory to calculate the shock radius, they believed that 90% of the energy absorbed by the plasma was used to drive the shock wave propagation. According to the research of Taylor and Phuoc, it was known that the energy absorbed by the plasma would not be used to drive the shock wave propagation.
Although previous articles had studied the factors that affect the morphology of LIPI flame kernels, most of the research had focused on the effects of pressure, laser energy and equivalence ratio; thus, the influence of oxygen content on the morphology of the flame kernel remains to be explored. In addition, for the process of the LIPI generation shock wave, although the research was early and a variety of models had been proposed, they were not suitable for the experimental conditions of this paper. Therefore, this paper also proposes a modified model based on the Jones model, and compares it with the Taylor model and the Jones model, in order to systematically study the impact of different factors on the characteristics of the shock wave and the flame kernel in the LIPI.

2. Experiment Methods

2.1. Experimental Setup

The experiment is carried out in a constant volume combustion chamber and the experimental system is shown in Figure 1, which is mainly composed of a constant volume combustion chamber, a high-speed shadow system, an ignition system, a gas distribution system, a data acquisition system, and a timing control system.
The constant volume combustion chamber cavity used was a stainless-steel cylinder with an inner diameter of 100 mm, a length of 100 mm, and a volume of 0.86 L (including a part of the pipeline), with two quartz glass windows on each of its axial and radial ends. The diameter of the axial window was 100 mm, which was used to shoot flame images; The diameter of the radial window was 25 mm, which was used for focusing the laser beam and collecting the remaining laser beam energy. A pressure transmitter (Rosemount 3051TA) was installed on the top of the constant volume combustion chamber to measure the static pressure of the air intake. In addition, the piezoelectric pressure sensor (Kistler 7001) was used to measure pressure fluctuations during the combustion process, which was installed at the bottom of the combustion chamber, and its pressure diaphragm was flush with the inner surface of the combustion chamber cavity.
The laser ignition system used a Nd:YAG lamp pumped nanosecond Q-switched laser (Grace NASOR800); three wavelengths of 355 nm/532 nm/1064 nm could be selected, and 532 nm was used in this article (RMS < 1.2%). The laser pulse energy was adjusted by an energy attenuator on the external optical path, which was composed of a half-wave plate and a polarizer. The measurement of laser pulse energy was mainly realized through a beam splitter and two energy meters (Coherent, model: J-50-MB-YAG and J-25-MB-LE). The laser beam could be divided into two beams of a certain proportion by the beam splitter, and by measuring one beam, the energy of the other laser beam could be known. The incident laser energy was measured by a spectroscope and an energy meter, while the remaining energy of the laser was measured by another energy meter.
The high-speed shadow camera system consisted of an LED light source, lens and high-speed camera (Phantom v2511). According to the difference of experimental conditions, the high-speed camera selected two frame rates and corresponding frames. When studying the impact of equivalence ratio on shock wave, the selected high-speed camera had a frame rate of 245,000 and a frame size of 768 × 112 pixels. In other experimental conditions, the high-speed camera had a frame rate of 138,000 and a frame of 768 × 208 pixels.

2.2. Experimental Conditions

According to the research on shock wave theory, the main factors affecting the development of shock wave are laser energy, environmental pressure and gas composition. Therefore, the influence of these parameters on the shock wave characteristics as well as the flame kernel will be studied by changing the initial pressure, laser energy, equivalence ratio, and oxygen content.
The initial temperature of the experiment is 298 K, and the experimental conditions are shown in Table 1. In order to ensure the accuracy and reliability of the experimental data, each working condition is repeated at least twice. The purity of methane, oxygen and nitrogen in the experiment are 99.99%. The laser energy values in Table 1 are obtained by averaging the laser energy multiple times by adjusting the settings of the laser, gradually making the average value of the multiple laser energy close to the laser energy required for the experiment. It should be noted that in order to ensure that the laser energy incident in the combustion chamber is consistent with the experimental conditions, the laser energy here takes into account the loss of quartz glass.

2.3. Data Processing

The radius and velocity of the shock wave could be calculated based on the pictures. Figure 2 shows the flame kernel and shock wave of the CH4/O2 mixture LIPI. The center area was the flame kernel and the outer contour was the shock wave. The cone-shaped structure on the left side of the flame kernel was determined by the characteristics of LIPI, which will be analyzed in detail below. In fact, the shape of the flame kernel directly affected the initial shape of the shock wave, which was the cause of the bulge on the left side of the shock wave, where the bulge direction of the shock wave was consistent with the incident direction of the laser. It was worth noting that the position and size of the plasma formed by laser ignition changed each time, there was a certain deviation between each other. Because the shock wave was generated by the instantaneous heat release of the plasma, its shape and center point position would be affected by the plasma energy. In order to facilitate the calculation of the shock wave velocity, a set of shock wave pictures were obtained in each experimental condition, and the position of the geometric center point needed to be determined separately. First, the position of the geometric center point of the shock wave in the first picture was determined, and the geometric center points in the other pictures were the same as the first picture, so that the geometric center in each experiment could be determined. As shown in Figure 2, the intersection of the vertical lines of the shock wave in the Y+ and X+ directions was determined as the geometric center point.
Due to the limitation of the high-speed camera frame, the number of shock waves obtained in the Y+ and Y− directions is lower. By comparing the radius and velocity of the shock wave of methane/oxygen LIPI in the four directions with the initial pressure of 2.5 atm, it could be found that the radius and velocity of the shock wave in the X−, X+, Y− and Y+ directions were basically the same, so this paper mainly studied the change of the radius and velocity of the shock wave in the X− and X+ direction. The radius could be obtained by measuring the distance from the boundary of the shock wave to the center point, and the shock wave velocity could be calculated according to the differentiation of the shock wave radius with respect to time.

2.4. Models

In order to simplify the calculations, the Taylor explosion wave model and the Jones explosion wave model had been sorted out and simplified, obtaining the following equations:
v a t ( 6.25 E 0 / ( B γ P 0 ) ) 1 / 3 = 0.543 [ ( 1 + 4.61 R ( 6.25 E 0 / ( B γ P 0 ) ) 1 / 3 2.5 ) 0.4 1 ] ,
where, v a —local speed of sound/mm·μs−1; R —shock wave propagation radius/mm; t —shock wave propagation time/μs; γ —medium specific heat ratio; E 0 —energy of driving shock wave propagation/mJ; P 0 —environmental pressure/atm; and B —relevant to γ .
In this paper, after linear fitting of and calculated by Peters [25], the following equations were obtained:
B = 25.525 14.434 γ ,
R = s ( γ ) ( E 0 ρ 0 ) 0.2 t 0.4 ,
where, R —shock wave propagation radius/mm; E 0 —energy of driving shock wave propagation/J; t —shock wave propagation time/μs; ρ 0 —density of surrounding gas/Kg·m−3; γ —medium specific heat ratio; and s ( γ ) —the dimensionless parameter obtained from the equation of motion, which is approximately 1.
Among them, Equation (1) was the Jones simplified model, and Equation (3) was the Taylor simplified model. E0 was calculated as 80% of the energy absorbed by the plasma, and which was 85% of the difference value between the incident laser energy and the remaining energy, so E0 was equal to 68% of the energy difference [8]. The basic parameter values of the model with different experimental conditions were shown in Table 2.
In this paper, by calculating the Jones explosion wave model and the Taylor explosion wave model, the curve of shock wave radius and velocity with time with different experimental conditions could be obtained; the results were compared with the experimental data. The calculation results of the Taylor model had great deviations from the experimental data, while the Jones model fits the experimental data relatively well, but there were still deviations. Therefore, this paper proposed a modified model Equation (4) based on the Jones model.
This model was derived by fitting the experimental data of CH4/O2 mixture with pressure of 2.5 atm. In this paper, the Jones model was used as the curve fitting function, and the two numbers 0.543 and 0.4 in Equation (1) were determined as variables. Then, by continuously adjusting these two variables, the curve fitting function with the best fitting effect was obtained:
v a t ( 6.25 E 0 / ( B γ P 0 ) ) 1 / 3 = 0.323 [ ( 1 + 4.61 R ( 6.25 E 0 / ( B γ P 0 ) ) 1 / 3 2.5 ) 0.45 1 ]

3. Results and Discussions

The LIPI shock wave and flame kernel of the CH4/O2 mixture at initial temperature 298 K is shown in Figure 3, with initial pressure of 2.5 atm, and an equivalence ratio of 1. The first picture is taken 7.4 μs after the camera is triggered; it can be seen that there is still a strong plasma signal in the central area. At the same time, a more obvious flame kernel and shock wave boundary can be observed, but the shape of the flame kernel is different from the usual one ring-shaped or third-lobe structure, and a cone-shaped structure appears. After that, flame kernel and shock wave gradually develop and spread outward. It can be noticed that after 44.4 μs, the brightness of the central area of the flame kernel gradually increases, while the boundary of the shock wave is gradually blurred and difficult to identify. In addition, in the initial stage of the formation of the flame kernel, the axial development speed is faster than the radial direction, but then the radial development speed is faster, which leads to the gradual disappearance of the cone-shaped structure. The following is the detailed analysis of the shock wave development and flame kernel morphology under different experimental conditions.

3.1. The Effect of the Pressure

Figure 4 shows the changes of the radius and velocity of the shock wave of CH4/O2 laser ignition with the initial pressures when the equivalence ratio is 1, and compares the experimental data with three models. Figure 4a,b compare the difference in radius and velocity in four directions with the initial pressure of 2.5 atm. The radius of the shock wave in the four directions gradually increases with the increase of time; the change of the velocity also gradually increases with time, and then stabilizes. In this paper, the radius of the shock wave is calculated by measuring the number of pixels in the pictures according to the scale. It is worth noting that the deviation of the radius in different directions in the figure is mainly caused by the deviation of the pixel point measurement, which directly leads to the large fluctuation of the velocity. In addition, the shock wave has thickness, and the boundary of the shock is not clear, which leads to a certain error in the number of pixels; the deviation of 1 pixel will cause an error of about 18 m/s.
It can be seen from the figure that there is a big deviation between the calculation results of the Taylor model and the experimental data. This is mainly because the model is only applicable to the strong explosion wave model and can predict the development and changes within a few microseconds after the shock wave is formed. However, after dozens of microseconds, the shock wave will gradually propagate to the sound speed due to energy loss and environmental pressure, so the Taylor model is not suitable for predicting the development of the shock wave in the later stage. Compared with the Taylor model, the Jones model fits the experimental data better. This is because the Jones model is used to predict the medium-intensity explosion wave, and the propagation of the shock wave in the later stage can be considered as the medium-intensity explosion wave propagation. Although the predicted results of the Jones model are relatively good, there are still deviations from the experimental data. For this reason, this paper has performed curve function fitting on the experimental data points in Figure 4a,b based on the Jones model, and obtained the modified model (Equation (4)) based on the Jones model. As can be seen from the figure, compared to the Jones and Taylor models, the modified model has the best fitting results.
The changes in the radius and velocity of the shock wave in the X– and X+ directions with different pressures are shown in Figure 4c–f. With the propagation of shock wave, its radius gradually increases, and its velocity gradually decreases and tends to be stable afterwards. It can be seen that, with the increase of pressure, there is no obvious change in the radius and velocity of the shock wave. The modified model fits well with the experimental results and can better show the development trend of the shock wave. Through the calculation results of the three models in the figure, it can be found that the initial pressure has a small effect on the velocity, but as the initial pressure increases, the shock wave velocity will gradually decrease, mainly because the increase in environmental pressure will hinder the propagation of shock wave. According to the calculation results of the modified model: the initial pressures are 0.4 atm and 2.5 atm, and the corresponding shock wave velocities are 428.96 m/s and 405.78 m/s at 20 μs respectively; the difference of the shock wave velocities corresponding to the two initial pressures are almost the same as the experimental data error, which may explain why the influence in the shock wave velocity of the initial pressure cannot be seen through the experimental data.
For the different initial pressures, the evolution of CH4/O2 flame kernel morphology, are shown in Figure 5. When the initial pressure is 0.4 atm and t = 14.8 μs, the initial flame kernel is close to a circle, but a cone-shaped structure can be seen. At t = 66.6 μs, a ring structure gradually forms on the left side of the flame kernel, and a vortex appeared in the middle. However, the intermediate vortex will not develop into a more obvious ring structure, and the size of the vortex gradually decreases. In comparison, when the initial pressure is 0.6 atm and t = 14.8 μs, a more obvious cone-shaped structure can be seen on the left side of the initial flame kernel. At t = 66.6 μs, a small vortex appears in the middle of the left side of the flame kernel, and will gradually disappear. As the initial pressure continues to increase, the cone-shaped structure becomes more obvious, and the axial length of the flame kernel increases; while the vortex structure gradually disappears, the flame kernel gets closer and closer to an ellipse.
From the analysis of the hot gas kernel morphology, there are two possible reasons for the formation of the ring and third-lobe structure: one is the interaction between the sparse wave and the hot gas kernel, and the other is the asymmetric deposition of laser energy. In this paper, the possible reason is that only the vortex structure appears at the initial pressure of 0.4 atm, the ring and third petal structures are not formed, for the chemical reaction rate is accelerated in the pure oxygen environment. The hot gas nucleus quickly forms a stable flame kernel, which forms the ring structure, and the vortex motion of the third-lobe structure is weakened. In addition, when the initial pressure is low, the absorbed energy of the focused plasma is reduced, and the influence of asymmetric deposition of laser energy is also weakened. As the initial pressure increases, the flame kernel gradually forms a conical structure on the side far from the incident direction of laser (hereinafter referred to as the distal end); the side close to the incident direction of laser (hereinafter referred to as the proximal end) forms a spherical structure. This may be due to the combined effects of the asymmetric deposition of laser energy and the acceleration of the chemical reaction rate, but the asymmetric deposition of laser energy may lead to more deposition energy at the proximal end, which makes the area where the surrounding gas heated larger and the chemical reaction rate faster. Therefore, the flame spreads rapidly in all directions in this area and the propagation rate is almost the same, thereby forming a spherical structure, inhibiting the formation of the third-lobe structure. In addition, the accelerated chemical reaction rate not only weakens the vortex motion, but also promotes the formation of a cone-shaped structure.

3.2. The Effect of the Energy

The changes of the radius and velocity of the shock wave with different laser energy in the X− and X+ direction of the CH4/O2 mixture at initial pressure of 1 atm, as shown in Figure 6a,b. The shock wave velocity gradually decreases, then tends to stabilize with the propagation of flame. It can be seen that there are certain fluctuations in the shock wave velocity in both directions, but there is no obvious relationship between the shock wave velocity and the laser energy. Furthermore, from Figure 6c,d, it can be found that compared with the Jones and Taylor models, the modified model can still better predict the development of shock wave velocity. The model results show that as the laser energy increases, the shock wave velocity gradually increases, because the increase in laser energy will increase the energy absorbed by the plasma, which in turn increases the energy used to drive the shock wave. When t = 20 μs and the laser energy is 21.42 mJ and 217.67 mJ, the corresponding shock wave velocities are 399.56 m/s and 435.87 m/s respectively, and the velocity difference value is 36.31 m/s.
Through the above analysis, it can be known that the laser energy has no obvious effect on the shock wave velocity, and a few little differences in the velocity may be caused by the experimental errors. Therefore, further experimental measurements and research are needed in the future.
The development of CH4/O2 flame kernel morphology using different laser energy is shown in Figure 7. When the laser energy is 21.42 mJ and t = 66.6 μs, a less obvious vortex structure appears on the left side of the flame kernel, while the vortex structure gradually disappears with the development of flame. With the increase of laser energy, a more obvious cone-shaped structure appears on the left side of the initial flame kernel, but the size of the flame kernel does not change significantly. This may be due to the fact that increasing the laser energy can speed up the chemical reaction rate, but when the energy is increased to a certain value, the chemical reaction will no longer be affected by the laser energy.

3.3. The Effect of Equivalence Ratio

When the initial pressure is 1 atm, the influence of CH4/O2/N2 mixture with different equivalence ratios on shock wave velocity is shown in Figure 8. From the experimental data in the figure, it can be seen that the change in equivalence ratio will hardly affect the development of the shock wave. Moreover, the modified model can better predict the change trend of the shock wave velocity, especially in the X+ direction. The simulation results of the modified model reveal that at t = 20 μs, when the equivalence ratios are 0.7 and 1.3, the corresponding shock wave velocities are 424.29 m/s and 423.63 m/s, respectively. This shows that when the laser energy and initial pressure are constant, the change of the component concentration has little effect on the shock wave velocity.
The flame kernel morphology of the CH4/O2/N2 mixture with different equivalence ratios is shown in Figure 9. Different from the shock wave, it can be seen that when the equivalence ratio is 0.7 and 1.3, the flame kernel has obvious ring-shaped and third-lobed structures, while when the equivalence ratio is 1.3, the characteristics of the ring-shaped and third-lobe structures are more obvious. This is mainly due to the fact that the flame development is more stable during rich combustion and the flame structure is easier to maintain.

3.4. The Effect of Oxygen Content

Figure 10 shows the relationship between the shock wave velocity and the oxygen content of the CH4/O2/N2 mixture when the equivalence ratio is 1, and the initial pressure is 1 atm. It can be seen from the figure that the shock wave velocity does not change significantly with the increase in oxygen content. The modified model can still better predict the trend of shock wave velocity. According to the results of the modified model, at t = 20 μs and an oxygen content of 21% and 100%, the corresponding shock wave velocities are 422.89 m/s and 417.44 m/s, respectively. However, it can be seen from Table 2 that γ decreases with the increase of oxygen content, so the change of component concentration has a certain effect on the shock wave velocity, and it decreases with the decrease of the ratio of the specific heat of the mixture. However, the change of velocity is not obvious, because the γ has only limited effect.
The change of the flame kernel morphology of the CH4/O2/N2 mixture with different oxygen content is shown in Figure 11. It can be seen that, as the oxygen content gradually increases, the plasma and flame kernel morphologies are obviously different. When t = 7.4 μs, as the oxygen content increases, the brightness of the plasma region gradually decreases, while the size of the flame kernel gradually increases. This may be due to the fact that the increase of oxygen content increases the rate of chemical reaction, then accelerates the absorption and dissipation of energy in the plasma region, so the brightness of this region decays rapidly. Meanwhile, a difference can be seen in the structure of the flame kernel with the increase of oxygen content, and the annular and third-lobed structures disappear completely when the oxygen content increases to 60%, but the size of the flame kernel also increases rapidly with the increase in oxygen content. Through comparative analysis of the development and changes of the flame kernel, the disappearance of the ring-shaped and third-lobe structures is mainly caused by the acceleration of the chemical reaction rate, so that the flame kernel rapidly develops and expands outwards, which will inhibit the airflow of the vortex, resulting in weakening of the vortex structure and the failure of the ring and third-lobe structures to form.

4. Conclusions

In this paper, the development characteristics of shock waves and flame kernels with different initial pressure, laser energy, equivalence ratio and oxygen content are studied in the constant volume combustion chamber. Furthermore, the Taylor model and the Jones model were used to simulate the development of shock wave with different conditions, and a new modified model is proposed. However, this experiment only studied the development of the flame kernel and the shock wave, and did not measure the products generated in the subsequent combustion process, such as water and carbon dioxide. This has important practical value for the application of laser ignition with some engineering working conditions, such as high pressure and oxygen enrichment. Future work can include measurement of the reaction products. The main conclusions are as follows:
  • The initial pressure and laser energy will affect the shock wave velocity from the theoretical analysis; however, the results are not obvious because of the experimental error and measurement devices. The shock wave velocity decreases with the increase of the initial pressure, while increases with the increase of the laser energy. Changes in the equivalence ratio and oxygen content have little effect on the shock wave.
  • The calculation results of the Taylor model have great deviations from the experimental data, while the Jones model fits the experimental data relatively well, but there are still deviations. This paper proposes a modified model based on the Jones model, and the modified model proposed can predict the development trend of shock waves well.
  • The chemical reaction rate of the CH4/O2 mixture is too fast, causing the flame kernel to expand rapidly, which will inhibit the motion of the vortex airflow, which makes it difficult for the flame kernel to form the ring and third-lobe structure. Reducing the initial pressure and laser energy is beneficial to the formation of the ring structure, but the third-lobe structure cannot be formed.
  • It is easier to maintain the shape of the flame kernel in a richer equivalence ratio than lower, and it is not easily affected by the movement of the hot air flow, and the boundary between the ring and the third-lobe structure is clearer. Increasing the oxygen content of the mixed gas will inhibit the formation of the ring and the third-lobe structure, which is also affected by the excessively fast chemical reaction rate. The increase in oxygen content will reduce the duration of the plasma, so that the energy in the plasma region is quickly absorbed and dissipated by the combustible mixture.

Author Contributions

Conceptualization, J.Z. and Q.G.; Methodology, J.Z. and Q.G.; Formal Analysis, J.Z.; Q.G. and G.Y.; Investigation, J.Z. and Q.G.; Original Draft Preparation, J.Z. and Q.G.; Data Curation, J.Z. and Q.G.; Writing, J.Z.; Validation, E.H.; Resources, E.H.; Supervision, E.H. and Z.H.; Project Administration, E.H. and Z.H.; Funding Acquisition, E.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (52176131).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, J.Q.; Li, Q.L.; Shen, C.B. Ignition characteristics and combustion performances of a LO2/GCH4 small thrust rocket engine. J. Cent. South Univ. 2018, 25, 646–652. [Google Scholar] [CrossRef]
  2. Figiel, A.; Klačková, I. Safety requirements for mining systems controlled in automatic mode. Acta Montan. Slovaca 2020, 25, 417–436. [Google Scholar]
  3. Kuric, I.; Klačková, I.; Nikitin, Y.R.; Zajačko, I.; Císar, M.; Tucki, K. Analysis of diagnostic methods and energy of production systems drives. Processes 2021, 9, 843. [Google Scholar] [CrossRef]
  4. Lu, L.; Han, X.; Li, J.; Hua, J.; Ouyang, M. A review on the key issues for lithium-ion battery management in electric vehicles. J. Power Sources 2013, 226, 272–288. [Google Scholar] [CrossRef]
  5. Phuoc, T. Laser-induced spark ignition fundamental and applications. Opt. Laser Eng. 2006, 44, 351–397. [Google Scholar] [CrossRef]
  6. Bradley, D.; Sheppard, C.; Suardjaja, I.M.; Woolley, R. Fundamentals of high-energy spark ignition with lasers. Combust. Flame 2004, 138, 55–77. [Google Scholar] [CrossRef]
  7. Bao, X.; Sahu, A.; Jiang, Y.; Badawy, T.; Xu, H. Flame kernel evolution and shock wave propagation with laser ignition in ethanol-air mixtures. Appl. Energ. 2019, 233–234, 86–98. [Google Scholar] [CrossRef]
  8. Bradley, D.; Lung, K.K. Spark ignition and the early stages of turbulent flame propagation. Combust. Flame 1987, 69, 71–93. [Google Scholar] [CrossRef]
  9. Taylor, G.I. The Formation of a Blast Wave by a Very Intense Explosion. Proc. R. Soc. Lond. A Math. Phys. Sci. 1950, 201, 159–174. [Google Scholar]
  10. Akram, M. Two-dimensional model for spark discharge simulation in air. AIAA J. 2012, 34, 1835–1842. [Google Scholar] [CrossRef]
  11. Kono, M.; Niu, K.; Tsukamoto, T.; Ujiie, Y. Mechanism of flame kernel formation produced by short duration sparks. Symp. Combust. 1989, 22, 1643–1649. [Google Scholar] [CrossRef]
  12. Kravchik, T.; Sher, E. Numerical modeling of spark ignition and flame initiation in a quiescent methane-air mixture. Combust. Flame 1994, 99, 635–643. [Google Scholar] [CrossRef]
  13. Spiglanin, T.A.; Mcilroy, A.; Fournier, E.W.; Cohen, R.B.; Syage, J.A. Time-resolved imaging of flame kernels: Laser spark ignition of H2/O2/Ar mixtures. Combust. Flame 1995, 102, 310–328. [Google Scholar] [CrossRef]
  14. Raizer, Y.P. Laser-induced discharge phenomena. Nucl. Fusion 1977, 18, 11. [Google Scholar]
  15. Wohlhüter, M.; Zhukov, V.P.; Sender, J.; Schlechtriem, S. Analysis of the laser ignition of methane/oxygen mixtures in a sub-scale rocket combustion chamber. CEAS Space J. 2017, 9, 211–225. [Google Scholar] [CrossRef]
  16. Dumitrache, C.; Baumgardner, M.; Boissiere, A.; Maria, A.; Roucis, J.; Marchese, A.J.; Yalin, A. A study of laser induced ignition of methane–air mixtures inside a Rapid Compression Machine. Proc. Combust. Inst. 2016, 36, 3431–3439. [Google Scholar] [CrossRef] [Green Version]
  17. Meier, W.; Gebel, G.; Mosbach, T. Analysis of Laser Ignition of Kerosene Spray by Optical Diagnostics. In Laser Ignition Conference; Optical Society of America: Argonne, IL, USA, 2015. [Google Scholar]
  18. Leela, C.; Bagchi, S.; Kumar, V.R.; Tewari, S.P.; Kiran, P.P. Dynamics of laser induced micro-shock waves and hot core plasma in quiescent air. Laser Part. Beams 2013, 31, 263–272. [Google Scholar] [CrossRef] [Green Version]
  19. Shen, Y.R. The principles of nonlinear optics. IOP Sci. 1984, 36, 178. [Google Scholar]
  20. Lin, S. Cylindrical Shock Waves Produced by Instantaneous Energy Release. J. Appl. Phys. 1954, 25, 54–57. [Google Scholar] [CrossRef]
  21. Sakurai, A. On the Propagation and Structure of a Blast Wave, II. J. Phys. Soc. Jpn. 1953, 8, 662. [Google Scholar] [CrossRef]
  22. Gatti, M.; Palleschi, V.; Salvetti, A.; Singh, D.P.; Vaselli, M. Spherical shock waves in laser produced plasmas in gas. Opt. Commun. 1988, 69, 141–146. [Google Scholar] [CrossRef]
  23. Sedov, L.I.; Friedman, M.; Holt, M.; Cole, J.D. Similarity and Dimensional Methods in Mechanics. J. Appl. Mech. 1961, 28, 159–160. [Google Scholar] [CrossRef] [Green Version]
  24. Jones, D.L. Intermediate Strength Blast Wave. Phys. Fluids 2003, 11, 1664–1667. [Google Scholar] [CrossRef]
  25. Peters, N.D. Investigation of the Multi-Physics of Laser-Induced Ignition of Transportation Fuels. Ph.D. Thesis, Syracuse University, Syracuse, NY, USA, 2017. [Google Scholar]
Figure 1. Diagram of the laser ignition experiment.
Figure 1. Diagram of the laser ignition experiment.
Energies 14 07976 g001
Figure 2. Schematic diagram of laser-induced plasma ignition (LIPI) shock wave radius measurement for CH4/O2 mixture.
Figure 2. Schematic diagram of laser-induced plasma ignition (LIPI) shock wave radius measurement for CH4/O2 mixture.
Energies 14 07976 g002
Figure 3. The image of the shock wave developed by CH4/O2 laser ignition over time (T = 298 K, P 0 = 2.5 atm, Φ = 1).
Figure 3. The image of the shock wave developed by CH4/O2 laser ignition over time (T = 298 K, P 0 = 2.5 atm, Φ = 1).
Energies 14 07976 g003
Figure 4. Changes of CH4/O2 shock wave radius and velocity with different initial pressures. (a) Radius in different directions; (b) velocity in different directions; (c) radius in X+ directions; (d) velocity in X+ directions; (e) radius in X– directions; (f) velocity in X– directions.
Figure 4. Changes of CH4/O2 shock wave radius and velocity with different initial pressures. (a) Radius in different directions; (b) velocity in different directions; (c) radius in X+ directions; (d) velocity in X+ directions; (e) radius in X– directions; (f) velocity in X– directions.
Energies 14 07976 g004
Figure 5. Comparison of CH4/O2 flame kernel morphology with different initial pressures.
Figure 5. Comparison of CH4/O2 flame kernel morphology with different initial pressures.
Energies 14 07976 g005aEnergies 14 07976 g005b
Figure 6. Changes of CH4/O2 shock wave velocity with different laser energy. (a) Velocity in X− directions; (b) velocity in X+ directions; (c) comparison of experimental data and model results in X− directions; (d) comparison of experimental data and model results in X+ directions.
Figure 6. Changes of CH4/O2 shock wave velocity with different laser energy. (a) Velocity in X− directions; (b) velocity in X+ directions; (c) comparison of experimental data and model results in X− directions; (d) comparison of experimental data and model results in X+ directions.
Energies 14 07976 g006
Figure 7. Comparison of CH4/O2 flame kernel morphology with different laser energy.
Figure 7. Comparison of CH4/O2 flame kernel morphology with different laser energy.
Energies 14 07976 g007
Figure 8. Changes of CH4/O2 shock wave velocity with different equivalence ratio. (a) Velocity in X− directions; (b) velocity in X+ directions.
Figure 8. Changes of CH4/O2 shock wave velocity with different equivalence ratio. (a) Velocity in X− directions; (b) velocity in X+ directions.
Energies 14 07976 g008
Figure 9. Comparison of CH4/O2 flame kernel morphology with different equivalence ratios.
Figure 9. Comparison of CH4/O2 flame kernel morphology with different equivalence ratios.
Energies 14 07976 g009aEnergies 14 07976 g009b
Figure 10. Changes of CH4/O2 shock wave velocity with different oxygen content. (a) Velocity in X− directions; (b) velocity in X+ directions.
Figure 10. Changes of CH4/O2 shock wave velocity with different oxygen content. (a) Velocity in X− directions; (b) velocity in X+ directions.
Energies 14 07976 g010
Figure 11. Comparison of CH4/O2 flame kernel morphology with different oxygen content.
Figure 11. Comparison of CH4/O2 flame kernel morphology with different oxygen content.
Energies 14 07976 g011
Table 1. Experimental conditions of CH4/O2/N2 mixed gas laser ignition.
Table 1. Experimental conditions of CH4/O2/N2 mixed gas laser ignition.
MixturePressure
P 0 / atm
Oxygen
Content
Equivalence
Ratio / Φ
Laser Energy
E 0 / mJ
CH4/O20.4–2.5100%1150
CH4/O2/N2121–100%1150
CH4/O2/N2121%0.7,1.3150
CH4/O21100%150–300
Table 2. Parameters in the Jones model and the Taylor model.
Table 2. Parameters in the Jones model and the Taylor model.
MixtureOxygen Content Φ E 0 / mJ P 0 / atm ρ 0 / Kg · m 3 v a / m · s 1 γ B
CH4/O2100%1.058.210.40.4364341.121.3625.86
59.570.60.6546
67.430.80.8728
71.331.01.0910
76.941.51.6366
82.762.02.1821
82.212.52.7276
CH4/O2/N221%1.073.091.01.1301344.491.3895.48
40%73.541.1184343.491.3815.59
60%72.751.1079342.501.3735.71
80%72.221.0989341.751.3675.79
100%72.331.0910341.121.3625.87
CH4/O2/N221%0.780.441.01.144344.861.3925.43
1.382.201.117344.111.3865.52
CH4/O2100%1.021.421.01.0910341.121.3625.86
55.42
97.17
159.44
189.31
217.67
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, J.; Hu, E.; Gao, Q.; Yin, G.; Huang, Z. Shock Wave Propagation and Flame Kernel Morphology in Laser-Induced Plasma Ignition of CH4/O2/N2 Mixture. Energies 2021, 14, 7976. https://doi.org/10.3390/en14237976

AMA Style

Zhang J, Hu E, Gao Q, Yin G, Huang Z. Shock Wave Propagation and Flame Kernel Morphology in Laser-Induced Plasma Ignition of CH4/O2/N2 Mixture. Energies. 2021; 14(23):7976. https://doi.org/10.3390/en14237976

Chicago/Turabian Style

Zhang, Junjie, Erjiang Hu, Qunfei Gao, Geyuan Yin, and Zuohua Huang. 2021. "Shock Wave Propagation and Flame Kernel Morphology in Laser-Induced Plasma Ignition of CH4/O2/N2 Mixture" Energies 14, no. 23: 7976. https://doi.org/10.3390/en14237976

APA Style

Zhang, J., Hu, E., Gao, Q., Yin, G., & Huang, Z. (2021). Shock Wave Propagation and Flame Kernel Morphology in Laser-Induced Plasma Ignition of CH4/O2/N2 Mixture. Energies, 14(23), 7976. https://doi.org/10.3390/en14237976

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