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

Exploring scenarios for enhanced fuel compression and performance on the National Ignition Facility with machine-learning-aided design techniques

Recent fusion experiments on the National Ignition Facility (NIF) have achieved ignition, producing multi-MJ fusion yields for input laser energies of roughly 2 MJ [Abu-Shawareb et al., Phys. Rev. Lett. 132, 065102 (2024)]. Building on the success of the target designs that have achieved ignition, we explore new implosion scenarios predicted to generate significantly more compression of the dense DT ice layer and correspondingly higher yields while preserving many of the key physics characteristics of present-day ignition designs. Our main result is a novel 3-shock implosion scheme that effectively minimizes the shock-induced entropy in the dense, accelerating DT shell and maximizes the resulting fuel compression subject to a fixed leading shock strength consistent with present-day ignition experiments, which is necessary to melt the crystalline high-density carbon ablator. Compared to the first NIF experiment to fulfill Lawson's ignition criterion, shot N210808 [Abu-Shawareb et al., Phys. Rev. Lett. 129, 075001 (2022)], our design exhibits a 40% increase in simulated peak areal density (ρR) and a 5× increase in 1D fusion yield using a 4% lighter ablator and identical DT payloads. We also present a complete integrated 2D hohlraum design and laser pulse specifications capable of generating the desired 3-shock drive and maintaining control of the low-mode capsule implosion symmetry, where the increase in simulated 2D yield relative to N210808 is > 10×. This new implosion regime was discovered with help from a machine-learning-enabled capsule design optimization framework. We outline the workflow this automated tool uses to identify improved design candidates by running several rounds of capsule simulations, constructing a surrogate model mapping input variations to key physics output quantities, and querying the resulting statistical model to propose adjustments to the x-ray drive and capsule to reach a set of physics objectives prescribed by the designer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Physics-informed machine learning of the correlation functions in bulk fluids

The Ornstein–Zernike (OZ) equation is the fundamental equation for pair correlation function computations in the modern integral equation theory for liquids. In this work, machine learning models, notably physics-informed neural networks and physics-informed neural operator networks, are explored to solve the OZ equation. The physics-informed machine learning models demonstrate great accuracy and high efficiency in solving the forward and inverse OZ problems of various bulk fluids. The results highlight the significant potential of physics-informed machine learning for applications in thermodynamic state theory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Merger Response of Halo Anisotropy Properties

Anisotropy properties—halo spin, shape, position offset, velocity offset, and orientation—are an important family of dark matter halo properties that indicate the level of directional variation of the internal structures of halos. These properties reflect the dynamical state of halos, which in turn depends on the mass assembly history. In this work, we study the evolution of anisotropy properties in response to merger activity using the IllustrisTNG simulations. We find that the response trajectories of the anisotropy properties significantly deviate from secular evolution. These trajectories have the same qualitative features and timescales across a wide range of merger and host properties. We propose explanations for the behavior of these properties and connect their evolution to the relevant stages of merger dynamics. We measure the relevant dynamical timescales. We also explore the dependence of the strength of the response on time of merger, merger ratio, and mass of the main halo. These results provide insight into the physics of halo mergers and their effects on the statistical behavior of halo properties. This study paves the way toward a physical understanding of scaling relations, particularly to how systematics in their scatter are connected to the mass assembly histories of halos.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1983

Statistical correlation of structural mode shapes from test measurements and NASTRAN analytical values

The software and procedures of a system of programs used to generate a report of the statistical correlation between NASTRAN modal analysis results and physical tests results from modal surveys are described. Topics discussed include: a mathematical description of statistical correlation, a user's guide for generating a statistical correlation report, a programmer's guide describing the organization and functions of individual programs leading to a statistical correlation report, and a set of examples including complete listings of programs, and input and output data.

Purves, L.↗