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

Parallel-Plate 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 CO 2 decomposition that works at atmospheric pressure with AC power. Results for pure CO 2 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 CO 2 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 CO 2 conversion and energy efﬁciency. 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 × 10 14 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 ﬂuid equations was developed for determining the main processes in the CO 2 splitting. The solution to this model, using the ﬁnite element method, gave the temporal and spatial behaviors of the formed species densities, the electric potential, and the temperatures of electrons.


Introduction
The mitigation of greenhouse gases is one of the most important challenges in the world today. In consequence, finding solutions for carbon dioxide (CO 2 ) reduction is a priority research topic. Technologies for the decomposition of these gases after their generation are being developed, together with solutions to reduce CO 2 production at source.
The applicability of high voltage electrical discharge plasmas to reform carbon dioxide (CO 2 ) as an environmental pollutant has been extensively studied [1][2][3][4][5]. Plasma discharge 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.

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). 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). 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 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 CO 2 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 (V diff > V drif >> V axi ).
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 CO 2 and the formation of CO and O 2 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 CO 2 , CO, O 2 , etc.

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  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]: 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 the reaction rate depends on the density of each species, nA and nB.
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]: The collisional cross section can be written as follows: 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

For a typical reaction between species
the reaction rate depends on the density of each species, n A and 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]: The collisional cross section can be written as follows: where I p 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), B i (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 CO 2 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 CO 2 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 nonasymmetric excited vibrational modes; V1 is the first level of the asymmetric vibrational mode, and the excited electronic states of the CO 2 are grouped in the effective level CO 2 *.  [24] 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, D i , i.e., 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: 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 CO 2 + and C + species were 0.0012 and 0.0009 m 2 /Vs, respectively. The diffusion coefficient of the electrons and ions were instead calculated from the Einstein relation: with T e(i) and q e(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]: where ω e is the electron energy density, R ω is the energy loss or gain due to inelastic collisions, the term eE·Γ 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: The electron energy loss or gain R ω is obtained by summing the collisional energy loss or gain over all reactions [14]: where x j is the mole fraction of the target species for reaction j, k j is the rate coefficient for reaction j, N n 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 T e are correlated with each other through [30]: For non-electron species, the following equation was solved for the mass fraction of each species [30]: 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]: 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]: 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, n k , and their charge, eZ k [31]: 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.

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: where n is the normal vector pointing toward the tube wall and, v th(i) is the thermal velocity of particles [32]: and the number a i is defined by: 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 th,e n e − ∑ p γ p Γ p ·n (20) 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]: 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: and the electric potential in the powered electrode was given by: Appl. Sci. 2021, 11, 10047 9 of 25 with f = 50 Hz and V 0 = 22 kV.

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 CO 2 conversion factor can be calculated from the areas of CO 2 peak, [CO 2 ] in and [CO 2 ] out by [33]: In a similar way, the selectivity of CO and O 2 species were calculated by [32]: Figure 4a shows the conversion factor of CO 2 and the selectivity of CO and O 2 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).

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, [ 2 ] and [ 2 ] by [33]: In a similar way, the selectivity of CO and O2 species were calculated by [32]: × 100 × 100 (25) 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]: where Vm is the molar volume, that can be considered 24.5 L/mol, ∆ 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: 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. The energy efficiency of this technique is given by [32]: where V m is the molar volume, that can be considered 24.5 L/mol, ∆H R is the reaction enthalpy for CO 2 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: 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 CO 2 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 CO 2 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. Table 4. Comparison of CO 2 conversion factor and energy efficiency for different plasma technologies [7].

Range of CO 2 Conversion (%) Range of Energy Efficiency (%)
Dielectric 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 CO 2 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 CO 2 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.

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 CO 2 plasma reactor. Vibro-rotational bands of CO 2 , CO 2 + , C 2 , 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. Table 5. Detected peaks for the AC-PPP reactor with CO 2 gas.

Species λ (nm) Transition Name Transition Symbol
338.3 Fox, Duffendack and Barker's System Fox, Duffendack and Barker's System A study of the dependencies of the normalized intensities of the headband of vibrorotational 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]. Appl Atomic Transition 3p 3 P → 3s 3 S -- A study of the dependencies of the normalized intensities of the headband of vibrorotational 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 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 ( 2 → 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.

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]: 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]: 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 cm 3 ), 〈R 2 〉 is the difference of the squared radius of the 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 CO 2 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]: 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]: ∆λ Pressure = ∆λ van der Waals = 8.18 × 10 −26 λ 2 α R 2 2/5 T g µ 3/10 P kT g (29) where µ is the reduced mass of the colliding particles (µ= 0.92308 for H/CO 2 collision), λ is the wavelength of the H β line (486.1 nm), α is the molecular polarizability of the CO 2 disturbing particles (α= 2.46 × 10 −24 cm 3 ), R 2 is the difference of the squared radius of the upper and lower levels of H β transition, T g 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]: 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]: where λ is the wavelength of β hydrogen line (486.1 nm), M is the molar weight of hydrogen, which is equal to 1, and T g is the gas temperature in Kelvin units. The total Gaussian profile is generally considered as a convolution of Doppler and instrumental profiles, then [43]: 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. 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]: 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]: 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]: 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 (30) and (28), 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 × 10 14 cm −3 . This high value supports the electron 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 × 10 14 cm −3 . This high value supports the electron kinetics for CO 2 decomposition assumed for this reactor.

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 CO 2 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.

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  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. 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).

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. 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).

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 CO 2 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.

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. 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.

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 exc > T v > T r ∼ T g [47]. Appl. Sci. 2021, 11, x FOR PEER REVIEW 18 of 25 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: ≳ > >~ [47].

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

CO 2 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 CO 2 ground state molecules. As predicted by the model, the most relevant mechanisms of CO 2 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 CO 2 + 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 CO 2 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: CO 2 + , CO, O, O 2 , and O − were studied to identify the reactions that dominated in each part of the reactor.
Density of CO 2 + Ion Figure 13 shows the spatial distribution of the density of CO 2 + 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 CO 2 + in the direction of this field. If a CO 2 + 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 CO 2 + ions are less effective than the electrons in the impacting ionization process. As shown in Figure 13, the CO 2 + ions accumulated near the ground electrode during the positive voltage applied, and they accumulated in the powered electrode with the negative voltages. 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 Figures 14 and 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 the CO 2 + 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 CO 2 + ion is an unstable species and will be quickly recombined with electrons or negative ions O 2 − to convert to CO and O, and/or O 2 [48].

Density of Neutral Species CO, O, and O 2
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 Figures 14 and 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 × 10 11 cm −3 and 2.3 × 10 11 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 CO 2 when the residence time is high [48]. However, the O atoms recombine quickly into oxygen molecules, O 2 and/or ozone, O 3 .
Oxygen molecules, O 2 , are formed mainly by the reaction of oxygen atoms with the CO 2 gas through neutral dissociation and 3-body recombination, reactions 17 and 20, respectively. Figure  The maximum values of CO and O densities are equal to 4.7 × 10 11 cm −3 and 2.3 × 10 11 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. The maximum values of CO and O densities are equal to 4.7 × 10 11 cm −3 and 2.3 × 10 11 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.

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. 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. 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 according to electron density. This negative oxygen ion density increased quickl became constant after the first millisecond.
The electron dissociative attachment that produced the O − ions was higher inlet/outlet pipes. The high densities and low temperatures in these positions favore 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 CO 2 Molecules
Reactions 7, 9, and 15 mainly describe the CO 2 conversion and the density distribution of this molecule in the reactor. The spatial distribution of CO 2 density between the two electrodes inside the AC-PPP reactor is presented in Figure 17. The CO 2 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 CO 2 was also discomposed inside these pipes. 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 CO 2 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 CO 2 conversion factor was made with the time averages of the CO 2 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).

Conclusions
In this study, a new AC-PPP reactor for CO 2 remediation was presented. We undertook experimental measurements by gas chromatography and optical emission spectroscopy to characterize and determine its capacity for CO 2 decomposition. The CO 2 conversion and energy efficiency were obtained by the GC method and compared with other plasma sources. It was found that CO 2 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 × 10 14 cm −3 .
A 2-D fluid model was developed to simulate the kinetics of CO 2 decomposition in this plasma. The temporal and spatial behaviors of the densities of CO 2 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 CO 2 decomposition by increasing collisions.
The model has demonstrated that the main mechanisms of CO 2 splitting were the electron impact dissociation, the electron ionization, and the electron dissociative attachment. The species O − and CO 2 + ions, CO and O 2 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 CO 2 conversion by this reactor.
The distribution of the CO 2 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 CO 2 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 CO 2 and its transformation into other species. This type of reactor could contribute to the optimization of working conditions and improvements in reactor design.

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