Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ensemble models”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quantifying Northern High Latitude Gross Primary Productivity (GPP) Using Carbonyl Sulfide (OCS)

The northern high latitude (NHL, 40°N to 90°N) is where the second peak region of gross primary productivity (GPP) other than the tropics. The summer NHL GPP is about 80% of the tropical peak, but both regions are still highly uncertain. Carbonyl sulfide (OCS) provides an important proxy for photosynthetic carbon uptake. Here we optimize the OCS plant uptake fluxes across the NHL by fitting atmospheric concentration simulation with the GEOS-CHEM global transport model to the aircraft profiles acquired over Alaska during NASA's Carbon in Arctic Reservoirs Vulnerability Experiment (2012–2015). We use the empirical biome-specific linear relationship between OCS plant uptake flux and GPP to derive the six plant uptake OCS fluxes from different GPP data. Such GPP-based fluxes are used to drive the concentration simulations. We evaluate the simulations against the independent observations at two ground sites of Alaska. The optimized OCS fluxes suggest the NHL plant uptake OCS flux of -247 Gg S year -1 , about 25% stronger than the ensemble mean of the six GPP-based OCS fluxes. GPP-based OCS fluxes systematically underestimate the peak growing season across the NHL, while a subset of models predict early start of season in Alaska, consistent with previous studies of net ecosystem exchange. The OCS optimized GPP of 34 PgC yr -1 for NHL is also about 25% more than the ensembles mean from six GPP data. Further work is needed to fully understand the environmental and biotic drivers and quantify their rate of photosynthetic carbon uptake in Arctic ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Accurate estimation of angular power spectra for maps with correlated masks

A common procedure when analyzing maps of the cosmic microwave background (CMB) or other cosmological signals is the need to remove ("mask") regions of the maps that are heavily contaminated, e.g., by non-cosmological foreground emission. After applying such a mask, one must account for its effect when inferring statistical properties of interest, such as the angular power spectrum of the field in the original map. A widely used approach to correct for such mask-induced effects was presented by Hivon et al. (2002), now widely known as the "MASTER" formalism. However, it is often the case that the map and mask are correlated in some way, such as point source masks used in CMB analyses, which have nonzero correlation with CMB secondary anisotropy fields and other mm-wave sky signals. In such situations, the MASTER approach gives biased results, as it assumes that the unmasked map and mask have zero correlation. While such effects have been discussed before with regard to specific physical models, here we derive a completely general formalism for any case where the map and mask are correlated. We show that our result ("reMASTERed") reconstructs ensemble-averaged angular power spectra to effectively exact precision, with significant improvements over traditional estimators for cases where the map and mask are correlated. An important consequence of our result is that for maps with correlated masks, it is no longer possible to invert a simple equation to obtain the true power spectrum from the observed (masked) power spectrum. Instead, our result necessitates the use of forward modeling from theory space into the observable domain of the masked power spectrum. We publicly release our software implementation of these results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Data-guided Multi-Map variables for ensemble refinement of molecular movies

Driving molecular dynamics simulations with data-guided collective variables offer a promising strategy to recover thermodynamic information from structure-centric experiments. In this study, the three-dimensional electron density of a protein, as it would be determined by cryo-EM or x-ray crystallography, is used to achieve simultaneously free-energy costs of conformational transitions and refined atomic structures. Unlike previous density-driven molecular dynamics methodologies that determine only the best map-model fits, our work employs the recently developed Multi-Map methodology to monitor concerted movements within equilibrium, non-equilibrium, and enhanced sampling simulations. Construction of all-atom ensembles along the chosen values of the Multi-Map variable enables simultaneous estimation of average properties, as well as real-space refinement of the structures contributing to such averages. Using three proteins of increasing size, we demonstrate that biased simulation along the reaction coordinates derived from electron densities can capture conformational transitions between known intermediates. The simulated pathways appear reversible with minimal hysteresis and require only low-resolution density information to guide the transition. The induced transitions also produce estimates for free energy differences that can be directly compared to experimental observables and population distributions. The refined model quality is superior compared to those found in the Protein Data Bank. We find that the best quantitative agreement with experimental free-energy differences is obtained using medium resolution density information coupled to comparatively large structural transitions. Practical considerations for probing the transitions between multiple intermediate density states are also discussed.

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

Predicting oxidation damage of ultra high-temperature carbide ceramics in extreme environments using machine learning

Determining the oxidation resistance of UHTC carbides in extreme environments is challenging theoretically and experimentally due to the high dimensional complexity of influencing variables and intricate testing setups. Herein we demonstrate the use of machine learning (ML) models trained with experimental literature data to predict the oxide thickness of UHTC carbides exposed to air based on composition, mean grain size, relative densification, holding time, and temperature. A multi-dimensional database with 76 occurrences is created containing experimental results of Hf, Zr, and Ta carbides plus additives. In this study, the preprocessed database is then used to train ML models to predict their oxidation behavior. The trained model predicts the oxidation damage in the form of an average oxide thickness in UHTC carbides with a Mean Absolute Error (MAE) of ±65.45 μm for samples in the testing set that developed thicknesses up to 1000 μm. The model successfully predicted oxidation damage for a recession rate lower than 60 μm/min. It is noticed that the ensemble method MAE is increased to ±134.34 μm while forecasting the oxidation of samples with a recession rate higher than the threshold. The unprecedented approach is a novel way to predict the damage through the oxidation of carbide compounds before processing for a smarter design with room for improvement.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Thermofield Theory for Finite-Temperature Electronic Structure

Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. Furthermore, they also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.

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

Studying phonon coherence with a quantum sensor

Nanomechanical oscillators offer numerous advantages for quantum technologies. Their integration with superconducting qubits shows promise for hardware-efficient quantum error-correction protocols involving superpositions of mechanical coherent states. Limitations of this approach include mechanical decoherence processes, particularly two-level system (TLS) defects, which have been widely studied using classical fields and detectors. In this manuscript, we use a superconducting qubit as a quantum sensor to perform phonon number-resolved measurements on a piezoelectrically coupled phononic crystal cavity. This enables a high-resolution study of mechanical dissipation and dephasing in coherent states of variable size ($\overline{n}$ ≃ 1 – 10 phonons). We observe nonexponential relaxation and state size-dependent reduction of the dephasing rate, which we attribute to TLS. Using a numerical model, we reproduce the dissipation signatures (and to a lesser extent, the dephasing signatures) via emission into a small ensemble (N = 5) of rapidly dephasing TLS. Our findings comprise a detailed examination of TLS-induced phonon decoherence in the quantum regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Spectral broadening from turbulence in multiscale lower hybrid current drive simulations

Abstract The scattering of lower hybrid (LH) waves due to scrape-off layer filaments is investigated. It is revealed that scattering can account for the LH spectral gap without any ad hoc modification to the wave-spectrum. This is shown using a multiscale simulation approach which allows, for the first time, the inclusion of full-wave scattering physics in ray-tracing/Fokker–Planck calculations. In this approach, full-wave scattering probabilities are calculated for a wave interacting with a statistical ensemble of filaments. These probabilities are coupled to ray-tracing equations using radiative transfer theory. This allows the modeling of scattering along the entire ray-trajectory, which can be important in the multi-pass regime. Simulations are conducted for LH current drive (LHCD) in Alcator C-Mod, resulting in excellent agreement with experimental current and hard x-ray profiles. A region in filament parameter space is identified in which the impact of scattering on LHCD is saturated. Such a state coincides with experimental LHCD measurements, suggesting saturation indeed occurs in C-Mod, and therefore the exact statistical properties of the filaments are not important.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Linking repeat lidar with Landsat products for large scale quantification of fire-induced permafrost thaw settlement in interior Alaska

The permafrost–fire–climate system has been a hotspot in research for decades under a warming climate scenario. Surface vegetation plays a dominant role in protecting permafrost from summer warmth, thus, any alteration of vegetation structure, particularly following severe wildfires, can cause dramatic top–down thaw. A challenge in understanding this is to quantify fire-induced thaw settlement at large scales (>1000 km 2 ). In this study, we explored the potential of using Landsat products for a large-scale estimation of fire-induced thaw settlement across a well-studied area representative of ice-rich lowland permafrost in interior Alaska. Six large fires have affected ~1250 km 2 of the area since 2000. We first identified the linkage of fires, burn severity, and land cover response, and then developed an object-based machine learning ensemble approach to estimate fire-induced thaw settlement by relating airborne repeat lidar data to Landsat products. The model delineated thaw settlement patterns across the six fire scars and explained ~65% of the variance in lidar-detected elevation change. Our results indicate a combined application of airborne repeat lidar and Landsat products is a valuable tool for large scale quantification of fire-induced thaw settlement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Statistical Treatment of Inverse Problems Constrained by Differential Equations-based Models with Stochastic Terms

This paper introduces a statistical treatment of inverse problems constrained by models with stochastic terms. The solution of the forward problem is given by a distribution represented numerically by an ensemble of simulations. The goal is to formulate the inverse problem, in particular the objective function, to find the closest forward distribution (i.e., the output of the stochastic forward problem) that best explains the distribution of the observations in a certain metric. We use proper scoring rules, a concept employed in statistical forecast verification, namely energy, variogram, and hybrid (i.e., combination of the two) scores. We study the performance of the proposed formulation in the context of two applications: a coefficient field inversion for subsurface flow governed by an elliptic partial differential equation with a stochastic source and a parameter inversion for power grid governed by differential-algebraic equations. In both cases we show that the variogram and the hybrid scores produce better parameter inversion results than does the energy score, whereas the energy score leads to better probabilistic predictions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines↗
DOE OSTI2022

Generative adversarial networks for ensemble projections of future urban morphology

As city planners design and adapt cities for future resilience and intelligence, interactions among neighborhood morphological development with respect to changes in population and resultant built infrastructure's impact on the natural environment must be considered. For deep understanding of these interactions, explicit representation of future neighborhoods is necessary for future city modeling. Generative Adversarial Networks (GANs) have been shown to produce spatially accurate urban forms at scales representing entire cities to those at neighborhood and single building scale. Here we demonstrate a GAN method for generating an ensemble of possible new neighborhoods given land use characteristics and designated neighborhood type.

Dumas, Melissa↗
DOE OSTI2023

Collisional-energy-cascade model for nonthermal velocity distributions of neutral atoms in plasmas

Nonthermal velocity distributions with much greater tails than a Maxwellian have been observed for radical atoms in plasmas for a long time. Historically, such velocity distributions have been modeled by a two-temperature Maxwell distribution. In this paper, I propose a model based on collisional energy cascade, which has been studied in the field of granular materials. In the collisional energy cascade, a particle ensemble undergoes energy input at the high-energy region, entropy production by elastic collisions among particles, and energy dissipation. For radical atoms, energy input may be caused by the Franck-Condon energy of molecular dissociation or charge-exchange collision with hot ions, and the input energy is eventually dissipated by collisions with the walls. Here, I show that the steady-state velocity distribution in the collisional energy cascade is approximated by the generalized Mittag-Leffler distribution, which is a one-parameter extension of the Maxwell distribution. This parameter indicates the degree of the nonthermality and is related to the relative importance of energy dissipation over entropy production. This model is compared with a direct molecular dynamics simulation for simplified gaseous systems with energy input, as well as some experimentally observed velocity distributions of light radicals in plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Understanding and control of Zener pinning via phase field and ensemble learning

Zener pinning refers to the dispersion of fine particles which influences grain size distribution via movement of grain boundaries in a polycrystalline material. Grain size distribution in polycrystals has a significant impact on their properties including physical, chemical, mechanical, and optical to name a few. We explore the use of Phase-field modeling and machine-learning techniques to understand and improve the control of grain size distribution via Zener pinning in polycrystalline materials. We develop a machine learning model that determines the relative importance of various parameters to exercise microstructure control via Zener pinning. Our workflow combines high-throughput phase-field simulations and machine learning to address the computational bottlenecks associated with large-scale simulations as well as identify features necessary for microstructure control in polycrystals. A random forest (RF) regression model was developed to predict grain sizes based on five Phase-field model parameters, achieving an average prediction error of 0.72 nm for the training data and 1.44 nm for the test data. The importance of the input parameters is analyzed using the SHapley Additive exPlanations (SHAP) approach which reveals that diffusivity, volume fraction, and particle diameter are the most important parameters in determining the final grain size. These findings will allow us to select the best second-phase particles, optimize grain size distributions and thus design microstructures with the desired properties. The developed method is a highly versatile and generalizable approach that can be used to assess the combined effects of individual features in the presence of multiple variables.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Uncertainty quantification in machine learning for engineering design and health prognostics: A tutorial

On top of machine learning (ML) models, uncertainty quantification (UQ) functions as an essential layer of safety assurance that could lead to more principled decision making by enabling sound risk assessment and management. The safety and reliability improvement of ML models empowered by UQ has the potential to significantly facilitate the broad adoption of ML solutions in high-stakes decision settings, such as healthcare, manufacturing, and aviation, to name a few. In this tutorial, we aim to provide a holistic lens on emerging UQ methods for ML models with a particular focus on neural networks and the applications of these UQ methods in tackling engineering design as well as prognostics and health management problems. Towards this goal, we start with a comprehensive classification of uncertainty types, sources, and causes pertaining to UQ of ML models. Next, we provide a tutorial-style description of several state-of-the-art UQ methods: Gaussian process regression, Bayesian neural network, neural network ensemble, and deterministic UQ methods focusing on spectral-normalized neural Gaussian process. Established upon the mathematical formulations, we subsequently examine the soundness of these UQ methods quantitatively and qualitatively (by a toy regression example) to examine their strengths and shortcomings from different dimensions. Then, we review quantitative metrics commonly used to assess the quality of predictive uncertainty in classification and regression problems. Afterward, we discuss the increasingly important role of UQ of ML models in solving challenging problems in engineering design and health prognostics. In conclusion, two case studies with source codes available on GitHub are used to demonstrate these UQ methods and compare their performance in the life prediction of lithium-ion batteries at the early stage (case study 1) and the remaining useful life prediction of turbofan engines (case study 2).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Solid-state platform for cooperative quantum dynamics driven by correlated emission

While traditionally regarded as an obstacle to quantum coherence, recent breakthroughs in quantum optics have shown that the dissipative interaction of a qubit with its environment can be leveraged to protect quantum states and synthesize many-body entanglement. Inspired by this progress, here we set the stage for the—as yet uncharted—exploration of analogous cooperative phenomena in hybrid solid-state platforms. We develop a comprehensive formalism for the quantum many-body dynamics of an ensemble of solid-state spin defects interacting with the magnetic field fluctuations of a common solid-state reservoir. Our framework applies to any solid-state reservoir whose fluctuating spin, pseudospin, or charge degrees of freedom generate magnetic fields. To understand whether correlations induced by dissipative processes can play a relevant role in a realistic experimental setup, we apply our model to a qubit array interacting via the spin fluctuations of a ferromagnetic bath. Our results show that the low-temperature collective relaxation rates of the qubit ensemble can display clear signatures of super- and subradiance, i.e., forms of cooperative dynamics traditionally achieved in atomic ensembles. We find that the solid-state analog of these cooperative phenomena is robust against spatial disorder in the qubit ensemble and thermal fluctuations of the magnetic reservoir, providing a route for their feasibility in near-term experiments. Finally, our work lays the foundation for a multiqubit approach to quantum sensing of solid-state systems and the direct generation of many-body entanglement in spin-defect ensembles. Furthermore, we discuss how the tunability of solid-state reservoirs opens up novel pathways for exploring cooperative phenomena in regimes beyond the reach of conventional quantum optics setups.

NV centers↗
DOE OSTI2020

Accelerated Climate Model for Energy (ACME): Ocean and Sea-Ice Processes

The primary goal of this project was to understand the role of the ocean and sea-ice components in the evolution of the climate in both a high-resolution fully-coupled 1850 preindustrial control and present day transient climate simulations run using the Energy Exascale Earth System Model (E3SM), formerly known as the Accelerated Climate Model for Energy (ACME). Two distinct approaches to the initialization of the present day transients were tested: in the first, an ensemble of three transients used atmospheric reanalysis forced 0.1-degree ocean/sea-ice states, while a fourth transient was branched off the preindustrial control. Biases from the control were inherited by its successor transient, in particular, sea surface temperatures (SSTs) at high latitudes were too warm and too much northward meridional heat transport occurred in the North Atlantic Ocean. SSTs in the Northern Hemisphere and Atlantic meridional heat transport increased in the transients initialized from the ocean/sea-ice simulation over a 40-50 year period to largely resemble the transient branched off the control simulation. Excessive absorbed solar radiation (ASR) is considered to be responsible for the SST biases as well as the formation of mode waters in lighter density classes than in reality both in the Southern Ocean and the North Atlantic Ocean subtropical gyre. These mode waters ventilate the thermocline waters to the north of their formation sites explaining the excessive Atlantic meridional heat transport and some of the warming in the high-latitude North Atlantic Ocean after multiple decades.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Modeling the Joint Effects of Vegetation Characteristics and Soil Properties on Ecosystem Dynamics in a Panama Tropical Forest

Abstract In tropical forests, both vegetation characteristics and soil properties are important not only for controlling energy, water, and gas exchanges directly but also determining the competition among species, successional dynamics, forest structure and composition. However, the joint effects of the two factors have received limited attention in Earth system model development. Here we use a vegetation demographic model, the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM), ELM‐FATES, to explore how plant traits and soil properties affect tropical forest growth and composition concurrently. A large ensemble of simulations with perturbed vegetation and soil hydrological parameters is conducted at the Barro Colorado Island, Panama. The simulations are compared against observed carbon, energy, and water fluxes. We find that soil hydrological parameters, particularly the scaling exponent of the soil retention curve ( B sw ), play crucial roles in controlling forest diversity, with higher B sw values (>7) favoring late successional species in competition, and lower B sw values (1 ∼ 7) promoting the coexistence of early and late successional plants. Considering the additional impact of soil properties resolves a systematic bias of FATES in simulating sensible/latent heat partitioning with repercussion on water budget and plant coexistence. A greater fraction of deeper tree roots can help maintain the dry‐season soil moisture and plant gas exchange. As soil properties are as important as vegetation parameters in predicting tropical forest dynamics, more efforts are needed to improve parameterizations of soil functions and belowground processes and their interactions with aboveground vegetation dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Enabling machine learning-ready HPC ensembles with Merlin

With the growing complexity of computational and experimental facilities, many scientific researchers are turning to machine learning (ML) techniques to analyze large scale ensemble data. With complexities such as multi-component workflows, heterogeneous machine architectures, parallel file systems, and batch scheduling, care must be taken to facilitate this analysis in a high performance computing (HPC) environment. Here, we present Merlin, a workflow framework to enable large ML-friendly ensembles of scientific HPC simulations. By augmenting traditional HPC with distributed compute technologies, Merlin aims to lower the barrier for scientific subject matter experts to incorporate ML into their analysis. As a producer–consumer workflow model, Merlin enables multi-machine, cross-batch job, dynamically allocated yet persistent workflows capable of utilizing surge-compute resources. Key features of Merlin are a flexible HPC-centric interface, low per-task overhead, multi-tiered fault recovery, and a hierarchical sampling algorithm that allows for $\mathscr{O}$(N) task execution and $\mathscr{O}$(N ln N) task queuing to ensembles of millions of tasks. In addition to Merlin’s design, we test the algorithm’s performance in an HPC center and demonstrate the ability to enqueue 40 million simulations in 100 s, with a 30 millisecond per-task overhead that is independent of ensemble size. Finally, we describe some example applications that Merlin has enabled on leadership-class HPC resources, such as the ML-augmented optimization of nuclear fusion experiments and the calibration of infectious disease models to study the progression of and possible mitigation strategies for COVID-19.

97 MATHEMATICS AND COMPUTING↗