Search NASA⌕ Search

SEARCH · Search NASA

Results for “SEPARATION PROCESSES”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A multi-step nucleation process determines the kinetics of prion-like domain phase separation

Compartmentalization by liquid-liquid phase separation (LLPS) has emerged as a ubiquitous mechanism underlying the organization of biomolecules in space and time. Here, we combine rapid-mixing time-resolved small-angle X-ray scattering (SAXS) approaches to characterize the assembly kinetics of a prototypical prion-like domain with equilibrium techniques that characterize its phase boundaries and the size distribution of clusters prior to phase separation. We find two kinetic regimes on the micro- to millisecond timescale that are distinguished by the size distribution of clusters. At the nanoscale, small complexes are formed with low affinity. After initial unfavorable complex assembly, additional monomers are added with higher affinity. At the mesoscale, assembly resembles classical homogeneous nucleation. Careful multi-pronged characterization is required for the understanding of condensate assembly mechanisms and will promote understanding of how the kinetics of biological phase separation is encoded in biomolecules.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Theoretical prediction and atomic kinetic Monte Carlo simulations of void superlattice self-organization under irradiation

Nano-structured superlattices may have novel physical properties and irradiation is a powerful mean to drive their self-organization. However, the formation mechanism of superlattice under irradiation is still open for debate. Here we use atomic kinetic Monte Carlo simulations in conjunction with a theoretical analysis to understand and predict the self-organization of nano-void superlattices under irradiation, which have been observed in various types of materials for more than 40 years but yet to be well understood. The superlattice is found to be a result of spontaneous precipitation of voids from the matrix, a process similar to phase separation in regular solid solution, with the symmetry dictated by anisotropic materials properties such as one-dimensional interstitial atom diffusion. This discovery challenges the widely accepted empirical rule of the coherency between the superlattice and host matrix crystal lattice. The atomic scale perspective has enabled a new theoretical analysis to successfully predict the superlattice parameters, which are in good agreement with independent experiments. The theory developed in this work can provide guidelines for designing target experiments to tailor desired microstructure under irradiation. Furthermore, it may also be generalized for situations beyond irradiation, such as spontaneous phase separation with reaction.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Carbon nanotube nanofluidics

Fluid flow under extreme spatial confinement exhibits unusual physical behaviors. This nanofluidic transport regime is relevant to a variety of mass transport, separation, and energy production processes in biological and industrial systems. Carbon nanotubes (CNTs) offer a nearly ideal platform for exploring nanofluidic transport because of their extremely narrow, smooth, hydrophobic inner pores, which enable very fast molecular flow while providing strong selectivity. In this review, we aim to provide a comprehensive understanding of nanofluidics in CNTs, focusing on the basic physics of mass transport in CNTs, various experimental platforms developed to investigate these phenomena, and key results on the permeation of water, protons, and ions. We focus on the critical factors that influence transport efficiency and selectivity, such as slip flow and charge regulation in CNTs, and the roles of entrance effects, dehydration processes and ion–charge interactions at the CNT entrances. We also explore the confinement effects, highlighting how the unique one-dimensional structure of CNTs imposes distinct constraints on fluid behavior and leads to novel single-file transport phenomena. Finally, we address current challenges and future directions of CNT nanofluidics.

Li, Zhongwu [Lawrence Livermore National Laborator↗
DOE OSTI2015

GaAs quantum dot solar cell under concentrated radiation

Effects of concentrated solar radiation on photovoltaic performance are investigated in well-developed GaAs quantum dot (QD) solar cells with 1-Sun efficiencies of 18%–19%. In these devices, the conversion processes are enhanced by nanoscale potential barriers and/or AlGaAs atomically thin barriers around QDs, which prevent photoelectron capture to QDs. Under concentrated radiation, the short circuit current increases proportionally to the concentration and the open circuit voltage shows the logarithmic increase. In the range up to hundred Suns, the contributions of QDs to the photocurrent are proportional to the light concentration. The ideality factors of 1.1–1.3 found from the V{sub OC}-Sun characteristics demonstrate effective suppression of recombination processes in barrier-separated QDs. The conversion efficiency shows the wide maximum in the range of 40–90 Suns and reaches 21.6%. Detailed analysis of I-V-Sun characteristics shows that at low intensities, the series resistance decreases inversely proportional to the concentration and, at ∼40 Suns, reaches the plateau determined mainly by the front contact resistance. Improvement of contact resistance would increase efficiency to above 24% at thousand Suns.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Growth far from equilibrium: Examples from III-V semiconductors

The development of new applications has driven the field of materials design and synthesis to investigate materials that are not thermodynamically stable phases. Materials which are not thermodynamically stable can be synthesized and used in many applications. These materials are kinetically stabilized during use. The formation of such metastable materials requires both an understanding of the associated thermochemistry and the key surface transport processes present during growth. Phase separation is most easily accomplished at the growth surface during synthesis where mass transport is most rapid. These surface transport processes are sensitive to the surface stoichiometry, reconstruction, and chemistry as well as the growth temperature. The formation of new metastable semiconducting alloys with compositions deep within a compositional miscibility gap serves as model systems for the understanding of the surface chemical and physical processes controlling their formation. The GaAs{sub 1−y}Bi{sub y} system is used here to elucidate the role of surface chemistry in the formation of a homogeneous metastable composition during the chemical vapor deposition of the alloy system.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Aqueous density fractionation of nanocarbons from colliding-wave PBX9502 detonation soot

The study of solid carbon (nanocarbon) nucleation and growth from detonating high explosives is limited because of the difficulty directly evaluating chemical reactions behind the shock front. Understanding the correlation between detonation conditions and nanocarbon formation requires isolation, purification, and analysis of the carbon products. Aqueous density fractionation of products recovered from a colliding wave TATB-based high explosive (PBX9502) detonation are evaluated as an environmentally benign, “green” chemical processing approach for the separation of carbon materials from metal and metal oxide byproducts. Recovery of material sufficient for X-ray scattering analysis was obtained in two fractions. X-ray scattering showed the less dense fraction (ρ < 1 g/cm 3 ) contains hollow spherical core (d avg = 27.8 nm) nanoparticles with a shell (thickness=3.6 nm) composed of predominately sp 2 hybridized carbon. Particles contained within the dense fraction (ρ > 1 g/cm 3 ) were a mixture of oblate spheroids (major axis length=43.8 nm and minor axis=2.9 nm) and small polydisperse solid spherical particles (d avg = 2.4 nm) composed of both sp 2 and sp 3 bonded carbon with residual contamination from the brass detonators.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

$\rm {H\,\small {I}}$ constraints from the cross-correlation of eBOSS galaxies and Green Bank Telescope intensity maps

We present the joint analysis of Neutral Hydrogen ($\rm {H\,\small {I}}$) Intensity Mapping observations with three galaxy samples: the Luminous Red Galaxy (LRG) and Emission Line Galaxy (ELG) samples from the eBOSS survey, and the WiggleZ Dark Energy Survey sample. The $\rm {H\,\small {I}}$ intensity maps are Green Bank Telescope observations of the redshifted $21\rm cm$ emission on $100 \, {\rm deg}^2$ covering the redshift range 0.6 < z < 1.0. We process the data by separating and removing the foregrounds present in the radio frequencies with FastI ICA. We verify the quality of the foreground separation with mock realizations, and construct a transfer function to correct for the effects of foreground removal on the $\rm {H\,\small {I}}$ signal. We cross-correlate the cleaned $\rm {H\,\small {I}}$ data with the galaxy samples and study the overall amplitude as well as the scale dependence of the power spectrum. We also qualitatively compare our findings with the predictions by a semianalytical galaxy evolution simulation. The cross-correlations constrain the quantity $\Omega _{\rm {H\,\small {I}}} b_{\rm {H\,\small {I}}} r_{\rm {H\,\small {I}},{\rm opt}}$ at an effective scale k eff , where $\Omega _\rm {H\,\small {I}}$ is the $\rm {H\,\small {I}}$ density fraction, $b_\rm {H\,\small {I}}$ is the $\rm {H\,\small {I}}$ bias, and $r_{\rm {H\,\small {I}},{\rm opt}}$ the galaxy–hydrogen correlation coefficient, which is dependent on the $\rm {H\,\small {I}}$ content of the optical galaxy sample. At $k_{\rm eff}=0.31 \, h\,{\rm Mpc^{-1}}$ we find $\Omega _{\rm {H\,\small {I}}} b_{\rm {H\,\small {I}}} r_{\rm {H\,\small {I}},{\rm Wig}} = [0.58 \pm 0.09 \, {\rm (stat) \pm 0.05 \, {\rm (sys)}}] \times 10^{-3}$ for GBT-WiggleZ, $\Omega _{\rm {H\,\small {I}}} b_{\rm {H\,\small {I}}} r_{\rm {H\,\small {I}},{\rm ELG}} = [0.40 \pm 0.09 \, {\rm (stat) \pm 0.04 \, {\rm (sys)}}] \times 10^{-3}$ for GBT-ELG, and $\Omega _{\rm {H\,\small {I}}} b_{\rm {H\,\small {I}}} r_{\rm {H\,\small {I}},{\rm LRG}} = [0.35 \pm 0.08 \, {\rm (stat) \pm 0.03 \, {\rm (sys)}}] \times 10^{-3}$ for GBT-LRG, at z ≃ 0.8. We also report results at $k_{\rm eff}=0.24$ and $k_{\rm eff}=0.48 \, h\,{\rm Mpc^{-1}}$. With little information on $\rm {H\,\small {I}}$ parameters beyond our local Universe, these are amongst the most precise constraints on neutral hydrogen density fluctuations in an underexplored redshift range.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Ab initio phonon dispersion in crystalline naphthalene using van der Waals density functionals

Acene molecular crystals are of current interest in organic optoelectronics, both as active materials and for exploring and understanding new phenomena. Phonon scattering can be an important facilitator and dissipation mechanism in charge separation and carrier transport processes. Here, we carry out density functional theory (DFT) calculations of the structure and the full phonon dispersion of crystalline naphthalene, a well-characterized acene crystal for which detailed neutron-diffraction measurements, as well as infrared and Raman spectroscopy, are available. We evaluate the performance, relative to experiments, of DFT within the local density approximation (LDA); the generalized gradient approximation of Perdew, Burke, and Ernzerhof (PBE); and a recent van der Waals-corrected nonlocal correlation (vdW-DF-cx) functional. We find that the vdW-DF-cx functional accurately predicts lattice parameters of naphthalene within 1%. Intermolecular and intramolecular phonon frequencies across the Brillouin zone are reproduced within 7.8% and 1%, respectively. As expected, LDA (PBE) underestimates (overestimates) the lattice parameters and overestimates (underestimates) phonon frequencies, demonstrating their shortcomings for predictive calculations of weakly bound materials. If the unit cell is fixed to the experimental lattice parameters, PBE is shown to lead to improved phonon frequencies. Our study provides a detailed understanding of the phonon spectrum of naphthalene, and highlights the importance of including van der Waals dispersion interactions in predictive calculations of lattice parameters and phonon frequencies of molecular crystals and related organic materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Measurement of the 230 Th( p ,2n)Pa229 and 230 Th( p ,3n)Pa228 reaction cross sections from 14.1 to 16.9 MeV

Actinium-225 is of interest for medical isotope production and there is on-going research into methods of producing Ac 225 , either directly or via the decay of its parent isotopes ( Th 229 , Pa 229 , and Ra 225 ). One method that has been suggested is the Th 230 ( p , 2 n ) Pa 229 reaction. However, there is no available cross-section data for this reaction in the literature. Purpose: Measure the Th 230 ( p , 2 n ) and Th 230 ( p , 3 n ) reaction cross sections in the energy range where the ( p , 2 n ) reaction is predicted to peak to determine the feasibility of Ac 225 production via the Th 230 ( p , 2 n ) reaction. Methods: Targets naturally enriched in Th 230 were irradiated at the Center for Accelerator Mass Spectrometry at Lawrence Livermore National Laboratory with energies ranging from 14.1 to 16.9 MeV. Furthermore, chemical processing was used to separate the protactinium activation products, followed by γ -ray spectroscopy to measure the activities of Pa 228 , 229 , 230 , 232 produced in the irradiation. Results: We find that excitation functions are reported for the first time in the literature for the Th 230 ( p , 2 n ) and Th 230 ( p , 3 n ) reactions in this energy range. The peak measured value of the Th 230 ( p , 2 n ) reaction was found to be 182 ± 12 mb at 14.4 ± 0.1 MeV. The Th 232 ( p , n ) Pa 232 reaction was used to verify the experimental conditions, the measured values are reported and are comparable to the existing literature values. From the γ -ray spectrometry data, the half-life of Pa 229 was measured as 1.5 ± 0.1 days, which is within the error of the half-life reported in the evaluated nuclear data as well as in the recent measurements, and the half-life of Pa 228 was measured as 19.5 ± 0.4 hours. Conclusions: Overall, the Th 230 ( p , 2 n ) Pa 229 reaction could reasonably be used for Ac 225 isotope production, although significant amounts of relatively isotopically pure Th 230 would be needed for significant production because the low alpha-decay branching ratio of Pa 229 and long half-life of Th 229 inhibit the in-growth of significant amounts of Ac 225 .

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2018

Evidence for the Associated Production of a Single Top Quark and a Photon in Proton-Proton Collisions at s = 13 TeV

The first evidence of events consistent with the production of a single top quark in association with a photon is reported. The analysis is based on proton-proton collisions at √s=13 TeV and recorded by the CMS experiment in 2016, corresponding to an integrated luminosity of 35.9 fb -1 . Events are selected by requiring the presence of a muon (μ), a photon (γ), an imbalance in transverse momentum from an undetected neutrino (ν), and at least two jets (j) of which exactly one is identified as associated with the hadronization of a b quark. A multivariate discriminant based on topological and kinematic event properties is employed to separate signal from background processes. An excess above the background-only hypothesis is observed, with a significance of 4.4 standard deviations. A fiducial cross section is measured for isolated photons with transverse momentum greater than 25 GeV in the central region of the detector. The measured product of the cross section and branching fraction is σ(pp→tγj)B(t→μνb)=115±17(stat)±30(syst) fb, which is consistent with the standard model prediction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

SU-E-T-423: Fast Photon Convolution Calculation with a 3D-Ideal Kernel On the GPU

Purpose: The calculation time is a trade-off for improving the accuracy of convolution dose calculation with fine calculation spacing of the KERMA kernel. We investigated to accelerate the convolution calculation using an ideal kernel on the Graphic Processing Units (GPU). Methods: The calculation was performed on the AMD graphics hardware of Dual FirePro D700 and our algorithm was implemented using the Aparapi that convert Java bytecode to OpenCL. The process of dose calculation was separated with the TERMA and KERMA steps. The dose deposited at the coordinate (x, y, z) was determined in the process. In the dose calculation running on the central processing unit (CPU) of Intel Xeon E5, the calculation loops were performed for all calculation points. On the GPU computation, all of the calculation processes for the points were sent to the GPU and the multi-thread computation was done. In this study, the dose calculation was performed in a water equivalent homogeneous phantom with 150{sup 3} voxels (2 mm calculation grid) and the calculation speed on the GPU to that on the CPU and the accuracy of PDD were compared. Results: The calculation time for the GPU and the CPU were 3.3 sec and 4.4 hour, respectively. The calculation speed for the GPU was 4800 times faster than that for the CPU. The PDD curve for the GPU was perfectly matched to that for the CPU. Conclusion: The convolution calculation with the ideal kernel on the GPU was clinically acceptable for time and may be more accurate in an inhomogeneous region. Intensity modulated arc therapy needs dose calculations for different gantry angles at many control points. Thus, it would be more practical that the kernel uses a coarse spacing technique if the calculation is faster while keeping the similar accuracy to a current treatment planning system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Recycling of multilayer plastic packaging materials by solvent-targeted recovery and precipitation

Many plastic packaging materials manufactured today are composites made of distinct polymer layers (i.e., multilayer films). Billions of pounds of these multilayer films are produced annually, but manufacturing inefficiencies result in large, corresponding postindustrial waste streams. Although relatively clean (as opposed to municipal wastes) and of near-constant composition, no commercially practiced technologies exist to fully deconstruct postindustrial multilayer film wastes into pure, recyclable polymers. Here, we demonstrate a unique strategy we call solvent-targeted recovery and precipitation (STRAP) to deconstruct multilayer films into their constituent resins using a series of solvent washes that are guided by thermodynamic calculations of polymer solubility. We show that the STRAP process is able to separate three representative polymers (polyethylene, ethylene vinyl alcohol, and polyethylene terephthalate) from a commercially available multilayer film with nearly 100% material efficiency, affording recyclable resins that are cost-competitive with the corresponding virgin materials.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Improved Electrochemical Performance of Carbon-Coated LiFeBO 3 Nanoparticles for Lithium-Ion Batteries

Carbon-coated LiFeBO 3 nanoparticles have been successfully prepared by surfactant-assisted ball milling and a size selection process based on centrifugal separation. We observed monodispersed LiFeBO 3 nanoparticles with dimensions of 10–20 nm by transmission electron microscope. The introduced surfactant acts as the dispersant as well as the carbon source for LiFeBO 3 nanoparticles. Greatly improved discharge capacities of 190.4 mA h g –1 at 0.1 C and 106.6 mA h g –1 at 1 C rate have been achieved in the LiFeBO 3 nanoparticles when cycling the cells between 1.0 V and 4.8 V. Meanwhile, the as-prepared micro-size LiFeBO 3 electrodes show lower discharge capacities of 142 mA h g –1 and 93.3 mA h g –1 at 0.1 C and 1 C rates. Moreover, the post-treated LiFeBO 3 nanostructure has drastically enhanced the electrochemical performance due to the short diffusion length and ameliorated electrical contract between LiFeBO 3 nano particles.

14 SOLAR ENERGY↗
DOE OSTI2020

Production of Cf-252 and other transplutonium isotopes at Oak Ridge National Laboratory

In 1957 Glenn T. Seaborg conceived and advocated for the construction of the High Flux Isotope Reactor (HFIR) and the Transuranium Processing Plant (since then renamed the Radiochemical Engineering Development Center, or REDC) at Oak Ridge National Laboratory. There, heavily shielded hot cells, glove boxes, and laboratories allow recovery of transuranium elements produced in substantial quantities. Seaborg’s vision of HFIR and REDC producing milligram quantities of berkelium, californium, and einsteinium has been fulfilled beginning in 1966 through May 2019 with 78 production campaigns yielding a cumulative totals of 1.2 g of 249 Bk, 10.2 g of 252 Cf, 39 mg of 253 Es, and 15 pg of 257 Fm. Notably, 252 Cf is a neutron source used in many industrial applications including oil exploration; process control systems for the cement industry, coal analysis, and power production; sources to start nuclear reactors and perform nondestructive materials analyses; homeland security and national defense detection devices; and medical research. Isotopes made available through transplutonium production at HFIR/REDC have enabled scientists to study the nuclear properties and reactions, chemical properties, optical properties, and solid-state properties of transplutonium elements. Long-lived isotopes have served as targets in heavy ion accelerators to produce heavier elements leading to the discovery of 104 Rf, 105 Db, 106 Sg, 113 Nh, 114 Fl, 115 Mc, 116 Lv, 117 Ts, and 118 Og. This paper reviews the evolution of the processing flowsheets to produce, separate, and purify transplutonium isotopes, which have evolved over 50 years of operation at HFIR and REDC, and summarizes directions of future work to improve the efficiency of the production operations.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Raw tensiometer soil matric potential data, from Oct 2017-Jun 2020 in Manaus, Brazil

This data package contains raw output from two Campbell Scientific CR1000 data loggers each connected to three UMS T4e tensiometer sensors in Manaus, Brazil. Tensiometers at each data logger location were installed to depths of 10 cm, 20 cm and 50 cm below the ground surface for the purpose of measuring soil matric potential at each depth. “Station 3” tensiometers were installed at approximately 2°36'31.68"S, 60°12'34.70"W. “Station 4” tensiometers were installed at approximately 2°36'32.46"S, 60°12'34.92"W. The files “tensio3.edit.202109.csv” and “tensio4.edit.202109.csv” contain raw data (with invalid data removed). A text editor will be required to view the raw files in .csv (comma-separated values) format. Data processing was limited to: 1) Editing column headers for clarity; 2) Applying a correction that converts voltage read by the data loggers to pressure units (kPa); 3) removing all data greater than 50 kPa and less than -115 kPa; 4) Adding columns entitled “Caution.10cm”, “Caution.20cm” and “Caution.50cm” that contains values of either “NA” or “Caution”, the latter of which is applied when data exceeds the manufacturer reported data validity range of less than -85 kPa to suggest use of these data with caution; 5) Adding columns entitled “Invalid.10cm”, “Invalid.20cm”, and “Invalid.50cm” that contains values of either “NA” or “Invalid”, the latter of which is applied when data accuracy is deemed compromised due to the introduction of air bubbles into the instrument. Data marked “Invalid” may still indicate general trends of wet/high matric potential and dry/low matric potential, but values are inaccurate and dry/low matric potential values are likely not low enough. Wet/high matric potential values are likely to be more accurate than dry/low matric potential values; 6) Data format changed from .dat to .csv. This research was performed as a part of the NGEE Tropics project, which aims to advance model predictions of tropical forest carbon cycle responses to a changing climate over the 21st Century.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Draft Rules for Alloy 617 Creep-Fatigue Design Using an EPP+SMT Approach

The report describes an improved new creep-fatigue design method for structural components in high temperature nuclear service. The new method uses an integrated elastic-perfectly plastic (EPP) analysis and Simplified Model Test (SMT) approach that greatly simplifies the design process by avoiding the separate evaluation of creep and fatigue damages and eliminating the requirement of stress classification. The report includes draft design rules, presented in a format compatible with an ASME nuclear Code Case, as well as a commentary on the rules and the validation of the rules using pressurized SMT (p-SMT) test data. Areas of potential improvement have been identified for further evaluation. The report also develops EPP+SMT design charts for Alloy 617 at temperatures between 800°C and 950°C by extrapolating high strain range, short hold test data to low strain range and long hold times, typical for structural components in high temperature nuclear service. The design charts are currently limited to 10,000 cycles to avoid excessive extrapolation outside the test database. More robust alternate extrapolation procedures will be considered in future work. The report includes several verification problems comparing the new EPP+SMT design method with current ASME creep-fatigue design methods. The comparison demonstrates that the EPP+SMT method significantly reduces the over-conservatism in current methods. Finally, the report includes a fully-documented sample problem with multiple load cases detailing the application of proposed EPP+SMT design rules.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Transient cycling of nitrogen, organic carbon and oxygen within the free-flowing Columbia River corridor: Linking exposure time dependent biogeochemical reactions to river stage fluctuations (Final project report)

The objective of this project was to develop next-generation techniques for representing the transformations of complex reaction networks like those impacting transient river-corridors, then apply them to develop upscaling tools. The project contributed significantly to the development of Lagrangian “mass transfer particle tracking” (MTPT) tools that are the first numerical framework that explicitly separates mixing and spreading processes; a crucial distinction for accurate representations of reaction rates across scales. The MTPT approach uses an operator splitting scheme where physical transport processes are simulated using classical random walk methods and each particle is treated as a “container” that carries with it any number of chemical species. Mixing between particles is simulated using a colocation-probability based mass transfer kernel and reactions are evaluated on each particle after mixing. The method is stable, accurate, and also allows the explicit modeling of residence time distributions of the individual containers of mass. The MTPT scheme allowed us to explore several research questions related to reactive transport in river corridors and complex reaction networks. We found that it is possible to quantify the exposure time distributions (ETDs) of reactants and that this does suggest some pathways to upscaling. However, comparable approximations can be obtained using simplified (reduced dimensionality) MTPT simulations in less time with similar uncertainty, so development of extensive ETD-based methods was not productive. The major findings are that i) the MTPT schemes are robust and highly accurate across a range of arbitrarily complex reactions, ii) the methods can be efficiently parallelized and the parallel performance characteristics of MTPT are predictable, and iii) the MTPT tools allow simultaneous tracking of residence time. The development of these research tools into efficient software packages continues. The methods are already available to the community because open-source, working examples have been included with all publications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Spark Channel Dynamics of Electrostatic Discharges

When two differently-charged objects are brought in close proximity to each other, the resulting high electric fields can cause electron avalanche breakdown of the air gap separating the objects, a process known as electrostatic discharge (ESD). If enough initial charge is stored on the objects, the electrical breakdown can proceed to ionize the air to such a degree that a highly conductive filament of plasma forms in the gap, known as a spark channel. The spark electrically bridges the air gap, resulting in a rapid pulse of current that neutralizes the charge difference. The current pulse produces significant heating of the gas in the spark, resulting in dissociation, ionization, thermal radiation, and hydrodynamic expansion. ESD presents a hazard to electrically-sensitive devices, with consequences such as economic losses (e.g. damaged electronics) or unsafe response (e.g. unintended ignition of flammable gas mixtures, initiation of detonators, etc.). For this thesis, the ESD spark is taken to occur between two conducting electrodes, with the spark channel being axisymmetric in a cylindrical coordinate system centered on the channel. An RLC-type circuit is used for the discharge model of the ESD event. The spark is treated as a time-dependent resistance that is in series with a capacitance, an inductance, and (optionally) a load resistance representing a “victim” component under threat from the ESD event. The primary motivation of this work is to use a numerical hydrodynamic model to understand the energy dissipation and transport processes in the spark. The model consists of the compressible Euler equations of mass, momentum, and energy conservation together with an Eddington/P1 approximation for thermal radiation transport. To close the hydrodynamic system, an equation of state (EOS) was fitted from tabular data for air that accounts for the dissociation and ionization of air species. The hydrodynamic equations are solved using a conservative Lagrangian finite volume method. These partial differential equations are coupled to the circuit equations by calculation of the spark resistance via numerical integration of the electrical conductivity of the channel. Computational results are compared against experimental measurements of discharge current and radial density of the spark channel.

42 ENGINEERING↗