Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “well data”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

The Cretaceous-Tertiary Impact Crater and the Cosmic Projectile that Produced it

Evidence gathered to date from topographic data, geophysical data, well logs, and drill-core samples indicates that the buried Chicxulub basin, the source crater for the approximately 65 Ma Cretaceous-Tertiary (K/T) boundary deposits, is approximately 300 km in diameter. A prominent topographic ridge and a ring of gravity anomalies mark the position of the basin rim at approximately 150 km from the center. Wells in this region recovered thick sequences of impact-generated breccias at 200-300 m below present sea level. Inside the rim, which has been severely modified by erosion following impact, the subsurface basin continues to deepen until near the center it is approximately 1 km deep. The best planetary analog for this crater appears to be the 270 km-diameter Mead basin on Venus. Seismic reflection data indicate that the central zone of downward displacement and excavation (the transient crater is approximately 130 km in diameter, consistent with previous studies of gravity anomaly data). Our analysis of projectile characteristics utilizes this information, coupled with conventional scaling relationships, and geochemical constraints on the mass of extraterrestrial material deposited within the K/T boundary layer. Results indicate that the Chicxulub crater would most likely be formed by a long-period comet composed primarily of nonsilicate materials (ice, hydrocarbons, etc.) and subordinate amounts (less than or equal to 50 percent) primitive chondritic material. This collision would have released the energy equivalent to between 4 x 10(exp 8) and 4 x 10(exp 9) megatons of TNT. Studies of terrestrial impact rates suggest that such an event would have a mean production rate of approximately 1.25 x 10(exp -9) y(exp -1). This rate is considerably lower than that of the major mass extinctions over the last 250 million years (approximately 5 x 10(exp -7) y(exp -1). Consequently, while there is substantial circumstantial evidence establishing the cause-effect link between the Chicxulub basin forming event and the K/T biological extinctions, the results of our analysis do not support models of impact as a common or singular causative agent of mass extinctions on Earth.

Sharpton, Virgil L.↗
NASA NTRS2021

An Evaluation of Remotely Sensed and In-Situ Data Sufficiency for SGMA-Scale Groundwater Studies in the Central Valley, California

California’s Central Valley aquifer is a critical freshwater resource for the state, providing drinking water to 6.5 million residents and irrigation water for more than half of the nation’s produce. Increasing demands have led to a steady decline of aquifer water levels, leading to irreversible compaction and land subsidence. The Sustainable Groundwater Management Act of California (SGMA), instituted in 2014, requires high-priority regions to establish sustainability plans by 2020-2022 that mitigate groundwater depletion and land subsidence. Many regions within the Central Valley lack groundwater monitoring wells or continuous well data, rendering in-situ monitoring challenging. When combined with other data, NASA’s Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow On (GRACE-FO) missions provide valuable information about groundwater storage changes at a sub-basin scale. Interferometric Surface Aperture Radar (InSAR) measurements can map land subsidence, and Global Positioning System (GPS) can be used to estimate crustal uplift. Yet, all of these measurements have differences in resolution, coverage, discretization, and record length. Here we assess the various datasets that are potentially useful for Central Valley sub-basin groundwater monitoring and provide an analysis of gaps and other issues in regard to their suitability for SGMA-related analysis. Finally, we offer next steps and recommendations on data treatment and integration.

Kyra H Kim↗
NASA NTRS2004

Progressive Aerodynamic Model Identification From Dynamic Water Tunnel Test of the F-16XL Aircraft

Development of a general aerodynamic model that is adequate for predicting the forces and moments in the nonlinear and unsteady portions of the flight envelope has not been accomplished to a satisfactory degree. Predicting aerodynamic response during arbitrary motion of an aircraft over the complete flight envelope requires further development of the mathematical model and the associated methods for ground-based testing in order to allow identification of the model. In this study, a general nonlinear unsteady aerodynamic model is presented, followed by a summary of a linear modeling methodology that includes test and identification methods, and then a progressive series of steps suggesting a roadmap to develop a general nonlinear methodology that defines modeling, testing, and identification methods. Initial steps of the general methodology were applied to static and oscillatory test data to identify rolling-moment coefficient. Static measurements uncovered complicated dependencies of the aerodynamic coefficient on angle of attack and sideslip in the stall region making it difficult to find a simple analytical expression for the measurement data. In order to assess the effect of sideslip on the damping and unsteady terms, oscillatory tests in roll were conducted at different values of an initial offset in sideslip. Candidate runs for analyses were selected where higher order harmonics were required for the model and where in-phase and out-of-phase components varied with frequency. From these results it was found that only data in the angle-of-attack range of 35 degrees to 37.5 degrees met these requirements. From the limited results it was observed that the identified models fit the data well and both the damping-in-roll and the unsteady term gain are decreasing with increasing sideslip and motion amplitude. Limited similarity between parameter values in the nonlinear model and the linear model suggest that identifiability of parameters in both terms may be a problem. However, the proposed methodology can still be used with careful experiment design and carefully selected values of angle of attack, sideslip, amplitude, and frequency of the oscillatory data.

Murphy, Patrick C.↗
NASA NTRS2012

Staring at 4U 1909+07 with Suzaku (Research Note)

We present an analysis of the neutron star High Mass X-ray Binary (HMXB) 4U 1909+07 mainly based on Suzaku data. We extend the pulse period evolution, which behaves in a random-walk like manner, indicative of direct wind accretion. Studying the spectral properties of 4U 1909+07 between 0.5 to 90keV we find that a power-law with an exponential cutoff can describe the data well, when additionally allowing for a blackbody or a partially covering absorber at low energies. We find no evidence for a cyclotron resonant scattering feature (CRSF), a feature seen in many other neutron star HMXBs sources. By performing pulse phase resolved spectroscopy we investigate the origin of the strong energy dependence of the pulse profile, which evolves from a broad two-peak profile at low energies to a profile with a single, narrow peak at energies above 20keV. Our data show that it is very likely that a higher folding energy in the high energy peak is responsible for this behavior. This in turn leads to the assumption that we observe the two magnetic poles and their respective accretion columns at different phases, and that these accretions column have slightly different physical conditions.

Fuerst, F.↗
NASA NTRS1999

Time course of sleep inertia dissipation in human performance and alertness

Alertness and performance on a wide variety of tasks are impaired immediately upon waking from sleep due to sleep inertia, which has been found to dissipate in an asymptotic manner following waketime. It has been suggested that behavioural or environmental factors, as well as sleep stage at awakening, may affect the severity of sleep inertia. In order to determine the time course of sleep inertia dissipation under normal entrained conditions, subjective alertness and cognitive throughput were measured during the first 4 h after habitual waketime from a full 8-h sleep episode on 3 consecutive days. We investigated whether this time course was affected by either sleep stage at awakening or behavioural/environmental factors. Sleep inertia dissipated in an asymptotic manner and took 2-4 h to near the asymptote. Saturating exponential functions fitted the sleep inertia data well, with time constants of 0.67 h for subjective alertness and 1.17 h for cognitive performance. Most awakenings occurred out of stage rapid eye movement (REM), 2 or 1 sleep, and no effect of sleep stage at awakening on either the severity of sleep inertia or the time course of its dissipation could be detected. Subjective alertness and cognitive throughput were significantly impaired upon awakening regardless of whether subjects got out of bed, ate breakfast, showered and were exposed to ordinary indoor room light (approximately 150 lux) or whether subjects participated in a constant routine (CR) protocol in which they remained in bed, ate small hourly snacks and were exposed to very dim light (10-15 lux). These findings allow for the refinement of models of alertness and performance, and have important implications for the scheduling of work immediately upon awakening in many occupational settings.

Non-NASA Center↗
DOE OSTI2023

Long-term sinking of nuclear waste canisters in salt formations by low-stress creep at high temperature

Abstract Rock salt has a self-sealing capacity, low permeability, and high thermal conductivity, making it a potential host for heat-generating nuclear waste. The feasibility of nuclear waste disposal within salt formations has been investigated mostly for small-sized canisters. Geologic disposal of larger-sized canisters originally designed for spent fuel storage and transportation has lately been examined as a cost-effective alternative. This raises questions about their long-term vertical movement due to their weight and high decay heat. Low-stress creep governs this movement; however, most salt constitutive models do not incorporate it. In this paper, the Norton and the WIPP creep models are compared with the Lux/Wolters/Lerche (LWL) model and a simpler model that combines linear and Norton creep laws (named combined creep model). The LWL and combined creep models consider pressure solution creep, though all incorporate dislocation creep. The models are first applied to creep tests under various stress levels. The LWL and the combined creep models results fit the experimental data well in both high and low stress ranges, whereas the Norton and WIPP models results only fit in higher-stress ranges. The different models are further applied for analyzing long-term canister movement. A sinking rate of $$-4.4\times 10^{-7}$$ - 4.4 × 10 - 7 mm/year was predicted using the Norton and WIPP models versus $$-2.1\times 10^{-2}$$ - 2.1 × 10 - 2 mm/year and $$-3.1\times 10^{-2}$$ - 3.1 × 10 - 2 mm/year using the LWL and the combined creep models, respectively. This comparative study confirms that creep models calibrated exclusively against high-deviatoric stress data might result in an inaccurate estimation of waste packages sinking rate in salt formations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

BayFlux: A Bayesian method to quantify metabolic Fluxes and their uncertainty at the genome scale

Metabolic fluxes, the number of metabolites traversing each biochemical reaction in a cell per unit time, are crucial for assessing and understanding cell function. 13 C Metabolic Flux Analysis ( 13 C MFA) is considered to be the gold standard for measuring metabolic fluxes. 13 C MFA typically works by leveraging extracellular exchange fluxes as well as data from 13 C labeling experiments to calculate the flux profile which best fit the data for a small, central carbon, metabolic model. However, the nonlinear nature of the 13 C MFA fitting procedure means that several flux profiles fit the experimental data within the experimental error, and traditional optimization methods offer only a partial or skewed picture, especially in “non-gaussian” situations where multiple very distinct flux regions fit the data equally well. Here, we present a method for flux space sampling through Bayesian inference (BayFlux), that identifies the full distribution of fluxes compatible with experimental data for a comprehensive genome-scale model. This Bayesian approach allows us to accurately quantify uncertainty in calculated fluxes. We also find that, surprisingly, the genome-scale model of metabolism produces narrower flux distributions (reduced uncertainty) than the small core metabolic models traditionally used in 13 C MFA. The different results for some reactions when using genome-scale models vs core metabolic models advise caution in assuming strong inferences from 13 C MFA since the results may depend significantly on the completeness of the model used. Based on BayFlux, we developed and evaluated novel methods (P- 13 C MOMA and P- 13 C ROOM) to predict the biological results of a gene knockout, that improve on the traditional MOMA and ROOM methods by quantifying prediction uncertainty.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1996

Diffusion lengths in irradiated N/P InP-on-Si solar cells

Indium phosphide (InP) solar cells were made on silicon (Si) wafers (InP/Si) by to take advantage of both the radiation-hardness properties of the InP solar cell and the light weight and low cost of Si wafers. The InP/Si cell application is for long duration and/or high radiation orbit space missions. Spire has made N/P InP/Si cells of sizes up to 2 cm by 4 cm with beginning-of-life (BOL) AM0 efficiencies over 13% (one-sun, 28C). These InP/Si cells have higher absolute efficiency and power density after a high radiation dose than gallium arsenide (GaAs) or silicon (Si) solar cells after a fluence of about 2e15 1 MeV electrons/sq. cm. In this work, we investigate the minority carrier (electron) base diffusion lengths in the N/P InP/Si cells. A quantum efficiency model was constructed for a 12% BOL AM0 N/P InP/Si cell which agreed well with the absolutely measured quantum efficiency and the sun-simulator measured AM0 photocurrent (30.1 mA/sq. cm). This model was then used to generate a table of AM0 photocurrents for a range of base diffusion lengths. AM0 photocurrents were then measured for irradiations up to 7.7e16 1 MeV electrons/sq. cm (the 12% BOL cell was 8% after the final irradiation). By comparing the measured photocurrents with the predicted photocurrents, base diffusion lengths were assigned at each fluence level. A damage coefficient K of 4e-8 and a starting (unirradiated) base electron diffusion length of 0.8 microns fits the data well. The quantum efficiency was measured again at the end of the experiment to verify that the photocurrent predicted by the model (25.5 mA/sq. cm) agreed with the simulator-measured photocurrent after irradiation (25.7 mA/sq. cm).

Wojtczuk, Steven↗
NASA NTRS1995

Determination of plasma temperatures and luminosities using multiple extreme-ultraviolet and X-ray filters

We carefully examine the techniques used to infer temperatures of stellar coronal plasmas from the count rates of several broadband instruments in the X-ray and extreme-ultraviolet spectral ranges. In particular, we determine to what extent temperatures can be constrained and the corresponding uncertainties in the luminosities and emission measures lowered by fitting simultaneously count rates from the Einstein imaging proportional counter (IPC), the ROSAT Position Sensitive Proportional Counter (PSPC), the ROSAT Wide Field Camera (WFC) (both filters), and the EXOSAT Low Energy Telescope (LET) with the 3-Lex filter. We use published plasma emissivities with solar photospheric abundances. Since it has been found that single-temperature plasmas do not fit IPC data well, we assume a two-temperature plasma model. We find that, even with count rates from all of the above filters and overly optimistic error estimates, it is still not possible to determine a unique two-temperature solution. However, since the use of count rates from many filters can reduce substantially the number of possible solutions, temperature solutions determined by other means can be tested. We carry out such an analysis on a set of 18 nearby late-type stars to determine possible two-temperature solutions using multifilter photometry, and we compare these results with the temperature solutions derived by Schmitt et al. (1990) using IPC spectral data. In general, the two-temperature fits derived from the IPC spectral data are inconsistent with our results, with our data implying that, for many stars, the two temperatures derived by the IPC may be too low by about a factor of 2. The EXOSAT transmission grating Spectrometer (TGS) spectra of capella and sigma(exp 2) CrB support this conclusion. For Procyon and 70 Oph, though, the presence of a temperature component cooler than a million degress (not detected by the IPC) is deduced. While our analysis suggests the existence of more than one temperature in the coronae of late-type stars, in many instances our WFC data appear to be inconsistent with the presence of significant emission measure over a broad temperature distribution. This, together with the success of two-temperature plasmas in fitting IPC and TGS data, implies that for many stars, the coronal emission measure distribution may in fact be dominated by two distinct temperature regimes.

Wood, Brian E.↗
DOE OSTI2021

Jet fragmentation transverse momentum distributions in pp and p-Pb collisions at $ \sqrt{s}$, $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV

Jet fragmentation transverse momentum ( j T ) distributions are measured in proton-proton (pp) and proton-lead (p-Pb) collisions at s NN = 5 . 02 TeV with the ALICE experiment at the LHC. Jets are reconstructed with the ALICE tracking detectors and electromagnetic calorimeter using the anti- k T algorithm with resolution parameter R = 0 . 4 in the pseudorapidity range |η| < 0 . 25. The j T values are calculated for charged particles inside a fixed cone with a radius R = 0 . 4 around the reconstructed jet axis. The measured j T distributions are compared with a variety of parton-shower models. Herwig and P ythia 8 based models describe the data well for the higher j T region, while they underestimate the lower j T region. The j T distributions are further characterised by fitting them with a function composed of an inverse gamma function for higher j T values (called the “wide component”), related to the perturbative component of the fragmentation process, and with a Gaussian for lower j T values (called the “narrow component”), predominantly connected to the hadronisation process. The width of the Gaussian has only a weak dependence on jet transverse momentum, while that of the inverse gamma function increases with increasing jet transverse momentum. For the narrow component, the measured trends are successfully described by all models except for Herwig. For the wide component, Herwig and PYTHIA 8 based models slightly underestimate the data for the higher jet transverse momentum region. These measurements set constraints on models of jet fragmentation and hadronisation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

X-Ray Polarization of the Magnetar 1E 1841−045

We report on IXPE and NuSTAR observations beginning 40 days after the 2024 outburst onset of magnetar 1E 1841−045, marking the first IXPE observation of a magnetar in an enhanced state. Our spectropolarimetric analysis indicates that both a blackbody (BB) plus double power-law (PL) and a double blackbody plus power-law spectral model fit the phase-averaged intensity data well, with a hard PL tail (Γ = 1.19 and 1.35, respectively) dominating above ≈5 keV. For the former model, we find the soft PL (the dominant component at soft energies) exhibits a polarization degree (PD) of ≈30% while the hard PL displays a PD of ≈40%. Similarly, the cool BB of the 2BB+PL model possesses a PD of ≈15% and a hard PL PD of ≈57%. For both models, each component has a polarization angle (PA) compatible with celestial north. Model-independent polarization analysis supports these results, wherein the PD increases from ≈15% to ≈70% in the 2–3 keV and 6–8 keV ranges, respectively, while the PA remains nearly constant. We find marginal evidence for phase-dependent variability of the polarization properties, namely a higher PD at phases coinciding with the hard X-ray pulse peak. We compare the hard X-ray PL to the expectation from resonant inverse Compton scattering (RICS) and secondary pair cascade synchrotron radiation from primary high-energy RICS photons; both present reasonable spectropolarimetric agreement with the data, albeit the latter does so more naturally. We suggest that the soft PL X-ray component may originate from a Comptonized corona in the inner magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRSDate not supplied

Experimental and Theoretical Needs for the JWST Early Release Science Program on Radiative Feedback from Massive Stars (PDRs4ALL): III. A Quantum Chemistry Perspective

Massive stars disrupt their natal molecular cloud material by dissociating molecules, ionizing atoms and molecules, and heating the gas and dust. These processes drive the evolution of interstellar matter in our Galaxy and throughout the Universe from the era of vigorous star formation at z=1-3 to the present day. Much of this interaction occurs in Photo- Dissociation Regions (PDRs) where far-ultraviolet photons of these stars create a largely neutral, but warm region of gas and dust. PDR emission dominates the IR spectra of starforming galaxies and provides a unique tool to study the physical and chemical processes that are relevant for inter- and circumstellar media including diffuse clouds, molecular cloud and protoplanetary disk surfaces, and starburst galaxies. The ERS program ID1288 is dedicated to provide template data as well as data processing and analysis tools for PDRs. To this end, it will observe the Orion Bar, the proto-typical PDR situated in the nearby Orion Nebula, using NIRSpec IFU, MIRI IFU, and NIRCAM and MIRI imaging. These observations will, for the first time, spatially resolve and perform a tomography of the PDR, revealing the individual IR spectral signatures from the key zones and sub-regions within the ionized gas, the PDR, and the molecular cloud. Quantum chemistry can be used to determine molecular data such as ionization potentials, chemical reaction rates, electronic spectra, dissociation energies, channels, and branching ratios, etc. These are often used in synergy with experiment to take advantage of the strengths of each approach. An area of active research is modeling the infrared (IR) emission spectra produced by PAH molecules, including the effects of anharmonicity. This will be important in the interpretation of JWST data which will yield higher-resolution spectral data compared to Spitzer or ISO. Quantum chemistry methods can compute accurate anharmonic IR absorption spectra of PAH molecules up to about 25 C atoms, and the results have been benchmarked against high-resolution experimental spectra. These feed into a quantum Monte Carlo (QMC) approach to generate temperature-dependent anharmonic IR emission spectra and then used to generate fully anharmonic cascade emission IR spectra of PAH molecules -- the very type of spectra that JWST will observe. Work continues to extend these methods to larger PAH molecules. We will discuss the current state of this work, as well as

Timothy Lee↗
DOE OSTI2022

Bayesian inference of fiber orientation and polymer properties in short fiber-reinforced polymer composites

Herein we present a Bayesian methodology to infer the elastic modulus of the constituent polymer and the fiber orientation state in a short-fiber reinforced polymer composite (SFRP). The properties are inversely determined using only a few experimental tests. Developing composite manufacturing digital twins for SFRP composite processes, including injection molding and extrusion deposition additive manufacturing (EDAM) requires extensive experimental material characterization. In particular, characterizing the composite mechanical properties is time consuming and therefore, micromechanics models are used to fully identify the elasticity tensor. Hence, the objective of this paper is to infer the fiber orientation and the effective polymer modulus and therefore, identify the elasticity tensor of the composite with minimal experimental tests. To that end, we develop a hierarchical Bayesian model coupled with a micromechanics model to infer the fiber orientation and the polymer elastic modulus simultaneously which we then use to estimate the composite elasticity tensor. We motivate and demonstrate the methodology for the EDAM process but the development is such that it is applicable to other SFRP composites processed via other methods. Our results demonstrate that the approach provides a reliable framework for the inference, with as few as three tensile tests, while accounting for epistemic and aleatory uncertainty. Posterior predictive checks show that the model is able to recreate the experimental data well. The ability of the Bayesian approach to calibrate the material properties and its associated uncertainties, make it a promising tool for enabling a probabilistic predictive framework for composites manufacturing digital twins.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Relationship between ductility and the porosity of additively manufactured AlSi10Mg

Additive manufacturing via selective laser melting can result in variable levels of internal porosity both between build plates and within components from the same build. In this work, sample porosity levels were compared to tensile properties for 176 samples spanning eight different build plates. Sample porosity was measured both by Archimedes density, which provided an estimation of overall porosity, and by observation of voids in the fracture surface, which provided an estimation of the porosity at the failure plane. The porosity observed at the fracture surface consistently demonstrated higher porosity than that suggested by Archimedes density. The porosity values obtained from both methods were compared against the mechanical results. Sample porosity appears to have some correlation to the ultimate tensile strength, yield strength, and modulus, but the strongest relationship is observed between porosity and ductility. Three different models were used to relate the fracture surface porosity to the ductility. The first method was a simple linear regression analysis, while the other two models have been used to relate porosity to ductility in cast alloys. It is shown that all three models fit the data well over the observed porosity ranges, suggesting that the models taken from casting theory can extend to additively manufactured metals. Finally, it is proposed that the non-destructive Archimedes method could be used to estimate an approximate sample ductility through the use of correlations realized here. Such a relationship could prove useful for design and for a deeper understanding of the impact of pores on tensile behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Oil in Alaska North Slope Gas Hydrate Reservoir: Micro-CT and Flow Simulation Insights into Permeability

Gas hydrate-bearing sands on the Alaska North Slope (ANS) host minor volumes of crude oil whose impact on formation permeability has never been quantified relative to a gas hydrate reservoir system. Here, we combine in situ pressure-core microcomputed-tomography (μ-CT), thermogravimetric analysis (TGA), gas-chromatography–mass-spectrometry (GC–MS), scanning-electron microscopy coupled with energy-dispersive spectroscopy (SEM-EDS), and Stokes-flow simulation to (i) quantify oil saturation, (ii) infer its source, and (iii) evaluate its influence on permeability after gas-hydrate dissociation. μ-CT and SEM-EDS imaging identify the host sediment as silty, with some grain-coating clay present. Eight representative 3003-voxel subvolumes extracted from a preserved pressure core (Core 15P-3c, from the B1 sand (unit B), HYDRATE 02 Geo Data Well (GDW), 886.96–887.29 m measured depth, MD) exhibit porosities of 0.34–0.37 and an average oil saturation of 0.09 ± 0.03, in comparison to an independent TGA-based estimate of 0.16. GC-MS fingerprints obtained from an adjacent depressurized core (Core 17P-1, 891.37–891.44 m MD), together with oil pore habit revealed by μ-CT, suggest that the oil is partly native but also includes contributions from mineral oil-based drilling fluid contamination. Further analysis of oil–porewater interactions and flushing tests confirms that native oil saturation in the analyzed samples ranges from 0.04 to 0.08 with an average of approximately 0.06, and the oil is highly immobile. Flow simulations demonstrate that the native oil saturation of ≈0.06 reduces permeability by approximately 50% due to both pore blockage and increased flow-path tortuosity. These findings reveal the potential presence of native oil within the ANS gas hydrate reservoir and suggest that future simulation models may need to account for its impact on permeability to improve long-term performance predictions of gas and water production.

02 PETROLEUM↗
DOE OSTI2026

CO 2 Sorption in Moisture Swing Anion Exchange Resins for Direct Air Capture: Experimental Isotherm Determination and Modeling

Moisture swing (MS) sorption is a promising direct air capture (DAC) technology to achieve negative CO 2 emissions and counter global warming. Here, in this study, MS CO 2 sorption in a model MS sorbent, IRA900, was investigated. IRA900 is a macroporous commercial strong-base anion exchange resin (AER) with quaternary ammonium functional groups. A rigorous CO 2 desorption process was developed to desorb all CO 2 from the sample and obtain CO 2 sorption isotherms as a function of relative humidity, CO 2 partial pressure, and temperature. CO 2 sorption increased with decreasing relative humidity, decreasing temperature, and increasing CO 2 partial pressure. Remarkably, the total CO 2 uptake from the sorption isotherms matched the ion exchange capacity (IEC) of the material, suggesting a stoichiometry of one CO 2 molecule reacting per active site. An isotherm model for CO 2 sorption starting with empty sorption sites (i.e., with the AER in the OH – form) described the experimental data well, supporting the hypothesis that the bicarbonate-loaded AER fully unloaded to the OH – state following rigorous CO 2 desorption. CO 2 sorption kinetics were studied, and carbon diffusion coefficients were estimated as a function of %RH and CO 2 partial pressure. Overall, this study provides a roadmap for the systematic evaluation of MS CO 2 DAC materials.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Magnetoelectric Coupling in a Mn 4 Na–Organic Complex under Pulsed Magnetic Fields up to 73 T

Research on the magnetoelectric (ME) effect (or spin-electric coupling) in molecule-based magnetic materials is a relatively nascent but promising topic. Molecule-based magnetic materials have diverse magnetic functionalities that can be coupled to electrical properties. Here, in this study, we investigate a realization of ME coupling that is fundamental but not heavily studied─the coupling of magnetic spin level crossings to changes in electric polarization. A mixed-valence Mn 4 Na complex with a total ground-state spin S = 5/2 under zero magnetic field and S = 17/2 under high magnetic field undergoes a cascade of ground-state level crossings of the S z states with increasing magnetic field. Magnetization and electrical polarization measurements under pulsed magnetic fields up to 73 T show that each spin level crossing is accompanied by a significant change in electric polarization that is an even function of the applied magnetic field. A molecular Hamiltonian describing antiferromagnetic exchange in a distorted tetrahedron of three Mn III and one Mn II ions matches the data well. We conclude that the ME coupling is caused by magnetostriction within the polar molecule as it distorts to lower its magnetic exchange energy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

On analysis and stochastic modeling of the particle kinetic energy equation in particle-laden isotropic turbulent flows

In this study, we analyze three-dimensional particle-laden, isotropic turbulence to develop an understanding of inertial particle dynamics from a kinetic energy perspective. Data trends implying inhomogeneous sampling of the flow by particles are identified and used to support a proposed particle behavior: particles appear to accumulate in regions of low flow kinetic energy over time because they lose kinetic energy and slow down in such regions, ultimately causing them to spend more time there. To elucidate this behavior, we derive a particle kinetic energy equation from the particle momentum equation, which incorporates inertial effects through the Schiller–Naumann drag correlation. Upon extracting fundamental physics from this equation, hypotheses regarding the role of the Stokes number in the temporal change of particle kinetic energy and the previously proposed particle behavior are evaluated using simulation data considering three Stokes numbers. Finally, a Fokker–Planck equation is used to derive the steady-state probability density function of the particle kinetic energy. The model fits the simulation data well and provides a tool for further investigation into understanding preferential concentration, as well as a reduced order model for predicting particle kinetic energy in turbulent flows.

42 ENGINEERING↗