Search NASA⌕ Search

SEARCH · Search NASA

Results for “Regional surrogate”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Neural-Network-Enhanced COTSIM: Advancing Predictive Capabilities for Fast DIII-D Simulations

Sustaining fusion reactions in tokamaks requires heating plasma to thermonuclear temperatures while maintaining confinement and stability. Neutral beam injection (NBI) provides heating, current drive, torque, and fueling, while electron cyclotron (EC) waves are widely used for heating and current drive; together, these actuators shape the plasma current, temperature, and density profiles. The control-oriented tokamak simulator (COTSIM), a predictive, control-oriented code, has been enhanced with neural-network surrogates for transport and sources. Turbulent transport is predicted by MMMnet—a neural-network version of the updated multimode model (MMM 9.0.10)—with significantly reduced computation time relative to MMM; neoclassical transport follows the Chang–Hinton model. NUBEAMnet, a surrogate of the Monte Carlo NUBEAM module, predicts beam-driven heating, current, and torque. EC heating and current drive use a control-oriented, empirically scaled source model; plasma resistivity follows the Spitzer formulation; bootstrap current uses the Sauter model. Equilibrium is computed using both prescribed and fixed-boundary solvers (FBSs), and the pedestal structure is modeled with an empirical pedestal model. For a representative DIII-D discharge, COTSIM predicts electron and ion temperature and safety-factor profiles in close agreement with TRANSP predictive and interpretive simulations while extending predictions through the pedestal region to the plasma edge (versus 80% of the minor radius in TRANSP). Furthermore, the equivalent COTSIM simulation runs in under 3 min compared to about 2 h for TRANSP, enabling rapid scenario planning, optimization of tokamak operation, and between-pulse control design.

Control-oriented tokamak simulator (COTSIM)↗
DOE OSTI2024

EPCAPE-PT-LANL Measurements: Wideband Integrated Bioaerosol Sensor

Coastal cities offer a unique environment for studying aerosol-cloud interactions and the effects of urban emissions on cloud properties. As part of the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE), the Partitioning Thrust by Los Alamos National Laboratory (EPCAPE-PT-LANL) was conducted. Our campaign focused on measuring the optical and chemical properties of aerosols and their interactions within marine stratocumulus clouds in La Jolla, California. EPCAPE-PT-LANL enhances the primary goals of EPCAPE through innovative observations of vapor-phase transitions between aerosols and cloud droplets, the impact of black carbon on aerosol-cloud dynamics, and the effects of cloud processing on aerosol optical properties. Instrument: Wideband Integrated Bioaerosol Sensor (Droplet Measurements Technology) Data Notes: The WIBS is an online single-particle measurement that detects FBAPs (within a size range of 0.5 - 30 microns in diameter) based on the excitation and emission wavelengths of the individual particles. Using two xenon lamps, the WIBS excites FBAPs at 280 nm and 370 nm. Their emission is detected across two wavebands of 310-400 nm and 420-650 nm. We classified the FBAPs into seven different categories (A, B, C, AB, BC, AC, and ABC) using the classification scheme in Perring et. al. (2015) [1]. Averaged number concentration of FBAPs (total and by category) and particles that non-fluorescent bioaerosols particles (NFBAPs). In separate files, we also present one-minute-averaged size distributions and the asymmetry factor (AF, a surrogate for shape) of all FBAPs and NFBAPs. The logarithmic bin width of the size bins are the same as the average bin width of the AOS's optical particle counter (OPC, Grimm) for the range of sizes in which they overlap (26 bins from 0.5 - 30 microns). AF of the particles ranges from 0-100 and is divided into five bins with a linear spacing at increments of 20. The smallest AF bin represents more spherical particles while the largest bin represents more rod-shaped particles. [1] Perring, A. E., et al. (2015), Airborne observations of regional variation in fluorescent aerosol across the United States, J. Geophys. Res. Atmos., 120, 1153–1170, doi:10.1002/2014JD022495. Abstract and description of the campaign can be found here : https://www.arm.gov/research/campaigns/amf2023epcape-pt-lanl. Files data_10min_WIBS_AFDist.csv Header: - FBAP_AFDist[/cm3]_Bin_1 to Bin_5: Concentration of fluorescent bioaerosol particles in the each of 5 AF bins, measured in particles per cubic centimeter. Each bin represents a specific range of particle AF, capturing the shapes of FBAPs detected during the measurement. - NFBAP_AFDist[/cm3]_Bin_1 to Bin_5: Concentration of non-fluorescent bioaerosol particles in the each of 5 AF bins, measured in particles per cubic centimeter. Each bin represents a specific range of particle AF, capturing the shapes of NFBAPs detected during the measurement. - CVI_Flag[bool]: A boolean flag indicating whether the Counterflow Virtual Impactor (CVI) was active (true) or inactive (false) during the measurement. AF Bins: • Bin 1: 0 – 20 [unitless] • Bin 2: 21 – 40 [unitless] • Bin 3: 41 – 60 [unitless] • Bin 4: 61 – 80 [unitless] • Bin 5: 81 – 100 [unitless] Files data_10min_WIBS_Conc.csv Header: - NumberConcentrationA[/cm3]: Number concentration of bioaerosol particles detected by fluorescence channel A, measured in particles per cubic centimeter. - NumberConcentrationB[/cm3]: Number concentration of bioaerosol particles detected by fluorescence channel B, measured in particles per cubic centimeter. - NumberConcentrationC[/cm3]: Number concentration of bioaerosol particles detected by fluorescence channel C, measured in particles per cubic centimeter. - NumberConcentrationAB[/cm3]: Combined number concentration of bioaerosol particles detected by both fluorescence channels A and B, measured in particles per cubic centimeter. - NumberConcentrationBC[/cm3]: Combined number concentration of bioaerosol particles detected by both fluorescence channels B and C, measured in particles per cubic centimeter. - NumberConcentrationAC[/cm3]: Combined number concentration of bioaerosol particles detected by both fluorescence channels A and C, measured in particles per cubic centimeter. - NumberConcentrationABC[/cm3]: Combined number concentration of bioaerosol particles detected by all three fluorescence channels A, B, and C, measured in particles per cubic centimeter. - NumberConcentrationNFBAP[/cm3]: Number concentration of non-fluorescent bioaerosol particles, measured in particles per cubic centimeter. - CVI_Flag[bool]: A boolean flag indicating whether the Counterflow Virtual Impactor (CVI) was active (true) or inactive (false) during the measurement. Files data_10min_WIBS_SizeDist.csv Header: - FBAP_SizeDist[/cm3]_Bin_1 to FBAP_SizeDist[/cm3]_Bin_26: Number concentrations of FBAP in each of 26 size bins, measured in particles per cubic centimeter. Each bin represents a specific range of particle sizes, capturing the size distribution of FBAPs detected during the measurement. - NFBAP_SizeDist[/cm3]_Bin_1 to NFBAP_SizeDist[/cm3]_Bin_26: Number concentrations of NFBAP in each of 26 size bins, measured in particles per cubic centimeter. Similar to FBAP, each bin covers a specific range of particle sizes, detailing the size distribution of NFBAPs detected. - CVI_Flag[bool]: A boolean flag indicating whether the Counterflow Virtual Impactor (CVI) was active (true) or inactive (false) during the measurement. Size Bins: • Bin 1: 0.48 to 0.57 μm • Bin 2: 0.57 to 0.67 μm • Bin 3: 0.67 to 0.79 μm • Bin 4: 0.79 to 0.93 μm • Bin 5: 0.93 to 1.1 μm • Bin 6: 1.1 to 1.29 μm • Bin 7: 1.29 to 1.52 μm • Bin 8: 1.52 to 1.8 μm • Bin 9: 1.8 to 2.11 μm • Bin 10: 2.11 to 2.5 μm • Bin 11: 2.5 to 2.94 μm • Bin 12: 2.94 to 3.46 μm • Bin 13: 3.46 to 4.08 μm • Bin 14: 4.08 to 4.81 μm • Bin 15: 4.81 to 5.67 μm • Bin 16: 5.67 to 6.68 μm • Bin 17: 6.68 to 7.88 μm • Bin 18: 7.88 to 9.29 μm • Bin 19: 9.29 to 10.96 μm • Bin 20: 10.96 to 12.92 μm • Bin 21: 12.92 to 15.23 μm • Bin 22: 15.23 to 17.96 μm • Bin 23: 17.96 to 21.17 μm • Bin 24: 21.17 to 24.96 μm • Bin 25: 24.96 to 29.43 μm • Bin 26: 29.43 to 34.70 μm

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Unexpected trends in the enhanced Ce 3+ surface concentration in ceria–zirconia catalyst materials

Despite the immense importance of ceria–zirconia solid solutions in heterogeneous catalysis, and the growing consensus that catalytic activity correlates with the concentration of reduced Ce 3+ species and accompanying oxygen vacancies, the extent of reduction at the surfaces of these materials, where catalysis occurs, is unknown. Using angle-resolved X-ray Absorption Near Edge Spectroscopy (XANES), we quantify under technologically relevant conditions the Ce 3+ concentration in the surface (2–3 nm) and bulk regions of ceria–zirconia films grown on single crystal yttria-stabilized zirconia, YSZ (001). In all circumstances, we observe substantial Ce 3+ enrichment at the surface relative to the bulk. Surprisingly, the degree of enhancement is highest in the absence of Zr. This behavior stands in direct contrast to that of the bulk in which the Ce 3+ concentration monotonically increases with increasing Zr content. Overall, these results suggest that while Zr enhances the oxygen storage capacity in ceria, undoped ceria may have higher surface catalytic activity. They further urge caution in the use of bulk properties as surrogate descriptors for surface characteristics and hence catalytic activity.

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

A scale-wise analysis of intermittent momentum transport in dense canopy flows

We investigate the intermittent dynamics of momentum transport and its underlying time scales in the near-wall region of the neutrally stratified atmospheric boundary layer in the presence of a vegetation canopy. This is achieved through an empirical analysis of the persistence time scales (periods between successive zero-crossings) of momentum flux events, and their connection to the ejection–sweep cycle. Using high-frequency measurements from the GoAmazon campaign, spanning multiple heights within and above a dense canopy, the analysis suggests that, when the persistence time scales ( $t_p$ ) of momentum flux events from four different quadrants are separately normalized by $\varGamma _{w}$ (integral time scale of the vertical velocity), their distributions $P(t_p/\varGamma _{w})$ remain height-invariant. This result points to a persistent memory imposed by canopy-induced coherent structures, and to their role as an efficient momentum-transporting mechanism between the canopy airspace and the region immediately above. Moreover, $P(t_p/\varGamma _{w})$ exhibits a power-law scaling at times $t_{p}<\varGamma _{w}$ , with an exponential tail appearing for $t_{p} \geq \varGamma _{w}$ . By separating the flux events based on $t_p$ , we discover that around 80 % of the momentum is transported through the long-lived events ( $t_{p} \geq \varGamma _{w}$ ) at heights immediately above the canopy, while the short-lived ones ( $t_{p} < \varGamma _{w}$ ) only contribute marginally ( $\approx 20\,\%$ ). To explain the role of instantaneous flux amplitudes in momentum transport, we compare the measurements with newly developed surrogate data and establish that the range of time scales involved with amplitude variations in the fluxes tends to increase as one transitions from within to above the canopy.

Mechanics↗
DOE OSTI2021

HYPPO: A Surrogate-Based Multi-Level Parallelism Tool for Hyperparameter Optimization

We present a new software, HYPPO, that enables the automatic tuning of hyperparameters of various deep learning (DL) models. Unlike other hyperparameter optimization (HPO) methods, HYPPO uses adaptive surrogate models and directly accounts for uncertainty in model predictions to find accurate and reliable models that make robust predictions. Using asynchronous nested parallelism, we are able to significantly alleviate the computational burden of training complex architectures and quantifying the uncertainty. HYPPO is implemented in Python and can be used with both TensorFlow and PyTorch libraries. We demonstrate various software features on time-series prediction and image classification problems as well as a scientific application in computed tomography image reconstruction. Finally, we show that (1) we can reduce by an order of magnitude the number of evaluations necessary to find the most optimal region in the hyperparameter space and (2) we can reduce by two orders of magnitude the throughput for such HPO process to complete.

adaptation models↗
DOE OSTI2026

Characterizing the impact of water uptake on light absorption by aerosol particles

Absorbing aerosols (AA) have an important impact on the global radiation budget. The composition and properties of AA can vary substantially throughout the atmosphere, depending on the particle source and the influence of chemical aging. Uncertainties associated with the radiative effects of AA remain substantial. A key contributor to this uncertainty is understanding the extent to which water uptake alters absorption by AA particles and how this depends on particle composition. We have used new and existing experimental tools to systematically characterize the relationship(s) between AA chemical and physical characteristics, relative humidity, and mixing of AA with other aerosol components through a combination of laboratory and field measurements. We have developed fundamental understanding through laboratory experiments that will enable interpretation of field observations at DOE ARM sites and facilitate process-based improvements in the simulation of absorption by carbonaceous aerosols in climate models. In the laboratory experiments we have considered a variety of different AA types, including fullerene soot (a surrogate for refractory black carbon, or BC) and nigrosine (a surrogate for partially soluble brown carbon, or BrC). By examining a range of AA types we will develop insight into how the chemical and physical properties of AA particles (e.g. whether they are refractory or soluble) impact the influence of water uptake on absorption. These AA particles will also be mixed—both internally and externally—with compounds having a wide range of hygroscopic properties to establish how the impact of water uptake on absorption will vary with photochemical processing. These laboratory observations have been compared with commonly used theoretical models, enabling development of empirical models that will facilitate improved representation of absorption by AA particles in regional and global climate models. Overall, our study has contributed to the mission of the Atmospheric System Research program by quantifying how interactions between aerosols and radiation depend on relative humidity. Through this, our work will improve understanding and model representation of aerosol processes as they affect the Earth’s radiation budget.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

sMF-BO-2CoGP: A sequential multi-fidelity constrained Bayesian optimization framework for design applications

Bayesian optimization (BO) is an effective surrogate-based method that has been widely used to optimize simulation-based applications. While the traditional Bayesian optimization approach only applies to single-fidelity models, many realistic applications provide multiple levels of fidelity with various levels of computational complexity and predictive capability. In this work, we propose a multi-fidelity Bayesian optimization method for design applications with both known and unknown constraints. The proposed framework, called sMF-BO-2CoGP, is built on a multi-level CoKriging method to predict the objective function. An external binary classifier, which we approximate using a separate CoKriging model, is used to distinguish between feasible and infeasible regions. Finally, the sMF-BO-2CoGP method is demonstrated using a series of analytical examples and a flip-chip application for design optimization to minimize the deformation due to warping under thermal loading conditions.

42 ENGINEERING↗
DOE OSTI2025

Predicting interface structure using the minima hopping method

Here, we adapt the minima hopping method (MHM) to the problem of interfacial structure prediction and apply it to study a canonical problem, the tilt grain boundaries in SrTiO 3 . Our method employs a hybrid approach by first exploring the potential energy surface (PES) of different grain boundary samplings with an empirical force field, among which the fifteen candidates with lower energies are then refined using ab initio density functional theory (DFT) calculations. During the exploratory stage, we bias the search using a local order parameter to primarily sample various reconstructions in the vicinity of the interface, while preserving the crystallinity of the bulk regions. We further enhance the search by incorporating initial structures with rigid body displacements to account for translational variations between bulk phases, enabling the MHM to effectively generate both stoichiometric and nonstoichiometric SrTiO 3 Σ⁢3(111)[110] and Σ⁢3(112)[110] grain boundaries. From an algorithmic standpoint, MHM outperforms earlier studies based on genetic algorithms (GA) by identifying more stable interfacial structures of several SrTiO 3 grain boundaries. The performance of the present implementation of the MHM approach is primarily limited by exploring an approximate description of the PES with a rather simple Buckingham potential. This limitation leads to variations in performance when compared to approaches utilizing more advanced surrogate PES models, such as direct DFT-PES sampling or GA with the embedded atom method (EAM). Despite the present limitations, the MHM approach is able to yield interfacial structures with comparable or lower interfacial energies in specific cases, such as Σ⁢3(111)[110] Γ=1, ±0.5 and Σ⁢3(112)[110] Γ= ±1, −2, underscoring the robustness of the MHM approach even with a simple approximation of the DFT PES. The MHM interfacial structure prediction method thus offers an efficient approach to understanding the grain boundaries and heterointerfaces at the atomic scale, providing an important prerequisite for effective materials design.

density functional theory↗
DOE OSTI2023

Black and gray box learning of amplitude equations: Application to phase field systems

Here, we present a data-driven approach to learning surrogate models for amplitude equations and illustrate its application to interfacial dynamics of phase field systems. In particular, we demonstrate learning effective partial differential equations describing the evolution of phase field interfaces from full phase field data. We illustrate this on a model phase field system, where analytical approximate equations for the dynamics of the phase field interface (a higher-order eikonal equation and its approximation, the Kardar-Parisi-Zhang equation) are known. For this system, we discuss data-driven approaches for the identification of equations that accurately describe the front interface dynamics. When the analytical approximate models mentioned above become inaccurate, as we move beyond the region of validity of the underlying assumptions, the data-driven equations outperform them. In these regimes, going beyond black box identification, we explore different approaches to learning data-driven corrections to the analytically approximate models, leading to effective gray box partial differential equations.

42 ENGINEERING↗
DOE OSTI2025

Rapid spatial analysis of surrogate TRISO fuel particles using laser-induced breakdown spectroscopy image analysis

Laser-induced breakdown spectroscopy (LIBS) mapping enables rapid elemental and spatial analysis of solid samples. In this work, surrogate TRi-structural ISOtropic (TRISO) particles with various layers consisting of Zr, W, C, and Si were used to demonstrate the use of a new thickness measurement tool developed to analyze elemental images generated from LIBS maps. Zr particles with varying outer-layer thicknesses ranging from 16 to 32 μm were measured with scanning electron microscopy with energy dispersive spectroscopy (SEM-EDS) and LIBS with both a complementary metal-oxide-semiconductor (CMOS) detector and an intensified charged coupled device (ICCD) detector. LIBS maps of particles were completed using CMOS and ICCD spectrometers with effective spatial resolutions of 4 and 2 μm, respectively. The novel thickness measurement tool identified layer regions within a LIBS map and then identified the locations of the boundaries between these layers. The tool then generated up to 1000 random profiles stemming radially from the center region, which were used to measure the layer thickness/radius. This image analysis tool demonstrated LIBS's ability to provide values comparable with SEM-EDS (3.7% relative difference) along with a 95% reduction in measurement time. Furthermore, the precision of these measurements was on par with the SEM-EDS measurements at <15% relative standard deviation. Following the analysis of the Zr particles, W particles with increased complexity (e.g., five versus three layers) were analyzed using both the CMOS and ICCD spectrometers simultaneously. The spatial dimensions were extracted with an average relative difference of 2.7% and an average relative standard deviation of 9%.

Andrews, Hunter B. [Oak Ridge National Laboratory ↗
DOE OSTI2023

Bioaerosol Properties using WIBS from SAIL on Mt. Crested Butte, CO

We deployed a Wideband Integrated Bioaerosol Sensor (WIBS-Neo, DMT) to measure different properties of fluorescent bioaerosol particles (FBAP) during SAIL. Our measurements were collected at the supplementary site on Crested Butte Mountain from inside the ARM Aerosol Observing System (AOS) sampling from the AOS inlet. We had two deployment periods at SAIL during biologically active seasons: summer 2022 (15 June 2022 - 13 September 2022) and spring 2023 (21 March 2023 - 19 June 2023). The WIBS is an online single-particle measurement that detects FBAPs (within a size range of 0.5 - 30 microns in diameter) based on the excitation and emission wavelengths of the individual particles. Using two xenon lamps, the WIBS excites FBAPs at 280 nm and 370 nm. Their emission is detected across two wavebands of 310-400 nm and 420-650 nm. We classified the FBAPs into seven different categories (A, B, C, AB, BC, AC, and ABC) using the classification scheme in Perring et. al. (2015) [1]. Here, we report one-minute averaged number concentration of FBAPs (total and by category) and NFBAPs. We also present one-minute-averaged size distributions and the asymmetry factor (AF, a surrogate for shape) of all FBAPs and NFBAPs. The logarithmic bin width of the size bins is the same as the average bin width of the AOS's optical particle counter (OPC, Grimm) for the range of sizes in which they overlap (26 bins from 0.5 - 30 microns). AF of the particles ranges from 0-100 and is divided into five bins with a linear spacing at increments of 20. The smallest AF bin represents more spherical particles while the largest bin represents more rod-shaped particles. We recommend referring to the readme file for more details of the variables present in the files. [1] Perring, A. E., et al. (2015), Airborne observations of regional variation in fluorescent aerosol across the United States, J. Geophys. Res. Atmos., 120, 1153-1170, doi:10.1002/2014JD022495. Abstract and description of the campaign can be found here: https://www.arm.gov/research/cam

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Spatial Optimization of Multiscale Biorefinery Deployment for a Diversified Bioeconomy in the United States

Strategic biorefinery siting is critical for a diversified bioeconomy, yet industry, policy, and research often focus on either large-scale biofuel plants or smaller-scale specialty bioproduct facilities, with limited coordination across scales. We address this gap by modeling biorefinery deployment spanning a 28-fold difference in capacity. We developed an open-source, spatially explicit framework integrating techno-economic analysis with logistics and refinery cost surrogate models to evaluate multiscale miscanthus-derived biorefineries across the rainfed U.S. for the production of ethanol, succinic acid, lactic acid, potassium sorbate, and acrylic acid. Overall costs change little as feedstock density increases, while transport distances decrease by ∼30 to 67% (∼100 km) and siting flexibility improves. Specifically, a 5-fold feedstock density increase (2% to 10% of suitable land) reduces minimum selling prices by <10% (e.g., 0.27 USD·gal –1 for ethanol). This limited economic sensitivity suggests dense planting is not required for competitive deployment, particularly for smaller-scale facilities. Representing collection areas as irregular rather than circular expands the feasible space under low-density scenarios. While large-scale refineries anchor regional supply chains, smaller facilities retain spatial flexibility even when large refineries are established. These findings highlight the importance of spatial representation and multiscale coordination for robust, regionally tailored biomanufacturing networks to advance renewable carbon integration without extensive land conversion.

biorefinery siting↗
DOE OSTI2023

Explainability and extrapolation of machine learning models for predicting the glass transition temperature of polymers

Abstract Machine learning (ML) offers promising tools to develop surrogate models for polymers' structure–property relations. Surrogate models can be built upon existing polymer data and are useful for rapidly predicting the properties of unknown polymers. The accuracy of such ML models appears to depend on the feature space representation of polymers, the range of training data, and learning algorithms. Here, we establish connections between these factors for predicting the glass transition temperature (T g ) of polymers. Our analysis suggests linear models with fewer fitting parameters are as accurate as nonlinear models with many hidden and unexplainable parameters. Also, the performance of a monomer topology‐based ML model is found to be qualitatively identical to that of a physicochemical descriptor‐based ML model. We find that the ML models's performance in the extrapolative region is enhanced as the property range of the training data increases. Moreover, we establish newT g – polymer chemistry correlations via ML. Our work illustrates how ML can advance the fundamental understanding of polymer structure–property correlations and its efficacy for extrapolation problems.

Polymer Science↗
DOE OSTI2023

Structure of a 28.5 kDa duplex-embedded G-quadruplex system resolved to 7.4 Å resolution with cryo-EM

Genomic regions with high guanine content can fold into non-B form DNA four-stranded structures known as G-quadruplexes (G4s). Extensive in vivo investigations have revealed that promoter G4s are transcriptional regulators. Little structural information exists for these G4s embedded within duplexes, their presumed genomic environment. Here, we report the 7.4 Å resolution structure and dynamics of a 28.5 kDa duplex-G4-duplex (DGD) model system using cryo-EM, molecular dynamics, and small-angle X-ray scattering (SAXS) studies. The DGD cryo-EM refined model features a 53° bend induced by a stacked duplex-G4 interaction at the 5’ G-tetrad interface with a persistently unstacked 3’ duplex. The surrogate complement poly dT loop preferably stacks onto the 3’ G-tetrad interface resulting in occlusion of both 5’ and 3’ tetrad interfaces. Structural analysis shows that the DGD model is quantifiably more druggable than the monomeric G4 structure alone and represents a new structural drug target. Our results illustrate how the integration of cryo-EM, MD, and SAXS can reveal complementary detailed static and dynamic structural information on DNA G4 systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Unveiling stimulation fluid-driven alterations in shale pore architecture through combined interpretation of TD-NMR and multi-component gas adsorption

The hydraulic fracturing stimulation process can alter petrophysical and flow properties of the shale matrix in unconventional reservoirs. These alterations can result from interactions with native and exogenous fluids, and can be detrimental to reservoir productivity. Different experimental procedures have been proposed to study reactive processes and thereby understand the aforementioned interactions. Most published works do not consider the rock initial saturation with formation fluid, or the effect of organic matter (OM). To gain insights into these effects, two reactive experimental sets, differing in experimental conditions, were conducted. In the first set, at 125°C and 45 MPa, a stimulation fluid surrogate was put in contact with Baxter Shale samples, previously saturated and equilibrated with synthetic formation fluid. A baseline consisted of exposure to formation fluid exclusively. In the second set, at 0.08 MPa and 25°C, contact with stimulation fluid took place as in the first set. Two additional assays in this second set either skipped the initial formation fluid saturation or removed the OM, respectively. Petrophysical properties analyses, focused on the pore architecture, relied on the combined interpretation of gas adsorption and Time-Domain Nuclear Magnetic Resonance (TD-NMR). Results show that the presence of the formation fluid significantly restricts pore space accessibility to the stimulation fluid and mitigates the pore-architecture alteration to the macro-pore region. Mineral dissolution and precipitation dominate over surface or wettability alteration. Next, results suggest that the OM restricts accessibility to some pores. Ultimately, the equilibration process with synthetic formation fluid drives mainly surface alteration, unhindered by OM.

58 GEOSCIENCES↗
DOE OSTI2026

Transit Rider/Travel Behavior Inventory Survey - Minneapolis-St. Paul Metro - 1990

The 1990 transit on-board survey aimed to update the 1988 survey, which was conducted as part of the Preliminary Engineering Study for the Hennepin County Light Rail Transit System. The 1988 study received Regional Transit Board funding, and survey results could be applied to a mode split model for projecting ridership on the proposed Light Rail Transit System. However, this survey was designed with the 1990 Travel Behavior Inventory in mind. The intention had been to update the 1988 survey in 1990 to be compatible with the 1990 Travel Behavior Inventory data. The results of the 1990 survey were used primarily to create a table of observed transit trips between each of the 1,200 traffic analysis zones in the region. This trip table was used to calibrate a new mode split model, which estimated future year travel by mode. The 1990 update survey focused on new routes and routes that had changed significantly since 1988. To preserve compatibility with the 1988 survey, the same survey questionnaire card was used, together with the same survey procedures for data collection. The procedures randomly sampled bus patrons during the transit trip, asking key questions about the patron and the transit trip. The survey card was intended for patrons to fill out quickly so it could be completed during the transit trip. The questions focused on conditions that have proven over time to significantly influence ridership. In all, a total of 20,126 valid survey records were processed. Adding surrogate trips, the survey data file is composed of a total of 27,159 trip records . About 10% of the records were filled out by persons who had answered more than one questionnaire.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Extended Modeling of DOE Sealed Canisters with Updated Chemistry Models

Road-ready and final disposition packaging configurations for the advanced test reactor (ATR) fuel currently specifies storage within helium backfilled DOE sealed standard canisters. The aluminum cladding of the ATR fuel contains an oxyhydroxide layer of boehmite/bayerite that generates hydrogen when subjected to irradiation. Understanding the effect of this hydrogen buildup over time to important for long term storage considerations. Previous modeling efforts have built a coupled CFD-chemical model to simulate the temperature gas phase concentrations within the DOE sealed standard canisters. These models have coupled the temperature conditions to both the gas phase radiolysis chemistry and the surface chemistry associated with the oxyhydroxide layer. The previous iteration of the model utilized constant G-values for the hydrogen generation a 50-year period. This new iteration of the model utilizes new experimental data to update the hydrogen generation rate, as well as increase the simulated time to a 200-year period. Given the half-life assumed for the primary Cs-137 isotope responsible for gamma radiation in the ATR spent fuel, the 200-year period decreases the decay heat and dose rate of the spent fuel by a factor of 100, which combined with the updated chemistry substantially decreases new hydrogen generation. Continued experimental work has identified trends for aluminum surrogate samples with oxyhydroxide layers for tests done at higher dose rates. At low initial doses a fast generation rate of hydrogen occurs which starts to roll over to a lower generation as the total dose applied increase. A small-scale chemical model was built to replicate as mini-canister surrogate system at SRNL as well as for the smaller capsule tests performed at INL. A variety of chemical models to capture this effect were tested, and a back reaction of H radical absorbing onto the surface, or inhibition of the reaction by significant H 2 cover gas were both able to fit both the mini-canister and the small capsule test data. Both kinetic fits were used to generate data from the new 200-year simulation. As additional experiments continue, the kinetic fits may be adjusted to adapt to new data. However, updated data shows that hydrogen atmosphere has little effect on actual generation data, so the model was reverted to use a star-stepped G-value for low-dose and high-dose regions. For the undried fuel case, the three models differ significantly with an end concentration of 1% for back reaction, 4.2% for inhibition, and 13.1% for constant G-value for the nominal scenario. For the dried fuel case, the nominal cases showed end concentrations of 0.23% for back reaction, 1.9% for inhibition, and 4.2% for const G-value for the nominal scenario. For the constant G-value case that is less conservative than the other two, the pressure for undried fuel increases to 1.94 atm over 200 years for the nominal case and 2.11 atm for the high decay heat case. The total hydrogen concentration after 200 years is 4.92% for the low decay heat case and 27.4% for the high decay heat case for undried fuel and is 1.4% and 9.6% for low and high decay heat fuel for the dried fuel case. If small amount of residual air is present, the potential for nitric acid formation of 595, 1568, and 2816 ppm for the lower, nominal, and upper fuel decay heat can occur. At long-term timeframes no significant shift in major species present occurs, so no appreciable amount of oxygen is present in the system. The continued rate of hydrogen generation in the canister occurs at an increase of 0.04% by mole over the final 10 years. Increasing the model range out to 1000 years continues to drop the hydrogen generation rate through decreasing dose rate.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Predicting nucleon-nucleus scattering observables using nuclear structure theory

Developing a predictive capability for inelastic scattering will find applications in multiple areas. Experimental data for neutron-nucleus inelastic scattering is limited and thus one needs a robust theoretical framework to complement it. Charged-particle inelastic scattering can be used as a surrogate for ( n, γ ) reactions to predict capture cross sections for unstable nuclei. Our work uses microscopic nuclear structure calculations for spherical nuclei to obtain nucleon-nucleus scattering potentials and calculate cross sections for these processes. We implement the Jeukenne, Lejeune, Mahaux (JLM) semi-microscopic folding approach, where the medium effects on nuclear interaction are parameterized in nuclear matter to obtain the nucleon-nucleon ( NN ) interaction in a medium at positive energies. We solve for the nuclear ground state using the Hartree-Fock-Bogliubov (HFB) many-body method, assuming the nucleons within the nucleus interact via the Gogny-D1M potential. The vibrational excited states of the target nucleus are calculated using the quasi-particle random phase approximation (QRPA). We demonstrate our approach for spherical nuclei in the medium-mass region, showing scattering results for the 90 Zr nucleus.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗