Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical physics”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

While legume-based rotation influences the chemical composition of mineral-associated organic matter, tillage has little effect on its persistence

Legume-based rotations are key to enhancing soil carbon (C) storage in stable forms such as mineral-associated organic matter (MAOM), improving soil health, and helping maintain long-term C stock. This study investigates how plant inputs and tillage practices influence the distribution, composition, and persistence of MAOM in wheat cropping systems when legumes are incorporated into the rotation. Soil samples were collected from Walla Walla silty loam of northeast Oregon, including conventional tillage (CT) and no-till (NT) treatments in a long-term dryland winter wheat (Triticum aestivum L.) and spring pea (Pisum sativum L.) rotation experiment, with an adjacent grass pasture (GP) plot serving as the baseline. Bulk soil, particulate organic matter (POM), and MAOM fractions were scanned using mid-infrared (Mid-IR) spectroscopy and analyzed for total C and nitrogen (N) as well as radiocarbon age, aggregate stability, and mineralizable C. In WP systems, NT promoted greater C accumulation and stabilization in MAOM than CT, reaching levels similar to those under native GP conditions. Furthermore, our findings indicate that decades of continuous legume-based rotation were insufficient to restore total bulk C content to levels observed in the undisturbed reference site (GP) in the 0–30 cm layer, with increases limited to the surface layer. The Mid-IR data provided insights into how cultivation may alter the chemical composition of MAOM. An increase in amide bands was observed in GP, while legume-based rotation systems favored an abundance of aliphatic (C-H) components in the MAOM fraction. In cultivated soils, C in MAOM exhibited longer residence times under arable conditions than in grassland, while POM remained an actively cycled pool across all treatments. Despite variations in the chemical composition of MAOM across treatments, mineralizable C emissions from MAOM fractions remained statistically unchanged during the 96-hour incubation period. Our findings suggest that soil physical protection, enhanced by reduced tillage, is the primary factor influencing MAOM cycling.

Carbon stabilization↗
DOE OSTI2022

Quantifying Errors in Effective Cluster Interactions of Lattice Gas Cluster Expansions

The promise of lattice gas (LG) cluster expansions (CEs) is that they can describe a given system property to any level of accuracy since the orthogonal “cluster basis functions” span the complete space of available configurations. Such an approach can be constructed to an arbitrarily large surface of a finite number of distinct adsorption sites. Unfortunately, this is only true for the case of an ideal, fixed lattice decorated with components at precise lattice points (the lattice “sites”) with no distortions or relaxations subsequently allowed. Since most systems, and surfaces specifically, do not conform to such an ideal set of constraints, errors in LG CEs must be expected or CE convergence severely hampered. Beyond this, numerical errors in the provided data can complicate the proper construction of a truly predictive and/or physically significant CE. In this work, we show here how reliance on typical statistical tools like confidence intervals cannot be expected to provide an accurate representation of the uncertainty of effective cluster interactions (ECIs) in the CE due to the nature of the target ab initio data and the nature of CEs themselves. We develop a method for estimating these errors that does not rely on statistical assumptions about the model or data. We then use these ECI errors to quantify fundamental consequences on the uncertainty of ECIs in CEs built from O/Fe(100) data whose surface and adsorbates have been allowed to relax in a typical manner and from O/Fe(100) data whose surface and adsorbates are fixed in ideal lattice positions. We also quantify the effect of using a different density functional theory exchange–correlation functional, using these ECI errors to assess the significance in any deviations. In both cases, our method is shown to have remarkable utility in the quantification of errors in the ECIs of CEs. While we stick to the lattice gas convention in this work, the method is equally applicable to the Ising convention or, in principle, any linear model of sufficient complexity.

08 HYDROGEN↗
DOE OSTI2020

Joint Estimation of Behind-the-Meter Solar Generation in a Community

Distribution grid planning, control, and optimization require accurate estimation of solar photovoltaic (PV) generation and electric load in the system. Most of the small residential solar PV systems are installed behind-the-meter making only the net load readings available to the utilities. This paper presents an unsupervised framework for joint disaggregation of the net load readings of a group of customers into the solar PV generation and electric load. Our algorithm synergistically combines a physical PV system performance model for individual solar PV generation estimation with a statistical model for joint load estimation. The electric loads for a group of customers are estimated jointly by a mixed hidden Markov model (MHMM) which enables modeling the general load consumption behavior present in all customers while acknowledging the individual differences. At the same time, the model can capture the change in load patterns over a time period by the hidden Markov states. The proposed algorithm is also capable of estimating the key technical parameters of the solar PV systems. Our proposed method is evaluated using the net load, electric load, and solar PV generation data gathered from residential customers located in Austin, Texas. Testing results show that our proposed method reduces the mean squared error of state-of-the-art net-load disaggregation algorithms by 67%.

behind-the-meter solar generation↗
DOE OSTI2022

Deep Generative Models for Materials Discovery and Machine Learning-Accelerated Innovation

Machine learning and artificial intelligence (AI/ML) methods are beginning to have significant impact in chemistry and condensed matter physics. For example, deep learning methods have demonstrated new capabilities for high-throughput virtual screening, and global optimization approaches for inverse design of materials. Recently, a relatively new branch of AI/ML, deep generative models (GMs), provide additional promise as they encode material structure and/or properties into a latent space, and through exploration and manipulation of the latent space can generate new materials. These approaches learn representations of a material structure and its corresponding chemistry or physics to accelerate materials discovery, which differs from traditional AI/ML methods that use statistical and combinatorial screening of existing materials via distinct structure-property relationships. However, application of GMs to inorganic materials has been notably harder than organic molecules because inorganic structure is often more complex to encode. In this work we review recent innovations that have enabled GMs to accelerate inorganic materials discovery. We focus on different representations of material structure, their impact on inverse design strategies using variational autoencoders or generative adversarial networks, and highlight the potential of these approaches for discovering materials with targeted properties needed for technological innovation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evaluating downscaled products with expected hydroclimatic co-variances

Abstract. There has been widespread adoption of downscaled products amongst practitioners and stakeholders to ascertain risk from climate hazards at the local scale (e.g., ∼ 5 km resolution). Such products must nevertheless be consistent with physical laws to be credible and of value to users. Here we evaluate statistically and dynamically downscaled products by examining local co-evolution of downscaled temperature and precipitation during convective and frontal precipitation events (two mechanisms testable with just temperature and precipitation). We find that two widely used statistical downscaling techniques (Localized Constructed Analogs version 2, LOCA2, and Seasonal Trends and Analysis of Residuals Empirical Statistical Downscaling Model, STAR-ESDM) generally preserve expected co-variances during convective precipitation events over the historical and future projected intervals as compared to European Centre for Medium-Range Weather Forecasts Reanalysis v5 (ERA5) and two observation-based data products (Livneh and nClimGrid-Daily). However, both techniques dampen future intensification of frontal precipitation that is otherwise robustly captured in global climate models (i.e., prior to downscaling) and with process-based dynamical downscaling across five different regional climate models. In the case of LOCA2, this leads to appreciable underestimation of future frontal precipitation event intensity. This study is one of the first to quantify a likely ramification of the stationarity assumption underlying statistical downscaling methods and identify a phenomenon where projections of future change diverge depending on data production method employed. Finally, our work proposes expected co-variances during convective and frontal precipitation as useful evaluation diagnostics that can be universally applied to a wide range of statistically downscaled products.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Tractable Learning in Underexcited Power Grids

Estimating the structure of physical flow networks, such as power grids, is critical to secure delivery of energy. This article discusses statistical structure estimation in power grids in the “underexcited” regime, where a subset of internal nodes has zero injection fluctuations. Prior estimation algorithms based on nodal voltages fail for such grids as the voltage covariance matrix is not invertible. Here, we propose a novel topology learning algorithm for learning underexcited general networks. Our algorithm uses physics-informed conservation laws to first identify the zero-injection buses and their neighbors, and then estimates the remaining edges in the grid. We prove the asymptotic correctness of our algorithm for grids with nonadjacent internal zero-injection nodes. More important, we theoretically analyze our algorithm’s efficacy under noisy measurements, and determine bounds on maximum noise under which asymptotically correct recovery is guaranteed. Our approach is validated through simulations with voltage samples generated on test distribution grids with real injection data and nonlinear power flow models.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Moment extraction using an unfolding protocol without binning

Deconvolving (“unfolding”) detector distortions is a critical step in the comparison of cross-section measurements with theoretical predictions in particle and nuclear physics. However, most existing approaches require histogram binning while many theoretical predictions are at the level of statistical moments. We develop a new approach to directly unfold distribution moments as a function of another observable without having to first discretize the data. Our moment unfolding technique uses machine learning and is inspired by Boltzmann weight factors and generative adversarial networks (GANs). We demonstrate the performance of this approach using jet substructure measurements in collider physics. With this illustrative example, we find that our moment unfolding protocol is more precise than bin-based approaches and is as or more precise than completely unbinned methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Effect of CLIQ on Training of HL-LHC Quadrupole Magnets

The high-luminosity LHC upgrade requires stronger than LHC low-beta quadrupole magnets to reach the luminosity goals of the project. The project is well advanced and HL-LHC quadrupole magnets are currently being commissioned in US Labs (MQXFA magnets) and CERN (MQXFB magnets). Those are the first Nb 3 Sn magnets to be used in any large particle accelerator. At development stages, many Nb 3 Sn accelerator sub-scale models showed relatively slow training and MQXFA magnets were projected to have low tens of quenches before reaching operational field. Recently it was shown that dedicated capacitor-based de-vices can affect Nb 3 Sn magnet training, and it was suggested that CLIQ, a capacitor-based device intended for quench protection, can do too. Here, the present paper investigates effects on training likely induced by CLIQ, using the base fact that only half the coils in a quadrupole experience upward current modulation at quench because of capacitor discharge. The study encompasses all MQXFA production magnets trained at BNL to date. No other high-statistics data from identical magnets (series) with CLIQ protection exist so far. Implications and opportunities stemming from data analysis are discussed and conclusions drawn.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

How data science methods can improve the quality and efficiency of ICF and HEDP research

Data Science methods (many that are Bayesian based) are widely used in the physical sciences to estimate model parameters from experimental data, synthesize heterogeneous data, calibrate models, design experiments, and determine statistical significance of data. These methods provide a wealth of advantages over traditional analysis techniques because: 1) uncertainties are rigorously defined and propagated naturally through complex systems including covariance, 2) prior information is captured within the analysis framework (including rad-MHD and rad-hydro simulations), 3) competing models can be selected and/or ruled out using quantitative criteria, and 4) complex, heterogeneous data can be incorporated simultaneously. While these methods have been widely adopted as the gold standard in fields such as particle physics, astronomy, and biology, they have been slow to catch on in Inertial Confinement Fusion (ICF) and High Energy Density Physics (HEDP) research. Recently, several teams at LLNL, SNL, LANL, and the LLE have been exploring the use of these tools in their research and have found success. Here we propose that a concerted effort to consolidate these independent research efforts by developing and deploying common tools for use across the complex can revolutionize the way we approach data analysis, assimilation of theory and experiment, and decision making. The Bayesian formalism provides a means to accomplish this, but we are lacking certain infrastructure to make it happen on a large scale. Furthermore, once adopted, these techniques can be used to develop standards by which discoveries can be judged, similar to the so-called 5σ rule in high energy particle physics. Such standards may be used in the future to address the issue of unknown reproducibility in ICF and HED experiments caused by low shot rate and high cost per experiment. Our goals as a group are to advance the state of the art in HED measurement science by enabling: 1) better inferences from data with well-defined uncertainties, 2) better use of the data we have and continue to collect, 3) intelligent synthesis of data, 4) evaluation of the statistical significance of our data, and 5) informed decision making regarding the design of new experiments and instruments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Energy dependence of light hypernuclei production in heavy-ion collisions from a coalescence and statistical-thermal model perspective

A comparison of light hypernuclei production, from Ultra-relativistic Quantum Molecular Dynamics (UrQMD)+coalescence and the thermal model, in heavy ion collisions over a wide range of beam energies and system sizes is presented. We find that both approaches provide generally similar results, with differences in specific details. Especially the ratios of hypertriton to Λ are affected by both the source radius Δr of the coalescence procedure as well as canonical effects. On the other hand, the double ratio S 3 is almost independent of canonical effects, which is in contrast to coalescence. Thus, both the beam energy dependence and centrality dependence of S 3 can be used to constrain the hypertriton source radius. To do so the currently available data are not yet sufficient. Elliptic flow is shown to be unaffected by the source size of the nuclei and an almost perfect mass scaling of the elliptic flow is observed. Furthermore, our predictions further suggest that the existence of the H-dibaryon(ΛΛ) seems ruled out by ALICE data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Optimizing time integration for accurate recovery of shockwave interface location in radiography

We present simulations and experiments of time integrated radiographic imaging of a moving 1D shock wave front and a quantitative method for determining the statistical error in locating the shock front as a function of integration time and noise in the radiograph. We discuss the trade-off between increasing motion blur, which leads to decreased shock front location certainty, and increasing signal-to-noise, which leads to improved image quality with increasing integration time. We find an optimum integration time between a short integration time, where noise limits the error, and a long integration time, where motion blurring limits the error. This methodology can be used to tune experimental configurations to obtain the highest quality radiograph for a given experimental configuration.

Bremsstrahlung↗
DOE OSTI2022

Molecular-based analysis of nanoparticle solvation: Classical density functional approach

Proper statistical mechanics understanding of nanoparticle solvation processes requires an accurate description of the molecular structure of the solvent. Achieving this goal with standard molecular dynamics (MD) simulation methods is challenging due to large length scales. An alternative approach to this problem can be formulated using classical density functional theory (cDFT), where a full configurational description of the positions of all the atoms is replaced by collective atomic site densities in the molecule. Using an example of the negatively charged silica-like system in an aqueous polar environment represented by a two-site water model, we demonstrate here that cDFT can reproduce MD data at a fraction of the computational cost. An important implication of this result is the ability to understand how the solvent molecular features may affect the system’s properties at the macroscopic scale. A concrete example highlighted in this work is the analysis of nanoparticle interactions with sizes of up to 100 nm in diameter.

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

Preface to special topic: The High Repetition Rate Frontier in High-Energy-Density Physics

High-repetition-rate (HRR) experiments can collect large datasets with high temporal, spatial, and/or parametric resolution or large numbers of repeat measurements for statistics. HRR experiments also enable new experimental designs, including active feedback control loops and novel diagnostics, that can improve the reproducibility as well as the quantity of measurements. Together, these attributes make HRR experiments ideal for performing high-quality repeatable science. Until recently, these techniques have not been applied to high-energy-density–physics (HEDP) experiments, which are typically restricted to repetition rates of a few per day. However, recent advancements in lasers, pulsed-power drivers, target fabrication, and diagnostics are starting to change this fact, opening an exciting new frontier of HRR HEDP experiments. A mini-conference on this subject at the 2021 meeting of the American Physical Society Division of Plasma Physics brought together members of this growing community. As a result, the “High Repetition Rate Frontier in High-Energy-Density Physics” special topic in Physics of Plasmas highlights current progress in this exciting area.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2024

A new Monte Carlo generator for BSM physics in B → K*ℓ+ℓ− decays with an application to lepton non-universality in angular distributions

Abstract Within the widely used EvtGen framework, we have added a new event generator model forB → K * ℓ + ℓ − with improved standard model (SM) decay amplitudes and possible BSM physics contributions, which are implemented in the operator product expansion in terms of Wilson coefficients. This event generator can then be used to estimate the statistical sensitivity of a simulated experiment to the most general BSM signal resulting from dimension-six operators. We describe the advantages and potential of the newly developed ‘Sibidanov Physics Generator’ in improving the experimental sensitivity of searches for lepton non-universal BSM physics and clarifying signatures. The new generator can properly simulate BSM scenarios, interference between SM and BSM amplitudes, and correlations between different BSM observables as well as acceptance bias. We show that exploiting such correlations substantially improves experimental sensitivity. As a demonstration of the utility of the MC generator, we examine the prospects for improved measurements of lepton non-universality in angular distributions forB→K * ℓ + ℓ − decays from the expected 50 ab −1 data set of the Belle II experiment, using a four-dimensional unbinned maximum likelihood fit. We describe promising experimental signatures and correlations between observables. The use of lepton-universality violating ∆-observables significantly reduces uncertainties in the SM expectations due to QCD and resonance effects and is ideally suited for Belle II with the large data sets expected in the next decade. Thanks to the clean experimental environment of ane + e − machine, Belle II should be able to probe BSM physics in the Wilson coefficientsC 7 and$$ {C}_7^{\prime } $$ C 7 ′ , which appear at lowq 2 in the di-electron channel.

Physics↗
DOE OSTI2021

Model-independent extraction of the proton charge radius from PRad data

The proton radius puzzle has motivated several new experiments that aim to extract the proton charge radius and resolve the puzzle. Recently, PRad, a new electron–proton scattering experiment at Jefferson Lab, reported a proton charge radius of [Formula: see text]. The value was obtained by using a rational function model for the proton electric form factor. We perform a model-independent extraction using [Formula: see text]-expansion of the proton charge radius from PRad data. We find that the model-independent statistical error is more than 50% larger compared to the statistical error reported by PRad.

Astronomy & Astrophysics↗
DOE OSTI2022

COG Beta-Effective Benchmarks

COG 11.3 calculations with the “new method” using ENDF/B-VIII.0 data are shown to produce results with |C/E – 1| < 3σ in all cases, whereas the other methods produce a few minor outliers. Furthermore, the “new method” using ENDF/B-VIII.0 data is observed to produce “over” and “under” results consistent with 50/50 expectations whereas the other methods appear to be slightly over or under predicting. JEFF-3.3 results are slightly larger overall than ENDF/B-VIII.0 results. In all but a few cases, it is also observed that $F_d/F_t ≥ (vF)_d/(vF)_t ≥ 1 – (k_{prompt}/k_{eff})$ and all estimators produce equivalent results for the most thermal systems. In those other cases, this relation holds within the given statistics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Progress Towards an Investigation of the MiniBooNE Low Energy Excess Using Neutral-Current Delta-Like Single Photons in MicroBooNE with Wire-Cell 3D Reconstruction Algorithms

This note presents progress towards a Neutral-Current (NC) Delta-like single photon search in MicroBooNE using Wire-Cell 3D reconstruction and pattern recognition algorithms. This analysis will help address the question of the MiniBooNE Low Energy Excess (LEE) [1]. The selections are similar to the previously published Pandora-based single photon selections [2] in many ways, but here we use independent reconstruction and selection tools and examine a larger phase space (including charged pions and multiple protons). In the future, we will significantly reduce systematic uncertainties via the use of constraining sideband observations in order to statistically test specific models of the MiniBooNE LEE.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗