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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The role of ProtoDUNE-SP in future oscillation physics

Long-baseline oscillation experiments are entering a new phase where neutrino beams can provide ultra-high fuxes over long distances. These beams will lead to an era of high statistics in large neutrino detectors, creating new systematic uncertainty-limited analyses. The Deep Underground Neutrino Experiment (DUNE) plans a high fux beam to travel 1300 km to 10-kT liquid argon detector modules. While neither the near detector nor the far detector has been constructed, prototypes of the DUNE Far Detector have been built and operated at CERN. One prototype, ProtoDUNE Single-Phase, is a full-scale mockup of the first 10-kT DUNE Far Detector module and can inform DUNE on detector performance and particle passage in argon. The proceedings discuss some published and in-progress results from ProtoDUNE-SP and how these results can impact current sensitivity studies and future analyses.

Diurba, Richard↗
DOE OSTI2020

Statistical Tests for Diagnosing Fission Source Convergence and Undersampling in Monte Carlo Criticality Calculations [Slides]

There is a very strong need for statistical testing to determine fission source convergence in Monte Carlo criticality calculations. Automation of such tests will greatly streamline and support the work carried out by NCS practitioners. Recent R&D work has shown that no single statistical test for convergence is sufficiently reliable, robust, and “guaranteed.” However, a combination of several standard statistical tests for the similarity of distributions, coupled with a high-fidelity estimate of the fission-matrix source is sufficiently robust, reliable, and repeatable that convergence can be “guaranteed.” During the course of the EG-AMCT studies, a number of statistical metrics and tests were proposed for diagnosing clustering and undersampling. None of these was robust and reliable enough for practical use in production codes. However, the expert group efforts came close. Some recent R&D work stemming from those past efforts has been very successful and promising. This Sub-Group will provide international input and collaboration on the development and implementation of statistical tests for convergence, with the primary goal of having the MC codes automatically detect convergence (or lack thereof). Newly proposed statistical tests to detect undersampling (after convergence) will also be reviewed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

(U) Correlated Sampling Using Batch Statistics to Reduce the Uncertainty of Combinations of KSEN Outputs with MCNP6

The relative sensitivity of k eff to the densities of nuclides in a material are combined to compute relative sensitivities to other inputs. When computed in a single Monte Carlo run, the nuclide density sensitivities are correlated, and the statistical uncertainties propagated to other inputs will be incorrect unless those correlations are accounted for. Equations are presented to apply correlated sampling using batch statistics for the sum of an arbitrary number of random tallies and the difference of two random tallies when each is multiplied by a different constant. When correlated sampling is used on a recent benchmark evaluation, the correct statistical uncertainties for certain combinations of sensitivities are dramatically smaller than the incorrect uncertainties. A user-controlled, regular output of MCNP6’s KSEN sensitivities that allows batch statistics to be applied to combinations would be extremely valuable.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Quantifying uncertainty in physics-based predictions of rare-isotope production cross sections via Bayesian-inspired model averaging across nuclear mass tables

Accurate prediction of fragmentation cross sections is essential for rare-isotope beam production, planning new-isotope searches, and designing experiments to study the most exotic regions of the nuclear chart. However, existing reaction models and phenomenological cross-section parametrizations often exhibit significant deviations over broad regions of mass and charge. In this work, a Bayesian-inspired model-averaging framework is developed to combine abrasion-ablation (AA) calculations based on multiple nuclear mass tables into a single statistically weighted estimate. For the calibrated systems, the model weights are assigned empirically according to the relative quality of fit to measured cross sections, thereby reducing systematic model bias while preserving the underlying physics content of the AA description. The weights are constrained using proton-rich fragmentation data for the 78 Kr and 124 Xe projectiles. The resulting parameter trends are then propagated to the 92 Mo and 144 Sm systems through a controlled scaling procedure. In the present implementation, the excitation-energy prescription is fixed, while the averaging is performed across nuclear-mass inputs; the framework provides both weighted cross sections and associated uncertainty estimates. Applied to proton-rich fragmentation, the present approach provides a practical basis for interpolation and limited extrapolation in regions relevant to rare-isotope production. The resulting predictions are used to assess the production of very proton-rich nuclei, and candidate new isotopes are discussed.

Bayesian methods↗
DOE OSTI2025

Origin of the ring ellipticity in the black hole images of M87*

We investigate the origin of the elliptical ring structure observed in the images of the supermassive black hole M87*, aiming to disentangle contributions from gravitational, astrophysical, and imaging effects. Leveraging the enhanced capabilities of the Event Horizon Telescope (EHT)'s 2018 array, including improved (u,v)-coverage from the Greenland Telescope, we measured the ring's ellipticity using five independent imaging methods, obtaining a consistent average value of τ = 0.08−0.02+0.03 with a position angle of ξ = 50.1−7.6+6.2 degrees. To interpret this measurement, we compared it to general relativistic magnetohydrodynamic (GRMHD) simulations spanning a wide range of physical parameters including the thermal or nonthermal electron distribution function, spins, and ion-to-electron temperature ratios in both low- and high-density regions. We find no statistically significant correlation between spin and ellipticity in GRMHD images. Instead, we identify a correlation between ellipticity and the fraction of non-ring emission, particularly in nonthermal models and models with higher jet emission. These results indicate that the ellipticity measured from the M87* emission structure is consistent with that expected from simulations of turbulent accretion flows around black holes, where it is dominated by astrophysical effects rather than gravitational ones. Future high-resolution imaging, including space very long baseline interferometry and long-term monitoring, will be essential to isolate gravitational signatures from astrophysical effects.Key words: accretion, accretion disks / black hole physics / gravitation / galaxies: active** These authors have contributed equally to this work.*** NASA Hubble Fellowship Program, Einstein Fellow.† Deceased.

Dahale, Rohan [Granada U., Theor. Phys. Astrophys.↗
DOE OSTI2023

Hierarchical reconstruction of 3D well-connected porous media from 2D exemplars using statistics-informed neural network

The relationships between porous microstructures and transport properties are of fundamental importance in various scientific and engineering applications. Due to the intricacy, stochasticity and heterogeneity of porous media, reliable characterization and modeling of transport properties often require a complete dataset of internal microstructure samples. However, it is often an unbearable cost to acquire sufficient 3D digital microstructures by purely using microscopic imaging systems. Herein this paper presents a machine learning-based technique to hierarchically reconstruct 3D well-connected porous microstructures from one isotropic or several anisotropic low-cost 2D exemplar(s). To compactly characterize the large-scale microstructural features, a Gaussian image pyramid is built for each 2D exemplar. Local morphology patterns are collected from the Gaussian image pyramids, and then they serve as the training data to embed the 2D morphological statistics into feed-forward neural networks at multiple length levels. By using a specially-developed morphology integration scheme, the 3D morphological statistics at different levels can be inferred from the statistics-informed neural networks. Gibbs sampling is adopted to hierarchically reconstruct 3D microstructures by using multi-level 3D morphological statistics, where the large-scale, regional and local morphological patterns are statistically generated and successively added to the same 3D random field. The proposed method is tested on a series of porous media with distinct morphologies, and the statistical equivalence between the reconstructed and the real microstructures is systematically evaluated by comparing morphological descriptors and transport properties. The results demonstrate that the proposed 2D-to-3D microstructure reconstruction method is a universal and efficient approach to generating morphologically and physically realistic samples of porous media.

42 ENGINEERING↗
DOE OSTI2023

zPerf: A Statistical Gray-Box Approach to Performance Modeling and Extrapolation for Scientific Lossy Compression

With the scaling up of simulation-based scientific discovery on high-performance computing systems, the disparity between compute and I/O has increased, forcing domain scientists to save only a small amount of simulation data to persistent storage. This can result in the loss of essential physics fields that are needed for data analysis. While error-bounded lossy compression has made tremendous progress in bridging the gap between compute and I/O, the lack of understanding of compression performance remains a key hurdle to its wide adoption. Here, in this work, we present zPerf, a statistical gray-box performance modeling approach for scientific lossy compression. Our contributions are threefold: 1) We develop zPerf to estimate the performance of lossy compression techniques, based on in-depth understanding and statistical modeling for data features and core compression metrics; 2) We demonstrate the in-detailed implementation of zPerf using two case studies, where we derive the performance modeling for SZ and ZFP, two leading lossy compressors; 3) We evaluate the effectiveness of zPerf on real-world datasets across various domains. Based on the evaluation, we demonstrate the efficacy of the zPerf performance model; 4) We further discuss three case studies where zPerf is applied to extrapolate the compression ratio of SZ and ZFP with alternative encoding schemes as well as ZFP with an alternative transform scheme. Through the case studies, we demonstrate the potential of zPerf for exploring the design space of lossy compression, which has hardly been studied in the literature.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Hypertriton puzzle in relativistic heavy-ion collisions

The yields of hadrons and light nuclei in relativistic collisions of heavy-nuclei at a center of mass energy of $\sqrt{s_{NN}}$ = 2.6 TeV can be described remarkably well by a thermal distribution of an ideal gas of hadrons and light nuclei interacting only via the decay of resonances. Given the particularly small binding energy of hypertritons relative to the temperature describing the yields (about 156 MeV), one might naturally expect hypertrions to dissociate in medium, making the agreement of hypertriton yields with thermal predictions highly puzzling. The puzzle is compounded by the fact that small binding energy is associated with the large size of the hypertriton. This size is on a similar scale to the overall size of the fireball and much larger than the length scale over which temperatures in the fireball vary over phenomenologically relevant amounts. Here, this paper quantifies the tension this effect causes and shows that it is sufficiently large to render the thermal model inconsistent: its natural assumptions are in conflict with its outputs. The possibility that hypertritons are formed at freeze out as compact objects, quark droplets, that subsequently evolve into hypertritons is considered as a way to resolve the puzzle. It is noted that beyond making the assumption that compact quark droplets form, additional detailed dynamical assumptions which have not been justified are needed to make the thermal model work. The issue of why, despite these issues, the hypertriton is well described by a simple statistical description at freeze out is unresolved. Resolving the hypertriton puzzle is important as it may clarify whether the phenomenological success of the simple thermal model for yields accurately reflects the simple picture of the underlying physics on which it is based.

hydrodynamic models↗
DOE OSTI2021

Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation

The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for γ astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of γ cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of γ absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift z=2 and to constrain or detect γ halos up to intergalactic-magnetic-field strengths of at least 0.3 pG . Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from γ astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of γ cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2013

Influence of Computational Drop Representation in LES of a Droplet-Laden Mixing Layer

Multiphase turbulent flows are encountered in many practical applications including turbine engines or natural phenomena involving particle dispersion. Numerical computations of multiphase turbulent flows are important because they provide a cheaper alternative to performing experiments during an engine design process or because they can provide predictions of pollutant dispersion, etc. Two-phase flows contain millions and sometimes billions of particles. For flows with volumetrically dilute particle loading, the most accurate method of numerically simulating the flow is based on direct numerical simulation (DNS) of the governing equations in which all scales of the flow including the small scales that are responsible for the overwhelming amount of dissipation are resolved. DNS, however, requires high computational cost and cannot be used in engineering design applications where iterations among several design conditions are necessary. Because of high computational cost, numerical simulations of such flows cannot track all these drops. The objective of this work is to quantify the influence of the number of computational drops and grid spacing on the accuracy of predicted flow statistics, and to possibly identify the minimum number, or, if not possible, the optimal number of computational drops that provide minimal error in flow prediction. For this purpose, several Large Eddy Simulation (LES) of a mixing layer with evaporating drops have been performed by using coarse, medium, and fine grid spacings and computational drops, rather than physical drops. To define computational drops, an integer NR is introduced that represents the ratio of the number of existing physical drops to the desired number of computational drops; for example, if NR=8, this means that a computational drop represents 8 physical drops in the flow field. The desired number of computational drops is determined by the available computational resources; the larger NR is, the less computationally intensive is the simulation. A set of first order and second order flow statistics, and of drop statistics are extracted from LES predictions and are compared to results obtained by filtering a DNS database. First order statistics such as Favre averaged stream-wise velocity, Favre averaged vapor mass fraction, and the drop stream-wise velocity, are predicted accurately independent of the number of computational drops and grid spacing. Second order flow statistics depend both on the number of computational drops and on grid spacing. The scalar variance and turbulent vapor flux are predicted accurately by the fine mesh LES only when NR is less than 32, and by the coarse mesh LES reasonably accurately for all NR values. This is attributed to the fact that when the grid spacing is coarsened, the number of drops in a computational cell must not be significantly lower than that in the DNS.

Bellan, Josette↗
DOE OSTI2024

Recent trends and variability of temperature and atmospheric water vapor over South Asia

This study examines the recent trends and variability of temperature and atmospheric water vapor over South Asia during the summer monsoon season from 1970 to 2022. We employ gridded reanalysis and observational datasets to assess the trends and variability in surface temperature and further investigate the associated physical processes. Our results show that the surface temperature over South Asia has increased by 0.5 °C per decade. Notably, we find a statistically significant major transition in the surface temperature series in 1995, such that 1970–1994 (1995–2022) is associated with anomalous cooling (anomalous warming). Since 1995, not only have the temporal variations in regional surface temperature, atmospheric water vapor, and clear-sky downward longwave radiation been synchronous at interannual scales, but their largest spatial variations have been located mainly over the same region. A diagnosis of the atmospheric moisture budget reveals that the observed changes in water vapor convergence over South Asia are dominated by changes in mean circulation dynamics, while the change in thermodynamics is relatively small. The enhancement of the Somali low-level jet may contribute to increased water vapor in South Asia by modulating low-level circulations.

Atmospheric circulations↗
DOE OSTI2021

Surface processing and discharge-conditioning of high voltage electrodes for the Ra EDM experiment

The Ra EDM experiment uses a pair of high voltage electrodes to search for the atomic electric dipole moment of 225 Ra. We use identical, plane–parallel electrodes with a primary high gradient surface of 200 mm 2 to generate reversible DC electric fields. Our statistical sensitivity is linearly proportional to the electric field strength in the electrode gap. We adapted surface decontamination and processing techniques from accelerator physics literature to chemical polish and clean a suite of newly fabricated large-grain niobium and grade-2 titanium electrodes. Three pairs of niobium electrodes and one pair of titanium electrodes were discharge-conditioned with a custom high voltage test station at electric field strengths as high as +52.5 kV/mm and –51.5 kV/mm over electrode gap sizes ranging from 0.4 mm to 2.5 mm. One pair of large-grain niobium electrodes was discharge-conditioned and validated to operate at ±20 kV/mm with steady-state leakage current ≤ 25 pA (1σ) and a polarity-averaged 98 ± 19 discharges per hour. These electrodes were installed in the Ra EDM experimental apparatus, replacing a copper electrode pair, and were revalidated to ±20 kV/mm. Here, the niobium electrodes perform at an electric field strength 3.1 times larger than the legacy copper electrodes and are ultimately limited by the maximum output of our 30 kV bipolar power supply.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Confinement-Driven Segregation Enables Glassy Polymer Hybrid Materials Featuring Disordered Hyperuniformity and Integrated Self-Healing

The blending of glassy copolymer-brush modified colloids with a viscoelastic linear copolymer featuring intrinsic self-healing enabled disordered hyperuniform hybrid materials that combined mechanical robustness with structural color, processability, environmental stability, and the ability to recover structure and properties after incurring physical damage via ‘integrated self-healing’. Symmetric linear n-butyl acrylate/methyl methacrylate (BA/MMA) were co-assembled with asymmetric glassy BA/MMA statistical copolymer brush (silica) particles. ‘Confinement-driven segregation’ resulted in a microphase-separated morphology in which the linear copolymer resided within the interstitial regions of a rigid (∼1 GPa) copolymer brush particle template with disordered hyperuniform microstructure. Diffusion of the self-healing copolymer additive into damage regions drove the recovery after damage, along with the restoration of structural color due to the materials hyperuniform microstructure. The synergistic action of intrinsic and extrinsic healing mechanisms could provide a versatile platform for the bottom-up fabrication of multifunctional hybrid materials with increased damage resistance and functional longevity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

One-shot gas detection with transformer paired neural networks in Mako collected longwave infrared hyperspectral imagery

To date, careful data treatment workflows and statistical detectors are used to perform hyperspectral image (HSI) detection of any gas contained in a spectral library, which is often expanded with physics models to incorporate different spectral characteristics. In general, surrounding evidence or known gas-release parameters are used to provide confidence in or confirm detection capability, respectively. This makes quantifying detection performance difficult as it is nearly impossible to develop an absolute ground truth for gas target pixel presence in collected HSI. Consequently, developing and comparing new detection methods, especially machine learning (ML)-based methods, is susceptible to subjectivity in derived detection map quality. Here, in this work, we demonstrate the first use of transformer-based paired neural networks (PNNs) for one-shot gas target detection for multiple gases while providing quantitative classification and detection metrics for their use on labeled data. Terabytes of training data are generated from a database of long-wave infrared HSI obtained from historical Mako sensor campaigns over Los Angeles. By incorporating labels, singular signature representations, and a model development pipeline, we can tune and select PNNs to detect multiple gas targets that are not seen in training on a quantitative basis. We additionally assess our test set detections using interpretability techniques widely employed with ML-based predictors, but less common with detection methods relying on learned latent spaces.

Hyperspectral imaging↗
DOE OSTI2026

Improving MicroBooNE's Inclusive Single-Photon Search with Low-Energy Hadronic Identification

This analysis aims to further investigate MicroBooNE’s inclusive single-photon search results, which reported a 2.2$\sigma$ excess below 600 MeV in shower energy for events with no reconstructed protons using roughly half of MicroBooNE's dataset. Taken together with MiniBooNE’s long-standing low-energy excess and MicroBooNE’s recent electron-like search results showing no observable excess with respect to Standard Model predictions, this result provides strong impetus for expanded exploration of the single-photon channel in Fermilab’s short-baseline liquid-argon time projection chamber (LArTPC) experiments. By identifying and classifying isolated MeV-scale energy depositions, or blips, in the vicinity of single-photon events selected by the Wire-Cell reconstruction framework, we establish a more comprehensive labeling scheme for nearby hadronic content. In particular, blips found backwards along the shower axis indicate the presence of previously-unidentified final-state protons, while elevated blip counts at wide angles signal the presence of final-state neutrons. By applying this new technique to its full dataset, MicroBooNE will perform a purer and higher-statistics test of the truly isolated nature of its modest photon-like excess, furthering its hunt for the presence of unexpected new physics beyond the Standard Model.

Andrade Aldana, Diego [Los Alamos; IIT, Chicago (m↗
NASA NTRS2001

Rotation of EOFs by the Independent Component Analysis: Towards A Solution of the Mixing Problem in the Decomposition of Geophysical Time Series

The Independent Component Analysis is a recently developed technique for component extraction. This new method requires the statistical independence of the extracted components, a stronger constraint that uses higher-order statistics, instead of the classical decorrelation, a weaker constraint that uses only second-order statistics. This technique has been used recently for the analysis of geophysical time series with the goal of investigating the causes of variability in observed data (i.e. exploratory approach). We demonstrate with a data simulation experiment that, if initialized with a Principal Component Analysis, the Independent Component Analysis performs a rotation of the classical PCA (or EOF) solution. This rotation uses no localization criterion like other Rotation Techniques (RT), only the global generalization of decorrelation by statistical independence is used. This rotation of the PCA solution seems to be able to solve the tendency of PCA to mix several physical phenomena, even when the signal is just their linear sum.

Aires, Filipe↗
NASA NTRS2000

Rotation of EOFs by the Independent Component Analysis: Towards a Solution of the Mixing Problem in the Decomposition of Geophysical Time Series

The Independent Component Analysis is a recently developed technique for component extraction. This new method requires the statistical independence of the extracted components, a stronger constraint that uses higher-order statistics, instead of the classical decorrelation, a weaker constraint that uses only second-order statistics. This technique has been used recently for the analysis of geophysical time series with the goal of investigating the causes of variability in observed data (i.e. exploratory approach). We demonstrate with a data simulation experiment that, if initialized with a Principal Component Analysis, the Independent Component Analysis performs a rotation of the classical PCA (or EOF) solution. This rotation uses no localization criterion like other Rotation Techniques (RT), only the global generalization of decorrelation by statistical independence is used. This rotation of the PCA solution seems to be able to solve the tendency of PCA to mix several physical phenomena, even when the signal is just their linear sum.

Aires, Filipe↗