Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical model”

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

Lagrangian and Eulerian Supersaturation Statistics in Turbulent Cloudy Rayleigh–Bénard Convection: Applications for LES Subgrid Modeling

Abstract Turbulent fluctuations of scalar and velocity fields are critical for cloud microphysical processes, e.g., droplet activation and size distribution evolution, and can therefore influence cloud radiative forcing and precipitation formation. Lagrangian and Eulerian water vapor, temperature, and supersaturation statistics are investigated in direct numerical simulations (DNS) of turbulent Rayleigh–Bénard convection in the Pi Convection Cloud Chamber to provide a foundation for parameterizing subgrid-scale fluctuations in atmospheric models. A subgrid model for water vapor and temperature variances and covariance and supersaturation variance is proposed, valid for both clear and cloudy conditions. Evaluation of phase change contributions through an a priori test using DNS data shows good performance of the model. Supersaturation is a nonlinear function of temperature and water vapor, and relative external fluxes of water vapor and heat (e.g., during entrainment-mixing and phase change) influence turbulent supersaturation fluctuations. Although supersaturation has autocorrelation and structure functions similar to the independent scalars (temperature and water vapor), the autocorrelation time scale of supersaturation differs. Relative scalar fluxes in DNS without cloud make supersaturation PDFs less skewed than the adiabatic case, where they are highly negatively skewed. However, droplet condensation changes the PDF shape response: it becomes positively skewed for the adiabatic case and negatively skewed when the sidewall relative fluxes are large. Condensation also increases correlations between water vapor and temperature in the presence of relative scalar fluxes but decreases correlations for the adiabatic case. These changes in correlation suppress supersaturation variability for the nonadiabatic cases and increase it for the adiabatic case. Implications of this work for subgrid microphysics modeling using a Lagrangian stochastic scheme are also discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Maximum Entropy Principle in Deep Thermalization and in Hilbert-Space Ergodicity

We report universal statistical properties displayed by ensembles of pure states that naturally emerge in quantum many-body systems. Specifically, two classes of state ensembles are considered: those formed by (i) the temporal trajectory of a quantum state under unitary evolution or (ii) the quantum states of small subsystems obtained by partial, local projective measurements performed on their complements. These cases, respectively, exemplify the phenomena of “Hilbert-space ergodicity” and “deep thermalization.” In both cases, the resultant ensembles are defined by a simple principle: The distributions of pure states have maximum entropy, subject to constraints such as energy conservation, and effective constraints imposed by thermalization. We present and numerically verify quantifiable signatures of this principle by deriving explicit formulas for all statistical moments of the ensembles, proving the necessary and sufficient conditions for such universality under widely accepted assumptions, and describing their measurable consequences in experiments. We further discuss information-theoretic implications of the universality: Our ensembles have maximal information content while being maximally difficult to interrogate, establishing that generic quantum state ensembles that occur in nature hide (scramble) information as strongly as possible. Our results generalize the notions of Hilbert-space ergodicity to time-independent Hamiltonian dynamics and deep thermalization from infinite to finite effective temperature. Our work presents new perspectives to characterize and understand universal behaviors of quantum dynamics using statistical and information-theoretic tools.

Eigenstate thermalization↗
DOE OSTI2022

Freely jointed chain models with extensible links

We report analytical relations for the mechanical response of single polymer chains are valuable for modeling purposes, on both the molecular and the continuum scale. These relations can be obtained using statistical thermodynamics and an idealized single-chain model, such as the freely jointed chain model. To include bond stretching, the rigid links in the freely jointed chain model can be made extensible, but this almost always renders the model analytically intractable. Here, an asymptotically correct statistical thermodynamic theory is used to develop analytic approximations for the single-chain mechanical response of this model. The accuracy of these approximations is demonstrated using several link potential energy functions. This approach can be applied to other single-chain models, and to molecular stretching in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Data Assimilation for Robust UQ Within Agent-Based Simulation on HPC Systems

Agent-based simulation provides a powerful tool for in silico system modeling. However, these simulations do not provide built-in methods for uncertainty quantification (UQ). Within these types of models a typical approach to UQ is to run multiple realizations of the model then compute aggregate statistics. This approach is limited due to the compute time required for a solution. When faced with an emerging biothreat, public health decisions need to be made quickly and solutions for integrating near real-time data with analytic tools are needed. We propose an integrated Bayesian UQ framework for agent-based models based on sequential Monte Carlo sampling. Given streaming or static data about the evolution of an emerging pathogen this Bayesian framework provides a distribution over the parameters governing the spread of a disease through a population. These estimates of the spread of a disease may be provided to public health agencies seeking to abate the spread. By coupling agent-based simulations with Bayesian modeling in a data assimilation, our proposed framework provides a powerful tool for modeling dynamical systems in silico. We propose a method which reduces model error and provides a range of realistic possible outcomes. Moreover, our method addresses two primary limitations of ABMs: the lack of UQ and an inability to assimilate data. Our proposed framework combines the flexibility of an agent-based model with UQ provided by the Bayesian paradigm in a workflow which scales well to HPC systems. We provide algorithmic details and results on a simulated outbreak with both static and streaming data.

Spannaus, Adam [ORNL] (ORCID:0000000225213657)↗
DOE OSTI2020

Reconstruction of effective potential from statistical analysis of dynamic trajectories

The broad incorporation of microscopic methods is yielding a wealth of information on the atomic and mesoscale dynamics of individual atoms, molecules, and particles on surfaces and in open volumes. Analysis of such data necessitates statistical frameworks to convert observed dynamic behaviors to effective properties of materials. Here, we develop a method for the stochastic reconstruction of effective local potentials solely from observed structural data collected from molecular dynamics simulations (i.e., data analogous to those obtained via atomically resolved microscopies). Using the silicon vacancy defect in graphene as a model, we apply the statistical framework presented herein to reconstruct the free energy landscape from the calculated atomic displacements. Evidence of consistency between the reconstructed local potential and the trajectory data from which it was produced is presented, along with a quantitative assessment of the uncertainty in the inferred parameters.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Mana and thermalization: Probing the feasibility of near-Clifford Hamiltonian simulation

Quantum hydrodynamics is the emergent classical dynamics governing transport of conserved quantities in generic strongly interacting quantum systems. Here, recent matrix product operator methods [1,2] have made simulations of quantum hydrodynamics in 1+1D tractable, but they do not naturally generalize to 2+1D or higher, and they offer limited guidance as to the difficulty of simulations on quantum computers. Near-Clifford simulation algorithms are not limited to one dimension, and future error-corrected quantum computers will likely be bottlenecked by non-Clifford operations. We therefore investigate the non-Clifford resource requirements for simulation of quantum hydrodynamics using mana, a resource theory of non-Clifford operations. For infinite-temperature starting states, we find that the mana of subsystems quickly approaches zero, while for starting states with energy above some threshold the mana approaches a nonzero value. Surprisingly, in each case the finite-time mana is governed by the subsystem entropy, not the thermal state mana; we argue that this is because mana is a sensitive diagnostic of finite-time deviations from canonical typicality.

1-dimensional spin chains↗
DOE OSTI2022

Rapid attribution analysis of the extraordinary heat wave on the Pacific coast of the US and Canada in June 2021

Towards the end of June 2021, temperature records were broken by several degrees Celsius in several cities in the Pacific Northwest areas of the US and Canada, leading to spikes in sudden deaths and sharp increases in emergency calls and hospital visits for heat-related illnesses. Here we present a multi-model, multi-method attribution analysis to investigate the extent to which human-induced climate change has influenced the probability and intensity of extreme heat waves in this region. Based on observations, modelling and a classical statistical approach, the occurrence of a heat wave defined as the maximum daily temperature (TXx) observed in the area 45–52 °N, 119–123 °W, was found to be virtually impossible without human-caused climate change. The observed temperatures were so extreme that they lay far outside the range of historical temperature observations. This makes it hard to state with confidence how rare the event was. Using a statistical analysis that assumes that the heat wave is part of the same distribution as previous heat waves in this region led to a first-order estimation of the event frequency of the order of once in 1000 years under current climate conditions. Using this assumption and combining the results from the analysis of climate models and weather observations, we found that such a heat wave event would be at least 150 times less common without human-induced climate change. Also, this heat wave was about 2 °C hotter than a 1-in-1000-year heat wave would have been in 1850–1900, when global mean temperatures were 1.2 °C cooler than today. Looking into the future, in a world with 2 °C of global warming (0.8 °C warmer than today), a 1000-year event would be another degree hotter. Our results provide a strong warning: our rapidly warming climate is bringing us into uncharted territory with significant consequences for health, well-being and livelihoods. Adaptation and mitigation are urgently needed to prepare societies for a very different future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Machine-Learning Assisted Identification of Battery Life Models

Predictive battery life models are commonly utilized to extrapolate degradation trends observed during accelerated aging tests for simulation of degradation in real-world applications. Thus, fitting accelerated aging data as accurately as possible and with low uncertainty is crucial for making believable projections of battery lifetime, but it is challenging to identify algebraic expressions that accurately fit multivariate degradation trends. A review of models published in literature reveal some common expressions for fitting calendar aging data, which is only dependent on temperature and state-of-charge, but no consistency across many models for fitting cycle aging data, indicating the need for a statistically rigorous data driven approach for developing empirical models. This talk will describe a machine-learning assisted method for identification of predictive battery life models utilizing bilevel optimization and symbolic regression. Bilevel optimization with cross-validation is used to statistically determine cell- and stress-dependent model parameters, while symbolic regression identifies both linear and multiplicative candidate expressions to predict stress-dependent degradation rates by selecting low-order subsets of features from a generated feature library. Because model expressions are identified empirically, it is crucial to ensure resulting models behave according to physical expectations, so the stability of models for interpolation or extrapolation is interrogated qualitatively through simulation and quantitatively through cross-validation and uncertainty quantification via bootstrap resampling. This model identification approach substantially improves upon models identified purely using expert judgement in terms of both accuracy and uncertainty. Model simulation and validation is then conducted by deriving a state-equation form of the predictive model, enabling simulation of battery aging under dynamic stresses. This enables validation of the predictive battery model on lab-based tests with varying conditions or on drive-cycle or application-cycle testing protocols. Parameter uncertainty can be carried forward into model simulation, giving lifetime estimates and confidence windows for cell- or system-level lifetime. The financial impact of battery model uncertainty can be estimated by incorporating uncertainty into a technoeconomic model.

battery↗
DOE OSTI2022

Ensemble approximate control variate estimators: Applications to multi-fidelity importance sampling.

The recent growth in multifidelity uncertainty quantification has given rise to a large set of variance reduction techniques that leverage information from model ensembles to provide variance reduction for estimates of the statistics of a high-fidelity model. In this paper we provide two contributions: (1) we utilize an ensemble estimator to account for uncertainties in the optimal weights of approximate control variate (ACV) approaches and derive lower bounds on the number of samples required to guarantee variance reduction; and (2) we extend an existing multifidelity importance sampling (MFIS) scheme to leverage control variates. Our approach directly addresses a limitation of many multifidelity sampling strategies that require the usage of pilot samples to estimate covariances. As such we make significant progress towards both increasing the practicality of approximate control variates—for instance, by accounting for the effect of pilot samples—and using multifidelity approaches more effectively for estimating low-probability events. The numerical results indicate our hybrid MFIS-ACV estimator achieves up to 50% improvement in variance reduction over the existing state-of-the-art MFIS estimator, which had already shown an outstanding convergence rate compared to the Monte Carlo method, on several problems of computational mechanics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Modeling of H2 Dispersion at ARIES

Hydrogen is a versatile and clean energy carrier that can be produced from various renewable sources such as wind, solar, and hydropower and help decarbonize electricity grids, industry, and transportation. Using the Hydrogen Research Facility under Advanced Research on Integrated Energy Systems (ARIES) at the National Renewable Energy Laboratory's (NREL) Flatirons campus as a test bench, the study examines the feasibility, useability, and value of using computational fluid dynamics (CFD) techniques to model hydrogen dispersion. The ARIES facility was chosen because controlled hydrogen releases can be performed at a rate of 27 kg-H2/hr. Site-specific atmospheric and weather condition data such as wind speed and temperature were used as inputs to the model. The results show statistical distributions and ranges of hydrogen concentrations at locations throughout the domain. Wind conditions are found to significantly impact the release behavior, including the hydrogen cloud's direction and concentrations. At low wind speeds (below 1 mph), hydrogen forms a cloud and at higher wind speeds (> 2-4 mph) hydrogen plume stretches in the direction of wind momentum. From >100 simulations for ARIES site-specific conditions, statistical quantities combined with a clustering algorithm were used to propose sensor location at various elevations from ground.

dispersion↗
DOE OSTI2024

Time Dilated Bundt Cake Analysis of PV Output [Poster]

We present a novel method for modeling time-dependent statistics in the power signal generated by a photovoltaic (PV) system. Our white-box machine learning method is interpretable and auditable, based on principles of multiperiodic basis functions and convex optimization. Our proposed method of time dilating the daily signal to remove night time values results in a novel representation of PV power signals, evocative of a ‘Bundt cake’. The proposed model describes the marginal distribution of power output as a function of date and time. The resulting probabilistic model of a PV system can be used to perform a variety of tasks, and here, we demonstrate the application of clear sky detection.

14 SOLAR ENERGY↗
DOE OSTI2026

First search for sterile neutrino oscillation leading to $ν_μ$ disappearance in the Booster Neutrino Beam at ICARUS

We present a search for muon neutrino disappearance in the Booster Neutrino Beam (BNB) at Fermilab using the ICARUS detector. Neutrino interactions identified as muon neutrinos interacting with argon nuclei via the charged current interaction and having only a muon and at least one proton in the final state (1$μ$Np) have been selected from data collected in 2022-2023 (ICARUS Run 2) and compared with a simulation-based expectation. In the context of a fit to a two-neutrino approximation of the sterile 3+1 model, including the impact of systematic uncertainty from the flux, neutrino interaction, and detector models, we find no statistically significant muon neutrino disappearance at the ICARUS baseline of 600 meters from the BNB target. Corresponding 90% C.L. exclusion contours in $Δm^2_{41}$ - sin$^22θ_{μμ}$ space are presented. This is the first oscillation analysis produced by ICARUS exposed to the BNB. We note that the analysis is systematics limited due to large unconstrained uncertainties from the flux and interaction models. In future joint analyses, data from ICARUS and the SBND detector, exposed to the BNB at 110 meters from target, will be combined to provide significant constraint of these uncertainties, enabling a robust, world-leading two-detector analysis.

Alrahman, F. Abd [Houston U.]↗
DOE OSTI2023

Soft matter roadmap

Soft materials are usually defined as materials made of mesoscopic entities, often self-organised, sensitive to thermal fluctuations and to weak perturbations. Archetypal examples are colloids, polymers, amphiphiles, liquid crystals, foams. The importance of soft materials in everyday commodity products, as well as in technological applications, is enormous, and controlling or improving their properties is the focus of many efforts. From a fundamental perspective, the possibility of manipulating soft material properties, by tuning interactions between constituents and by applying external perturbations, gives rise to an almost unlimited variety in physical properties. Together with the relative ease to observe and characterise them, this renders soft matter systems powerful model systems to investigate statistical physics phenomena, many of them relevant as well to hard condensed matter systems. Understanding the emerging properties from mesoscale constituents still poses enormous challenges, which have stimulated a wealth of new experimental approaches, including the synthesis of new systems with, e.g. tailored self-assembling properties, or novel experimental techniques in imaging, scattering or rheology. Theoretical and numerical methods, and coarse-grained models, have become central to predict physical properties of soft materials, while computational approaches that also use machine learning tools are playing a progressively major role in many investigations. This Roadmap intends to give a broad overview of recent and possible future activities in the field of soft materials, with experts covering various developments and challenges in material synthesis and characterisation, instrumental, simulation and theoretical methods as well as general concepts.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Calibration and Commissioning of the LSST

Understanding the nature of dark energy and dark matter remains one of the fundamental questions in physics today; impacting our understanding of particle physics, cosmology, and possibly theories of gravity. Given the scale and complexity of the next generation of cosmology experiments (e.g., the Rubin Observatory, the Euclid satellite mission, and the Roman space telescope) we are entering an era where statistical noise no longer determines the accuracy to which we can measure cosmological parameters. Our ability to control and correct for systematics will ultimately determine the scientific impact of these experiments. This award addressed the challenge of how we determine what limits the accuracy of our cosmological measures, what techniques are appropriate for measuring and calibrating the properties of galaxies to best constrain cosmological models, how to develop statistical techniques that are insensitive to systematic errors, and how to optimize survey strategies in order to minimize systematics while maximizing the speed at which an experiment can achieve its science objectives. In this final technical report for award DE-SC0011635 we describe a set of open-source frameworks that simulate the characteristics and properties of current and planned cosmology surveys and the application of these frameworks to the development of new methodologies for estimating the properties and distances to galaxies that are robust to noisy and incomplete data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Measurement of the total and differential Higgs boson production cross-sections at $ \sqrt{s} $ = 13 TeV with the ATLAS detector by combining the H → ZZ* → 4ℓ and H → γγ decay channels

The total and differential Higgs boson production cross-sections are measured through a combined statistical analysis of the H → ZZ* → 4ℓ and H → γγ decay channels. The results are based on a dataset of 139 fb –1 of proton–proton collisions at a centre-of-mass energy of 13 TeV, recorded by the ATLAS detector at the Large Hadron Collider. The measured total Higgs boson production cross-section is $55.5$$^{+4.0}_{–3.8}$ pb, consistent with the Standard Model prediction of 55.6 ± 2.5 pb. All results from the two decay channels are compatible with each other, and their combination agrees with the Standard Model predictions. A combined statistical interpretation of the measured fiducial cross-sections as a function of the Higgs boson transverse momentum is performed in order to probe the Yukawa couplings to the bottom and charm quarks. A similar interpretation is performed by including also the constraints from the measurements of Higgs boson production in association with a W or Z boson in the H → bb¯ and cc¯ decay channels.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Separate the Role of Southern and Northern Extra‐Tropical Pacific in Tropical Pacific Climate Variability

Abstract Observational and modeling studies have elucidated the influential role played by the southern and northern extratropical Pacific (SEP and NEP) forcing in shaping dynamics of tropical Pacific climate variability. However, the relative importance of the NEP and SEP and the timescale on which they impact the tropics remain unclear. Using a linear inverse model (LIM) that selectively incorporates or excludes tropical‐extratropical coupling, we find a reduction in tropical interannual variability (∼40%) and low‐frequency (sub‐decadal to decadal) variability in the southeastern tropical Pacific region (∼70%) in the absence of SEP. Conversely, the absence of NEP yields no significant impact on tropical interannual variability but markedly diminishes low‐frequency variability in the central tropical Pacific region (∼70%). LIM and statistic diagnostics on CMIP6 models show the low‐frequency to total variability ratio in the tropical Pacific depending on their NEP and SEP representation. Models with more (less) low‐frequency power tend to show stronger NEP (SEP) dynamics.

Geology↗
DOE OSTI2021

A Comprehensive Open-Source R Software For Statistical Metrology Calculations: From Uncertainty Evaluation To Risk Analysis

Whether calibrating equipment or inspecting products on the factory floor, metrology requires many complicated statistical calculations to achieve a full understanding and evaluation of measurement uncertainty and quality. In order to assist its workforce in performing these calculations in a consistent and rigorous way, the Primary Standards Lab at Sandia National Laboratories (SNL) has developed a free and open-source software package for computing various metrology calculations from uncertainty propagation to risk analysis. In addition to propagating uncertainty through a measurement model using the well-known Guide to Expression of Uncertainty in Measurement or Monte Carlo approaches, evaluating the individual Type A and Type B uncertainty components that go into the measurement model often requires other statistical methods such as analysis of variance or determining uncertainty in a fitted curve. Once the uncertainty in a measurement has been calculated, it is usually evaluated from a risk perspective to ensure the measurement is suitable for making a particular conformance decision. Finally, SNL’s software can perform all these calculations in a single application via an easy-to-use graphical interface, where the different functions are integrated so the results of one calculation can be used as inputs to another calculation.

97 MATHEMATICS AND COMPUTING↗