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: Packed bed plasma reactors (PBPRs) are gaining increasing interest for use in environmental applications, such as greenhouse gas conversion into value-added chemicals or renewable fuels and volatile pollutant removal (e.g. NOx, VOC, K), as they enhance the conversion and energy efficiency of the process compared to a non-packed reactor. However, the plasma behaviour in a PBPR is not well understood. In this paper we demonstrate, by means of a fluid model, that the discharge behaviour changes considerably when changing the size of the packing beads and their dielectric constant, while keeping the interelectrode spacing constant. At low dielectric constant, the plasma is spread out over the full discharge gap, showing significant density in the voids as well as in the connecting void channels. The electric current profile shows a strong peak during each half cycle. When the dielectric constant increases, the plasma becomes localised in the voids, with a current profile consisting of many smaller peaks during each half cycle. For large bead sizes, the shift from full gap discharge to localised discharges takes place at a higher dielectric constant than for smaller beads. Furthermore, smaller beads or beads with a lower dielectric constant require a higher breakdown voltage to cause plasma formation.

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 studied the evolution of an isolated pulsed plasma in a vortex flow stabilised microwave (MW) discharge in N2 at 25 mbar via the combination of 0D kinetics modelling, iCCD imaging and laser scattering diagnostics. Quenching of electronically excited N2 results in fast gas heating and the onset of a thermal-ionisation instability, contracting the discharge volume. The onset of a thermal-ionisation instability driven by vibrational excitation pathways is found to facilitate significantly higher N2 conversion (i.e. dissociation to atomic N2 ) compared to pre-instability conditions, emphasizing the potential utility of this dynamic in future fixation applications. The instability onset is found to be instigated by super-elastic heating of the electron energy distribution tail via vibrationally excited N2 . Radial contraction of the discharge to the skin depth is found to occur post instability, while the axial elongation is found to be temporarily contracted during the thermal instability onset. An increase in power reflection during the thermal instability onset eventually limits the destabilising effects of exothermic electronically excited N2 quenching. Translational and vibrational temperature reach a quasi-non-equilibrium after the discharge contraction, with translational temperatures reaching ∼1200 K at the pulse end, while vibrational temperatures are found in near equilibrium with the electron energy (1 eV, or ∼11 600 K). This first description of the importance of electronically excited N2 quenching in thermal instabilities gives an additional fundamental understanding of N2 plasma behaviour in pulsed MW context, and thereby brings the eventual implementation of this novel N2 fixation method one step closer.

Abstract: A combined Monte Carlo (MC)/analytical surface model is developed to study the plasma processes occurring during the reactive sputter deposition of TiOx thin films. This model describes the important plasma species with a MC approach (i.e. electrons, Ar+ ions, {\rm O}_2

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 behaviour in an Ar/CF4 inductively coupled plasma is investigated by a Langmuir probe and a hybrid model. The simulated and measured results include electron density, temperature and electron energy distribution function for different values of Ar/CF4 ratio, coil power and gas pressure. The hybrid plasma equipment model simulations show qualitative agreement with experiment. The effect of F2 electron attachment on the electron behaviour is explored by comparing two sets of data based on different F atom boundary conditions. It is demonstrated that electron attachment at F2 molecules is responsible for the depletion of low-energy electrons, causing a density decrease as well as a temperature increase when CF4 is added to an Ar plasma.