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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Some statistical analysis techniques for the effect of systematic instrument errors on the results of measurements in a space experiment

Mathematical models of certain systematic telemetry signal errors arising during the passage of the signal through the recording and transmitting blocks of the physical instruments aboard the spacecraft are classified and constructed. Statistical analysis methods are used to determine the background function and relative sensitivity coefficients of the channels of a charged particle detector. Results are presented. A block diagram for a method taking into account the effect of distortions is given.

Kondakov, S. P.↗
DOE OSTI2018

A two-phase mixing layer between parallel gas and liquid streams: multiphase turbulence statistics and influence of interfacial instability

The two-phase mixing layer formed between parallel gas and liquid streams is an important fundamental problem in turbulent multiphase flows. The problem is relevant to many industrial applications and natural phenomena, such as air-blast atomizers in fuel injection systems and breaking waves in the ocean. The velocity difference between the gas and liquid streams triggers an interfacial instability which can be convective or absolute depending on the stream properties and injection parameters. In the present study, a direct numerical simulation of a two-phase gas–liquid mixing layer that lie in the absolute instability regime is conducted. A dominant frequency is observed in the simulation and the numerical result agrees well with the prediction from viscous stability theory. As the interfacial wave plays a critical role in turbulence transition and development, the temporal evolution of turbulent fluctuations (such as the enstrophy) also exhibits a similar frequency. To investigate the statistical response of the multiphase turbulence flow, the simulation has been run for a long physical time so that time-averaging can be performed to yield the statistically converged results for Reynolds stresses and the turbulent kinetic energy (TKE) budget. An extensive mesh refinement study using from 8 million to about 4 billions cells has been performed. The turbulent dissipation is shown to be highly demanding on mesh resolution compared with other terms in TKE budget. The results obtained with the finest mesh are shown to be close to converged results of turbulent dissipation which allow us to obtain estimations of the Kolmogorov and Hinze scales. The estimated Kolmogorov scale is found to be similar to the cell size of the finest mesh used here. The computed Hinze scale is significantly larger than the size of droplets observed and does not seem to be a relevant length scale to describe the smallest size of droplets formed in atomization.

Mechanics↗
DOE OSTI2021

Physics‐based iterative reconstruction for dual‐source and flying focal spot computed tomography

Purpose For single‐source helical Computed Tomography (CT), both Filtered‐Back Projection (FBP) and statistical iterative reconstruction have been investigated. However, for dual‐source CT with flying focal spot (DS‐FFS CT), a statistical iterative reconstruction that accurately models the scanner geometry and acquisition physics remains unknown to researchers. Therefore, our purpose is to present a novel physics‐based iterative reconstruction method for DS‐FFS CT and assess its image quality. Methods Our algorithm uses precise physics models to reconstruct from the native cone‐beam geometry and interleaved dual‐source helical trajectory of a DS‐FFS CT. To do so, we construct a noise physics model to represent data acquisition noise and a prior image model to represent image noise and texture. In addition, we design forward system models to compute the locations of deflected focal spots, the dimension, and sensitivity of voxels and detector units, as well as the length of intersection between x‐rays and voxels. The forward system models further represent the coordinated movement between the dual sources by computing their x‐ray coverage gaps and overlaps at an arbitrary helical pitch. With the above models, we reconstruct images by an advanced Consensus Equilibrium (CE) numerical method to compute the maximum a posteriori estimate to a joint optimization problem that simultaneously fits all models. Results We compared our reconstruction with Siemens ADMIRE, which is the clinical standard hybrid iterative reconstruction (IR) method for DS‐FFS CT, in terms of spatial resolution, noise profile, and image artifacts through both phantoms and clinical scan datasets. Experiments show that our reconstruction has a higher spatial resolution, with a Task‐Based Modulation Transfer Function (MTF task ) consistently higher than the clinical standard hybrid IR. In addition, our reconstruction shows a reduced magnitude of image undersampling artifacts than the clinical standard. Conclusions By modeling a precise geometry and avoiding data rebinning or interpolation, our physics‐based reconstruction achieves a higher spatial resolution and fewer image artifacts with smaller magnitude than the clinical standard hybrid IR.

Wang, Xiao↗
DOE OSTI2022

HL-LHC Analysis With ROOT: ROOT Project Input to the HL-LHC Computing Review (Stage 2)

ROOT is high energy physics' software for storing and mining data in a statistically sound way, to publish results with scientific graphics. It is evolving since 25 years, now providing the storage format for more than one exabyte of data; virtually all high energy physics experiments use ROOT. With another significant increase in the amount of data to be handled scheduled to arrive in 2027, ROOT is preparing for a massive upgrade of its core ingredients. As part of a review of crucial software for high energy physics, the ROOT team has documented its R&D plans for the coming years.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Computational and analytical comparison of flux discretizations for the semiconductor device equations beyond Boltzmann statistics

We compare three thermodynamically consistent numerical fluxes known in the literature, appearing in a Voronoï finite volume discretization of the van Roosbroeck system with general charge carrier statistics. Our discussion includes an extension of the Scharfetter–Gummel scheme to non-Boltzmann (e.g. Fermi–Dirac) statistics. It is based on the analytical solution of a two-point boundary value problem obtained by projecting the continuous differential equation onto the interval between neighboring collocation points. Hence, it serves as a reference flux. The exact solution of the boundary value problem can be approximated by computationally cheaper fluxes which modify certain physical quantities. One alternative scheme averages the nonlinear diffusion (caused by the non-Boltzmann nature of the problem), another one modifies the effective density of states. To study the differences between these three schemes, we analyze the Taylor expansions, derive an error estimate, visualize the flux error and show how the schemes perform for a carefully designed p-i-n benchmark simulation. We present strong evidence that the flux discretization based on averaging the nonlinear diffusion has an edge over the scheme based on modifying the effective density of states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Energy considerations in the Community Atmosphere Model (CAM)

An error in the energy formulation in the Community Atmosphere Model (CAM) is identified and corrected. Ten year AMIP simulations are compared using the correct and incorrect energy formulations. Statistics of selected primary variables all indicate physically insignificant differences between the simulations, comparable to differences with simulations initialized with rounding sized perturbations. The two simulations are so similar mainly because of an inconsistency in the application of the incorrect energy formulation in the original CAM. CAM used the erroneous energy form to determine the states passed between the parameterizations, but used a form related to the correct formulation for the state passed from the parameterizations to the dynamical core. If the incorrect form is also used to determine the state passed to the dynamical core the simulations are significantly different. In addition, CAM uses the incorrect form for the global energy fixer, but that seems to be less important. The difference of the magnitude of the fixers using the correct and incorrect energy definitions is very small.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Assessing concurrent effects of climate change on hydropower supply, electricity demand, and greenhouse gas emissions in the Upper Yangtze River Basin of China

Hydropower importantly provides flexible low-carbon electricity, however, climate change will affect the hydropower system through altering hydrologic regimes while also affecting electricity demands for heating and cooling that hydropower resources serve. This study assesses the effect of climate change on hydropower and electricity demand in the Upper Yangtze River Basin (UYRB) in China on the regional net electric load and greenhouse gas (GHG) emissions. This is accomplished by using climate projections from five global climate models (GCMs) to simultaneously force (1) a physically-based hydrological model and a statistically-based hydropower model to estimate the future generating capacity of 21 large hydropower plants in the UYRB and (2) an empirical electricity demand model accounting for socioeconomic and climatic factors. Under climate change, the projected hydropower generation in the UYRB tends to increase in the 21st century but is far less than the increase in electricity demand, increasing the gap between demand and supply. Future increases in overall electricity demand are driven by GDP growth, but climate change will alter the distribution of the seasonal electricity demand. Climate warming decreases electricity demand for heating in winter and increases electricity demand for cooling in summer, but ultimately increases demand. Meanwhile, there is an increasing mismatch between electricity demand and hydropower supply associated with inter- and intra-annual variations, owing to the temporal climate change and increase in compound climate extremes (droughts and heatwaves). Finally, meeting the gap between supply and demand due to climate change is estimated to contribute 79.0–184.6 and 50.6–316.2 MMT CO 2e /yr of additional GHG emissions by the mid and end of 21st century, respectively.

13 HYDRO ENERGY↗
DOE OSTI2021

Statistical description of coalescing magnetic islands via magnetic reconnection

The physical picture of interacting magnetic islands provides a useful paradigm for certain plasma dynamics in a variety of physical environments, such as the solar corona, the heliosheath and the Earth's magnetosphere. In this work, we derive an island kinetic equation to describe the evolution of the island distribution function (in area and in flux of islands) subject to a collisional integral designed to account for the role of magnetic reconnection during island mergers. This equation is used to study the inverse transfer of magnetic energy through the coalescence of magnetic islands in two dimensions. We solve our island kinetic equation numerically for three different types of initial distribution: Dirac delta, Gaussian and power-law distributions. The time evolution of several key quantities is found to agree well with our analytical predictions: magnetic energy decays as $\tilde {t}^{-1}$ , the number of islands decreases as $\tilde {t}^{-1}$ and the averaged area of islands grows as $\tilde {t}$ , where $\tilde {t}$ is the time normalised to the characteristic reconnection time scale of islands. General properties of the distribution function and the magnetic energy spectrum are also studied. Finally, we discuss the underlying connection of our island-merger models to the (self-similar) decay of magnetohydrodynamic turbulence.

Physics↗
DOE OSTI2025

Spinon excitations and spin correlations in the one-dimensional quantum magnet β -VOSO 4 probed by Raman spectroscopy

Fractionalized excitations such as spinons and anyons have emerged as a central theme in condensed matter physics with broad implications for superconductivity, quantum statistics, and quantum computation. The nearly ideal one-dimensional S = 1/2 system β-VOSO 4 without long-range order down to 85 mK provides a promising platform to experimentally explore such fractionalized excitations. Here, we employ Raman spectroscopy to probe magnetic excitations and the evolution of spin correlations in β-VOSO 4 . Spinon signatures are found along the chain direction, evidenced by a broad, gapless scattering continuum at low temperatures. The temperature dependence of the spinon spectral weight aligns considerably with numerical density matrix renormalization group calculations. By comparing the experimental spinon spectral weight with calculated results and evaluating the associated quantum Fisher information (QFI) therefrom, we observe a steep increase in QFI upon cooling, indicating rapidly growing correlation lengths. Our study showcases QFI as a probe of spin correlations in quantum magnets.

Wulferding, Dirk [Sejong Univ., Seoul (Korea, Repu↗
DOE OSTI2014

Modeling the neutral hydrogen distribution in the post-reionization Universe: intensity mapping

We model the distribution of neutral hydrogen (HI) in the post-reionization era and investigate its detectability in 21 cm intensity mapping with future radio telescopes like the Square Kilometer array (SKA). We rely on high resolution hydrodynamical N-body simulations that have a state-of-the-art treatment of the low density photoionized gas in the inter-galactic medium (IGM). The HI is assigned a-posteriori to the gas particles following two different approaches: a halo-based method in which HI is assigned only to gas particles residing within dark matter halos; a particle-based method that assigns HI to all gas particles using a prescription based on the physical properties of the particles. The HI statistical properties are then compared to the observational properties of Damped Lyman-α Absorbers (DLAs) and of lower column density systems and reasonable good agreement is found for all the cases. Among the halo-based method, we further consider two different schemes that aim at reproducing the observed properties of DLAs by distributing HI inside halos: one of this results in a much higher bias for DLAs, in agreement with recent observations, which boosts the 21 cm power spectrum by a factor ∼ 4 with respect to the other recipe. Furthermore, we quantify the contribution of HI in the diffuse IGM to both Ω{sub HI} and the HI power spectrum finding to be subdominant in both cases. We compute the 21 cm power spectrum from the simulated HI distribution and calculate the expected signal for both SKA1-mid and SKA1-low configurations at 2.4 ≤ z ≤ 4. We find that SKA will be able to detect the 21 cm power spectrum, in the non-linear regime, up to k ∼ 1 h/Mpc for SKA1-mid and k ∼ 5 h/Mpc for SKA1-low with 100 hours of observations. We also investigate the perspective of imaging the HI distribution. Our findings indicate that SKA1-low could detect the most massive HI peaks with a signal to noise ratio (SNR) higher than 5 for an observation time of about 1000 hours at z = 4, for a synthesized beam width of 2'. Detection at redshifts z≥2.4 with SKA1-mid would instead require a much longer observation time to achieve a comparable SNR level.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Controlling Strain Bursts and Avalanches at the Nano- to Micrometer Scale

We demonstrate, through three-dimensional discrete dislocation dynamics simulations, that the complex dynamical response of nano- and microcrystals to external constraints can be tuned. Under load rate control, strain bursts are shown to exhibit scale-free avalanche statistics, similar to critical phenomena in many physical systems. For the other extreme of displacement rate control, strain burst response transitions to quasiperiodic oscillations, similar to stick-slip earthquakes. External load mode control is shown to enable a qualitative transition in the complex collective dynamics of dislocations from self-organized criticality to quasiperiodic oscillations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Sci—Fri PM: Dosimetry—05: Megavoltage electron backscatter: EGSnrc results versus 21 experiments

The accuracy of electron backscatter calculations at megavoltage energies is important for many medical physics applications. In this study, EGSnrc calculations of megavoltage electron backscatter (1–22 MeV) are performed and compared to the data from 21 experiments published between 1954 and 1993 for 25 single elements with atomic numbers from 3 to 92. Typical experimental uncertainties are 15%. For EGSnrc simulations, an ideal detector is assumed, and the most accurate electron physics options are employed, for a combined statistical and systematic uncertainty of 3%. The quantities compared are the backscatter coefficient and the energy spectra (in the backward hemisphere and at specific detector locations). For the backscatter coefficient, the overall agreement is within ±2% in the absolute value of the backscatter coefficient (in per cent), and within 11% of the individual backscatter values. EGSnrc results are systematically on the higher end of the spread of the experimental data, which could be partially from systematic experimental errors discussed in the literature. For the energy spectra, reasonable agreement between simulations and experiments is observed, although there are significant variations in the experimental data. At the lower end of the spectra, simulations are higher than some experimental data, which could be due to reduced experimental sensitivity to lower energy electrons and/or over-estimation by EGSnrc for backscattered secondary electrons. In conclusion, overall good agreement is observed between EGSnrc backscatter calculations and experimental measurements for megavoltage electrons. There is a need for high quality experimental data for the energy spectra of backscattered electrons.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Features of the application of the Monte Carlo method to calculations for large RBMK reactors and to model correction on the basis of data from in-core detectors

Statistical errors in sampling neutron fields in physically large systems like an RBMK are analyzed both qualitatively and quantitatively. Recommendations concerning the choice of parameters for calculations are given. A new procedure for Monte Carlo RBMK calculations with model corrections on the basis of data from in-core detectors is proposed. Dedicated software based on the CUDA software and hardware platform is developed for computational research. Results of testing the procedure and software in question via calculations for real RBMK reactors are discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

A Histogram-Free Multicanonical Monte Carlo Algorithm for the Basis Expansion of 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.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Deep Learning Parameterization of Vertical Wind Velocity Variability via Constrained Adversarial Training

Atmospheric models with typical resolution in the tenths of kilometers cannot resolve the dynamics of air parcel ascent, which varies on scales ranging from tens to hundreds of meters. Small-scale wind fluctuations are thus characterized by a subgrid distribution of vertical wind velocity W with standard deviation σ W . The parameterization of σ W is fundamental to the representation of aerosol–cloud interactions, yet it is poorly constrained. Using a novel deep learning technique, this work develops a new parameterization for σ W merging data from global storm-resolving model simulations, high-frequency retrievals of W , and climate reanalysis products. The parameterization reproduces the observed statistics of σ W and leverages learned physical relations from the model simulations to guide extrapolation beyond the observed domain. Incorporating observational data during the training phase was found to be critical for its performance. The parameterization can be applied online within large-scale atmospheric models, or offline using output from weather forecasting and reanalysis products.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Energy Audit of MC-1

The goal of the Muon g-2 project is to measure the muon magnetic anomalous moment to a precision of 140 ppb. This is an unprecedented endeavor that will allow for the comparison of the standard model’s prediction and experimental value of the muon magnetic anomaly. If the difference is statistically significant this could mean there is physics beyond the standard model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Threshold photoproduction of J/ψ with the GlueX experiment

The GlueX experiment uses a linearly-polarized tagged-photon beam produced by electrons from the 12 GeV CEBAF machine. The detector system is based on a 2 T solenoid and includes e.m. calorimeters and drift chambers providing full acceptance coverage. The experiment allows us, for the first time, to study J/ψ photoproduction from the threshold, at 8.2 GeV, up to a photon energy of 12 GeV. Results from the recent commissioning run, where we took 25 billion triggers, will be presented. The precise knowledge of the beam energy, combined with the recoil proton reconstruction and the identification of the J/ψ decay products, e + /e - , allows us to achieve a high mass resolution. With the data from the upcoming physics runs, we expect to have enough statistics to be able to distinguish between different production mechanisms.

Pentchev, Lubomir↗
DOE OSTI2022

Advanced Laboratory and Field Arrays: Debris Modeling, Detection,& Mitigation (Task 1)

The statement of project objectives for this task was: develop tools, methods and models to assess, and mitigate the risk of damage to MHK infrastructure from woody debris. Develop the capability to detect woody debris using sonar and/or physical methods for purposes of characterizing debris statistics in river (at UAF’s Tanana River Test Site) and near-shore wave (at Yakutat, AK) environments and to activate debris mitigation measures. Develop debris impact risk maps and tables using statistics on debris size, geometry, type, prevalence, mobility and location. Improve and apply the COUPi discrete element method (DEM) to develop models of debris movement and impact on MHK infrastructure to evaluate risk of damage, and interference, to operations from debris. The proposed final deliverable for the task was a set of tools or techniques for providing estimates of the probability of debris impact, and resulting impact forces, on MHK infrastructure as a function of debris size, type, wave regime, and current velocity. Such estimates are required to assess damage risks to operational MHK infrastructure.

13 HYDRO ENERGY↗