Calculation method of spherically expanding flame propagation radius to consider ignition electrode effects | Scientific Reports
Scientific Reports volume 14, Article number: 8435 (2024) Cite this article
305 Accesses
Metrics details
Ignition electrodes have an immense impact on the accurate measurement of the flame propagation spherical radius. In this study, a flame-radius calculation method is designed. The method is able to eliminate effects due to the ignition electrodes. The adaptability and optimization effects of the proposed method are analyzed. The results show that the ratio of the angle is affected by the ignition electrodes under the Han II method. There are three obvious divisions include a high-value area, a sharp-variation area, and a mild-variation area. The ratio of the angle affected by the ignition electrodes is only applicable to the mild-variation region when the flame presents respective convex and concave distributions. For these distributions, the increment rate of the mean radius is 0.4–0.85% and 0.42–3.19%. The reduced rate of the standard deviation of the radius extraction value is 11.91–22.1% and 5.13–17.99%, and the reduced rate of the radius extraction value range is 20.32–39.51% and 0.32–8.09%.
The laminar flame is the most basic combustion phenomenon that comprehensively embodies the characteristics of the gas flow property, chemical reaction, and molecular transport. It is an important means of developing and verifying the kinetic mechanism of a fuel chemical reaction1,2. It is also an important basis for the study of turbulent flames3. Moreover, it provides basic data for engine combustion simulation and optimization research4. To fully elucidate the laminar flame, an effective method is to obtain the flame propagation radius using a constant volume incendiary bomb and high-speed camera. Accordingly, Schlieren images are obtained, processed, and calculated. Therefore, the calculation method for the high-precision flame propagation radius is the primary input condition that is used to understand the characteristics of laminar flame combustion. Meanwhile, it is the key to determining the accuracy of laminar flame propagation velocity.
Many studies have been conducted on flame radius calculation methods. Some examples of common methods of this type are the equivalent area method5,6, random circle detection algorithm7, circle fitting8,9,10, and center measurement method11. Gu12, Milton13, and Zhang14 determines the flame radius by calculating the equivalent area. Chen et al.6 proposes a random circle detection algorithm to measure the flame radius by determining the true circle according to the distance threshold requirements. Broustail8, Han9, Bouvet10, DuttaRoy15, Gong16, and Tahtouh17 calculates the flame propagation radius by circle fitting of the flame front-edge points. Zuo et al.11 extracted the flame propagation radius through center measurement method. Wu et al.18uses operators to detect flame edges and improve the traditional spherical flame arc fitting method to calculate the flame radius. Li et al.5 applies the inter-frame difference method to process the images. They use the Random Circle Detection (RCD) algorithm, equivalent area method to identify flame rings, calculate the flame radius and center position.
By conducting in-depth research on the flame radius calculation method, some researchers determine that the flame radius near the ignition electrodes appeared as flame bumps and non-standard circle phenomena. Liang et al.19 use the Canny edge detection operator to detect the front edge of the flame in Schlieren images and then compare the difference between the arc fitting method and the flame area method in calculating the flame radius. They selected the area method for the flame radius calculation. They find out that the shape of the spherical flame is affected by electrode factors in the propagation process, which leads to a certain protrusion near the ignition electrodes. Wang et al.20 applies the Schlieren method to scale the recorded flame images and measure the flame radius. It is determined that the cooling effects of the ignition electrodes caused by the flame to propagate slowly along the ignition electrode direction. Han et al.21 confirmed that the existence of ignition electrodes is one of the main reasons for the distortion of the flame propagation radius. The existence of the ignition electrode leadS to an increase in the dispersion of the radius extraction value. The flame profile deformation caused by the ignition electrode is shown in Fig. 1.
Deformation of flame profile caused by the ignition electrode.
Therefore, some scholars have gradually begun to study the phenomenon that the flame propagation radius near the ignition electrode produces a non-standard circle. Mevel and Lafosse22 writes a software program that takes the radius of the equivalent circle as the flame propagation radius. Although the nonstandard circle is considered in this method, the effects of ignition electrodes on accurate flame radius measurement are not separately considered. Liang19 and Wang20 only proves the effects of the ignition electrodes on the flame radius, but they do not address the effects on the accurate measurement of the flame radius. Varea et al.23 proposes a method of flame radius post-processing based on the effects of ignition electrodes on the flame radius. They apply this method to the extraction of the laminar combustion parameters of ethanol/air and methanol/ethanol/isooctane/air mixture. Although they propose a corresponding calculation method for the effects of ignition electrodes on the flame propagation radius, they only solve the problem of the non-standard circle of the flame surface engendered by the electrodes. The influence range of the ignition electrode is not further studied. Han et al.21 establishes an optimization method for the influence of ignition electrodes. However, the optimization method can only solve the problem of electrode interference with the protrusion of the flame propagation radius near the ignition electrodes in certain conditions. It is not suitable for the non-standard flame radius near most ignition electrodes.
In summary, ignition electrodes that cause flame bumps and non-standard circle phenomena significantly affects the accurate measurement of the spherical radius of the flame propagation. Nevertheless, many researchers have only focused on the measurement of the laminar flame propagation radius. To date, few calculation methods have been devised that eliminate the effects of ignition electrodes on the laminar flame propagation radius, which would enable the spherical flame propagation radius to be accurately measured. Therefore, a flame radius calculation method is herein proposed to remove the effects of ignition electrodes based on the above flame radius calculation method. Figure 2 shows the flame profile without the influence of ignition electrode. Compared with the equivalent circle area method, the circle fitting method and the center measurement method, the maximum difference of the calculated flame radius is 4.54%, 5.19% and 9.17%, respectively. This method not only provides a means of supporting the accurate measurement of the spherical flame propagation radius, but also sets a foundation for scholars to conduct high-precision laminar combustion studies and contribute to the advancement of basic combustion research.
Flame profile without ignition electrode effect.
The experiments of this study are conducted in a constant volume combustion vessel (CVCV), which is a sphere with an inner diameter of 260 mm. The CVCV has two 100 mm quartz windows on each side to provide visual access for the Schlieren photography system. Electrodes extend into the center of the CVCV to realize ignition in Fig. 3. The initial temperature in the CVCV is controlled within 600 K by using a heating control system. A high-speed camera is used to record the flame propagation process at an image resolution of 512 × 512. The parameters of the sensors are used in this test are shown in Table 1. Moreover, all experimental protocols are approved by the Xihua University.
Schematic of the CVCV experimental system.
Extract the residual exhaust gas from the CVCV after the previous test and bring it to a vacuum state. Then, the gas of each component is slowly entered into the CVCV by using the corresponding partial-pressure method, and the CVCV is simultaneously heated until the initial temperature is reached. Only the flame images with 6 mm < R < 25 mm are used to avoid the effects of ignition energy and combustion pressure24. Moreover, flame images with obvious cellular structures are omitted25.
The test conditions are shown in Table 2. The diluent gas fraction \(({\Phi }_{r})\) represents the volume fraction of diluted gas (CO2) in all the gas mixtures. The diluting gas fraction is defined as the volume fraction of the diluting gas in the mixed gas, which is used to describe the relative content of the diluting gas in the mixed gas. The equivalence ratio is defined as the ratio of the theoretically required amount of air for complete combustion to the actual amount of air supplied during fuel combustion.
The flame image processing is shown in Fig. 4. A program using MATLAB code that is independently developed by the present research group is employed to process the Schlieren images through six steps: removing the background, transforming it into a gray image, increasing the gray level, removing the island, removing isolated highlights, and detecting image edges.
Flame image processing by MATLAB.
The two-dimensional coordinate system and measuring center Q of the Schlieren images are obtained based on the Han I method. The Schlieren images are divided into four regions. Then, the fitted edge contour points and extracted value of the fitting radius (\({r}_{\theta }\)) are obtained after removing the ignition electrodes and noise points of the images. The characteristic radius (\({r}_{k}\)) is the extracted value of the fitting radius in the direction toward the ignition electrodes in each of the four areas. The affected angle caused by the ignition electrodes (\({\theta }_{b}\)) is determined based on the values of the evaluation parameters of the affected angle caused by the ignition electrodes (\({\Delta r}_{max}\)). The above method for calculating the flame radius with greater accuracy by removing the ignition electrode effects supplements and improves the Han I method. In this paper, the optimized method is named the Han II method. The flame edge contour points are obtained using the output results of the method, measurement center Q', and the flame radius after removal of the ignition electrode effects (\({R}_{\theta }^{\prime}\)) are obtained by removing the angle of the ignition electrode effects, as shown in Fig. 5.
Frame diagram of the Han II method.
In Fig. 6, a two-dimensional coordinate system is established based on measurement center point Q, which is determined by the Han I method9. The Schlieren images are divided into four regions using the horizontal and vertical coordinate system. The corresponding angle regions are defined as region I (0°–90°), region II (90°–180°), region III (180°–270°), and region IV (270°–360°), respectively. Passing through measurement center point Q in the direction of the ignition electrode connection is taken as the vertical axis, and the straight line perpendicular to the vertical axis is the horizontal axis. The intersection point between one side of the horizontal axis and the edge of the flame front surface is M1, which is the starting point. Then, QM1 is a fixed edge of the measuring angle in this coordinate system. In determining any flame front surface edge contour point P in the clockwise direction, angle M1QP is measurement angle \(\theta \) (\(\theta \in [0^\circ , 360^\circ ]\)) in this coordinate system.
Dividing the regions.
In the process of brightening gray-scale images through processing, numerous contour points and nearby noise points will be generated near the ignition electrodes. This will greatly affect the recognition of contour points by the MATLAB code. It is thus necessary to remove the contour points and nearby noise points of the ignition electrodes. The intersection points between the outline of ignition electrodes and the flame edge are A, D and B, C in Fig. 7. The ignition electrode pixel edge is defined by the two straight lines that pass-through points A, D and B, C, respectively. The contour and noise points within the width of the ignition electrode pixel edge must be removed.
Removing the contour and noise points within the width of the ignition electrode pixel edge.
The corresponding random errors are generated since the original edge contour points are discrete points with an uneven distribution and randomness in Fig. 8. This makes it difficult to identify the ignition electrode effect boundary. To effectively eliminate random errors, the original contour points must be fitted because the correlation coefficients of the fitting above the ninth order are very close. Since the polynomial fitting coefficients are almost equal when the polynomial order is greater than 9, the least square method is used to fit the original contour points with nine-order polynomial.
Discrete points of edge contour.
After obtaining the fitting edge contour points, a two-dimensional coordinate system is established using Schlieren images with measurement center point Q as the center of the circle, as determined by the Han I method. As shown in Fig. 9, the coordinates of point Q and point P in the coordinate system are defined as (\({x}_{0}\),\({y}_{0}\)) and (\({x}_{i}\),\({y}_{i}\)), respectively. In region I, it can be observed that the relation between line PQ and the corresponding angle θ is
Calculating the extracted value of fitting radius (\({r}_{\theta }\)).
In region II and III, the relation between line segment PQ and corresponding angle θ is
In region IV, the relation between line segment PQ and corresponding angle θ is
The pixel difference between two points of PQ can be obtained by the circle radius function formula:
The extracted value of the fitting flame radius (\({r}_{\theta }\)) can be calculated from calibration ratio K21, \(\theta \in [0^\circ , 84^\circ ] \cup [96^\circ , 264^\circ ] \cup [276^\circ , 360^\circ ]\). The equation is as follows:
Here, calibration ratio \(K=S/{S}_{0}\), where S is the actual size of the scale, and \({S}_{0}\) is the pixel difference corresponding to the actual size of the scale in the Schlieren images.
The extracted value of the fitting radius (\({r}_{\theta }\)), \(\theta \in [0^\circ , 84^\circ ]\cup [96^\circ , 264^\circ ]\cup [276^\circ , 360^\circ ]\) is obtained. The extracted radius value \({r}_{84^\circ }\) in region I is defined as the characteristic radius \({{\text{r}}}_{{\text{I}}}\), the extracted radius value \({r}_{96^\circ }\) in region II is defined as the characteristic radius \({{\text{r}}}_{{\text{II}}}\), the extracted radius value \({r}_{264^\circ }\) in region III is defined as the characteristic radius \({{\text{r}}}_{{\text{III}}}\), the extracted radius value \({r}_{276^\circ }\) in region IV is defined as the characteristic radius \({{\text{r}}}_{{\text{IV}}}\) in Fig. 10.
Definition of characteristic radius.
The effects of the ignition electrodes on the flame edge contour are manifested as the mutation region, which includes the convex and concave shapes of the flame. As shown in Fig. 11, the maximum difference between the extracted value of the fitting radius (\({r}_{\theta }\)) and the characteristic radius (\({r}_{k}\)) can properly reflect the effects of the mutation region caused by the ignition electrodes on the flame propagation radius. Thus, the evaluation parameters of the angle affected by the ignition electrodes (\(\Delta {r}_{max}\)) for the four regions are defined as
Determining the parameters.
In Eq. (6), extracted radius value \({r}_{\theta }\) is calculated first from the flame radius near the direction of the ignition electrode connection. The abrupt boundary of the mutation region caused by the ignition electrodes can be obtained from Eq. (6). Thus, the boundary angle (\({\theta }_{bk}\)) of each of the four regions is obtained from \(\Delta {r}_{max}\) in Fig. 12. It is calculated by
Removing the angle affected by the ignition electrodes.
The new measurement center Q′ is determined again according to the Han I method9 after removing the angle affected by the ignition electrodes in Fig. 13. Here, \(\theta ^{\prime}\) is defined as the measurement angle based on measurement center Q' in the two-dimensional coordinate system, and \({R}_{\theta }\)′ is flame radius after removal of ignition electrode effects corresponding to angle \(\theta ^{\prime}\). Thus, the mean radius calculated through the Han II method (\({Ra}_{2}\)) can be obtained.
Re-determining the measurement center.
To evaluate the optimization effect of the Han II method compared with the Han I method, the optimization rate of the mean radius (\(\Delta Ra /{Ra}_{1}\)) is used to describe them. The standard deviation optimization rate of the extracted radius value (\(\Delta \sigma /{\sigma }_{1}\)) and the range optimization rate of the extracted radius value (\(\Delta R/{R}_{1}\)) are used to describe the distribution rule. The difference of \(\Delta Ra /{Ra}_{1}\) with time is described by the optimization rate of the flame mean radius per unit angle (\({R}_{d}\)). The difference of \(\Delta \sigma /{\sigma }_{1}\) with time is described by the standard deviation optimization rate of the flame radius per unit angle (\({\sigma }_{d}\)). The specific evaluation of characteristic parameters is defined as follows:
Optimization rate of mean radius (\(\Delta Ra /{Ra}_{1}\)):
where \({Ra}_{1}\) and \({Ra}_{2}\) are the mean radius calculated by the Han I method and Han II method, respectively. This evaluation index is used to describe the optimization rate of the flame propagation radius.
where \({R}_{\theta }\) and \({R}_{\theta }^{\prime}\) are the extracted radius values calculated by the Han I method and Han II method, respectively. \(n\) is the total number of radius extracted values.
Optimization rate of the flame mean radius per unit angle (\({R}_{d}\)):
where \(\Delta Ra\) is the difference of the mean radius calculated by the Han I method and Han II method, respectively, and \({\theta }_{b}\) is the angle affected by the ignition electrodes. This evaluation index is used to describe the difference of the mean radius per unit angle.
Standard deviation optimization rate of the extracted radius value (\(\Delta \sigma /{\sigma }_{1}\)):
where \({\sigma }_{1 }\) and \({\sigma }_{2}\) are the standard deviation of extracted radius values calculated by the Han I method and Han II method, respectively. This evaluation index quantifies the improvement of the non-uniform flame propagation state.
Standard deviation of the extracted radius value (\(\sigma \)):
Standard deviation optimization rate of the flame radius per unit angle (\({\sigma }_{d}\)):
where \(\Delta\upsigma \) is the standard deviation of the extracted radius value calculated by the Han I method and Han II method, respectively, and \({\theta }_{b}\) is the angle affected by the ignition electrodes. This evaluation index is used to describe the standard deviation per unit angle.
Range optimization rate of the extracted radius value (\(\Delta R/{R}_{1}\)):
where \({R}_{1}\) and \({R}_{2}\) are the range of the extracted radius values calculated by the Han I method and Han II method, respectively. This evaluation index is used to describe the optimization effects when the flame presents both convex and concave distributions.
Range calculated through the Han I method (\({R}_{1}\)):
Range calculated through the Han II method (\({R}_{2}\)):
where \({R}_{\theta max}\mathrm{ and}\) \({R}_{\theta min}\) are the maximum and minimum extracted radius values calculated by the Han I method respectively. \({{{R}_{\theta }}^{\mathrm{^{\prime}}}}_{max}\) and \({{{R}_{\theta }}^{\mathrm{^{\prime}}}}_{min}\) are the maximum and minimum extracted radius values calculated by the Han II method, respectively.
Two typical working conditions in which the flame presents a convex or concave distribution are selected to conduct the adaptability analysis of the Han II method, as well as to investigate the affected rule of ignition electrodes on the flame propagation radius. Figure 14 shows the Schlieren images of the two typical working conditions.
Schlieren images of the two typical working conditions.
The ratio of the angle affected by ignition electrodes \(P\) is defined as follows:
where \({\theta }_{b}\) is the angle affected by the ignition electrodes.
Figure 15 displays the ratio of the angle affected by the ignition electrodes of the two typical working conditions. The flame propagation radius is 6–25 mm. The flame schlieren images with flame development time between 1.5–12 ms and 7.8–42.9 ms under convex or concave flame distributions were selected for analysis. Figure 15a shows the ratio of the angle affected by the ignition electrodes when the flame presented a convex distribution. It can be observed that the maximum ratio of the angle is 84–85% when the flame development time is 1.5 to 4.2 ms. As the flame continues to spread, the ratio of the angle decreases dramatically from 85 to 21% when the flame development time is 4.2 to 6.0 ms. However, the ratio eventually begins to slightly increase again from 21 to 30% when the flame development time is 4.2 to 6.0 ms. The ratio of the angle can be divided into three regions including a high value region, sharp variation region, and mild variation region for the flame that presents a convex distribution.
Ratio of the angle affected by the ignition electrodes of two typical working conditions. (a) Schlieren images of the flame presenting a convex distribution; (b) Schlieren images of the flame presenting a concave distribution.
With in-depth analysis of the Schlieren images showing the flame presenting a convex distribution, the boundary of the flame is affected by the ignition electrodes may have been misjudged for the different flame propagation stages, resulting in a large ratio of the angle affected by the ignition electrodes. The flame development time is 1.5 to 4.2 ms, and the flame front remains smooth for the high-value region, which may lead to a misjudgment of the boundary in Fig. 15a. Therefore, the ratio of the angle affected by the ignition electrodes in the high-value region is not reliable. For the sharp-variation region, the flame development time is 4.2 to 6.0 ms. The flame front surface near the ignition electrodes gradually appears slightly convex, which causes a decrease of the ratio of the angle affected by the ignition electrodes. The error of the boundary also gradually decreases. However, owing to the small value of \({\Delta r}_{max}\) (1.52 to 1.64 mm), the ratio of the angle affected by the ignition electrodes in this region is also not reliable. For the mild-variation region, the flame development time is 1.64 to 2.65 ms, and the slightly convex shape of the flame front surface near the ignition electrodes develops into an obviously convex shape. Here, \({\Delta r}_{max}\) is increased from 1.64 to 2.65 mm, which is sufficient to identify the boundary. Thus, the ratio of the angle affected by the ignition electrodes in the mild-variation region can be considered highly reliable. It can be concluded that the ratio of the angle affected by the ignition electrodes is only applicable to the mild-variation region when the flame presents a convex distribution.
Figure 15b shows the ratio of the angle affected by the ignition electrodes when the flame presents a concave distribution. The ratio of the angle can be divided into three regions include a mild-variation region, sharp-variation region, and high-variation region. The ratio of the angle affected by the ignition electrodes shows a relatively small reduction from 37 to 11% when the flame development time is 7.8 to 27.3 ms. This region is defined as a mild-variation region. As the flame continues to spread, the ratio of the angle increases dramatically from 11 to 98% when the flame development time is 27.3 to 29.9 ms. This region of the sudden change is defined as a sharp-variation region. The ratio of the angle affected by the ignition electrodes eventually increases to more than 95% when the flame development time is 29.9 to 42.9 ms. Hence, this region is defined as a high-value region. It can be concluded that the ratio of the angle is only applicable to the mild-variation region when the flame presents a concave distribution.
To systematically evaluate the accuracy of the flame radius after the removal of the ignition electrode effects in the convex distribution, the above evaluation of characteristic parameters are used to describe the flame propagation radius and the difference of optimization.
As the flame propagates, the optimization rate of the mean radius (\(\Delta Ra/{\Delta Ra}_{1}\)) tends to linearly decrease in Fig. 16a. The Schlieren images at five different moments within the flame time of 6 to 12 ms are selected for analysis. Figure 16b shows the mean radius \((Ra)\) and the optimization rate of the mean radius versus time when the flame presents a convex distribution. After removal of the ignition electrode effects, the mean radius increases, and the increment rate is between 0.4 and 0.85%. As the flame propagates, the optimization rate of the mean radius gradually decreases. The optimization rate of the flame mean radius per unit angle (\({R}_{d}\)) gradually decreases, which is why \(\Delta Ra/{\Delta Ra}_{1}\) decreases as the flame propagates in Fig. 16c.
Optimization rate of the mean radius when the flame presents a convex distribution. (a) Optimization rate of the mean radius; (b) Mean radius and the optimization rate of the mean radius; (c) Optimization rate of the flame mean radius per unit angle.
The standard deviation optimization rate of the extracted radius value (Δσ/\({\sigma }_{1}\)) tends to linearly increase as the flame propagates in Fig. 17a. In this study, Schlieren images at five different moments within the flame time of 6 to 12 ms are selected for analysis. Figure 17b shows the standard deviation of extracted radius value \(\sigma \) and the standard deviation optimization rate of the extracted radius value versus time when the flame presents a convex distribution. After removal of the ignition electrode effects, the standard deviation decreases, which implies that the inhomogeneity of the flame-radius extracted values decreases, and the reduced rate is between 11.91 and 22.1%. As the flame propagates, the standard deviation optimization rate of the extracted radius value gradually increases.
Optimization rate of standard deviation when flame presents convex distribution. (a) Standard deviation optimization rate of extracted radius value; (b) Standard deviation of extracted radius value and standard deviation optimization rate of extracted radius value; (c) Standard deviation optimization rate of flame radius per unit angle.
In Fig. 17c, the standard deviation optimization rate of the flame radius per unit angle (\({\sigma }_{d}\)) gradually increases as the flame propagates, which why \(\Delta \sigma /{\sigma }_{1}\) increases. This implies that the extracted radius values calculated by the Han II method are closer to the actual flame propagation radius. The convex distribution near the ignition electrodes is caused by the irregular shape of the initial flame core. However, the resistance caused by the flame-retarding effect near the ignition electrode surface inhibits the convex distribution near the ignition electrodes. Meanwhile, the stretching effect during flame development also inhibits a convex distribution near the ignition electrodes. With the flame development, the self-accelerating phenomenon aggravates the convex tendency near the ignition electrodes to a certain extent. In the process of the combined effects of those four factors, the self-accelerating effect of the flame is dominant and its effect is greater than the ignition electrode surface-retarding effect and the flame-stretching effect. The flame near the ignition electrodes presents a convex distribution and tends to change more intensely. The specific mechanism of the influence of those four factors on the flame shape near the ignition electrodes needs further study.
In Fig. 18a, the range optimization rate of the extracted radius value (\(\Delta R/{R}_{1}\)) tends to linearly increase as the flame propagates. Thus, Schlieren images at five different moments within the flame time of 6 to 12 ms are selected for analysis. Figure 18b shows the distribution of the maximum and minimum extracted radius values respectively calculated by the Han I method and Han II method when the flame presents a convex distribution. It can be observed that the maximum and minimum extracted radius values calculated by the Han I method are distributed within the angle affected by the ignition electrodes. However, the maximum and minimum extracted radius values calculated by the Han II method are distributed outside the angle affected by the ignition electrodes. Figure 18c shows the range of the extracted radius value (\(R\)) and the range optimization rate of the extracted radius value versus time when the flame presents a convex distribution. After removal of the ignition electrode effects, the range of the extracted radius value decreases, and the reduced rate of \(R\) is between 20.32 and 39.51%. As the flame propagates, the range optimization rate of the extracted radius value (\(\Delta R/{R}_{1}\)) gradually increases. This implies that the Han II method can effectively remove the ignition electrode effects.
Range of the extracted radius value when the flame presents a convex distribution.
As shown in Fig. 19a, as the flame propagates, the optimization rate of the mean radius (\(\Delta Ra/{\Delta Ra}_{1}\)) tends to linearly decrease. Thus, Schlieren images at five different moments within the flame time of 7.8 to 27.3 ms are selected for analysis. Figure 19b shows the mean radius \((Ra)\) and optimization rate of the mean radius versus time when the flame presents a concave distribution. After removal of the ignition electrode effects, the mean radius increases, and the increment rate is between 0.42 and 3.19%. As the flame propagates, the optimization rate of the mean radius gradually decreases. As shown in Fig. 19c, the optimization rate of the flame mean radius per unit angle (\({R}_{d}\)) gradually decreases, which is why \(\Delta Ra/{\Delta Ra}_{1}\) decreases as the flame propagates.
Optimization rate of the mean radius when the flame presents a concave distribution. (a) Optimization rate of mean radius; (b) Mean radius and optimization rate of mean radius; (c) Optimization rate of flame mean radius per unit angle.
As shown in Fig. 20a, the standard deviation optimization rate of the extracted radius value (Δσ/\({\sigma }_{1}\)) tends to linearly decrease as the flame propagates. In this study, Schlieren images at five different moments within the flame time of 7.8 to 27.3 ms are selected for analysis. Figure 20b shows the standard deviation of the extracted radius value \((\sigma \)) and the standard deviation optimization rate of the extracted radius value versus time when the flame presents a concave distribution. After removal of the ignition electrode effects, the standard deviation decreases, which implies that the inhomogeneity of the flame-radius extracted values decreases, and the reduce rate is between 5.13 and 17.99%. As the flame propagates, the standard deviation optimization rate of the extracted radius value gradually decreases.
Optimization rate of the standard deviation when the flame presents a concave distribution. (a) Standard deviation optimization rate of extracted radius value; (b) Standard deviation of extracted radius value and standard deviation optimization rate of extracted radius value; (c) Standard deviation optimization rate of flame radius per unit angle.
As shown in Fig. 20c, the standard deviation optimization rate of the flame radius per unit angle (\({\sigma }_{d}\)) gradually decreases as the flame propagates, which why \(\Delta \sigma /{\sigma }_{1}\) decreases. This implies that the extracted radius value calculated by the Han II method is closer to the actual flame propagation radius; however, it is slightly different from the variation trends that occur when the flame presents a convex distribution. The convex distribution near the ignition electrodes is caused by the irregular shape of the initial flame core. However, the resistance caused by the flame-retarding effect near the ignition electrode surface inhibits the convex distribution near the ignition electrodes. Meanwhile, the stretching effect during flame development also inhibits the convex distribution near the ignition electrodes. With the flame development, the self-accelerating phenomenon aggravates the convex tendency near the ignition electrodes to a certain extent. In the process of the combined effects of those four factors, the ignition electrode surface retarding and flame stretching are dominant. The convex tendency of the initial flame fire core is gradually inhibited until it reverses. The flame near the ignition electrodes displays a concave distribution, and it tends to slow under this working condition with the flame development. However, the specific mechanism of those four factors on the flame shape near the ignition electrodes needs further study.
In Fig. 21a, the range optimization rate of the extracted radius value (\(\Delta R/{R}_{1}\)) tends to linearly decrease as the flame propagates. Schlieren images at five different moments within the flame time of 7.8 to 27.3 ms are selected for analysis. Figure 21b shows the distribution of the maximum and minimum extracted radius values respectively calculated by the Han I method and Han II method when the flame presents a concave distribution. The maximum and minimum extracted radius values calculated by the Han I method are distributed within the angle affected by the ignition electrodes; nevertheless, the maximum and minimum extracted radius values calculated by the Han II method are distributed outside the angle affected by the ignition electrodes. Figure 21c shows the range of the extracted radius value (\(R\)) and range optimization rate of the extracted radius value (\(\Delta R/{R}_{1}\)) versus time when the flame presents a concave distribution. After removal of the ignition electrode effects, the range of the extracted radius value (\(R\)) decreases, and the reduced rate of \(R\) is between 0.32 and 8.09%. As the flame propagates, the range optimization rate of the extracted radius value (\(\Delta R/{R}_{1}\)) gradually decreases. This implies that the Han II method can effectively remove the ignition electrode effects, however, this trend is slightly different from the variation trends that occur when the flame presents a convex distribution. This is due to the different flame shapes resulting from the ignition electrode effects. The effect of the ignition electrodes is weakened as the flame propagates when presenting a concave distribution.
Range of extracted radius values when flame presents a concave distribution.
In this paper, a flame-radius calculation method is proposed to remove the effects of ignition electrodes. The adaptability and optimization effects of this method are analyzed. The main conclusions are as follows:
The ratio of the angle affected by the ignition electrodes calculated by the Han II method showed three obvious divisions: a high-value area, sharp-variation area, and mild-variation area. The ratio of the angle affected by the ignition electrodes is only applicable to the mild-variation region when the flame presented respective convex and concave distributions.
In the case of the convex distribution, after removing the ignition electrode effects, the increment rate of the mean radius is 0.4–0.85%. The reduced rate of the standard deviation of the extracted radius value is 11.91–22.1%. The decreased rate of the range of the extracted radius value is 20.32–39.51%.
In the case of a concave distribution, after removing the ignition electrode effects, the increment rate of the mean radius is 0.42–3.19%. The reduced rate of the standard deviation of the extracted radius value is 5.13–17.99%. Moreover, the decreased rate of the range of the extracted radius value is 0.32–8.09%.
The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.
Extracted Value of Fitting Radius
Characteristic Radius
Evaluation Parameters of Affected Angle Due to Ignition Electrodes
Boundary Angle
Affected Angle Due to Ignition Electrodes
Measurement Angle of Flame Radius
Initial Temperature
Initial Pressure
Equivalence Ratio
Diluent Gas Fraction
Constant Volume Combustion Vessel
Flame Radius after Removal of Ignition Electrode Effects
Mean Radius
Mean Radius Calculated through Han I Method
Mean Radius Calculated through Han II Method
Range of Extracted Radius Value
Range Calculated through Han I Method
Range Calculated through Han II Method
Standard Deviation of Extracted Radius Value
Standard Deviation Calculated through Han I Method
Standard Deviation Calculated through Han II Method
Ratio of Angle Affected by Ignition Electrodes
Optimization Rate of Mean Radius
Range Optimization Rate of Extracted Radius Value
Standard Deviation Optimization Rate of Extracted Radius Value
Optimization Rate of Flame Mean Radius Per Unit Angle
Standard Deviation Optimization Rate of Flame Radius Per Unit Angle
Huang, Z. et al. Measurements of laminar burning velocities for natural gas–hydrogen–air mixtures. Combust Flame 146(1–2), 302–311. https://doi.org/10.1016/j.combustflame.2006.03.003 (2006).
Article ADS CAS Google Scholar
Mohapatra, S. et al. Adaptability of different mechanisms and kinetic study of methane combustion in steam diluted environments. Sci. Rep. 12, 4577. https://doi.org/10.1038/s41598-022-08648-5 (2022).
Article ADS CAS PubMed PubMed Central Google Scholar
Cheng, Q. et al. Hyperspectral image reconstruction from colored natural flame luminosity imaging in a tri-fuel optical engine. Sci. Rep. 13, 2445. https://doi.org/10.1038/s41598-023-29673-y (2023).
Article ADS CAS PubMed PubMed Central Google Scholar
Wang, Q. et al. Suppression of deflagration flame propagation of methane–air in tube by argon gas and explosion-eliminating chamber. Sci. Rep. 12, 4965. https://doi.org/10.1038/s41598-022-09086-z (2022).
Article ADS CAS PubMed PubMed Central Google Scholar
Li, L. Research on determination of radii of spherical flames and computer realization. Ph.D. Thesis; Wuhan University of Technology (2014). https://doi.org/10.7666/d.D617237.
Li, H. et al. Effect of dilution on laminar burning characteristics of H2/CO/CO2/air premixed flames with various hydrogen fractions. Exp. Therm. Fluid Sci. 74, 160–168. https://doi.org/10.1016/j.expthermflusci.2015.11.026 (2016).
Article CAS Google Scholar
Chen, T. C. & Chung, K. An efficient randomized algorithm for detecting circles. Comput. Vis. Image Undergr. 83(2), 172–191. https://doi.org/10.1006/cviu.2001.0923 (2001).
Article Google Scholar
Broustail, G., Seers, P., Halter, F., Moreac, G. & Mounaim-Rousselle, C. Experimental determination of laminar burning velocity for butanol and ethanol iso-octane blends. Fuel 90(1), 1–6. https://doi.org/10.1016/j.fuel.2010.09.021 (2011).
Article CAS Google Scholar
Han, Z., Xiao, B., Tian, W., Li, J. & Yu, W. Study on difference and adaptability of calculation method of spherical flame radius. Fuel 259, e116216. https://doi.org/10.1016/j.fuel.2019.116216 (2020).
Article CAS Google Scholar
Bouvet, N., Chauveau, C., Gökalp, I. & Halter, F. Experimental studies of the fundamental flame speeds of syngas (H2/CO)/air mixtures. Proc. Combust. Inst. 33(1), 913–920. https://doi.org/10.1016/j.proci.2010.05.088 (2011).
Article CAS Google Scholar
Zuo, Z., Pei, Y., Qin, J., Xu, H. & Lu, L. Laminar burning characteristics of premixed methane-dissociated methanol–air mixtures under lean burn conditions. Appl. Therm. Eng. 140, 304–312. https://doi.org/10.1016/j.applthermaleng.2018.05.040 (2018).
Article ADS CAS Google Scholar
Gu, X. J., Haq, M. Z., Lawes, M. & Woolley, R. Laminar burning velocity and Markstein lengths of methane-air mixtures. Combust. Flame 121(1–2), 41–58. https://doi.org/10.1016/S0010-2180(99)00142-X (2000).
Article ADS CAS Google Scholar
Milton, B. E. & Keck, J. C. Laminar burning velocities in stoichiometric hydrogen and hydrogen hydrocarbon gas mixtures. Combust. Flame 58(1), 13–22. https://doi.org/10.1016/0010-2180(84)90074-9 (1984).
Article ADS CAS Google Scholar
Zhang, X., Hou, X., Wang, Y. & Zhang, J. Study on flame characteristics of low heat value gas. Energy Convers. Manag. 196, 344–353. https://doi.org/10.1016/j.enconman.2019.05.024 (2019).
Article CAS Google Scholar
DuttaRoy, R., Chakravarthy, S. R. & Sen, A. K. Experimental investigation of flame propagation in a meso-combustor. Proc. Inst. Mech. Eng. A J. Power Energy 234(8), 1131–1146. https://doi.org/10.1177/0957650919897755 (2020).
Article CAS Google Scholar
Gong, X., Huo, J., Ren, Z. & Law, C. K. Extrapolation and DNS-mapping in determining laminar flame speeds of syngas/air mixtures. Combust. Flame 200, 365–373. https://doi.org/10.1016/j.combustflame.2018.11.033 (2019).
Article ADS CAS Google Scholar
Tahtouh, T., Halter, F. & Mounaïm-Rousselle, C. Measurement of laminar burning speeds and Markstein lengths using a novel methodology. Combust. Flame 156(9), 1735–1743. https://doi.org/10.1016/j.combustflame.2009.03.013 (2009).
Article ADS CAS Google Scholar
Wu, F., Jomaas, G. & Law, C. K. An experimental investigation on self-acceleration of cellular spherical flames. Proc. Combust. Inst. 34(1), 937–945. https://doi.org/10.1016/j.proci.2012.05.068 (2013).
Article CAS Google Scholar
Li, G. et al. Experimental investigation on laminar burning velocities and markstein lengths of premixed methane–n-heptane–air mixtures. Energy Fuels 29, 4549–4556. https://doi.org/10.1021/acs.energyfuels.5b00355 (2015).
Article CAS Google Scholar
Wang, X. et al. Premixed combustion of spherically propagating methanol-air-nitrogen flames. Trans. CSICE 27(3), 207–214. https://doi.org/10.16236/j.cnki.nrjxb.2009.03.005 (2009).
Article Google Scholar
Han, Z., Xiao, B., Su, Q., Wu, X. & Qian, Y. Study on accuracy of measuring propagation radius of constant volume combustion flame. Trans. CSICE 37(5), 417–425. https://doi.org/10.16236/j.cnki.nrjxb.201905054 (2019).
Article Google Scholar
Mevel, R., Lafosse, F., Chaumeix, N., Dupre, G. & Paillard, C. E. Spherical expanding flames in H2–N2O-Ar mixtures: Flame speed measurements and kinetic modeling. Int. J. Hydrog. Energy 34(21), 9007–9018. https://doi.org/10.1016/j.ijhydene.2009.08.054 (2009).
Article CAS Google Scholar
Varea, E., Modica, V., Vandel, A. & Renou, B. Measurement of laminar burning velocity and Markstein length relative to fresh gases using a new postprocessing procedure: Application to laminar spherical flames for methane, ethanol and isooctane/air mixtures. Combust. Flame 159(2), 577–590. https://doi.org/10.1016/j.combustflame.2011.09.002 (2012).
Article ADS CAS Google Scholar
Xie, Y., Wang, J., Cai, X. & Huang, Z. Self-acceleration of cellular flames and laminar flame speed of syngas/air mixtures at elevated pressures. Int. J. Hydrog. Energy 41(40), 18250–18258. https://doi.org/10.1016/j.ijhydene.2016.07.239 (2016).
Article CAS Google Scholar
Law, C. K. & Sung, C. Structure, aerodynamics, and geometry of premixed flamelets. Prog. Energy Combust. 26(4–6), 459–505. https://doi.org/10.1016/S0360-1285(00)00018-6 (2000).
Article CAS Google Scholar
Download references
This work is supported by the National Defense Science and Technology Key Laboratory Fund Project [Grant No. 2022-JCJQ-LB-062-0102]; the National Natural Science Foundation of China [Grant No. 51776177]; the Sichuan Science and Technology Program [Grant No. 2023NSFSC0836]; and the Chunhui Plan of the Ministry of Education of the People's Republic of China [Grant Nos. HZKY20220586, HZKY20220569].
Key Laboratory of Fluid and Power Machinery of Ministry of Education, Xihua University, Chengdu, 610039, China
Likang Fan, Yan Yan, Zinong Zuo, Zhiqiang Han & Jia Fang
Key Laboratory of Vehicle Measurement, Control and Safety of Sichuan Province, Xihua University, Chengdu, 610039, China
Likang Fan, Xin Fu, Mingjie Hu, Yan Yan, Zinong Zuo, Zhiqiang Han, Jia Fang & Bang Xiao
Engineering Research Center of Ministry of Education for Intelligent Air-Ground Fusion Vehicles and Control, Xihua University, Chengdu, 610039, China
Zhiqiang Han
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
The authors all contributed to the deployment of the paper. Conceptualization, L.F.; Methodology, L.F. and Z.H.; Data curation, L.F., B.X., Z.Z. and X.F.; Formal analysis, L.F.; Project administration L.F. and Z.H.; Investigation, L.F., X.F. and B.X.; Resources, L.F., Z.H. and J.F.; Writing—original draft preparation, L.F.; Writing—review and editing, L.F., X.F., M.H., Y.Y., Z.Z., Z.H. and J.F.; All authors reviewed the manuscript.
Correspondence to Zhiqiang Han.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and permissions
Fan, L., Fu, X., Hu, M. et al. Calculation method of spherically expanding flame propagation radius to consider ignition electrode effects. Sci Rep 14, 8435 (2024). https://doi.org/10.1038/s41598-024-58940-9
Download citation
Received: 17 November 2023
Accepted: 04 April 2024
Published: 10 April 2024
DOI: https://doi.org/10.1038/s41598-024-58940-9
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative