Abstract: CO2 conversion into value-added products has gained significant interest over the few last years, as the greenhouse gas concentrations constantly increase due to anthropogenic activities. Here we report on experiments for CO2 conversion by means of a cold atmospheric plasma using a cylindrical flowing dielectric barrier discharge (DBD) reactor. A detailed comparison of this DBD ignited in a so-called burst mode (i.e. where an AC voltage is applied during a limited amount of time) and pure AC mode is carried out to evaluate their effect on the conversion of CO2 as well as on the energy efficiency. Decreasing the duty cycle in the burst mode from 100% (i.e. corresponding to pure AC mode) to 40% leads to a rise in the
conversion from 16–26% and to a rise in the energy efficiency from 15 to 23%. Based on a detailed electrical analysis, we show that the conversion correlates with the features of the microfilaments. Moreover, the root-mean-square voltage in the burst mode remains constant as a function of the process time for the duty cycles <70%, while a higher duty cycle or the usual pure AC mode leads to a clear voltage decay by more than 500 V, over approximately 90 s, before reaching a steady state regime. The higher plasma voltage in the burst mode yields a higher electric field. This causes the increasing the electron energy, and therefore their
involvement in the CO2 dissociation process, which is an additional explanation for the higher CO2 conversion and energy efficiency in the burst mode.

Abstract: Electrical breakdown by the application of an electric field occurs more easily in hot gases than in cold gases because of the extra electron-species interactions that occur as a result of dissociation, ionization and excitation at higher temperature. This paper discusses some overlooked physics and clarifies inaccuracies in the evaluation of the effective ionization coefficients and the critical reduced breakdown electric field of CO2 at elevated temperature, considering the influence of excited states and ion kinetics. The critical reduced breakdown electric field is obtained by balancing electron generation and loss mechanisms using the electron energy distribution function (EEDF) derived from the Boltzmann transport equation under the two-term approximation. The equilibrium compositions of the hot gas mixtures are determined based on Gibbs free energy minimization considering the ground states as well as vibrationally and electronically excited states as independent species, which follow a Boltzmann distribution with a fixed excitation temperature. The interaction cross sections between electrons and the excited species, not reported previously, are properly taken into account. Furthermore, the ion kinetics, including electron–ion recombination, associative electron detachment, charge transfer and ion conversion into stable negative ion clusters, are also considered. Our results indicate that the excited species lead to a greater population of high-energy electrons at higher gas temperature and this affects the Townsend rate coefficients (i.e. of electron impact ionization and attachment), but the critical reduced breakdown electric field strength of CO2 is only affected when also properly accounting for the ion kinetics. Indeed, the latter greatly influences the effective ionization coefficients and hence the critical reduced breakdown electric field at temperatures above 1500 K. The rapid increase of the dissociative electron attachment cross-section of molecular oxygen with rising vibrational quantum number leads to a larger electron loss rate and this enhances the critical reduced breakdown electric field strength in the temperature range where the concentration of molecular oxygen is relatively high. The results obtained in this work show reasonable agreement with experimental results from literature, and are important for the evaluation of the dielectric strength of CO2 in a highly reactive environment at elevated temperature.

Abstract: In this study we report on a 2D fluid model of a gliding arc discharge in argon. Despite the 3D nature of the discharge, 2D models are found to be capable of providing very useful information about the operation of the discharge. We employ two modelsan axisymmetric and a Cartesian one. We show that for the considered experiment and the conditions of a low current arc (around 30 mA) in argon, there is no significant heating of the cathode surface and the discharge is sustained by field electron emission from the cathode accompanied by the formation of a cathode spot. The obtained discharge power and voltage are relatively sensitive to the surface properties and particularly to the surface roughness, causing effectively an amplification of the normal electric field. The arc body and anode region are not influenced by this and depend mainly on the current value. The gliding of the arc is modelled by means of a 2D Cartesian model. The arcelectrode contact points are analysed and the gliding mechanism along the electrode surface is discussed. Following experimental observations, the cathode spot is simulated as jumping from one point to another. A complete arc cycle is modelled from initial ignition to arc decay. The results show that there is no interaction between the successive gliding arcs.

Abstract: We use a zero-dimensional reaction kinetics model to simulate CO2 conversion in microwave discharges where the excitation of the vibrational levels plays a significant role in the dissociation kinetics. The model includes a description of the CO2 vibrational kinetics, taking into account state-specific VT and VV relaxation reactions and the effect of vibrational excitation on other chemical reactions. The model is used to simulate a general tubular microwave reactor, where a stream of CO2 flows through a plasma column generated by microwave radiation. We study the effects of the internal plasma parameters, namely the reduced electric field, electron density and the total specific energy input, on the CO2 conversion and its energy efficiency. We report the highest energy efficiency (up to 30%) for a specific energy input in the range 0.41.0 eV/molecule and a reduced electric field in the range 50100 Td and for high values of the electron density (an ionization degree greater than 10−5). The energy efficiency is mainly limited by the VT relaxation which contributes dominantly to the vibrational energy losses and also contributes significantly to the heating of the reacting gas. The model analysis provides useful insight into the potential and limitations of CO2 conversion in microwave discharges.

Abstract: The conversion of methane to value-added chemicals and fuels is considered to be one of the challenges of the 21st century. In this paper we study, by means of fluid modeling, the conversion of methane to higher hydrocarbons or oxygenates by partial oxidation with CO2 or O2 in a dielectric barrier discharge. Sixty-nine different plasma species (electrons, ions, molecules, radicals) are included in the model, as well as a comprehensive set of chemical reactions. The calculation results presented in this paper include the conversion of the reactants and the yields of the reaction products as a function of residence time in the reactor, for different gas mixing ratios. Syngas (i.e. H2 + CO) and higher hydrocarbons (C2Hx) are typically found to be important reaction products.

Abstract: The electron bounce resonance heating (BRH) in dual-frequency capacitively coupled plasmas operated in oxygen is studied by different experimental methods and a particle-in-cell/Monte Carlo collision (PIC/MCC) simulation, and compared with the electropositive argon discharge. In comparison with argon, the experimental results show that in an oxygen discharge the resonance peaks in positive-ion density and light intensity tend to occur at larger electrode gaps. Moreover, at electrode gaps L > 2.5 cm, the positive-ion (and electron) density and the light emission drop monotonically in the oxygen discharge upon increasing L, whereas they rise (after an initial drop) in the argon case. At resonance gap the electronegativity reaches its maximum due to the BRH. All these experimental observations are explained by PIC/MCC simulations, which show that in the oxygen discharge the bulk electric field becomes quite strong and is out of phase with the sheath field. Therefore, it retards the resonance electrons when traversing the bulk, resulting in a suppressed BRH. Both experiment and simulation results show that this effect becomes more pronounced at lower high-frequency power, when the discharge mode changes from electropositive to electronegative. In a pure oxygen discharge, the BRH is suppressed with increasing pressure and almost diminishes at 12 Pa. Finally, the driving frequency significantly affects the BRH, because it determines the phase relation between bulk electric field and sheath electric field.

Abstract: In this work particle-in-cell/Monte Carlo collision simulations are performed to study the heating mechanism and plasma characteristics in direct current (dc) superposed radio-frequency (RF) capacitively coupled plasmas, operated both in single-frequency (SF) and dual-frequency (DF) regimes. An RF (60/2 MHz) source is applied on the bottom electrode to sustain the discharge, and a dc source is fixed on the top electrode. The heating mechanism appears to be very different in dc superposed SF and DF discharges. When only a single source of 60 MHz is applied, the plasma bulk region is reduced by the dc source, thus the ionization rate and hence the electron density decrease with rising dc voltage. However, when a DF source of 60 and 2 MHz is applied, the electron density can increase upon addition of a dc voltage, depending on the gap length and applied dc voltage. This is explained from the spatiotemporal ionization rates in the DF discharge. In fact, a completely different behavior is observed for the ionization rate in the two half-periods of the LF source. In the first LF half-period, the situation resembles the dc superposed SF discharge, and the reduced plasma bulk region due to the negative dc bias results in a very small effective discharge area and a low ionization rate. On the other hand, in the second half-period, the negative dc bias is to some extent counteracted by the LF voltage, and the sheath close to the dc electrode becomes particularly thin. Consequently, the amplitude of the high-frequency sheath oscillations at the top electrode is largely enhanced, while the LF sheath at the bottom electrode is in its expanding phase and can thus well confine the high-energy electrons. Therefore, the ionization rate increases considerably in this second LF half-period. Furthermore, in addition to the comparison between SF and DF discharges and the effect of gap length and dc voltage, the effect of secondary electrons is examined.

Abstract: Understanding the nature and effect of the multitude of plasma–surface interactions in plasma catalysis is a crucial requirement for further process development and improvement. A particularly intriguing and rather unique property of a plasma-catalytic setup is the ability of the plasma to modify the electronic structure, and hence chemical properties, of the catalyst through charging, i.e. the absorption of excess electrons. In this work, we develop a quantum chemical model based on density functional theory to study excess negative surface charges in a heterogeneous catalyst exposed to a plasma. This method is specifically applied to investigate plasma-catalytic CO2 activation on supported M/Al2O3 (M=Ti, Ni, Cu) single atom catalysts. We find that (1) the presence of a negative surface charge dramatically improves the reductive power of the catalyst, strongly promoting the splitting of CO2 to CO and oxygen, and (2) the relative activity of the investigated transition metals is also changed upon charging, suggesting that controlled surface charging is a powerful additional parameter to tune catalyst activity and selectivity. These results strongly point to plasma-induced surface charging of the catalyst as an important factor contributing to the plasma-catalyst synergistic effects frequently reported for plasma catalysis.

Abstract: A gliding arc plasmatron (GAP) is very promising for CO2 conversion into value-added chemicals, but to further improve this important application, a better understanding of the arc behavior is indispensable. Therefore, we study here for the first time the dynamic arc behavior of the GAP by means of a high-speed camera, for different reactor configurations and in a wide range of operating conditions. This allows us to provide a complete image of the behavior of the gliding arc. More specifically, the arc body shape, diameter, movement and rotation speed are analyzed and discussed. Clearly, the arc movement and shape relies on a number of factors, such as gas turbulence, outlet diameter, electrode surface, gas contraction and buoyance force. Furthermore, we also compare the experimentally measured arc movement to a state-of-the-art 3D-plasma model, which predicts the plasma movement and rotation speed with very good accuracy, to gain further insight in the underlying mechanisms. Finally, we correlate the arc dynamics with the CO2 conversion and energy efficiency, at exactly the same conditions, to explain the effect of these parameters on the CO2 conversion process. This work is important for understanding and optimizing the GAP for CO2 conversion.

Abstract: Low-temperature plasma chemical kinetic models are particularly important to the plasma community. These models typically require dozens of inputs, especially rate coefficients. The latter are not always precisely known and it is not surprising that the error on the rate coefficient data can propagate to the model output. In this paper, we present a model that uses N = 400 different combinations of rate coefficients based on the uncertainty attributed to each rate coefficient, giving a good estimation of the uncertainty on the model output due to the rate coefficients. We demonstrate that the uncertainty varies a lot with the conditions and the type of output. Relatively low uncertainties (about 15%) are found for electron density and temperature, while the uncertainty can reach more than an order of magnitude for the population of the vibrational levels in some cases and it can rise up to 100% for the CO2 conversion. The reactions that are mostly responsible for the largest uncertainties are identified. We show that the conditions of pressure, gas temperature and power density have a great effect on the uncertainty and on which reactions lead to this uncertainty. In all the cases tested here, while the absolute values may suffer from large uncertainties, the trends observed in previous modeling work are still valid. Finally, in accordance with the work of Turner, a number of ‘good practices’ is recommended.

Abstract: Plasma catalysis is gaining increasing interest for various environmental applications. Catalytic
material can be inserted in different shapes in the plasma, e.g., as pellets, (coated) beads, but also
as honeycomb monolith and 3DFD structures, also called ‘structured catalysts’, which have high
mass and heat transfer properties. In this work, we examine the streamer discharge propagation
and the interaction between plasma and catalysts, inside the channels of such structured catalysts,
by means of a two-dimensional particle-in-cell/Monte Carlo collision model. Our results reveal
that plasma streamers behave differently in various structured catalysts. In case of a honeycomb
structure, the streamers are limited to only one channel, with low or high plasma density when
the channels are parallel or perpendicular to the electrodes, respectively. In contrast, in case of a
3DFD structure, the streamers can distribute to different channels, causing discharge
enhancement due to surface charging on the dielectric walls of the structured catalyst, and
especially giving rise to a broader plasma distribution. The latter should be beneficial for plasma
catalysis applications, as it allows a larger catalyst surface area to be exposed to the plasma.

Abstract: In this work, a hybrid model is used to investigate the effect of different gas ratios on the Si etching and polymer film deposition characteristics in an Ar/CF4 inductively coupled plasma. The influence of the surface processes on the bulk plasma properties is studied, and also the spatial characteristics of important gas phase and etched species. The densities of F and CF2 decrease when the surface module is included in the simulations, due to the species consumption caused by etching and polymer deposition. The influence of the surface processes on the bulk plasma depends on the Ar/CF4 gas ratio. The deposited polymer becomes thicker at high CF4 content because of more abundant CFx radicals. As a result of the competition between the polymer thickness and the F flux, the etch rate first increases and then decreases upon increasing the CF4 content. The electron properties, more specifically the electron density profile, affect the Si etch characteristics substantially by determining the radical density and flux profiles. In fact, the radial profile of the etch rate is more uniform at low CF4 content since the electron density has a smooth distribution. At high CF4 content, the etch rate is less uniform with a minimum halfway along the wafer radius, because the electron density distribution is more localized. Therefore, our calculations predict that it is better to work at relatively high Ar/CF4 gas ratios, in order to obtain high etch rate and good profile uniformity for etch applications. This, in fact, corresponds to the typical experimental etch conditions in Ar/CF4 gas mixtures as found in the literature, where Ar is typically present at a much higher concentration than CF4.