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: Plasma catalysis is gaining increasing interest, but the underlying mechanisms are far from understood. Different catalyst materials will have different chemical effects, but in addition, they might also have different dielectric constants, which will affect surface charging, and thus the plasma behavior. In this work, we demonstrate that surface charging plays an important role in the streamer propagation and discharge enhancement inside catalyst pores, and in the plasma distribution along the dielectric surface, and this role greatly depends on the dielectric constant of the material. For εr50, surface charging causes the plasma to spread along the dielectric surface and inside the pores, leading to deeper plasma streamer penetration, while for εr>50 or for metallic coatings, the discharge is more localized, due to very weak surface charging. In addition, at εr=50, the significant surface charge density near the pore entrance causes a large potential drop at the sharp pore edges, which induces a strong electric field and results in most pronounced plasma enhancement near the pore entrance.

Abstract: In this work we study energy transport in a gliding arc discharge with two diverging flat
electrodes in argon gas at atmospheric pressure. The discharge is ignited at the shortest electrode
gap and it is pushed downstream by a forced gas flow. The current values considered are
relatively low and therefore a non-equilibrium plasma is produced. We consider two cases, i.e.
with high and low discharge current (28 mA and 2.8mA), and a constant gas flow of 10 lmin −1 ,
with a significant turbulent component to the velocity. The study presents an analysis of the
various energy transport mechanisms responsible for the redistribution of Joule heating to the
plasma species and the moving background gas. The objective of this work is to provide a
general understanding of the role of the different energy transport mechanisms in arc formation
and sustainment, which can be used to improve existing or new discharge designs. The work is
based on a three-dimensional numerical model, combining a fluid plasma model, the shear stress
transport Reynolds averaged Navier–Stokes turbulent gas flow model, and a model for gas
thermal balance. The obtained results show that at higher current the discharge is constricted
within a thin plasma column several hundred kelvin above room temperature, while in the low-
current discharge the combination of intense convective cooling and low Joule heating prevents
discharge contraction and the plasma column evolves to a static non-moving diffusive plasma,
continuously cooled by the flowing gas. As a result, the energy transport in the two cases is
determined by different mechanisms. At higher current and a constricted plasma column, the
plasma column is cooled mainly by turbulent transport, while at low current and an unconstricted
plasma, the major cooling mechanism is energy transport due to non-turbulent gas convection. In
general, the study also demonstrates the importance of turbulent energy transport in
redistributing the Joule heating in the arc and its significant role in arc cooling and the formation
of the gas temperature profile. In general, the turbulent energy transport lowers the average gas
temperature in the arc, thus allowing additional control of thermal non-equilibrium in the
discharge.

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: The vibrational level populations of the electronic ground state of the nitrogen molecule have been calculated for typical glow discharge conditions in argonnitrogen mixtures with nitrogen concentrations between 0.1 and 1%. Stationary solutions of the master equations of the vibrational levels have been obtained using numerical methods. The main mechanisms responsible for the population and depopulation of the vibrational levels, and for the overall shape of the vibrational distribution function are pointed out. It has been found that vibrationvibration collisions play only a minor role and therefore the population of the vibrational levels is basically determined by the electron temperature.

Abstract: In this study, we present a computational model of a cylindrical electric probe in atmospheric pressure argon plasma. The plasma properties are varied in terms of density and electron temperature. Furthermore, results for plasmas with Maxwellian and non-Maxwellian electron energy distribution functions are also obtained and compared. The model is based on the fluid description of plasma within the COMSOL software package. The results for the ion saturation current are compared and show good agreement with existing analytical Langmuir probe theories. A strong dependence between the ion saturation current and electron transport properties was observed, and attributed to the effects of ambipolar diffusion.

Abstract: In this study we quantitatively investigate for the first time the plasma characteristics of an argon gliding arc with a 3D model. The model is validated by comparison with available experimental data from literature and a reasonable agreement is obtained for the calculated gas temperature and electron density. A complete arc cycle is modeled from initial ignition to arc decay. We investigate how the plasma characteristics, i.e., the electron temperature, gas temperature,
reduced electric field, and the densities of electrons, Ar+ and Ar2+ ions and Ar(4s) excited states, vary over one complete arc cycle, including their behavior in the discharge and post-discharge. These plasma characteristics exhibit a different evolution over one arc cycle, indicating that either the active discharge stage or the post-discharge stage can be beneficial for certain applications.

Abstract: In recent years there has been growing interest in the use of plasma technology for CO2 conversion. To improve this application, a good insight into the underlying mechanisms is of great importance. This can be obtained from modeling the detailed plasma chemistry in order to understand the chemical reaction pathways leading to CO2 conversion (either in pure form or mixed with another gas). Moreover, in practice, several plasma reactor types are being investigated for CO2 conversion, so in addition it is essential to be able to model these reactor geometries so that their design can be improved, and the most energy efficient CO2 conversion can be achieved. Modeling the detailed plasma chemistry of CO2 conversion in complex reactors is, however, very time-consuming. This problem can be overcome by using a combination of two different types of model: 0D chemical reaction kinetics models are very suitable for describing the detailed plasma chemistry, while the characteristic features of different reactor geometries can be studied by 2D or 3D fluid models. In the first instance the latter can be developed in argon or helium with a simple chemistry to limit the calculation time; however, the ultimate aim is to implement the more complex CO2 chemistry in these models. In the present paper, examples will be given of both the 0D plasma chemistry models and the 2D and 3D fluid models for the most common plasma reactors used for CO2 conversion in order to emphasize the complementarity of both approaches. Furthermore, based on the modeling insights, the paper discusses the possibilities and limitations of plasma-based CO2 conversion in different types of plasma reactors, as well as what is needed to make further progress in this field.

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: Particle-in-cell/Monte Carlo collision (PIC/MCC) simulations are performed to investigate the asymmetric secondary electron emission (SEE) effects when electrons strike two different material electrodes in low pressure capacitively coupled plasmas (CCPs). To describe the electron-surface interactions, a realistic model, considering the primary electron impact energy and angle, as well as the corresponding surface property-dependent secondary electron yields, is employed in PIC/MCC simulations. In this model, three kinds of electrons emitted from the surface are considered: (i) elastically reflected electrons, (ii) inelastically backscattered electrons, and (iii) electron induced secondary electrons (SEs, i.e., delta-electrons). Here, we examined the effects of electron-surface interactions on the ionization dynamics and plasma characteristics of an argon discharge. The discharge is driven by a voltage source of 13.56MHz with amplitudes in the range of 200-2000V. The grounded electrode material is copper (Cu) for all cases, while the powered electrode material is either Cu or silicon dioxide (SiO2). The simulations reveal that the electron impact-induced SEE is an essential process at low pressures, especially at high voltages. Different electrode materials result in an asymmetric response of SEE. Depending on the instantaneous local sheath potential and the phase of the SEE, these SEs either are reflected by the opposite sheath or strike the electrode surface, where they can induce delta-electrons upon their residual energies. It is shown that highly energetic delta-electrons contribute significantly to the ionization rate and a self-bias forms when the powered electrode material is assumed to be made of SiO2. Complex dynamics is observed due to the multiple electron-surface interaction processes and asymmetric yields of SEs in CCPs.

Abstract: The present study analyses a specific laser ablation cell, the High Efficiency Aerosol Dispersion (HEAD) cell (see J. Pisonero et al., J. Anal. At. Spectrom., 2006, 21, 922931), by means of computational fluid dynamics (CFD) simulations. However, this cell consists of different modular parts, therefore, the results are probably of interest for the further development of other ablation cells. In the HEAD cell, the ablation spot is positioned below an orifice in the ceiling of the sample chamber. The particle transport through this orifice has been analysed for a ceiling height of 0.8 mm. The critical velocity for the onset of particle losses was found to be independent on the ejection angle at the crater spot. The deceleration of the particles can be described as the stopping in an effectively steady gas. Particle losses were negligible in this modular part of the cell at the evaluated laser ablation conditions. The transport efficiency through the Venturi chamber was investigated for different sample gas flow rates. In this case, it was found that small particles were predominantly lost at low flow rates, the large particles at higher flow rates. Making use of the simulation results, it was possible to design a modification of the HEAD cell that results in extremely short calculated washout times. The simulations yielded a signal of less than 10 ms, which was produced by more than 99% of the introduced mass.

Abstract: A packed bed dielectric barrier discharge plasma reactor is computationally studied with a fluid model. Two different complementary axisymmetric 2D geometries are used to mimic the intrinsic 3D problem. It is found that a packing enhances the electric field strength and electron temperature at the contact points of the dielectric material due to polarization of the beads by the applied potential. As a result, these contact points prove to be of direct importance to initiate the plasma. At low applied potential, the discharge stays at the contact points, and shows the properties of a Townsend discharge. When a high enough potential is applied, the plasma will be able to travel through the gaps in between the beads from wall to wall, forming a kind of glow discharge. Therefore, the inclusion of a so-called ‘channel of voids’ is indispensable in any type of packed bed modelling.

Abstract: In this experimental study, a flowing dielectric barrier discharge operating at atmospheric pressure is used for the splitting of CO2 into O2 and CO. The influence of the applied frequency and plasma power on the microdischarge properties is investigated to understand their role on the CO2 conversion. Electrical measurements are carried out to explain the conversion trends and to characterize the microdischarges through their number, their lifetime,
their intensity and the induced electrical charge. Their influence on the gas and electrode temperatures is also evidenced through optical emission spectroscopy and infrared imaging. It is shown that, in our configuration, the conversion depends mostly on the charge delivered in the plasma and not on the effective plasma voltage when the applied power is modified. Similarly, at constant total current, a better conversion is observed at low frequencies, where a less filamentary discharge regime with a higher effective plasma voltage than that at a higher
frequency is obtained.

Abstract: Dielectric barrier discharges (DBDs) are commonly used to generate cold plasmas at
atmospheric pressure. Whatever their configuration (tubular or planar), the presence of a dielectric barrier is mandatory to prevent too much charge build up in the plasma and the formation of a thermal arc. In this article, the role of the barrier thickness (2.0, 2.4 and 2.8 mm) and of the kind of dielectric material (alumina, mullite, pyrex, quartz) is investigated on the filamentary behavior in the plasma and on the CO2 conversion in a tubular flowing DBD, by means of mass spectrometry measurements correlated with electrical characterization and IR imaging. Increasing the barrier thickness decreases the capacitance, while preserving the electrical charge. As a result, the voltage over the dielectric increases and a larger number of microdischarges is generated, which enhances the CO2 conversion. Furthermore, changing the dielectric material of the barrier, while keeping the same geometry and dimensions, also affects the CO2 conversion. The highest CO2 conversion and energy efficiency are obtained for quartz and alumina, thus not following the trend of the relative permittivity. From the
electrical characterization, we clearly demonstrate that the most important parameters are the somewhat higher effective plasma voltage (yielding a somewhat higher electric field and electron energy in the plasma) for quartz, as well as the higher plasma current (and thus larger electron density) and the larger number of microdischarge filaments (mainly for alumina, but also for quartz). The latter could be correlated to the higher surface roughness for alumina and to the higher voltage over the dielectric for quartz.