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: 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: Nanosecond pulsed discharge plasma shows a high degree of non-equilibrium, and exhibits relatively high conversions in the dry reforming of methane. To further improve the application, a good insight of the underlying mechanisms is desired. We developed a chemical kinetics model to explore the underlying plasma chemistry in nanosecond pulsed discharge. We compared the calculated conversions and product selectivities with experimental results, and found reasonable agreement in a wide range of specific energy input. Hence, the chemical kinetics model is able to provide insight in the underlying plasma chemistry. The modeling results predict that the most important dissociation reaction of CO₂and CH₄is electron impact dissociation. C₂H₂is the most abundant hydrocarbon product, and it is mainly formed upon reaction of two CH₂radicals. Furthermore, the vibrational excitation levels of CO₂contribute for 85% to the total dissociation of CO₂.

Abstract: The catalytic oxidative dehydrogenation of hydrocarbons is an industrially important process, in which selectivity is a key issue. We here investigate the conversion of methanol to formaldehyde on a vanadia surface employing long timescale simulations, reaching a time scale of seconds. In particular, we compare the thermal process to the case where an additional external electric field is applied, as would be the case in a direct plasma-catalysis setup. We find that the electric field influences the retention time of the molecules at the catalyst surface. These simulations provide an atomic scale insight in the thermal catalytic oxidative dehydrogenation process, and in how an external electric field may affect this process.

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: 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: 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, we investigate the etch process of so-called low-k organic material by means of a N2/H2 capacitively coupled plasma, as applied in the micro-electronics industry for the manufacturing of computer chips. In recent years, such an organic material has emerged as a possible alternative for replacing bulk SiO2 as a dielectric material in the back-end-of-line, because of the smaller parasitic capacity between adjacent conducting lines, and thus a faster propagation of the electrical signals throughout the chip. Numerical simulations with a hybrid plasma model, using an extensive plasma and surface chemistry set, as well as experiments are performed, focusing on the plasma properties as well as the actual etching process, to obtain a better insight into the underlying mechanisms. Furthermore, the effects of gas pressure, applied power and gas composition are investigated to try to optimize the etch process. In general, the plasma density reaches a maximum near the wafer edge due to the so-called 'edge effect'. As a result, the etch rate is not uniform but will also reach its maximum near the wafer edge. The pressure seems not to have a big effect. A higher power increases the etch rate, but the uniformity becomes (slightly) worse. The gas mixing ratio has no significant effect on the etch process, except when a pure H2 or N2 plasma is used, illustrating the synergistic effects of a N2/H2 plasma. In fact, our calculations reveal that the N2/H2 plasma entails an ion-enhanced etch process. The simulation results are in reasonable agreement with the experimental values. The microscopic etch profile shows the desired anisotropic shape under all conditions under study.

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: A fluid model is self-consistently established to investigate the harmonic effects in an inductively coupled plasma, where the electromagnetic field is solved by the finite difference time domain technique. The spatiotemporal distribution of harmonic current density, harmonic potential, and other plasma quantities, such as radio frequency power deposition, plasma density, and electron temperature, have been investigated. Distinct differences in current density have been observed when calculated with and without Lorentz force, which indicates that the nonlinear Lorentz force plays an important role in the harmonic effects, especially at low frequencies. Moreover, the even harmonics are larger than the odd harmonics both in the current density and the potential. Finally, the dependence of various plasma quantities with and without the Lorentz force on various driving frequencies is also examined. It is shown that the deposited power density decreases and the depth of penetration increases slightly because of the Lorentz force. The electron density increases distinctly while the electron temperature remains almost the same when the Lorentz force is taken into account.

Abstract: The formation process of a microdischarge (MD) in both μm- and nm-sized catalyst pores is simulated by a two-dimensional particle-in-cell/Monte Carlo collision model. A parallel-plate dielectric barrier discharge configuration in filamentary mode is considered in ambient air. The discharge is powered by a high voltage pulse. Our calculations reveal that a streamer can penetrate into the surface features of a porous catalyst and MDs can be formed inside both μm- and nm-sized pores, yielding ionization inside the pore. For the μm-sized pores, the ionization mainly occurs inside the pore, while for the nm-sized pores the ionization is strongest near and inside the pore. Thus, enhanced discharges near and inside the mesoporous catalyst are observed. Indeed, the maximum values of the electric field, ionization rate and electron density occur near and inside the pore. The maximum electric field and electron density inside the pore first increase when the pore size rises from 4 nm to 10 nm, and then they decrease for the 100 nm pore, due to
a more pronounced surface discharge for the smaller pores. However, the ionization rate is highest for the 100 nm pore due to the largest effective ionization region.

Abstract: Plasma catalysis is gaining increasing interest for various applications, but the underlying mechanisms are still far from understood. Hence, more fundamental research is needed to understand these mechanisms. This can be obtained by both modelling and experiments. This foundations paper describes the fundamental insights in plasma catalysis, as well as efforts to gain more insights by modelling and experiments. Furthermore, it discusses the state-of-the-art of the major plasma catalysis applications, as well as successes and challenges of technology transfer of these applications.

Abstract: We studied hydrogen sulfide (H2S) decomposition into hydrogen (H2) and sulfur (S2) in a gliding arc plasmatron (GAP) and microwave (MW) plasma by a combination of 0D and 2D models. The conversion, energy efficiency, and plasma distribution are examined for different discharge conditions, and validated with available experiments from literature. Furthermore, a comparison is made between GAP and MW plasma. The GAP operates at atmospheric pressure, while the MW plasma experiments to which comparison is made were performed at reduced pressure. Indeed, the MW discharge region becomes very much contracted near atmospheric pressure, at the conditions under study, as revealed by our 2D model. The models predict that thermal reactions play the most important role in H2S decomposition in both plasma types. The GAP has a higher energy efficiency but lower conversion than the MW plasma at their typical conditions. When compared at the same conversion, the GAP exhibits a higher energy efficiency and lower energy cost than the MW plasma.