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 271 records · Page 15

Virtual refrigerant charge sensing algorithm for residential CO₂ heat pumps

Natural refrigerants are increasingly adopted in next-generation heat pump systems, among which CO₂ heat pumps have attracted significant attention. However, due to their high operating pressures, the leakage risk is higher, resulting in undercharge conditions and degraded heat pump performance. Thus, developing an accurate refrigerant charge level detection technique is necessary to guarantee safe and efficient operation. Although virtual refrigerant charge (VRC) level calculation algorithms for CO₂ heat pumps exist, they typically rely on empirically selected features without a systematic selection framework, leading to multicollinearity and potential overfitting, which limit their prediction accuracy and generalizability. To address these issues, this study proposes a VRC algorithm framework with a systematic feature selection method that identifies physically meaningful and statistically significant features, and is applied using a residential CO₂ heat pump as a case study. The method is extended from previous work on conventional refrigerants to account for charge behavior in CO₂ gas coolers. The selected features include gas cooler outlet density, evaporator pressure, and superheat temperature. The results demonstrate that the proposed feature selection method significantly improves prediction accuracy compared to existing VRC approaches. A relatively small training dataset (∼30 samples) is sufficient for feature identification and model development. The developed algorithm achieves less than 3% prediction error under both undercharge and overcharge conditions, representing reductions of 46.7% and 35.3% compared to two recent reference VRC algorithms for transcritical CO₂ heat pumps reported in the literature. The proposed algorithm and feature selection method enhance leakage detection capability, facilitate the deployment of CO₂ heat pump systems, and contribute to reduced energy waste and maintenance costs.

Guo, Fangzhou [Lawrence Berkeley National Laborato↗

A better way to define dark matter haloes

ABSTRACT Dark matter haloes have long been recognized as one of the fundamental building blocks of large-scale structure formation models. Despite their importance – or perhaps because of it! – halo definitions continue to evolve towards more physically motivated criteria. Here, we propose a new definition that is physically motivated, effectively unique, and parameter-free: ‘A dark matter halo is comprised of the collection of particles orbiting in their own self-generated potential’. This definition is enabled by the fact that, even with as few as ≈300 particles per halo, nearly every particle in the vicinity of a halo can be uniquely classified as either orbiting or infalling based on its dynamical history. For brevity, we refer to haloes selected in this way as physical haloes. We demonstrate that (1) the mass function of physical haloes is Press–Schechter, provided the critical threshold for collapse is allowed to vary slowly with peak height; and (2) the peak-background split prediction of the clustering amplitude of physical haloes is statistically consistent with the simulation data, with accuracy no worse than ≈5 per cent.

Astronomy & Astrophysics↗

Toward Robust and Routine Determination of Mw for Small Earthquakes: Application to the 2020 Mw 5.7 Magna, Utah, Seismic Sequence

To better characterize seismic hazard, particularly, for induced seismicity, there is an increasing interest in methods to estimate moment magnitude (M w ) for small earthquakes. M w is generally preferred over other magnitude types, but, it is difficult to estimate M w for earthquakes with local magnitude (M L ) <3-3.5, using conventional moment tensor (MT) inversion. The 2020 M ww 5.7 Magna, Utah, seismic sequence provides an opportunity to illustrate and evaluate the value of spectral methods for this purpose. Starting with a high-quality seismic catalog of 2103 earthquakes (M L <5.6), we estimate M w using two independent spectral methods—one based on direct waves, yielding M w,direct , and the other based on coda waves, yielding M w,coda . For the direct-wave method, we present a non-parametric (NP) inversion scheme that solves for apparent geometrical spreading, G(R), and site effects (S), similar to other NP procedures that have been used to calibrate regional M L scales. The NP inversion is constrained using M ws derived from MTs for nine events in the Magna sequence. We recover statistically robust and physically reasonable G(R) and S and compute M w,direct for 635 Magna earthquakes down to M L 0.7. For the coda-wave method, we consider two separate calibration schemes involving previous MT solutions and compute M w,coda for 311 earthquakes down to M L 1.0. For 280 of the events that were processed with both methods—M w,direct and M w,coda —are strongly correlated (r = 0.98), with a mean difference of only 0.05. We compare M w,direct and M w,coda with M L and find reasonably good agreement for M L <3.6 with the theoretically predicted relationship of M w =(2/3)M L +C, in which C is a regional constant. Our results imply that seismic network operators can use spectral-based M w estimates to replace M L estimates for events with M L ≥1.0, and possibly smaller. The main requirement is the existence of a small number of MT solutions for calibration purposes.

58 GEOSCIENCES↗

Phasor-Measurement-Unit-Based Data Analytics Using Digital Twin and PhasorAnalytics Software

A major objective of this project was to apply GE’s commercial machine learning and data analytics toolsets to large-scale, real-world, anonymized Phasor Measurement Unit (PMU) datasets in order to extract signatures, correlated and/or causal factors, and precursor patterns associated with significant power system phenomena. The project had a particular emphasis on extraction of insights relevant to asset health monitoring, real-time load modeling and cybersecurity monitoring. Additionally, the team was directed to undertake a comprehensive data quality analysis for the provided datasets and encouraged to estimate the ‘machine-learning readiness’ of the datasets by documenting any major obstacles to the application of commercial machine learning algorithms. To accomplish the aforementioned objectives, the project team’s work centered around the identification of key event signatures and application of the identified event signatures for event detection and event classification. The industry-validated, semi-supervised machine learning strategy employed for event signature identification involved several major tasks, including data-preprocessing, generation of an overabundance of features, normal data identification, normality modeling, and event signature identification through a methodical, quantitative ranking of features in order of relevance to each studied event type. Throughout the project, data quality issues and mitigation techniques were investigated. In this report, insights are provided regarding the readiness of the provided synchrophasor datasets for application of machine learning and data analytics. The methodologies employed for this technical strategy are summarized in this report. With regards to data preprocessing and feature generation, the provided Training and Test Datasets were ingested into GE’s big data environment. Subsequently, the team applied bad data cleansing and data imputation scripts, event detection scripts, and application programming interfaces (APIs) to the datasets for convenient data access. The project team completed development and validation of dozens of physics-based, statistics-based and transformation-based feature functions used for the extraction of over 60 synchrophasor features. Using a new parallel feature generation technology developed on this project, over 60 features have been rapidly generated for the full two years’ worth of Training and Test Dataset data associated with both the Eastern and Western interconnects. Even accommodating for temporal down-sampling inherent to the feature extraction procedure, this parallel feature generation activity resulted in a massive feature set with a storage requirement approximately equal to that of the raw training dataset itself. With regards to normal data identification and normality modeling, a normality model was built using the feature data extracted from the Training Dataset and iteratively refined subsequent to incremental adjustments and expansions of the Training Dataset feature data. With respect to event characterization and signature identification, an event signature identification pipeline was developed and used in conjunction with the normality model to identify over 15 event signatures for key event categories within the Training Dataset. The identified event signatures were used to characterize hundreds of key events in terms of relative severity, duration, and location of the event. An investigation was undertaken to identify correlated and causal factors involved in transformer events. A separate investigation into temporal trends in ring-down analysis results was undertaken to determine possible associations between system dynamics and various other factors such as loading, season or year. To validate the identified event signatures, additional work was undertaken to develop signature-based anomaly detection and classification tools suitable for convenient application to the synchrophasor datasets. The anomaly detection and classification tools, suitable for online application, were then applied to the entirety of the Eastern Interconnect Training and Test Datasets. Performance of the event detection and classification tools was evaluated upon receipt of the Test Dataset event logs (i.e., the labels for events contained in the Test Dataset), and promising results were obtained despite several challenges (documented herein) associated with application of supervised or semi-supervised machine learning methods to large-scale, anonymized datasets. Finally, the detection and classification tools were used to detect, classify, and characterize thousands of new events not included in the original event logs provided by the DOE within both the Training and Test Datasets.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Highly Resolved Reference Projections of Building Energy Use for the Contiguous United States: Building Sector Energy Baselines, Projection Methods, and Results

This report describes one methodology of projecting energy consumption of the US residential and commercial building sectors using NREL's ResStock™ and ComStock™ as well as growth rates derived from EIA's Annual Energy Outlook (AEO). The impetus for this work is to provide an intermediate method for compiling demand-side sectoral energy projections that is suitable for grid-scale analysis, such as NREL's Standard Scenarios. ResStock and ComStock are physics-based and statistically representative building stock models of the US residential and commercial sector, respectively. Using the 2012 actual meteorological year (AMY) weather data, the sectoral energy baselines are simulated and then segmented along key dimensions (e.g., geography, dwelling/building type). The segmented results are then scaled using the corresponding annual growth rates derived from the 2021 AEO reference case to produce energy projections out to 2050. The compiled result is a demand-side grid model (dsgrid) data set suitable for use in NREL's large-scale grid models, such as the Regional Energy Deployment System (ReEDS). This simple projection method does not endogenously represent how the building stock could evolve through time. Most notably, it does not reflect large-scale electrification, for example, the conversion of space heating, water heating, clothes drying, and cooking from primary fossil fuels to electricity, as this is not part of AEO's reference case assumptions. Nonetheless this approach is more resolved and potentially extensible compared to the current method used by Standard Scenarios's reference case, which augments a sector's total load based on a single growth rate from AEO.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Improving Protein–Ligand Interaction Modeling with cryo-EM Data, Templates, and Deep Learning in 2021 Ligand Model Challenge

Elucidating protein–ligand interaction is crucial for studying the function of proteins and compounds in an organism and critical for drug discovery and design. The problem of protein–ligand interaction is traditionally tackled by molecular docking and simulation, which is based on physical forces and statistical potentials and cannot effectively leverage cryo-EM data and existing protein structural information in the protein–ligand modeling process. In this work, we developed a deep learning bioinformatics pipeline (DeepProLigand) to predict protein–ligand interactions from cryo-EM density maps of proteins and ligands. DeepProLigand first uses a deep learning method to predict the structure of proteins from cryo-EM maps, which is averaged with a reference (template) structure of the proteins to produce a combined structure to add ligands. The ligands are then identified and added into the structure to generate a protein–ligand complex structure, which is further refined. The method based on the deep learning prediction and template-based modeling was blindly tested in the 2021 EMDataResource Ligand Challenge and was ranked first in fitting ligands to cryo-EM density maps. These results demonstrate that the deep learning bioinformatics approach is a promising direction for modeling protein–ligand interactions on cryo-EM data using prior structural information.

59 BASIC BIOLOGICAL SCIENCES↗

A Hybrid Deep Learning Approach to Cosmological Constraints from Galaxy Redshift Surveys

We present a deep machine learning (ML)–based technique for accurately determining σ g and Ω m from mock 3D galaxy surveys. The mock surveys are built from the AbacusCosmos suite of N -body simulations, which comprises 40 cosmological volume simulations spanning a range of cosmological parameter values, and we account for uncertainties in galaxy formation scenarios through the use of generalized halo occupation distributions (HODs). We explore a trio of ML models: a 3D convolutional neural network (CNN), a power spectrum–based fully connected network, and a hybrid approach that merges the two to combine physically motivated summary statistics with flexible CNNs. We describe best practices for training a deep model on a suite of matched-phase simulations, and we test our model on a completely independent sample that uses previously unseen initial conditions, cosmological parameters, and HOD parameters. Despite the fact that the mock observations are quite small (~0.07 h -3 Gpc 3 ) and the training data span a large parameter space (six cosmological and six HOD parameters), the CNN and hybrid CNN can constrain estimates of σ g and Ω m to ~3% and ~4%, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

Comets and the formation of planets

Morphological study of the physical and dynamical processes of planet formation, with emphasis on the role of comet nuclei. A consistent model proposes the formation of comets and planets in preplanetary rings of the residual solar nebula, with subsequent ejection, chiefly by Jupiter, of the comets to Oort's (1950) sphere. Physically, dynamically, or statistically evaluated items include: (1) the total number and mass of comets in Oort's cloud; (2) reevaluation of the diameters and masses of comet nuclei; (3) the processes of nucleation from gravitational and 'Boltzmann' instabilities of gaseous media to agglomerations of particulate matter; and (4) the statistical-dynamical conditions and time scales of orbital interaction of comets with the planets and the consequences of disintegration.

Opik, E. J.↗

Statistical model for premixed turbulent flames

Statistical concepts and physical arguments are used to develop closure models for the mean energy equation governing one-dimensional, premixed, intense turbulent flames. It is proposed that the turbulent transport terms can be modeled by a turbulent conductivity based on turbulence kinetic energy (k) and energy dissipation rate (epsilon). The models obtained are used to derive correlations for the turbulent flame speed. Results show that there are three different flame regimes: a low-intensity turbulence, wrinkled flame; a high-intensity, fast-chemistry flame; and a high-intensity, finite-rate-chemistry flame.

Gouldin, F. C.↗

Analysis of turbofan propulsion system weight and dimensions

Weight and dimensional relationships that are used in aircraft preliminary design studies are analyzed. These relationships are relatively simple to prove useful to the preliminary designer, but they are sufficiently detailed to provide meaningful design tradeoffs. All weight and dimensional relationships are developed from data bases of existing and conceptual turbofan engines. The total propulsion system is considered including both engine and nacelle, and all estimating relations stem from physical principles, not statistical correlations.

Waters, M. H.↗

Lagrange thermodynamic potential and intrinsic variables for He-3 He-4 dilute solutions

For a two-fluid model of dilute solutions of He-3 in liquid He-4, a thermodynamic potential is constructed that provides a Lagrangian for deriving equations of motion by a variational procedure. This Lagrangian is defined for uniform velocity fields as a (negative) Legendre transform of total internal energy, and its primary independent variables, together with their thermodynamic conjugates, are identified. Here, similarities between relations in classical physics and quantum statistical mechanics serve as a guide for developing an alternate expression for this function that reveals its character as the difference between apparent kinetic energy and intrinsic internal energy. When the He-3 concentration in the mixtures tends to zero, this expression reduces to Zilsel's formula for the Lagrangian for pure liquid He-4. An investigation of properties of the intrinsic internal energy leads to the introduction of intrinsic chemical potentials along with other intrinsic variables for the mixtures. Explicit formulas for these variables are derived for a noninteracting elementary excitation model of the fluid. Using these formulas and others also derived from quantum statistical mechanics, another equivalent expression for the Lagrangian is generated.

Jackson, H. W.↗

Classification, seasonality and persistence of low-frequency atmospheric circulation patterns

The seasonality and persistence of the major modes of interannual variability in the low-frequency atmospheric circulation were studied using orthogonally rotated principle component analysis (RPCA) of Northern Hemisphere 1-month mean 700 mb heights. The twice-daily data for the 1950-1984 period were used. Winter results are similar to those of other recent RPCA and teleconnection studies. The strongest summer pattern is the North Atlantic Oscillation, which is also the strongest winter pattern; it systematically contracts northward in summer and expands southward in winter, being the only pattern found for every month of the year. The robustness of the RPCA results was examined through consistency with results of other studies and of adjacent month solutions within this study, as well as by replicating the results using 3-month and 10-day means of 700-mb height. It is concluded that the RPCA method provides a physically meaningful and statistically stable product with the simplicity of teleconnection patterns but with superior pattern choice and depiction.

Barnston, Anthony G.↗

Relativistic beaming and the nuclei of double-lobed quasars

If the simple relativistic beaming model for extragalactic radio sources is at all correct, Doppler boosting may bias us toward choosing those objects pointing nearly at us in flux-density limited samples. In order to minimize any such effects in statistical tests of physical theories of these objects, an attempt is made to define a complete sample of sources with random orientations. For a sample of 26 double-lobed quasars, the dependence of a number of observable features on the angle to the line of sight and on the superluminal velocity is examined, and the distribution of superluminal velocities for the sources is measured. Two new superluminals are reported, and their apparent velocities and nuclear strengths are discussed.

Hough, David H.↗

Structures Of Turbulent Boundary Layers

Report presents second part of two-part comprehensive study of coherent motions in, and structures of, turbulent boundary-layer flows. Study critically reviews current fragmented knowledge with view toward eventual unification and understanding of physical causes and statistical properties of turbulence phenomena.

Robinson, Stephen K.↗

The GEOS Ozone Data Assimilation System: Design and Validation

An ozone data assimilation system has been developed at the Data Assimilation Office of the NASA/Goddard Space Flight Center to provide global three-dimensional analyzed ozone mixing ratio and total column ozone. The Total Ozone Mapping Spectrometer (TOMS) total column ozone and the Solar Backscatter Ultraviolet (SBUV) or SBUV/2 partial ozone profile data are assimilated. The analyzed winds from the Goddard Earth Observing System Data Assimilation System (GEOS-DAS) drive the ozone transport. Following every transport model timestep, the model prediction is combined with the observations using a global, physical-space based, statistical analysis scheme. Due to the smaller size of the ozone system than that of a global meteorological data assimilation system, new statistical analysis methodology, including anisotropic and flow-dependent forecast error correlation models, can be implemented and tested in the ozone system more easily. Sample results from the winter 1992 validation period are presented. There is a close agreement between the analyzed fields and the independent observations from ozone sondes and the Halogen Occultation Experiment (HALOE).

Stajner, I.↗

Atlas of Seasonal Means Simulated by the NSIPP 1 Atmospheric GCM

This atlas documents the climate characteristics of version 1 of the NASA Seasonal-to-Interannual Prediction Project (NSIPP) Atmospheric General Circulation Model (AGCM). The AGCM includes an interactive land model (the Mosaic scheme), and is part of the NSIPP coupled atmosphere-land-ocean model. The results presented here are based on a 20-year (December 1979-November 1999) "ANIIP-style" integration of the AGCM in which the monthly-mean sea-surface temperature and sea ice are specified from observations. The climate characteristics of the AGCM are compared with the National Centers for Environmental Prediction (NCEP) and the European Center for Medium-Range Weather Forecasting (ECMWF) reanalyses. Other verification data include Special Sensor Microwave/Imager (SSNM) total precipitable water, the Xie-Arkin estimates of precipitation, and Earth Radiation Budget Experiment (ERBE) measurements of short and long wave radiation. The atlas is organized by season. The basic quantities include seasonal mean global maps and zonal and vertical averages of circulation, variance/covariance statistics, and selected physics quantities.

Suarez, Max J.↗

Studying the Iron Line Complex in the Bright Seyfert Galaxy NGC 5506

This grant was to support the reduction and analysis of our approved XMM observation of the nearby Seyfert 2 galaxy NGC 5506. The observation has been carried out simultaneously with a BeppoSAX observation of the same source. The proposal was aimed to study in detail the Compton reflection component and the complex Iron K line of this source, combining the still unique capability of BeppoSAX in hard X-rays (to strongly constrain the reflection component, and then the intrinsic nuclear continuum), and the sensitivity of XMM at the energy of the Iron Line complex. NGC 5506 is one of the brightest AGN in hard X-rays and has been intensively studied in the past. GINGA detected the complex iron line as well as the reflection component. Both ASCA (spectroscopically) and Rossi-XTE (through variability analysis) suggested that the FeK line is complex, possibly made up of several distinct components. The centroid of the FeK complex in a subsequent BeppoSAX observation was bluer than the 6.4 keV energy of the relatively low-ionization iron Kalpha transition. NGC 5506 has been observed simultaneously by NewtonXMM and BeppoSAX on February 2-3 2001. we have reduced and analyzed both the NewtonXMM and the BeppoSAX data, and have written and published a paper on our results (appeared in Volume 377 (page 31) of A&A-Letters). Our main results can be summarized as follows: (a) we confirm that the FeK line is complex, and for the first time disentangle its components: we find that at least two components made up the FeK complex, one neutral and narrow, at 6.4 keV (rest energy), and another one either broader and highly ionized, at about 6.7 keV (rest frame), or, in turn, made up of two narrow and unresolved components from the He-like and the H-like ions of Fe; (b) the two possible solution for the high-ionization Fe-K component, are statistically indistinguishable. However, physically, a blend of two narrow lines from photoionized matter seems to be preferable to emission of a relativistically broadened line from an ionized accretion disk; (c) The bulk of the Compton reflection continuum originates in neutral matter, and is therefore associated with the narrow FeK line at 6.4 keV: they are most likely emitted by distant matter.

Nicastro, Fabrizio↗

Introduction to Data Assimilation

Atmospheric data assimilation is a class of techniques used for producing descriptions of fields of air temperature, pressure, humidity, wind, etc. on a spatial grid or in terms of a finite functional representation. These are then used to initialize numerical weather forecasts or to analyze the atmosphere for other purposes. The techniques combine past, present, and even future observations in an approximate statistically optimal way. Various types of statistical or physically-based models and their corresponding adjoints are employed to relate diverse fields in both time and space and to relate what is observed to what is being analyzed. Computationally, the problem is very demanding and onstraining on the techniques that can be employed on a routine basis.

Atmospheric data assimilation↗