Data set of fully kinetic simulations of capacitively-coupled radio-frequency plasma discharges in argon
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Arsenic contamination in groundwater threatens public health, particularly in small, low-income communities lacking affordable treatment solutions. This study investigated the field implementation of novel air cathode assisted iron electrocoagulation (ACAIE) technology for arsenic removal in Allensworth, California, where groundwater arsenic concentrations exceeded 250 µg/L. Over four months, a pilot-scale ACAIE system, operating at 600 L/h, consistently reduced arsenic levels to below the EPA’s maximum contaminant level of 10 µg/L. Laboratory experiments informed the optimization of charge dosage and flow rates, which were validated during field testing of the ACAIE 600 L/h system. The in-situ generation of hydrogen peroxide at the cathode speeded up the reaction kinetics, ensuring high arsenic removal efficiency while allowing high throughput, even with a compact reactor size. An economic analysis demonstrated a treatment cost of USD 0.02/L excluding labor, highlighting the system’s affordability compared to conventional methods. Adding labor costs increased the treatment cost to USD 0.09/L. The regeneration of air cathodes extended their operational life, addressing a key maintenance challenge, thus reducing the costs slightly. Intermittent challenges were encountered with filtration and secondary contaminant removal; these issues highlight opportunities for further operational improvements. Despite these challenges, ACAIE’s low operational complexity, scalability, and cost-effectiveness make it a promising solution for underserved small communities. These findings provide critical insights into deploying sustainable arsenic remediation technologies that are tailored to the needs of rural, low-resource communities.
Explore the source record for details and available documents.
Electrification and use of renewable hydrogen is currently a necessity for decarbonizing the iron-and-steel industry. In this regard, hydrogen plasma smelting reduction (HPSR) is a novel pathway that is being explored for reduction of iron ore. HPSR provides several decarbonization merits compared to conventional blast furnaces. Firstly, the use of renewable hydrogen drastically reduces the CO2 emissions compared to the use of coke. Secondly, renewable electricity in the form of a thermal plasma for making reactive hydrogen species (radicals, ions) are more efficient at reducing iron ore compared to neutral H2. Thirdly, a molten product compatible with downstream processes is obtained from the intense heat transfer from the plasma. However, the scale-up of this technology requires fundamental exploration of hydrogen plasma dynamics and its interaction with complex solid material that include phase changing iron-ore and slag. In this work, we present a first principles continuum scale model for thermal plasmas in Ar/H2 gas mixtures typically used for HPSR. The thermal plasma governing equations for mass, momentum and energy with Lorentz force and Joule heating source terms are solved along with electromagnetic equations for electrostatic and magnetic vector potential. Our solver will be based on Pele, a suite of reacting flow solvers designed for advanced scientific computing architectures (Henry De Frahan et al., Proceedings of SIAM Parallel Processing, 13-25, 2024), and will utilize adaptive mesh generation for enhanced resolutions at locations of intense physicochemical interactions. This study will present the impact of Ar to H2 ratios on excited/dissociated hydrogen species concentrations, plasma temperature and conductivity along with the impact of outgassed species (water, metal vapor, O, OH radicals) from ore surface on gas phase chemistry. Furthermore, the heat and species flux to the surface will be quantified as a function of applied voltages in a transferred arc configuration.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
A digital-to-analog converter circuit that creates an analog waveform from an input digital waveform. Operating the circuit comprises using the input digital waveform to 1) operate a charge control switch to set a charge time period, 2) operate a discharge control switch to set a discharge time period, 3) set a charge current magnitude using a charge gain, and 4) set a discharge current magnitude using a discharge gain. A charge source electrically charges a load capacitor during the charge time period (i.e., the charge mode). A discharge source electrically discharges the load capacitor during the discharge time period (i.e., the discharge mode). A circuit output transmits the analog waveform defined by the charge mode and the discharge mode. A charge current magnitude greater than the discharge current magnitude produces an upward-sloping analog waveform. A charge current magnitude less than the discharge current magnitude produces a downward-sloping analog waveform.
The frequency of type-I ELMs decreases as the initiation of the neutral beam injection (NBI) heating is delayed with respect to the time when plasma current (I p ) reaches flat-top in the ITER Baseline Scenario discharges in DIII-D. Henceforth, the time gap between the NBI initiation and I p flat-top will be referred to as “heating delay.” As the heating delay is modified, pedestal formation follows different trajectories in the edge current density–pedestal pressure gradient (j edge -∇p e ped ) space from the L-H transition toward the first ELM event. During the stationary phase after the first ELM, the ELM frequency (f ELM ) decreases by a factor of ~2 as the heating delay is increased. A longer pedestal recovery time in the inter-ELM period is observed for the low f ELM discharges as compared to the high f ELM discharges. Both low and high f ELM discharges show nearly identical profiles of electron density and temperature and have a similar MHD stability just before an ELM crash. However, a marked difference is observed in the magnetic spectrogram of the high and low f ELM discharges in response to the variation in the heating delay. The main difference is in the 200–400 kHz range of the magnetic spectra. A quasi-coherent mode (QCM) at 220 kHz and weaker broadband fluctuations are observed in the high f ELM discharges, while only strong broadband fluctuations are prevalent in the low f ELM discharges. ELM-synchronized analysis shows that the time evolution of these modes is different for the high and low f ELM discharges. The localization of both these modes is confirmed at the maximum gradient region of the pedestal. We hypothesize that these modes cause important pedestal transport and that the difference in the pedestal recovery of the high and low f ELM discharges is a result of the difference in transport driven by these modes, as they change with changes in the heating delay. It is demonstrated experimentally for the first time that discharges with similar pedestal parameters can carry the history of the heating delay into the stationary phase and that changes in turbulent-driven transport are a likely cause of changes in f ELM observed with variations of heating delay.
Abstract The present work improves a phenomenological plasma-assisted combustion model by integrating the spatiotemporal distribution of plasma power density, thereby considering the evolution of plasma streamers in the modeling, and subsequently, better predicting the ignition kernel evolution. The improved phenomenological model is validated against experiments representing the plasma discharge and post-discharge ignition kernel evolution. Specifically, the new model demonstrates a more accurate prediction of ultrafast gas heating and O 2 dissociation during the plasma discharge, compared to the original model. In addition, the new model is found to closely match the experimental pressure wave and heated channel profiles post-discharge without the need for tuning the energy deposition (unlike the original model), highlighting its accuracy of post-discharge ignition kernel dynamics. The improved phenomenological model is then employed to investigate ignition kernel evolution for a stoichiometric methane-air discharge across various discharge gap configurations. Simulations reveal a non-uniform temperature and streamer distribution progressing from the electrode tips toward the center, contrasting uniform cylindrical discharges previously described in the original model. Streamer propagation is observed to be faster for larger gaps when maintained at the same average electric field for different discharge gaps. The tendency of smaller gaps to produce detached toroidal ignition kernels is observed, while larger gaps promote cylindrical and attached ignition kernels. Interactions between successive ignition kernels from consecutive discharges varied significantly, with the smallest gap (1 mm) promoting the quenching of the preceding ignition kernel due to the initial kernel–kernel separation. The intermediate gap (2 mm) promotes detached kernel growth. In contrast, in the largest gap (4 mm), kernels consistently combine and expand attached to electrodes. The impact of homogeneous isotropic turbulence is also explored, showing the persistence of ignition kernels early on but eventually quenching due to enhanced radical and heat losses with pronounced turbulence intensity.
Describing flow resistance using the physical properties of an underlying surface is a recalcitrant problem in overland flow models. If discharge measurements are available, an equivalent roughness (e.g., Manning’s n) can be calibrated to represent the effects of surface properties within the domain with a single numerical value. Alternatively, the flow resistance can be estimated from discharge and velocity measured at a point, typically a runoff plot outlet. However, such experimental estimates are often inconsistent with the equivalent roughness determined from calibration to discharge, even if both derive from the same dataset. For example, if Manning’s equation is used to parameterize flow resistance, the Manning’s n obtained by calibrating a model to discharge differs from the value of n calculated from measured flow and velocity at the hillslope outlet. Here, this discrepancy is resolved by deriving a correction factor relating experimentally-determined flow resistance to the equivalent roughness. The derived correction factor is tested for four commonly-used resistance formulations using 129 rainfall simulator experiments. The correction factor is necessary to reproduce measured velocities, and yields minor improvements in discharge prediction. Plain Language Summary: Accurate runoff prediction is needed for land and water management in dryland regions, where sporadic and limited rainfall necessitate efficient water use and drought mitigation strategies. The skill of runoff models is known to be hindered by out ability to estimate flow resistance, which is the quantity that describes how energy is lost from flowing water to the underlying surface. Typically, models represent flow resistance with an equivalent roughness, e.g., Manning’s n, that is adjusted until the model can reproduce available discharge observations at watershed scale. However, the flow resistance measured in plot-scale experiments (1–10 m) often exceeds equivalent roughness coefficients by a factor of 10. This means that the direct use of plot-scale experimental data to parameterize runoff models could cause errors in discharge and runoff velocity predictions. Here, we resolve these differences by deriving an analytic correction factor that relates flow resistance to the equivalent roughness required for models to reproduce experimental velocity and discharge data. This correction factor is tested using rainfall simulator data from 129 experiments performed in the US Southwest covering a wide range of precipitation intensities, soil textures and vegetation types. Use of the correction factor substantially improves model prediction of flow velocity, which is needed for reproducing the timing of flood events and the estimation of erosion.
Previous 2D radio frequency interferometric observations showed that lightning is commonly started with a positive fast discharge (+FD) before a normal negative leader continues from the origin of the +FD to the ensuing lightning flash. However, the inception and development of the +FD cannot be convincingly explained by existing discharge theories. With our new 3D broadband interferometric mapping and polarization system, we observed that the +FD was sometimes followed by an even faster and more extensive negative discharge (−FD) that propagated backward and overshoot the origin by a few hundred meters, as reported by an earlier 2D study. Surprisingly, the signal polarization, which measures the orientation of the discharge current, systematically slanted from the discharge propagation direction and rotated between the two opposite discharges, showing the +FD and −FD were driven by other storm-independent factors in addition to the storm electric field, or otherwise all would align in the same direction. Assuming a cosmic-ray shower (CRS) piercing through the cloud immediately before the +FD/−FD discharge, we found that their path is consistent with a pre-ionized path by the CRS, and the polarizations for the two opposite discharges are consistent with the respective deflected trajectories of high-energy positrons and electrons in the geomagnetic and an electric field. We further analyzed the more commonly observed +FD and showed that it is consistent with the CRS interpretation, suggesting these flashes in thunderstorms were ignited by cosmic-ray showers.
Here, this work numerically studies the plasma assisted deflagration to detonation transition (DDT) of H 2 /O 2 mixtures in a microscale channel with detailed chemistry and transport. The results show that the DDT onset time is non-monotonically dependent on the discharge pulse number. The DDT is accelerated with small pulse numbers, whereas retarded with large ones. Two different DDT regimes, respectively at a small and large plasma discharge number, via acoustic choking of the burned gas and plasma-enhanced reactivity gradient without acoustic choking, are observed. Without plasma discharge, pronounced pressure and temperature gradients in front of the flame are generated by acoustic compression after the choking of the burned gas, triggering DDT via autoignition. With small plasma pulse numbers, the plasma-generated species enhance the ignition kinetics and lead to an increased reactivity in the boundary layer. After the choking of the burned gas, the plasma-enhanced reactivity advances the sequence of autoignition near the wall, strengthens ignition-shock wave coupling, and accelerates DDT. However, with a large discharge pulse number, a direct autoignition initiating DDT can occur without the acoustic choking of the burned gas due to the strongly accelerated reactivity and elevated temperature. In this case, DDT onset is retarded because the elevated temperature increases sonic velocity and the increased reactivity accelerates fuel oxidation in front of the flame, decelerating the formation of a leading shock and subsequent pressure buildup ahead of the flame. The present modeling reveals that no matter with or without plasma discharge, DDT is initiated by autoignition in thermal, pressure, and reactivity gradient fields via the Zel'dovich gradient mechanism. The acoustic choking of the burned gas may not be the necessary condition of DDT with strong plasma-enhanced reactivity gradient. This work provides an answer to the experimentally observed non-monotonic DDT onset time by plasma, which provides guidance to control DDT in advanced detonation engines and fire safety of hydrogen-fueled catalytic reactors in microchannels by non-equilibrium plasma discharge.
Groundwater dependent ecosystems in areas with industrial and military land use are at risk of direct exposure to a wide range of contaminants, including PFAS chemicals. Glaciated terrain often has mixed high and low permeability sediments coupled with groundwater flow-through lake features. These hydrogeologic attributes create highly complex ‘source to seep’ dynamics that make spatiotemporal contaminant transport patterns difficult to predict. We investigated one such system in detail using a suite of heat-tracing and chemical methods. Numerous (n=57) preferential groundwater discharge zones (vertical flux rates ranging 0.2 to 3.2 m/d) were identified across the upper Quashnet River stream-wetland system in Mashpee, MA, USA, adjacent to an Air Force Base with several known PFAS source areas. Surface-water and groundwater samples were collected and analyzed (n=145) for precursors and terminal PFAS compounds between March and September 2022. Samples were collected at identified seeps along the Quashnet River (n=59), from wells upgradient from the stream-wetland system (n= 44), from contributing flow-through kettle lakes (n=8), and at multiple locations along the Quashnet River (n=34). Samples from seeps and wells had measured PFAS concentrations ranging from non-detect to approximate 3,500 ng/L (mean= 1,650 ng/L), and a range of deuterium excess values (3.2 to 15.9 per mil) indicative of varying degrees of groundwater-lake interaction prior to emergence at the discharge zones. Groundwater-lake interaction along flowpaths that sourced the sampled seeps was farther supported by significant correlations (p < 0.01) between deuterium excess and %PFAS precursors, and between %PFAS precursors and multiple terminal PFAS compounds (e.g., PFPeS, PFBS, PFHxS). However, some sampled seeps contributing groundwater to the stream-wetland system had much higher total PFAS concentrations (>1000 ng/L) than the upgradient kettle lakes, despite showing lake (evaporative) isotopic signatures, indicating the potential for groundwater flowpath convergence at wetland discharge zones and the influence of lakebed PFAS precursor reactions. PFAS compounds and water isotopic composition at sampled multilevel groundwater wells, rivers, lakes, and seeps suggest that a complex mixture of source groundwater and flowpath characteristics are responsible for diverse observed PFAS mixtures at preferential discharge zones across the stream-wetland system. Further, total PFAS loading patterns to the Quashnet River via groundwater discharge remained remarkably similar from winter to summer to fall conditions, despite a regional dry period in late summer 2022 with the upper river channel completely drying. This work addresses gaps in the existing PFAS literature by demonstrating the importance of subsurface fate and transport on PFAS compound concentrations in controlling contaminant mass loading in preferential groundwater discharge zones and presents a transferrable field toolkit for efficient characterization of spatially preferential PFAS transport dynamics.
Density functional theory (DFT) is utilized to explore the effects of increasing concentrations of vanadium (V) and chromium (Cr) substitution on the discharge of α-MnO 2 cathode in hydrated zinc-ion batteries (ZIBs). During H + -intercalation and Zn 2+ -intercalation, Cr-substitution proves to be more effective than V-substitution in promoting discharge behaviors. Transitioning from Mn 0.875 Cr 0.125 O 2 , Mn 0.75 Cr 0.25 O 2 to Mn 0.625 Cr 0.375 O 2 is found to consistently enhance the discharge voltage, along with improved tunnel structure retention and volume expansion suppression. In comparison, the promoting effect of increasing V-substitution is relatively small at initial discharge stages and leads to degradation at later stages, primarily due to an increased concentration of unstable Mn 2+ ion. The superior effect of Cr-substitution is attributed to the unique atomic and electronic structures of substituted Cr 4+ and reduced Cr 3+ ions during discharge. These ions serve as active electron acceptors to limit the formation of Mn 3+ and Mn 2+ ions, and as anchors to stabilize the α-MnO 2 framework and intercalated H + /Zn 2+ ions, respectively. Finally, our study highlights fine-tuning through substitution to enhance the performance of α-MnO 2 -based cathode materials in ZIBs.
This study investigates historical and future coastal runoff across the Arctic and sub-Arctic subregions of Alaska, USA, emphasising changes in terrestrial hydrology and their implications for coastal discharge. Using the physically based SnowModel-HydroFlow framework, historical (1991–2020) and future (2071–2100) coastal discharge climatologies were simulated, incorporating meteorological, topographic, and land cover inputs. Results reveal distinct seasonal and spatial patterns in temperature, precipitation, and snow dynamics across the Beaufort, Chukchi, Yukon, and Bristol subregions. Future projections suggest significant warming and shifts in hydrological regimes, such as increased annual precipitation and coastal discharge totals. Following a high emissions, business as usual pathway (SSP5-8.5), coastal discharge hydrographs in the Beaufort and Yukon subregions remain snow dominated with increases in discharge in nearly all months, while Chukchi and Bristol subregions shift from snow dominated to rain dominated. These changes are linked to rising temperatures, decreasing snow precipitation fractions, and altered freezing levels, with implications for water storage, timing, and availability. This research provides insights into the drivers and consequences of hydrological changes in a rapidly warming Arctic, supporting efforts to predict and mitigate the impacts on coastal and marine ecosystems.
Self-discharge and chemically induced mechanical effects degrade calendar and cycle life in intercalation-based electrochromic and electrochemical energy storage devices. In rechargeable lithium-ion batteries, self-discharge in cathodes causes voltage and capacity loss over time. The prevailing self-discharge model centers on the diffusion of lithium ions from the electrolyte into the cathode. Here, we demonstrate an alternative pathway, where hydrogenation of layered transition metal oxide cathodes induces self-discharge through hydrogen transfer from carbonate solvents to delithiated oxides. In self-discharged cathodes, we further observe opposing proton and lithium ion concentration gradients, which contribute to chemical and structural heterogeneities within delithiated cathodes, accelerating degradation. Hydrogenation occurring in delithiated cathodes may affect the chemo-mechanical coupling of layered cathodes as well as the calendar life of lithium-ion batteries.
The lithium vapor divertor concept is being developed as a method to achieve detached divertor conditions in a tokamak while minimizing impurity radiation losses from the core plasma. SOLPS-ITER modeling has previously been used to identify some of the geometric constraints and required lithium evaporation rate of a lithium vapor divertor in a medium-sized tokamak during steady-state operation. Here an updated conceptual design based on these operating requirements is introduced and the thermal response of the system is modeled during cyclical operation, consistent with operation in a short-pulse tokamak. Controllability of the temperature of the lithium capillary porous system (CPS) is achieved by adopting a design where there is no line-of-sight for radiation from the plasma to reach the heated CPS surface. Operational strategies to minimize the amount of lithium evaporated between plasma discharges while achieving steady evaporation rates during plasma discharges are discussed and modeled here. The optimal feedforward control strategy demonstrated in this work is to ramp up the temperature of the evaporator as quickly as possible immediately before a plasma discharge and then reduce the heating to match the desired steady-state net evaporation rate just before the plasma discharge begins, allowing the thermal inertia of the system to stabilize the evaporation rate during the first second of the plasma discharge.