1. Introduction
Studying a combustion process is a very complex task
[1]. However, the characterization of a flame is essential to improve the combustion process’s efficiency and reduce pollutant emissions. Measuring several variables in a flame is used to monitor the efficiency of combustion processes. Generally, some measured variables are flame temperature, soot concentration, radical intensities, flame size, luminosity, and flickering
[2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22]. Full-field optical techniques are usually used to make this kind of measurement because they give a global description of the process and do not disturb the sample. On the other hand, simultaneously measuring several variables is preferable because this avoids the lack of correlation between them
[9][16][17][18][19]. However, a complex optical setup is usually employed for measurement in such cases, making the task costly and challenging
[16][18].
The flame equivalence ratio (
𝛷) is an adimensional variable used to determine if a fuel-oxidizer mixture is rich, lean, or stoichiometric
[1]. When
𝛷>1, it denotes a fuel-rich mixing process; on the other hand,
𝛷<1 is related to fuel-lean mixtures, and the last
𝛷=1 in the ratio relationship is the value for stoichiometric mixing. This relation is intrinsically linked to determining a system’s performance. This study, as in others, corroborates that the maximum temperature value is obtained at a value of
𝛷=1.05 for LPG fuel. Other variables that can be connected directly to the flame equivalence ratio are the temperature and radical intensities (chemiluminescence)
[10][11][12][13][14][15][16][17][18][19]. In a combustion process, some released species are OH, CH, CH
2O, and C
2, which have a specific emission spectrum
[10][11][12][13][14][15][16][17][18][19]. In these combustion flames, the presence of radicals
C∗2 and
CH∗ are inherent in the reaction zone. These two chemical species are most abundant within the flame. Furthermore, these radicals emit in the visible region of the electromagnetic spectrum, centered at wavelengths 430 nm and 525 nm, respectively
[14]. For this reason, digital color cameras play an essential role in chemiluminescence measurements
[11][12][13][14][15][16][17][18]. In this analysis, each color channel contributes specific information about the object under study. However, most digital color cameras suffer from the crosstalk effect; this occurs when sensors have overlapping sensibilities, contaminating the recorded data and needing correction
[15][17].
There are few works in which chemiluminescence and flame temperature are measured
[3][16][17][18]. In one of the works, the optical system is simple; they only use a camera, and the temperature is measured using two-color pyrometry. However, temperature measurement is applicable only for sooty flames
[3]. In other research, two-color pyrometry is also used to determine the flame’s temperature; however, the optical arrangement is complex
[16][18]. In
[17], the optical configuration is simple to implement, and Digital Laser Speckle Displacement (DLSD) is used to determine the temperature fields. However, the chemiluminescence measurements are contaminated with speckle noise, making determining these values difficult.
2. Theoretical Development of Intensity Profiles on Flames
When a frame of information is registered on the camera plane, it corresponds to the projected intensity profile of the flame. In some studies of combustion analysis, it is a common topic of interest to determine the reconstructed profile
[23]. Therefore, homolog Abel deconvolution techniques can be used from speckle displacement to intensity. The direct and inverse Abel transforms are linked to analyze the desirable plane:
where
𝜌(𝑥,𝑧) is represented in rectangular reference coordinates and (
x,
z) is the profile intensity on the camera plane (see
Figure 1).
𝜌 can be the chemiluminescence and temperature values. Moreover,
𝐼(𝑟,𝑧) is the original profile intensity of the object in three dimensions in cylindrical coordinates; this is due to the nature of the axisymmetric flames. The main interest is obtaining the flame’s reconstructed profile intensity
𝐼(𝑟,𝑧). These profiles are contaminated with speckle phenomena and the crosstalk effect, which makes signal clean-up difficult.
Figure 1. Intensity profile projected on the camera plane from the phase object.
3. Algorithm for Profile Intensity Noise Reduction
Temperature and chemiluminescence values are obtained using DLSD. The data corresponding to chemiluminescence values require speckle noise reduction. Noise reduction is achieved using a four-order PDE to create a later curve fitting, utilizing a set of Gaussian functions.
1. Elimination of background noise: Firstly, it is necessary to perform a subtraction from the object image minus the reference image:
In Equation (3), the prefix
S refers to the subtracted image obtained from the object (
O) minus the reference image (
R); the superscript
rgb means that these are color images and all the channels are involved. In order to reduce the residual speckle presence in the neighborhood of the flame, a median filter is used. Afterward, a simple segmentation is performed to isolate the flame from the background:
where
N and
F are the resultant matrices with the background and flame information, respectively, and the superscript
𝑆𝑐 means this operation needs to be applied to each channel independently; here,
𝑡
is found by an iterative method
[20].
2. Filter using a four-order (PDE): The residual existence of speckles is characteristic as the median filter is insufficient to eliminate it.
Let the Laplacian of the intensity data be
∇𝜌𝑛𝑖,𝑗 and the subscripts
i,
j be the discretized values of the coordinates (
x,
z) from the data obtained on the camera, then:
where the superscript
𝑛 means that the target function
𝜌(𝑖,𝑗) is calculated by an iterative method, and
ℎ is the stepsize, in the first iteration, a seed value is necessary to obtain the next value on the iteration
𝑛+1; therefore, the recursive relation is:
where
𝜔 is a step size linked to the velocity of convergence in the solution, with a suggested value of
𝜔= 0.25.
3. Curve fitting by a set of Gaussian bases and deconvolution: Using a finite differences method implies a stagger effect on the pixel values. Therefore, a curve fitment adjustment improves integration performance over the intensity data
[24].
where
𝜌(𝑥,𝑧) is the profile to adjust for a fixed
z; in this case, for axisymmetric flames, a set of Gaussian functions can be used appropriately to fit the intensity profiles. Here,
𝑓𝑖(𝑥,𝑧)=𝑒𝑥𝑝(−(𝑥−𝜇 )22𝜎2), where
μ is the mean value of the data from the respective profile, and
𝜎 is the standard deviation of the intensity values. Equation (7) can be applied to all the profiles of the image and be represented as a matrix system:
The matrix system in (8) is a system of the type
Aw=B, where m is the profile position and can be solved by a matrix solution process to obtain the weights.
3. Abel Reconstruction
After profile noise reduction and curve fitting, 3D Abel transform reconstruction is applied to chemiluminescence and temperature values
[17]. Equation (3) is used to fulfill such a purpose. A numerical approximation proposal in N. A. Fomin is used to obtain an approximate reconstruction of the object intensity
[25].
4. Crosstalk Correction
Given the nature of the optical setup to perform temperature and chemiluminescence measurements, the construction of the color camera that uses a Bayer-type filter responsible for the sensibility is more prominent on the green channel than the other two in this device. Therefore, crosstalk correction is needed to obtain the trust intensity values on the object plane. Furthermore, due to the regulation of parameters such as exposure time, saturation, and others, an approximation of intensity values in the pixels can be proposed in the following way:
In this notation, the letter
I means that this correction applies to the reconstructed object intensity profiles. Moreover,
𝐼𝑖̂(𝑟,𝑧) in the equation is the corrected profile against the crosstalk effect for each RGB color, which means there is a proportionality for each pixel to denote the influence from the overlapped spectral response in the camera, the
𝐴𝑖,𝑗 coefficient denotes this; as long as
𝐼𝑗 is the uncorrected intensity registered on the camera sensor; here, the indices ‘
i’ and ‘
j’ refer to the color of the illumination source and each pixel, respectively. This equation can be represented in matrix form as follows.
The values of the
A coefficients can be obtained from the camera’s graph spectral response. It has been established that certain chemical species are related to efficiency in a combustion process. Two of them are analyzed due to the nature of the premixed flames. These radicals are CH* and
C∗2, which emit on the visible region of the electromagnetic spectrum around
~430 nm and
~515 nm, respectively
[14]. Green and blue channels of the color camera are used to measure them. On the other hand, the illumination source, the He-Ne laser, has a wavelength of 633 nm, and is related to calculating temperature profiles related to speckle displacements. Therefore, the effect or red illumination source is depreciated for crosstalk analysis, so subscripts 1 and 2 can be considered for chemiluminescence analysis, substituting them for G and B, respectively; the set of equations
[12] can be approximated to a system of two equations with two unknowns.
The previous equations are solved to find the solutions for intensity values from each color channel, as follows:
Therefore, the values of intensities
𝐼𝐺 and
𝐼𝐵 correspond to the chemiluminescence measurements.