Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied statistics”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Develop and Connect TRISO Failure Analysis and Uncertainty Quantification to Fission Product Release Calculation Capability

The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering-scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. Stress-based failure probability has been developed and analyzed to assess the integrity of tri-structural isotropic (TRISO) fuel particles during fuel life cycles. While simple, stress-based approaches to failure probability leveraging the Weibull statistical distribution entails a number of drawbacks when stress concentration occurs near crack tips, including finite element mesh size dependency. In this report, we use an interaction integral approach to the computation of stress intensity factors in functionally graded materials (FGM) for axisymmetric models. The inner pyrolytic carbon (IPyC) cracking induced silicon carbide (SiC) failure is one of the dominated failure modes in TRISO failure analysis. In this study, we consider a crack in the IPyC layer perpendicular to the SiC layer. The interface between these two TRISO layers is considered to be porous, which we simulate considering a transition of mechanical properties over the porous length. These aspects are considered in the computation of stress intensity factor (SIF) from a fracture mechanics approach and compared with the known stress-based failure probability approach.

42 ENGINEERING↗
DOE OSTI2026

First constraints from marked angular power spectra with Subaru Hyper Suprime-Cam Survey First-Year Data

We present the first application of marked power spectra to weak lensing data, using maps from the Subaru Hyper Suprime-Cam Year 1 (HSC-Y1) survey. Marked convergence fields, constructed by weighting the convergence field with non-linear functions of its smoothed version, are designed to encode higher-order information while remaining computationally tractable. Using simulations tailored to the HSC-Y1 data, we test three mark functions that up- or down-weight different density environments. Our results show that combining multiple types of marked auto and cross-spectra improves constraints on the clustering amplitude parameter S8≡σ8Ωm/0.3 by ≈43 per cent compared to standard two-point power spectra. When applied to the HSC-Y1 data, this translates into a constraint on S8=0.807±0.024⁠. We assess the sensitivity of the marked power spectra to systematics, including baryonic effects, intrinsic alignment, photometric redshifts, and multiplicative shear bias. We note that some of the additional information introduced by the marked field originates from scales smaller than the scale cut, and is partly Gaussian in nature. This does not invalidate our systematic tests. These results demonstrate the promise of marked statistics as a practical and powerful tool for extracting non-Gaussian information from weak lensing surveys.

Cowell, Jessica A. [Oxford U.; Tokyo U., IPMU] (OR↗
DOE OSTI2023

A Machine Learning Framework to Deconstruct the Primary Drivers for Electricity Market Price Events

As the electricity grid is moving towards a 100% Renewable Energy Source Bulk Power Grid, the overall operations of the power system operations and electricity markets are changing. The electricity markets are not only dispatching resources economically but also taking into account various controllable actions like renewable curtailment, transmission congestion mitigation, and energy storage optimization to make sure the grid is operating reliably. As a result, price formations in electricity markets have become quite complex. Traditional root cause analysis and statistical approaches are rendered inapplicable to analyze and infer the main drivers behind price formation in the modern grid and markets with variable renewable energy (VRE). In this paper, we propose a machine learning analysis framework to deconstruct some primary drivers for price formation in modern electricity markets with high renewable energy and the outcomes can be utilized for various critical aspects of market design, renewable dispatch and curtailment, operations, and cyber-security applications. The framework can be applied to any ISO or market data and in this paper it is applied to open-source publicly available datasets from California Independent System Operator (CAISO) and ISO New England.

machine learning (ML), electricity markets, Renewa↗
DOE OSTI2025

Muon momentum reconstruction in ICARUS-T600 LArTPC via Multiple Scattering

The SBN program at Fermilab investigates the possible existence of eV-scale sterile neutrinos. This work focuses on reconstructing muon momentum in $\nu$-Ar interactions within the ICARUS-T600 LArTPC detector using Multiple Coulomb Scattering (MCS), the only viable method for non-contained muons. Two algorithms are evaluated: the Gran Sasso algorithm (2D, Collection view) and the MicroBooNE algorithm (3D angles), applied to 2391 simulated and 2255 real stopping muon tracks (0.4-1 GeV/c). Performance is assessed via bias and resolution, using range-based momentum as a reference. Two enhanced versions of the Gran Sasso algorithm are also developed: one combining all three wire plane views, and one using 3D angles. Results show that the MicroBooNE algorithm performs slightly better, with resolution down to 14\% in simulation and 16\% in real data. These improvements significantly boost $\nu_\mu$ CC statistics, enhancing sensitivity to $\nu_\mu$ disappearance and 3+1 oscillation signals.

Chiello, Giovanni [Pisa U.] (ORCID:000900064571071↗
DOE OSTI2021

Prototype equation of state for phase transition of confined fluids based on the generalized van der Waals partition function

A simple self-consistent prototype equation of state (EOS) based on the Generalized van der Waals (vdW) partition function has been demonstrated to describe the phase transition of simple fluids in nanopores with uniform size. Different from those commonly presented in the literature, the new EOS does not need an auxiliary equation that is conventionally applied to provide the capillary pressure derived from surface tension. Lastly, the encouraging performance of the EOS calls for further extension to applications with more complex fluids and porous media.

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

Effect of Groundwater Age and Recharge Source on Nitrate Concentrations in Domestic Wells in the San Joaquin Valley

Nitrate is one of the most abundant contaminants in groundwater globally, in the United States, and in California (CA). We studied well construction information, water chemistry, stable isotopes, and noble gases to understand how groundwater travel time and recharge source and mechanism control nitrate concentrations in domestic wells in the San Joaquin Valley (SJV), CA, a large semiarid, irrigated agricultural region. Using nonparametric statistics, we find a decreasing trend in nitrates with groundwater travel time and well depth. Samples collected from wells that are closer to rivers and that show indications of river water recharge, either low recharge temperature or low δ 18 O signature, have lower concentrations of nitrates than samples with isotopic signatures indicating mixed source or local precipitation recharge. The curbing effect of river water recharge on nitrate concentrations in domestic wells is similar for direct river recharge and water applied as irrigation. This suggests that irrigation with river water also has a diluting effect that reduces the concentration of nitrate found in groundwater. This conclusion supports the idea that flood-managed aquifer recharge may be considered for remediation of groundwater nitrate when designing replenishment of aquifers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Predicting Flow in Fracture Networks With Quantum Algorithms

Uncertainty quantification plays a crucial role in the modeling of subsurface flow. For instance, uncertainties in the properties of geologic fracture networks significantly impact flow, requiring numerous simulations to accurately estimate quantities of interest. However, each simulation is computationally expensive because it requires solving a large linear system to capture features that involve both small and large fractures. An example is in percolation, where the interaction of many small fractures (which cumulatively can have a large surface area) with the rock matrix must be modeled precisely. Quantum computing is an emerging tool with the potential to address this issue. Quantum algorithms offer a significant speedup in solving linear systems, achieving efficiencies that are challenging to match with classical approaches. These classical approaches include direct solvers, such as LU decomposition, and iterative methods, notably preconditioned conjugate gradient, commonly used in subsurface modeling to solve large sparse systems. However, applying quantum algorithms to geologic fracture flow requires careful attention to algorithmic and problem-specific constraints to fully realize this quantum advantage. In this work we describe a quantum algorithm for generalized Monte Carlo applications with a quadratic speedup over the classical approaches which can be combined with the quantum speedup, currently under investigation, for solving quantum linear systems for subsurface flow. We show that for quantum algorithms the computational cost of estimating a quantity of interest for a statistical ensemble of networks is roughly the same as that of a single realization, essentially implying that one can get uncertainty quantification for free.

58 GEOSCIENCES↗
DOE OSTI2022

Using kernel-based statistical distance to study the dynamics of charged particle beams in particle-based simulation codes

Measures of discrepancy between probability distributions (statistical distance) are widely used in the fields of artificial intelligence and machine learning. We describe how certain measures of statistical distance can be implemented as numerical diagnostics for simulations involving charged-particle beams. Related measures of statistical dependence are also described. The resulting diagnostics provide sensitive measures of dynamical processes important for beams in nonlinear or high-intensity systems, which are otherwise difficult to characterize. Here, the focus is on kernel-based methods such as maximum mean discrepancy, which have a well-developed mathematical foundation and reasonable computational complexity. Several benchmark problems and examples involving intense beams are discussed. While the focus is on charged-particle beams, these methods may also be applied to other many-body systems such as plasmas or gravitational systems.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Emissions and Flame Structure Assessment of Aeroengine Micromixing Injectors for Lean Direct Injection of Hydrogen and Hydrogen/Natural Gas Blends

Abstract The present study focuses on whether adapting low emissions liquid fuel injection concepts can be applied to gaseous mixtures. The goal is to develop an injector with specific geometric characteristics that will achieve the lowest emissions. This is conducted by evaluating emissions and flame diagnostics on a set of test hardware that has been developed for combustion of pure hydrogen and hydrogen/methane mixtures. The current test hardware adapts lean direct injection concepts from aeroengine injectors to a ground-based engine based on target flows associated with a Solar Turbines Industrial Engine gas turbine. A statistical test plan was developed to assess these geometric characteristics using a three level Box Behnken design. The test plan investigates the effect of injector flow splits, air swirl, fuel swirl, pressure drop, preheat temperature, fuel composition, and the adiabatic flame temperature on the flame structure and emissions performance. A total of 16 injectors were considered for a total of 688 test points for emission and imaging data. Emissions data were reported on a mass basis in ng/J and chemical reactor network analysis were conducted to further interpret emissions behavior. Images were also captured to analyze the flame structure features of interest including reaction volume, flame length, flame spread, and flame width. With these results, a predictive model for emissions was made to conclude the most optimal injector configuration is to have the highest level of air split (more volume of air through the injector plate versus the injector) and air swirl angle for both methane and hydrogen gas mixtures. The predictive emissions level for all of these cases were below 10 ng/J.

08 HYDROGEN↗
DOE OSTI2020

Computing Potential of the Mean Force Profiles for Ion Permeation Through Channelrhodopsin Chimera, C1C2

Umbrella sampling, coupled with a weighted histogram analysis method (US-WHAM), can be used to construct potentials of mean force (PMFs) for studying the complex ion permeation pathways of membrane transport proteins. Despite the widespread use of US-WHAM, obtaining a physically meaningful PMF can be challenging. Here, we provide a protocol to resolve that issue. Then, we apply that protocol to compute a meaningful PMF for sodium ion permeation through channelrhodopsin chimera, C1C2, for illustration.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A best-practice guide to predicting plant traits from leaf-level hyperspectral data using partial least squares regression

Partial least squares regression (PLSR) modelling is a statistical technique for correlating datasets, and involves the fitting of a linear regression between two matrices. One application of PLSR enables leaf traits to be estimated from hyperspectral optical reflectance data, facilitating rapid, high-throughput, non-destructive plant phenotyping. This technique is of interest and importance in a wide range of contexts including crop breeding and ecosystem monitoring. The lack of a consensus in the literature on how to perform PLSR means that interpreting model results can be challenging, applying existing models to novel datasets can be impossible, and unknown or undisclosed assumptions can lead to incorrect or spurious predictions. We address this lack of consensus by proposing best practices for using PLSR to predict plant traits from leaf-level hyperspectral data, including a discussion of when PLSR is applicable, and recommendations for data collection. Further, we provide a tutorial to demonstrate how to develop a PLSR model, in the form of an R script accompanying this manuscript. This practical guide will assist all those interpreting and using PLSR models to predict leaf traits from spectral data, and advocates for a unified approach to using PLSR for predicting traits from spectra in the plant sciences.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Determination of jet calibration and energy resolution in proton–proton collisions at $\sqrt{s} = 8~\hbox {TeV}$ using the ATLAS detector

The jet energy scale, jet energy resolution, and their systematic uncertainties are measured for jets reconstructed with the ATLAS detector in 2012 using proton–proton data produced at a centre-of-mass energy of 8 TeV with an integrated luminosity of 20fb -1 . Jets are reconstructed from clusters of energy depositions in the ATLAS calorimeters using the anti-k t algorithm. A jet calibration scheme is applied in multiple steps, each addressing specific effects including mitigation of contributions from additional proton–proton collisions, loss of energy in dead material, calorimeter non-compensation, angular biases and other global jet effects. The final calibration step uses several in situ techniques and corrects for residual effects not captured by the initial calibration. Furthermore, these analyses measure both the jet energy scale and resolution by exploiting the transverse momentum balance in γ + jet, Z + jet, dijet, and multijet events. A statistical combination of these measurements is performed. In the central detector region, the derived calibration has a precision better than 1% for jets with transverse momentum 150 GeV < p T < 1500 GeV, and the relative energy resolution is (8.4 ± 0.6)% for p T =100 GeV and (23 ± 2)% for p T =20 GeV. The calibration scheme for jets with radius parameter R = 1.0, for which jets receive a dedicated calibration of the jet mass, is also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Data Evaluation of Actinide Cross Sections: 230 Th and 231 Th

A new evaluation of the ENDL cross section set for Thorium (Z=90) is developed using the TALYS statistical model cross section code. The primary goal of this effort is to produce an evaluation that attempts to match as closely as possible fission cross sections developed through surrogate reaction techniques on actinide targets 230 Th and 231 Th. This evaluation effort and the processing needed to render its results into data libraries is a necessary step in making the efforts of nuclear experimentalists useful to the broad community of researchers engaged in simulations of nuclear fusion for basic and applied science. Another aspect, verification and validation against various AGEX experiments, is also presented. The end-product is an updated library that includes the latest measurements of fission cross sections for comparison against those measured via traditional techniques. All the steps in the evaluation, processing, validation and verification, and library release are described in the following sections. For completeness, the appendix contains all the parameters used in the TALYS cross section evaluation for neutrons incident on 230 Th and 231 Th.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Extreme Wind and Waves in U.S. East Coast Offshore Wind Energy Lease Areas

The Outer Continental Shelf along the U.S. east coast exhibits abundant wind resources and is now a geographic focus for offshore wind deployments. This analysis derives and presents expected extreme wind and wave conditions for the sixteen lease areas that are currently being developed. Using the homogeneous ERA5 reanalysis dataset it is shown that the fifty-year return period wind speed (U50) at 100 m a.s.l. in the lease areas ranges from 29.2 to 39.7 ms -1 . After applying corrections to account for spectral smoothing and averaging period, the associated pseudo-point U50 estimates are 34 to 46 ms -1 . The derived uncertainty in U50 estimates due to different distributional fitting is smaller than the uncertainty associated with under-sampling of the interannual variability in annual maximum wind speeds. It is shown that, in the northern lease areas, annual maximum wind speeds are generally associated with intense extratropical cyclones rather than cyclones of tropical origin. Extreme wave statistics are also presented and indicate that the 50-year return period maximum wave height may substantially exceed 15 m. From this analysis, there is evidence that annual maximum wind speeds and waves frequently derive from the same cyclone source and often occur within a 6 h time interval.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Full‐dimensional coupled‐channel statistical approach to atom‐triatom systems and applications to H/D + O 3 reaction

Abstract The statistical quantum model (SQM), which assumes that the reactivity is controlled by entrance/exit channel quantum capture probabilities, is well suited for chemical reactions with a long‐lived intermediate complex. In this work, a time‐independent coupled‐channel implementation of the SQM approach is developed for atom‐triatom systems in full dimensionality. As SQM treats the capture dynamics quantum mechanically, it is capable of handling quantum effects such as tunneling. A detailed study of the H/D + O 3 capture dynamics was performed by applying the newly developed SQM method on an accurate global potential energy surface. Agreement with previous ring polymer molecular dynamics (RPMD) results on the same potential energy surface is excellent except for very low temperatures. The SQM results are also in reasonably good agreement with available experimental rate coefficients. The strong H/D kinetic isotope effect underscores the dominant role of quantum tunneling under an entrance channel barrier at low temperatures.

Yang, Dongzheng↗
DOE OSTI2025

Navigating the Noise: Bringing Clarity to ML Parameterization Design With O $\boldsymbol{\mathcal{O}}$(100) Ensembles

Abstract Machine‐learning (ML) parameterizations of subgrid processes (here of turbulence, convection, and radiation) may one day replace conventional parameterizations by emulating high‐resolution physics without the cost of explicit simulation. However, uncertainty about the relationship between offline and online performance (i.e., when integrated with a large‐scale general circulation model) hinders their development. Much of this uncertainty stems from limited sampling of the noisy, emergent effects of upstream ML design decisions on downstream online hybrid simulation. Our work rectifies the sampling issue via the construction of a semi‐automated, end‐to‐end pipeline for size ensembles of hybrid simulations, revealing important nuances in how systematic reductions in offline error manifest in changes to online error and online stability. For example, removing dropout and switching from a Mean Squared Error to a Mean Absolute Error loss both reduce offline error, but they have opposite effects on online error and online stability. Other design decisions, like incorporating memory, converting moisture input from specific humidity to relative humidity, using batch normalization, and training on multiple climates do not come with any such compromises. Finally, we show that ensemble sizes of may be necessary to reliably detect causally relevant differences online. By enabling rapid online experimentation at scale, we can empirically settle debates regarding subgrid ML parameterization design that would have otherwise remained unresolved in the noise.

Lin, Jerry [Department of Earth System Sciences Un↗
DOE OSTI2022

Probabilistic inversion of circular phase spectra: application to two-station phase-velocity dispersion estimation in western Canada

SUMMARY Periodic directional and temporal measurements are common in seismology, and necessitate specific statistical analyses that are appropriate for circular quantities. In this work, we explore the use of a von Mises distribution as a representation of errors on circular seismological observations. Specifically, we automate the estimation of surface-wave phase-velocity dispersion for the teleseismic two-station method, which generally suffers from a 2π phase ambiguity. The use of Bayesian inverse techniques, which aim to rigorously quantify model parameter uncertainty, have become widespread throughout seismology over the last decade. Here, we apply Bayesian inversion to measurements of surface-wave phase spectra in order to estimate 1-D, path-averaged Earth structure between station pairs. The dispersion curve and associated uncertainties are additional results of the inversion, which can then be used as input for subsequent analyses (e.g. tomography). We demonstrate this technique through application to surface-wave recordings from long-running seismic stations throughout western Canada. Our results for over 10 000 station pairs reveal first-order tectonic features consistent with previous studies, which provides confidence in our approach as well as an initial step towards resolving a full 3-D seismic velocity model for the region. This work also presents a foundation for the inversion of surface-wave phase spectra to estimate 3-D Earth structure directly. Finally, the ideas presented in this work are not limited to the inversion of surface-wave phase spectra, but can also be considered for Bayesian geophysical inversion of any circular quantities.

Gosselin, Jeremy M. (ORCID:0000000203754102)↗
DOE OSTI2020

AI for Science: Report on the Department of Energy (DOE) Town Halls on Artificial Intelligence (AI) for Science

The report documents the DOE Town Halls held during 2019 at Argonne National Laboratory, Oak Ridge National Laboratory, Lawrence Berkeley National Laboratory, and in Washington, DC. From July to October 2019, the Argonne, Oak Ridge, and Berkeley National Laboratories hosted a series of four town hall meetings attended by more than 1,000 U.S. scientists and engineers. The goal of the town hall series was to examine scientific opportunities in the areas of artificial intelligence (AI), Big Data, and high-performance computing (HPC) in the next decade, and to capture the big ideas, grand challenges, and next steps to realizing these opportunities. In this report and in the Department of Energy (DOE) laboratory community, we use the term “AI for Science” to broadly represent the next generation of methods and scientific opportunities in computing, including the development and application of AI methods (e.g., machine learning, deep learning, statistical methods, data analytics, automated control, and related areas) to build models from data and to use these models alone or in conjunction with simulation and scalable computing to advance scientific research. The AI for Science town hall discussions focused on capturing the transformational uses of AI that employ HPC and/or data analysis, leveraging data sets from HPC simulations or instruments and user facilities, and addressing scientific challenges unique to DOE user facilities and the agency’s wide-ranging fundamental and applied science enterprise.

36 MATERIALS SCIENCE↗