Search NASA⌕ Search

SEARCH · Search NASA

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

A histogram-free multicanonical Monte Carlo algorithm for the construction of analytical density of states

We report a new multicanonical Monte Carlo (MC) algorithm to obtain the density of states (DOS) for physical systems with continuous state variables in statistical mechanics. Our algorithm is able to obtain an analytical form for the DOS expressed in a chosen basis set, instead of a numerical array of finite resolution as in previous variants of this class of MC methods such as the multicanonical (MUCA) sampling and Wang-Landau (WL) sampling. This is enabled by storing the visited states directly in a data set and avoiding the explicit collection of a histogram. This practice also has the advantage of avoiding undesirable artificial errors caused by the discretization and binning of continuous state variables. Our results show that this scheme is capable of obtaining converged results with a much reduced number of Monte Carlo steps, leading to a significant speedup over existing algorithms.

Li, Ying Wai↗
DOE OSTI2018

The Muon g-2 Straw Tracking Detectors

• Built to measure the muon anomalous magnetic moment to an uncertainty of 140 ppb and search for a muon electric dipole moment. • Determine whether discrepancy from the Standard Model in the previous experiment $\Delta 𝒂_{\mu}$ 𝑬𝟖𝟐𝟏 − 𝑺𝑴 = 𝟐𝟓𝟓 ± 𝟖𝟎 ×$𝟏𝟎^{-𝟏𝟏}$ is a statistical fluctuation or indicates non-SM physics. • A muon storage ring located at Fermilab, Illinois. • 15m diameter storage ring. • 1.45T uniform magnetic field. • 24 calorimeters around the ring. • 2 tracking stations. • Storing 100s of billions of muons.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

An Energy Calculator for Simple Commercial Buildings

According to the EIA, simple commercial buildings account for 97% of total commercial building stock. However, most simple commercial buildings for example small- to mid-sized offices, retail, schools and warehouses do not benefit from the data-driven decision-making capabilities of whole-building energy modeling. The high cost of custom modeling limits the use of energy modeling of simple buildings for new construction or retrofit measures. Lack of tools providing helpful information on interactive savings estimates creates difficulties in meeting aggressive decarbonization and energy efficiency goals for simple building designers and utility program managers. This paper reviews a beta phase Simple Building Calculator with the ability to generate relatively accurate and interactive modeling results based on a limited but robust set of inputs. It can evaluate whole-building or single measure savings in new or existing buildings, compare measure package choices, or provide simplified performance modeling for energy codes and utility incentives. The tool combines physical (annual whole building prototype simulation) and statistical modeling techniques to predict annual energy performance. It supports a variety of building characteristics for envelope, HVAC, and lighting with parameters ranging from vintage to max tech configurations, as well as support for single-zone and simple multi-zone HVAC systems. The Simple Building Calculator was designed to provide immediate feedback for otherwise computationally intensive tasks like measure comparison, development of multiple measure package combinations, or verification that measures meet efficiency targets—all with the goal of providing a tool for quick annual energy simulation of simple commercial buildings.

Hart, Reid↗
DOE OSTI2020

Realizable hyperuniform and nonhyperuniform particle configurations with targeted spectral functions via effective pair interactions

The capacity to identify realizable many-body configurations associated with targeted functional forms for the pair correlation function $g_2(r)$ or its corresponding structure factor $\textit{S(k)}$ is of great fundamental and practical importance. While there are obvious necessary conditions that a prescribed structure factor at number density ρ must satisfy to be configurationally realizable, sufficient conditions are generally not known due to the infinite degeneracy of configurations with different higher-order correlation functions. A major aim of this paper is to expand our theoretical knowledge of the class of pair correlation functions or structure factors that are realizable by classical disordered ensembles of particle configurations, including exotic “hyperuniform” varieties. We first introduce a theoretical formalism that provides a means to draw classical particle configurations from canonical ensembles with certain pairwise-additive potentials that could correspond to targeted analytical functional forms for the structure factor. This formulation enables us to devise an improved algorithm to construct systematically canonical-ensemble particle configurations with such targeted pair statistics, whenever realizable. As a proof of concept, we test the algorithm by targeting several different structure factors across dimensions that are known to be realizable and one hyperuniform target that is known to be nontrivially unrealizable. Our algorithm succeeds for all realizable targets and appropriately fails for the unrealizable target, demonstrating the accuracy and power of the method to numerically investigate the realizability problem. Subsequently, we also target several families of structure-factor functions that meet the known necessary realizability conditions but are not known to be realizable by disordered hyperuniform point configurations, including $\textit{d}$-dimensional Gaussian structure factors, $\textit{d}$-dimensional generalizations of the two-dimensional one-component plasma (OCP), and the $\textit{d}$-dimensional Fourier duals of the previous OCP cases. Moreover, we also explore unusual nonhyperuniform targets, including “hyposurficial” and “antihyperuniform” examples. In all of these instances, the targeted structure factors are achieved with high accuracy, suggesting that they are indeed realizable by equilibrium configurations with pairwise interactions at positive temperatures. Remarkably, we also show that the structure factor of nonequilibrium perfect glass, specified by two-, three-, and four-body interactions, can also be realized by equilibrium pair interactions at positive temperatures. In this work, our findings lead us to the conjecture that any realizable structure factor corresponding to either a translationally invariant equilibrium or nonequilibrium system can be attained by an equilibrium ensemble involving only effective pair interactions. Our investigation not only broadens our knowledge of analytical functional forms for $g_2(r)$ and $\textit{S(k)}$ associated with disordered point configurations across dimensions but also deepens our understanding of many-body physics. Moreover, our work can be applied to the design of materials with desirable physical properties that can be tuned by their pair statistics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Big data health physics: a case study in operational health physics

This paper, written as a two part project as required by Illinois Institute of Technology’s Phys597 Reading & Special Problems class and health physics professional master’s program, focuses on the results of a prototype data collection process implemented at a Los Alamos National Laboratory facility. In particular, this paper examines the data in a statistical context in an effort to optimize operational health physics processes at the Los Alamos National Laboratory.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Modeling Flow and Transport in Fractured Rock using Machine Learning [Slides]

Fractured systems in the subsurface play a role in many natural and engineered applications such as geologic carbon sequestration, hydraulic fracturing and underground nuclear test detection. Structural information (fracture size, orientation, etc.) plays a key role in governing the dominant physics for these systems but can only be known statistically. Traditional modeling approaches either ignore or idealize structural information at these larger scales because we lack a computational framework that utilizes it in its entirety. The work presented here integrates computational physics, machine learning and graph theory to make a paradigm shift from computationally intensive high-fidelity models to coarse-scale graphs without loss of critical structural information. We exploit the underlying discrete structure of fracture networks in systems considering flow through fractures and fracture propagation. We demonstrate that compact graph representations require significantly fewer degrees of freedom (dof) to capture micro-fracture information and further accelerate these models with Machine Learning.

58 GEOSCIENCES↗
NASA NTRS1986

Physics through the 1990s: Elementary-particle physics

The volume begins with a non-mathematical discussion of the motivation behind, and basic ideas of, elementary-particle physics theory and experiment. The progress over the past two decades with the quark model and unification of the electromagnetic and weak interactions is reviewed. Existing theoretical problems in the field, such as the origin of mass and the unification of the fundamental forces, are detailed, along with experimental programs to test the new theories. Accelerators, instrumentation, and detectors are described for both current and future facilities. Interactions with other areas of both theoretical and applied physics are presented. The sociology of the field is examined regarding the education of graduate students, the organization necessary in large-scale experiments, and the decision-making process involved in high-cost experiments. Finally, conclusions and recommendations for maintaining US excellence in theory and experiment are given. Appendices list both current and planned accelerators, and present statistical data on the US elementary-particle physics program. A glossary is included.

Source record↗
NASA NTRS2012

Likelihood-Based Climate Model Evaluation

Climate models are deterministic, mathematical descriptions of the physics of climate. Confidence in predictions of future climate is increased if the physics are verifiably correct. A necessary, (but not sufficient) condition is that past and present climate be simulated well. Quantify the likelihood that a (summary statistic computed from a) set of observations arises from a physical system with the characteristics captured by a model generated time series. Given a prior on models, we can go further: posterior distribution of model given observations.

climate models↗
NASA NTRS2011

Applications of the BIOPHYS Algorithm for Physically-Based Retrieval of Biophysical, Structural and Forest Disturbance Information

Canopy reflectance model inversion using look-up table approaches provides powerful and flexible options for deriving improved forest biophysical structural information (BSI) compared with traditional statistical empirical methods. The BIOPHYS algorithm is an improved, physically-based inversion approach for deriving BSI for independent use and validation and for monitoring, inventory and quantifying forest disturbance as well as input to ecosystem, climate and carbon models. Based on the multiple-forward mode (MFM) inversion approach, BIOPHYS results were summarized from different studies (Minnesota/NASA COVER; Virginia/LEDAPS; Saskatchewan/BOREAS), sensors (airborne MMR; Landsat; MODIS) and models (GeoSail; GOMS). Applications output included forest density, height, crown dimension, branch and green leaf area, canopy cover, disturbance estimates based on multi-temporal chronosequences, and structural change following recovery from forest fires over the last century. Good correspondences with validation field data were obtained. Integrated analyses of multiple solar and view angle imagery further improved retrievals compared with single pass data. Quantifying ecosystem dynamics such as the area and percent of forest disturbance, early regrowth and succession provide essential inputs to process-driven models of carbon flux. BIOPHYS is well suited for large-area, multi-temporal applications involving multiple image sets and mosaics for assessing vegetation disturbance and quantifying biophysical structural dynamics and change. It is also suitable for integration with forest inventory, monitoring, updating, and other programs.

Peddle, Derek R.↗
DOE OSTI2023

Engineering Super–Poissonian Photon Statistics of Spatial Light Modes

The nature of light sources is defined by the statistical fluctuations of the electromagnetic field. As such, the photon statistics of light sources are typically associated with distinct emitters. Here, the possibility of producing light beams with various photon statistics through the spatial modulation of coherent light is demonstrated. This is achieved by the sequential encoding of controllable Kolmogorov phase screens in a digital micromirror device. Interestingly, the flexibility of this scheme allows for the shaping of spatial light modes with engineered photon statistics at different spatial positions. The performance of this scheme is assessed through the photon-number-resolving characterization of different families of spatial light modes with engineered photon statistics. Furthermore, it is believed that the possibility of controlling the photon fluctuations of the light field at arbitrary spatial locations has important implications for quantum spectroscopy, sensing, and imaging.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

Parameters Impacting Columnated Granular Soil Pneumatic Seal Performance

The ability of a column of loose granular soil to form a pneumatic seal was investigated by varying the diameter of the soil column, the effective column height, and the level of compaction in the soil. The soil column diameter was tested at three levels using pipes with inner diameters measuring 5.08, 10.16, and 15.24 centimeters (2, 4, and 6 inches). Soil was filled in each pipe to form 15, 30, and 45 centimeter (6, 12, and 18 inch) tall soil columns. Each diameter/height configuration was also tested at three levels of soil compaction, compared by calculating the bulk density of the soil with mass and volume measurements. GRC-1a simulant was used, with approximate low/medium/high bulk densities of 1.6, 1.75, and 1.9 g/cc achieved with a combination of vibration and tamping. The order of tests for a given column diameter was randomized and repeated three times. With the top of the soil open to atmosphere at room temperature, compressed air was injected through a small diffuser at the column base with several small downward-facing holes. The number and size of these holes was scaled such that a constant total inlet orifice area to column cross sectional area ratio was maintained for each column diameter. Inlet air pressure was slowly increased via a precision regulator to preserve quasi-static equilibrium in the soil column to minimize the impact of dynamics. Air pressure was increased all the way through the static and bubbling regimes until slugging or turbulent behavior was observed in the soil to ensure that the entire static regime had been captured during data collection. A three-factor, three-level analysis of variance (ANOVA) statistical analysis was performed on the resulting data to determine the extent to which each physical parameter impacted the soil column seal performance. It was concluded that there is statistically significant evidence that column height, and the interaction between column height and diameter impact soil seal performance. In all other cases there was insufficient data to identify a statistically significant causal relationship. Additionally, plots were generated comparing experimental data to the predictive formula developed by Ogino et al. for fluidized beds in 1993. Because this model was developed for industrial spouted fluidized beds, the accuracy of its output prior to fluidization in the static seal ‘edge case’ is unknown, especially considering in this application the working gas was diffused across the column base rather than being injected through a spout. Further, the fidelity of the model had not yet been tested with lunar soil simulants. Plotting the Ogino et al. model alongside test data allows for a more intuitive sense of the impact of test parameters on soil seal performance, as well as providing a quick means to further tune this predictive model for more accurate use with static seals across granular lunar soil simulants. The model provided by Ogino et al. was further tuned using test data to determine the degree to which a spouted bed model could be applied to a slightly modified set of testing conditions: lunar soil simulants and a more diffuse, homogeneous application of pneumatic pressure. A Matlab script was created to test different values for the leading coefficient and exponents in Ogino’s formula. The script swept preset ranges, then iterated with higher resolutions across narrower ranges to converge on the optimal value for each parameter. The resulting adjusted model was compared to the original, as well as test data with noticeable improvements across the entire test domain.

Jack Stewart↗
NASA NTRS1987

Is solar neutrino capture rate correlated with sunspot number?

The statistical significance of the apparent correlation between sunspots and the observed neutrino rate is quantified. It is shown that the correlation depends almost entirely upon four low neutrino capture rates near the beginning of 1980. A calculation based on standard electroweak theory and neutrino production processes demonstrates that a correlation, if real, would be extremely puzzling on energetic grounds alone. It is concluded that measurements with the Cl-37 detector during the next sunspot cycle will be needed to show that there is a physical correlation, since the existing data are not statistically significant at a definitive level.

Bahcall, J. N.↗
DOE OSTI2023

Modeling single-molecule stretching experiments using statistical thermodynamics

Single-molecule stretching experiments are widely utilized within the fields of physics and chemistry to characterize the mechanics of individual bonds or molecules, as well as chemical reactions. Analytic relations describing these experiments are valuable, and these relations can be obtained through the statistical thermodynamics of idealized model systems representing the experiments. Since the specific thermodynamic ensembles manifested by the experiments affect the outcome, primarily for small molecules, the stretching device must be included in the idealized model system. Though the model for the stretched molecule might be exactly solvable, including the device in the model often prevents analytic solutions. In the limit of large or small device stiffness, the isometric or isotensional ensembles can provide effective approximations, but the device effects are missing. Here a dual set of asymptotically correct statistical thermodynamic theories are applied to develop accurate approximations for the full model system that includes both the molecule and the device. In conclusion, the asymptotic theories are first demonstrated to be accurate using the freely jointed chain model and then using molecular dynamics calculations of a single polyethylene chain.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

To have value, comparisons of high-throughput phenotyping methods need statistical tests of bias and variance

The gap between genomics and phenomics is narrowing. The rate at which it is narrowing, however, is being slowed by improper statistical comparison of methods. Quantification using Pearson’s correlation coefficient ( r ) is commonly used to assess method quality, but it is an often misleading statistic for this purpose as it is unable to provide information about the relative quality of two methods. Using r can both erroneously discount methods that are inherently more precise and validate methods that are less accurate. These errors occur because of logical flaws inherent in the use of r when comparing methods, not as a problem of limited sample size or the unavoidable possibility of a type I error. A popular alternative to using r is to measure the limits of agreement (LOA). However both r and LOA fail to identify which instrument is more or less variable than the other and can lead to incorrect conclusions about method quality. An alternative approach, comparing variances of methods, requires repeated measurements of the same subject, but avoids incorrect conclusions. Variance comparison is arguably the most important component of method validation and, thus, when repeated measurements are possible, variance comparison provides considerable value to these studies. Statistical tests to compare variances presented here are well established, easy to interpret and ubiquitously available. The widespread use of r has potentially led to numerous incorrect conclusions about method quality, hampering development, and the approach described here would be useful to advance high throughput phenotyping methods but can also extend into any branch of science. The adoption of the statistical techniques outlined in this paper will help speed the adoption of new high throughput phenotyping techniques by indicating when one should reject a new method, outright replace an old method or conditionally use a new method.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Proton imaging of stochastic magnetic fields

Recent laser-plasma experiments report the existence of dynamically significant magnetic fields, whose statistical characterization is essential for a complete understanding of the physical processes these experiments are attempting to investigate. In this paper, we show how a proton-imaging diagnostic can be used to determine a range of relevant magnetic-field statistics, including the magnetic-energy spectrum. To achieve this goal, we explore the properties of an analytic relation between a stochastic magnetic field and the image-flux distribution created upon imaging that field. This ‘Kugland image-flux relation’ was previously derived under simplifying assumptions typically valid in actual proton-imaging set-ups. We conclude that, as with regular electromagnetic fields, features of the beam’s final image-flux distribution often display a universal character determined by a single, field-scale dependent parameter – the contrast parameter $\unicode[STIX]{x1D707}\equiv d_{s}/{\mathcal{M}}l_{B}$ – which quantifies the relative size of the correlation length $l_{B}$ of the stochastic field, proton displacements $d_{s}$ due to magnetic deflections and the image magnification ${\mathcal{M}}$ . For stochastic magnetic fields, we establish the existence of four contrast regimes, under which proton-flux images relate to their parent fields in a qualitatively distinct manner. These are linear, nonlinear injective, caustic and diffusive. The diffusive regime is newly identified and characterized. The nonlinear injective regime is distinguished from the caustic regime in manifesting nonlinear behaviour, but as in the linear regime, the path-integrated magnetic field experienced by the beam can be extracted uniquely. Thus, in the linear and nonlinear injective regimes we show that the magnetic-energy spectrum can be obtained under a further statistical assumption of isotropy. This is not the case in the caustic or diffusive regimes. We discuss complications to the contrast-regime characterization arising for inhomogeneous, multi-scale stochastic fields, which can encompass many contrast regimes, as well as limitations currently placed by experimental capabilities on one’s ability to extract magnetic-field statistics. Furthermore, the results presented in this paper are of consequence in providing a comprehensive description of proton images of stochastic magnetic fields, with applications for improved analysis of proton-flux images.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A variance deconvolution estimator for efficient uncertainty quantification in Monte Carlo radiation transport applications

Monte Carlo simulations are at the heart of many high-fidelity simulations and analyses for radiation transport systems. As is the case with any complex computational model, it is important to propagate sources of input uncertainty and characterize how they affect model output. Unfortunately, uncertainty quantification (UQ) is made difficult by the stochastic variability that Monte Carlo transport solvers introduce. The standard method to avoid corrupting the UQ statistics with the transport solver noise is to increase the number of particle histories, resulting in very high computational costs. In this contribution, we propose and analyze a sampling estimator based on the law of total variance to compute UQ variance even in the presence of residual noise from Monte Carlo transport calculations. Here, we rigorously derive the statistical properties of the new variance estimator, compare its performance to that of the standard method, and demonstrate its use on neutral particle transport model problems involving both attenuation and scattering physics. We illustrate, both analytically and numerically, the estimator’s statistical performance as a function of available computational budget and the distribution of that budget between UQ samples and particle histories. We show analytically and corroborate numerically that the new estimator is unbiased, unlike the standard approach, and is more accurate and precise than the standard estimator for the same computational budget.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

A Qualitative Strategy for Fusion of Physics into Empirical Models for Process Anomaly Detection

To facilitate the automated online monitoring of power plants, a systematic and qualitative strategy for anomaly detection is presented. This strategy is essential to provide credible reasoning on why and when an empirical versus hybrid (i.e., physics-supported) approach should be used and to determine the ideal mix of these two approaches for a defined anomaly detection scope. Empirical methods are usually based on pattern, statistical, and causal inference. Hybrid methods include the use of physics models to train and test data methods, reduce data dimensionality, reduce data-model complexity, augment data, and reduce empirical uncertainty; hybrid methods also include the use of data to tune physics models. The presented strategy is driven by key decision points related to data relevance, simple modeling feasibility, data inference, physics-modeling value, data dimensionality, physics knowledge, method of validation, performance, data availability, and suitability for training and testing, cause-effect, entropy inference, and model fitting. The strategy is demonstrated through a pilot use case for the application of anomaly detection to capture a valve packing leak at the high-pressure coolant injection system of a nuclear power plant.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Search for dijet resonances with data scouting in proton-proton collisions at $\sqrt{s}=13$ TeV

A search is presented for narrow resonances, with a mass between 0.6 and 1.8 TeV, decaying to pairs of jets, in proton-proton collisions at $\sqrt{s}=13$ TeV. The search is performed using dijets that are reconstructed, selected, and recorded in a compact form by the high-level trigger in a technique referred to as “data scouting”, from data collected in 2016–2018 corresponding to an integrated luminosity of 117 fb −1 . The dijet mass spectra are well described by a smooth parameterization, and no significant evidence for the production of new particles is observed. Model-independent upper limits are presented on the product of the cross section, branching fraction, and acceptance for the individual cases of narrow quark-quark, quark-gluon, and gluon-gluon resonances, and are compared to the predictions from a variety of models of narrow dijet resonance production. The upper limit on the coupling of a dark matter mediator to quarks is presented as a function of the mediator mass. The sensitivity of this search goes beyond what is expected from statistical scaling with the integrated luminosity alone, as a consequence of the use of fewer parameters in the background function within a more robust statistical procedure.

beyond Standard Model↗