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: 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: We present a zero-dimensional plasma kinetics model, including both surface and gas phase kinetics, to determine the role of vibrationally excited states in plasma-catalytic ammonia synthesis. We defined a new method to systematically capture the conditions of dielectric barrier discharges (DBDs), including those found in packed bed DBDs. We included the spatial and temporal nature of such discharges by special consideration of the number of micro-discharges in the model. We introduce a parameter that assigns only a part of the plasma power to the microdischarges, to scale the model conditions from filamentary to uniform plasma. Because of the spatial and temporal behaviour of the micro-discharges, not all micro-discharges occurring in the plasma reactor during a certain gas residence time are affecting the molecules. The fraction of power considered in the model ranges from 0.005 %, for filamentary plasma, to 100 %, for uniform plasma. If vibrational excitation is included in the plasma chemistry, these different conditions, however, yield an ammonia density that is only varying within one order of magnitude. At only 0.05 % of the power put into the uniform plasma component, a model neglecting vibrational excitation clearly does not result in adequate amounts of ammonia. Thus, our new model, which accounts for the concept in which not all the power is deposited by the micro-discharges, but some part may also be distributed in between them, suggests that vibrational kinetic processes are really important in (packed bed) DBDs. Indeed, vibrational excitation takes place in both the uniform plasma between the micro-discharges and in the strong micro-discharges, and is responsible for an increased N2 dissociation rate. This is shown here for plasma-catalytic ammonia synthesis, but might also be valid for other gas conversion applications.

Abstract: Plasma has been studied for several years to convert CO2 into value-added products. If CO2 could be converted in the presence of H2O as a cheap H-source for making syngas and oxygenates, it would mimic natural photosynthesis. However, CO2/H2O plasmas have not yet been extensively studied, not by experiments, and certainly not computationally. Therefore, we present here a kinetic modelling study to obtain a greater understanding of the vibrational kinetics of a CO2/H2O microwave plasma. For this purpose, we first created an electron impact cross section set for H2O, using a swarm-derived method. We added the new cross section set and CO2/H2O-related chemistry to a pure CO2 model. While it was expected that H2O addition mainly causes quenching of the CO2 asymmetric mode vibrational levels due to the additional CO2/H2O vibrational-translational relaxation, our model shows that the modifications in the vibrational kinetics are mainly induced by the strong electron dissociative attachment to H2O molecules, causing a reduction in electron density, and the corresponding changes in the input of energy into the CO2 vibrational levels by electron impact processes.

Abstract: We derive dispersion relations for a system of identical particles confined in a two-dimensional annular harmonic well and which interact through a Yukawa potential, e.g., a dusty plasma ring. When the particles are in a single chain (i.e., a one-dimensional ring), we find a longitudinal acoustic mode and a transverse optical mode which show approximate agreement with the dispersion relation for a straight configuration for large radii of the ring. When the radius decreases, the dispersion relations modify: there appears an anticrossing of the modes near the crossing point resulting in a frequency gap between the lower and upper branches of the modified dispersion relations. For the double chain (i.e., a two-dimensional zigzag configuration), the dispersion relation has four branches: longitudinal acoustic and optical and transverse acoustic and optical.

Abstract: The physics of plasmas containing positive and negative ions is discussed with special attention to the recently produced pair-ion plasma containing ions of equal mass and opposite charge. The effects of the density gradient in the direction perpendicular to the ambient magnetic field vector are discussed. The possible presence of electrons is discussed in the context of plasma modes propagating at an angle with respect to the magnetic field vector. It is shown that the electron plasma mode may become a backward mode in the presence of a density gradient, and this behavior may be controlled either by the electron number density or the mode number in the perpendicular direction. In plasmas with hot electrons an instability may develop, driven by the combination of electron collisions and the density gradient, and in the regime of a sound ions' response. In the case of a pure pair-ion plasma, for lower frequencies and for parameters close to those used in the recent experiments, the perturbed ions may feel the effects of the magnetic field. In this case the plasma mode also becomes backward, resembling features of an experimentally observed but yet unexplained backward mode. (C) 2008 American Institute of Physics.

Abstract: In this paper, a self-consistent fluid model is developed focusing on the plasma parameters in capacitively coupled 20% N 2-80% Ar discharges. Measurements of ion density are performed with the help of a floating double probe, and the emission intensities from Ar(4p) and N 2 ( B ) transitions are detected by an optical emission spectroscopy to estimate their relative densities. The consistency between the numerical and experimental results confirms the reliability of the simulation. Then the plasma characteristics, specifically the reaction mechanisms of active species, are analyzed under various voltages. The increasing voltage leads to a monotonous increase in species density, whereas a less homogeneous radial distribution is observed at a higher voltage. Due to the high concentration of Ar gas, Ar + becomes the main ion, followed by the N 2 +</mml:msubsup> ion. Besides the electron impact ionization of neutrals, the charge transfer processes of Ar +/ N 2 and N 2 +</mml:msubsup>/Ar are found to have an impact on the ionic species. The results indicate that adopting the lower charge transfer reaction rate coefficients weakens the Ar + ion density and yields a higher N 2 +</mml:msubsup> ion density. However, the effect on the species spatial distributions and other species densities is limited. As for the excited-state species, the electron impact excitation of background gases remains overwhelming in the formation of Ar(4p), N 2 ( B ), and N 2 ( a ' ), whereas the <mml:msub> N 2 ( A ) molecules are mainly formed by the decay of <mml:msub> N 2 ( B ). In addition, the dissociation of <mml:msub> N 2 collided by excited-state Ar atoms dominates the N generation, which are mostly depleted to produce N + ions.

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.