Abstract: Research in plasma reactor designs is developing rapidly as plasma technology is gaining increasing interest for sustainable gas conversion applications, like the conversion of greenhouse gases into value-added chemicals and renewable fuels, and fixation of N₂from air into precursors of mineral fertilizer. As plasma is generated by electric power and can easily be switched on/off, these applications allows for efficient conversion and energy storage of intermittent renewable electricity. In this paper, we present a new comprehensive modelling approach for the design and development of gliding arc plasma reactors, which reveals the fluid dynamics, the arc behaviour and the plasma chemistry by solving a unique combination of five complementary models. This results in a complete description of the plasma process, which allows one to efficiently evaluate the performance of a reactor and indicate possible design improvements before actually building it. We demonstrate the capabilities of this method for an experimentally validated study of plasma-based NO_xformation in a rotating gliding arc reactor, which is gaining increasing interest as a flexible, electricity-driven alternative for the Haber–Bosch process. The model demonstrates the importance of the vortex flow and the presence of a recirculation zone in the reactor, as well as the formation of hot spots in the plasma near the cathode pin and the anode wall that are responsible for most of the NO_xformation. The model also reveals the underlying plasma chemistry and the vibrational non-equilibrium that exists due to the fast cooling during each arc rotation. Good agreement with experimental measurements on the studied reactor design proves the predictive capabilities of our modelling approach.

Abstract: Given the ecological problems associated to the CO2 emissions of fossil fuels, plasma technology has gained
interest for conversion of CO2 into value-added products. Microwave plasmas operating at atmospheric pressure
have proven to be especially interesting, due to the high gas temperatures inside the reactor (i.e. up to 6000 K)
allowing for efficient thermal dissociation of CO2 into CO and O2. However, the performance of these high
temperature plasmas is limited by recombination of CO back into CO2 once the gas cools down in the afterglow.
In this work, we computationally investigated several quenching nozzles, developed and experimentally tested
by Hecimovic et al., [1] for their ability to quickly cool the gas after the plasma, thereby quenching the CO
recombination reactions. Using a 3D computational fluid dynamics model and a quasi-1D chemical kinetics
model, we reveal that a reactor without nozzle lacks gas mixing between hot gas in the center and cold gas near
the reactor walls. Especially at low flow rates, where there is an inherent lack of convective cooling due to the
low gas flow velocity, the temperature in the afterglow remains high (between 2000 and 3000 K) for a relatively
long time (in the 0.1 s range). As shown by our quasi-1D chemical kinetics model, this results in a important loss
of CO due to recombination reactions. Attaching a nozzle in the effluent of the reactor induces fast gas quenching
right after the plasma. Indeed, it introduces (i) more convective cooling by forcing cool gas near the walls to mix
with hot gas in the center of the reactor, as well as (ii) more conductive cooling through the water-cooled walls of
the nozzle. Our model shows that gas quenching and the suppression of recombination reactions have more
impact at low flow rates, where recombination is the most limiting factor in the conversion process.

Abstract: 200 years ago, humanity started the industrial revolution by discovering fossil fuels, which lead to unprecedented technological advancements. However it has become alarmingly clear that the major environmental concerns associated with fossil fuels require a short-term transition from a carbon-based energy economy to a sustainable one based on green electricity. A key step concerning this transition exists in developing electricity-driven alternatives for chemical processes that rely on fossil fuels as a raw material. A technology that is gaining increasing interest to achieve this, is plasma technology. Using plasmas to induce chemical reactions by selectively heating electrons in a gas has already delivered promising results for gas conversion applications like CO2 conversion and N2 fixation, but plasma reactors still require optimization to be considered industrially competitive to existing fossil fuel-based processes and emerging other electricity-based technologies. In this thesis I develop computational models to describe plasma reactors and identify key mechanisms in three different plasma reactors for three different gas conversion applications, i.e. N2 fixation, combined CO2-CH4 conversion and CO2 splitting. I first developed models to describe a new rotating gliding arc (GA) reactor operating in two arc modes, which, as revealed by my model, are characterized by distinct plasma chemistry pathways. Subsequently, my colleague and I study the quenching effect of an effusion nozzle to this rotating GA reactor, reaching the best results to date for N2 fixation into NOx at atmospheric pressure, i.e., NOx concentrations up to 5.9%, at an energy cost down to 2.1 MJ/mol. Afterwards, I investigate the possible improvement of N2 admixtures in plasma-based CO2 and CH4 conversion, as significant amounts of N2 are often found in industrial CO2 waste streams, and gas separations are financially costly. Through combining my models with the experiment from a fellow PhD student, we reveal that moderate amounts of N2 (i.e. around 20%) increase both the electron density and the gas temperature to yield an overall energy cost reduction of 21%. Finally, I model quenching nozzles for plasma-based CO2 conversion in a microwave reactor, to explain the enhancements in CO2 conversion that were demonstrated in experiments. Through computational modelling I reveal that the nozzle introduces fast gas quenching resulting in the suppression of recombination reactions, which have more impact at low flow rates, where recombination is the most limiting factor in the conversion process.

Abstract: This work proposes to use core-shell structured spheres to evaluate whether it allows to individually optimize bulk and surface effects of a packing material, in order to optimize conversion and energy efficiency. Different core-shell materials have been prepared by spray coating, using dense spheres (as core) and powders (as shell) of SiO2, Al2O3, and BaTiO3. The materials are investigated for their performance in CO2 dissociation and compared against a benchmark consisting of a packed-bed reactor with the pure dense spheres, as well as an empty reactor. The results in terms of CO2 conversion and energy efficiency show various interactions between the core and shell material, depending on their combination. Al2O3 was found as the best core material under the applied conditions here, followed by BaTiO3 and SiO2, in agreement with their behaviour for the pure spheres. Applying a thin shell layer on the cores showed equal performance between the different shell materials. Increasing the layer thickness shifts this behaviour, and strong combination effects were observed depending on the specific material. Therefore, this method of core-shell spheres has the potential to allow tuning of the packing properties more closely to the application by designing an optimal combination of core and shell.

Abstract: Plasma reactors are interesting for gas-based chemical conversion but the fundamental relation between the plasma chemistry and selected conditions remains poorly understood. Apparent kinetic parameters for the loss and formation processes of individual components of gas conversion processes, can however be extracted by performing experiments in an extended residence time range (2–75 s) and fitting the gas composition to a firstorder kinetic model of the evolution towards partial chemical equilibrium (PCE). We specifically investigated the differences in kinetic characteristics and PCE state of the CO2 dissociation and CH4 reforming reactions in a dielectric barrier discharge reactor (DBD), how these are mutually affected when combining both gases in the dry reforming of methane (DRM) reaction, and how they change when a packing material (non-porous SiO2) is added to the reactor. We find that CO2 dissociation is characterized by a comparatively high reaction rate of 0.120 s−1 compared to CH4 reforming at 0.041 s−1; whereas CH4 reforming reaches higher equilibrium conversions, 82% compared to 53.6% for CO2 dissociation. Combining both feed gases makes the DRM reaction to proceed at a relatively high rate (0.088 s−1), and high conversion (75.4%) compared to CO2 dissociation, through accessing new chemical pathways between the products of CO2 and CH4. The addition of the packing material can also distinctly influence the conversion rate and position of the equilibrium, but its precise effect depends strongly on the gas composition. Comparing different CO2:CH4 ratios reveals the delicate balance of the combined chemistry. CO2 drives the loss reactions in DRM, whereas CH4 in the mixture suppresses back reactions. As a result, our methodology provides some of the insight necessary to systematically tune the conversion process.

Abstract: Experimental evidence in the literature has shown that low-current direct current nitrogen discharges can exist in both glow and arc regimes at atmospheric pressure. However, modelling investigations of the positive column that include the influence of the cathode phenomena are scarce. In this work we developed a 2D axisymmetric model of a plasma discharge in flowing nitrogen gas, studying the influence of the two cathode emission mechanisms—thermionic field emission and secondary electron emission—on the cathode region and the positive column. We show for an inlet gas flow velocity of 1 m s⁻¹in the current range of 80–160 mA, that the electron emission mechanism from the cathode greatly affects the size and temperature of the cathode region, but does not significantly influence the discharge column at atmospheric pressure. We also demonstrate that in the discharge column the electron density balance is local and the electron production and destruction is dominated by volume processes. With increasing flow velocity, the discharge contraction is enhanced due to the increased convective heat loss. The cross sectional area of the conductive region is strongly dependent on the gas velocity and heat conductivity of the gas.

Abstract: Promising biomedical uses for nonthermal plasma (NTP) in the fields of regenerative medicine, cancer therapy, and vaccine delivery involve the noninvasive application of uniform nonequilibrium plasma (including dielectric barrier discharge plasma) to living skin. Whereas most investigations have focused on achieving desired therapeutic outcomes, fewer studies have examined the mechanisms and pathways by which epithelial cells respond to NTP exposure. Using a transwell apical-basolateral-chambered system to culture the human keratinocyte HaCaT cell line, in vitro experiments were performed to demonstrate the effects of nanosecond-pulsed dielectric barrier discharge (nsDBD) plasma on polarized epithelial cell viability, monolayer permeability, intracellular oxidative stress, and the release of adenosine triphosphate (ATP). Application of nsDBD plasma at 60 Hz or below had minimal or no effect on HaCaT monolayer viability or permeability. nsDBD plasma exposure did, however, result in frequency-dependent reductions in intracellular glutathione (indicating direct induction of oxidative stress by nsDBD plasma) and increased extracellular ATP concentrations in the ba-solateral (subepithelial) media, which are indicators of cellular stress and an NTP-induced inflammatory response. These studies provide new insights into nsDBD plasma-induced inflammation and local innate immune responses initiated by polarized epithelial tissues.

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: In this computational study, a gliding arc plasma reactor with a reverse-vortex flow stabilization is modelled for the first time by a fluid plasma description. The plasma reactor operates with argon gas at atmospheric pressure. The gas flow is simulated using the k-ε Reynolds-averaged Navier–Stokes turbulent model. A quasi-neutral fluid plasma model is used for computing the plasma properties. The plasma arc movement in the reactor is observed, and the results for the gas flow, electrical characteristics, plasma density, electron temperature, and gas temperature are analyzed.

Abstract: Cryogenic plasma etching is a promising technique for high-control wafer development with limited plasma induced damage. Cryogenic wafer temperatures effectively reduce surface damage during etching, but the fundamental mechanism is not well understood. In this study, the influences of wafer temperature, gas mixture and substrate bias on the (cryogenic) etch rates of Si with SF6/O2 inductively coupled plasmas are experimentally and computationally investigated. The etch rates are measured in situ with double-point reflectometry and a hybrid computational Monte Carlo – fluid model is applied to calculate plasma properties. This work allows the reader to obtain a better insight in the effects of wafer temperature on the etch rate and to find operating conditions for successful anisotropic (cryo)etching.