SEARCH · Engineering Papers

Results for “large deviations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Towards testing the theory of gravity with DESI: summary statistics, model predictions and future simulation requirements

Shortly after its discovery, General Relativity (GR) was applied to predict the behavior of our Universe on the largest scales, and later became the foundation of modern cosmology. Its validity has been verified on a range of scales and environments from the Solar system to merging black holes. However, experimental confirmations of GR on cosmological scales have so far lacked the accuracy one would hope for — its applications on those scales being largely based on extrapolation and its validity there sometimes questioned in the shadow of the discovery of the unexpected cosmic acceleration. Future astronomical instruments surveying the distribution and evolution of galaxies over substantial portions of the observable Universe, such as the Dark Energy Spectroscopic Instrument (DESI), will be able to measure the fingerprints of gravity and their statistical power will allow strong constraints on alternatives to GR. In this paper, based on a set of N-body simulations and mock galaxy catalogs, we study the predictions of a number of traditional and novel summary statistics beyond linear redshift distortions in two well-studied modified gravity models — chameleon f(R) gravity and a braneworld model — and the potential of testing these deviations from GR using DESI. These summary statistics employ a wide array of statistical properties of the galaxy and the underlying dark matter field, including two-point and higher-order statistics, environmental dependence, redshift space distortions and weak lensing. We find that they hold promising power for testing GR to unprecedented precision. The major future challenge is to make realistic, simulation-based mock galaxy catalogs for both GR and alternative models to fully exploit the statistic power of the DESI survey (by matching the volumes and galaxy number densities of the mocks to those in the real survey) and to better understand the impact of key systematic effects. Using these, we identify future simulation and analysis needs for gravity tests using DESI.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Non-Fermi liquid behaviour in a correlated flat-band pyrochlore lattice

Electronic correlation effects are manifested in quantum materials when either the on-site Coulomb repulsion is large or the electron kinetic energy is small. The former is the dominant effect in cuprate superconductors and heavy-fermion systems whereas it is the latter in twisted bilayer graphene and geometrically frustrated metals. However, the simultaneous cooperation of both effects in the same quantum material remains rare. The design aim is to produce correlated topological flat bands pinned at the Fermi level. Here, we observe a flat band at the Fermi level in a 3d pyrochlore metal CuV 2 S 4 . Our angle-resolved photoemission spectroscopy data reveal that destructive quantum interference associated with the V pyrochlore sublattice and further renormalization to the Fermi level by electron interactions induce this flat band. Consequently, we discover transport signatures that evidence a deviation from Fermi liquid behaviour as well as an enhanced Sommerfeld coefficient. Our work illustrates the combined cooperation of local Coulomb interactions and geometric frustration in a pyrochlore lattice system to induce correlated topology by constructing and pinning correlated flat bands near the Fermi level.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Redshift space power spectrum beyond Einstein-de Sitter kernels

Abstract We develop a framework to compute the redshift space power spectrum (PS), with kernels beyond Einstein-de Sitter (EdS), that can be applied to a wide variety of generalized cosmologies. We build upon a formalism that was recently employed for standard cosmology in Chen, Vlah & White (2020), and utilize an expansion of the density-weighted velocity moment generating function that explicitly separates the magnitude of the k-modes and their angle to the line-of-sight direction dependencies. We compute the PS for matter and biased tracers to 1-loop Perturbation Theory (PT) and show that the expansion has a correct infrared and ultraviolet behavior, free of unwanted divergences. We also add Effective Field Theory (EFT) counterterms, necessary to account for small-scale contributions to PT, and employ an IR-resummation prescription to properly model the smearing of the BAO due to large scale bulk flows within Standard-PT. To demonstrate the applicability of our formalism, we apply it on the ΛCDM and the Hu-Sawicki f(R) models, and compare our numerical results against the elephant suite of N-body simulations, finding very good agreement up to k = 0.27 Mpc -1 h at z = 0.5 for the first three non-vanishing Legendre multipoles of the PS. To our knowledge, the model presented in this work is the most accurate theoretical EFT-PT for modified gravity to date, being the only one that accounts for beyond linear local biasing in redshift-space. Hence, we argue our RSD modeling is a promising tool to construct theoretical templates in order to test deviations from ΛCDM using real data obtained from the next stage of cosmological surveys such as DESI and LSST.

Astronomy & Astrophysics↗

DOE OSTI2025

New precious metal containing normal spinels: LiRhRu 1-x Ir x O 4 , LiFeIr 1-x Ru x O 4 , and LiCoIr 1-x Ru x O 4

Large spin-orbit-coupled cations in geometrically frustrated crystal structures have the most suitable setting for exploring novel exotic states of matter. Spinel oxides (AM 2 O 4 ) are well-known examples of geometrically frustrated systems. In this study, we report for the first time the synthesis of compositions LiRhRu 1-x Ir x O 4 (x = 0–0.5), LiFeIr 1-x Ru x O 4 (x = 0–0.5), and LiCoIr 1-x Ru x O 4 (x = 0–0.3) containing precious metal cations on edge-sharing octahedral M-sites, and systematically investigate their magnetic and electrical properties. 57 Fe Mössbauer spectroscopy revealed that iron is trivalent in all LiFeIr 1-x Ru x O 4 solid solutions. Magnetic measurements indicate deviations from theoretical spin-only magnetic moment values, indicating the influence of spin-orbit coupling owing to the presence of 4d and 5d block elements. The LiFeIr 1-x Ru x O 4 series shows spin-glass-like freezing behavior with T g ≈ 20 K, and a small frustration index (f ≈ 1-2), indicating that the frustration originates from site disorder. LiRhRu 1-x Ir x O 4 and LiCoIr 1-x Ru x O 4 exhibit strongly geometrically frustrated magnetism. Electrical resistivity measurements as a function of temperature indicate that all phases are semiconducting. Seebeck coefficient measurements show that LiRhRu 1-x Ir x O 4 and LiFeIr 1-x Ru x O 4 are p-type semiconductors with holes as the major charge carriers. A sign reversal of the Seebeck coefficient indicates both holes and electrons as carriers for LiCoIr 1-x RuxO 4 (x = 0–0.2), but only holes as major carriers for x = 0.3. Here, the Seebeck coefficient and power factor increase drastically in the LiRhRu 1-x Ir x O 4 solid solution with Ir substitution, reaching a maximum of ≈ +125 μV/K and ≈2.3×10 -6 W/mK 2 at ∼650 K for x = 0.5.

Charge carriers↗

DOE OSTI2024

An idealized model for the crystal structure of intermetallic compounds isostructural with Mg 3 Cr 2 Al 18

Mg 3 Cr 2 Al 18 (abbreviated in this report as MCA) is the parent phase for a large class of intermetallic compounds that belong to the cubic crystal space group, $Fd\overline{3}m$. The purpose of this paper is to introduce an ideal, unrelaxed crystal structure for compounds isostructural with MCA. There are five distinct atomic sublattices in MCA compounds, which can be denoted, $A, B, C, D,$ and $E$. With this, a general description for MCA structures can be written as $A^{8a}_{1}B^{16c}_{2}C^{16d}_{2}D^{48f}_{6}E^{96g}_{12}$, where the superscripts represent the Wyckoff special equipoints associated with the various sublattices in MCA, and the subscripts indicate the contributions of each sublattice to the stoichiometry of one formula unit in an any given MCA structured compound. Sublattices D and E are where deviations from ideality occur in real, MCA-like compounds. This paper examines MCA bond lengths, nearest-neighbour polyhedral arrangements, 3-D sublattice crystal structures, 2-D atom tessellation patterns, and crystal chemical effects associated with atomic relaxations on the $D$ and $E$ sublattices. The ideal MCA crystal structure developed in this report provides an appropriate initial structure for use as input to crystal structure refinements of diffraction data for MCA-like phases being examined experimentally, or as input for computational, atomistic simulations of the structures of such compounds.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Probing the γ -Ray Emission Origin of Two Star-forming Galaxies NGC 2403 and NGC 3424 with the Fermi-LAT

Star-forming galaxies (SFGs) are a subclass of γ-ray emitters, and a correlation between their γ-ray luminosity (L γ ) and the total infrared (IR) luminosity (L IR ) has been established based on the Fermi Large Area Telescope (LAT) data. NGC 2403 and NGC 3424 have been reported as outliers in the L γ –L IR correlation with light curves showing significant variability, which contrasts with the temporally stable γ-ray emission in other SFGs, originating primarily from cosmic rays interacting with interstellar medium. In this study, we reanalyze the γ-ray emission in the directions of NGC 2403 and NGC 3424 using more than 16.5 yr Fermi-LAT data. NGC 3424 is found to be spatially coincident with the detected γ-ray source, while NGC 2403 is significantly offset from the nearest γ-ray source, suggesting an implausible association. We confirm the previously reported variability of both γ-ray sources and the significant deviation from the L γ –L IR correlation when assuming an association of both γ-ray sources with the two galaxies. Our findings lend further support to the interpretation that their γ-ray emission is driven primarily by alternative radiative processes—rather than by star formation activity—such as the ejecta of the Type IIP supernova SN 2004dj in NGC 2403 interacting with a surrounding high-density shell and an obscured active galactic nucleus in NGC 3424.

Liu, Linjie [Chinese Academy of Sciences (CAS), Ku↗

DOE OSTI2026

Physics-aware adaptive checkpointing with shadow systems for nonlinear PDE simulations

Large-scale simulations of nonlinear partial differential equations (PDEs) that exhibit strongly transient behavior and pattern-forming dynamics produce enormous amounts of data, which, even with modern storage systems, cannot be stored for later curation. Current I/O strategies either write dense time series of snapshots, which is often prohibitive in I/O and storage, or store a few checkpoints that enable restart but incur expensive recomputation cost and provide no control over post-restart error growth, especially when lossy compression is used. Moreover, most, if not all, existing strategies take no account of the actual physical state of the system. Here, we present a simple physics-aware I/O framework in which a low-cost shadow system adaptively triggers lossy checkpoints when the shadow system deviates from the fine-scale simulation. The shadow system can be a coarsened replica of the fine-scale simulation that evolves concurrently. This means that checkpoints are taken based on the physical state of the system: fewer checkpoints are triggered when the system is quiescent while more are taken when the system undergoes a rapid change. This type of behavior is observed in many systems such as Brusselator and FitzHugh–Nagumo. We illustrate that our framework maintains stable restarts, keeps fine-scale restart errors bounded by shadow errors, and reconstructs the time history with significantly lower error and storage than interpolating fixed-interval snapshots, with low-cost shadow replay and modest online synchronization overhead.

Gong, Qian [ORNL] (ORCID:0000000235704142)↗

DOE OSTI2020

Reaction Profiles and Kinetics for Radical–Radical Hydrogen Abstraction via Multireference Coupled Cluster Theory

Radical-radical abstractions in hydrocarbon oxidation chemistry are disproportionation reactions that are generally exothermic with little or no barrier yet are underappreciated and poorly studied. Such challenging multireference electronic structure problems are tackled here using the recently developed state-specific multireference coupled cluster methods Mk-MRCCSD and Mk-MRCCSD(T), as well as the companion perturbation theory Mk-MRPT2 and the established MRCISD, MRCISD+Q, and CASPT2 approaches. Reaction paths are investigated for five prototypes involving radical-radical hydrogen abstraction: H + BeH → H 2 + Be, H + NH 2 → H 2 + NH, CH 3 + C 2 H 5 → CH 4 + C 2 H 4 , H + C2H 5 → H 2 + C 2 H 4 , and H + HCO → H 2 + CO. Full configuration interaction (FCI) benchmark computations for the H + BeH, H + NH 2 , and H + HCO reactions prove that Mk-MRCCSD(T) provides superior accuracy for the interaction energies in the entrance channel, with mean absolute errors less than 0.3 kcal mol(-1) and percentage deviations less than 10% over the fragment separations of relevance to kinetics. To facilitate combustion studies, energetics for the CH 3 + C 2 H 5 , H + C 2 H 5 , and H + HCO reactions were computed at each level of theory with correlation-consistent basis sets (cc-pVXZ, X = T, Q 5) and extrapolated to the complete basis set (CBS) limit. These CBS energies were coupled with CASPT2 projected vibrational frequencies along a minimum energy path to obtain rate constants for these three reactions. The rigorous Mk-MRCCSD(T)/CBS results demonstrate unequivocally that these three reactions proceed with no barrier in the entrance channel, contrary to some earlier predictions. Mk-MRCCSD(T) also reveals that the economical CASPT2 method performs well for large interfragment separations but may deteriorate substantially at shorter distances.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Insights into Rotational and Translational Dynamics in Mixtures of Ethylene Glycol and Choline Chloride Using Nuclear Magnetic Resonance Techniques

This work examines molecular dynamics and interactions in ethylene glycol–choline chloride (EG–ChCl) mixtures across 0–33 mol % ChCl, spanning the true eutectic region near 17–20 mol % and the commonly used 1:2 formulation. We combine pulsed-field-gradient (PFG) diffusion, fast-field-cycling (FFC) relaxometry, temperature-dependent 13 C T 1 , and nuclear Overhauser effect spectroscopy (NOESY) to disentangle local from macroscopic dynamics. PFG and FFC show that both translational and average rotational motions largely track the strong increase in viscosity with ChCl content, with ethylene glycol consistently diffusing faster than the choline cation and no global dynamical anomaly at the eutectic composition. More subtle, site-specific composition effects nevertheless emerge. The ratio of the diffusion coefficient of the hydroxyl group of choline to the diffusion coefficient of the methyl group of choline displays a shallow minimum in the 17–25 mol % region, indicating a modest change in how the hydroxyl-bearing end of choline samples the underlying translational motion relative to the methyl groups. 13 C T 1 analysis shows that rotational correlation times at 25 °C generally increase with ChCl, reflecting viscosity-coupled slowing, while the CH 2 –N α site exhibits a small but reproducible deviation from this monotonic trend near the eutectic. NOESY spectra at similar compositions reveal enhanced cross-relaxation between EG and choline protons, consistent with increased headgroup–solvent contact density rather than a wholesale structural rearrangement. Overall, our multitechnique study demonstrates that EG–ChCl dynamics are predominantly viscosity-dominated, with the eutectic region acting as a subtle dynamical crossover where specific choline segments become maximally coupled to the hydrogen-bond network. These insights refine the structure–dynamics picture of choline-chloride DESs and provide practical guidance for tuning composition in electrochemical, separation, and catalytic applications.

diffusion↗

DOE OSTI2021

Microbial Community Field Surveys Reveal Abundant Pseudomonas Population in Sorghum Rhizosphere Composed of Many Closely Related Phylotypes

While the root-associated microbiome is typically less diverse than the surrounding soil due to both plant selection and microbial competition for plant derived resources, it typically retains considerable complexity, harboring many hundreds of distinct bacterial species. Here, we report a time-dependent deviation from this trend in the rhizospheres of field grown sorghum. In this study, 16S rRNA amplicon sequencing was used to determine the impact of nitrogen fertilization on the development of the root-associated microbiomes of 10 sorghum genotypes grown in eastern Nebraska. We observed that early rhizosphere samples exhibit a significant reduction in overall diversity due to a high abundance of the bacterial genus Pseudomonas that occurred independent of host genotype in both high and low nitrogen fields and was not observed in the surrounding soil or associated root endosphere samples. When clustered at 97% identity, nearly all the Pseudomonas reads in this dataset were assigned to a single operational taxonomic unit (OTU); however, exact sequence variant (ESV)-level resolution demonstrated that this population comprised a large number of distinct Pseudomonas lineages. Furthermore, single-molecule long-read sequencing enabled high-resolution taxonomic profiling revealing further heterogeneity in the Pseudomonas lineages that was further confirmed using shotgun metagenomic sequencing. Finally, field soil enriched with specific carbon compounds recapitulated the increase in Pseudomonas, suggesting a possible connection between the enrichment of these Pseudomonas species and a plant-driven exudate profile.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Extremely Metal-poor Galaxies in DESI DR1: Connections to Galaxies in the Early Universe

Extremely metal-poor galaxies (XMPGs), defined as having metallicities below 10% of the solar value, are considered possible local analogs to primordial systems and offer a unique window into early galaxy evolution. This study presents a large-scale search for XMPGs using data from the Dark Energy Spectroscopic Instrument Data Release 1, systematically evaluating their resemblance to high-redshift galaxies. From a parent sample of more than 14 million galaxies, we identify 662 (556 new) confirmed XMPGs and 763 (666 new) high-quality candidates via the direct T$_{e}$ method. Our results reveal that XMPGs follow a distinct star-forming main sequence (SFMS) that is elevated and shallower than that of the comparing star-forming galaxies. Notably, at higher stellar masses (M$_{⋆}$ > 10$^{7.5}$M$_{⊙}$), the XMPG SFMS converges with the sequence observed in high-redshift galaxies by the James Webb Space Telescope, indicating that mature XMPGs sustain star formation rates comparable to their primordial counterparts. Furthermore, XMPGs consistently deviate below the local fundamental metallicity relation, mirroring high-redshift galaxy behavior. These findings demonstrate that XMPGs not only exhibit low metallicities but also preserve scaling relations characteristic of the early Universe, confirming their potential value as local laboratories for studying early galaxy formation processes.

Sui, Jipeng [Beijing Observ.; Beijing, GUCAS] (ORC↗

DOE OSTI2025

Commercial building HVAC demand flexibility with model predictive control: Field demonstration and literature insights

Model Predictive Control (MPC) for building Heating Ventilation and Air Conditioning (HVAC) systems is beginning to gain traction in the market, with a few controls companies incorporating it into their product offerings. However, it remains difficult to assess whether the energy cost savings are enough to justify the cost of MPC implementation for a particular building, given the limited number of reported demonstrations. For small commercial and residential buildings with relatively uniform systems, standardized approaches can help lower implementation costs. In contrast, for large buildings or district systems, the potential magnitude of cost savings could justify more customized solutions. Estimating the cost-effectiveness of MPC becomes more challenging for medium and large commercial buildings, where a one-size-fits-all solution may not be suitable, and the potential energy cost savings may be insufficient to justify a customized solution. To make MPC technology more appealing, incorporating additional value streams beyond energy efficiency alone can significantly increase its attractiveness. One such revenue stream is demand flexibility, in response to dynamic electricity prices, where MPC can leverage the thermal mass of the building to shift the load and support the grid. Building on an extensive literature review of MPC field studies focused on cost savings and demand flexibility, this paper presents the results of implementing MPC control in a large office building HVAC system in Berkeley, CA. Four different dynamic electricity price profiles were integrated into the MPC objective function to shift building demand while maintaining comfort, and field testing was performed with each price profile across four seasons. The results show potential for 40–65 % demand decrease percentage and up to 61 % annual cost savings compared to the existing rule-based control strategy, under the tested dynamic price scenarios. This paper also presents a sensitivity analysis on the cost savings with respect to the price profile variability, discusses the implementation effort for the price-responsive MPC, and compares the cost savings found in this study to those found in literature on the basis of dynamic price variability, or so-called Electricity Price Relative Standard Deviation.

Zanetti, Ettore↗

DOE OSTI2023

Emulation of seismic-phase traveltimes with machine learning

SUMMARY We present a machine learning (ML) method for emulating seismic-phase traveltimes that are computed using a global-scale 3-D earth model and physics-based ray tracing. Accurate traveltime predictions based on 3-D earth models are known to reduce the bias of event location estimates, increase our ability to assign phase labels to seismic detections and associate detections to events. However, practical use of 3-D models is challenged by slow computational speed and the unwieldiness of pre-computed lookup tables that are often large and have prescribed computational grids. In this work, we train a ML emulator using pre-computed traveltimes, resulting in a compact and computationally fast way to approximate traveltimes that are based on a 3-D earth model. Our model is trained using approximately 850 million P-wave traveltimes that are based on the global LLNL-G3D-JPS model, which was developed for more accurate event location. The training-set consists of traveltimes between 10 393 global seismic stations and randomly sampled event locations that provide a prescribed, distance-dependent geographic sample density for each station. Prediction accuracy is dependent on event-station distance and whether the station was included in the training set. For stations included in the training set the mean absolute deviation (MAD) of the difference between traveltimes computed using ray tracing through the 3-D model and the ML emulator for local, regional, and teleseismic distances are 0.090, 0.125 and 0.121 s, respectively. For tested station locations not included in the training set, MAD values for the three distance ranges increase to 0.173, 0.219 and 0.210 s, respectively. Empirical traveltime residuals for a global reference data are indistinguishable when ML emulation or the 3-D model is used to compute traveltimes. This result holds regardless of whether the recording station is used in ML training or not.

58 GEOSCIENCES↗

DOE OSTI2025

Assessing Heterogeneity of Surface Water Temperature Following Stream Restoration and a High-Intensity Fire from Thermal Imagery

Thermal heterogeneity of rivers is essential to support freshwater biodiversity. Salmon behaviorally thermoregulate by moving from patches of warm water to cold water. When implementing river restoration projects, it is essential to monitor changes in temperature and thermal heterogeneity through time to assess the impacts to a river’s thermal regime. Lightweight sensors that record both thermal infrared (TIR) and multispectral data carried via unoccupied aircraft systems (UASs) present an opportunity to monitor temperature variations at high spatial (<0.5 m) and temporal resolution, facilitating the detection of the small patches of varying temperatures salmon require. Here, we present methods to classify and filter visible wetted area, including a novel procedure to measure canopy cover, and extract and correct radiant surface water temperature to evaluate changes in the variability of stream temperature pre- and post-restoration followed by a high-intensity fire in a section of the river corridor of the South Fork McKenzie River, Oregon. We used a simple linear model to correct the TIR data by imaging a water bath where the temperature increased from 9.5 to 33.4 °C. The resulting model reduced the mean absolute error from 1.62 to 0.35 °C. We applied this correction to TIR-measured temperatures of wetted cells classified using NDWI imagery acquired in the field. We found warmer conditions (+2.6 °C) after restoration (p < 0.001) and median absolute deviation for pre-restoration (0.30) to be less than both that of post-restoration (0.85) and post-fire (0.79) orthomosaics. In addition, there was statistically significant evidence to support the hypothesis of shifts in temperature distributions pre- and post-restoration (KS test 2009 vs. 2019, p < 0.001, D = 0.99; KS test 2019 vs. 2021, p < 0.001, D = 0.10). Moreover, we used a Generalized Additive Model (GAM) that included spatial and environmental predictors (i.e., canopy cover calculated from multispectral NDVI and photogrammetrically derived digital elevation model) to model TIR temperature from a transect along the main river channel. This model explained 89% of the deviance, and the predictor variables showed statistical significance. Collectively, our study underscored the potential of a multispectral/TIR sensor to assess thermal heterogeneity in large and complex river systems.

Barker, Matthew I. (ORCID:0000000252864930)↗

DOE OSTI2022

Accumulator Sizing for the Hybrid Hydraulic Electric Architecture (HHEA) Using Dynamic Programming

Abstract The Hybrid Hydraulic Electric Architecture (HHEA) is an approach to electrifying and improving that does not require very large electric machines. The majority power is provided hydraulically via a set of common pressure rails while small electric machines are used to modulate that power to meet exact demand. Each common pressure rail has an accumulator, which keeps the pressure of each rail near constant as long as the volume in the accumulator does not change much. A main pump adds or removes fluid from each accumulator as needed to reduce volume changes in the accumulators. It is the operation of this main pump which is investigated here. In previous analysis it was assumed that the size of the accumulators was large enough to account for any difference in fluid volume without significantly changing the pressure. In this paper, the required accumulator sizes are studied using a dynamic programming approach. At each time step the inlet and outlet of the main pump could be connected to any of the common pressure rails. These decisions for the main pump were optimized to minimize the required accumulator size. Case studies are conducted using drive cycles from a 22-ton excavator, a 20-ton wheel loader and a 5-ton excavator. The machines could be operated with total accumulator sizes of less than 20, 10, and 5 liters respectively. Required accumulator size was found to vary with pump flow rate, the frequency at which the pump/motor was able to change operating conditions, the tolerance level of pressure deviation and the accumulator’s polytropic index.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Transport in Stochastic Media with Random Chord Length Distributions

Thermal radiation transport computations in binary Markovian random mixtures rely almost exclusively on the Levermore-Pomraning (LP) model which is obtained by applying a heuristic closure to the ensemble averaged random medium transport equation. The validity of this model has been extensively tested by comparing numerical results over a broad parameter range (material types and mixing parameters) against benchmark solutions in planar geometry. The conditions under which the LP-model provides useful results and when it breaks down are now well established, but work to date has been largely restricted to homogeneous mixing statistics, i.e., the mean chord lengths of both materials are taken to be spatially constant. In recent work, this limitation was relaxed by allowing the mean chord lengths and, in a consistent fashion, the volume fractions in the LP-model to vary spatially and in a follow-up investigation benchmark solutions were obtained by ensemble averaging results over material realizations sampled from a nonhomogeneous Poisson process (NHPP). Numerical experiments in rod geometry with specifically linear and quadratic spatial dependence of chord lengths showed that the material averaged radiation intensities vary nonmonotonically with depth into the medium, in stark contrast to solutions obtained assuming uniform chord lengths. Moreover, depending on the local optical thickness and strength of scattering, the LP-model results showed locally more nuanced deviations from the benchmark solutions than was the case with constant chord lengths. These limited numerical investigations highlight the nontrivial qualitative and quantitative consequences of nonhomogeneous mixing statistics, in particular that closure approximations may not be uniformly valid or invalid over the problem domain.

42 ENGINEERING↗

DOE OSTI2020

Application of entropy and signal energy for ultrasound-based classification of three-dimensional printed polyetherketoneketone components

This paper describes a preliminary method for the classification of annealed and unannealed polyetherketoneketone (PEKK) components manufactured using a material extrusion three-dimensional (3D) printing process. PEKK is representative of a class of high-performance thermoplastics that are increasingly employed as feedstocks for use in 3D printing. PEKK components may be used continuously at elevated temperatures, are chemically resistant, and able to withstand large mechanical loads. These properties render PEKK suitable as a metal component replacement in aerospace applications, high-temperature industrial applications, and surgical implants. The structure of PEKK is semi-crystalline with the specific crystallinity correlating to the final properties during application, making determination of this property crucial. This study compares three different signal processing techniques intended to distinguish annealed (high crystallinity) from unannealed (low crystallinity) components using backscattered ultrasound. The first is energy-based and is unable to detect annealing. The second two are based on different entropies of the backscattered signal: a limiting form of Renyi's entropy and a limiting form of joint entropy. The joint entropy values for the annealed and unannealed specimens fall into two non-overlapping intervals and have a statistical separation of two standard deviations.

36 MATERIALS SCIENCE↗

DOE OSTI2026

3D Deep Learning Joint Inversion of Active Seismic Full Waveform and Passive Seismic Traveltime Data for Reservoir Imaging and Uncertainty Quantification

Here, we present deep learning (DL) networks for three-dimensional (3D) joint inversion of active seismic full waveform and passive seismic traveltime data to image reservoirs and their properties and quantify imaging uncertainties. Active seismic full-waveform data can provide high-resolution monitoring images but are collected only intermittently because of their high acquisition cost. In contrast, passive seismic data can be gathered at relatively low cost between regular active surveys, although their imaging quality can be compromised by factors such as low signal-to-noise ratios and limited ray coverage of the target. Although these datasets are routinely acquired together at CO 2 storage sites, their combined inversion within a 3D DL framework has not been previously demonstrated. To our knowledge, this is the first study to address this gap, combining the strength of both data types. For efficient data storage and DL training with large 3D seismic datasets, we use a 3D data matrix in which a random number of passive seismic traveltime data are stored as parabolic envelopes using one-hot encoding and a 3D full-waveform data matrix in which multiple shot gathers are summed. Two network architectures are evaluated: a single-encoder U-Net for single-data type inversion and a dual-encoder U-Net for joint inversion of active and passive seismic data. We also evaluate the single-encoder U-Net for joint inversion by concatenating full-waveform data and traveltime data. We propose a systematic approach for selecting an optimal dropout rate that balances regularization during training and Monte Carlo dropout-based uncertainty quantification during prediction by examining the correlation coefficient between standard deviation and prediction error, along with the training misfit, across a range of dropout rates. 3D DL inversion experiments include five different network configurations, with evaluations under ideal, noisy and dropout-enabled conditions. Both model and data uncertainties are assessed, as well as their combined effects. Across all conditions, the networks consistently predict accurate CO 2 saturation models with low prediction errors, such as a structural similarity index of 0.993 and CO 2 difference of 1.1%. Uncertainty estimates show strong spatial correlation with prediction errors, confirming the effectiveness of the proposed dropout selection approach. The results demonstrate that our DL approach, utilizing compact data representations and appropriate uncertainty quantification, yields accurate subsurface images under various inversion conditions and provides valuable insights into the reliability of predictions.

Um, Evan Schankee [Lawrence Berkeley National Labo↗