Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mass diffusivity”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

DiffLense: a conditional diffusion model for super-resolution of gravitational lensing data

Abstract Gravitational lensing data is frequently collected at low resolution due to instrumental limitations and observing conditions. Machine learning-based super-resolution techniques offer a method to enhance the resolution of these images, enabling more precise measurements of lensing effects and a better understanding of the matter distribution in the lensing system. This enhancement can significantly improve our knowledge of the distribution of mass within the lensing galaxy and its environment, as well as the properties of the background source being lensed. Traditional super-resolution techniques typically learn a mapping function from lower-resolution to higher-resolution samples. However, these methods are often constrained by their dependence on optimizing a fixed distance function, which can result in the loss of intricate details crucial for astrophysical analysis. In this work, we introduce DiffLense , a novel super-resolution pipeline based on a conditional diffusion model specifically designed to enhance the resolution of gravitational lensing images obtained from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Our approach adopts a generative model, leveraging the detailed structural information present in Hubble space telescope (HST) counterparts. The diffusion model, trained to generate HST data, is conditioned on HSC data pre-processed with denoising techniques and thresholding to significantly reduce noise and background interference. This process leads to a more distinct and less overlapping conditional distribution during the model’s training phase. We demonstrate that DiffLense outperforms existing state-of-the-art single-image super-resolution techniques, particularly in retaining the fine details necessary for astrophysical analyses.

Computer Science↗
DOE OSTI2021

Composite fermion nonlinear sigma models

We study particle-hole symmetry at the integer quantum Hall plateau transition using composite fermion mean-field theory. Because this theory implicitly includes some electron-electron interactions, it also has applications to certain fractional quantum Hall plateau transitions. Previous work [P. Kumar et al., Phys. Rev. B 100, 235124 (2019)] using this approach showed that the diffusive quantum criticality of this transition is described by a nonlinear sigma model with topological θ = π term. This result, which holds for both the Dirac and Halperin, Lee, and Read composite fermion theories, signifies an emergent particle-hole (reflection) symmetry of the integer (fractional) quantum Hall transition. Here we consider the stability of this result to various particle-hole symmetry-violating perturbations. In the presence of quenched disorder that preserves particle-hole symmetry, we find that finite longitudinal conductivity at this transition requires the vanishing of a symmetry-violating composite fermion effective mass, which if present would generally lead to θ ≠ π and a corresponding violation of particle-hole symmetric electrical transport σ xy ≠ $\frac1{2}$$\frac{e^2}{h}$ . When the disorder does not preserve particle-hole symmetry, we find that θ can vary continuously within the diffusive regime. Our results call for further study of the universality of the quantum Hall plateau transition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Automatically adaptive stabilized finite elements and continuation analysis for compaction banding in geomaterials

Under compressive creep, viscoplastic solids experiencing internal mass transfer processes can accommodate singular cnoidal wave solutions as material instabilities at the stationary wave limit. These instabilities appear when the loading rate is significantly faster than the material's capacity to diffusive internal perturbations, leading to localized failure features (e.g., cracks and compaction bands). These cnoidal waves, generally found in fluids, have strong nonlinearities that produce periodic patterns. Due to the singular nature of the solutions, the applicability of the theory is currently limited. Additionally, practical simulation tools require proper regularization to overcome the challenges that singularity induces. We focus on the numerical treatment of the governing equation using a nonlinear approach building on a recent adaptive stabilized finite element method. This automatic refinement method provides an error estimate that drives mesh adaptivity, a crucial feature for the problem at hand. We compare the performance of this adaptive strategy against analytical and standard finite element solutions. We then investigate the sensitivity of the diffusivity ratio, the parameter controlling the process, and identify multiple possible solutions with several stress peaks. We also show the evolution of the spacing between peaks for all solutions as a function of that parameter.

42 ENGINEERING↗
DOE OSTI2024

Dynamics of dark matter misalignment through the Higgs portal

A light singlet scalar field feebly coupled through the super-renormalizable Higgs portal provides a minimal and well-motivated realization of ultra-light bosonic dark matter. We study the cosmological production of dark matter in this model by elucidating the dynamics of two sources of scalar field misalignment generated during the radiation era. For large scalar masses (above $\mathcal{O}$ (10 −3 eV)), dark matter is produced through thermal misalignment, by which the scalar field is driven towards large field values as a result of the finite-temperature effective potential. The dominance of thermal misalignment in this mass range leads to a sharp relic abundance prediction which is, to a significant extent, insensitive to the initial conditions of the scalar field. On the other hand, for low mass scalars (below $\mathcal{O}$ (10 −5 eV)), dark matter is produced via VEV misalignment, which is caused by the induced scalar field vacuum expectation value triggered by the electroweak phase transition. We show that the relic abundance in this low mass range is sensitive to the scalar field initial conditions. In the intermediate mass range, the relic abundance is a consequence of a competition between thermal misalignment and VEV misalignment, which can potentially lead to novel forced resonance effects which cause a recurring enhancement and suppression in the late time oscillation amplitude as a function of the scalar mass. We compare our relic abundance predictions with constraints and projections from equivalence principle and inverse square law tests, stellar cooling, resonant molecular absorption, and observations of extra-galactic background light and diffuse X-ray backgrounds. New experimental ideas are needed to probe most of the cosmologically motivated regions of parameter space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

In situ catalyst activation and regeneration enable energy-efficient high-current CO 2 reduction to ethanol-rich C 2+ mixtures

Electrochemical conversion of dissolved CO 2 in bicarbonate electrolytes, i.e., bicarbonate electrolysis, offers distinct advantages over gas diffusion electrode systems by enabling direct utilization of the CO 2 capture electrolyte while bypassing the energy-intensive CO 2 release step. However, bicarbonate electrolysis faces challenges such as CO 2 mass-transfer limitation, local pH-driven CO 2 depletion, and high cathodic potentials. The higher potential often causes catalyst surface reorganization, leading to a gradual loss of active sites and variations in selectivity during CO 2 reduction. Here, we report a directed, in situ activation and regeneration method that allows precatalysts to equilibrate under dynamic (pulsed) electrolysis conditions. We demonstrate in situ activation of a scalable Cu 2 O/Cu mesh that, under short-width (t = 4 s) pulsed electrolysis, provides stable mixed oxidation states of Cu, favoring the formation of an ethanol-rich crude mixture. The pulsed electrolysis waveform, consisting of six distinct segments, is tuned to form Cu + oxides, which are then reduced to generate local alkaline conditions favoring C–C coupling. This synergistic effect results in FEs of 73% for C2+ products and 39% for ethanol at an applied current density of −150 mA cm −2 and a cathodic potential of −1.45 V (vs. RHE). The overall half-cell energy efficiency is ∼30% for C 2+ products. The in situ Raman experiments confirm the role of pCO 2 R in dynamically regenerating Cu+-containing surface species during pulsed operation, thereby steering selectivity towards C 2+ products. A comprehensive multiscale, multiphysics model is developed to investigate the dynamic behavior of copper surface species (Cu, Cu + , and Cu 2+ ) and local microenvironmental conditions during the pCO 2 R. The results reveal that the coexistence of different copper oxidation states, especially the Cu+ intermediate, is critical in steering selectivity towards multicarbon (C 2+ ) products. The dynamic modulation of surface redox states via tailored pulsing strategies favors C–C coupling pathways by inducing localized alkaline conditions and stabilizing reactive intermediates. This work establishes a predictive modeling platform that links pulse waveform design with mechanistic insights into catalyst state evolution and product selectivity. Overall, this study provides valuable insights into the synergistic effect of in situ activation of pre-catalysts and pulsed electrolysis for higher selectivity towards C 2+ products.

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

Kinetic simulations of anomalous resistivity in high-temperature current carrying plasmas

Anomalous resistivity caused by lower hybrid drift instabilities (LHDIs) has been extensively studied in the literature, and has been invoked to explain the rates of magnetic diffusion and plasma profile evolution in the low-density plasma periphery of theta and screw pinch plasmas. In this article, we present a fully kinetic simulation study of LHDIs with a realistic mass ratio using the particle-in-cell code Chicago over a large range of drift speeds, 0.5 < v$^{2}_{di}$/v$^{2}_{ti}$ < 14.5, plasma β, β ≥ 2×10 –5 , and temperature ratios, 0.1 < T e /T i < 10. The resistivity quantified from the simulations is compared with the analytic estimates from the literature and a generalized resistivity expression is presented, which can be applied over the entire range of plasma environments that were simulated in this study. This expression uses the amplitude of the LHDI fluctuations resulting from instability saturation by electron resonance broadening rather than plateau formation, current relaxation, or ion trapping. The generalized resistivity expression is appropriate for full-scale fluid simulations that encounter a large range of plasma conditions, and has better agreement with the resistivity quantified from kinetic simulations than previous expressions in the literature, particularly at low plasma β. Anomalous resistivity has also been considered as an explanation of the rapid magnetic field dissipation during magnetic reconnection events at large plasma β. At a large β, the magnetic field fluctuations of LHDIs can dominate the electrostatic fluctuations and this article investigates enhanced collisionality due to the nonlinear magnetic force, δJ×δB/c.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Stellar splashback: the edge of the intracluster light

We examine the outskirts of galaxy clusters in the C-EAGLE simulations to quantify the ‘edges’ of the stellar and dark matter distribution. The radius of the steepest slope in the dark matter, commonly used as a proxy for the splashback radius, is located at $\sim \, r_{200 \rm m}$; the strength and location of this feature depends on the recent mass accretion rate, in good agreement with previous work. Interestingly, the stellar distribution (or intracluster light, ICL) also has a well-defined edge, which is directly related to the splashback radius of the halo. Thus, detecting the edge of the ICL can provide an independent measure of the physical boundary of the halo, and the recent mass accretion rate. We show that these caustics can also be seen in the projected density profiles, but care must be taken to account for the influence of substructures and other non-diffuse material, which can bias and/or weaken the signal of the steepest slope. This is particularly important for the stellar material, which has a higher fraction bound in subhaloes than the dark matter. Finally, we show that the ‘stellar splashback’ feature is located beyond current observational constraints on the ICL, but these large projected distances (>> 1 Mpc) and low surface brightnesses (μ>>32 mag arcsec -2 ) can be reached with upcoming observational facilities such as the Vera C. Rubin Observatory, the Nancy Grace Roman Space Telescope, and Euclid.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Numerical Simulations of Orographic Convection across Multiple Gray Zones

Idealized simulations are used to determine the sensitivity of moist orographic convection to horizontal grid spacing Δ h . In simulated mechanically (MECH) and thermally (THERM) forced convection over an isolated ridge, Δ h is varied systematically over both the deep-convection (Δ h ~ 10 –1 km) and turbulence (Δ h ~ 1 km–100 m) gray zones. To aid physical interpretation, a new parcel-based bulk entrainment/detrainment diagnosis for horizontally heterogeneous flows is developed. Within the deep-convection gray zone, the Δ h sensitivity is dominated by differences in parameterized versus explicit convection; the former initiates convection too far upstream of the ridge (MECH) and too early in the diurnal heating cycle (THERM). These errors stem in part from a large underprediction of parameterized entrainment and detrainment. Within the turbulence gray zone, sensitivities to Δh arise from the representation of both subcloud- and cloud-layer turbulence. As Δ h is decreased, MECH exhibits stronger cloud-layer entrainment to enhance the convective mass flux M co , while THERM exhibits stronger detrainment to suppress Mco and delay convection initiation. The latter is reinforced by increased subcloud turbulence at smaller Δ h , which leads to drying and diffusion of the central updraft responsible for initiating moist convection. Numerical convergence to a robust solution occurs only in THERM, which develops a fully turbulent flow with a resolved inertial subrange (for Δ h ≤ 250 m). In MECH, by contrast, turbulent transition occurs within the orographic cloud, the details of which depend on both physical location and Δ h .

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

Searching for High-energy Neutrino Emission from Galaxy Clusters with IceCube

Galaxy clusters have the potential to accelerate cosmic rays (CRs) to ultrahigh energies via accretion shocks or embedded CR acceleration sites. The CRs with energies below the Hillas condition will be confined within the cluster and eventually interact with the intracluster medium gas to produce secondary neutrinos and gamma rays. Using 9.5 yr of muon neutrino track events from the IceCube Neutrino Observatory, we report the results of a stacking analysis of 1094 galaxy clusters with masses ≳10 14 M ⊙ and redshifts between 0.01 and ~1 detected by the Planck mission via the Sunyaev–Zel’dovich effect. We find no evidence for significant neutrino emission and report upper limits on the cumulative unresolved neutrino flux from massive galaxy clusters after accounting for the completeness of the catalog up to a redshift of 2, assuming three different weighting scenarios for the stacking and three different power-law spectra. Weighting the sources according to mass and distance, we set upper limits at a 90% confidence level that constrain the flux of neutrinos from massive galaxy clusters (≳10 14 M ⊙ ) to be no more than 4.6% of the diffuse IceCube observations at 100 TeV, assuming an unbroken E -2.5 power-law spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Soil porous microstructure control over soil organic matter mobility: A multimethod workflow for understanding chemistry-dependent organic matter binding in soil

Soil organic matter (SOM) has attracted a great deal of interest; particularly for its potential to mitigate human derived CO 2 emissions. Studies have demonstrated that SOM plays a critical role in carbon storage and CO 2 sequestration. However, the sorption properties of SOM, which influence its transport in pore water and stabilization within the soil, remain poorly understood. This study develops a workflow to: (1) examine compound-specific advective and diffusive transport and desorption behaviors, (2) quantify desorption rates through stop-flow and continuous-flow column experiments, and (3) evaluate the impact of soil microporosity on SOM mobility using high-resolution imaging and extractions. Intact core column experiments were conducted on Uncultivated (Natural) and Cultivated soil samples, both were arid soils, collected in Washington State. X-ray computed tomography was employed to measure porosity and pore connectivity, while Fourier-transform ion cyclotron resonance mass spectrometry was used to analyze SOM composition. The findings revealed that cultivation increased total carbon and nitrogen levels due to irrigation and fertilization, enhancing carbon capture potential in arid soils. In contrast, the Natural soil, characterized by higher porosity and connectivity, contained more oxidized carbon. Pore network analysis indicated that soil compaction in the Cultivated soil may lead to longer diffusion pathways, significantly influencing SOM transport and stability.

Hydraulic Properties↗
DOE OSTI2021

Engineered Thin Diffusion Layers for Anion-Exchange Membrane Electrolyzer Cells with Outstanding Performance

Anion-exchange membrane electrolyzer cells (AEMECs) are one of the most promising technologies for carbon-neutral hydrogen production. Over the past few years, the performance and durability of AEMECs have substantially improved. Herein, we report an engineered liquid/gas diffusion layer (LGDL) with tunable pore morphologies that enables the high performance of AEMECs. The comparison with a commercial titanium foam in the electrolyzer indicated that the engineered LGDL with thin-flat and straight-pore structures significantly improved the interfacial contacts, mass transport, and activation of more reaction sites, leading to outstanding performance. We obtained a current density of 2.0 A/cm 2 at 1.80 V with an efficiency of up to 81.9% at 60 °C under 0.1 M NaOH-fed conditions. The as-achieved high performance in this study provides insight to design advanced LGDLs for the production of low-cost and high-efficiency AEMECs.

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

Porous transport electrodes for oxygen evolution reaction in proton exchange membrane water electrolysis -cells: Materials, designs, and diagnoses

H 2 production using proton exchange membrane (PEM) water electrolysis (PEMWE) cells has received considerable attention because of the high efficiencies of these cells and no harmful emissions from the related process. In PEMWE cells, porous transport electrodes (PTEs) composed of a catalyst layer (CL) comprising O 2 evolution reaction (OER) catalysts, porous transport layer (PTL), and PEM play key roles in the stack performance and lifetime. Herein, Ir-based and non-precious-metal OER catalysts that are highly active and stable at low pH values and high anodic potentials are reviewed to understand their OER mechanisms. Various strategies are proposed for engineering CLs and PTLs to improve the interfacial properties and mass transfers of reactants and products to and from the active sites. Additionally, diagnoses of PTEs is significantly crucial for interpreting electrochemical processes and addressing their current challenges. Therefore, half-cell analyses, including diffusion electrode (DE), floating electrode (FE), and modified rotating disk electrode (MRDE) techniques, are explored, and membrane electrode assembly (MEA)-based analyses, such as the polarization technique, electrochemical impedance spectroscopy, and magnetic field analysis, are established. In conclusion, this study aims to provide an overview of recent technologies used for the engineering and diagnostic tools of PEMWE cells and insights into the advanced components and systems to be developed in this field.

Diagnosis of PEMWE Cells↗
DOE OSTI2023

Gas nitriding behavior of refractory metals and implications for multi-principal element alloy design

Multi-principal element alloys (MPEAs) comprise a large, flexible compositional space that enables tuning of their chemistry, structure, and properties. Here, to facilitate the development of nitriding-based surface-enhancement strategies that harness a broad compositional space, this study examined the gas nitriding behavior of Hf, Mo, Nb, Ta, Ti, and Zr as a function of time, temperature (750 and 1000 °C), and nitriding potential (i.e. ammonia-to-hydrogen ratio). These metals were selected because they have a strong driving force to form nitrides, and appear in many promising refractory MPEA compositions. The nitriding temperatures were selected based on the phase transformation temperature of Ti and Zr, and the nitriding potentials were chosen such that all elements are expected to form nitrides. Mass gain measurements indicate that all six elements follow parabolic kinetics. The microstructure observations and quantitative microchemical analysis show formation of dense and well-adhered compound layers for Mo, Nb, and Ta. Thick diffusion zones appear in Hf, Ta, Ti, and Zr, and diffusion coefficients were fit to the composition profiles. Partial delamination of the compound layer occurred for Ti and Zr. Peak hardness values above 30 GPa are obtained in the dense compound layers, and the solute hardening of the underlying alloy is correlated with the nitrogen content. The results provide insight into the dynamics of nitride compound formation relative to interstitial dissolution of nitrogen, and are discussed in the context of MPEA composition and processing design.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Moisture-driven CO 2 pump for direct air capture

Direct Air Capture, the chemical separation of carbon dioxide (CO 2 ) from air, is considered a necessary approach to generate negative carbon emissions and in turn limit global warming, though energetic and monetary costs must continue to decline to meet the needed gigaton scale. The promise of energy-efficient and continuous carbon dioxide membrane separation has only recently been investigated with the development of facilitated chemical transport and novel electrochemically driven methods. This work demonstrates the potential for continuous pumping of CO 2 from air using commercial anion exchange membranes driven solely by a water vapor gradient. A reactive transport model that couples ionic transport with moisture sensitive bicarbonate chemistry in charged polymers, also termed moisture swing sorption, was developed alongside experimentally measured equilibrium and transport properties of CO 2 and water vapor in Fumasep FAA-3 to determine the mass transport limits of such a system. The results show that the commercial membrane is kinetically limited by the moisture sensitive chemistry, enabling a CO 2 current of 1.1 μmol/m 2 • s. In diffusion limited materials, a CO 2 pumping flux can reach 11 μmol/m 2 • s in and peak performance is expected with a dry air relative humidity around 50% pumping into a saturated water vapor permeate. In all configurations, water loss becomes the major cost in such a separation and must be managed through optimal polymer design or water recovery mechanisms. The work demonstrates the applicability of a new driving force, the evaporation of water, to meet the energetic demand of CO 2 separation from air.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tailoring the Physicochemical Properties of Nb Thin Films via Surface Engineering Methods

The modification of surface oxide layers formed on niobium (Nb) thin films via chemical mechanical planarization (CMP) and accelerated neutral atom beam (ANAB) processing provides a promising route toward tailoring their emergent properties and performance when used as superconducting qubits. Here, in this study, we show that CMP- and ANAB-formed Nb oxides are significantly thinner and smoother than the native oxide, as revealed by transmission electron microscopy (TEM) and atomic force microscopy. Scanning TEM and energy-dispersive X-ray spectroscopy along with X-ray photoelectron spectroscopy identified an oxidation gradient within the native and surface-engineered oxides. The topside layer is dominated by Nb 5+ (Nb 2 O 5 ), with various Nb suboxides present closer to the oxide/metal interface. Time-of-flight secondary ion mass spectrometry (ToF-SIMS) depth profiling confirmed the presence of an oxygen content gradient and demonstrated the enhanced resistance of the CMP- and ANAB-formed oxides to oxygen surface exchange and subsequent diffusion via 18 O 2 isotopic labeling experiments. ToF-SIMS also identified an interfacial layer containing trapped hydrogen (H)-containing species at the Nb oxide/metal interface. In situ ToF-SIMS and TEM revealed migration of the H/OH interfacial layer coinciding with decomposition of the surface oxide. Furthermore, our density functional theory calculations indicated that both H from moisture present in ambient air and bulk H in Nb films tend to segregate at the interface. These findings underscore the importance of understanding surface oxidation mechanisms, hydrogen incorporation, and their impact on the designed functionalities of Nb-based devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Observationally Constrained Modeling of the Reactive Uptake of Isoprene-Derived Epoxydiols under Elevated Relative Humidity and Varying Acidity of Seed Aerosol Conditions

Isoprene is the non-methane volatile organic compound (VOC) emitted in largest amounts to the atmosphere, and it is a significant source of secondary organic aerosol (SOA) mass. The uptake of isoprene oxidation products followed by multiphase chemistry in fine particles is the key pathway to form isoprene epoxydiol-derived SOA (IEPOX-SOA). However, many parameters that relate to diffusion and reaction of IEPOX in the particle phase remain uncertain, since reaction kinetics previously measured in bulk aqueous phase solutions might be different from atmospheric aerosols. Here, we use simultaneous environmental chamber measurements of multiple parameters governing IEPOX-SOA formation at timescales of ~hours: particle size distribution, composition, and volatility of IEPOX-SOA to constrain the key parameters governing IE-POX-SOA formation under humid (i.e., 50% relative humidity, RH) and varying seed aerosol acidity conditions. Reducing the 2-methyltetrol (tetrol) reaction rate constants by a factor of 4 brings the model predictions in agreement with the IEPOX-SOA measurements with acidified ammonium bisulfate seed aerosols. For less acidic ammonium sulfate aerosols both the organo-sulfate (OS) and tetrol reaction rate constants need to be reduced to bring model predictions closer to chamber observations. Using the measured non-volatile content of IEPOX-SOA we constrain the oligomerization timescale of 2-methyltetrols. We find that the oligomerization timescale is 4 hours with acidified seed aerosols, but a much longer time scale of 24 hours is needed for non-acidified seed aerosols, indicating that the aerosol acidity greatly affects the oligomerization rate of tetrols. Here we show that the actual kinetics of IEPOX-SOA formation rate on aerosol seeds consisting of both ammonium bisulfate and ammonium sulfate are a factor of 4~5 slower under 50-60% RH conditions compared to their application in previous models, which were based on bulk aqueous solution measurements.

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

Versatile soil gas concentration and isotope monitoring: optimization and integration of novel soil gas probes with online trace gas detection

Abstract. Gas concentrations and isotopic signatures can unveil microbial metabolisms and their responses to environmental changes in soil. Currently, few methods measure in situ soil trace gases such as the products of nitrogen and carbon cycling or volatile organic compounds (VOCs) that constrain microbial biochemical processes like nitrification, methanogenesis, respiration, and microbial communication. Versatile trace gas sampling systems that integrate soil probes with sensitive trace gas analyzers could fill this gap with in situ soil gas measurements that resolve spatial (centimeters) and temporal (minutes) patterns. We developed a system that integrates new porous and hydrophobic sintered polytetrafluoroethylene (sPTFE) diffusive soil gas probes that non-disruptively collect soil gas samples with a transfer system to direct gas from multiple probes to one or more central gas analyzer(s) such as laser and mass spectrometers. Here, we demonstrate the feasibility and versatility of this automated multiprobe system for soil gas measurements of isotopic ratios of nitrous oxide (δ18O, δ15N, and the 15N site preference of N2O), methane, carbon dioxide (δ13C), and VOCs. First, we used an inert silica matrix to challenge probe measurements under controlled gas conditions. By changing and controlling system flow parameters, including the probe flow rate, we optimized recovery of representative soil gas samples while reducing sampling artifacts on subsurface concentrations. Second, we used this system to provide a real-time window into the impact of environmental manipulation of irrigation and soil redox conditions on in situ N2O and VOC concentrations. Moreover, to reveal the dynamics in the stable isotope ratios of N2O (i.e., 14N14N16O, 14N15N16O, 15N14N16O, and 14N14N18O), we developed a new high-precision laser spectrometer with a reduced sample volume demand. Our integrated system – a tunable infrared laser direct absorption spectrometry (TILDAS) in parallel with Vocus proton transfer reaction mass spectrometry (PTR-MS), in line with sPTFE soil gas probes – successfully quantified isotopic signatures for N2O, CO2, and VOCs in real time as responses to changes in the dry–wetting cycle and redox conditions. Broadening the collection of trace gases that can be monitored in the subsurface is critical for monitoring biogeochemical cycles, ecosystem health, and management practices at scales relevant to the soil system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Calculation of self-diffusion coefficients in supercritical carbon dioxide using mean force kinetic theory

This paper presents an application of mean force kinetic theory (MFT) to the calculation of the self-diffusivity of CO 2 in the supercritical fluid regime. Two modifications to the typical application of MFT are employed to allow its application to a system of molecular species. Therefore, the first is the assumption that the inter-particle potential of mean force can be obtained from the molecule center-of-mass pair correlation function, which in the case of CO 2 is the C–C pair correlation function. The second is a new definition of the Enskog factor that describes the effect of correlations at the surface of the collision volume. The new definition retains the physical picture that this quantity represents a local density increase, resulting from particle correlations, relative to that in the zero density homogeneous fluid limit. These calculations are facilitated by the calculation of pair correlation functions from molecular dynamics (MD) simulations using the FEPM2 molecular CO2 model. The self-diffusivity calculated from theory is in good agreement with that from MD simulations up to and slightly beyond the density at the location of the Frenkel line. The calculation is compared with and is found to perform similarly well to other commonly used models but has a greater potential for application to systems of mixed species and to systems of particles with long range interatomic potentials due to electrostatic interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗