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 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

A42F-03: Small-Scale Variability in Tropical Tropopause Layer Humidity

Recent advances in statistical parameterizations of cirrus cloud processes for use in global models are highlighting the need for information about small-scale fluctuations in upper tropospheric humidity and the physical processes that control the humidity variability. To address these issues, we have analyzed high-resolution airborne water vapor measurements obtained in the Airborne Tropical TRopopause EXperiment (ATTREX) over the tropical Pacific between 14 and 20 km. Using accurate and precise 1-Hz water vapor measurements along approximately-level aircraft flight legs, we calculate structure functions spanning horizontal scales ranging from about 0.2 to 50 km, and we compare the water vapor variability in the lower (about 14 km) and upper (16-19 km) Tropical Tropopause Layer (TTL). We also compare the magnitudes and scales of variability inside TTL cirrus versus in clear-sky regions. The measurements show that in the upper TTL, water vapor concentration variance is stronger inside cirrus than in clear-sky regions. Using simulations of TTL cirrus formation, we show that small variability in clear-sky humidity is amplified by the strong sensitivity of ice nucleation rate to supersaturation, which results in highly-structured clouds that subsequently drive variability in the water vapor field. In the lower TTL, humidity variability is correlated with recent detrainment from deep convection. The structure functions indicate approximately power-law scaling with spectral slopes ranging from about minus 5 divided by 3, to minus 2.

water vapor↗
NASA NTRSDate not supplied

Prospects for Cross-Correlations of UHECR Events with Astrophysical Sources with Upcoming Space-Based Experiments

Ultra-high energy cosmic rays (UHECRs) are the messengers of the most extreme physics in the cosmos; however, efforts to identify their origins have thus far been thwarted by the fact that they don’t point back to their sources. Using statistical studies cross-correlating UHECR arrival directions with astrophysical catalogs, the ground-based Pierre Auger Observatory has reported hints of a correlation with nearby starburst galaxies, as well as lower-significance correlations with other classes of astrophysical sources. Space-based UHECR experiments, such as POEMMA and ZAP, will monitor large interaction volumes on the Earth or the Moon. Within a few years of mission operation time, both missions will achieve unprecedented exposures at energies above40EeV across the entire sky. In this contribution, we present studies of the cross-correlation between UHECR event arrival directions and astrophysical catalogs as motivated by expectations for the detector performance for POEMMA and ZAP. We find that both POEMMA and ZAP will achieve 5σ discovery reach for many plausible astrophysical scenarios.

Tonia M Venters↗
NASA NTRS2022

The Earth Model Column Collaboratory (EMC2) v1.1: An Open-Source Ground-Based Lidar and Radar Instrument Simulator and Subcolumn Generator for Large-Scale Models

Climate models are essential for our comprehensive understanding of Earth's atmosphere and can provide critical insights on future changes decades ahead. Because of these critical roles, today's climate models are continuously being developed and evaluated using constraining observations and measurements obtained by satellites, airborne, and ground-based instruments. Instrument simulators can provide a bridge between the measured or retrieved quantities and their sampling in models and field observations while considering instrument sensitivity limitations. Here we present the Earth Model Column Collaboratory (EMC2), an open-source ground-based lidar and radar instrument simulator and subcolumn generator, specifically designed for large-scale models, in particular climate models, but also applicable to high-resolution model output. EMC2 provides a flexible framework enabling direct comparison of model output with ground-based observations, including generation of subcolumns that may statistically represent finer model spatial resolutions. In addition, EMC2 emulates ground-based (and air- or space-borne) measurements while remaining faithful to large-scale models' physical assumptions implemented in their cloud or radiation schemes. The simulator uses either single particle or bulk particle size distribution lookup tables, depending on the selected scheme approach, to perform the forward calculations. To facilitate model evaluation, EMC2 also includes three hydrometeor classification methods, namely, radar- and sounding-based cloud and precipitation detection and classification, lidar-based phase classification, and a Cloud Feedback Model Intercomparison Project Observational Simulator Package (COSP) lidar simulator emulator. The software is written in Python, is easy to use, and can be straightforwardly customized for different models, radars, and lidars. Following the description of the logic, functionality, features, and software structure of EMC2, we present a case study of highly supercooled mixed-phase cloud based on measurements from the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment (AWARE). We compare observations with the application of EMC2 to outputs from four configurations of the NASA Goddard Institute for Space Studies (GISS) climate model (ModelE3) in single-column model (SCM) mode and from a large-eddy simulation (LES) model. We show that two of the four ModelE3 configurations can form and maintain highly supercooled precipitating cloud for several hours, consistent with observations and LES. While our focus is on one of these ModelE3 configurations, which performed slightly better in this case study, both of these configurations and the LES results post-processed with EMC2 generally provide reasonable agreement with observed lidar and radar variables. As briefly demonstrated here, EMC2 can provide a lightweight and flexible framework for comparing the results of both large-scale and high-resolution models directly with observations, with relatively little overhead and multiple options for achieving consistency with model microphysical or radiation scheme physics.

Earth Model Column Collaboratory↗
DOE OSTI2022

Grassmannian Shape Representations for Aerodynamic Applications: Preprint

Airfoil shape design is a classical problem in engineering, science, and manufacturing. Our motivation is to combine principled physics-based considerations for the shape design problem with modern computational techniques informed by a data-driven approach. Traditional analyses of airfoil shapes emphasize a flow-based sensitivity to deformations which can be represented generally by affine transformations (rotation, scaling, shearing, shifting). We present a novel representation of shapes which decouples affine-style deformations from a rich set of data-driven deformations over a submanifold of the Grassmannian. The Grassmannian representation, informed by a database of physically relevant airfoils, offers (i) a rich set of novel 2D airfoil deformations not previously captured in the data, (ii) improved low-dimensional parameter domain for inferential statistics, and (iii) consistent 3D blade representation and perturbation over a sequence of nominal shapes.

blade representation↗
DOE OSTI2022

Grassmannian Shape Representations for Aerodynamic Applications

Airfoil shape design is a classical problem in engineering, science, and manufacturing. Our motivation is to combine principled physics-based considerations for the shape design problem with modern computational techniques informed by a data-driven approach. Traditional analyses of airfoil shapes emphasize a flow-based sensitivity to deformations which can be represented generally by affine transformations (rotation, scaling, shearing, shifting). We present a novel representation of shapes which decouples affine-style deformations from a rich set of data-driven deformations over a submanifold of the Grassmannian. The Grassmannian representation, informed by a database of physically relevant airfoils, offers (i) a rich set of novel 2D airfoil deformations not previously captured in the data, (ii) improved low-dimensional parameter domain for inferential statistics, and (iii) consistent 3D blade representation and perturbation over a sequence of nominal shapes.

blade representation↗
DOE OSTI2026

Measuring the neutrino-oxygen neutral current quasielastic cross section using the accelerator neutrino neutron interaction experiment

The Accelerator Neutrino Neutron Interaction Experiment (ANNIE) is a 26-ton gadolinium-doped water Cherenkov detector located on-axis to Fermilab’s Booster Neutrino Beam (BNB). ANNIE is uniquely positioned to perform high-statistics measurements of neutrino-nucleus interactions in water, benefiting from a large neutrino flux due to a short (100-meter) baseline. A central focus of ANNIE’s physics program is the measurement of both charged current (CC) and neutral current (NC) cross sections on water, including neutral current quasielastic (NCQE) and CC-inclusive channels. The NCQE measurement is particularly critical for constraining uncertainties in rare-event searches such as the Diffuse Supernova Neutrino Background (DSNB), where atmospheric $\nu$NCQE interactions constitute a significant and poorly constrained background. This dissertation presents a measurement of the flux-averaged neutrino-oxygen neutral current quasielastic ($\nu$NCQE) cross section using $2.573 \times 10^{20}$~POT of BNB exposure from the 2022 and 2023 beam years. The $\nu$NCQE interaction is identified through the primary $\gamma$-rays produced by nuclear de-excitation of the residual $^{15}$N$^*$ or $^{15}$O$^*$ nucleus following nucleon knockout from $^{16}$O. A dedicated Monte Carlo (MC) re-tuning campaign was conducted using an americium-beryllium (AmBe) calibration source, Michel electrons from stopped muons, and throughgoing dirt muons originating upstream of the detector. This multi-sample approach provided a wide-ranging $\mathcal{O}(\text{MeV})$--$\mathcal{O}(\text{GeV})$ dataset for tuning the simulated detector response, which was subsequently validated against AmBe neutron and Michel electron data for use in the $\nu$NCQE analysis. A dedicated laser calibration campaign was carried out to reduce timing uncertainties across the PMT system, enabling reconstruction of the BNB bunch substructure with sufficient resolution to serve as a background rejection tool. By selecting events in-time with individual neutrino bunches, beam-correlated $\nu$NCQE events are separated from diffuse and accelerator-induced backgrounds, notably skyshine neutrons and externally-originating events, that would otherwise dominate traditional charge-based selections within a small-scale, surface-level, short-baseline detector. A data-driven estimation of the skyshine neutron and external background rates was performed and incorporated into the systematic uncertainty budget. The flux-averaged $\nu$NCQE cross section on oxygen is measured to be $1.57 \pm 0.06\,(\text{stat.})$ $^{+0.91}_{-0.67}\,(\text{syst.})$ $\times 10^{-38}\ \text{cm}^{2}$. A full systematic budget is constructed by propagating uncertainties in the secondary hadronic interaction modeling, background cross section normalizations, detector response, neutrino flux, and the primary $\gamma$-ray emission probabilities from oxygen nuclear de-excitation. An idealized de-excitation model, constructed from existing measurements in the literature is developed to benchmark the predictions of the \textsc{GENIE} event generator. A comparison reveals that \textsc{GENIE} systematically overpredicts the primary $\gamma$-ray emission probability from oxygen de-excitation by a factor of $1.49\times$ for $E_\gamma > 6$~MeV and $3.07\times$ in the $3$--$6$~MeV band. This comparison motivates the dominant systematic uncertainty in this analysis, where a conservative uncertainty of $^{+39.9\%}_{-0\%}$ on the primary $\gamma$-ray signal prediction is assigned. The ANNIE result is consistent with and complementary to existing flux-averaged $\nu$NCQE cross section measurements from T2K and Super-Kamiokande, providing an independent measurement with a different detector, neutrino beam, and analysis methodology. Looking ahead, an upgrade to the ANNIE DAQ infrastructure enabling continuous extended readout will allow a complementary $\nu$NCQE neutron multiplicity measurement, directly relevant to constraining the NCQE background in DSNB searches, competitive with the recent T2K measurement at SK-Gd. The planned Super-SANDI upgrade, deploying a large Water-based Liquid Scintillator (WbLS) volume, will further extend ANNIE's reach to hadronic final states and exclusive NC channels, and enable joint measurements with liquid argon detectors sharing the BNB beamline ahead of DUNE and Hyper-Kamiokande.

Doran, Steven [Iowa State U.]↗
DOE OSTI2025

Hadronic light-by-light contribution to the muon anomaly from lattice QCD with infinite volume QED at physical pion mass

The hadronic light-by-light scattering contribution to the muon anomalous magnetic moment, (g–2)⁢/2, is computed in the infinite volume QED framework with lattice QCD. We report $a^{HLbL}_μ$ = 12.47⁢(1.15)⁢(0.95) ×10 –10 where the first error is statistical and the second systematic. The result is mainly based on the 2+1 flavor Möbius domain wall fermion ensemble with inverse lattice spacing a –1 = 1.73 GeV, lattice size L = 5.5 fm, and m π = 139 MeV, generated by the RBC-UKQCD collaborations. The leading systematic error of this result comes from the lattice discretization. This result is consistent with previous determinations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Interpreting the spatial distribution of soil properties with a physically-based distributed hydrological model

Digital soil maps are commonly data-driven as the development of physically-based models for soil mapping is difficult due to the complexity of soils. However, physically-based hydrologic models have been successful in simulating water dynamics. Since water movement is a major driver of pedogenesis, the physical rules that govern water movement might help explain and predict the spatial variation of soil properties. Here, we demonstrate the novel use of a physically-based, distributed hydrologic model to inform the spatial distribution of soil properties. The Distributed Hydrology Soil Vegetation Model (DHSVM) was utilized to simulate soil moisture content (SM) and water table depth (WTD) in two hillslope catchments under pasture and forest management wherein hydrologic model outputs were then compared with soil properties measured in situ. SM sensors and wells were installed in both catchments to validate simulations of soil water movement via Nash-Sutcliffe Efficiency (E). In-situ observations were made at 87 sites within both catchments to study the connection between simulated water movement (SM and WTD) and observed soil properties, namely the depth and thickness of the argillic (Bt), fragic (Btx), and C horizons, and the depth of redoximorphic features. The simulated time series of SM and WTD were also clustered per season using Dynamic Time Warping (DTW), which identified similarity among time series at varying timescales. Model validation suggested that simulations of surficial SM (0–20 cm) were reasonable (E = 0.45), however, simulated subsurface SM (45–60 cm) and WTD were not sufficiently accurate. The thickness of Btx horizons were spatially grouped into different populations by SM clusters from every season except spring. For the other properties, only SM dynamics of specific seasons grouped into significantly different populations, suggesting that the explanatory power of simulated water movement varies seasonally and was greater during winter. Here, we show clusters of simulated SM separated soil properties into statistically different populations, showing that hydrologic models could inform areas that followed different water dynamics related to pedogenic trajectories and related biogeochemical processes not necessarily simulated by the model. As such, physically-based modeling of water dynamics can, therefore, inform and advance digital soil mapping by linking water movement patterns stemming from hydrologic model outputs to spatial patterns of soil properties and pedogenesis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

State densities of heavy nuclei in the static-path plus random-phase approximation

We report that nuclear state densities are important inputs to statistical models of compound-nucleus reactions. State densities are often calculated with self-consistent mean-field approximations that do not include important correlations and must be augmented with empirical collective enhancement factors. Here, we benchmark the static-path plus random-phase approximation (SPA + RPA) to the state density in a chain of samarium isotopes 148–155 Sm against exact results (up to statistical errors) obtained with the shell-model Monte Carlo (SMMC) method. The SPA + RPA method incorporates all static fluctuations beyond the mean field together with small-amplitude quantal fluctuations around each static fluctuation. Using a pairing plus quadrupole interaction, we show that the SPA + RPA state densities agree well with the exact SMMC densities for both the even- and odd-mass isotopes. For the even-mass isotopes, we also compare our results with mean-field state densities calculated with the finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation. We find that the SPA + RPA repairs the deficiencies of the mean-field approximation associated with broken rotational symmetry in deformed nuclei and with the violation of particle-number conservation in the pairing condensate. In particular, in deformed nuclei the SPA + RPA reproduces the rotational enhancement of the state density relative to the mean-field state density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Exact quantum scars in the chiral nonlinear Luttinger liquid

While the chiral linear Luttinger liquid is integrable via bosonization, its nonlinear counterpart does not admit for an analytic solution. In this work, we find a subextensive number of exact eigenstates for a large family of density-density interaction terms. These states are embedded in a continuum of strongly correlated excited states. The real-space entanglement entropy of some exact states scales logarithmically with system size while that of others has volume-law scaling. We introduce momentum-space entanglement as an unambiguous differentiator between these exact states and the remaining excited states. With regard to momentum space, the exact states behave as bona fide quantum many-body scars: they exhibit identically zero momentum-space entanglement, while typical eigenstates behave thermally. We corroborate this finding by a level statistics analysis. Here we detail the general formalism for systematically finding all interaction terms and associated exact states, and present a number of infinite exact state sequences extending to arbitrarily high energies. Unlike many previous examples of quantum many-body scars, the exact states uncovered here do not lie at equidistant energies and do not follow from a special operator algebra. Instead, they are uniquely enabled by the interplay of Fermi statistics and chirality.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Quantal diffusion description of isotope production via the multinucleon transfer mechanism in 48 Ca+ 238 U collisions

As an extension of previous work, we calculate the production cross section of heavy neutron-rich isotopes by employing the quantal diffusion description to 48 Ca+ 238 U collisions. The quantal diffusion is deduced from stochastic mean-field approach, and transport properties are determined in terms of time-dependent single-particle wave functions of the time-dependent Hartree-Fock theory. As a result, the approach allows for prediction of production cross sections without any adjustable parameters. Furthermore, the secondary cross sections by particle emission are calculated with the help of the statistical gemini++ code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Observation of 𝑡⁢𝑊⁢𝑍 Production at the CMS Experiment

The first observation of single top quark production in association with a 𝑊 and a 𝑍 boson in proton-proton collisions is reported. The analysis uses data at center-of-mass energies of 13 and 13.6 TeV recorded with the CMS detector at the CERN LHC, corresponding to a total integrated luminosity of 200 fb −1 . Events with three or four charged leptons, which can be electrons or muons, are selected. Advanced machine-learning algorithms and improved reconstruction methods, compared to an earlier analysis, result in an unprecedented sensitivity to 𝑡⁢𝑊⁢𝑍 production. The measured cross sections for 𝑡⁢𝑊⁢𝑍 production are 248 ± 52 fb and 242 ± 77 fb for $\sqrt{s}$ =13 and 13.6 TeV, respectively. The signal is established with a statistical significance of 5.8 standard deviations, with 3.5 expected, compared to the background-only hypothesis.

Hayrapetyan, Aram [Yerevan Physics Institute]↗
DOE OSTI2022

Constraints from isoscaling on the source size in energetic heavy ion collisions

In the framework of the statistical multifragmentation model, the nuclear isoscaling analysis is extended to constrain the ratio between the sizes of the decaying sources formed in a collision between two heavy ions. It is found that the ratio between the probabilities of observing n fragments in each event, for each of the sources, follows a scaling law, similar to the traditional nuclear isoscaling. However, the corresponding slope is sensitive to the source sizes. This property is explained analytically using the grand-canonical ensemble. Furthermore, the extent to which our findings are affected by finite size effects and by the deexcitation of the hot primary fragments is investigated. The scaling turns out to be robust and weakly affected by effects implied by these two aspects. We also find that the Poisson distribution is a fairly good approximation to the above mentioned probabilities, associated with both the primordial fragments, produced at the breakup stage, and the final ones, found at the end of the fragment deexcitation process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Learning-based approach to plasticity in athermal sheared amorphous packings: Improving softness

The plasticity of amorphous solids undergoing shear is characterized by quasi-localized rearrangements of particles. While many models of plasticity exist, the precise relationship between the plastic dynamics and the structure of a particle’s local environment remains an open question. Previously, machine learning was used to identify a structural predictor of rearrangements called “softness.” Although softness has been shown to predict which particles will rearrange with high accuracy, the method can be difficult to implement in experiments where data are limited and the combinations of descriptors it identifies are often difficult to interpret physically. Here, we address both of these weaknesses, presenting two major improvements to the standard softness method. First, we present a natural representation of each particle’s observed mobility, allowing for the use of statistical models that are both simpler and provide greater accuracy in limited datasets. Second, we employ persistent homology as a systematic means of identifying simple, topologically informed, structural quantities that are easy to interpret and measure experimentally. We test our methods on two-dimensional athermal packings of soft spheres under quasi-static shear. We find that the same structural information that predicts small variations in the response is also predictive of where plastic events will localize. We also find that an excellent accuracy is achieved in athermal sheared packings using simply a particle’s species and the number of nearest neighbor contacts.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Accelerating high-strain continuum-scale brittle fracture simulations with machine learning

Failure in brittle materials under dynamic loading conditions is a result of the propagation and coalescence of microcracks. Simulating this discrete crack evolution at the continuum level is computationally expensive or, in some cases, intractable, resulting in the need to make broad assumptions or neglect key physics. In this work, we have developed an approach using machine learning that overcomes the current inability to represent meso-scale physics at the macro-scale. Our approach leverages damage and stress data from a computationally expensive high-fidelity model that explicitly resolves microcrack behavior to build an inexpensive machine learning emulator. Once trained, the machine learning emulator is used to predict the evolution of crack length statistics, which then informs a continuum-scale constitutive model. This results in a significant speed-up of the workflow by four orders of magnitude. Both the machine learning emulator and the continuum-scale model are validated against the high-fidelity model and experimental data, respectively, showing excellent agreement. There are two key findings. The first is that we can reduce the dimensionality of the problem, establishing that the machine learning emulator only needs the length of the longest crack and one of the maximum stress components to capture the necessary physics. Another compelling finding is that the emulator can be trained in one experimental setting and transferred successfully to predict behavior in a different setting.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An Analysis of FPGA LUT Bias and Entropy for Physical Unclonable Functions

Process variations within Field Programmable Gate Arrays (FPGAs) provide a rich source of entropy and are therefore well-suited for the implementation of Physical Unclonable Functions (PUFs). However, careful considerations must be given to the design of the PUF architecture as a means of avoiding undesirable localized bias effects that adversely impact randomness, an important statistical quality characteristic of a PUF. Here in this paper, we investigate a ring-oscillator (RO) PUF that leverages localized entropy from individual look-up table (LUT) primitives. A novel RO construction is presented that enables the individual paths through the LUT primitive to be measured and isolated at high precision, and an analysis is presented that demonstrates significant levels of localized design bias. The analysis demonstrates that delay-based PUFs that utilize LUTs as a source of entropy should avoid using FPGA primitives that are localized to specific regions of the FPGA, and instead, a more robust PUF architecture can be constructed by distributing path delay components over a wider region of the FPGA fabric. Compact RO PUF architectures that utilize multiple configurations within a small group of LUTs are particularly susceptible to these types of design-level bias effects. The analysis is carried out on data collected from a set of identically designed, hard macro instantiations of the RO implemented on 30 copies of a Zynq 7010 SoC.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Recurrent neural networks for short-term and long-term prediction of geothermal reservoirs

Accurate prediction of geothermal reservoir responses to alternative energy production scenarios is critical for optimizing the development of the underlying resources. While the conventional physics-based models offer a comprehensive prediction tool, data-driven models provide an efficient alternative to build fit-for-purpose predictive models by extracting and using the statistical patterns in the collected data to make predictions. The recurrent neural network (RNN) is a data-driven model that is commonly applied to predict time series sequences. This paper presents a variant of RNN that also utilizes the efficiency of convolutional neural networks (CNN) for the prediction of energy production from geothermal reservoirs. Specifically, a CNN–RNN architecture is developed that takes historical well controls as input (features) and their corresponding production response data as output (labels) to learn an input-output mapping that can predict the future well production responses/performance for any given future well control inputs. The model is paired with a labeling scheme to handle real field disturbances that create data gaps. In addition to the model structure, we introduce a thorough workflow for applying the model, which includes data pre-processing, feature selection, as well as different training strategies for short-term and long-term prediction. Finally, the performance and accuracy of the model are evaluated by applying it to multiple datasets, including a field reservoir model.

15 GEOTHERMAL ENERGY↗