Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical modeling”

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

Thermodynamics of order and randomness in dopant distributions inferred from atomically resolved imaging

Abstract Exploration of structure-property relationships as a function of dopant concentration is commonly based on mean field theories for solid solutions. However, such theories that work well for semiconductors tend to fail in materials with strong correlations, either in electronic behavior or chemical segregation. In these cases, the details of atomic arrangements are generally not explored and analyzed. The knowledge of the generative physics and chemistry of the material can obviate this problem, since defect configuration libraries as stochastic representation of atomic level structures can be generated, or parameters of mesoscopic thermodynamic models can be derived. To obtain such information for improved predictions, we use data from atomically resolved microscopic images that visualize complex structural correlations within the system and translate them into statistical mechanical models of structure formation. Given the significant uncertainties about the microscopic aspects of the material’s processing history along with the limited number of available images, we combine model optimization techniques with the principles of statistical hypothesis testing. We demonstrate the approach on data from a series of atomically-resolved scanning transmission electron microscopy images of Mo x Re 1- x S 2 at varying ratios of Mo/Re stoichiometries, for which we propose an effective interaction model that is then used to generate atomic configurations and make testable predictions at a range of concentrations and formation temperatures.

25 ENERGY STORAGE↗
DOE OSTI2022

Solvation Thermodynamics of Solutes in Water and Ionic Liquids Using the Multiscale Solvation-Layer Interface Condition Continuum Model

We report molecular assembly processes are generally driven by thermodynamic properties in solutions. Atomistic modeling can be very helpful in designing and understanding complex systems, except that bulk solvent is very inefficient to treat explicitly as discrete molecules. The SLIC/CDC multiscale model combines continuum solvent electrostatics based on the solvent layer interface condition (SLIC) with statistical thermodynamic models for hydrogen bonding and nonpolar modes: cavity formation, dispersion interactions, combinatorial mixing (CDC). The SLIC/CDC model predicts Gibbs energies of solvation for a database of 500 solutes in water with average accuracy better than 1 kcal/mol both for experimental measurements and for explicit-solvent molecular dynamics simulations. The separate SLIC/CDC energy mode values agree quantitatively with those computed from explicit-solvent molecular dynamics. The SLIC/SASA multiscale model combines the SLIC continuum electrostatic model with the solvent accessible surface area (SASA) nonpolar energy mode. The SLIC/SASA model predicts Gibbs energies of solvation with better than 1.4 kcal/mol average accuracy in aqueous systems and better than 1.6 kcal/mol average accuracy in ionic liquids. Both models predict solvation entropies, and are the first implicit-solvation models capable of predicting solvation heat capacities.

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

Random insights into the complexity of two-dimensional tensor network calculations

Projected entangled pair states (PEPS) offer memory-efficient representations of some quantum many-body states that obey an entanglement area law and are the basis for classical simulations of ground states in two-dimensional (2d) condensed matter systems. However, rigorous results show that exactly computing observables from a 2d PEPS state is generically a computationally hard problem. Yet approximation schemes for computing properties of 2d PEPS are regularly used, and empirically seen to succeed, for a large subclass of (“not too entangled”) condensed matter ground states. Adopting the philosophy of random matrix theory, in this work, we analyze the complexity of approximately contracting a 2d random PEPS by exploiting an analytic mapping to an effective replicated statistical mechanics model that permits a controlled analysis at a large bond dimension. Through this statistical-mechanics lens, we argue that (i) although approximately sampling wave-function amplitudes of random PEPS faces a computational-complexity phase transition above a critical bond dimension, and (ii) one can generically efficiently estimate the norm and correlation functions for any finite bond dimension. Furthermore, these results are supported numerically for various bond-dimension regimes. It is an important open question whether the above results for random PEPS apply more generally also to PEPS representing physically relevant ground states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

β — -delayed fission in the coupled quasiparticle random-phase approximation plus Hauser-Feshbach approach

Beta-delayed neutron emission and β-delayed fission (βdf) probabilities were calculated for heavy, neutronrich nuclei using the Los Alamos coupled quasiparticle random-phase approximation plus Hauser-Feshbach (QRPA + HF) approach. In this model, the compound nucleus is initially populated by β decay and is followed through subsequent statistical decays taking into account competition among neutrons, γ rays, and fission. The primary output of these calculations includes branching ratios along with neutron and γ -ray spectra. Here, we find a relatively large region of heavy nuclides where the probability of βdf is near 100%. For a subset of nuclei near the neutron drip line, delayed neutron emission and the probability of fission are both large, which leads to the possibility of multichance βdf (mc-βdf). We comment on prospective neutron-rich nuclei that could be probed by future experimental campaigns and provide a full table of branching ratios in ASCII format in the Supplemental Material for use in various applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

System-level performance and degradation of 21 GWDC of utility-scale PV plants in the United States

We assess the performance of a fleet of 411 utility-scale (i.e., >5 MWAC and ground-mounted) photovoltaic (PV) projects totaling 21.1 GWDC (16.3 GWAC) of capacity, which achieved commercial operations in the United States from 2007 to 2016. This fleet of projects contributed more than 50% of all solar electricity generated in the United States in 2017. Using detailed information on individual project characteristics, in conjunction with modeled irradiance data, we assess the extent to which actual first-year performance has lived up to both modeled and stated expectations. We then analyze system-level performance degradation in subsequent years by employing a “fixed effects” regression model to statistically isolate the impact of age on system performance. We find that this fleet of utility-scale PV projects has generally lived up to ex ante expectations for first-year performance but that subsequent system-level degradation—found to be −1.3%/year (±0.2%)—has, on average, been worse than both ex ante expectations (commonly −0.5%/year) and results from past studies (ranging from −0.8%/year to −1.0%/year). We emphasize that −1.3%/year is a system-level estimate that captures more than just module degradation (e.g., including soiling, balance of plant degradation, and downtime for maintenance and/or other events). A side analysis of a variety of project characteristics suggests that system-level degradation rates tend to be of lower magnitude among newer projects and larger projects and at sites with lower long-term average temperatures.

14 SOLAR ENERGY↗
DOE OSTI2022

An approximate inertial manifold (AIM) based closure for turbulent flows

A closure model for turbulent flows is developed based on a dynamical system theory. An appropriately discretized formulation of the governing equations is considered for this process. The key ingredient is an approximation of the system’s attractor, where all the trajectories in phase space are confined. This approximate inertial manifold based approach provides a path to track trajectories of the system in a lower-dimensional subspace. Unlike conventional coarse-graining approaches, the turbulent field is decomposed into resolved and unresolved dynamics using the properties of the governing equations. The novelty of the approach relies on the reconstruction of the unresolved field constrained by the governing equations. A posteriori tests for homogeneous isotropic turbulence and the Kuramoto–Sivashinsky equation show promising results for considerable dimension reduction with strong convergence properties. The proposed model outperforms the dynamic Smagorinsky model, and the computational overhead is competitive with similar approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

First simultaneous measurement of the γ-ray and neutron emission probabilities in inverse kinematics at a heavy-ion storage ring

The probabilities for γ-ray and particle emission as a function of the excitation energy of a decaying nucleus are valuable observables for constraining the ingredients of the models that describe the deexcitation of nuclei near the particle emission threshold. These models are essential in nuclear astrophysics and applications. In this paper, we have for the first time simultaneously measured the γ-ray and neutron emission probabilities of 208 Pb. The measurement was performed in inverse kinematics at the Experimental Storage Ring (ESR) of the GSI/FAIR facility, where a 208 Pb beam interacted through the 208 Pb (p,p') reaction with a hydrogen gas jet target. Instead of detecting the γ rays and neutrons emitted by 208 Pb, we detected the heavy beamlike residues produced after γ and neutron emission. These heavy residues were fully separated by a dipole magnet of the ESR and were detected with outstanding efficiencies. The comparison of the measured probabilities with model calculations has allowed us to test and select different descriptions of the γ-ray strength function and the nuclear level density available in the literature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Multiphenomenology explosion monitoring (MultiPEM): a general framework for data interpretation and yield estimation

SUMMARY An underground nuclear explosion (UNE) couples mechanical energy into crustal rock, which propagates as seismic and acoustic waves. These different physical phenomena transport, by different pathways, to standoff detectors at varying distances. The transport pathways attenuate the original signal but in different ways. Enabled by correct statistical weighting, signal attenuation models can be used to combine these disparate sensor data to estimate the yield of an UNE. Contemporaneous statistical models, used in yield estimation, can be improved with an advanced partition of error for these physical signal propagation models. We present an advanced multivariate approach to error modelling of multiphenomenology physical signatures. In addition to measurement error, our error model represents physical model biases as random with a physics-based covariance structure. To illustrate this proposed framework, we demonstrate the estimation of explosion yield using openly available seismic and acoustic data from chemical single-point explosions.

Williams, Brian J.↗
DOE OSTI2022

Pseudogauge dependence of the spin polarization and of the axial vortical effect

The mean spin polarization vector of spin 1/2 particles in a relativistic fluid at local thermal equilibrium (LTE) in different pseudogauges (PGs), i.e., with different choices for the decomposition of orbital and spin angular momentum, is obtained. The spin polarization obtained in the canonical PG differs from the one obtained in the Belinfante PG. It is found that this difference can not be written by replacing the thermal vorticity with the spin potential in the usual polarization formula. Other PG choices affect the contribution of thermal shear to the spin polarization. So, the choice of a PG affects the predictions for the polarization measured in heavy-ion collisions. In general, it is shown that the Wigner function of a noninteracting Dirac field at LTE is affected by the PG transformations. Explicit expressions of the mean axial current are also calculated and it is found that even the axial vortical effect conductivity depends on the PG.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Discussion on “Saving Storage in Climate Ensembles: A Model-Based Stochastic Approach”

Abstract We thank the authors for this interesting paper that highlights important ideas and concepts for the future of climate model ensembles and their storage, as well as future uses of stochastic emulators. Stochastic emulators are particularly relevant because of the statistical nature of climate model ensembles, as discussed in previous work of the authors (Castruccio et al. in J Clim 32:8511–8522, 2019; Hu and Castruccio in J Clim 34:8409–8418, 2021). We thank the authors for sharing of some of their data with us in order to illustrate this discussion. In the following, in Sect. 1 we discuss alternative techniques currently used and studied, namely lossy compression and ideas emerging from the climate modeling community, that could feed the discussion on ensemble and storage. In that section, we also present numerical results of compression performed on the data shared by the authors. In Sect. 2, we discuss the current statistical model proposed by the authors and its context. We discuss other potential uses of stochastic emulators in climate and Earth modeling.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Foundational Research for Fusion Systems Safety Assessment Overview

Presentation detailing work package to develop comprehensive fusion device safety models over the course of the next three years, including references to a blanket, divertor and magnet design, and statistical reliability models

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Computational and Experimental Mechanistic Insights into the Ethanol-to-Butanol Upgrading Reaction over MgO

The mechanism of ethanol upgrading to higher products is still under debate, especially regarding intermediate species and hydrogenation and dehydrogenation steps. In this work, we conducted a combined theoretical and experimental approach to contribute to this discussion. For such, detailed electronic structure density functional theory calculations (aiming at probing density of states, infrared spectra, geometric parameters, charge densities, and reaction energetics) and diffuse reflectance infrared Fourier transform spectroscopy experiments were carried out revealing the relevance of an appropriate combination of reactive surface sites to support the formation of several intermediates that are formed in the C-C coupling over MgO. The roles of Mg and O sites were also studied under an electronic perspective and different geometrical arrangements. We found that a kink configuration was the most adequate for ethanol to 1-butanol upgrading. Our calculations also gave us arguments to propose distinct reaction routes, whose mutual predominance would depend upon reaction temperature. At temperatures up to 573 K, the so-called β-route, which goes through scission of a Cβ-H bond and formation of an oxametallacycle-like intermediate, would dominate the coupling, whereas at higher temperatures, up to 673 K, a more usual Guerbet mechanism, via an aldol coupling step and then consecutive hydrogenations, would be expected. Here, the theoretical conclusions were followed by a careful experimental strategy using sequential experimental planning techniques in order to estimate accurate parameters with the lowest possible experimental load. Information from these different sources were coupled to develop a mathematical model for the rate of the ethanol upgrading reaction, using a Langmuir-Hinshelwood-Hougen-Watson approach. The developed and statistically validated model adequately described the experimental data at 673 K and 1.1 bar total pressure for ethanol partial pressures in the range from 0 to 20 kPa.

09 BIOMASS FUELS↗
DOE OSTI2026

Tula: Optimizing Time, Cost, and Generalization in Distributed Large-Batch Training

Distributed training increases the number of batches processed per iteration either by scaling-out (adding more nodes) or scaling-up (increasing the batch-size). However, the largest configuration does not necessarily yield the best performance. Horizontal scaling introduces additional communication overhead, while vertical scaling is constrained by computation cost and device memory limits. Thus, simply increasing the batch-size leads to diminishing returns: training time and cost decrease initially but eventually plateaus, creating a knee-point in the time/cost vs. batch-size pareto curve. The optimal batch-size therefore depends on the underlying model, data and available compute resources. Large batches also suffer from worse model quality due to the well-known “generalization gap”. In this paper, we present Tula, an online service that automatically optimizes time, cost, and convergence quality for large-batch training of convolutional models. It combines parallel-systems modeling with statistical performance prediction to identify the optimal batchsize. Tula predicts training time and cost within 7.5−14% error across multiple models, and achieves up to 20× overall speedup and improves test accuracy by ≈9% on average over standard large-batch training on various vision tasks, thus successfully mitigating the generalization gap and accelerating training at the same time.

Tyagi, Sahil [ORNL] (ORCID:0009000783144745)↗
DOE OSTI2023

Quantifying model prediction sensitivity to model-form uncertainty

Computational and mathematical models are essential to understanding complex systems and phenomena. However, when developing such models, limited knowledge and/or resources necessitates the use of simplifying assumptions. It is therefore crucial to quantify the impact of such simplifying assumptions on the reliability and accuracy of resulting model predictions. This work develops a first-of-its-kind approach to quantify the impact of physics modeling assumptions on predictions. Here, we leverage the emerging field of model-form uncertainty (MFU) representations, which are parameterized modifications to modeling assumptions, in combination with grouped Sobol’ indices to quantitatively measure an assumption’s importance. Specifically, we compute the grouped Sobol’ index for the MFU representation’s parameters as a single importance measure of the assumption for which the MFU representation characterizes uncertainty. To ensure this approach is robust to the subjective choice of how to parameterize a MFU representation, we establish bounds for the difference between sensitivity results for two different MFU representations based on differences in model prediction statistics. The capabilities associated with this approach are demonstrated on three exemplar problems: an upscaled subsurface contaminant transport problem, ablation modeling for hypersonic flight, and nuclear waste repository modeling. We found that our grouped approach is able to assess the impact of modeling assumptions on predictions and offers computational advantages over classical Sobol’ index computation while providing more interpretable results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Probabilistic projections of the Amery Ice Shelf catchment, Antarctica, under conditions of high ice-shelf basal melt

Abstract. Antarctica's Lambert Glacier drains about one-sixth of the ice from the East Antarctic Ice Sheet and is considered stable due to the strong buttressing provided by the Amery Ice Shelf. While previous projections of the sea-level contribution from this sector of the ice sheet have predicted significant mass loss only with near-complete removal of the ice shelf, the ocean warming necessary for this was deemed unlikely. Recent climate projections through 2300 indicate that sufficient ocean warming is a distinct possibility after 2100. This work explores the impact of parametric uncertainty on projections of the response of the Lambert–Amery system (hereafter “the Amery sector”) to abrupt ocean warming through Bayesian calibration of a perturbed-parameter ice-sheet model ensemble. We address the computational cost of uncertainty quantification for ice-sheet model projections via statistical emulation, which employs surrogate models for fast and inexpensive parameter space exploration while retaining critical features of the high-fidelity simulations. To this end, we build Gaussian process (GP) emulators from simulations of the Amery sector at a medium resolution (4–20 km mesh) using the Model for Prediction Across Scales (MPAS)-Albany Land Ice (MALI) model. We consider six input parameters that control basal friction, ice stiffness, calving, and ice-shelf basal melting. From these, we generate 200 perturbed input parameter initializations using space filling Sobol sampling. For our end-to-end probabilistic modeling workflow, we first train emulators on the simulation ensemble and then calibrate the input parameters using observations of the mass balance, grounding line movement, and calving front movement with priors assigned via expert knowledge. Next, we use MALI to project a subset of simulations to 2300 using ocean and atmosphere forcings from a climate model for both low- and high-greenhouse-gas-emission scenarios. From these simulation outputs, we build multivariate emulators by combining GP regression with principal component dimension reduction to emulate multivariate sea-level contribution time series data from the MALI simulations. We then use these emulators to propagate uncertainty from model input parameters to predictions of glacier mass loss through 2300, demonstrating that the calibrated posterior distributions have both greater mass loss and reduced variance compared to the uncalibrated prior distributions. Parametric uncertainty is large enough through about 2130 that the two projections under different emission scenarios are indistinguishable from one another. However, after rapid ocean warming in the first half of the 22nd century, the projections become statistically distinct within decades. Overall, this study demonstrates an efficient Bayesian calibration and uncertainty propagation workflow for ice-sheet model projections and identifies the potential for large sea-level rise contributions from the Amery sector of the Antarctic Ice Sheet after 2100 under high-greenhouse-gas-emission scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Intermittency and Dissipative Structures Arising from Relativistic Magnetized Turbulence

Abstract Kinetic simulations of relativistic turbulence have significantly advanced our understanding of turbulent particle acceleration. Recent progress has highlighted the need for an updated acceleration theory that can account for particle acceleration within the plasma’s coherent structures. Here, we investigate how intermittency modeling connects statistical fluctuations in turbulence to regions of high-energy dissipation. This connection is established by employing a generalized She–Leveque model to characterize the exponents ζ p for the structure functions S p ∝ l ζ p . The fitting of the scaling exponents provides us with a measure of the codimension of the dissipative structures, for which we subsequently determine the filling fraction. We perform our analysis for a range of magnetizations σ and relative fluctuation amplitudes δ B 0 / B 0 . We find that increasing values of σ and δ B 0 / B 0 allow the turbulent cascade to break sheetlike structures into smaller regions of dissipation that resemble chains of flux ropes. However, as their dissipation measure increases, the dissipative regions become less volume filling. With this work, we aim to inform future turbulent acceleration theories that incorporate particle energization from interactions with coherent structures within relativistic turbulence.

Astronomy & Astrophysics↗