Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atmospheric physics parameters”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Detection of the keto-enol tautomerization in acetaldehyde, acetone, cyclohexanone, and methyl vinyl ketone with a novel VUV light source

The discovery of enols in combustion environments and our atmosphere has garnered increasing attention to the many unanswered questions surrounding enol chemistry. The scarcity of experimental data concerning these enols renders combustion and atmospheric models with a lack of constraining parameters, leading to varying computational predictions. Experimental detection is difficult because mass spectrometry, a powerful tool for probing a wide variety of species, cannot distinguish between enols and their thermodynamically favorable ketone isomers. A solution to this ambiguity is to use tunable vacuum ultraviolet (VUV) light from a synchrotron to identify the presence of the enol by its lower ionization energy compared to the isomer. We present a tabletop-scale VUV light source that implements highly cascaded harmonic generation, a new regime of cascaded nonlinear optics, to provide a set of spectral lines spaced by 1.2 eV. We demonstrate that the variety of photon energies available allows us to detect the keto-enol tautomerization of four aldehydes and ketones. By combining this novel VUV light source with an established microreactor, we first revisit the formation of vinyl alcohol from acetaldehyde and confirm that the observed isomerization is indeed unimolecular. Secondly, we observe the thermal tautomerization of acetone to propen-2-ol for the first time. Finally, we observe the thermal tautomerization of cyclohexanone to 1-cyclohexenol and methyl vinyl ketone to 2-hydroxybutadiene, where the results are in good agreement with those reported at a synchrotron. Furthermore, our measurements can be used to constrain models, inform future experimental studies of enol reactivity, and potentially enhance current understanding of combustion and environmental chemistry.

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

ATOMIC Simulations and Experimental Data for CaCO3 Mixtures

This data consists of simulations and experimental measurements of laser-induced breakdown spectroscopy (LIBS). The simulations are produced by ATOMIC, a general purpose plasma modeling and kinetics code that has been designed to compute emission (or absorption) spectra from plasmas [1] and are used to develop a statistical characterization of matrix effects. Our overall suite of simulations includes contains several sets of simulations: training and validation sets of simulations for three and four element mixtures of calcium, carbon, oxygen, and nitrogen (included to account for atmosphere) along with simulations of the individual elements. The 4-element simulations include the mixture of all four elements mentioned and for each of the four individual elements. The 3-element simulations include output for the mixture of calcium, carbon, oxygen and for these three individual elements. The training data were produced using a 600-run design, shown in Figure 1, that varies input parameters temperature (T), electron density (Ne), and proportion of the elements calcium, carbon, oxygen, and nitrogen (Ca; C; O; N) for the 4 element output. The 3-element output includes all parameters except for the proportion of nitrogen. The element proportions (all the variables but T and Ne) sum to one and are unused in the single-element simulations. The validation data was produced with a 80-run design shown in Figure 2. The training and validation simulation outputs for the 4-element simulations for the mixture and for the single element calcium are shown as sample simulations in Figures 3 and 4 respectively. The simulations produce spectra over a range of 190nm - 950nm that roughly mimics the range collected by the SciAps Z-300 LIBS instrument that was used for the experimental data. The measured spectra for a CaCO3 (which may include contribution from Earth's atmosphere) in the experiment is shown in in Figure 5. All files are kept in directories whose names indicate the elemental composition (CaCO3, Ca, C, O, or N), number of elements (3 or 4), and purpose (training, which is not labeled in the file name, or validation) with file names numbered to indicate the line in the design files used to produce the simulation. The designs are provided as text files with names indicating their purpose. The experimental data is provided as a CSV file. [1] J Colgan, EJ Judge, DP Kilcrease, and JE Barefield II. Ab-initio modeling of an iron laser-induced plasma: Comparison between theoretical and experimental atomic emission spectra. Spectrochimica Acta Part B: Atomic Spectroscopy, 97:65{73}, 2014.

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

Multi-scale dynamics of atmospheric-pressure discharges ignited over liquid electrodes

A multiphase computational model of atmospheric-pressure gas discharges ignited over liquid electrodes is developed. The model takes into account both the liquid electrode motion under applied electric fields and plasma generation in the gas phase. The influence of the applied voltage and liquid properties (density, viscosity, and surface tension) on the liquid surface deformation is analyzed. It is shown that the shape of dynamic cones formed on the liquid surface differs from the static Taylor's cones. The influence of the liquid surface protrusions on gas breakdown dynamics is demonstrated. It is shown that the breakdown develops in two stages: first, a fast ionization wave propagates from the cathode to the anode, and then the cathode sheath collapses due to secondary electron emission from the liquid surface. It is shown that dynamics of the multiphase system containing gas, plasma, and liquid states is characterized by three disparate time scales: the fast electron time scale, the intermediate ion time scale, and the slow liquid dynamics time scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Single-Particle Isotopic Analysis Using the Liquid Sampling-Atmospheric Pressure Glow Discharge Microplasma Coupled to an Orbitrap Mass Spectrometer

Isotope ratio (IR) determinations on individual particles can provide valuable information relative to the sample, including processing and dating, which could be valuable to geological and nuclear material analysis communities. In comparison to “bulk” measurements, in which all particle information is lost and homogenized within a sample, “particle” measurements allow for fingerprinting and can provide unique insights related to the sample’s nano- and micro- compositions. Furthermore, the complex nature of particles, with potential isobaric interferences from either the matrix or other elements within the same particle, poses a significant analytical challenge for conventional mass spectrometry platforms employed for elemental analysis due to their limited mass resolution. Presented here is the first demonstration of an ultrahigh resolution method for single particle (SP) isotopic analysis using the liquid sampling-atmospheric pressure glow discharge (LS-APGD) microplasma ionization source coupled to an Orbitrap mass spectrometer and FTMS Booster X2T acquisition/processing system. For proof of concept, well-characterized CeO 2 microparticles with a nominal diameter of 1.0 μm were analyzed at a mass resolution of ∼330,000. Important instrumental parameters were optimized to enable the detection of single particles. The 142 Ce/ 140 Ce ratio of individual particles and the population average were determined and compared to an SP inductively coupled plasma time-of-flight mass spectrometry (ICP-TOF-MS) analysis. The isotopic accuracy and precision were determined by two methods and found to be in good agreement with the SP-ICP-TOF-MS method, with the linear regression slope method providing a more accurate ratio, showing a ∼1.4% relative difference from the ratio determined from a particle digest. The encouraging performance of the method was further supported by a determined detection limit of 2.9 fg ( 142 Ce). The developed method is anticipated to overcome many of the challenges posed by isobaric interferences encountered in conventional mass spectrometry techniques when analyzing real-world samples for particle isotopic composition.

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

Impact of metal diffusion, lattice distortions, native defects, and ambient on dielectric breakdown in Ni–Ga 2 O 3 Schottky diodes

Ga 2 O 3 unipolar devices are of high interest due to their ~8 MV/cm predicted breakdown fields, which have not yet been achieved due to premature device failure. Pre- and post-failure defect analysis of Ni–Ga 2 O 3 Schottky diodes in ultrahigh vacuum (UHV) and air were performed using depth-resolved cathodoluminescence, high angle annular dark field scanning transmission electron microscopy, and energy dispersive x-ray analysis to understand the physical mechanisms that precede premature breakdown. The breakdown voltage in UHV was dramatically reduced by nearly 40% compared with the breakdown in air. This reduction in the breakdown voltage correlated with post-breakdown differences in Ni distribution, indicating that the coordination and bonding of Ni contribute strongly to electrical behavior in Ni–Ga 2 O 3 Schottky diodes. Breakdown studies in UHV revealed that Ni diffuses away more from the metal–semiconductor interface than with air breakdown, where Ni localizes more near the interface, indicative of the preferential formation of a Ni oxide under O-poor conditions. These measurements also identified the formation of divacancy-interstitial complexes and their characteristic luminescence signature ~150 nm from the interface, the densities of which correlated with breakdown fields. These findings show that electric-field-induced degradation occurs via the rearrangement of native point defects, which act as an additional precursor to device failure. Macroscopically, they show the impact of both vacuum conditions and metal reactivity on Ga 2 O 3 device fabrication.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Scale‐Adaptive Urban Hydrologic Framework: Incorporating Network‐Level Storm Drainage Pipes Representation

Abstract Below‐ground urban stormwater networks (BUSNs) significantly influence urban flood dynamics, yet their representation at the watershed or larger scales remains challenging. We introduce a scalable urban hydrologic framework that centers on a novel network‐level BUSN representation, balancing the needs for physical basis, parameter parsimony, and computational efficiency. Our framework conceptualizes an urban watershed into four interacting zones: hillslopes (natural), storm‐sewersheds (urban), a sub‐network channel (tributaries), and a main channel. We develop an innovative Graph Theory‐based algorithm to derive network‐level BUSN parameters from publicly available datasets, enabling efficient, scalable parameterization. We demonstrate this framework's applicability at nine representative watersheds in the Houston metropolitan region, USA, with urban imperviousness ranging from 0% to 64% and drainage areas ranging from 24 to 302 . Our model achieves satisfying computational efficiency, completing hourly time step simulations for 18 years in less than 5 sec per watershed on a standard PC. Validation against observed daily streamflow confirms that the model can capture small‐to‐large flood peaks and seasonal and annual water balance over these watersheds. Comparisons with the National Water Model show better performance in predicting flood peaks and overall water balance, underscoring the promises of our new framework for urban hydrologic modeling at large scales. Furthermore, analysis reveals nonlinear relationships between BUSNs' designed capacities and flood reduction effects. Our approach bridges the gap between detailed hydraulic and large‐scale hydrologic models, providing a valuable tool for urban flood prediction and management across broader spatial and temporal scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Estimating and Evaluating Roughness Length and Displacement Height in Heterogeneous Urban Environments

The roughness length (z 0 ) and displacement height (z d ) are essential surface-layer parameters in numerical models (e.g., weather, climate, wall-modeled LES, etc.). This work evaluates the consistency of z 0 and z d estimates from morphometric and anemometric methods using data from two eddy-covariance flux towers (AmeriFlux US-INg and US-INc) in Indianapolis, IN. Results show inconsistencies in estimated z 0 and z d values depending on the chosen method. The two evaluated anemometric methods estimate non-physical values of z d when compared to roughness elements surrounding both towers. Additionally, predictions of mean wind speed using surface-layer similarity theory with morphometric estimates exhibit a bias during near-neutral and stable conditions relative to observations. The overestimation of mean wind speed by surface layer similarity theory is consistent with previous observational and modeling studies in urban areas, suggesting that the application of similarity theories to urban environments may have limitations. Differentiation of vegetation from built structures appears to impact morphometric z 0 and z d estimates, particularly where vegetation is abundant; however, it has little impact on correcting biases in the similarity theory. Specifically, we find that existing similarity theories using morphometric estimates underestimate integral velocity and length scales, and the degree of underestimation depends on the stability conditions. Accounting for the degree of anisotropy in surface-layer turbulence helps reduce the biases between similarity theories and observations during unstable conditions, but not in near-neutral cases. Future work is needed to identify the cause of such biases for near-neutral conditions.

Aerodynamic roughness length↗
DOE OSTI2022

Isotopic discrimination of natural and anthropogenic perchlorate sources in groundwater in a semi-arid region of northeastern Oregon (USA)

Perchlorate (ClO 4 - ) has synthetic and natural sources. Synthetic ClO 4 - is released to the environment from its use as an oxidant in military and aerospace applications, and from its presence in a variety of common commercial products, such as safety flares, chlorate herbicides, and fireworks. Natural sources of ClO 4 - in the environment include imported nitrate fertilizers derived from salt deposits in the Atacama Desert of Chile and indigenous natural ClO 4 - that accumulates in unsaturated soils and groundwaters in other arid and semi-arid environments, largely from atmospheric deposition. The stable isotope ratios of chlorine ( 37 Cl/ 35 Cl) and oxygen ( 18 O/ 16 O, 17 O/ 16 O) and the isotopic abundance of radioactive 36 Cl in ClO 4 - can be used to discriminate these different sources. Perchlorate was previously detected at relatively high concentrations (3.8–34.7 μg/L) in groundwater from many wells in the Boardman-Umatilla area near the Columbia River in northeastern Oregon, which is a semi-arid, highly agricultural, heavily irrigated area that includes several past and current military installations. Eight representative groundwater wells were sampled throughout this region and isotopic characteristics of ClO 4 - collected from each well were measured along with other chemical and isotopic parameters including tritium and other groundwater age indicators. We report isotopic data indicate that indigenous natural ClO 4 - was present in groundwater from all sampled wells and was the predominant source in five of the wells. Synthetic ClO 4 - was present in the three remaining wells with natural ClO 4 - , and a minor fraction of Atacama-fertilizer-derived ClO4- was indicated in one of the wells. Data from this study expand the geographic area of the USA in which indigenous natural ClO 4 - has been detected to include the semi-arid northwest. This study also illustrates the role of irrigation recharge as a mechanism for producing relatively high concentrations of indigenous natural ClO 4 - in groundwater by flushing accumulated salts from the unsaturated zone.

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

Modification of the Loss Cone for Energetic Particles in the Earth's Inner Magnetosphere

The behavior of energetic particles in the Earth's dipole and stretched magnetic-field models is explored with a focus on the shift of the atmospheric loss cone away from the magnetic-field direction and on non-adiabatic behavior occurring when the particle gyroradius becomes comparable with the gradient scale length of the local magnetic field. It is shown that the equatorial loss cone is aberrated away from the magnetic-field direction (pitch angle of 0°) by the perpendicular drift of a charged particle (which is referred to as the “Mozer transform” described in Mozer (1966, https://doi.org/10.1029/JZ071i011p02701). It is found that the Mozer coordinate transformation in pitch-angle/gyrophase-angle space better organizes the behavior of bounce times, mirror altitudes, drift speeds, and first adiabatic invariant. It also describes the loss-cone shift accurately for a certain range of values of the adiabaticity parameter ϵ, defined by the ratio of the particle gyroradius to the local radius of curvature of the magnetic field. For particles with larger ϵ, the Mozer transform (evaluated with the standard gradient-curvature drift) breaks down and the “central trajectory” theory can be used to calculate the angular shift of the loss cone, which now includes an Earthward shift. When ϵ becomes relatively large, stochastic field line curvature (FLC) scattering occurs: we show the intimate connection between the strength of FLC scattering and the loss-cone shift. Finally, by comparing ion orbits in the dipole and Tsyganenko (Ts89 and Ts04) magnetic fields, we conclude that the loss cone of ring-current ions is significantly modified during geomagnetically active times.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A theoretical investigation of the hydrolysis of uranium hexafluoride: the initiation mechanism and vibrational spectroscopy

Depleted uranium hexafluoride (UF 6 ), a stockpiled byproduct of the nuclear fuel cycle, reacts readily with atmospheric humidity, but the mechanism is poorly understood. Here we compare several potential initiation steps at a consistent level of theory, generating underlying structures and vibrational modes using hybrid density functional theory (DFT) and computing relative energies of stationary points with double-hybrid (DH) DFT. A benchmark comparison is performed to assess the quality of DH-DFT data using reference energy differences obtained using a complete-basis-limit coupled-cluster (CC) composite method. The associated large-basis CC computations were enabled by a new general-purpose pseudopotential capability implemented as part of this work. Dispersion-corrected parameter-free DH-DFT methods, namely PBE0-DH-D3(BJ) and PBE-QIDH-D3(BJ), provided mean unsigned errors within chemical accuracy (1 kcal mol -1 ) for a set of barrier heights corresponding to the most energetically favorable initiation steps. The hydrolysis mechanism is found to proceed via intermolecular hydrogen transfer within van der Waals complexes involving UF 6 , UF 5 OH, and UOF 4 , in agreement with previous studies, followed by the formation of a previously unappreciated dihydroxide intermediate, UF 4 (OH) 2 . The dihydroxide is predicted to form under both kinetic and thermodynamic control, and, unlike the alternate pathway leading to the UO 2 F 2 monomer, its reaction energy is exothermic, in agreement with observation. Finally, harmonic and anharmonic vibrational simulations are performed to reinterpret literature infrared spectroscopy in light of this newly identified species.

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

The Curious Case of [AnH(NR 2 ) 3 ] (An = Th, U; R = SiMe 3 ): Two Monomeric Actinide Hydrides Revisited

The reaction of AnCl 4 (DME) x (An = Th, x = 2; An = U, x = 0) with 4 equiv of NaNR 2 (R = SiMe 3 ) in THF at 65 °C results in the formation of [An{N(R)(SiMe 2 CH 2 )}(NR 2 ) 2 ] (An = U, 1; An = Th, 2), and not the reported monomeric actinide hydrides, [AnH(NR 2 ) 3 ], as expected. Furthermore, both complexes 1 and 2 were characterized by X-ray crystallography. Surprisingly, their unit cell parameters are remarkably close to those reported for [AnH(NR 2 ) 3 ], suggesting that the original crystals of [AnH(NR 2 ) 3 ] were, in fact, [An{N(R)(SiMe 2 CH 2 )}(NR 2 ) 2 ], but were misidentified. Reduction of 1 with 1.1 equiv of KC 8 in THF, in the presence of 1 equiv of 2.2.2-cryptand, results in the formation of [K(2.2.2-cryptand)][U{N(R)(SiMe 2 CH 2 )}(NR 2 ) 2 ] (3) in good yield. Likewise, the reaction of 1 with 1 equiv of bis(diisopropylamino)cyclopropenylidene (BAC) results in the formation of the BAC adduct, [(BAC)U{N(R)(SiMe 2 CH 2 )}(NR 2 ) 2 ] (4), in moderate yield. Finally, the addition of H 2 (10 bar) to 2 in C 6 D 6 at room temperature results in the formation of the targeted monomeric hydride, [ThH(NR 2 ) 3 ], in 32% yield, according to integrations against an internal standard. However, removal of the H 2 atmosphere results in rapid reformation of 2. In contrast, the addition of H 2 (10 bar) to 1 in C 6 D 6 at room temperature results in no apparent reaction.

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

Search for Heavy Neutral Leptons at the MINER$\nu$A detector

Heavy Neutral Leptons (HNL) are particles hypothesised to provide a mass generation mechanism for the active (observed) neutrino species, which are known to have nonzero mass from the definitive observation of neutrino oscillations. HNL are eigenstates of mass of the order $\mathcal{O}(0.1 − 1\,\,\mathrm{GeV}/c^{2})$, which mix into the active flavour eigenstates through the extended leptonic mixing matrix. Apart from neutrino mass, they could provide a natural dark matter candidate and a mechanism for matter-antimatter asymmetry in the early Universe, giving rise to today’s matter-dominated cosmos. Searches for HNL typically attempt to either confirm the existence of HNL through an excess in data that is most compatible with an HNL hypothesis, or by setting limits on the HNL parameter space $\left(M_{N4}, \left|U_{\alpha 4}\right|^{2}\right)$ in the case of statistically insignificant excess. Such searches are intensifying as part of a global research programme at both colliders and accelerator / atmospheric neutrino experiments. One such setting is the MINER$\nu$A experiment, located in the NuMI beamline at Fermilab. MINER$\nu$A has collected a large amount of data over seven years of operation to measure the cross-sections of neutrino-nucleus interactions necessary to drive systematic uncertainties down, in order for neutrino oscillation experiments to achieve sensitivity to CP violation in the neutrino sector. With the high-energy, high-intensity NuMI beam, and with good timing, position, and angular resolution leading to sensitivity to HNL decays, MINER$\nu$A is fertile ground for an HNL search, which is performed in this thesis. Specifically, a novel, experiment-agnostic and general simulation of HNL production and decay is presented and deployed; this simulation has been incorporated in the ubiquitous GENIE neutrino event generator for use with neutrino experiments in the future. The event selection and background characterisation is discussed in detail, including the primary background coming from charged-current coherent and diffractive pion production from neutrino-nucleus interactions, and constraints on the background using control regions from MINER$\nu$A data are derived. Finally, by means of a fake-data study given a background of $\mathcal{O}(500)$ events, the discovery potential and limit-setting capacity of MINER$\nu$A is demonstrated, and ways to improve this capacity are expounded upon.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗
DOE OSTI2023

Techno-economic analysis of carbon dioxide capture from low concentration sources using membranes

Rising carbon dioxide (CO 2 ) levels in the atmosphere lead to global warming, causing climate change. As such, carbon capture has become necessary to slow the increase and reduce CO 2 levels in the atmosphere. Point source emissions have a wide range of CO 2 concentrations, but emissions below 3% CO 2 have mostly been ignored because Carbon capture from these sources has been viewed as costly and economically unsustainable. Membrane technologies are considered the most viable solution by virtue of more energy-efficient operation. Our group at Idaho National Laboratory (INL) has developed poly[bis((2-methoxyethoxy)ethoxy)phosphazene] (MEEP)-based carbon dioxide selective membranes with CO 2 /N 2 selectivity greater than 40 and CO 2 permeability greater than 450 Barrer. To understand the economics of carbon capture, a spreadsheet-based techno-economic analysis (TEA) model was developed to consider multiple parameters, including selectivity and permeability of the membranes, performance conditions such as the number of stages, module material, electricity price, membrane price, and capital financing. The cost of carbon capture in US $\$$/metric ton was calculated at various purities and compared with other membrane processes, cryogenic capture, solvent-based capture, and pressure swing adsorption-based capture. It was determined that a MEEP-based three-stage process had a capture cost of US $\$$ 50.1/metric ton for 99.8% purity CO 2 from a 1% CO 2 feed source in nitrogen (N 2 ). In conclusion, the capture cost using the best performing Pebax-based membrane was 464% higher, cryogenic capture was 60%–140% higher, pressure swing adsorption was 55%–165% higher, and chemical absorption was -10%–110% higher than MEEP-based membrane capture, respectively.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Combined Meteorological and Hydrologic Uncertainties Shape Projections of Future Soil Moisture in the Eastern United States

Physical hazards pose risks to many critical systems. Designing adaptive measures to mitigate these risks is challenging due to large uncertainties in modeling future hazards and the associated sectoral responses. Here, we help address this challenge in a hydrologic context by examining the combined role of meteorological forcing and hydrologic parameter uncertainties in shaping projections of future soil moisture. By encoding a simple conceptual water balance model in a differentiable programming framework, we facilitate fast runtimes and an efficient calibration, enabling an improved uncertainty analysis. We characterize uncertainty in model parameters by calibrating against different target data sets and by using several loss functions. We then convolve the resulting parameter ensemble with a set of Earth system model projections to produce a large ensemble (2,340 members) of daily soil moisture simulations. Focusing on the eastern United States, we find that most ensemble members project a drying of soils across the region, although some simulate wetter conditions throughout this century. Our ensemble shows an increase in the frequency and intensity of dry extremes while there is less agreement for wet extremes. We conduct sensitivity analyses on several soil moisture signatures to measure the relative influence of meteorological and hydrologic uncertainties across space and time. Both meteorological and hydrologic factors contribute consistently to uncertainty surrounding long-term trends, while changes to both wet and dry soil extremes are typically more sensitive to hydrologic parameter uncertainty. Our results underscore the need to account for varied sources of uncertainty when developing long-term hydrometeorological projections.

Lafferty, David C. [University of Illinois Urbana‐↗
DOE OSTI2020

Characterization of an RF-driven argon plasma at atmospheric pressure using broadband absorption and optical emission spectroscopy

Atmospheric pressure plasmas in argon are of particular interest due to the production of highly excited and reactive species enabling numerous plasma-aided applications. In this contribution, we report on absolute optical emission and absorption spectroscopy of a radio frequency (RF) driven capacitively coupled argon glow discharge operated in a parallel-plate configuration. This enabled the study of all key parameters including electron density and temperature, gas temperature, and absolute densities of atoms in highly electronically excited states. The space and time-averaged electron density and temperature were determined from the measurement of the absolute intensity of the electron-atom bremsstrahlung in the visible range. Considering the non-Maxwellian electron energy distribution function, an electron temperature (T e ) of 2.1 eV and an electron density (n e ) of 1.1 × 10 19 m –3 were obtained. The time-averaged and spatially resolved absolute densities of atoms in the metastable (1s 5 and 1s 3 ) and resonant (1s 4 and 1s 2 ) states of argon in pure Ar and Ar/He mixture were obtained by broadband absorption spectroscopy. The 1s 5 metastable atoms had the largest density near the sheath region with a maximum value of 8 × 10 17 m –3 , while all other 1s states had densities of at most 2 × 10 17 m –3 . Here, the dominant production and loss mechanisms of these atoms were discussed, in particular, the role of radiation trapping. We conclude with a comparison of the plasma properties of the argon RF glow discharges with the more common He equivalent and highlight their differences.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗
DOE OSTI2025

Combining molecular beam epitaxy and low-energy electron microscopy with in situ magnetic susceptibility measurements within an integrated ultrahigh vacuum system

Quantum two-dimensional materials, including ultrathin superconducting films, are of great current research interest. These films are typically fabricated under ultra-high vacuum (UHV) conditions and are sensitive to the environment—prone to oxidation and contamination when exposed to the atmosphere. This hampers the study of their intrinsic properties by standard ex situ techniques. Here, we present a variable-temperature mutual inductance probe system integrated under UHV with molecular beam epitaxy (MBE) synthesis and low-energy electron microscopy, enabling nondestructive in situ characterization of superconducting thin films. The system employs a reflection-type configuration and reaches a low temperature (∼4 K) using a high-cooling-power, vibration-isolated cryocooler. In conclusion, we demonstrate the system performance by measuring the superconducting critical temperature in a copper-oxide thin film.

2D materials↗