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: 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.

Abstract: Plasma catalysis is promising for various environmental, energy and chemical synthesis applications, but the underlying mechanisms are far from understood. Modeling can help to obtain a better insight in these mechanisms. Some burning questions relate to the plasma behavior inside packed bed reactors and whether plasma can penetrate into catalyst pores. In this paper, we try to provide answers to these questions, by means of both fluid modeling and particle-in-cell/Monte Carlo collision simulations. We present a short overview of recent findings obtained in our group by means of modeling, i.e., the enhanced electric field near the contact points and the streamer propagation through the packing in packed bed reactors, as well as the plasma behavior in catalyst pores, to determine the minimum pore size in which plasma streamers can penetrate.

Abstract: The laser-modulated magnetopolaron effect on shallow donors in semiconductors is investigated in the presence of a magnetic field in the Voigt configuration. A nonperturbative approach is used to describe the electron-photon interaction by including the radiation field in an exact way via a laser-dressed interaction potential. Through a variational approach we evaluate the donor binding energy. We find that the interaction strength of the laser-dressed Coulomb potential in the z direction cannot only be enhanced but also weakened by the radiation field, while that in the x-y plane is only weakened. In this way, the binding energy of the states with odd z parity, like 2p(z) can be decreased or increased with respect to its static binding energy by the radiation field, while that of the other states can be only decreased. Furthermore, all binding energies become insensitive to the magnetic field if the radiation field is strong. The magnetopolaron effect on these energies is studied within second-order time-dependent perturbation theory. In the nonresonant region, a laser-modulated magnetopolaron correction, including the effect of single-photon processes, is observed. In the resonant region, a laser-modulated magnetopolaron effect, accompanied by the emission and absorption of a single photon, is found. Moreover, the 1s -> 2p(+) transition, accompanied by the emission of a single photon, is tuned by the radiation field into resonance with the longitudinal-optical phonon branch. This is electrically analogous to the magnetopolaron effect, and therefore we name it the dynamical magnetopolaron effect. Finally, by changing the frequency of the radiation field, these interesting effects can be tuned to be far away from the reststrahlen band and, therefore, can be detected experimentally. This in turn provides a direct measure of the electron-phonon interaction.

Abstract: This chapter discusses modeling efforts for plasma-based CO2 conversion, which are needed to obtain better insight in the underlying mechanisms, in order to improve this application. We will discuss two types of (complementary) modeling efforts that are most relevant, that is, (i) modeling of the detailed plasma chemistry by zero-dimensional (0D) chemical kinetic models and (ii) modeling of reactor design, by 2D or 3D fluid dynamics models. By showing some characteristic calculation results of both models, for CO2 splitting and in combination with a H-source, and for packed bed DBD and gliding arc plasma, we can illustrate the type of information they can provide.

Abstract: Nonthermal plasma is a promising alternative for ammonia synthesis at gentle conditions. Metal meshes of Fe, Cu, Pd, Ag, and Au were employed as catalysts in radio frequency plasma for ammonia synthesis. The energy yield for all these transition metal catalysts ranged between 0.12 and 0.19 g-NH3/kWh at 300 W and, thus, needs further improvement. In addition, a semimetal, pure gallium, was used for the first time as catalyst for ammonia synthesis, with energy yield of 0.22 g-NH3/kWh and with a maximum yield of ∼10% at 150 W. The emission spectra, as well as computer simulations, revealed hydrogen recombination as a primary governing parameter, which depends on the concentration or flux of H atoms in the plasma and on the catalyst surface. The simulations helped to elucidate the underlying mechanism, implicating the dominance of surface reactions and surface adsorbed species. The rate limiting step appears to be NH2 formation on the surface of the reactor wall and on the catalyst surface, which is different from classical catalysis.

Abstract: Landau levels (LLs) are modified by the Frohlich interaction which we investigate within the improved Wigner-Brillouin theory for energies both below and above the longitudinal-optical-continuum in monolayer MoS2.., WS2, MoSe2, and WSe2. Polaron corrections to the LLs are enhanced in monolayer MoS2 as compared to WS2. A series of levels are found at h omega(LO) + lh omega(c), and in addition, the Frohlich interaction lifts the degeneracy between the levels nh omega(c) and h omega(LO) + lh omega(c) resulting in an anticrossing. The screening effect due to the environment plays an important role in the polaron energy corrections, which are also affected by the effective thickness r(eff) parameter. The polaron anticrossing energy gap E-gap decreases with increasing effective thickness r(eff). Published by AIP Publishing.

Abstract: One of the main issues in plasma chemistry modeling is that the cross sections and rate coefficients are subject to uncertainties, which yields uncertainties in the modeling results and hence hinders the predictive capabilities. In this paper, we reveal the impact of these uncertainties on the model predictions of plasma-based dry reforming in a dielectric barrier discharge. For this purpose, we performed a detailed uncertainty analysis and sensitivity study. 2000 different combinations of rate coefficients, based on the uncertainty from a log-normal distribution, are used to predict the uncertainties in the model output. The uncertainties in the electron density and electron temperature are around 11% and 8% at the maximum of the power deposition for a 70% confidence level. Still, this can have a major effect on the electron impact rates and hence on the calculated conversions of CO2 and CH4, as well as on the selectivities of CO and H2. For the CO2 and CH4 conversion, we obtain uncertainties of 24% and 33%, respectively. For the CO and H2 selectivity, the corresponding uncertainties are 28% and 14%, respectively. We also identify which reactions contribute most to the uncertainty in the model predictions. In order to improve the accuracy and reliability of plasma chemistry models, we recommend using only verified rate coefficients, and we point out the need for dedicated verification experiments.

Abstract: In this work, a chemical kinetics study on methane activation for hydrogen production in a warm plasma, i.e., N-2 rotating gliding arc (RGA), was performed for the first time to get new insights into the underlying reaction mechanisms and pathways. A zero-dimensional chemical kinetics model was developed, which showed a good agreement with the experimental results in terms of the conversion of CH4 and product selectivities, allowing us to get a better understanding of the relative significance of various important species and their related reactions to the formation and loss of CH4, H-2, and C2H2 etc. An overall reaction scheme was obtained to provide a realistic picture of the plasma chemistry. The results reveal that the electrons and excited nitrogen species (mainly N-2(A)) play a dominant role in the initial dissociation of CH4. However, the H atom induced reaction CH4+ H -> CH3+ H-2, which has an enhanced reaction rate due to the high gas temperature (over 1200 K), is the major contributor to both the conversion of CH4 and H-2 production, with its relative contributions of > 90% and > 85%, respectively, when only considering the forward reactions. The coexistence and interaction of thermochemical and plasma chemical processes in the rotating gliding arc warm plasma significantly enhance the process performance. The formation of C-2 hydrocarbons follows a nearly one-way path of C2H6 -> C2H4 -> C2H2, explaining why the selectivities of C-2 products decreased in the order of C2H2 > C2H4 > C2H6.

Abstract: A packed bed dielectric barrier discharge (DBD) is widely used for plasma catalysis applications, but the exact plasma characteristics in between the packing beads are far from understood. Therefore, we study here these plasma characteristics by means of fluid modelling and experimental observations using ICCD imaging, for packing materials with different dielectric constants. Our study reveals that a packed bed DBD reactor in dry air at atmospheric pressure may show three types of discharges, i.e. positive restrikes, filamentary microdischarges, which can also be localized between two packing beads, and surface discharges (so-called surface ionization
waves). Restrikes between the dielectric surfaces result in the formation of filamentary microdischarges, while surface charging creates electric field components parallel to the dielectric surfaces, leading to the formation of surface discharges. A transition in discharge mode occurs from surface discharges to local filamentary discharges between the packing beads when the dielectric constant of the packing rises from 5 to 1000. This may have implications for the efficiency of plasma catalytic gas treatment, because the catalyst activation may be limited by constraining the discharge to the contact points of the beads. The production of reactive species occurs most in the positive restrikes, the surface discharges and the local microdischarges in between the beads, and is less significant in the longer filamentary microdischarges. The faster streamer propagation and discharge development with higher dielectric constant of the packing beads leads to a faster production of reactive species. This study is of great interest for plasma catalysis, where packing beads with different dielectric constants are often used as supports for the catalytic materials. It allows us to better understand how different packing materials can influence the performance of packed bed plasma reactors for environmental applications.

Abstract: The gliding arc plasmatron (GAP) is a highly efficient atmospheric plasma source, which is very promising for CO2 conversion applications. To understand its operation principles and to improve its application, we present here comprehensive modeling results, obtained by means of computational fluid dynamics simulations and plasma modeling. Because of the complexity of the CO2 plasma, a full 3D plasma model would be computationally impractical. Therefore, we combine a 3D turbulent gas flow model with a 2D plasma and gas heating model in order to calculate the plasma parameters and CO2 conversion characteristics. In addition, a complete 3D gas flow and plasma model with simplified argon chemistry is used to evaluate the gliding arc evolution in space and time. The calculated values are compared with experimental data from literature as much as possible in order to validate the model. The insights obtained in this study are very helpful for improving the application of CO2 conversion, as they allow us to identify the limiting factors in the performance, based on which solutions can be provided on how to further improve the capabilities of CO2 conversion in the GAP.

Abstract: Ablation-controlled plasmas have been used in a range of technical applications where local thermodynamic equilibrium (LTE) is often violated near the wall due to the strong cooling effect caused by the ablation of wall materials. The thermodynamic and transport properties of ablated polytetrafluoroethylene (PTFE) vapor, which determine the flowing plasma behavior in such applications, are calculated based on a two-temperature model at atmospheric pressure. To our knowledge, no data for PTFE have been reported in the literature. The species composition and thermodynamic properties are numerically determined using the two-temperature Saha equation and the Guldberg-Waage equation according to van de Sanden et al's derivation. The transport coefficients, including viscosity, thermal conductivity and electrical conductivity, are calculated with the most recent collision interaction potentials using Devoto's electron and heavy-particle decoupling approach but expanded to the third-order approximation (second-order for viscosity) in the frame of the Chapman-Enskog method. Results are computed for different degrees of thermal non-equilibrium, i.e. the ratio of electron to heavy-particle temperatures, from 1 to 10, with electron temperature ranging from 300 to 40 000 K. Plasma transport properties in the LTE state obtained from the present work are compared with existing published results and the causes for the discrepancy analyzed. The two-temperature plasma properties calculated in the present work enable the modeling of wall ablation-controlled plasma processes.

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: The conversion of atmospheric nitrogen into valuable compounds, that is, so-called nitrogen fixation, is gaining increased interest, owing to the essential role in the nitrogen cycle of the biosphere. Plasma technology, and more specifically gliding arc plasma, has great potential in this area, but little is known about the underlying mechanisms. Therefore, we developed a detailed chemical kinetics model for a pulsed-power gliding-arc reactor operating at atmospheric pressure for nitrogen oxide synthesis. Experiments are performed to validate the model and reasonable agreement is reached between the calculated and measured NO and NO2 yields and the corresponding energy efficiency for NOx formation for different N2/O2 ratios, indicating that the model can provide a realistic picture of the plasma chemistry. Therefore, we can use the model to investigate the reaction pathways for the formation and loss of NOx. The results indicate that vibrational excitation of N2 in the gliding arc contributes significantly to activating the N2 molecules, and leads to an energy efficient way of NOx production, compared to the thermal process. Based on the underlying chemistry, the model allows us to propose solutions on how to further improve the NOx formation by gliding arc technology. Although the energy efficiency of the gliding-arc-based nitrogen fixation process at the present stage is not comparable to the world-scale HaberBosch process, we believe our study helps us to come up with more realistic scenarios of entering a cutting-edge innovation in new business cases for the decentralised production of fertilisers for agriculture, in which low-temperature plasma technology might play an important role.

Abstract: The discharge properties of the plasma bulk flow near the surface of heated compound-materials strongly affects the kinetic layer parameters modeled and manifested in the Knudsen layer. This paper extends the widely used two-layer kinetic ablation model to the ablation controlled non-equilibrium discharge due to the fact that the local thermodynamic equilibrium (LTE) approximation is often violated as a result of the interaction between the plasma and solid walls. Modifications to the governing set of equations, to account for this effect, are derived and presented by assuming that the temperature of the electrons deviates from that of the heavy particles. The ablation characteristics of one typical material, polytetrafluoroethylene (PTFE) are calculated with this improved model. The internal degrees of freedom as well as the average particle mass and specific heat ratio of the polyatomic vapor, which strongly depends on the temperature, pressure and plasma non-equilibrium degree and plays a crucial role in the accurate determination of the ablation behavior by this model, are also taken into account. Our assessment showed the significance of including such modifications related to the non-equilibrium effect in the study of vaporization of heated compound materials in ablation controlled arcs. Additionally, a two-temperature magneto-hydrodynamic (MHD) model accounting for the thermal non-equilibrium occurring near the wall surface is developed and applied into an ablation-dominated discharge for an electro-thermal chemical launch device. Special attention is paid to the interaction between the non-equilibrium plasma and the solid propellant surface. Both the mass exchange process caused by the wall ablation and plasma species deposition as well as the associated momentum and energy exchange processes are taken into account. A detailed comparison of the results of the non-equilibrium model with those of an equilibrium model is presented. The non-equilibrium results show a non-equilibrium region near the plasma-wall interaction region and this indicates the need for the consideration of the influence of the possible departure from LTE in the plasma bulk on the determination of ablation rate.