Abstract

A coaxial dielectric barrier discharge plasma-assisted combustion actuator (DBD-PACA) system was set up to study its discharge and optical emission spectrum (OES) characteristics in space in this paper. Results showed that each discharge cycle can be divided into four stages: a, b, c, and d. Discharge-on only occurred in stages b and d. Comparatively, the discharge intensity was larger in stage d due to the memory effect of excited electrons. Moreover, Lissajous figure and current-voltage methods were utilized to calculate the power of the coaxial DBD-PACA, and both methods produced roughly similar results. The power presented an upward trend with increasing input voltage and airflow rate. In addition, numerous second positive system (SPS) excited nitrogen molecules were detected from the OES signals. The intensity of the spectral lines (297.54 nm, 315.76 nm, 336.96 nm, and 357.56 nm) first increased, then maintained, and then increased rapidly with the increased radius; however, the intensity of the spectral lines (380.34 nm, 405.80 nm, and 434.30 nm) basically remained unchanged, then increased, and finally decreased with the increased radius. The vibrational temperature first decreased quickly and then increased and reached the minimum at r = 18 mm with the increased radius. The vibrational temperatures at all collection points decreased with the increased input voltage. However, within the range of 0–280 L/min, when r was lower than 15 mm, the vibrational temperatures first increased rapidly and then decreased slowly; when r was greater than 15 mm, the vibrational temperatures first increased and then basically remained stable.

1. Introduction

Plasma, which is an equal density of negative and positive particles ionized gas, belongs to the fourth state of matter and is composed of charged particles (electrons and ions), atoms, radicals, molecules in ground and excited states, and photons [1, 2]. Plasma has been widely used in many fields, such as plasma-assisted combustion (PAC) [3], plasma surface modification [4], plasma sterilization [5], plasma active flow control [6], and plasma mutation breeding [7]. Among these, PAC is a remarkable solution to enlarging stable combustion boundary, improving combustion stability and combustion efficiency, and lowering nitric oxide emissions; it addresses these issues at a low energy cost, as put forward and studied by numerous relevant experts [811]. The superior performance of the PAC actuator (PACA) is the key of PAC technology, and numerous plasma discharge devices, such as dielectric barrier discharge (DBD) [1215], pulsed corona [16], microwave discharges [17], and gliding arc [18], were investigated to improve combustion. Retter et al. [19] designed a DBD actuator to enhance a hydrogen diffusion flame and studied the fundamental aspects of DBD plasmas such as the nonthermal self-limiting, and steamer was found with the presence of hydrogen. Mintoussov et al. [20] studied the characteristics of the DBD actuator at high temperature and found that the percentage of absorbed energy by DBD actuator grew with temperature and the energy input determined the gap infill. Tang et al. [21] studied the effects of airflows on a multicathode DBD and found that the discharge intensity increased along the direction of airflow.

Since DBD plasma can produce heat, electrons, excited molecules, long-lifetime intermediate species, radicals, ions, ionic wind, and Lorentz force, it improves combustion mainly through three pathways: thermal, kinetic, and transport [22]. Depending on the Arrhenius law, plasma produces heat and then improves temperature to accelerate combustion chemical reactions in the thermal pathway. Plasma also produces active radicals via direct electron impact dissociation, ion impact and recombination dissociation, and collisional dissociations of reactants with electronically excited and vibrationally excited molecules to lower the energy barrier of combustion chemical reactions and shorten the combustion chemical reactions chains in the kinetic pathway. Plasma can break down the large fuel molecules and produce ionic changing the local flow field and facilitating flow turbulization and mixing of the reactants in the transport pathway. Although lots of efforts have been done on PAC, qualitative or even quantitative understanding of the three pathways has not been clear. Meanwhile, few of PACAs have been used to enhance combustion performance through injecting ionized airflows and gaseous fuel into combustor. The coaxial DBD-PACA we designed is used in this way, and its characteristics need to be studied before being installed in combustor.

The purpose of PAC is to consume little energy to achieve combustion improvement, so studying discharge characteristic of the PACA is essential, and the common diagnostic method is to obtain the voltage and current signals which can deduce current density, discharge pulse width, microdischarge, and so on. However, the breakdown and extinction voltages and the power consumption are more important discharge parameters for PACA applied to aero-engine combustor [23]. Breakdown and extinction voltages are used to evaluate the stable operating range of the PACA, and power consumption is used to evaluate the performance of the PACA. Therefore, the influences of input voltage and airflow rate on the breakdown and extinction voltages and power consumption were studied in this paper.

The optical emission spectrum (OES) is a common and effective method to diagnose the kinetic effect of PACA. Radiative transitions relative to OH, CN, H, N, O, N2, C, and NO were identified in the PAC via the OES, and when the equivalence ratio increased from 0.4 to 1.4, the emissions respective to NO and C decreased, while the OH band intensity increased [10]. Lin et al. [24] applied the OES to study gliding arc plasma in gliding arc fuel injector and found that the significant thermal effect was caused by the arc column. To the best of the authors' knowledge, most researchers only study the spectra of a single point in the discharge area; thus, the results obtained were one-sided and incomplete. In this paper, the OES signals were captured at ten points with equal spacing along the radial direction of the coaxial DBD-PACA, and the results were more reliable and complete.

In the current study, a coaxial DBD-PACA experimental facility was established, and an oscilloscope was used to capture the discharge signals of the coaxial DBD-PACA. Meanwhile, a spectrograph was used to capture the OES signals of the plasma produced by the coaxial DBD-PACA. The major research focus of this work is to study the discharge and OES characteristics of the coaxial DBD-PACA and the effect of airflow rate and input voltage on the power calculated by Lissajous figure and current-voltage methods, analyze the spatial distribution of the spectrum, and obtain the trend of the vibrational temperature of the whole discharge area with increased airflow rate and input power.

2. Experimental Setup and Theory Analysis

2.1. Experimental Setup

Figure 1(a) shows a schematic of the experimental system. The experimental apparatus consisted of four primary parts: coaxial DBD-PACA, gas supply system, power supply, and measurement system. The inner electrode of coaxial DBD-PACA was a stainless steel tube with a 20 mm outer diameter and functioned in transporting gaseous fuel. The inner diameter of the quarter glass tube was 38 mm, whereas quarter glass thickness was 2 mm. The copper foil electrode adhered to the quarter glass tube and five layers of Kapton covering to prevent the glow discharge. For uniform discharge, the four parts of coaxial DBD-PACA should maintain concentricity, and the top view of the DBD-PACA is illustrated in Figure 1(b). In order to obtain homogeneous discharge plasma, a swirler was installed between the quartz glass tube and the inner electrode, and it was located at the upstream of the discharge area. Figures 1(c) and 1(d) show the discharge diagrams without and with swirler. The coaxial DBD-PACA was driven by an alternating current single high-voltage differential power (CTP-2000K), and this power supply was manufactured by Suman Electronics. The high-voltage end of the power supply was attached to the outer copper foil electrode and controlled by the voltage regulator, whereas the other end was connected to ground lead. The working frequency of the power supply was set to 9.5 kHz. The essence of the power supply is to increase the input voltage by thousands times to the applied voltage, and the applied voltage is controlled by the input voltage. The gas supply system mainly included a piston air compressor, and the flowmeter produced by Sevenstar was utilized to control the airflow rate flowing into the coaxial DBD-PACA.

The measurement system mainly consisted of a spectrograph, oscilloscope, and Canon digital signal lens reflex (DSLR) camera. The spectrograph (AvaSpec-ULS2048-4-USB2) produced by Avantes contained four channels to capture spectrum signals with the wavelength ranging from 200 nm to 950 nm and a slit width of 10 μm. The probe to collect the spectrum of the spectrograph was fixed on the slideway. And the spectra were collected at ten locations (r = 10 mm, 11 mm, 12 mm, 13 mm, 14 mm, 15 mm, 16 mm, 17 mm, 18 mm, and 19 mm), as shown in Figure 1(b). Tektronix P6015A (1000 : 1) high-voltage probe was used to measure the applied voltage of the coaxial DBD-PACA; Tektronix TPP1000 (10 : 1) low-voltage probe was used to measure the applied voltage of the monitor capacitor; Tektronix TCP0030 current probe (1 : 1 for measurement at 5 A and 10 : 1 for measurement at 30 A; bandwidth, 120 MHz) was used to measure the current of the series circuit; all discharge signals were recorded by the oscilloscope (DPO4104B; Tektronix). Direct discharge plasma images were captured by a DSLR camera (Canon EOS 80D).

2.2. Power Measurement Methods

Two methods were used to calculate the power consumption of the coaxial DBD-PACA: Lissajous figure method and current-voltage method [25].

2.2.1. Lissajous Figure Method

The power was estimated through calculating the area of charge-voltage curve of the coaxial DBD-PACA. The key to this method was to connect a monitor capacitor () with a value much larger compared with the capacitance of the coaxial DBD-PACA without plasma discharge to the series circuit (Figure 2) [26]. And the value of is 1000 pF in this experiment. Therefore, the current passing through the monitor capacitor was computed as follows:

Instantaneous charge on the monitor capacitor was computed as follows:

The capacitance of the monitor capacitor was assumed to be constant over the applied voltage and frequency range of the coaxial DBD-PACA, thus obtaining the following:

As the current through the monitor capacitor must be identical to that passing through the coaxial DBD-PACA (); the instantaneous power dissipated by the coaxial DBD-PACA can be written as follows:

So, the average power P over one discharge cycle period T was calculated with the following formula [27]:

Therefore, the instantaneous transport charge Qm and the instantaneous applied voltage VPACA plotted against each other generated a Lissajous curve in the Qm-VPACA plane. And the area inside the closed Lissajous curve divided by one discharge cycle period T was equal to the power of the coaxial DBD-PACA.

2.2.2. Current-Voltage Method

According to electrical properties, instantaneous power consumption can be calculated by multiplying the current-voltage signals on the coaxial DBD-PACA:where denoted the instantaneous applied voltage on the coaxial DBD-PACA, represented the coaxial DBD-PACA instantaneous current, and referred to the instantaneous power.

Current-voltage signals consist of numerous discontinuous points. Thus, the difference and summation method was introduced to the instantaneous power. The power was then computed as follows:where N specified the total signal points in one discharge cycle and Δt was the time interval between two adjacent signal points.

2.3. Vibrational Temperature Calculation Method

The particles in the excited vibrational state were easy to absorb the chemical reactants, which can overcome the reaction barrier, reduce the energy consumption, and promote the chemical reaction. The vibrational temperature described the intensity of molecular vibration excitation, which can be calculated by the multispectral line slope method. The spectral line intensity of the vibration system of diatomic molecular emission spectrum waswhere and represent the vibration quantum number of the upper and lower state, respectively, was Planck’s constant (h = 6.62606975 × 10−34 J.s), was the speed of light, was the emission spectral frequency, was the transition probability between the upper and lower energy levels, and was the molecular number of the upper state.

In the same electronic excited state, the number of molecules with different vibration-excited states waswhere N0 indicated the number of the excited vibrational ground state in this electronic excited state, k was a constant (k = 8.617689 × 10−5 eV/K), and was the vibrational temperature.

According to spectroscopy, the vibrational energy of the upper state can be expressed as follows:where was the vibrational frequency of the harmonic oscillator, was the first-order correction of the anharmonic vibration, and was the second-order correction of the anharmonic vibration [28]. The following equation can be obtained by combining the above three equations:

Therefore, was a linear equation of the (), so the vibrational temperature of the plasma was the reciprocal of the slope of this equation. However, the fitting error exists in the vibrational temperature calculated by this method, and the error of spectral intensity also leads to the systematic error of vibrational temperature. In this paper, three groups of vibration-excited state bands, Δν = −1 (0-1, 1-2, 2-3), Δν = −2 (0–2, 1–3, 2–4), and Δν = −3 (0–3, 1–4, 2–5), were substituted into the above equation to calculate the vibrational temperature, respectively, and the average value was taken as the vibrational temperature of the plasma to reduce any discrepancies. Table 1 shows some spectral parameters of the SPS of nitrogen molecules used for calculating the vibrational temperature.

3. Results and Discussion

3.1. Discharge Characteristic

Figure 3 shows typical applied voltage, current, and charge waveforms with airflow rate of 200 L/min and input voltage of 130 V. Based on waveforms, the applied voltage and charge of the coaxial DBD-PACA discharge presented sinusoidal shapes, which were related to power supply characteristics. However, the current contained a series of intermittent high-amplitude spikes. The coaxial DBD-PACA discharge contained numerous filament channels (micro-discharge) with continuous formation and extinguishment, which resulted in the appearance of spikes [30]. The phase difference between the applied voltage and charge was approximated 6 μs. Each discharge cycle can be roughly divided into four stages: a, b, c, and d. No spike (discharge-off) was observed in stages a and c. On the other hand, numerous spikes were noted in stages b and d, and the quantity and intensity were both considerably larger in stage d. The reason was attributed to the massive electrons attached to the quartz dielectric inner surface and formed a built-in electric field [31] with the direction opposite to the applied electric field in stage b, thus weakening the intensity of total electric field; however, both directions were same in stage d, thus strengthening the intensity of total electric field.

Figure 4 showed quarter discharge cycle waveform, which was the partial enlarged version of the waveforms in the green oval in Figure 3. The appearance of every spike caused a decrease in applied voltage across the coaxial DBD-PACA and increase in charge. When one discharge channel was broken down, resistance of the coaxial DBD-PACA decreased. However, the resistance of the monitor capacitor remained constant. Thus, the applied voltage across the coaxial DBD-PACA reduced. The circuit charge increased with increasing current. However, the number of spikes barely increased with increasing applied voltage amplitude.

Extinction and breakdown input voltages were essential characteristics of the coaxial DBD-PACA [32]. Figure 5 shows the extinction input voltage and breakdown input voltage versus the airflow rate. The data were the average of five experimental results under the same experimental conditions. The breakdown input voltage was notably greater than the extinction input voltage due to the energy of the coaxial DBD-PACA originated entirely from the power supply before breakdown of air; however, during the coaxial DBD-PACA discharge, a low voltage was needed to maintain discharge under the action of electrons memory effect and built-in electric field. Moreover, the breakdown input voltage and extinction input voltage both increased with airflow rate as the aerodynamic force increased the relative minimum breakdown distance between two electrodes. Thus, a higher input voltage was needed to break down or maintain the microdischarge channels. The minimum extinction input voltage and breakdown input voltage reached 72.7 V and 77.0 V, respectively, in still air. The maximum extinction input voltage and breakdown input voltage measured 87.8 V and 91.0 V, respectively, when the airflow rate was 280 L/min.

Figure 6 shows the five cycles of Lissajous figures corresponding to the waveforms of the applied voltage and charge in Figure 3. Lissajous figures presented an almost-parallelogram shape, and the four edges represent the four stages of the coaxial DBD-PACA discharge. Stages a and c were in discharge-off state, but, stages b and d were in the state of discharge-on. Smooth Lissajous figures were obtained for stages a and c. However, Lissajous figures fluctuated in stages b and d, and fluctuation became more severe in stage d. When the coaxial DBD-PACA was discharging, negatively charged electrons were produced and moved along the direction of electric field, and wall charges appeared on the quarter dielectric inner side. When the applied voltage was positive, the electric field induced by wall charges was contrary to the external electric field, thus weakening the discharges. By contrast, when the applied voltage was negative, the direction of the electric field induced by wall charges was the same as that of the external electric field, thus strengthening the discharges.

The red loop was the average value of five cycles of Qm-VPACA Lissajous figures in Figure 6. In this study, to avoid systematic and experimental errors, the red loop was used as standard to calculate the area of Lissajous figure, which was used to obtain the power of the coaxial DBD-PACA.

Figure 7 shows the Lissajous figures and the power with the input voltage ranging from 90 V to 190 V at the condition of the airflow rate being 200 L/min. Figure 7(a) showed that the area of the Lissajous figure and maximum transport charge distinctly increased with increasing input voltage; however, the maximum applied voltage of the coaxial DBD-PACA increased slightly with increasing input voltage. So the input voltage has little effect on the maximum applied voltage, but the input voltage has a significant effect on the maximum transport charge. This phenomenon occurred as the current increased with increasing input voltage in the circuit; then more microdischarge channel breakdowns were observed in the coaxial DBD-PACA, and its equivalent resistance decreased. Thus, the assigned applied voltage of the monitor capacitor evidently increased. And depending on equation (2), the maximum transported charge also deservedly increased. Figure 7(b) demonstrated that the power calculated by Lissajous figures method and current-voltage method both increased with increasing input voltage. The power calculated by the two methods presented differences caused by the theoretical errors.

Figure 8 shows the Lissajous figures and the power with the airflow rate from 0 L/min to 280 L/min at the condition of the input voltage being 130 V. Based on Figure 8(a), the area of the Lissajous figure fluctuated slightly with increasing airflow rate. As shown in Figure 8(b), the power followed a macroscopic upward trend along with fluctuation within a narrow range. The turbulence of airflow possibly caused the fluctuation of the power. As the excited particles and heat energy produced by the coaxial DBD-PACA discharge were blown away and given the presence of aerodynamic force, thus the microdischarge channels lengthened and the plasma temperature declined. Therefore, the applied voltage increased, and the quantity of the microdischarge channels rose and thus the power increased.

3.2. OES Characteristic
3.2.1. Emission Spectra of the Coaxial DBD-PACA

As a nonintrusive diagnostic method, the OES was widely used to diagnose plasma [3337]; it can obtain important plasma data, such as the kinds and concentrations of excited particles, electron temperature, rotational temperature, and vibrational temperature. In this part, the integration time of the four channels of the spectrograph was set to 200 ms and the acquisition spectrum denoted the average of 20 cycles.

Figure 9 shows the typical coaxial DBD-PACA OES under the experimental conditions of 130 V input voltage and 200 L/min airflow rate. As depicted in Figure 9, the spectral lines were mainly composed of numerous nitrogen molecules SPS (), several nitrogen ions first negative system (FNS, ), and nitric oxide γ system (). The SPS of nitrogen was produced by the N2 (C3Πu) and N2 (B3Πg) transitions, which was the strongest and most stable system in the emission spectrum of air. Since the ionization excitation potentials of N2 (B3Πg) and N2 (C3Πu) were 7.35 eV and 11.03 eV, respectively, as shown in equations (12) and (13), there were free electrons with energy higher than 11.03 eV in the plasma generated by the coaxial DBD-PACA discharge. As the excited nitrogen atoms had a short lifetime and fleetly recombined into nitrogen molecules, no spectral lines of excited nitrogen atoms were detected [38]. As ionization excitation potential of particle was as high as 18.76 eV and the energy level transition was mainly due to the collision of free electrons, thus the weak relative intensity of FNS spectral lines indicated that a notably low number of electrons exhibited an energy higher than 18.76 eV in the plasma. The most common excited particles in oxygen were O2 (a1Δg), O2 (b1Σg+), and O2 (A3Σu+), whose excitation energies were 0.977 eV, 1.627 eV, and 4.43 eV, respectively, and can be produced by the collision between electrons and oxygen molecules, as shown in (14)–(16). These particles, however, have very low dissociation energies of 4.54 eV, 3.89 eV, and 1.15 eV, respectively, so they can easily collide with high-energy electrons to form ground state oxygen O (3P) and metastable oxygen O (1D), O (1S). Therefore, no oxygen molecules and atoms in the excited state were observed in the plasma generated by the coaxial DBD-PACA discharge. Next, the spectral lines of the SPS of nitrogen molecules with wavelength between 290 nm and 440 nm were studied.

In the experiment, the discharge plasma spectra were collected at ten points with equal spacing along the radial direction of the coaxial DBD-PACA, as shown in Figure 10. The spectral intensity along the radial direction varied obviously in the discharge area. As a whole, the intensity of the spectra increased in the direction away from the inner electrode. Next, the trends of the seven representative spectral lines (297.54 nm (2, 0), 315.76 nm (1, 0), 336.96 nm (0, 0), 357.56 nm (0, 1), 380.34 nm (0, 2), 405.80 nm (0, 3), 434.30 nm (0, 4)) of the SPS of nitrogen molecules were selected to study, as shown in Figure 11. As the radius increased, the intensity of spectral lines (297.54 nm, 315.76 nm, 336.96 nm, and 357.56 nm) first increased, then remained, and then increased rapidly and reached the maximum value at r = 19 mm. However, as the radius increased, the intensity of spectral lines (380.34 nm, 405.80 nm, and 434.30 nm) basically remained unchanged, then increased, and finally decreased and reached its maximum value at r = 16 mm. And the intensity of these three spectral lines almost disappeared at the position of r = 19 mm. In the SPS of nitrogen molecules, the spectral line intensity of the 336.96 nm was the strongest at almost all discharge regions of the coaxial DBD-PACA. This was mainly because the ground state was the easiest to obtain and the most abundant of all vibrational level states of the excited particles N2 (C3Πu); meanwhile, the transitions between the ground state of N2 (C3Πu) and N2 (B3Πg) were the easiest transition path. However, the spectral line intensity of the 380.34 nm was the strongest at radiuses of 10 mm and 16 mm, respectively.

3.2.2. Effects of Input Voltage on Vibrational Temperature

The OES were collected at ten equidistant points in the radial direction of the coaxial DBD-PACA with different input voltages. And according to the multispectral line slope method mentioned above, when the airflow rate was 200 L/min, the vibrational temperatures at all collection points were calculated, as shown in Figure 12. It can be observed from Figure 12(a) that the vibrational temperature had the same tendency in the radial direction of the coaxial DBD-PACA when the input voltage changed, which first decreased quickly and then increased, and reached the minimum at r = 18 mm. Figure 12(b) shows that the vibrational temperature at each collection point reduced but slightly fluctuated, with the increase of the input voltage. It was found that vibrational excitation dominated the input energy coupling when the reduced electric field was small, while electronic excitation dominated when the reduced electric field was great by solving the Boltzmann equation [39]. And when the coaxial DBD-PACA driven by a single-high-voltage power supply discharged, the closer to the inner electrode connected to the ground, the smaller the reduced electric field; and the closer to the outer electrode connected to the high voltage, the greater the reduced electric field [40]. Thereby the vibrational temperature decreased with the increase of the input voltage.

3.2.3. Effects of Airflow Rate on Vibrational Temperature

The OES were collected at ten equidistant points in the radial direction of the coaxial DBD-PACA at different airflow rates. And Figure 13 shows the vibrational temperature at all collection positions at the condition of the input voltage being 130 V. As the radius increased, the vibrational temperature first decreased and then increased and reached the minimum value at r = 18 mm, as shown in Figure 13(a). As we can see from Figure 13(b), within the range of 0–280 L/min, when r was less than 15 mm, the vibrational temperature first increased rapidly and then decreased slowly; when r was greater than 15 mm, the vibrational temperature first increased and then basically remained stable. As the discharge was nonuniform in still air, the vibrational temperature reduced obviously when the airflow rate was 0 L/min. The smaller the radius was, the more the vibrational temperature was influenced by the airflow rate.

4. Conclusions

In this paper, a coaxial DBD-PACA system was established to study its discharge and OES characteristics. Lissajous figure method and current-voltage method were used to calculated the power consumption of the coaxial DBD-PACA. And the vibrational temperature was obtained by means of fitting spectral intensity of nitrogen molecule SPS; meanwhile, the influences of input voltage and airflow rate on vibrational temperature were studied. The main conclusions were as follows.

It was found that the breakdown input voltage and extinction input voltage both increased with increasing airflow rate, and breakdown input voltage was notably greater than extinction input voltage. The coaxial DBD-PACA only discharged in the stages b and d, and the discharge became more intense in stage d. The powers measured by Lissajous figure method and current-voltage method were similar, and the power increased with increasing input voltage and airflow rate.

A number of SPS excited nitrogen molecules were detected in the discharge plasma. And the trends of seven representative spectral intensities of the SPS nitrogen molecules were different in different discharge region. Meanwhile, the vibrational temperature varied significantly in different discharge region. The vibrational temperature decreased with the increase of the input voltage; however, the vibrational temperature first increased, then remained, and lastly decreased with the increase of the airflow rate.

Furthermore, the results obtained revealed the performance of the coaxial DBD-PACA. However, the ultimate goal is to install the coaxial DBD-PACA on the aero-engine to assist combustion. Therefore, the coaxial DBD-PACA will be installed on the combustion chamber to study the influence of the coaxial DBD-PACA on combustion in the future.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was funded by the Key Project of National Natural Science Foundation of China (Grant no. 51436008) and the Youths Science Foundation Project of National Natural Science Foundation of China (Grant no. 51806245).