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: Plasma (catalysis) reactors are increasingly being used for gas-based chemical conversions, providing an alternative method of energy delivery to the molecules. In this work we explore whether classical concepts such as
equilibrium constants, (overall) rate coefficients, and catalysis exist under plasma conditions. We specifically
investigate the existence of a so-called partial chemical equilibrium (PCE), and how process parameters and
packing properties influence this equilibrium, as well as the overall apparent rate coefficient, for CO2 splitting in
a DBD plasma reactor. The results show that a PCE can be reached, and that the position of the equilibrium, in
combination with the rate coefficient, greatly depends on the reactor parameters and operating conditions (i.e.,
power, pressure, and gap size). A higher power, higher pressure, or smaller gap size enhance both the equilibrium constant and the rate coefficient, although they cannot be independently tuned. Inserting a packing
material (non-porous SiO2 and ZrO2 spheres) in the reactor reveals interesting gap/material effects, where the
type of material dictates the position of the equilibrium and the rate (inhibition) independently. As a result, no
apparent synergistic effect or plasma-catalytic behaviour was observed for the non-porous packing materials
studied in this reaction. Within the investigated parameters, equilibrium conversions were obtained between 23
and 71%, while the rate coefficient varied between 0.027 s−1 and 0.17 s−1. This method of analysis can provide
a more fundamental insight in the overall reaction kinetics of (catalytic) plasma-based gas conversion, in order
to be able to distinguish plasma effects from true catalytic enhancement.

Abstract: DBD plasma reactors are commonly used in a static ‘one inlet – one outlet’ design that goes against reactor design principles for multi-component reactions, such as dry reforming of methane (DRM). Therefore, in this paper we have developed a novel reactor design, and investigated how the shape and size of the reaction zone, as well as gradual gas addition, and the method of mixing CO2 and CH4 can influence the conversion and product com­ position of DRM. Even in the standard ‘one inlet – one outlet’ design, the direction of the gas flow (i.e. short or long path through the reactor, which defines the gas velocity at fixed residence time), as well as the dimensions of the reaction zone and the power delivery to the reactor, largely affect the performance. Using gradual gas addition and separate plasma activation zones for the individual gases give increased conversions within the same operational parameters, by optimising mixing ratios and kinetics. The choice of the main (pre-activated) gas and the direction of gas flow largely affect the conversion and energy cost, while the gas inlet position during separate addition only influences the product distribution.

Abstract: DBD plasma reactors are of great interest for environmental and energy applications, such as CO2 conversion, but they suffer from limited conversion and especially energy efficiency. The introduction of packing materials has been a popular subject of investigation in order to increase the reactor performance. Reducing the discharge gap of the reactor below one millimetre can enhance the plasma performance as well. In this work, we combine both effects and use a packed-bed DBD micro plasma reactor to investigate the influence of gap size reduction, in combination with a packing material, on the conversion and efficiency of CO2 dissociation. Packing materials used in this work were SiO2, ZrO2, and Al2O3 spheres as well as glass wool. The results are compared to a regular size reactor as a benchmark. Reducing the discharge gap can greatly increase the CO2 conversion, although at a lower energy efficiency. Adding a packing material further increases the conversion when keeping a constant residence time, but is greatly dependent on the material composition, gap and sphere size used. Maximum conversions of 50–55% are obtained for very long residence times (30 s and higher) in an empty reactor or with certain packing material combinations, suggesting a balance in CO2 dissociation and recombination reactions. The maximum energy efficiency achieved is 4.3%, but this is for the regular sized reactor at a short residence time (7.5 s). Electrical characterization is performed to reveal some trends in the electrical behaviour of the plasma upon reduction of the discharge gap and addition of a packing material.

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: 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: Plasma is gaining increasing interest for N2 fixation, being a flexible, electricity-driven alternative for the current conventional fossil fuel-based N2 fixation processes. As the vibrational-induced dissociation of N2 is found to be an energy-efficient pathway to acquire atomic N for the fixation processes, plasmas that are in vibrational nonequilibrium seem promising for this application. However, an important challenge in using nonequilibrium plasmas lies in preventing vibrational−translational (VT) relaxation processes, in which vibrational energy crucial for N2 dissociation is lost to gas heating. We present here both experimental and modeling results for the vibrational and gas temperature in a microsecond-pulsed microwave (MW) N2 plasma, showing how power pulsing can suppress this unfavorable VT relaxation and achieve a maximal vibrational nonequilibrium. By means of our kinetic model, we demonstrate that pulsed plasmas take advantage of the long time scale on which VT processes occur, yielding a very pronounced nonequilibrium over the whole N2 vibrational ladder. Additionally, the effect of pulse parameters like the pulse frequency and pulse width are investigated, demonstrating that the advantage of pulsing to inhibit VT relaxation diminishes for high pulse frequencies (around 7000 kHz) and long power pulses (above 400 μs). Nevertheless, all regimes studied here demonstrate a clear vibrational nonequilibrium while only requiring a limited power-on time, and thus, we may conclude that a pulsed plasma seems very interesting for energyefficient vibrational excitation.

Abstract: Plasma medicine has been claimed to provide a novel route to heal wounds and regenerate skin, although very little is currently known about the elementary processes taking place. We carried out a series of ReaxFF-based reactive molecular dynamics simulations to investigate the interaction of O and OH radicals with lipids, more specifically with α-linolenic acid as a model for the free fatty acids present in the upper skin layer. Our calculations predict that the O and OH radicals most typically abstract a H atom from the fatty acids, which can lead to the formation of a conjugated double bond, but also to the incorporation of alcohol or aldehyde groups, thereby increasing the hydrophilic character of the fatty acids and changing the general lipid composition of the skin. Within the limitations of the investigated model, no formation of possibly toxic products was observed.

Abstract: The mechanisms of plasma in medicine are broadly attributed to plasma-derived reactive oxygen and nitrogen species (RONS). In order to exert any intracellular effects, these plasma-derived RONS must first traverse a major barrier in the cell membrane. The cell membrane lipid composition, and thereby the magnitude of this barrier, is highly variable between cells depending on type and state (e.g. it is widely accepted that healthy and cancerous cells have different membrane lipid compositions). In this study, we investigate how plasma-derived RONS interactions with lipid membrane components can potentially be exploited in the future for treatment of diseases. We couple phospholipid vesicle experiments, used as simple cell models, with molecular dynamics (MD) simulations of the lipid membrane to provide new insights into how the interplay between phospholipids and cholesterol may influence the response of healthy and diseased cell membranes to plasma-derived RONS. We focus on the (i) lipid tail saturation degree, (ii) lipid head group type, and (iii) membrane cholesterol fraction. Using encapsulated molecular probes, we study the influence of the above membrane components on the ingress of RONS into the vesicles, and subsequent DNA damage. Our results indicate that all of the above membrane components can enhance or suppress RONS uptake, depending on their relative concentration within the membrane. Further, we show that higher RONS uptake into the vesicles does not always correlate with increased DNA damage, which is attributed to ROS reactivity and lifetime. The MD simulations indicate the multifactorial chemical and physical processes at play, including (i) lipid oxidation, (ii) lipid packing, and (iii) lipid rafts formation. The methods and findings presented here provide a platform of knowledge that could be leveraged in the development of therapies relying on the action of plasma, in which the cell membrane and oxidative stress response in cells is targeted.

Abstract: In recent years, the ability of cold atmospheric pressure plasmas (CAPS) to selectively induce cell death in cancer cells has been widely established. This selectivity has been assigned to the reactive oxygen and nitrogen species (RONS) created in CAPs. To provide new insights in the search for an explanation
for the observed selectivity, we calculate the transfer free energy of multiple ROS across membranes containing a varying amount of cholesterol. The cholesterol fraction is investigated as a selectivity parameter because membranes of cancer cells are known to contain lower fractions of cholesterol compared to healthy cells. We find that cholesterol has a significant effect on the permeation of
reactive species across a membrane. Indeed, depending on the specific reactive species, an increasing cholesterol fraction can lead to (i) an increase of the transfer free energy barrier height and width, (ii) the formation of a local free energy minimum in the center of the membrane and (iii) the creation of extra free energy barriers due to the bulky sterol rings. In the context of plasma oncology, these observations suggest that the increased ingress of RONS in cancer cells can be explained by the decreased cholesterol fraction of their cell membrane.

Abstract: A zero-dimensional, semi-empirical model is used to describe the plasma chemistry in an argon plasma jet flowing into humid air, mimicking the experimental conditions of a setup from the Eindhoven University of Technology. The model provides species density profiles as a function of the position in the plasma jet device and effluent. A reaction chemistry set for an argon/humid air mixture is developed, which considers 84 different species and 1880 reactions. Additionally, we present a reduced chemistry set, useful for higher level computational models. Calculated species density profiles along the plasma jet are shown and the chemical pathways are explained in detail. It is demonstrated that chemically reactive H, N, O and OH radicals are formed in large quantities after the nozzle exit and H2, O2(1Δg), O3, H2O2, NO2, N2O, HNO2 and HNO3 are predominantly formed as 'long living' species. The simulations show that water clustering of positive ions is very important under these conditions. The influence of vibrational excitation on the calculated electron temperature is studied. Finally, the effect of varying gas temperature, flow speed, power density and air humidity on the chemistry is investigated.

Abstract: In this paper we study two cold atmospheric pressure plasma jets, operating in Ar + 2% air, with a different electrode geometry but with the same power dissipated in the plasma. The density profiles of the biomedically active NO and O species throughout the plasma jet, previously obtained by laser diagnostics, are calculated by means of a zero-dimensional semi-empirical reaction kinetics model. A good agreement between the calculated and measured data is demonstrated. Furthermore, the most probable spatial power distribution in an RF driven plasma jet is obtained for the first time by comparing measured and calculated species density profiles. This was possible due to the strong effect of the power distribution on the NO and O density profiles. In addition the dominant reaction pathways for both the NO and the O species are identified. The model allows us to obtain key information on the reactive species production inside the jet, which is difficult to access by laser diagnostics in a coaxial geometry. Finally, we demonstrate that water impurities in the order of 100 ppm in the gas feed can have a significant effect on the spatial distribution of the NO and O density.

Abstract: The use of plasma technology for CO2 splitting is gaining increasing interest, but one of the major obstacles to date for industrial implementation is the considerable energy cost. We demonstrate that the introduction of a packing of dielectric zirconia (ZrO2) beads into a dielectric barrier discharge (DBD) plasma reactor can enhance the CO2 conversion and energy efficiency up to a factor 1.9 and 2.2, respectively, compared to that in a normal (unpacked) DBD reactor. We obtained a maximum conversion of 42 % and a maximum energy efficiency of 9.6 %. However, it is the ability of the packing to almost double both the conversion and the energy efficiency simultaneously at certain input parameters that makes it very promising. The improved conversion and energy efficiency can be explained by the higher values of the local electric field and electron energy near the contact points of the beads and the lower breakdown voltage, demonstrated by 2 D fluid modeling.