Next Article in Journal
A Thermal Discomfort Index for Demand Response Control in Residential Water Heaters
Next Article in Special Issue
Applications of Plasma Produced with Electrical Discharges in Gases for Agriculture and Biomedicine
Previous Article in Journal
Predicting ICD-9 Codes Using Self-Report of Patients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Carbon Dioxide Decomposition by a Parallel-Plate Plasma Reactor: Experiments and 2-D Modelling

1
Faculty of Physics, Shahid Bahonar University of Kerman, Kerman 7616913439, Iran
2
Department of Physics, Vali-e-Asr University of Rafsanjan, Rafsanjan 7718897111, Iran
3
Department of Physics, University of Cordoba, E-14071 Cordoba, Spain
4
Department of HVAC Engineering, Bialystok University of Technology, 15-351 Bialystok, Poland
5
Department of Plasma Engineering, Graduate University of Advanced Technology, Kerman 7616913439, Iran
6
Department of Pilotage, Mahanair UAST University of Kerman, Kerman 7616913439, Iran
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2021, 11(21), 10047; https://doi.org/10.3390/app112110047
Submission received: 2 October 2021 / Revised: 21 October 2021 / Accepted: 23 October 2021 / Published: 27 October 2021
(This article belongs to the Special Issue Novel Applications of Plasma Techniques for the Environment)

Abstract

:

Featured Application

This research is part of a project which aims to investigate the applicability of plasma technologies to air quality. This paper describes the development of a new plasma reactor for CO2 decomposition that works at atmospheric pressure with AC power. Results for pure CO2 are presented that elucidate the main mechanisms causing splitting of this molecule in the developed reactor.

Abstract

The applicability of high voltage electrical discharges for the decomposition of CO2 has been extensively demonstrated. In this study, a new AC parallel-plate plasma reactor is presented which was designed for this purpose. Detailed experimental characterization and simulation of this reactor were performed. Gas chromatography of the exhaust gases enabled calculation of the CO2 conversion and energy efficiency. A conversion factor approximating 25% was obtained which is higher in comparison to existing plasma sources. Optical emission spectroscopy enabled the determination of the emission intensities of atoms and molecules inside the plasma and characterization of the discharge. The Stark broadening of the Balmer hydrogen line Hβ was used for the estimation of the electron density. The obtained densities were of the order of 5 × 1014 cm−3 which indicates that the electron kinetic energy dominated the discharge. The rotational, vibrational, and excitation temperatures were determined from the vibro-rotational band of the OH radical. A 2-temperature plasma was found where the estimated electron temperatures (~18,000 K) were higher than the gas temperatures (~2000 K). Finally, a 2-D model using the fluid equations was developed for determining the main processes in the CO2 splitting. The solution to this model, using the finite element method, gave the temporal and spatial behaviors of the formed species densities, the electric potential, and the temperatures of electrons.

1. Introduction

The mitigation of greenhouse gases is one of the most important challenges in the world today. In consequence, finding solutions for carbon dioxide (CO2) reduction is a priority research topic. Technologies for the decomposition of these gases after their generation are being developed, together with solutions to reduce CO2 production at source.
The applicability of high voltage electrical discharge plasmas to reform carbon dioxide (CO2) as an environmental pollutant has been extensively studied [1,2,3,4,5]. Plasma discharge reactors have shown improved efficiency and fewer side effects compared with other existing methods, such as steam methane reforming, steam methane/oxygen reforming, partial oxidation, catalytic partial oxidation, auto-thermal reforming, etc., [5].
Carbon dioxide is a stable and inert molecule that needs a large amount of energy for its splitting which means low energy conversion efficiency. Plasma (ionized gas) is an active medium containing different species including neutral species, ions, electrons, excited species, etc. The ionized species are accelerated by an applied electromagnetic field. Electrons with small mass can reach high energy and produce CO2 conversion by different impact processes including ionization and vibrational or electronic excitation. New species appear that can favor this conversion. In this way, the energy efficiency is substantially improved [6]. The efficiency of this conversion can be improved with the synergy of a co-reactant that has a higher (less negative) Gibbs free energy; CH4 (ΔGº = −50.7 kJ mol−1) and H2 (ΔGº = 0 kJ mol−1) are the candidates most commonly selected for this purpose [7].
Li et al. [8] studied the variation of CO2 and CH4 conversion with the CH4/CO2 ratio using an atmospheric pressure DC corona discharge. A CO2 conversion of 70% was found for a 2:1 ratio.
A conversion factor approximating 90% was obtained by Li et al. [9] using an atmospheric pressure glow discharge plasma (APGD). In this case, the CH4/CO2 ratio was 4:6. The proposed reactor presented the advantage of large scale treatment and high conversion ability.
In addition, a toroidal transformer-coupled plasma (TCTP) source was used for CO2 and CH4 conversion by Li et al. [10]. Optical and mass spectrometric measurements of this source was performed.
Yu et al. [11] evaluated a dielectric packed-bed plasma reactor for CO2 conversion. The authors showed that the dielectric properties and morphology of packing dielectric pellets notably influenced the electron energy distribution in the formed plasma discharge and the reactions inside the plasma reactor. The effects of the reverse reactions in the CO2 decomposition and the oxidation of CO were examined. Wang et al. [12] developed a plasma reactor for CO2 reforming based on dielectric barrier discharge (DBD) with CH4 and a catalyst. The influences of the different species formed inside the reactor on CO2 reforming were studied. It was found that the catalyst could substantially improve reduction in CO2 concentration. The conversion of CO2 into more valuable chemical products using catalytic plasmas was studied experimentally by Liu et al. [13]. Their results showed that the CO2 plasma discharges could generate oxygen and other active plasma species for further reaction. In some cases, these reactions lead to the formation of more valuable chemicals such as ethylene, propylene and oxygenates. The experiments also confirmed that the CO2 plasma was an excellent “catalyst” for the conversion of low alkanes to alkenes. This method was found to be an effective method for the utilization of CO2 and low alkanes.
This study aims to contribute to the development of new plasma technologies for CO2 conversion. An AC parallel-plate plasma reactor (AC-PPP) for CO2 remediation is presented that works at atmospheric pressure and uses alternating current (AC). This reactor has a very simple and low-cost design that can operate at atmospheric pressure, allowing straightforward scaling up for industrial applications. The new design is based on a high voltage (HV) discharge between two parallel electrodes where inlet and outlet metal pipes have been added. This makes it possible to extend the electromagnetic field inside these pipes and expand the treatment area, consequently increasing the conversion efficiency.
A complete experimental and theoretical study of this new reactor was performed to determine its CO2 remediation ability. The CO2 conversion factor, CO and O2 selectivity, and energy efficiency were determined by analyzing the exhaust gases using gas chromatography (GC).
Optical emission spectroscopy confirmed the CO2 decomposition in this reactor. The former species formed inside the discharge were detected. This technique also made it possible to determine the characteristic parameters of the plasma including electron density and rotational, vibrational, and excitation temperatures. In consequence, it was possible to determine if the plasma conditions were adequate for the application.
As indicated, the performance of this device was influenced by the spatial distribution of discharge inside the metallic inlet and outlet pipes. A 2-D fluid model was developed to simulate the spatial and temporal behavior of the plasma in all positions of the new AC-PPP reactor. This model allowed determination of the main kinetic mechanisms of CO2 decomposition in this reactor. The spatial distributions and temporal variations of the different species formed from the CO2 decomposition were calculated.

2. Materials and Methods

2.1. Experimental Set-Up

A scheme of the experimental arrangements with the AC-PPP reactor used in this study is shown in Figure 1.
For the manufacture of this reactor, two copper pipes (of 10 mm diameter) were placed axially on both sides of the reactor structure and were connected to the two ring disks (of 50 mm diameter). The pipes performed the role of both gas inlet/outlets and electrodes which were connected to the AC power supply (0–22 kV). The two dielectric cylinders (of 50 mm diameter) were applied to fix these copper pipes. These dielectric cylinders were joined by a Pyrex tube. The axial distance between electrodes was kept at 10 mm (see Figure 2).
A strong electric field was created between the copper electrodes when applying the AC high voltage. The applied voltage varied between 10–22 kV. In this way, a plasma was generated between the electrodes at atmospheric pressure.
The gas entered via one pipe and exited by a separate one. The flow rate of CO2 gas (purity > 99.5%) was low and constant (~0.15 SLM) to obtain improved control of the applied field. Therefore, the derived axial velocity could be neglected compared to other velocities, i.e., drift and diffusion velocities (Vdiff > Vdrif >> Vaxi).
The application of an electric potential between the two electrodes led to the formation of an axial electric field. Subsequent electron acceleration led to the ionization of background carbon dioxide. During the ionization process, free radicals, ions, and neutral atoms were created inside the AC-PPP reactor. The operating situation of the AC-PPP reactor for a flow rate of 0.15 SLM and an applied voltage fixed of 22 kV is shown in the photographic images in Figure 2.
The consumed power in the discharge was measured with a current probe (Tektronix. Beaverton, OR, USA, TCP202) and a high voltage probe (Tektronix, P6015A) which were connected to an oscilloscope (HM1508).
An HR4000 Ocean Optic spectrometer was used to detect the spectra and perform optical emission spectroscopy (OES). This spectrometer is a compact and flexible device with no moving parts and a linear CCD-array of 3648 pixels. It stores a full spectrum every millisecond with a wavelength range of 200–1100 nm. The optical features of this spectrometer include an input and output focal distance 101.6 mm, spatial groove density grating 2400 line/mm, slit width 5µm, and optical resolution of width ~0.025 nm.
A confocal quartz lens, with a diameter of 50 mm and focal length of 75 mm, was used to focus the emitted light of the plasma onto the optical fiber probe, which was placed on the focal point of the lens. This lens was focused on the highest intensity positions between electrodes (middle of radial positions, r = 15 mm). This configuration allowed the determination of species present in the plasma and enabled measurement of the plasma parameters, such as electron density, and the rotational, vibrational, and excitation temperatures in this zone.
Gas chromatography was used to study the decomposition of CO2 and the formation of CO and O2 compounds. The feed and exhaust gases were analyzed using a compact-gas chromatograph (CGC) type GC, Agilent 6890 N, equipped with a flame ionization detector (FID) and the packed GC columns Molecular Sieve 139 (MS-139) and HayeSep type Q and N. The FID can evaluate hydrocarbons such as propane, acetylene, ethylene, ethane, and others. Moreover, a thermal detector connected by columns, was used to analyze the gas components such as CO2, CO, O2, etc.

2.2. Two-Dimensional Fluid Model

2.2.1. Model Equations

For modeling purposes, half of the AC-PPP reactor was considered and azimuthal symmetry around the reactor axis was assumed. Thus, the spatial description of the problem was mathematically two-dimensional (with only axial and radial directions). The simulated domain was the discharge gap between the high-voltage (HV) and ground electrodes. This domain was extended into the conductive inlet/outlet pipes that can affect the electric field distribution (see Figure 3). The grid size was 4.5 µm.
The spatial and temporal macroscopic description of the gas discharge inside the reactor was determined by solving the fluid continuity equations for different species coupled with Poisson’s equation. These equations were solved using the finite element method (FEM).
The continuity equation for all the formed species inside the AC reactor is expressed as follows [14]:
n i t + · Γ i = m R i , m  
where ni is the number density, Γi expresses the flux for the species i, and Ri,m are the reaction rates between species i and species m.
For a typical reaction between species
A + B C + D  
the reaction rate depends on the density of each species, nA and nB.
R = k n A n B
with k, the reaction constant [14,15].
In this study, two different approaches were considered to obtain the reaction constants. For some reactions, the experimental data for these reaction rates were available in the literature [16]. In other cases, the reaction rate constants were calculated using the total collision cross sections in terms of the collisional energy, σ(ε), by the following relationship [17]:
k ( T ) = 1 k B T ( 8 π μ k B T ) 1 / 2 0 ε σ ( ε ) e ε / k B T d ε
The collisional cross section can be written as follows:
σ ( ε ) = 10 13 E I p [ B 1 ln ( ε I p ) + i = 1 N 1 B i + 1 ( 1 I p ε ) i ]  
where Ip is a parameter close (but not always equal) to the ionization or appearance potential for a given ionization channel (expressed in electron-Volts (eV)), ε is the collision energy (eV), Bi (i = 1,…, N) are fitting coefficients, and N is determined from those conditions for achieving a standard deviation of the fit from data smaller than 3–5% [18].
In the model, only the reactions that showed the greatest influence on the plasma CO2 kinetics were considered [19,20,21,22,23,24]. These reactions were based mainly on the chemistry of Wang [20] and Kozak [23] who provided useful results for gliding arc and microwave plasma, respectively. Table 1 shows the selected reactions inside the AC-PPP reactor. A total of 16 electron collisions and six atomic processes were considered.
The excited CO2 states were simplified according to the scheme explained in a previous paper from our group [25]. The levels Va, …, and Vd of Table 1 represent the two non-asymmetric excited vibrational modes; V1 is the first level of the asymmetric vibrational mode, and the excited electronic states of the CO2 are grouped in the effective level CO2*.
The flux term in the continuity fluid equations for all species (Equation (1)) was based on the momentum conservation of each species. The corresponding flux term for each species i depends on its electrical mobility, μi, and spatial diffusivity, Di, i.e.,
Γ i = ± μ i n i E D i n i  
The plus or minus sign in this equation corresponds to the sign of the charged particles [1,26].
Only the mobility coefficients for ions and electrons were included [27]. The mobility of ions was calculated according to the Langevin equation:
μ i = 0.514 m i 1 / 2 T g P t o t α i 1 / 2  
where αi is the polarization of background gas per unit of cubic angstroms; its value for various gases is presented in the existing literature on gaseous discharges [27]. In this work, the mobilities for CO2+ and C+ species were 0.0012 and 0.0009 m2/Vs, respectively.
The diffusion coefficient of the electrons and ions were instead calculated from the Einstein relation:
D e ( i ) = k B T e ( i ) q e ( i ) μ e ( i )  
with Te(i) and qe(i) being the temperature and charge of electrons and ions [28].
For neutral species, the diffusion coefficients were calculated using the distribution coefficients of Lennard–Jones [29].
The rate of change of the electron energy density is described by [1]:
ω e t + · Γ ω + e E · Γ e = R ω  
where ωe is the electron energy density, R ω is the energy loss or gain due to inelastic collisions, the term e E · Γ e accounts for the ohmic or joule heating of the electrons in the electric field, and Γ ω is the electron flux energy, that is described by:
Γ ω = 5 3 ( μ e ω e E D ε ω e )  
The electron energy loss or gain R ω is obtained by summing the collisional energy loss or gain over all reactions [14]:
R ω = j = 1 P x j k j N n n e Δ ε j  
where xj is the mole fraction of the target species for reaction j, kj is the rate coefficient for reaction j, Nn is the total neutral number density and Δ ε j is the energy loss from reaction j.
The electron energy density ωe, the mean electron energy ε, and the electron temperature Te are correlated with each other through [30]:
ω e = n e ε = 3 2 k B n e T e  
For non-electron species, the following equation was solved for the mass fraction of each species [30]:
ρ ξ k t + ρ ( u · ) ξ k = · j k + R k  
where j k is the diffusive flux vector, R k is the rate expression for species k, u is the mass averaged fluid velocity vector, ρ denotes the density of the mixture and ξ k is the mass fraction of the kth species.
The diffusive flux vector is defined as [30]:
j k = ρ ξ k V k  
with V k , being the multicomponent diffusion velocity for species k.
To initiate discharge in the reactor, electric potential should be applied between the electrodes, thus Poisson’s equation must also be considered in the model [14]:
ε 0 2 φ = σ  
where φ is the electric potential, ε 0 is the vacuum permittivity and σ is the charge density, that can be written in terms of density of the charged species, nk, and their charge, eZk [31]:
σ = e ( k = 1 k Z k n k n e )  
In this work, 16 different neutral and ionized species were considered in the model (Table 2). Thus, 16 continuity equations together with Poisson’s equation were solved with the employment of a stabilized FEM.

2.2.2. Boundary Conditions

To obtain a unique solution for the system of coupled equations with the geometry presented in Figure 3, the boundary conditions (Dirichlet and Neumann boundary conditions) must be imposed. The boundary conditions applied for the AC plasma reactor corresponded to those found in the existing literature [32]. The following boundary condition was used to account for the particle flux in walls:
Γ i · n = ( a i   s g n ( q i ) μ i E · n + 1 4 v t h ( i ) ) n i
where n is the normal vector pointing toward the tube wall and, v t h ( i ) is the thermal velocity of particles [32]:
v t h ( i ) = 8 k B π T i m i
and the number ai is defined by:
a i = { 1   s g n ( q i ) μ i E · n 0 0   s g n ( q i ) μ i E · n 0  
For electrons, as a special case, the particle flux due to secondary electron emission (SEE) was added to the system and is defined as follows [32]:
Γ e · n = ( a e μ e E · n + 1 4 v t h , e ) n e p γ p Γ p · n
where γ p are the SEE coefficients, which defines the average number of electrons emitted per impact of ions p on the tube wall.
Similarly, a boundary condition for electron energy was [32]:
Γ ω · n = 5 6 v t h , e ω e p γ p ε p ( Γ p · n )
Here, the second term is the SEE energy flux, being ε p the mean energy of the secondary electrons.
The discharge was driven by a sinusoidal electric potential applied to an electrode and the other electrode was grounded. Then the boundary condition of electric potential in the grounded electrode was:
φ = 0
and the electric potential in the powered electrode was given by:
φ = V 0 s i n ( 2 π f t )
with f = 50 Hz and V0 = 22 kV.

3. Results and Discussion

3.1. Experimental Results

3.1.1. Gas Chromatography

The inlet and outlet gases were analyzed with the gas chromatograph GC, Agilent 6890 N. The areas of chromatogram peaks are proportional to the concentrations of compounds in these gases. Then, the CO2 conversion factor can be calculated from the areas of CO2 peak, [ C O 2 ] i n   and   [ C O 2 ] o u t by [33]:
X C O 2 ( % ) = [ C O 2 ] i n [ C O 2 ] o u t [ C O 2 ] i n × 100
In a similar way, the selectivity of CO and O2 species were calculated by [32]:
S C O ( % ) = 0.5   [ C O ] o u t [ C O 2 ] i n [ C O 2 ] o u t × 100 S O 2 ( % ) = [ O 2 ] o u t [ C O 2 ] i n [ C O 2 ] o u t × 100
Figure 4a shows the conversion factor of CO2 and the selectivity of CO and O2 obtained by gas chromatography for an applied voltage between 10 and 22 kV. A maximum value of the conversion factor of 25% was found at 22 kV voltage. This factor was also compared with that obtained by the 2-D model (see Section 3.2).
The energy efficiency of this technique is given by [32]:
η ( % ) = X C O 2 ( % )   H R S E I * V m
where Vm is the molar volume, that can be considered 24.5 L/mol, H R is the reaction enthalpy for CO2 splitting ( CO 2 CO + 1 2 O 2 ) , i.e., 279.8 kJ/mol, and SEI is the specific energy input, that is calculated from discharge power and gas flow rate:
S E I = P o w e r F l o w  
The power consumed in discharge was estimated from measurements of the electrical current and HV voltage between the electrodes of the AC-PPP reactor. Then the energy efficiency could be calculated. Table 3 shows the values of this efficiency for different applied voltages. In this case, the maximum efficiency was 6.1% for an applied voltage of 22 kV. This experimental energy efficiency of CO2 conversion is plotted in Figure 4b.
In the review paper by Snoeckx et al. [7], a comparison of the conversion factor and the energy efficiency of different plasmas technologies for CO2 remediation was performed. Typical plasma technologies used for this purpose are dielectric barrier discharges (DBDs), microwave discharges (MWs), and gliding arc discharges (GAs). These systems are compared with the proposed AC-PPP reactor in Table 4.
In general, the DBD plasma has a high conversion factor, but it is not energy efficient. Its maximum energy efficiency values were approximately 6%. The MW is the technology with the highest CO2 conversion with systems that have up to 80% energy efficiency. Some systems can also achieve a high energy efficiency of 40%, but the typical values for this technology are less than 10%. The power consumption of this system is usually very high. The worst CO2 conversion is obtained for the GA discharge. The reactor introduced in this study had a conversion factor value greater than GA and close to the maximum for DBD. Energy efficiency is low compared to the other technologies, but it has the advantages of its simple design and low power consumption, which are important factors for its industrial implementation.

3.1.2. Optical Emission Spectroscopy

In the plasma discharge media, the collision process excites several plasma species to upper states decaying and emitting photons with various specific wavelengths. By recording the emission spectra of the plasma discharge medium, these photons can be detected and analyzed.
In this work, the OES technique was used to identify the different reactive species produced between the electrodes of the CO2 plasma reactor. Vibro-rotational bands of CO2, CO2+, C2, and CO radicals and atomic lines of oxygen were detected (see Figure 5 and Table 5). Vibro-rotational bands of radical ·OH and the Hβ emission line of hydrogen, coming from small water impurities in the inlet gas, were also identified. The amount of this water was so small that it had no relevant influence on the plasma.
A study of the dependencies of the normalized intensities of the headband of vibro-rotational bands of these formed species and the oxygen atomic line 771.1 nm on the plasma position was carried out. Figure 6 plots these dependencies for an applied voltage of 22 kV and a frequency of 50 Hz. The intensity of CO 2 species was higher than that of other species although it was reduced close to the electrodes. Furthermore, the intensities of CO , O, OH, and C 2 species increased near the electrodes but they had a minimum value at the middle of discharge. The relative intensities of Figure 6 are not directly related to the populations of these species because these are affected by many parameters, such as electron density, the cross sections of particles, etc.; however, it is indicative of the presence of species between the reactor electrodes.
The density and temperature of particles in the plasma are the most fundamental parameters in the physical characteristics of gas discharges and play a very important role in understanding the discharge physics and optimization of the operating parameters of the plasma reactor. From an experimental point of view, these parameters are commonly measured using different techniques including Langmuir probes, microwave interferometry, laser Thomson scattering, and optical emission spectroscopy (OES) techniques [34].

Rotational, Vibrational, and Excitation Temperatures

As indicated, the OH emission band ( A 2 Σ X 2 Π ) was detected in the wavelength ranges of 306–312 nm. The formation of the ·OH radical is generally caused by water vapor molecules that are present as impurities in the gas [35,36]. The rotational, vibrational, and electronic states of this radical are populated by inelastic collisions with particles within the plasma (heavy particles and electrons) which produce variations in their energies.
The rotational temperature gives the population distribution of rotational states. For these states, which have a small separation in energy, the effect of collisions with heavy particles are predominant; the population of rotational states is given by the energy of these particles. In equilibrium conditions, the rotational temperature is considered a good approximation of the gas temperature (mean kinetic temperature of heavy particles).
On the other hand, the population of vibrational and electronic states, with higher energy separation, is given mainly by collisions with electrons. In this case, the obtained vibrational and excitation temperatures can be considered a closer approximation to the electron temperature.
In this work, the emission spectrum of the OH A 2 Σ X 2 Π band was used for the determination of the rotational, vibrational, and excitation temperatures in the AC plasma reactor. Figure 7a shows an example of this spectrum in the AC plasma reactor at the applied AC voltage of 22 kV. Rotational, vibrational, and excitation temperatures were calculated using SPECAIR software that fits a simulated spectrum to experimental data to estimate these temperatures (see Figure 7a) [26]. For this simulation work, all the factors affecting the line shape, such as the instrumental resolution or the collisional broadenings, were considered.
By applying this fitting to vibro-rotational bands, the rotational, vibrational, and excitation temperatures were obtained in terms of position with an error of 7% (Figure 7b). For the three temperatures, their values were constant along the reactor, due to the parameters oscillating between the electrodes during the cycle of AC voltage (see next section). These results correspond to time average values during this cycle.
Figure 7b shows that the experimental rotational temperature was about 2000 K for all positions. The vibrational and excitation temperatures were about 5000 K and 18,000 K, respectively, which means the plasma was in 2-T conditions, where the electron temperature was higher than the gas temperature. The energy of the heavy particles and electrons were enough to produce the conversion of the CO2 molecules.

Electron Number Density

To find out whether the electron collisions are the main cause of molecule dissociation in the formed discharges, the electron number density was experimentally calculated in the plasma positions focused on by the lens.
The Stark broadening analysis of the spectral profile of the Hβ emission line (486.1 nm) is the most usual procedure for the experimental determination of electron density in a plasma discharge [37]. The Stark broadening of this line depends on electron density according to the expression [38]:
Δ λ s t a r k = 2 × 10 11 ( n e 2 / 3 )
where density is in cm−3 and Stark broadening is in nm.
The pressure broadening occurs when the energy states of the emitting species are disturbed by the neutral species in the plasma discharge. This broadening depends on resonance and van der Waals effects. In this experiment, the hydrogen atom density was very low, and the resonance effect can be neglected. Thus, the van der Waals broadening was the only contribution to the pressure broadening, which can be described as follows [39,40,41]:
Δ λ P r e s s u r e = Δ λ v a n   d e r   W a a l s = 8.18 × 10 26   λ 2 ( α R 2 ) 2 / 5 ( T g μ ) 3 / 10 P k T g
where µ is the reduced mass of the colliding particles (µ= 0.92308 for H/CO2 collision), λ is the wavelength of the Hβ line (486.1 nm), α is the molecular polarizability of the CO2 disturbing particles (α= 2.46 × 10−24 cm3), 〈R2〉 is the difference of the squared radius of the upper and lower levels of Hβ transition, Tg is the gas temperature in K, and P is the pressure (1 atm for atmospheric pressure).
The Stark and pressure broadenings are the main contributions to the Lorentz shape of a line. So, the full-width at half-maximum (FWHM) value of the Lorentz profile can be obtained from these broadenings [42]:
λ L o r e n t z = λ S t a r k + λ P r e s s u r e  
Another contribution to the line broadening is the Doppler effect through particle movement. This broadening has a Gaussian profile and can be written as follows [42]:
λ D ( n m ) = 7.2 × 10 7 T g M λ
where λ is the wavelength of β hydrogen line (486.1 nm), M is the molar weight of hydrogen, which is equal to 1, and Tg is the gas temperature in Kelvin units.
The total Gaussian profile is generally considered as a convolution of Doppler and instrumental profiles, then [43]:
λ G = λ D 2 + λ I 2
where ΔλG and ΔλI are the Gaussian and instrumental broadenings, respectively. The optical components resolution directly affects the instrumental broadening.
A convolution of these Gaussian and Lorentzian profiles results in a Voigt shape. Fitting of experimental data to this analytical curve allowed us to separate Lorentz and Gaussian contributions (Figure 8a). In this study, Origin software was used for this purpose. In this way, the FWHM of the Lorentzian broadening could be distinguished.
The Stark broadening and, consequently, the plasma electron density could be calculated via Equations (28) and (30), respectively.
Figure 8b shows the experimental values of the electron density obtained in this way. A constant distribution along the discharge was again obtained. The experimental value of electron density was of the order of 5 × 1014 cm−3. This high value supports the electron kinetics for CO2 decomposition assumed for this reactor.

3.2. Model Results

The considered continuity model for such a plasma reactor contributes to understanding the chemical and dynamical behaviors of the charged particles inside the reactor. Using the FEM, the equation system described in Section 2.2.1 was solved. This system corresponds to the continuity equation for each species considered in the model (Equation (1)), the energy equation for electrons (Equation (9)), and Poisson’s equation for the electric field (Equation (15)). For studying the non-electron species (neutral species), the mass fraction equation should be solved (Equation (13)). Therefore, the temporal and spatial behaviors of these species inside the CO2 plasma reactor were calculated. The simulation parameters used in these calculations are presented in Table 6.
To show the time evolution of the discharge, the results of the plasma simulation for five different times during the AC cycle are presented in the following subsections. Figure 9 plots the selected time against the applied HV voltage.

3.2.1. Electric Potential

Figure 10 shows variations of the electric potential as a function of position at times 2, 6, 10, 14, and 18 ms. The electric potential presents the highest gradient near the electrodes, which is due to electron accumulation at the anode shielding the electric field. This sheath formation has also been found in other contact-electrode plasmas, such as DBD [44,45,46]. Higher electric potential gradients at the electrodes gave a higher electric field in these positions. The electric field was able to reach a value of 4000 kV/m for 22 kV near the electrodes, decreasing to 1000 kV/m in the zone between them. Inside the inlet and outlet metal pipes, the potential was almost constant, which means the electric field was zero in these zones.
When the applied HV voltage was changed, the potential oscillated between electrodes according to the voltage variation, maintaining similar spatial distribution behavior (see different times in Figure 10).

3.2.2. Electron Density

After applying the electric field, ionization of species was started inside the AC-PPP reactor, producing electrons, which are responsible for most of the mechanisms of the CO2 splitting. These electrons start moving to the electrodes, where they concentrate, under the effect of the electric field. During their movement, they can collide with other species present in the reactor, producing new reactions, such as ionization, dissociation, etc. The determination of the spatial distribution of electrons is essential to find the kinetic energy of the discharge.
The electron density inside the plasma reactor for times 2, 6, 10, 14, and 18 ms are shown in Figure 11. During this time, when a positive voltage was applied, the electrons were concentrated in the powered electrode. They moved inside the plasma reactor to the surface electrode, where they were lost by recombination with its positive charge. When the voltage applied on the electrodes changed, these electrons moved to the ground electrode.
In the model results, high values of electron density inside the inlet/outlet pipes were observed, which means that these pipes contained high levels of chemical activity even though the electric field was low.
The experimental results were time-averaged values over the measuring period. Then, to compare with the experimental results, the time averages of the simulated electron density over a period (T = 20 ms) were performed for the radial position r = 15 mm, where the lens was focused. The results of these calculations are plotted in Figure 8b. A close agreement between the experimental results and this averaged result was found, which supports the soundness of the model.

3.2.3. Electron Temperature

In plasma dominated by electron interactions, its chemistry depends on the kinetic energy of these particles. The electrons take the energy from the electric field and share their energy with other particles by collisions. The electron temperature is a measure of the kinetic energy of the electrons.
Figure 12 plots the dependence of the electron temperature on the positions for times 2, 6, 10, 14, and 18 ms. The highest values of the electron temperature corresponded to positions where this field was highest, i.e., the powered electrode for positive voltage and the ground electrode for negative voltage. Otherwise, the electron temperature in the inlet/outlet pipes was lower by the low electric field in these zones.
The model results of electron temperature have also been averaged over the period (T = 20 ms). The averaged values for r = 15 mm are compared with the experimental measures of rotational, vibrational, and excitation temperatures in Figure 7b. The values of the theoretical electron temperatures were higher than these temperatures, which is typical of a 2-temperature plasma, characterized by the following relation between the temperatures: T e T e x c > T v > T r ~   T g [47].

3.2.4. CO2 Splitting

The high values of electron density and temperature indicate that the chemistry in the AC-PPP reactor was dominated by the collisions of electrons with CO2 ground state molecules. As predicted by the model, the most relevant mechanisms of CO2 splitting were:
-
the electron impact dissociation (or dissociative de-excitation collision) producing CO and O (see reactions 7 and 9 in Table 1),
-
the electron ionization producing CO2+ ions (see reaction 12 in Table 1),
-
the electron dissociative attachment producing CO molecules and O ions (see reaction 15 in Table 1).
These mechanisms align with those proposed by Bogaerts et al. [48] for pure CO2 splitting. The energy required by these three processes is higher than that needed for splitting by C=O bound breaking ( CO 2 CO + 1 2 O 2 ) , i.e., 279.8 kJ/mol. This explains the low energy efficiency found for this plasma, i.e., lower than 6%.
The distribution of the four formed species: CO2+, CO, O, O2, and O were studied to identify the reactions that dominated in each part of the reactor.

Density of CO2+ Ion

Figure 13 shows the spatial distribution of the density of CO2+ ions formed for times, 2, 6, 10, 14, and 18 ms. In contrast with the electrons which moved toward the power electrode during the positive voltage applied (2, 6, and 10 ms) and were aggregated there, the electric field accelerated the produced positive ions CO2+ in the direction of this field. If a CO2+ ion obtains enough acceleration from the electric field, it collides with the background gas molecules and ionizes them. However, this phenomenon rarely happens in high electron density plasmas, where CO2+ ions are less effective than the electrons in the impacting ionization process. As shown in Figure 13, the CO2+ ions accumulated near the ground electrode during the positive voltage applied, and they accumulated in the powered electrode with the negative voltages.
The maximum values of the CO2+ density were found inside the metallic inlet/outlet pipes. The positive ions accumulated in this zone because they did not move toward the electrodes due to the low electric field inside the pipes. However, the CO2+ ion is an unstable species and will be quickly recombined with electrons or negative ions O2 to convert to CO and O, and/or O2 [48].

Density of Neutral Species CO, O, and O2

The densities of the neutral species CO and O generated inside the plasma discharge medium of the developed AC-PPP reactor at times 2, 6, 10, 14, and 18 ms are presented in Figure 14 and Figure 15. These neutral species were produced near the powered electrode during the positive voltage and the ground electrode during the negative voltage. This was due to the high electron density near the electrodes, which produced the dissociation and attachment processes of the neutral species in these regions.
The maximum values of CO and O densities are equal to 4.7 × 1011 cm−3 and 2.3 × 1011 cm−3, respectively. The maximum CO density is higher than O because this molecule is produced by two mechanisms, electron impact dissociation (reaction 7 and 9) and electron dissociative attachment (reaction 15). CO is a relatively stable molecule that can recombine with O and O to form new CO2 when the residence time is high [48]. However, the O atoms recombine quickly into oxygen molecules, O2 and/or ozone, O3.
Oxygen molecules, O2, are formed mainly by the reaction of oxygen atoms with the CO2 gas through neutral dissociation and 3-body recombination, reactions 17 and 20, respectively. Figure 16 shows the model results of the spatial distribution of the O2 density at times 2, 6, 10, 14, and 18 ms. The maximum values of this density correspond to positions where O and CO2 densities are higher.

Density of O Ion

The electrons traveling to the electrodes can produce negative ions by attachment with the species present in the plasma. This was the case for the negative ion O (reaction 15).
The spatial distribution of the density of O ions produced by this process for times 2, 6, 10, 14, and 18 ms are shown in Figure 17.
The negative O ions moved toward the electrodes and were aggregated there according to electron density. This negative oxygen ion density increased quickly and became constant after the first millisecond.
The electron dissociative attachment that produced the O ions was higher in the inlet/outlet pipes. The high densities and low temperatures in these positions favored the electron attachment in these zones.

Density of CO2 Molecules

Reactions 7, 9, and 15 mainly describe the CO2 conversion and the density distribution of this molecule in the reactor. The spatial distribution of CO2 density between the two electrodes inside the AC-PPP reactor is presented in Figure 17. The CO2 density decreased near the electrodes due to the high values of electron density and temperature in this position. A zoom of the zone near the electrodes was plotted in Figure 18. Due to the plasma is still active inside the conductor inlet/outlet pipes, the CO2 was also discomposed inside these pipes.
As is shown, during the positive voltage (2, 6, and 10 ms), the maximum CO2 decomposition was near the powered electrode, where the electron density was highest. Conversely, the maximum value was near the ground electrode during the negative voltage (14 and 18 ms). A theoretical estimation of the CO2 conversion factor was made with the time averages of the CO2 density over a period, T = 20 ms, in the outlet tube position, which could be compared with the experimental results of gas chromatography. A good agreement between experimental and theoretical results was found (see Figure 4a).

4. Conclusions

In this study, a new AC-PPP reactor for CO2 remediation was presented. We undertook experimental measurements by gas chromatography and optical emission spectroscopy to characterize and determine its capacity for CO2 decomposition. The CO2 conversion and energy efficiency were obtained by the GC method and compared with other plasma sources. It was found that CO2 conversion in the developed plasma reactor was close to 25%, higher than GA, and similar to DBD. Although the energy efficiency of the AC-PPP reactor (near 6.1%) is lower than the maximum values of the other technologies, this AC-PPP reactor could be very useful and suitable for industrial applications due to its simple design and capacity to work with AC at atmospheric pressure.
The experimental electron density and the rotational, vibrational, and excitation temperatures were calculated by OES spectroscopy. The high values of vibration and excitation temperatures, around 5100 K and 18,300 K, respectively, compared with the rotational temperature (~2000 K), show that it is a 2-temperature plasma, where the electrons have higher kinetic energy and are mainly responsible for the splitting reactions. This result was confirmed by the high electron number density found experimentally in this discharge with a maximum value of 5 × 1014 cm−3.
A 2-D fluid model was developed to simulate the kinetics of CO2 decomposition in this plasma. The temporal and spatial behaviors of the densities of CO2 and the formed species inside the plasma reactor were calculated according to this model. According to the theoretical results, the charged species oscillated between the electrodes at each cycle of the AC potential with a frequency of 50 Hz. This movement can cause more ionization in the plasma reactor and a better CO2 decomposition by increasing collisions.
The model has demonstrated that the main mechanisms of CO2 splitting were the electron impact dissociation, the electron ionization, and the electron dissociative attachment. The species O and CO2+ ions, CO and O2 molecules, and O atoms were formed by these mechanisms. The spatial distributions of these species were calculated. The results of the model showed that the discharge was also formed inside the inlet and outlet pipes. This extension of the discharge influences the improvement of CO2 conversion by this reactor.
The distribution of the CO2 density was also determined. A final reduction of this molecule of about 22% could be achieved after the passing of gas along the reactor. This CO2 conversion agreed with experimental values of this parameter obtained by GC.
The electron density and temperature averaged over a period was calculated. The value and behavior of this electron density under study conditions also agreed with the experimental results.
We consider that the developed 2-D model has been shown to be valid for describing and understanding the main processes that occur in the reactor and which lead to the decomposition of CO2 and its transformation into other species. This type of reactor could contribute to the optimization of working conditions and improvements in reactor design.

Author Contributions

A.B. designed the research; A.B., S.I.M., A.F. and S.K. performed the experiments and the numerical calculations; A.B. and A.R. supervised the research (experiments and calculations); A.B., A.R. and D.A.K. were involved in writing and editing the manuscript; A.R. and D.A.K. were involved in funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

The research and APC was funded by the project BUT InterAcademic Partnerships (PPI/APM/2018/1/00033/U/001) of The Polish National Agency for Academic Exchange.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors are grateful for the part- financing of the study by the Andalusian Regional Government (Research Groups FQM-136). This publication was financed within the project BUT InterAcademic Partnerships (PPI/APM/2018/1/00033/U/001) of The Polish National Agency for Academic Exchange. The study was conducted in the Scientific Cooperation Agreement, “The possibility of the renewable energy sources usage in the context of improving energy efficiency and air quality in the buildings and civil constructions” between BUT and UCO supported by science research funds from the Bialystok University of Technology WZ/WB-IIS/9/2019.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liebermanm, M.A.; Lichtenberg, A.J. Principles of Plasma Discharges and Materials Processing, 2nd ed.; John Wiley & Sons: New York, NY, USA, 2005; pp. 122–250. [Google Scholar]
  2. Kogelschatz, U. Dielectric-barrier discharges: Their history, discharge physics, and industrial applications. Plasma Chem. Plasma Process. 2003, 23, 1–46. [Google Scholar] [CrossRef]
  3. Aresta, M. Carbon Dioxide as Chemical Feedstock, 1st ed.; John Wiley & Sons: New York, NY, USA, 2010; pp. 152–180. [Google Scholar]
  4. Kraus, M.; Eliasson, B.; Kogelschatz, U.; Wokaun, A. CO2 reforming of methane by the combination of dielectric-barrier discharges and catalysis. Phys. Chem. Chem. Phys. 2001, 3, 294–300. [Google Scholar] [CrossRef]
  5. Fridman, A. Plasma Chemistry, 1st ed.; Cambridge University Press: Cambridge, UK, 2008; pp. 231–254. [Google Scholar]
  6. Trencheva, G.; Nikiforov, A.; Wanga, W.; Kolevc, S.; Bogaerts, A. Atmospheric pressure glow discharge for CO2 conversion: Model-based exploration of the optimum reactor configuration. Chem. Eng. J. 2019, 362, 830–841. [Google Scholar] [CrossRef]
  7. Snoeckx, R.; Bogaerts, A. Plasma technology—A novel solution for CO2 conversion? Chem. Soc. Rev. 2017, 46, 5805–5863. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Li, D.; Li, X.; Bai, M.; Tao, X.; Shang, S.; Dai, X.; Yin, Y. CO2 Reforming of CH4 by Atmospheric Pressure Glow Discharge Plasma: A High Conversion Ability. Int. J. Hydrogen Energy 2009, 34, 308–313. [Google Scholar] [CrossRef]
  9. Li, H.; Zhou, Y.; Donnelly, V.M. Optical and Mass Spectrometric Measurements of the CH4–CO2 Dry Reforming Process in a Low Pressure, Very High Density, and Purely Inductive Plasma. J. Phys. Chem. A 2020, 124, 7271–7282. [Google Scholar] [CrossRef] [PubMed]
  10. Li, M.; Xu, G.; Tian, Y.; Chen, L.; Fu, H. Carbon Dioxide Reforming of Methane Using DC Corona Discharge Plasma Reaction. J. Phys. Chem. A 2004, 108, 1687–1693. [Google Scholar] [CrossRef]
  11. Yu, Q.; Kong, M.; Liu, T.; Fei, J.; Zheng, X. Characteristics of the decomposition of CO2 in a dielectric packed-bed plasma reactor. Plasma Chem. Plasma Process. 2012, 32, 153–163. [Google Scholar] [CrossRef]
  12. Wang, Q.; Shi, H.; Yan, B.; Jin, Y.; Cheng, Y. Steam enhanced carbon dioxide reforming of methane in DBD plasma reactor. Int. J. Hydrogen Energy 2011, 36, 8301–8306. [Google Scholar] [CrossRef]
  13. Liu, C.J.; Xu, G.H.; Wang, T. Non-thermal plasma approaches in CO2 utilization. Fuel Process. Technol. 1999, 58, 119–134. [Google Scholar] [CrossRef]
  14. Shrestha, R.; Tyata, R.; Subedi, D. Estimation of electron temperature in atmospheric pressure dielectric barrier discharge using line intensity ratio method. Kathmandu Univ. J. Sci. Eng. Technol. 2013, 8, 37–42. [Google Scholar] [CrossRef] [Green Version]
  15. Raju, G.G. Gaseous Electronics Theory and Practice, 1st ed.; CRC Press Taylor & Francis Group: New York, NY, USA, 2006; pp. 351–423. [Google Scholar]
  16. Torres, J.; Van De Sande, M.J.; Van Der Mullen, J.; Gamero, A.; Sola, A. Stark broadening for simultaneous diagnostics of the electron density and temperature in atmospheric microwave discharges. Spectrochim. Acta B At. Spectrosc. 2006, 61, 58–68. [Google Scholar] [CrossRef]
  17. Laux, C.O.; Spence, T.G.; Kruger, C.H.; Zare, R.N. Optical diagnostics of atmospheric pressure air plasmas. Plasma Sources Sci. Technol. 2003, 12, 125. [Google Scholar] [CrossRef]
  18. Dorai, R. Modeling of Plasma Remediation of NOx Using Global Kinetic Models Accounting for Hydrocarbons. Master’s Thesis, University of Illinois at Urbana-Champaign, Champaign, IL, USA, 2000. [Google Scholar]
  19. Hayashi Database. Available online: www.lxcat.net/Hayashi (accessed on 19 March 2020).
  20. Wang, W.; Berthelot, A.; Kolev, S.; Tu, X.; Bogaerts, A. CO2 conversion in a gliding arc plasma: 1D cylindrical discharge model. Plasma Sources Sci. Technol. 2016, 25, 065008. [Google Scholar] [CrossRef]
  21. PLASIMO Software. Plasma Matters. Eindhoven University of Technology. Available online: http://plasimo.phys.tue.nl/ (accessed on 19 March 2020).
  22. Moravej, M.; Yang, X.; Hicks, R.F.; Penelon, J.; Babayan, S.E. A radio-frequency nonequilibrium atmospheric pressure plasma operating with argon and oxygen. J. Appl. Phys. 2006, 99, 093305. [Google Scholar] [CrossRef]
  23. Kozák, T.; Bogaerts, A. Splitting of CO2 by vibrational excitation in non-equilibrium plasmas: A reaction kinetics model. Plasma Sources Sci. Technol. 2014, 23, 045004. [Google Scholar] [CrossRef]
  24. Mikoviny, T.; Kocan, M.; Matejcik, S.; Mason, N.J.; Skalny, J.D. Experimental study of negative corona discharge in pure carbon dioxide and its mixtures with oxygen. J. Phys. Appl. Phys. 2004, 37, 64–73. [Google Scholar] [CrossRef]
  25. Barkhordari, A.; Mirzaei, S.I.; Falahat, A.; Rodero, A. Numerical and experimental study of an Ar/CO2 plasma in a point-to-plane reactor at atmospheric pressure. Spectrochim. Acta Part B At. Spectrosc. 2021, 177, 106048. [Google Scholar]
  26. Chiper, A.S.; Chen, W.; Mejlholm, O.; Dalgaard, P.; Stamate, E. Atmospheric pressure plasma produced inside a closed package by a dielectric barrier discharge in Ar/CO2 for bacterial inactivation of biological samples. Plasma Sources Sci. Technol. 2011, 20, 025008. [Google Scholar] [CrossRef]
  27. Prigogine, I. Chemical kinetics and dynamics. Ann. N. Y. Acad. Sci. 2003, 988, 128–132. [Google Scholar] [CrossRef] [PubMed]
  28. Marrero, T.; Mason, E.A. Gaseous diffusion coefficients. J. Phys. Chem. Ref. Data 1972, 1, 3–118. [Google Scholar] [CrossRef]
  29. Jasper, A.W.; Miller, J.A. Lennard–Jones parameters for combustion and chemical kinetics modeling from full-dimensional intermolecular potentials. Combust. Flame 2014, 161, 101–110. [Google Scholar] [CrossRef]
  30. Raizer, Y.P.; Kisin, V.I.; Allen, J.E. Gas Discharge Physics, 1st ed.; Springer: Berlin/Heidelberg, Germany, 1991; pp. 283–301. [Google Scholar]
  31. Becker, M.; Loffhagen, D.; Schmidt, W. A stabilized finite element method for modeling of gas discharges. Comput. Phys. Commun. 2009, 180, 1230–1241. [Google Scholar] [CrossRef]
  32. Hagelaar, G.; De Hoog, F.; Kroesen, G. Boundary conditions in fluid models of gas discharges. Phys. Rev. E 2000, 62, 1452. [Google Scholar] [CrossRef] [PubMed]
  33. Aerts, R.; Somers, W.; Bogaerts, A. Carbon dioxide splitting in a dielectric barrier discharge plasma: A combined experimental and computational study. ChemSusChem 2015, 8, 702–716. [Google Scholar] [CrossRef] [PubMed]
  34. Reyes, P.G.; Gomez, A.; Vergara, J.; Martinez, H.; Torres, C. Plasma diagnostics of glow discharges in mixtures of CO2 with noble gases. Rev. Mex. Fis. 2017, 63, 363–371. [Google Scholar]
  35. Shimizu, K.; Oda, T. Emission spectrometry for discharge plasma diagnosis. Sci. Technol. Adv. Mater. 2001, 2, 577. [Google Scholar] [CrossRef] [Green Version]
  36. Maouhoub, E. Determination of the real experimental composition of the Ar-CO2 mixture from calculation of the atomic excitation temperature. EPJ Appl. Phys. 1998, 3, 81–89. [Google Scholar] [CrossRef]
  37. Griem, H.R. Plasma Spectroscopy, 2nd ed.; Cambridge University Press: Cambridge, UK, 2005; pp. 342–390. [Google Scholar]
  38. Gigosos, M.A.; Cardenoso, V. New plasma diagnosis tables of hydrogen Stark broadening including ion dynamics. J. Physical B-At. Mol. Opt. 1996, 29, 4795. [Google Scholar] [CrossRef]
  39. Belostotskiy, S.G.; Ouk, T.; Donnelly, V.M.; Economou, D.J.; Sadeghi, N. Gas temperature and electron density profiles in an argon dc microdischarge measured by optical emission spectroscopy. J. Phys. D Appl. Phys. 2010, 107, 053305. [Google Scholar] [CrossRef]
  40. Nikiforov, A.Y.; Leys, C.; Gonzalez, M.A.; Walsh, J.L. Electron density measurement in atmospheric pressure plasma jets: Stark broadening of hydrogenated and non-hydrogenated lines. Plasma Sources Sci. Technol. 2015, 24, 034001. [Google Scholar] [CrossRef]
  41. Boffard, J.B.; Lin, C.C.; DeJoseph, C.A., Jr. Application of excitation cross sections to optical plasma diagnostics. J. Phys. D Appl. Phys. 2004, 37, R143. [Google Scholar] [CrossRef]
  42. Mehr, F.J.; Biondi, M.A. Electron-temperature dependence of electron-ion recombination in Argon. Phys. Rev. 1968, 176, 322. [Google Scholar] [CrossRef]
  43. Pandhija, S.; Rai, A.K. In situ multielemental monitoring in coral skeleton by CF-LIBS. Appl. Physical B 2009, 94, 545–552. [Google Scholar] [CrossRef]
  44. Braun, D.; Gibalov, V.; Pietsch, G. Two-dimensional modelling of the dielectric barrier discharge in air. Plasma Sources Sci. Technol. 1992, 1, 166. [Google Scholar] [CrossRef]
  45. Moreno Wandurraga, S.H. Reduced Reaction Kinetics Model for CO2 Dissociation in Non-Thermal Microwave Discharges. Master’s Thesis, Delft University of Technology, Delft, The Netherlands, 2015. [Google Scholar]
  46. Khodja, K.; Belasri, A. Sheath formation study in Ne-Xe DBD discharge. Radiat. Eff. Defects Solids 2012, 167, 734–742. [Google Scholar] [CrossRef]
  47. Staack, D.; Farouk, B.; Gutsol, A.; Fridman, A. Characterization of a dc atmospheric pressure normal glow discharge. Plasma Sources Sci. Technol. 2005, 14, 700. [Google Scholar] [CrossRef]
  48. Bogaerts, A.; Snoeckx, R.; Trenchev, G.; Wang, W. Plasma Chemistry and Gas Conversion; BoD–Books on Deman: Norderstedt, Germany, 2018. [Google Scholar]
Figure 1. Experimental arrangement.
Figure 1. Experimental arrangement.
Applsci 11 10047 g001
Figure 2. Photographic image and the 3D scheme of the AC-PPP reactor: (1) high voltage dielectric, (2) Pyrex tube, (3) copper disk, (4) copper pipe.
Figure 2. Photographic image and the 3D scheme of the AC-PPP reactor: (1) high voltage dielectric, (2) Pyrex tube, (3) copper disk, (4) copper pipe.
Applsci 11 10047 g002
Figure 3. The simulated domain for the AC-PPP reactor in the 2-D model.
Figure 3. The simulated domain for the AC-PPP reactor in the 2-D model.
Applsci 11 10047 g003
Figure 4. (a) Dependence of CO2 conversion and CO/O2 selectivity on applied AC voltage. (b) Experimental dependence of energy efficiency on applied AC voltage.
Figure 4. (a) Dependence of CO2 conversion and CO/O2 selectivity on applied AC voltage. (b) Experimental dependence of energy efficiency on applied AC voltage.
Applsci 11 10047 g004
Figure 5. Emission spectrum of the AC-PPP at the AC voltage of 22 kV and the position z = 5 mm.
Figure 5. Emission spectrum of the AC-PPP at the AC voltage of 22 kV and the position z = 5 mm.
Applsci 11 10047 g005
Figure 6. Relative intensities for emission lines of various species as a function of position at an AC voltage of 22 kV and 1 cm distance between electrodes.
Figure 6. Relative intensities for emission lines of various species as a function of position at an AC voltage of 22 kV and 1 cm distance between electrodes.
Applsci 11 10047 g006
Figure 7. (a) Experimental emission spectrum of the OH A 2 Σ X 2 Π band (dots) with their SPECAIR fitting (line) for the determination of the rotational, vibrational, and excitation temperatures. (b) Variations of rotational, vibrational, and excitation temperatures as a function of position at the AC voltage of 22 kV and 1 cm distance between electrodes.
Figure 7. (a) Experimental emission spectrum of the OH A 2 Σ X 2 Π band (dots) with their SPECAIR fitting (line) for the determination of the rotational, vibrational, and excitation temperatures. (b) Variations of rotational, vibrational, and excitation temperatures as a function of position at the AC voltage of 22 kV and 1 cm distance between electrodes.
Applsci 11 10047 g007
Figure 8. (a) Detected Hβ emission line with its fitting to a Voigt profile. (b) Variations of electron density as a function of position at the AC voltage of 22 kV and 1 cm distance between electrodes.
Figure 8. (a) Detected Hβ emission line with its fitting to a Voigt profile. (b) Variations of electron density as a function of position at the AC voltage of 22 kV and 1 cm distance between electrodes.
Applsci 11 10047 g008
Figure 9. Applied HV voltage with the selected times for the presentation of the model results.
Figure 9. Applied HV voltage with the selected times for the presentation of the model results.
Applsci 11 10047 g009
Figure 10. Spatial distribution of the electric potential for times (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms and (e) 18 ms for the reactor at 22 kV AC voltage and 1 cm distance between electrodes.
Figure 10. Spatial distribution of the electric potential for times (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms and (e) 18 ms for the reactor at 22 kV AC voltage and 1 cm distance between electrodes.
Applsci 11 10047 g010
Figure 11. Spatial distribution of the electron density at time (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms, and (e) 18 ms for the reactor at 22 kV AC voltage and 1 cm distance between electrodes.
Figure 11. Spatial distribution of the electron density at time (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms, and (e) 18 ms for the reactor at 22 kV AC voltage and 1 cm distance between electrodes.
Applsci 11 10047 g011
Figure 12. Spatial distribution of the electron temperature for times (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms, and (e) 18 ms for the reactor at 22 kV AC voltage and 1 cm distance between electrodes.
Figure 12. Spatial distribution of the electron temperature for times (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms, and (e) 18 ms for the reactor at 22 kV AC voltage and 1 cm distance between electrodes.
Applsci 11 10047 g012
Figure 13. Spatial distribution of the CO2+ density for times (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms, and (e) 18 ms for the reactor at 22 kV AC voltage and 1 cm distance between electrodes.
Figure 13. Spatial distribution of the CO2+ density for times (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms, and (e) 18 ms for the reactor at 22 kV AC voltage and 1 cm distance between electrodes.
Applsci 11 10047 g013
Figure 14. Spatial distribution of the CO density for times (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms, and (e) 18 ms for the reactor at 22 kV AC voltage and 1 cm distance between electrodes.
Figure 14. Spatial distribution of the CO density for times (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms, and (e) 18 ms for the reactor at 22 kV AC voltage and 1 cm distance between electrodes.
Applsci 11 10047 g014
Figure 15. Spatial distribution of the O density for times (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms, and (e) 18 ms for the reactor at AC 22 kV voltage and 1 cm distance between electrodes.
Figure 15. Spatial distribution of the O density for times (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms, and (e) 18 ms for the reactor at AC 22 kV voltage and 1 cm distance between electrodes.
Applsci 11 10047 g015
Figure 16. Spatial distribution of the O2 density for times (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms, and (e) 18 ms for the reactor at AC 22 kV voltage and 1 cm distance between electrodes.
Figure 16. Spatial distribution of the O2 density for times (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms, and (e) 18 ms for the reactor at AC 22 kV voltage and 1 cm distance between electrodes.
Applsci 11 10047 g016
Figure 17. Spatial distribution of the O- ion density for times (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms, and (e) 18 ms for the reactor at 22 kV AC voltage and 1 cm distance between electrodes.
Figure 17. Spatial distribution of the O- ion density for times (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms, and (e) 18 ms for the reactor at 22 kV AC voltage and 1 cm distance between electrodes.
Applsci 11 10047 g017
Figure 18. Spatial distribution of CO2 density at time (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms, and (e) 18 ms for the reactor at AC 22 kV voltage and 1 cm distance between electrodes.
Figure 18. Spatial distribution of CO2 density at time (a) 2 ms, (b) 6 ms, (c) 10 ms, (d) 14 ms, and (e) 18 ms for the reactor at AC 22 kV voltage and 1 cm distance between electrodes.
Applsci 11 10047 g018
Table 1. Selected reactions for different formed species in the AC-PPP reactor.
Table 1. Selected reactions for different formed species in the AC-PPP reactor.
Num.Process NameReactionsReaction Rate ConstantsRef.
Electron Collisions
1Vibrational Excitation e + C O 2 e + C O 2 ( V 1 ) k 1 [19]
2Vibrational Excitation e + C O 2 e + C O 2 ( V a ) k 2 [19]
3Vibrational Excitation e + C O 2 e + C O 2 ( V b ) k 3 [19]
4Vibrational Excitation e + C O 2 e + C O 2 ( V c ) k 4 [19]
5Vibrational Excitation e + C O 2 e + C O 2 ( V d ) k 5 [19]
6Electronic Excitation e + C O 2 e + C O 2 * k 6 [21]
7Dissociative De-excitation Collision e + C O 2 ( V x ) e + C O + O k 7 [19]
8Electron Ionization e + C O 2 2 e + O 2 + C + k 8 [19]
9Electron Dissociation e + C O 2 e + O + C O k 9 [19]
10Electronic Excitation e + C O 2 e + C O 2 * k 10 [21]
11Elastic Collision e + C O 2 e + C O 2 k 11 [21]
12Electron Ionization e + C O 2 2 e + C O 2 + k 12 [22]
13Electron Ionization e + C O 2 2 e + C O + O + k 13 [22]
14Dissociative Electron-Ion Recombination e + C O 2 + C O + O k 14 [22]
15Dissociative Attachment e + C O 2 C O + O k 15 [23]
16Electron Dissociation e + O 2 e + 2 O k 16 [23]
Atomic Collisions
17Neutral Dissociation O + C O 2 C O + O 2 k 17 [20]
18Charge Exchange C O 2 + C O 2 + C O 2 + + C O 2 k 18 [20]
19Scattering O + O 2 O + O 2 k 19 [20]
203-body Recombination O + O + C O 2 O 2 + C O 2 k 20 [23]
213-body Neutral-Neutral Collision O + C + C O 2 C O + C O 2 k 21 [24]
223-body Neutral-Neutral Collision O + C O + C O 2 C O 2 + C O 2 k 22 [24]
Table 2. Species considered in the model.
Table 2. Species considered in the model.
Neutrals CO , CO 2 , O 2 ,   O , C
Pos. ions CO 2 + , C + , O + , O 2 +
Neg. ions O
Elec. excited CO 2 *
Vib. excited CO 2 ( V a d ) , CO 2 ( V 1 )
Table 3. The average consumed power and energy efficiency for different applied voltages.
Table 3. The average consumed power and energy efficiency for different applied voltages.
Applied Voltage (kV)Power (W)η (%)
1050.951.69
1261.602.89
1471.353.58
1682.144.23
1892.644.74
20105.495.11
22113.056.1
Table 4. Comparison of CO2 conversion factor and energy efficiency for different plasma technologies [7].
Table 4. Comparison of CO2 conversion factor and energy efficiency for different plasma technologies [7].
TechnologyRange of CO2 Conversion (%)Range of Energy Efficiency (%)
Dielectric Barrier Discharges0–300–15
Microwave Discharges0–800–40
Gliding Arc Discharges0–100–30
AC-PPP Reactor *2.5–251–6
* Results of this work.
Table 5. Detected peaks for the AC-PPP reactor with CO2 gas.
Table 5. Detected peaks for the AC-PPP reactor with CO2 gas.
Species λ (nm) Transition NameTransition Symbol v v
OH(Q)306–312Molecular Transition A 2 Σ X 2 Π 00
C O 2 326.1Fox, Duffendack and Barker’s System X 2 Π g A 2 Π u 31a
338.3Fox, Duffendack and Barker’s System X 2 Π g A 2 Π u 20
368.2Fox, Duffendack and Barker’s System X 2 Π g A 2 Π u 12
C O 2 + 387.5Fox, Duffendack and Barker’s System X 2 Π g A 2 Π u 24
407.3Fox, Duffendack and Barker’s System X 2 Π g A 2 Π u 14
H486.1Hβ Balmer series linen = 4 n = 2--
656.3Hα Balmer series linen = 3 n = 2--
C 2 471.6Swan System A 3 Π g B 3 Π u 01
516.4Swan System A 3 Π g B 3 Π u 00
560.1Swan System A 3 Π g B 3 Π u 10
C O 589.5High Pressure Bands A 3 Π g B 3 Π u 68
458.5The Triplet Bands d 3 Π a 3 Π 60
564.2The Triplet Bands d 3 Π a 3 Π 20
601.4The Triplet Bands d 3 Π a 3 Π 10
643.2The Triplet Bands d 3 Π a 3 Π 00
C O + 427.3Comet-tail system (First Negative System) X 2 Π A 2 Σ 20
O 2 385.3Chamberlains Airglow System a 1 Δ g C 3 Δ u 01
O 2 + 550–600Molecular emission Band B 4 Σ g A 4 Π u 01
O777.1Atomic Transition 3 p 5 P 3 s 5 S --
843.2Atomic Transition 3 p 3 P 3 s 3 S --
Table 6. The simulation parameters.
Table 6. The simulation parameters.
SymbolValue
Mobility, μe, m2/Vs (for electrons)1.28
Initial temperatures, T, eV (for electron and heavy particles)0.54, 0.026
Potential amplitude, KV22
Diffusion coefficient (for electrons)0.033
Initial density, cm−3 (charged species)1010
Initial density, cm−3 (background species)1019
Total time, ms20
Time steps, ns1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Barkhordari, A.; Karimian, S.; Rodero, A.; Krawczyk, D.A.; Mirzaei, S.I.; Falahat, A. Carbon Dioxide Decomposition by a Parallel-Plate Plasma Reactor: Experiments and 2-D Modelling. Appl. Sci. 2021, 11, 10047. https://doi.org/10.3390/app112110047

AMA Style

Barkhordari A, Karimian S, Rodero A, Krawczyk DA, Mirzaei SI, Falahat A. Carbon Dioxide Decomposition by a Parallel-Plate Plasma Reactor: Experiments and 2-D Modelling. Applied Sciences. 2021; 11(21):10047. https://doi.org/10.3390/app112110047

Chicago/Turabian Style

Barkhordari, Ali, Saeed Karimian, Antonio Rodero, Dorota Anna Krawczyk, Seyed Iman Mirzaei, and Amir Falahat. 2021. "Carbon Dioxide Decomposition by a Parallel-Plate Plasma Reactor: Experiments and 2-D Modelling" Applied Sciences 11, no. 21: 10047. https://doi.org/10.3390/app112110047

APA Style

Barkhordari, A., Karimian, S., Rodero, A., Krawczyk, D. A., Mirzaei, S. I., & Falahat, A. (2021). Carbon Dioxide Decomposition by a Parallel-Plate Plasma Reactor: Experiments and 2-D Modelling. Applied Sciences, 11(21), 10047. https://doi.org/10.3390/app112110047

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