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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Heat pumps for all? Distributions of the costs and benefits of residential air-source heat pumps in the United States

Electrification of fossil-fuel combustion in buildings is a key component of achieving global greenhouse gas emissions targets. We use physics simulations of 550,000 statistically representative households to analyze distributions of the costs and benefits of three air-to-air heat pump performance levels, with and without insulation upgrades, across the diversity of the US housing stock. We find positive greenhouse gas reductions in every US state for all performance levels across five 2022-2038 electric grid scenarios, with full adoption reducing national emissions by 5%-9%. We find that air-to-air heat pumps could be cost effective without subsidies in 59% of households (65 million). However, efficiency is key: whereas minimum-efficiency equipment could increase energy bills in 39% of households, this fraction is only 19% when also upgrading insulation or 5% when using higher-efficiency equipment, though both strategies have higher upfront costs. Such affordability challenges could be addressed through supportive incentives, policy, and innovation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Information-theoretic astrophysical uncertainties in the effective theory of dark matter direct detection

The impact of astrophysical uncertainties in direct detection searches can vary significantly across particle dark matter models and detector targets, due to the different velocity and momentum dependencies of the scattering cross section. We address these uncertainties for all operators of the nonrelativistic effective field theory of dark-matter/nucleon interactions, making use of the Kullback-Leibler (KL) information divergence to measure the deviation of the true dark matter velocity distribution from the Maxwell-Boltzmann form. This approach quantifies how astrophysical uncertainties affect each operator in the effective theory, without assuming any specific functional form for the velocity distribution. While for some operators the uncertainties are smaller than 1 order of magnitude for entropically motivated deviations from the Maxwell-Boltzmann form, for other operators, these uncertainties can be as large as three orders of magnitude near threshold. Furthermore, we identify the dependence of the scattering rate for various operators of the effective theory with different velocity-weighted moments of the velocity distribution, functionally analogous to the mean, variance, or skewness. This provides new analytic insight into which features of the velocity distribution are most relevant to detect a given particle dark matter model. Our technique is general and could be applied to a broader class of physics problems where a physical observable depends on the statistical moments of an uncertain theoretical distribution.

Herrera, Gonzalo [MIT, MKI; Harvard U.; Virginia T↗
DOE OSTI2023

Fast estimation of the look-elsewhere effect using Gaussian random fields

Abstract We discuss the use of Gaussian random fields to estimate the look-elsewhere effect correction. We show that Gaussian random fields can be used to model the null-hypothesis significance maps from a large set of statistical problems commonly encountered in physics, such as template matching and likelihood ratio tests. Some specific examples are searches for dark matter using pixel arrays, searches for astronomical transients, and searches for fast-radio bursts. Gaussian random fields can be sampled efficiently in the frequency domain, and the excursion probability can be fitted with these samples to extend any estimation of the look-elsewhere effect to lower p values. In addition, in cases where the Gaussian random field is stationary and the parameter space is Euclidean, the look-elsewhere effect correction can be computed analytically. We demonstrate these methods using two example template matching problems. Finally, we apply these methods to estimate the trial factor of a $$4^3$$ 4 3 accelerometer array for the detection of dark matter tracks in the Windchime project. When a global significance of $$3\sigma $$ 3 σ is required, the estimated trial factor for such an accelerometer array is $$10^{14}$$ 10 14 for a one-second search, and $$10^{22}$$ 10 22 for a 1-year search.

Qin, Juehang (ORCID:0000000182288949)↗
DOE OSTI2023

Code for the manuscript "Lagrangian Large Eddy Simulations vis Physics-informed Machine Learning"

We would like to disclose a code based on Python/Pytorch, in which we implement the data-driven physics-informed machine learning of Lagrangian large eddy simulation that is described in the manuscript "Lagrangian Large Eddy Simulations via Physics-Informed Machine Learning ". High-Reynolds number Homogeneous Isotropic Turbulence is fully described within the Navier-Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as Large Eddy Simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about the effects of sub-grid scales on the resolved scales. Here, we take an alternative approach and design novel LES heuristics stated in terms of Lagrangian particles moving with the flow. Our \emph{Lagrangian LES}, thus L-LES, is described by equations generalizing the weakly compressible Smoothed Particle Hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from Direct Numerical Simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The sub-grid scale contributions are modeled separately with physical constraints to account for the effects from un-resolved scales. We build the resulting model under the Differentiable Programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our Lagrangian LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

Tian, Yifeng↗
DOE OSTI2023

Secondary Pion Production using the LINAC Beam on Graphite and Tungsten Targets

The Fermilab 400-MeV Linac primary proton beam can produce precision, single-species secondary beams using a production target in the MeV Test Area experimental hall; specifically, pions from 4 up to 120 MeV kinetic energy. (These low energy pions rapidly decay into a muon beam in a secondary collection and transport beamline.) Although graphite is a standard production target material given its ultra-high temperature tolerance and low impact on the post-target primary beam, secondary pion/muon total production cross sections increase approximately as Z1/3 for positive muons and N2/3 for negative muons. Higher mass targets such as tungsten can potentially increase low-energy m+ and m- rates by factors of 3 and 8, respectively. Initial production studies have confirmed higher muon yield from heavy targets (tungsten vs carbon), but also, surprisingly, show significant differences in pion production between modern hadronic models (GENIEhad) and between GEANT and MCNP, state of the art production modeling codes. This project will further investigate the production discrepancies between these two target materials and between different physics models with high-statistics runs using the GEANT-based code, G4Beamline to better quantify the results and compare with production data which will be obtained using the MTA secondary beamline.

Ahmed, Shiza↗
NASA NTRS1977

Modeling of spacecraft charging

Three types of modeling of spacecraft charging are discussed: statistical models, parametric models, and physical models. Local time dependence of circuit upset for DoD and communication satellites, and electron current to a sphere with an assumed Debye potential distribution are presented. Four regions were involved in spacecraft charging: (1) undisturbed plasma, (2) plasma sheath region, (3) spacecraft surface, and (4) spacecraft equivalent circuit.

Whipple, E. C., Jr.↗
NASA NTRS1984

183 GHz water line variation: An energetic outburst in orion KL

Observations of the 3(13)-2(20) transition of water vapor in the direction of Ori MC1 in 1980 February show a 50% flux increase and an apparent additional red shift of approximately 2 km/s relative to the line observed in 1977 December. From a detailed examination of the amplitude and frequency calibration, it appears unlikely that the effect is due to systematic error. The increase is attributed to the appearance of a new component at a velocity of 12 km/s with respect to the local standard of rest. The new component also has broad wings. Increased emission from a region in the high-velocity core of Ori MC1 can be due either to additional far-IR radiation to pump the 1983 GHz transition or to a change in the physical conditions in the gas. Statistical equilibrium calculations using the large-velocity-gradient formalism were carried out to develop a model for the emission. The calculations support a model in which the gas in the region of enhanced emission is hotter than the dust. The temporal coincidence between the 183 GHZ increase and the 22 GH1 water maser outburst suggests a common, impulsive cause, which has heated the gas in a part of the HV source, enhancing the emission in both transitions.

Kuiper, T. B. H.↗
NASA NTRS1984

Convex profiles from asteroid lightcurves

A lightcurve inversion method that yields a two-dimensional convex profile is introduced. The number of parameters that characterize the profile is limited only by the number of Fourier harmonics used to represent the parent lightcurve. The implementation of the method is outlined by a recursive quadratic programming algorithm, and its application to photoelectric lightcurves and radar measurements is discussed. Special properties of the lightcurves of geometrically scattering ellipsoids are pointed out, and those properties are used to test the inversion method and obtain a criterion for judging whether any lightcurve could actually be due to such an object. Convex profiles for several asteroids are shown, and the method's validity is discussed from a physical as well as purely statistical point of view.

Ostro, S. J.↗
NASA NTRS1984

183 GHz water line variation - An energetic outburst in Orion KL

Observations of the 3(13)-2(20) transition of water vapor in the direction of Ori MC1 in 1980 February show a 50 percent flux increase and an apparent additional red shift of approximately 2 km/s relative to the line observed in 1977 December. From a detailed examination of the amplitude and frequency calibration, it appears unlikely that the effect is due to systematic error. The increase is attributed to the appearance of a new component at a velocity of 12 km/s with respect to the local standard of rest. The new component also has broad wings. Increased emission from a region in the high-velocity core of Ori MC1 can be due either to additional far-IR radiation to pump the 1983 GHz transition or to a change in the physical conditions in the gas. Statistical equilibrium calculations using the large-velocity-gradient formalism were carried out to develop a model for the emission. The calculations support a model in which the gas in the region of enhanced emission is hotter than the dust. The temporal coincidence between the 183 GHz increase and the 22 GH1 water maser outburst suggests a common, impulsive cause, which has heated the gas in a part of the HV source, enhancing the emission in both transitions.

Kuiper, T. B. H.↗
NASA NTRS1985

Convex-profile Inversion of Asteroid Lightcurves

A lightcurve inversion method that yields a two-dimensional convex profile is introduced. The number of parameters that characterize the profile is limited only by the number of Fourier harmonics used to represent the parent lightcurve. The implementation of the method is outlined by a recursive quadratic programming algorithm, and its application to photoelectric lightcurves and radar measurements is discussed. Special properties of the lightcurves of geometrically scattering ellipsoids are pointed out, and those properties are used to test the inversion method and obtained a criterion for judging whether any lightcurve could actually be due to such an object. Convex profiles for several asteroids are shown, and the method's validity is discussed from a physical as well as purely statistical point of view.

Ostro, S. J.↗
NASA NTRS1987

Composite load spectra for select space propulsion structural components

The objective of this program is to develop generic load models to simulate the composite load spectra (CLS) that are induced in space propulsion system components representative of the space shuttle main engines (SSME). These models are being developed through describing individual component loads with an appropriate mix of deterministic and state-of-the-art probabilistic models that are related to key generic variables. Combinations of the individual loads are used to synthesize the composite loads spectra. A second approach for developing the composite loads spectra load model simulation, the option portion of the contract will develop coupled models which combine the individual load models. Statistically varying coefficients of the physical models will be used to obtain the composite load spectra.

Newell, J. F.↗
NASA NTRS1986

Discrete Dirac equation on a finite half-integer lattice

The formulation of the Dirac equation on a discrete lattice with half-integer spacing and periodic boundary conditions is investigated analytically. The importance of lattice formulations for problems in field theory and quantum mechanics is explained; the concept of half-integer Fourier representation is introduced; the discrete Dirac equation for the two-dimensional case is derived; dispersion relations for the four-dimensional case are developed; and the spinor formulation for the Dirac fields on the half-integer lattice and the discrete time variable for the four-dimensional time-dependent Dirac equation are obtained. It is argued that the half-integer lattice, because it takes the Dirac Lagrangian into account, is more than a mere relabeling of the integer lattice and may have fundamental physical meaning (e.g., for the statistics of fermions). It is noted that the present formulation does not lead to species doubling, except in the continuum limit.

Smalley, L. L.↗
NASA NTRS1990

Direct simulations of wall-bounded compressible turbulence

A study has been initiated on the effects of compressibility (Mach number) on turbulent boundary layers. An understanding of both qualitative (turbulence structures, physical processes) and quantitative (turbulence statistics) effects are desired. This understanding should lead to better turbulence models for applications involving supersonic wall-bounded flows. Direct numerical simulations of an idealized problem will be used to accomplish these objectives. Among several possibilities we chose plane Couette flow with constant-temperature walls as the first problem to be studied. The lack of mean streamwise pressure gradient plus isothermal walls implies that both horizontal directions can be assumed to be homogeneous and that the flow can reach a statistically steady state. Together, these features greatly simplify the calculations and analyses of the results. To date, an algorithm has been developed and implemented (but not yet fully tested) for the accurate solution of the Navier-Stokes equations with the assumptions noted above.

Buell, J. C.↗
NASA NTRS1994

Multivariate space - time analysis of PRE-STORM precipitation

This paper presents the methodologies and results of the multivariate modeling and two-dimensional spectral and correlation analysis of PRE-STORM rainfall gauge data. Estimated parameters of the models for the specific spatial averages clearly indicate the eastward and southeastward wave propagation of rainfall fluctuations. A relationship between the coefficients of the diffusion equation and the parameters of the stochastic model of rainfall fluctuations is derived that leads directly to the exclusive use of rainfall data to estimate advection speed (about 12 m/s) as well as other coefficients of the diffusion equation of the corresponding fields. The statistical methodology developed here can be used for confirmation of physical models by comparison of the corresponding second-moment statistics of the observed and simulated data, for generating multiple samples of any size, for solving the inverse problem of the hydrodynamic equations, and for application in some other areas of meteorological and climatological data analysis and modeling.

Polyak, Ilya↗
NASA NTRS2011

Uncertainty Assessment of Hypersonic Aerothermodynamics Prediction Capability

The present paper provides the background of a focused effort to assess uncertainties in predictions of heat flux and pressure in hypersonic flight (airbreathing or atmospheric entry) using state-of-the-art aerothermodynamics codes. The assessment is performed for four mission relevant problems: (1) shock turbulent boundary layer interaction on a compression corner, (2) shock turbulent boundary layer interaction due a impinging shock, (3) high-mass Mars entry and aerocapture, and (4) high speed return to Earth. A validation based uncertainty assessment approach with reliance on subject matter expertise is used. A code verification exercise with code-to-code comparisons and comparisons against well established correlations is also included in this effort. A thorough review of the literature in search of validation experiments is performed, which identified a scarcity of ground based validation experiments at hypersonic conditions. In particular, a shortage of useable experimental data at flight like enthalpies and Reynolds numbers is found. The uncertainty was quantified using metrics that measured discrepancy between model predictions and experimental data. The discrepancy data is statistically analyzed and investigated for physics based trends in order to define a meaningful quantified uncertainty. The detailed uncertainty assessment of each mission relevant problem is found in the four companion papers.

Bose, Deepak↗
NASA NTRS2004

The Elegance of Disordered Granular Packings: A Validation of Edwards' Hypothesis

We have found a way to analyze Edwards' density of states for static granular packings in the special case of round, rigid, frictionless grains assuming constant coordination number. It obtains the most entropic density of single grain states, which predicts several observables including the distribution of contact forces. We compare these results against empirical data obtained in dynamic simulations of granular packings. The agreement between theory and the empirics is quite good, helping validate the use of statistical mechanics methods in granular physics. The differences between theory and empirics are mainly due to the variable coordination number, and when the empirical data are sorted by that number we obtain several insights that suggest an underlying elegance in the density of states

Metzger, Philip T.↗
NASA NTRS2019

Lightning Strike Distance Distribution Beyond a Preexisting Lightning Area

The 45th Weather Squadron (45 WS) asked the Applied Meteorology Unit (AMU) to review the 30-year-old, lightning stand-off distances of 5 nautical miles (nmi) for applicability to today's operations. This was based on the realization that previous lightning strike distance studies did not match how 45 WS issues lightning warnings (Roeder, 2008). The previous lightning distance studies were from the point of origin of the lightning or the average starting location that would tend to be in the core of the thunderstorm. However, the 45 WS issues lightning warnings based on the edge of a preexisting lightning area. Before beginning the AMU project, it took several years to develop a method to calculate a distance distribution beyond a preexisting area (Roeder, 2015). The AMU pulled Lightning Detection and Ranging (LDAR) sensor data from 1/1/2013 to 12/31/2013. This dataset consisted of 37 million individual source data points from the LDAR sensors. Only sources within 50 km north, south, east or west of the LDAR grid center were included in the dataset. This limited the use of LDAR data to that with the greatest accuracy of source detection and increased data processing. Points were grouped into flashes based on spatial and temporal criteria. Based on the sensitivity analysis the AMU performed on the flash clustering algorithm, a time value of 0.3 seconds was found to model flashes adequately. Distance parameters were tested from 1,500 to 7,500 meters (m) in 500 m increments. Distance parameters of both 3,000 m and 4,000 m produced results in the plotting tool that were most representative of the physical behavior of lightning. Thus statistics were gathered for the most representative of these spatial and temporal criteria on the flash size and the polygon expansion distance in order to find the correct data distributions. The best fit curves for the LDAR polygon expansion frequency vs. distances for both the 3 kilometer (km) and 4 km distance threshold values were exponential decay functions and had R2 values of > 0.998, indicating good model fits. The equations of the best fit curves were then used to calculate a desired safety radius of 4 nmi for either 3 km or 4 km distance threshold criteria. The AMU analysis concludes the safe reduction of the 5 nmi lightning warning circles to 4 nmi should improve the operational impact by 36% if based on distance from the center of the property area being protected. If based on the edge of the property being protected, then the reduction is 4.5 nmi to 4 nmi and the operational impact is decreased by 21%. For the 6 nmi lightning warning circles, the recommended 4 nmi stand-off distance will result in a safe reduction of operational impact of 31% if based on the center of the area being protected, or 16% if based on the edge of the property.

Flash clustering algorithm↗
DOE OSTI2024

Size Up or Size Down? National Analysis of Heat Pump Sizing and Impacts

Electrification of on-site fossil fuel combustion in buildings is recognized as a key component of achieving global greenhouse gas emissions targets. Air-source heat pumps are an efficient approach for space heating electrification. However, a major barrier to heat pump adoption is the high installation cost relative to furnaces and air conditioners. Because the installation cost of heat pumps - especially cold climate models - is dependent on their size, it is important to understand the distribution of heat pump sizes using different sizing methods and factors that impact heat pump sizing. In this study, we use sub-hourly physics simulations of 550,000 statistically representative dwelling units to analyze the distribution of heat pump size in the U.S. for three different efficiency levels of heat pumps, with and without building insulation and air sealing upgrades. We will present results from the national analysis of different factors affecting heat pump sizing along with the impact of sizing decisions on upgrade costs and operating costs.

electrification↗