Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ar”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Improved sampling technique to collect natural gas from hydrate-bearing pressure cores

High quality gas compositional data are an important factor in interpreting the genetic source of natural gas hosted in hydrate-bearing sediments and other subsurface systems. In order to accurately characterize the composition of gas samples degassed from hydrate-bearing pressure cores, one must use a reproducible sampling technique that minimizes artifacts of the sampling process. Herein, we review sediment core degassing techniques and compare data obtained from a commonly used degassing approach, which we term the standard quantitative degassing (SQD) technique, to our newly developed modified quantitative degassing (MQD) method designed to minimize atmospheric contamination and gas-water interactions. The SQD method allows sample gas to interact with water in a bubbling chamber, which we hypothesize could alter the gas composition following mixing with water or dissolved gases in the bubbling chamber. Whereas, the MQD method allows for the collection of sample gas prior to the bubbling chamber. To compare the SQD and MQD methods, we performed a side-by-side comparison of noble (He, Ne, Ar, Kr, and Xe), major (H 2 , N 2 , O 2 , and CO 2 ), and hydrocarbon (CH 4 , C 2 H 6 , C 3 H 8 , i-C 4 H 10 , C 4 H 10 , i-C 5 H 12 , C 5 H 12 ) gas concentrations and select isotopic compositions obtained using both sample collection techniques. Gas samples were collected from hydrate-bearing pressure cores recovered and maintained under hydrate stable conditions from the northern Gulf of Mexico during the UT-GOM 2 -1 Expedition. The MQD method displayed significantly lower concentrations of atmospheric gases, higher proportions of hydrocarbon gases, lower ratios of C 1 /C 2 + , and heavier stable carbon and hydrogen isotopes of methane than the SQD method. These results demonstrate that the MQD method reduced air contamination and minimized alteration of the hydrocarbon gases. Finally, we conclude this method may be important for future work that seeks to determine the composition of natural gas from pressure cores using quantitative degassing experiments, especially those seeking to measure major (e.g., N 2 ) and noble gases.

58 GEOSCIENCES↗
DOE OSTI2020

Elucidating the Influence of Residual Polymer and Gas Environment on the Electronic Structure of a Graphene Layer Using in Situ APXPS

We use in situ ambient pressure X-ray photoelectron spectroscopy (APXPS) and ultraviolet photoelectron spectroscopy (UPS) to develop an effective method for studying changes in the graphene (Gr) electronic structure according to certain circumstance. The amount of polymethyl methacrylate (PMMA) residual polymer (RP), inevitably generated during the Gr transfer process, is significantly reduced from Gr surface by thermal annealing. This processed Gr is then sequentially exposed to specific gas environments (Ar, N 2 , O 2 , and CO 2 ), and APXPS or UPS is carried out to investigate the variations in the Gr electronic structure including the work function. When the amount of PMMA RP on Gr is reduced, the position of the main carbon peak shifts by more than 0.4 eV to a higher binding energy (in XPS spectra), and the secondary electron cutoff moves by about 0.2 eV to a lower binding energy (in UPS spectra). These changes are generally caused by a decrease in the Gr work function. On the other hand, exposure to the gas environments at different temperatures that we investigated did not produce significant changes in the work function and chemical states of Gr. These results confirm that the material in contact with Gr should be considered to achieve the desired Gr performance in electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Irradiation-induced amorphization of UO 2 films prepared by spraying-assisted combustion synthesis

Spraying-assisted combustion synthesis with uranyl nitrate – acetylacetone – 2-methoxyethanol solutions was used to prepare UO 2 films on an aluminum alloy substrate. The tuning of the spraying parameters and annealing temperatures allowed the preparation of UO 2 films with thicknesses varying from 10–300 nm and 5–10 nm UO 2 grain size. High-resolution electron microscopy and X-ray photoelectron spectroscopy showed that increasing the annealing temperature promotes Mg diffusion from the substrate into the films. The incorporation of Mg reduced the overall crystallinity of the films. The irradiation with Ar 2+ ions (1.7 MeV energy and a fluence of 2 × 10 16 ions/cm 2 ) did not degrade the quality of the films. However, the Mg content significantly influenced the irradiation-induced restructuring of the UO 2 films. Irradiated films with low or no Mg content exhibit high crystallinity, and the UO 2 /Al interfacial layer becomes highly porous. Films with higher Mg content are mostly amorphized after irradiation. The origin of irradiation-induced amorphization was related to the formation of Mg y U 1-y O 2±x solid solutions. Chemically complex, pore-free, and amorphous Mg-Al-O film/substrate interfacial layers enable continuous Mg diffusion during irradiation. As a result, the gradual increase in Mg amounts triggers irradiation-induced precipitation of a crystalline MgO-rich phase within the amorphous films.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2027

Low thermal budget dopant activation in shallow junctions of implanted Si via Tunable plasma enhanced annealing

Low thermal budget (LTB) annealing has become increasingly critical for shallow junctions as semiconductor devices are scaled down. Plasma Enhanced Annealing (PEA) has emerged as a promising LTB annealing process, however its mechanisms remain unclear. In this study, arsenic and boron implanted silicon wafers were subjected to helium or argon ion bombardment generated by annealing plasmas under controlled ion energies and doses in a home-built annealing reactor. The structural and electrical changes were characterized by four-point probe measurements, secondary ion mass spectrometry, Raman spectroscopy and spectroscopic ellipsometry. A pronounced reduction in sheet resistance, accompanied by a decrease in the damage layer thickness demonstrates the surface selective annealing capability of PEA process. Comparative studies using He and Ar plasmas and related ion bombardment reveal a dependence of activation efficacy on both ion species and the dopant implantation profile. Furthermore, a multi-step annealing mechanism is proposed, consisting of an initial (low dose) structural recovery stage followed by the generation of “extended” defects and concluded at large doses by an enhanced dopant activation. The proposed annealing kinetics are thought to be controlled by ion flux, dose and energy. Finally, these results highlight PEA as an effective, surface-selective, LTB approach for shallow-dopant activation and near surface annealing, and provide mechanistic insights into PEA processes.

dopant activation↗
DOE OSTI2020

Water vapor oxidation behaviors of nuclear graphite IG-110 for a postulated accident scenario in high temperature gas-cooled reactors

Water leakage in accidental conditions of high temperature gas-cooled reactors is one of the most critical problems that can compromise the integrity of different nuclear components. In this study, oxidation behaviors of nuclear graphite IG-110 in water ingress accidental conditions were investigated. Mass loss and oxidation rates were evaluated after oxidation tests at temperatures up to 1400 °C in an Ar-20 vol% H 2 O mixed atmosphere. The activation energy decreased from 318.6 to 148.9 kJ/mol with temperature, indicating two different oxidation regimes. The cross-sections of the oxidized samples were systematically characterized. The corresponding logarithmic porosity profiles showed a temperature dependency. Pore formation moved toward near-surface regions with increasing temperature and preferential binder oxidation, with filler particle degradation. Furthermore, oxidant concentration profiles and oxidation depths were estimated using a theoretical model and compared with the experimental results. This work provides important benchmark data and safety analysis guidance for the accident scenario in high temperature gas-cooled reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Effects of transition metals on the evolution of polymer-derived SiOC ceramics

This work focuses on catalytic graphitization from different metals, iron (Fe), cobalt (Co), and nickel (Ni), for inducing carbon growth in silicon oxycarbides (SiOCs) during the pyrolysis. Fe, Co, and Ni-modified SiOCs were synthesized by pyrolysis to 900, 1100, and 1300 °C respectively in Ar. The transition metals induced the formation of the corresponding metal silicides, β-SiC, and graphitic carbon with the catalytic activity in the order of Fe > Co > Ni, in agreement with the activation energy calculation based on the carbon types and amounts. Lateral growth of turbostratic carbon followed a 2D grain growth process and its point-like defect density decreased based on the catalytic order of the transition metals, with SiOC/Fe having the lowest defect density. Quantitative analysis of the XPS results with Rietveld refinement illustrated that the phase separation of SiOC is more dominant than local carbothermal reduction between SiO 2 and C in the SiOC/M (M = Ni, Co, or Fe) systems. Understanding the catalytic graphitization effect of Ni, Co, and Fe on polymer derived ceramics offers new strategies in increasing high-temperature phase amounts and thus creating novel materials for heterogeneous catalysis, magnetic, and other applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A macro-scale ruck and tuck mechanism for deformation in ion-irradiated polycrystalline graphite

A vein structure, which becomes more pronounced with increasing ion dose, was found on the surface of polycrystalline HOPG (highly oriented pyrolytic graphite) implanted by ex situ C + (up to 1.8 × 10 17 ions/cm 2 ), and in situ Ar + in a transmission electron microscope (TEM). These veins are found to be independent of the crystallographic orientations and are associated with the formation of pores. Underneath the veins, a triangular-shaped core was formed with the graphite platelet inside the core displaced up towards the surface. A macro-scale ‘ruck&tuck’ geometry was thus generated at these triangle structure boundaries. Progressive movement of dislocations along basal planes during irradiation was observed, and a mechanistic model was proposed on this basis to explain the vein formation. A small increase of c-spacing was observed with irradiation but it is believed that macro-scale vein formation plays a more vital role in the dimensional and property changes in polycrystalline graphite, especially when a stress gradient is present. Here, the model proposed also explains the change of thermal expansion in HOPG with irradiation. Together with Heggie’s ‘ruck&tuck’ and Barsoum’s ‘ripplocations’ models, the present model is considered to have provided an additional experimentally proven mechanism responsible for irradiation behaviour in graphite materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nanoparticle dynamics in the spatial afterglows of nonthermal plasma synthesis reactors

Nonthermal plasma flow tube reactors are industrially scalable systems for the production of nanocrystal (NC) based materials and coatings. One key advantage of nonthermal plasma synthesis is the ability to both synthesize NCs and deposit films in a single reactor, as at the reactor outlet, NCs can be inertially deposited onto a target substrate. The size and morphology of deposited particles can substantially influence the film structure and function. Though NCs are typically near-spherical and monodispersed as-produced in plasma synthesis reactors, NC charge and growth dynamics can be altered substantially when NCs are sampled out of the plasma and through the spatial afterglow region, affecting deposition. Experiments have demonstrated changes of NC size and charge in the spatial afterglow; however, these dynamics remain unexplored and unexplained via theory and simulation. To address this, we developed a constant number Monte Carlo (CNMC) simulation model to examine the mechanisms of NC decharging and growth in the spatial afterglow of plasma flow tube reactors. Collisions between NC and plasma species, diffusive deposition, and electron desorption from NCs are incorporated in the CNMC simulation. The simulation results are specifically compared with previous experiments on Si NCs synthesized from a low pressure Ar-SiH 4 nonthermal plasma reactor. Furthermore, the experiment-model comparison shows that CNMC models can be implemented which accurately model NC size distribution evolution in a spatial afterglow. Simultaneously, results show that improved collision models, energetic species diffusion models, and electron desorption models will be necessary to accurately depict NC dynamics in spatial afterglows.

42 ENGINEERING↗
DOE OSTI2022

Quantifying the anion effect of gas solubility within ionic liquids using the solvation affinity index

The solubility of gas molecules within ionic liquids (ILs) is strongly dictated by the anion effect and the free volume effect, and of these two influences, the anion effect can be challenging to predict. We evaluate a wide range of different gas solutes (O 2 , NH 3 , CH 4 , H 2 S, SO 2 , Ar, H 2 , N 2 , CO 2 ) in different ILs, and we find that the solubility is strongly correlated to one of three different IL properties, which can be directly calculated from the anion electrostatic potential (ESP) surfaces. Thus, the solutes can be grouped into three fundamental types, leading to a general descriptor, the solvation affinity index (SAI), which can identify the affinity between a wide range of gas solutes and IL solvents. The SAI of nine common industrial gases with 47 common anions are generated. Furthermore, this can provide a fast approach for pre-screening IL solvents or even for ILs embedded within porous supports.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Position-specific isotope effects during alkaline hydrolysis of 2,4-dinitroanisole resolved by compound-specific isotope analysis, 13 C NMR, and density-functional theory

Compound-specific isotope analysis (CSIA), position-specific isotope analysis (PSIA), and computational modeling (e.g., quantum mechanical models; reactive-transport models) are increasingly being used to monitor and predict biotic and abiotic transformations of organic contaminants in the field. However, identifying the isotope effect(s) associated with a specific transformation remains challenging in many cases. We describe and interpret the position-specific isotope effects of C and N associated with a SN 2 Ar reaction mechanism by a combination of CSIA and PSIA using quantitative 13 C nuclear magnetic resonance spectrometry, and density-functional theory, using 2,4-dinitroanisole (DNAN) as a model compound. The position-specific 13 C enrichment factor of O–C 1 bond at the methoxy group attachment site (ε C1 ) was found to be approximately -41‰, a diagnostic value for transformation of DNAN to its reaction products 2,4-dinitrophenol and methanol. Theoretical kinetic isotope effects calculated for DNAN isotopologues agreed well with the position-specific isotope effects measured by CSIA and PSIA. This combination of measurements and theoretical predictions demonstrates a useful tool for evaluating degradation efficiencies and/or mechanisms of organic contaminants and may promote new and improved applications of isotope analysis in laboratory and field investigations.

13C NMR↗
DOE OSTI2020

Spherically expanding flame in silane–hydrogen–nitrous oxide–argon mixtures

The effect of silane addition on the laminar flame speed ($S^0_u$) of flames propagating in hydrogen–nitrous oxide–argon mixtures has been investigated experimentally for the first time using the spherically expanding flame technique in a constant volume combustion chamber. Replacing hydrogen by silane and maintaining the equivalence ratio constant, much higher flame speeds, explosion peak pressures, and pressure rise coefficients were measured. A previously developed detailed reaction model has been updated based on ab initio thermodynamic properties calculations and collision limit violation analysis. The improved reaction model demonstrates encouraging performance in predicting the flame speed, with a mean absolute error below 11%. To explain the effect of silane addition on the flame dynamics, a number of parameters have been calculated including OH and H rate of production, heat release rate per reaction, and sensitivity coefficient on . The dynamics of freely propagating flames in SiH 4 –H 2 –N 2 O–Ar mixtures is essentially controlled by reactions of the H–O–N chemical system: N 2 O+H=N 2 +OH, OH+H 2 =H 2 O+H, and N 2 O(+M)=N 2 +O(+M). Whereas silane addition does not influence much the rate of production of OH, it significantly modifies that of H with a number of pyrolytic chemical pathways of silicon hydrides, such as SiH+H 2 =SiH 2 +H and Si+H 2 =SiH+H, which act as sink of H atom as they proceed in the backward direction. The reactions forming SiO(s) and SiO 2 (s), such as SiO+OH=SiO 2 (s)+H and 2SiO=2SiO(s), are exothermic and significantly contribute to the temperature increase. The adiabatic, constant pressure flame temperature for mixture containing silane is significantly higher, up to several 100’s K. The increase of induced by silane addition seems to be mostly related to the large increase of the flame temperature which leads to higher energy release rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Time-domain modelling and thermometry of the CH 4 v 1 Q-branch using hybrid femtosecond/picosecond coherent anti-Stokes Raman scattering

Here, we present time-domain measurements and modelling of the CH 4 v 1 vibrational Q-branch, as well as simultaneous temperature acquisition from the CH 4 and N 2 vibrational Q-branches using broadband hybrid femtosecond/picosecond coherent anti-Stokes Raman scattering (fs/ps CARS). Accurate measurements with fs/ps CARS require knowledge of the time-domain response of the probed molecule at different temperatures and pressures. In this work, a fs/ps CARS model was developed for theCH 4 v 1 Q-branch and validated using time-resolved probe delay scans at pressures between 70 and 600 Torr and temperatures between 295 and 1000 K. A simplified modified exponential energy gap model was used to determine the linewidths of the Raman transitions within the Q-branch. Total Q-branch linewidths for the entire explored parameter space were also determined from the probe delay scans. This enabled quantitative measurements of CH 4 with fs/ps CARS as well as the temperature through the hot bands of the CH 4 v 1 Q-branch. With a broadband fs laser pulse in a two-beam phase matching scheme, the N 2 Q-branch and CH 4 v 1 Q-branch were simultaneously measured in a heated cell, and the temperatures from each molecule's spectra was evaluated. Here, the fitted temperatures corresponded well with differences ranging from 4 to 17%. It was also shown that exchanging 90% of CH 4 for N 2 or Ar at 295 K and 500 Torr did not significantly impact the normalized time-domain signal prior to a 100 picosecond probe delay. Spatially resolved one-dimensional (1-D) thermometry using the CH 4 Q-branch was also demonstrated with an imaged spatial resolution of 40 µm. With the developed model, 1-D quantitative measurements of CH 4 concentration and temperature can be measured with time resolution of tens of picoseconds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Time-resolved measurements of OH during auto-ignition of syngas with trimethylsilanol and hexamethyldisiloxane

The effects of trimethylsilanol (TMSO) and hexamethyldisiloxane (HMDSO) addition on OH time histories during syngas (H 2 and CO) ignition were investigated using the University of Michigan rapid compression facility. Experiments spanned temperatures of 1010–1080 K, at a pressure of approximately 5 atm. Syngas mixtures of 1.2 % H 2 /2.8 % CO/20 % O 2 by volume (balance N 2 and Ar) provided a baseline for comparison with mixtures that included 100, 200, and 1000 ppm of the TMSO and 100 ppm of HMDSO. Narrow-line ultraviolet laser-absorption was used to measure OH mole-fraction during ignition. The addition of TMSO and HMDSO significantly shifted the OH time-histories earlier in time, by up to 51 %, compared with the baseline syngas mixture. The value of the maximum OH mole fraction was consistent between the 100 and 200 ppm TMSO mixtures and the 100 ppm HMDSO mixtures, but the maximum OH increased significantly with the 1000 ppm TMSO mixtures. Here, the OH data indicate TMSO and HMDSO were not direct sources of OH radicals. Analysis further indicates the TMSO and HMDSO decompose rapidly followed by reactions that enhance the production of H atoms, and the increased reactivity observed is via the H + O 2 = OH + O reaction.

Hexamethyldisiloxane↗
DOE OSTI2020

High temperature oxidation behaviors of bulk SiC with low partial pressures of air and water vapor in argon

Oxidation behaviors of bulk SiC at 1000–1400 °C in Ar-20 vol% H 2 O-10 vol% air were investigated. The thickness of the passivation oxide layer showed a parabolic increase with time. Dense and uniform layers formed below 1200 °C, whereas pores and cracks were observed at higher temperatures. An analytical model was proposed to estimate the oxide layer growth and SiC mass loss rate. The SiC layer should be effective in retaining the fission products after oxidation by excessive water and air. This work is the first detailed study of SiC layer oxidation during loss of coolant accidents (LOCA) for HTGRs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Protective role of silicon in the high-temperature oxidation of titanium

We investigated the oxidation of pure titanium and Ti-0.8 at.% Si at 800 ° C in Ar-O 2 and O 2 to clarify how dilute silicon additions enhance the oxidation resistance of titanium. The effect of silicon stemmed from differences in oxide growth rather than oxygen dissolution in the metal. Initially, silicon reduced outward oxide growth and completely eliminated inward growth, which we attributed to the continuous Ti 5 Si 3 layer formed at the oxide-metal interface. Furthermore, after breakdown of the Ti 5 Si 3 layer, the oxide — containing silicon — became lamellar but remained more compact and possibly slower growing than the oxide on pure titanium.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Role of niobium in the high-temperature oxidation of titanium

Here, we compared the oxidation of Ti and Ti-2 at.% Nb at 800°C in Ar-20% O 2 and in sequential 16 O 2 - 18 O 2 atmospheres to clarify the role of niobium in the high-temperature oxidation of titanium. Niobium had no effect on the amount of oxygen dissolved in the metal or on oxygen diffusion in the metal. Instead, niobium enhanced oxidation resistance by segregating to TiO 2 grain boundaries, where it likely slowed outward Ti diffusion. These findings spotlight doping of TiO 2 grain boundaries as a strategy for improving the oxidation resistance of Ti alloys and, more generally, TiO 2 -forming alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

PotLib 2023: New version of a potential energy surface library for chemical systems

POTLIB is a library of global and semiglobal potential energy surface subprograms. The library currently features 410 entries, including both single-state entries and multi-state entries. When one calls the routine of a single-state entry, it returns the ground-electronic-state adiabatic potential energy surface at the input geometry. In addition, some entries also return the gradient of the surface. Further, when one calls a multi-state entry, it returns a diabatic potential energy matrix (DPEM). If the entry also has the gradient of the DPEM, one can compute adiabatic surfaces, their gradients, and the nonadiabatic coupling vectors (NACs) from the DPEM and its gradient by diagonalization. Some but not all the routines conform to one of a set of standard interfaces. The goal is to facilitate chemical dynamics research by collecting and disseminating a comprehensive collection of state-of-the-art potential energy routines (developed by a wide, international group of researchers) with systematic and well-defined interfaces for use with chemical dynamics programs. Systems in the library include CHArO 2 , CHN 2 O + , CH 2 O, CH 2 O 2 , CH 3 N 2 , CH 3 O, CH 4 , CH 4 Br, CH 4 Cl, CH 4 F, CH 4 O, CH 4 OCl, CH 4 OF, CH 4 OH, CH 5 , CH$^{+}_{5}$, CH 5 N, CH 5 O 2 , CH 8 O 2 , C 2 H 2 N 2 O, C 2 H 2 O, C 2 H 4 N, C 2 H 4 O 4 , C 2 H 6 Cl, C 2 H 6 F, C 2 H 6 H, C 2 H 6 O, C 2 H 6 OH, C 2 H 7 , C 2 O 2 , C 3 H 4 O 2 , C 3 H 7 NO, C 6 H 6 O, C 6 H 6 S, C 7 H8S, HBrCl, HCl 2 , HF 2 , HI 2 , HLiF, HNaF, HO 2 , HO 3 , HOBr, HSiO, H 2 Br, H 2 ClO, H 2 F, H 2 F 2 , H 2 FO, H 2 Na, H 2 O, H 2 O 2 , H 2 OBr, H 3 , H 3 Cl, H 3 ClN, H 3 ClO, H 3 N, H 3 O, H 3 O 2 , H$^{+}_{3}$, H 3 S, H 4 ClSi, H 4 N, H 4 NO, H 4 O 2 , H 5 GeO, H 5 Si, H$^{+}_{7}$, Al m H n , Al n , Ar N O, K 2 Rb 2 , NO 2 , N 2 O, N 2 O 2 , N 3 , N 4 , O 3 , and O 4 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Vidyut3d: A GPU accelerated fluid solver for non-equilibrium plasmas on adaptive grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure three-electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate ~ 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗