Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sparse”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Simulations of plasmas and fluids using anti-symmetric models

ALMA (Anti-symmetric, Large-Moment, Accelerated) is a fast, flexible, and scalable toolkit designed to solve hyperbolic conservation law systems in hybrid supercomputers. Here this manuscript describes the theoretical background and implementation of ALMA, which uses the anti-symmetric formulation of fluids to obtain simple, robust, and easily paralellizable code. Practical GPU acceleration is realized on entire applications with an overall gain factor of 2 to 4. ALMA also provides a parallel, GPU accelerated sparse solver based on geometric multigrid, capable of diagonalizing linear systems with 239 unknowns. Here we demonstrate ALMA's scaling and performance in petascale supercomputers and use standard fluid models to verify the overall approach with canonical benchmark problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Permutation-adapted complete and independent basis for atomic cluster expansion descriptors

Atomic cluster expansion (ACE) methods provide a systematic way to describe particle local environments of arbitrary body order. For practical applications it is often required that the basis of cluster functions be symmetrized with respect to rotations and permutations. Existing methodologies yield sets of symmetrized functions that are over-complete. These methodologies thus require an additional numerical procedure, such as singular value decomposition (SVD), to eliminate redundant functions. In this work, it is shown that analytical linear relationships for subsets of cluster functions may be derived using recursion and permutation properties of generalized Wigner symbols. From these relationships, subsets (blocks) of cluster functions can be selected such that, within each block, functions are guaranteed to be linearly independent. It is conjectured that this block-wise independent set of permutation-adapted rotation and permutation invariant (PA-RPI) functions forms a complete, independent basis for ACE. Along with the first analytical proofs of block-wise linear dependence of ACE cluster functions and other theoretical arguments, numerical results are offered to demonstrate this. The utility of the method is demonstrated in the development of an ACE interatomic potential for tantalum. Using the new basis functions in combination with Bayesian compressive sensing sparse regression, some high degree descriptors are observed to persist and help achieve high-accuracy models.

Angular momentum↗
DOE OSTI2025

A dynamic likelihood approach to filtering transport processes: advection-diffusion dynamics

A Bayesian data assimilation scheme is formulated for advection-dominated advective and diffusive evolutionary problems, based upon the Dynamic Likelihood (DLF) approach to filtering. The DLF was developed specifically for hyperbolic problems –waves–, and in this paper, it is extended via a split step formulation, to handle advection-diffusion problems. In the dynamic likelihood approach, observations and their statistics are used to propagate probabilities along characteristics, evolving the likelihood in time. The estimate posterior thus inherits phase information. For advection-diffusion the advective part of the time evolution is handled on the basis of observations alone, while the diffusive part is informed through the model as well as observations. We expect, and indeed show here, that in advection-dominated problems, the DLF approach produces better estimates than other assimilation approaches, particularly when the observations are sparse and have low uncertainty. The added computational expense of the method is cubic in the total number of observations over time, which is on the same order of magnitude as a standard Kalman filter and can be mitigated by bounding the number of forward propagated observations, discarding the least informative data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Joint state-parameter estimation for the reduced fracture model via the united filter

Here, in this paper, we introduce an effective United Filter method for jointly estimating the solution state and physical parameters in flow and transport problems within fractured porous media. Fluid flow and transport in fractured porous media are critical in subsurface hydrology, geophysics, and reservoir geomechanics. Reduced fracture models, which represent fractures as lower-dimensional interfaces, enable efficient multi-scale simulations. However, reduced fracture models also face accuracy challenges due to modeling errors and uncertainties in physical parameters such as permeability and fracture geometry. To address these challenges, we propose a United Filter method, which integrates the Ensemble Score Filter (EnSF) for state estimation with the Direct Filter for parameter estimation. EnSF, based on a score-based diffusion model framework, produces ensemble representations of the state distribution without deep learning. Meanwhile, the Direct Filter, a recursive Bayesian inference method, estimates parameters directly from state observations. The United Filter combines these methods iteratively: EnSF estimates are used to refine parameter values, which are then fed back to improve state estimation. Numerical experiments demonstrate that the United Filter method surpasses the state-of-the-art Augmented Ensemble Kalman Filter, delivering more accurate state and parameter estimation for reduced fracture models. This framework also provides a robust and efficient solution for PDE-constrained inverse problems with uncertainties and sparse observations.

Bayesian inference↗
DOE OSTI2026

Active operator learning with predictive uncertainty quantification for partial differential equations

With the increased prevalence of neural operators being used to provide rapid solutions to partial differential equations (PDEs), understanding the accuracy of model predictions and the associated error levels is necessary for deploying reliable surrogate models in scientific applications. Existing uncertainty quantification (UQ) frameworks employ ensembles or Bayesian methods, which can incur substantial computational costs during both training and inference. Here, we propose a lightweight predictive UQ method tailored for Deep operator networks (DeepONets) that also generalizes to other operator networks. Numerical experiments on linear and nonlinear PDEs demonstrate that the framework’s uncertainty estimates are unbiased and provide accurate out-of-distribution uncertainty predictions with a sufficiently large training dataset. Our framework provides fast inference and uncertainty estimates that can efficiently drive outer-loop analyses that would be prohibitively expensive with conventional solvers. We demonstrate how predictive uncertainties can be used in the context of Bayesian optimization and active learning problems to yield improvements in accuracy and data-efficiency for outer-loop optimization procedures. In the active learning setup, we extend the framework to Fourier Neural Operators (FNO) and describe a generalized method for other operator networks. To enable real-time deployment, we introduce an inference strategy based on precomputed trunk outputs and a sparse placement matrix, reducing evaluation time by more than a factor of five. Our method provides a practical route to uncertainty-aware operator learning in time-sensitive settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A score-based diffusion model approach for adaptive learning of stochastic partial differential equation solutions

In this paper, we propose a novel framework for adaptively learning the time-evolving solutions of stochastic partial differential equations (SPDEs) using score-based diffusion models within a recursive Bayesian inference setting. SPDEs play a central role in modeling complex physical systems under uncertainty, but their numerical solutions often suffer from model errors and reduced accuracy due to incomplete physical knowledge and environmental variability. To address these challenges, we encode the governing physics into the score function of a diffusion model using simulation data and incorporate observational information via a likelihood-based correction in a reverse-time stochastic differential equation. This enables adaptive learning through iterative refinement of the solution as new data becomes available. To improve computational efficiency in high-dimensional settings, we introduce the ensemble score filter, a training-free approximation of the score function designed for real-time inference. Numerical experiments on benchmark SPDEs demonstrate the accuracy and robustness of the proposed method under sparse and noisy observations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Influence of background sources and topographic resolution in the Weather Research and Forecasting Model on xenon plume characteristics at monitoring stations

For many atmospheric monitoring applications, networks of measurement sites—such as the radionuclide stations of the International Monitoring System—can be sparse. With measurement locations potentially hundreds to thousands of kilometers from a release it is important to quantify the effects of physical processes on transport and dispersion of plumes between source and measurement locations. This study addresses the effects of background sources and topography resolution near the release location of radionuclides. We use the Weather Research and Forecasting (WRF) model with inline chemistry to investigate (1) how an additional, time-varying source of 133 Xe, such as an operational medical isotope production facility, contributes to activity concentration measurements at monitoring sites, and (2) how complex topography influences on atmospheric conditions near emission sources impact plume concentrations at varying distances from the source. Two 133 Xe emission sources, including (1) a high flux rate of short duration representing an explosive event, and (2) a variable and continuous background source, are simulated. The continuous background source contributes significantly to total 133 Xe concentrations at several monitoring stations. Further, a WRF simulation at 9 km horizontal resolution is compared with a nested grid simulation, where the innermost domain has a resolution of 1 km. Increased topographic resolution leads to an improved representation of plume responses to local winds, with topographic influences greatest at locations closest to the sources. Differences between the two domain resolutions decrease at greater distances from the sources, as plumes have time to spread and mix and are influenced by synoptic scale circulation patterns that are represented similarly in both simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

High-resolution modeling of indoor radon exposure with uncertainty quantification in Utah

Indoor radon accounts for 37% of population-level exposure to ionizing radiation in the United States. However, radon metrics are typically reported at coarse spatial scales, potentially obscuring meaningful local variation. We developed a high-resolution modeling framework to estimate indoor radon concentrations across Utah while explicitly quantifying predictive uncertainty. A total of 19,497 residential radon measurements collected between 2006 and 2017 were combined with environmental and housing characteristics and analyzed using a geospatial neural network that accommodates spatial dependence and nonlinear associations. Predictions were generated on a uniform hexagonal grid at 0.73 km2 resolution (H3 level 8). Out-of-sample predictions aggregated to the H3 level 8 grid showed good agreement with observed concentrations (Pearson r=0.64), while household-level predictions exhibited more moderate agreement (r=0.45). The model produced well-calibrated uncertainty estimates, with 24.1% of held-out observations exceeding the predicted 75th-percentile threshold. Maps of predicted radon concentrations and the probability of exceeding the U.S. EPA action level of 148 Bq/m3 (4 pCi/L) revealed substantial fine-scale spatial heterogeneity that was not apparent in conventional coarse-resolution summaries, with greater local variability observed in densely monitored urban counties than in sparsely sampled regions. High-resolution radon models that explicitly quantify uncertainty provide a useful framework for characterizing the spatial distribution of indoor radon and identifying areas of elevated exceedance risk. These findings highlight the value of fine-scale monitoring data and uncertainty-aware modeling approaches for radon exposure assessment, environmental risk characterization, and radon-related health research.

Wu, Yunhan [ORNL] (ORCID:0000000178842994)↗
DOE OSTI2026

High pressure effects on EBC systems in high temperature environments

Environmental barrier coatings (EBCs) are designed to protect SiCfiber/SiCmatrix ceramic matrix composites (CMCs) in turbine engines by mitigating wear in high-temperature water vapor environments. The failure of EBCs is frequently attributed to the accelerated oxidation of the silicon bond coating layer when exposed to high-temperature steam, leading to the formation of a thickened thermally grown oxide (TGO). TGO growth increases interfacial stress, weakens adhesion, and results in coating spallation. Understanding the impact of high-temperature oxidation of EBC systems is essential for developing accurate lifespan models for turbine components, although pressurized oxidation testing is extremely sparse in the open literature. In this work, oxidation tests were performed on rare earth silicate EBCs coated onto SiC substrates under increased pressure conditions. The coated specimens were tested at 1100°C, 1200°C, and 1300°C at both 1 atm and 10 atm total pressure in steam environments. Subsequent characterization focused on the microstructural evolution of the EBC/Si/SiC system. The experimental findings indicated that TGO behavior is dependent on high-pressure conditions, with elevated pressure leading to an increase in oxide scale thickness and modifications in its morphology.

Ardrey, Kristyn [ORNL] (ORCID:0000000184407796)↗
DOE OSTI2024

Hydrologic applicability of satellite-based precipitation estimates for irrigation water management in the data-scarce region

Reliable precipitation estimates are crucial for planning and managing water resources, monitoring hydrologic extremes, and fulfilling irrigation water requirements. Accurate precipitation estimates are particularly challenging in complex mountain terrains, where monitoring gauges are often sparsely distributed due to their remote locations, and high installation and long-term operation costs. Recent advances in satellite-based precipitation estimates offer promising opportunities to improve our understanding of hydrologic processes and their applications for irrigation water management. Several datasets are available varying considerably in terms of their data sources, quality control methods, estimation procedure, and spatiotemporal resolutions. Choosing the most suitable dataset for a particular application is a complex task. In this study, we (1) evaluate the performance of six satellite-based precipitation estimates (SPEs): i) CHIRPS v2.0, ii) CMORPH v1.0, iii) ERA5, iv) IMERG v6, v) MSWEP v2.8, and vi) PERSIANN-CDR against the gauge precipitation using continuous statistical and categorical indices, (2) integrate SPEs with a calibrated semi-distributed hydrologic model to predict streamflow, and (3) demonstrate practical implications of improved streamflow prediction for irrigation water management in the central Himalayan region, Nepal. Our results illustrate that satellite-based precipitation estimates have competitive performance in capturing a wide range of rainfall characteristics, with demonstrated variability across river basins and time scales. Further, there are no significant discrepancies observed in satellite-based precipitation estimates for estimating irrigation water requirements for the three major crops (maize, wheat, and paddy) during the cropping period across the selected river basins, showing a greater promise for irrigation water management planning and decision making.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Inference of phase field fracture models

The phase field approach to modeling fracture uses a diffuse damage field to represent cracks. This representation mollifies singularities that arise in computations with sharp interface models and some of the resultant difficulties in the mathematical and numerical treatment of fracture. Phase field fracture models have proven effective at representing crack propagation, branching, and merging. Specific formulations, beginning with brittle fracture, have also been shown to converge to classical solutions. Extensions to cover the range of material failure, including ductile and cohesive fracture, lead to an array of possible models. There exists a large body of literature focusing on this class of models and on the impact of model form on the predicted crack evolution. However, there have not been systematic studies into how optimal models may be chosen. Here, we take a first step in this direction by developing formal methods for identification of the best parsimonious model of phase field fracture given full-field data on the damage and deformation fields. We consider some of the main models that have been used for the degradation of elastic response due to damage and its propagation. Our approach builds upon Variational System Identification (VSI), a weak form variant of the Sparse Identification of Nonlinear Dynamics (SINDy). Furthermore, in this first communication we focus on synthetically generated data but we also consider central issues associated with the use of experimental full-field data, such as data sparsity and noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Simple through-plane spatial filter for 2D MRI projections

A simple spatial filter for 2D projection MR imaging is introduced. It works in the third (unresolved) direction to eliminate uniform or slowly varying interfering background signals. A constant amplitude gradient pulse in the unresolved direction is applied at the same time as the usual phase encode gradient during 2D acquisition. The filter is demonstrated for root imaging in soil, where background soil water signals can be troublesome. The filter suppresses the soil water signal while preserving the desired signal of plant roots. Fundamental to the operation of the filter is that the roots are sparse in the image domain, meaning there are relatively few pixels with multiple roots present. Here, the performance of the through-plane filter is demonstrated and compares favorably to more conventional in-plane spatial filtering.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Role of interaction between hydraulic and natural fractures on production

One of the main purposes of hydraulic fracturing in unconventional resources such as shale is to improve the connectivity between existing natural fractures and the production well. Since the inherent permeability of shale is extremely low, this newly formed fracture connectivity provides faster paths for the hydrocarbon stored in preexisting natural fractures. In this study, we performed numerical simulations of the free gas production at large scale fractured reservoir with different pre-existing natural fracture network intensities and with variable hydraulic fracture sizes and staging. We explored production characteristics as a function of the properties of both, natural fractures (intensity) and hydraulic fractures (size and staging). Our hypothesis is that once maximum fracture connectivity is achieved, additional stimulation does not increase production significantly. In this numerical study we use an advanced modeling tool to study fractures connectivity, Discrete Fracture Network (DFN) approach, which can represent fracture networks similar to those observed on individual reservoir sites. In our model the horizontal production well and vertical hydraulic fractures were deterministically defined and surrounded by stochastically generated natural fracture networks. The performed numerical experiments show how fracture connectivity depends on hydraulic fracture settings in fractured reservoir with different preexisting fracture densities. We numerically showed that cumulative free gas production is not affected significantly by hydraulic fractures size (25 m–150 m) or spacing (15 m–45 m) for dense natural fracture networks. However, the hydraulic fracture size is important in reservoirs with sparse natural fractures, although spacing (number of stages) does not play a significant role.

58 GEOSCIENCES↗
DOE OSTI2020

The Sporadic History of Rubidium and its Role in Corrosion of Steel Related to Nuclear Material Storage

The influence of rubidium (Rb), both in its metallic and oxidized states, on the corrosion resistance of steels has direct relevance to the storage of radioactive Kr-85 containing nuclear materials. Kr-85 undergoes a ß-decay to generate stable, metallic Rb. The literature to date is sparse, contained primarily in company reports and divisive as to whether Rb metal, or its oxides, can lead to corrosion of steel materials. In both the direct storage of Kr-85 and used nuclear fuel that will contain Kr-85, if Rb is corrosive then continual decay of the Kr-85 will generate an increasingly corrosive environment for the steel canister it is stored in. This review was written to consolidate the known data to date on Rb corrosion of steels to identify consistencies, contradictions and remaining gaps in our understanding of the Rb : steel system.

krypton, corrosion, steel, Rubidium↗
DOE OSTI2020

Atom probe tomography characterization of ion and neutron irradiated Alloy 800H

Alloy 800H is a high temperature and creep resistant alloy, and is considered as a candidate alloy for use in Generation IV nuclear reactor systems. To clarify the alloy's behavior under irradiation, samples were subjected to single ion irradiations up to 20 dpa at 440°C, a dual ion/He irradiation to 17 dpa at 460°C, and a neutron irradiation to 17 dpa at 385°C. The irradiated microstructures were characterized using atom probe tomography to complement previously published transmission electron microscopy studies. Additionally, after single ion irradiation, sparse fine Al and Ti clusters were observed after 1dpa, while high number densities of nanoscale Ni-Al-Ti clusters, Cr-Ti rich carbides, and Si-decorated dislocation loops developed after 10 and 20 dpa. The microstructure formed during single ion irradiation exhibited a strong depth dependence. No significant differences were observed between the single and dual ion irradiations. While Ni-Al-Ti clusters, Cr-Ti rich carbides, and Si-decorated dislocation loops were also observed after neutron irradiation, the neutron-irradiated microstructure differed from those found in the ion-irradiated samples.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Characterization of the molten salt FMgNaK through ab initio molecular dynamics and experimental density measurements

The selection of a salt species for use in a molten salt reactor MSR is a key part of any MSR design. However, many salts have sparse or no thermophysical property data sets, especially those with higher melting points. One such salt is the eutectic mixture of (NaF)0.345(KF)0.59(MgF 2 )0.065, or FMgNaK, and is explored here through ab initio molecular dynamics (AIMD) simulation (1023-1273K) as well as experimental measurement of the liquid density using the Archimedean method (973-1223K). Predicted densities from AIMD are similar to experimentally-measured densities, though with more noise, suggesting that the small scale necessitated by AIMD simulations may be problematic for simulating FMgNaK in particular. Here the coefficient of thermal expansion is predicted from simulation, and salt structure is characterized. Mg-F-Mg chaining is observed in the salt network, though the low concentration of Mg inhibits chaining on the scale that is observed in (LiF) 0.67 (BeF 2 ) 0.33 , or FLiBe.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Determining the effects of U/Pu ratio on subsolidus phase transitions in U-Pu-Zr metallic fuel alloys

Here, ternary alloys consisting primarily of uranium, plutonium, and zirconium (U-Pu-Zr) are among the leading candidate fuel systems considered for fast spectrum nuclear reactors. Despite historical operation data from the testing of U-Pu-Zr rods in the Experimental Breeder Reactor-II, considerable uncertainty about the evolution of phases and microstructure across the ternary composition space exists. Due to sluggish kinetics and other difficulties in handling metal actinide specimens, quantitative measurements of phase-transitions in U-Pu-Zr alloys remain sparse in scientific literature, with most investigators reporting either phase-transition temperatures or phase identification data, but not both from the same specimens. The purpose of this paper is to critically compare experimental and calculated phase transition data and correlate with the microstructure and phase characterization data of as-cast and annealed U-Pu-Zr alloys. Phase transition peaks were measured using differential scanning calorimetry in the subsolidus regions (723-948 K) of three ternary U-Pu-Zr alloys with the same zirconium concentration but various U/Pu ratios. Overlapping peaks were deconvoluted using a Frazier-Suzuki peak fitting algorithm, and the critical peak temperatures and enthalpies were calculated. In general, increasing concentrations of Pu were associated with enhanced thermal stability of the body-centered cubic γ phase upon both heating and cooling. Experimental phase transition temperatures in this study tended to agree well with the predictions of the established ternary phase diagrams and other reported phase transition temperatures in literature. Additionally, the TAF-ID thermodynamic database was used to calculate a U-Pu-40 at.% Zr pseudobinary diagram as well as ternary diagrams from 773 to 973 K. The equilibrium phase transition temperatures tended to be considerably lower than measured peak temperatures upon both heating and cooling. Recommendations for improving the quality of data in future U-Pu-Zr characterization studies are also discussed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Radiation induced athermal diffusivity in uranium mononitride

Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuel considered for light water reactors and advanced reactor designs. Properties like self- and fission gas diffusivity need to be better understood, given that they influence key fuel performance phenomena such as fission gas swelling and release. In particular, the radiation induced athermal (D 3 ) diffusivity remains challenging to accurately predict and has only been sparsely characterized in UN, despite its importance as it likely governs diffusion at the low temperatures this high-thermal-conductivity fuel form may operate. Molecular Dynamics simulations are used to estimate the mean square displacement induced by a primary knock-on atom (PKA) with a given kinetic energy. These results are combined with the PKA energy distributions obtained from binary collision approximation calculations to obtain the displacement due to a particular fission fragment. Finally, this is combined with experimental fission fragment yields to determine the displacement due to an average fission event and, thus, express the athermal diffusivity as a function of the fission rate density. These results are in excellent agreement with available experimental data. In conclusion, a particular importance is given to the understanding and the quantification of the variability of these results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗