Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surrogates”

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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Uncertainty quantification for Multiphase-CFD simulations of bubbly flows: a machine learning-based Bayesian approach supported by high-resolution experiments

In this paper, we developed a machine learning-based Bayesian approach to inversely quantify and reduce the uncertainties of multiphase computational fluid dynamics (MCFD) simulations for bubbly flows. The proposed approach is supported by high-resolution two-phase flow measurements, including those by double-sensor conductivity probes, high-speed imaging, and particle image velocimetry. Local distributions of key physical quantities of interest (QoIs), including the void fraction and phasic velocities, are obtained to support the Bayesian inference. In the process, the epistemic uncertainties of the closure relations are inversely quantified while the aleatory uncertainties from stochastic fluctuations of the system are evaluated based on experimental uncertainty analysis. The combined uncertainties are then propagated through the MCFD solver to obtain uncertainties of the QoIs, based on which probability-boxes are constructed for validation. The proposed approach relies on three machine learning methods: feedforward neural networks and principal component analysis for surrogate modeling, and Gaussian processes for model form uncertainty modeling. The whole process is implemented within the framework of an open-source deep learning library PyTorch with graphics processing unit (GPU) acceleration, thus ensuring the efficiency of the computation. The results demonstrate that with the support of high-resolution data, the uncertainties of MCFD simulations can be significantly reduced. The proposed approach has the potential for other applications that involve numerical models with empirical parameters.

42 ENGINEERING↗
DOE OSTI2022

Physics-informed machine learning assisted uncertainty quantification for the corrosion of dissimilar material joints

Jointing techniques like the Self-Piercing Riveting (SPR), Resistance Spot Welding (RSW) and Rivet-Weld (RW) joints are used for mass production of dissimilar material joints due to their high performance, short cycle time, and adaptability. However, the service life and safety usage of these joints can be largely impacted by the galvanic corrosion due to the difference in equilibrium potentials between the metals with the presence of electrolyte. Here, in this paper, we focus on Al-Fe galvanic corrosion and develop physics-informed machine learning based surrogate model for statistical corrosion analysis, which enables the reliability analysis of dissimilar material joints under corrosion environment. In this study, a physics-based finite element (FE) corrosion model has been developed to simulate the galvanic corrosion between a Fe cathode and an Al anode. Geometric and environmental factors including crevice gap, roughness of anode, conductivity, and the temperature of the electrolyte are investigated. Further, a thorough Uncertainty Quantification (UQ) analysis is conducted for the overall corrosion behavior of the Fe-Al joints. It is found that the electrolyte conductivity has the largest effects on the material loss and needs to be managed closely for better corrosion control. This will help in designing and manufacturing joints with improved corrosion performance.

42 ENGINEERING↗
DOE OSTI2025

Bayesian Entropy Neural Networks for physics-aware prediction

This article addresses the need for deep learning models to integrate well-defined constraints into their outputs, driven by their application in surrogate models, learning with limited data and partial information, and scenarios requiring flexible model behavior to incorporate non-data sample information. We introduce Bayesian Entropy Neural Networks (BENN), a framework grounded in Maximum Entropy (MaxEnt) principles, designed to impose constraints on Bayesian Neural Network (BNN) predictions. BENN is capable of constraining not only the predicted values but also their derivatives and variances, ensuring a more robust and reliable model output. To achieve simultaneous uncertainty quantification and constraint satisfaction, we employ the method of multipliers approach. This allows for the concurrent estimation of neural network parameters and the Lagrangian multipliers associated with the constraints. Our experiments, spanning diverse applications such as beam deflection modeling and microstructure generation, demonstrate the effectiveness of BENN. The results highlight significant improvements over traditional BNNs and showcase competitive performance relative to contemporary constrained deep learning methods.

14 SOLAR ENERGY↗
DOE OSTI2021

MULTI-MODAL global surveillance methodology for predictive and on-demand characterization of localized processes using cube satellite platforms and deep learning techniques

This paper presents the work completed towards the development of a multi-modal global surveillance methodology using cube satellite (CubeSat) platforms and novel data analysis techniques. A CubeSat system equipped with adequate sensors and data analytics capabilities can autonomously characterize various phenomena of interest on the Earth’s surface. CubeSats are advantageous over conventional satellites in certain remote monitoring applications because of their reduced construction costs (due to the availability of commercially-off-the-shelf components) and are easier to launch. The CubeSat surveillance system developed in this paper focused on phenomena of interest surrounding the nuclear fuel cycle in support of nuclear non-proliferation and emergency response. To observe the phenomena, a constellation of 3U and 6U CubeSats deployed from the ISS with adequate components was chosen. Four different sensor configurations were identified for remote sensing: panchromatic/multispectral in the visible and near-infrared spectrum, multispectral in infrared spectrum, hyperspectral in infrared spectrum, and multispectral in ultraviolet spectrum. While a panchromatic/multispectral sensor configuration has CubeSat flight heritage at the required spatial resolutions, the other three sensor types need future 3 development to meet signature and system requirements. Once each sensor onboard the CubeSat system collects data on a target of interest, the onboard computers would then apply the deep learning-based characterization methodology developed in this paper to identify phenomena. Four surrogate datasets containing representative simplified “images” were created for each sensor type to train the characterization methodology. A convolutional neural network was applied to each dataset and produced recall rates for the phenomena between 89.7% - 99.3% and precision rates between 92.3% - 99.9%. Each phenomenon’s presence probability from each network is then combined into a final characterization solution for a target area. This paper covers multiple interdisciplinary areas to develop the foundation for a CubeSat surveillance system focused on phenomena surrounding the nuclear fuel cycle.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2022

Spectroscopic signatures and oxidation characteristics of nanosecond laser-induced cerium plasmas

Improving technologies related to the wide area environmental sampling of nuclear materials supports the nuclear nonproliferation mission of preventing the proliferation of nuclear weapons by monitoring nuclear weapons tests and detecting undeclared nuclear fuel cycle activities. Standoff, laser-based detection techniques such as laser-induced breakdown spectroscopy have the potential to offer robust, field-deployable methods that provide rapid element-specific and phase identifiable measurements over a wide range of materials. This work aims to elucidate the effects of atmospheric conditions and oxidation reactions on the highly complex and transient spectroscopic signatures of laser-induced plutonium surrogate plasmas. Time-resolved spectra of nanosecond laser ablation cerium plasmas were measured using laser-induced breakdown spectroscopy in a range of atmospheres containing low to high concentrations of oxygen. Here, the growth of strong CeO molecular emission bands was observed in the visible spectrum, where it was shown that the persistence of CeO is reduced from around 60 μs to 50 μs in oxygen rich atmospheric environments. To further investigate the growth and depletion of CeO in the laser-produced plasma, ratios of CeO-to-Ce emission were generated using integrated intensities corresponding to the Q-branch of the CeO D 1 -X 1 transitions and numerous strong atomic Ce peaks. It was determined that the fastest rate of formation of CeO in argon occurred for moderate oxygen mass fractions between 0.10 and 0.15 while the ratios were reduced at higher oxygen mass fractions (i.e., Y O 2 = 0.20) due to competing oxidation reactions and lower plasma temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Elemental fractionation in aerosol laser-induced breakdown spectroscopy with nanosecond and femtosecond laser ablation

For the purpose of realizing atmospheric radiation plume tracking, laser-induced breakdown spectroscopy (LIBS) was performed to characterize Cs atomic emissions from airborne nanoparticles in a binary particle matrix using nanosecond and femtosecond laser ablation. As a surrogate to multi-element nuclear fallout particulates, both Cs and Na solutions were mixed with Cu at a 1:1.45 M concentration ratio and aerosolized to study elemental fractionation in atomic emissions, revealing Na I resonance emissions are enhanced whereas Cs I resonance emissions are reduced with Cu present despite sharing similar electronic and chemical properties. Femtosecond laser ablation reduced elemental fractionation effects by more than 50% while retaining the same temporal fractionation trends for Cs and Na. The origin of elemental fractionation in aerosol LIBS was further investigated using Al and Ti aerosols which revealed that the particle matrix induces changes to the local plasma temperature, population of atomic states, and subsequent molecular association and emission.

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

The confluence of machine learning and multiscale simulations

Multiscale modeling has a long history of use in structural biology, as computational biologists strive to overcome the time- and length-scale limits of atomistic molecular dynamics. Contemporary machine learning techniques, such as deep learning, have promoted advances in virtually every field of science and engineering and are revitalizing the traditional notions of multiscale modeling. Deep learning has found success in various approaches for distilling information from fine-scale models, such as building surrogate models and guiding the development of coarse-grained potentials. However, perhaps its most powerful use in multiscale modeling is in defining latent spaces that enable efficient exploration of conformational space. In conclusion, this confluence of machine learning and multiscale simulation with modern high-performance computing promises a new era of discovery and innovation in structural biology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Neural network representation of electronic structure from ab initio molecular dynamics

Despite their rich information content, electronic structure data amassed at high volumes in ab initio molecular dynamics simulations are generally under-utilized. Here we introduce a transferable high-fidelity neural network representation of such data in the form of tight-binding Hamiltonians for crystalline materials. This predictive representation of ab initio electronic structure, combined with machine-learning boosted molecular dynamics, enables efficient and accurate electronic evolution and sampling. When it is applied to a one-dimension charge-density wave material, carbyne, we are able to compute the spectral function and optical conductivity in the canonical ensemble. The spectral functions evaluated during soliton-antisoliton pair annihilation process reveal significant renormalization of low-energy edge modes due to retarded electron-lattice coupling beyond the Born–Oppenheimer limit. The availability of an efficient and reusable surrogate model for the electronic structure dynamical system will enable calculating many interesting physical properties, paving the way to previously inaccessible or challenging avenues in materials modeling.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Predictive modeling of Néel temperature in austenitic alloys using CALPHAD and data analytics

The Néel temperature is a crucial yet often overlooked parameter in calculating the stacking fault energy (SFE) of austenitic alloys. Several empirical equations have been proposed to estimate the Néel temperature of austenitic alloys, which are then used to calculate the SFE and explain deformation mechanisms. However, these empirical equations, typically derived using linear regression algorithms, are often simplistic and may fail to capture the complex interactions among multiple alloying elements that influence the Néel temperature. Moreover, their applicability is usually limited to specific compositional ranges. In this study, we propose a CALPHAD based approach and develop a surrogate decision tree based regression model capable of capturing the interactions among multiple alloying elements to predict the Néel temperature. Predictions from both the CALPHAD approach and the regression model show close agreement with experimental measurements reported in the literature. In conclusion, the implications of accurate Néel temperature predictions on the calculated SFE and deformation mechanisms are also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Impact of large-scale effects on mass transfer and concentration polarization in Reverse Osmosis membrane systems

Here, we present well-resolved computational fluid dynamics simulations of a large-scale reverse osmosis membrane-spacer configuration (1 m). Our computational model solves the flow and transport equations with variable solute-dependent properties. We utilize a high resolution computational mesh to resolve all relevant length scales associated with spacer-induced mixing and thin concentration boundary layers. An important contribution of this work is the development of a modified mass-transfer correlation that accounts for the development of the concentration boundary layer along the channel. A set of 2D axisymmetric simulations were performed for a spiral wound module layer with varying cross-flow conditions and spacer diameters which indicate a significant entrance length effect for concentration profile development at lower flow rates while mixing effects dominate at higher flow rates. The mass-transfer correlations at higher flow rates compare well with published correlations while a surrogate model for Sherwood number was obtained that depends on an additional similarity variable that accounted for entrance length effects at lower flow rates. Finally, a large-scale membrane-spacer design relevant to high-pressure reverse osmosis is studied with a non-uniform arrangement of spacers, which indicate a substantial saving in pressure drop (~40%) compared to traditional uniformly spaced pattern with minor variations (~2%) in concentration polarization, product water quality (~1%) and water recovery (~7%) compared to a uniform spacer pattern.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Enabling scientific machine learning in MOOSE using Libtorch

A neural-network-based machine learning interface has been developed for the Multiphysics Object-Oriented Simulation Environment (MOOSE). The interface relies on Libtorch, the C++ front-end of PyTorch, and enables an online interaction between modern machine learning algorithms and all the existing simulation, modeling, and analysis processes available in MOOSE. New capabilities in MOOSE include the native generation and training of artificial neural networks together with options to load pretrained neural networks in TorchScript format. Furthermore, the MOOSE stochastic tools module (MOOSE-STM) has been enhanced with neural network-based surrogate and reduced-order model generation options for efficient stochastic analyses. Lastly, a reinforcement learning capability has been added to MOOSE-STM for the interactive control and optimization of complex multiphysics problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Influence of the air-cooled supercritical carbon dioxide recompression cycle design-point on molten salt power tower techno-economic performance

Here, this paper investigates the techno-economic performance of a molten salt power tower integrated with an air-cooled supercritical carbon dioxide recompression cycle. Prior research into the cycle design point showed trade-offs between efficiency, cost, and temperature difference across the heat input, resulting in trade-offs between the design-point generation of the molten salt power tower system and the cost of sensible-heat thermal energy storage. Prior research into cycle off-design showed that colder ambient temperature designs have higher efficiencies on cold days but lower efficiency and capacity degradation on hot days. We present a framework to evaluate how cycle design, cost, and corresponding off-design performance affect the techno-economic performance of the molten salt power tower system. We use a steady-state off-design cycle model assuming fixed shaft speeds, inventory control, and air-cooler fan power control and integrate the cycle model with the molten salt power tower model from the National Renewable Energy Laboratory's System Advisor Model. First, we present a surrogate model that uses a heuristic approach to correct the interpolation when the cycle reaches capacity limitations at ambient temperatures hotter than design. Next, we evaluate how the recompression cycle design-point ambient temperature, air-cooler approach temperature, air-cooler fan power, and recuperator conductance influence cycle and system design metrics. Finally, we complete annual simulations using an hourly "duck curve" pricing schedule and present annual performance metrics to evaluate the cumulative effects of the cycle design parameters and corresponding off-design performance on molten salt power tower techno-economics.

14 SOLAR ENERGY↗
DOE OSTI2026

Heliostat sizing methodology for concentrating solar thermal industrial process heat projects

This study presents a method to obtain a heliostat size that minimizes the levelized cost of heat (LCOH) of a heliostat-based concentrating solar thermal system for applications of solar heating for industrial processes at operating temperatures from 565 to 1550°C. The method extends prior work by embedding a routine for system design that obtains near-optimal subsystem sizes, increasing the fidelity of drive cost functions, and adding an optical performance model to supplement the previously developed cost models, which we update to reflect current pricing trends. An illustrative business case is developed for Daggett, California, targeting specified annual thermal energy outputs of 50 to 400 GWh th . Optical performance is modeled using verified estimates from the literature. A surrogate heliostat cost model, derived from commercial heliostat designs and scaled for production volume, installation, and operations and maintenance costs, is used to develop cost functions. Results show that heliostat size strongly affects the LCOH, producing a characteristic U-shaped trend with a robust near-optimal window of 7-20 m 2 ; the heliostat size producing the lowest project cost in our study grows slightly as the project size increases, and is reduced as the operating temperature increases. The findings in this study are consistent with the general trend of smaller heliostats being deployed at existing projects for high-temperature industrial process heat and reflect the significant reduction in power electronics and other per-heliostat costs. The methodology we propose is general and can be tailored to revised cost curves as the technology continues to evolve.

14 SOLAR ENERGY↗
DOE OSTI2020

Reactions during conversion of simplified low-activity waste glass feeds

The mechanisms that affect the incorporation of 99Tc, a volatile radioactive component of concern, into glass melt during vitrification of low-activity waste (LAW) are being investigated to develop the method to increase the retention of 99Tc in glass waste form. Previous studies with simulated LAW glass feeds (slurry mixture of liquid waste and chemical/mineral additives) demonstrated that the early stage feed-to-glass conversion reactions below 800°C are critical for the Re (used as a nonradioactive surrogate of 99Tc) retention in glass. To examine the effect of feed composition on the feed-to-glass conversion reactions, simplified systems containing major LAW components (NaNO3 and NaOH) and representative additive components (SiO2 and H3BO3) were designed and tested. The ratio of H3BO3 to NaNO3 was varied in three-component system without NaOH and that of NaOH to NaNO3 was varied in the four-component system at a fixed H3BO3 to NaNO3 ratio. As a first step of testing with simplified feeds, this study applied thermal analyses and phase characterization of the reacting feeds, which were performed without the addition of Re, to investigate the evolution of salt phases during slurry drying process and upon heating of dried feeds.

Jin, Tongan↗
DOE OSTI2021

Genotoxicity evaluation using primary hepatocytes isolated from rhesus macaque ( Macaca mulatta )

Non-human primates (NHPs) have played a vital role in fundamental, pre-clinical, and translational studies because of their high physiological and genetic similarity to humans. We report a method to isolate primary hepatocytes from the livers of rhesus macaques (Macaca mulatta) after in situ whole liver perfusion. Isolated primary macaque hepatocytes (PMHs) were treated with various compounds known to have different pathways of genotoxicity/carcinogenicity and the resulting DNA damage was evaluated using the high-throughput CometChip assay. The comet data were quantified using benchmark dose (BMD) modeling and the BMD50 values for treatments of PMHs were compared with those generated from primary human hepatocytes (PHHs) in our previous study. The results showed that despite varying CYP450 enzyme activities, PMHs had the same sensitivity and specificity as PHHs in detecting four indirect-acting (i.e., requiring metabolic activation) and seven direct-acting genotoxicants/carcinogens, as well as five non-carcinogens that are negative or equivocal for genotoxicity in vivo. The BMD50 estimates and their confidence intervals revealed species differences for DNA damage potency, especially for direct-acting compounds. The present study provides a practical method for maximizing the use of animal tissues by isolating primary hepatocytes from NHPs. Our data support the use of PMHs as a reliable surrogate of PHHs for evaluating the genotoxic hazards of chemical substances for humans.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Implicit neural representations for experimental steering of advanced experiments

Scattering measurements using electrons, neutrons, or photons are essential for obtaining microscopic insights into materials. However, limited facility availability and high-dimensional scattering data necessitate more efficient experimental steering techniques. Here, we report a machine learning method that guides scattering data collection and facilitates real-time estimation of model parameters, given a reliable forward model to simulate experimental signals. We employ implicit neural representations as efficient surrogates that link model parameters with simulated spectroscopies. This enables a Bayesian optimal experimental design framework to estimate the probability distributions of parameters from high-dimensional scattering data. We demonstrate the proposed method using inelastic neutron scattering with simulated and real experimental data, highlighting the method’s ability to provide real-time parameter estimation with quantified uncertainties and to deliver informed experimental guidance that reduces experimental time while maximizing scientific output. This approach paves the way for accelerated discoveries in condensed matter through scattering measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Reduced Order Modeling conditioned on monitored features for response and error bounds estimation in engineered systems

Reduced Order Models (ROMs) form essential tools across engineering domains by virtue of their function as surrogates for computationally intensive digital twinning simulators. Although purely data-driven methods are available for ROM construction, schemes that allow to retain a portion of the physics tend to enhance the interpretability and generalization of ROMs. However, physics-based techniques can adversely scale when dealing with nonlinear systems that feature parametric dependencies. This study introduces a generative physics-based ROM that is suited for nonlinear systems with parametric dependencies and is additionally able to provide numerical error bounds associated with the respective estimates. A main contribution of this work is the conditioning of these parametric ROMs to features that can be derived from monitoring measurements, feasibly in an online fashion. This is contrary to most existing ROM schemes, which remain restricted to the prescription of the physics-based, and usually a priori unknown, system parameters. Our work utilizes conditional Variational Autoencoders to continuously map the required reduction bases to a feature vector extracted from limited output measurements, while additionally allowing for a probabilistic assessment of the ROM-estimated Quantities of Interest. An auxiliary task using a neural network-based parametrization of suitable probability distributions is introduced to re-establish the link with physical model parameters. We verify the proposed scheme on a series of simulated case studies incorporating effects of geometric and material nonlinearity under parametric dependencies related to system properties and input load characteristics.

Conditional VAEs↗
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↗