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,567 records · Page 87

An engineering model of woven composites based on micromechanics

Composites with three-dimensional woven architectures exhibit large strains to failure when compared to composites made up of the same materials but not with three-dimensional interlocking tows. The fracture mechanics of such three-dimensional architectures is a subject requiring substantial investigation and experimental testing. Classical fracture mechanics concepts (for instance, an isolated defect in a homogeneous body) will not be applicable to the woven fracture test specimen. The use of an isolated singularity to characterize an entire specimen is inadequate when the density of defects is considerable and the material is heterogeneous. Modelling of such a complex system requires a great deal of insight and consideration as well as prudent choices of model sizes to make numerical schemes feasible. The purpose of this manuscript is to review our recently acquired knowledge of damage accumulation in woven composites and to describe a practicable model of the macroscopic behavior in these and other complex composite architectures based on such knowledge. In this manuscript, discussion will be limited to uniaxial compressive loading; considerations of general loading (monotonic and cyclic) will appear in a subsequent manuscript. Our modelling efforts may be briefly described as follows: the composite is subdivided into microstructural elements (microelements) in which the micromechanical modelling is either understood rigorously or can be represented adequately by statistical parameters. There can be microstructural elements for many different types of composite components, such as the various types of warp and weft and matrix for three-dimensional woven composites. The physical dimensions of microelements are made as large as possible while the response within the element can still be represented by a single micromechatlical calculation. The various elements are linked together(sometimes by associating distinct corners and edges, sometimes by superposition) in a pattern which resembles a particular weave architecture. The model can then be loaded in any manner and the linear and nonlinear elastic responses of representative weaves can be calculated. After the elastic regime, the fracture response is determined by monitoring the damage accumulation.

Carter, W. C.↗

Quantum Cohesion Oscillation of Electron Ground State in Low Temperature Laser Plasma

The development of radically new technological and economically efficient methods for obtaining chemical products and for producing new materials with specific properties requires the study of physical and chemical processes proceeding at temperature of 10(exp 3) to 10(exp 4) K, temperature range of low temperature plasma. In our paper, by means of Wigner matrix of quantum statistical theory, a formula is derived for the energy of quantum coherent oscillation of electron ground state in laser plasma at low temperature. The collective behavior would be important in ion and ion-molecule reactions.

Zhao, Qingxun↗

Utility of AIRS Retrievals for Climate Studies

Satellites provide an ideal platform to study the Earth-atmosphere system on practically all spatial and temporal scales. Thus, one may expect that their rapidly growing datasets could provide crucial insights not only for short-term weather processes/predictions but into ongoing and future climate change processes as well. Though Earth-observing satellites have been around for decades, extracting climatically reliable information from their widely varying datasets faces rather formidable challenges. AIRS/AMSU is a state of the art infrared/microwave sounding system that was launched on the EOS Aqua platform on May 4, 2002, and has been providing operational quality measurements since September 2002. In addition to temperature and atmospheric constituent profiles, outgoing longwave radiation and basic cloud parameters are also derived from the AIRS/AMSU observations. However, so far the AIRS products have not been rigorously evaluated and/or validated on a large scale. Here we present preliminary assessments of monthly and 8-day mean AIRS "Version 4.0" retrieved products (available to the public through the DAAC at NASA/GSFC) to assess their utility for climate studies. First we present "consistency checks" by evaluating the time series of means, and "anomalies" (relative to the first 4 full years' worth of AIRS "climate statistics") of several climatically important retrieved parameters. Finally, we also present preliminary results regarding interrelationships of some of these geophysical variables, to assess to what extent they are consistent with the known physics of climate variability/change. In particular, we find at least one observed relationship which contradicts current general circulation climate (GCM) model results: the global water vapor climate feedback which is expected to be strongly positive is deduced to be slightly negative (shades of the "Lindzen effect"?). Though the current AIRS climatology covers only -4.5 years, it will hopefully extend much further into the future.

Molnar, Guyla I.↗

A Hybrid Transfer Function Procedure for Broadband Auralization Within Small Flight Vehicle Interiors

Auralization of the sound fields within flight vehicle interiors is of interest to aircraft designers considering crew and passenger ride quality. Computer models can be used to predict the filtering effects that the fuselage structure and interior volume have on the exterior sound field in order to simulate the resulting cabin interior noise. However, this can become expensive when a wide frequency range of analysis suitable for passenger auralization is desired. This is particularly the case when using deterministic modeling methods such as finite element analysis (FEA) wherein prohibitively high levels of detail and finite element discretizations are required to resolve the wavelengths at high frequencies. On the contrary, Statistical Energy Analysis (SEA), while not capturing important modal behaviors at low frequencies, is better suited for higher frequency analyses where structural and interior dynamics exhibit high modal overlap and spatially averaged physical quantities suffice. For this reason, a hybrid procedure is considered here that utilizes FE modal solutions at lower frequencies and SEA results at higher frequencies. A transfer function dataset containing the exterior to interior vibroacoustic filtering effect is calculated accordingly and subsequently stored to be used for a variety of exterior dynamic load cases. At low frequencies, the transfer function dataset detail is related to the level of FE model refinement whereas at high frequencies it functions as an equalizer array with a level of detail corresponding to the SEA subsystem partitioning scheme. The above procedure is described herein and demonstrated on a six passenger flight vehicle with eight propulsors in a hover flight condition.

auralization↗

A Hybrid Transfer Function Procedure for Broadband Auralization Within Small Flight Vehicle Interiors

Auralization of the sound fields within flight vehicle interiors is of interest to aircraft designers considering crew and passenger ride quality. Computer models can be used to predict the filtering effects that the fuselage structure and interior volume have on the exterior sound field in order to simulate the resulting cabin interior noise. However, this can become expensive when a wide frequency range of analysis suitable for passenger auralization is desired. This is particularly the case when using deterministic modeling methods such as finite element analysis (FEA) wherein prohibitively high levels of detail and finite element discretizations are required to resolve the wavelengths at high frequencies. On the contrary, Statistical Energy Analysis (SEA), while not capturing important modal behaviors at low frequencies, is better suited for higher frequency analyses where structural and interior dynamics exhibit high modal overlap and spatially averaged physical quantities suffice. For this reason, a hybrid procedure is considered here that utilizes FE modal solutions at lower frequencies and SEA results at higher frequencies. A transfer function dataset containing the exterior to interior vibroacoustic filtering effect is calculated accordingly and subsequently stored to be used for a variety of exterior dynamic load cases. At low frequencies, the transfer function dataset detail is related to the level of FE model refinement whereas at high frequencies it functions as an equalizer array with a level of detail corresponding to the SEA subsystem partitioning scheme. The above procedure is described herein and demonstrated on a six passenger flight vehicle with eight propulsors in a hover flight condition.

auralization↗

Geometry, Disorder and Phase Transitions in Topological States of Matter

The quantum Hall effect is the birthplace of topological states of matter, a major theme at the forefront of condensed matter physics in the past two decades. The fractional quantum Hall (FQH) effect revolutionized our understanding of phases of electronic matter. FQH states support exotic fractionally charged excitations that obey Abelian or non-Abelian fractional statistics, which are topological excitations that result from the underlying topological order. During this project, our group discovered a previously unrecognized geometric degree of freedom of incompressible FQH states and studied that for a variety of gapped FQH states. We brought this new concept into direct contact with experiments for the first time by generalizing it to Fermi-liquid states of composite fermions. Using the newly formulated powerful infinite Density Matrix Renormalization Group method, our numerical calculations yielded a parameter free prediction that was found to be in excellent agreement with experimental findings on electron systems in semiconductor heterostructures. In parallel, we performed extensive numerical studies on different, competing phases at various Landau level filling factors, and quantum phase transitions that result from such a competition, e.g. Abelian-non-Abelian phase transitions in bilayer systems. We studied geometrical excitations dubbed “gravitons” (because of their analogy with excitations in the theory of gravitation) and ways to excite and detect them, and explored how they couple with topological excitations. In graphene-based chiral materials, we realized the ability to tune through different incompressible and compressible states in a single Landau level, and found appropriate experimental parameters for the exploration of universal Luttinger liquid behavior not obtained in semiconductor-based electron systems. We showed that topological systems had a very different response from nontopological systems to strong disorder (many-body localization) as well as periodic drives.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measurement of the anomalous spin precession frequency $\omega_a$ in the Muon $g-2$ experiment at Fermilab

The muon anomaly, $a_\mu=(g_{\mu}-2)/2$, is a low-energy observable which can be both measured and computed to high precision, making it a sensitive test of the Standard Model (SM) and a probe for new physics. The current discrepancy between the experimental value and the Standard Model calculation from the Muon $g-2$ Theory Initiative is $a_{\mu}^{exp}-a_{\mu}^{SM}=(251\pm59)\cdot10^{-11}$, with a significance of $4.2\,\sigma$.The Fermilab E989 experiment aims, with the full statistical power, to improve by a factor of four the precision of the measurement.In April 2021 the collaboration published the first measurement, based on the first year of data taking.This paper will present the status of the measurement of the anomalous muon spin precession frequency, $\omega_a$, performed on the datasets collected during Run-2 and Run-3 (2019 and 2020 campaigns), with a preliminary projection of the systematic uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

On principles of emergent organization

After more than a century of concerted effort, physics still lacks basic principles of spontaneous organization. To appreciate why, we first state the problem, outline historical approaches, and survey the present state of the physics of self-organization. This frames the particular challenges arising from mathematical intractability and the resulting need for computational approaches, as well as those arising from a chronic failure to define structure. Then an overview of two modern mathematical formulations of organization—intrinsic computation and evolution operators—lays out a way to overcome these challenges. Together, the vantage point they afford shows how to account for the emergence of structured states via a statistical mechanics of systems arbitrarily far from equilibrium. The result is a constructive path forward to principles of organization that builds on mathematical identification of structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A statistical analysis of simulated data from potential stratospheric satellite missions

A simplified one-dimensional model is used to analyze simulated data from typical proposed satellite measurements of trace species in the stratosphere. Uncertainties are introduced in simulated species concentration measurements and applied to the basic physical and chemical model data. The results of the calculations are then used to predict uncertainties in the species concentrations computed by the model. Candidate satellite missions can be evaluated for net information acquired by applying the error analysis. An abbreviated depiction of the computational procedure is given, and sample histograms based on an accumulation of 200 calculations are shown for H2O2 and ClNO3, demonstrating that emergent species distributions are roughly Gaussian in the logarithm of number density. Calculated vertical concentration profiles for these two species are presented.

Nealy, J. E.↗

Coronal Heating as Determined by the Solar Flare Frequency Distribution Obtained by Aggregating Case Studies

Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counterintuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms that could explain it: nanoflares or Alfvén waves. To date, neither can be directly observed. Nanoflares are, by definition, extremely small, but their aggregate energy release could represent a substantial heating mechanism, presuming they are sufficiently abundant. One way to test this presumption is via the flare frequency distribution, which describes how often flares of various energies occur. If the slope of the power law fitting the flare frequency distribution is above a critical threshold, α = 2 as established in prior literature, then there should be a sufficient abundance of nanoflares to explain coronal heating. We performed >600 case studies of solar flares, made possible by an unprecedented number of data analysts via three semesters of an undergraduate physics laboratory course. This allowed us to include two crucial, but nontrivial, analysis methods: preflare baseline subtraction and computation of the flare energy, which requires determining flare start and stop times. We aggregated the results of these analyses into a statistical study to determine that α = 1.63 ± 0.03. This is below the critical threshold, suggesting that Alfvén waves are an important driver of coronal heating.

Astrostatistics distributions↗

The Jefferson Lab Eta Factory Experiment and Applications of PbWO4 Calorimeters in Future Experimental Facilities

The goal of the new JLab Eta Factory (JEF) experiment, conducted with the GlueX detector in Hall D at Jefferson Lab, is to perform measurements of various ¿(') decays with a primary focus on rare neutral modes. The experiment’s physics program ranges from precision tests of low-energy QCD to searches for gauge bosons with masses below 1 GeV that could couple the Standard Model (SM) sector to the dark sector. The experiment will collect a high-statistics data sample of ¿(') mesons produced via a beam of tagged photons. The GlueX detector features a large, nearly uniform acceptance for both neutral and charged particles, enabling efficient identification of complex multi-particle final states. To meet the requirements of the JEF experiment, the inner section of the forward lead-glass calorimeter in the GlueX detector has been upgraded with lead tungstate (PbWO4) scintillating crystals. PbWO4 offers exceptional characteristics, such as a small radiation length and Molire radius, and large light yield, that make it ideal for constructing high- granularity, high-resolution, radiation-hard detectors. These properties enable excellent spatial separation and energy resolution of reconstructed electromagnetic showers, establishing PbWO4 as the material of choice for many high-precision experiments. The JEF experiment began data collection in April 2025 and will operate concurrently with the GlueX experiment, whose primary objective is the search for gluonic excitations in the meson spectrum. I will give an overview of the JEF experiment, the GlueX detector, and the feasibility of further upgrades to support future ¿ physics studies. Special attention will be given to the newly constructed PbWO4 scintillating calorimeter and recent advancements in calorimeter instrumentation.

Somov, Alexander↗

Uncertainty Quantification for a Multi-phase Carbon Equation of State Model

Many physics models have tunable parameters that are calibrated by matching the model output to experimental or calculated data. However, given that calibration data often contain uncertainty and that different model parameter sets might result in a very similar simulated output for a finite calibration data set, it is advantageous to provide an ensemble of parameter sets that are consistent with the calibration data. Uncertainty quantification (UQ) provides a means to generate such an ensemble in a statistically rigorous fashion. In this work, we perform UQ for a multi-phase equation of state (EOS) model for carbon containing the diamond, graphite, and liquid phases. We use a Bayesian framework for the UQ and introduce a novel strategy for including phase diagram information in the calibration. The method is highly general and accurately reproduces the calibration data without any material-specific prior knowledge of the EOS model parameters.

36 MATERIALS SCIENCE↗

The PENELOPE Physics Models and Transport Mechanics. Implementation into Geant4

A translation of the PENELOPE physics subroutines to C++, designed as an extension of the GEANT4 toolkit, is presented. The Fortran code system PENELOPE performs Monte Carlo simulation of coupled electron-photon transport in arbitrary materials for a wide energy range, nominally from 50 eV up to 1 GeV. PENELOPE implements the most reliable interaction models that are currently available, limited only by the required generality of the code. In addition, the transport of electrons and positrons is simulated by means of an elaborate class II scheme in which hard interactions (involving deflection angles or energy transfers larger than pre-defined cutoffs) are simulated from the associated restricted differential cross sections. After a brief description of the interaction models adopted for photons and electrons/positrons, we describe the details of the class-II algorithm used for tracking electrons and positrons. The C++ classes are adapted to the specific code structure of GEANT4. They provide a complete description of the interactions and transport mechanics of electrons/positrons and photons in arbitrary materials, which can be activated from the G4ProcessManager to produce simulation results equivalent to those from the original PENELOPE programs. The combined code, named PENG4, benefits from the multi-threading capabilities and advanced geometry and statistical tools of GEANT4.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the Search for Low W 0

The magnitude of the vacuum expectation value of the Gukov-Vafa-Witten superpotential |W o | plays a central role in the phenomenology of type IIB flux compactifications. Recent analytical constructions have shown that perturbatively flat vacua can be used to obtain very low values of |W o |. We present systematic algorithms to carry out exhaustive numerical searches for such vacua. Here we also analyse them in the statistical context, as part of the entire ensemble of type IIB flux vacua at low |W o |. Our preliminary analysis indicates that these perturbatively flat vacua are statistically sparse in the whole set of vacua at low |W o | as calculated by Denef and Douglas. Two-moduli examples are used to illustrate these more general findings in specific settings. We find that these simple (two moduli) cases are good examples for existence proofs but they do not feature a large statistical tuning freedom for phenomenological applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Comparison between mirror Langmuir probe and gas-puff imaging measurements of intermittent fluctuations in the Alcator C-Mod scrape-off layer

Statistical properties of the scrape-off layer plasma fluctuations are studied in ohmically heated plasmas in the Alcator C-Mod tokamak. For the first time, plasma fluctuations as well as parameters that describe the fluctuations are compared across measurements from a mirror Langmuir probe (MLP) and from gas-puff imaging (GPI) that sample the same plasma discharge. This comparison is complemented by an analysis of line emission time-series data, synthesized from the MLP electron density and temperature measurements. The fluctuations observed by the MLP and GPI typically display relative fluctuation amplitudes of order unity together with positively skewed and flattened probability density functions. Such data time series are well described by an established stochastic framework that models the data as a superposition of uncorrelated, two-sided exponential pulses. The most important parameter of the process is the intermittency parameter, $γ=τ_d/τ_w$, where τd denotes the duration time of a single pulse and $τ_w$ gives the average waiting time between consecutive pulses. In this study we show, using a new deconvolution method, that these parameters can be consistently estimated from different statistics of the data. We also show that the statistical properties of the data sampled by the MLP and GPI diagnostic are very similar. Finally, a synthetic GPI signal using only plasma parameters sampled by the MLP shows qualitatively different fluctuation statistics from the measured GPI signal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Methane Adsorption on Heteroatom-Modified Maquettes of Porous Carbon Surfaces

Experimental and theoretical studies disagree on the energetics of methane adsorption on carbon materials. However, this information is critical for the rational design and optimization of the structure and composition of adsorbents for natural gas storage. The delicate nature of dispersion interactions, polarization of both the adsorbent and the adsorbate, interplay between H-bonding and tetrel bonding, and induced dipole/Coulomb interactions inherent to methane physisorption require computational treatment at the highest possible level of theory. In this study, we employed the smallest reasonable computational model, a maquette of porous carbon surfaces with a central site for substitution and methane binding. The most accurate predictions of methane adsorption energetics were achieved by electron-correlated molecular orbital theory CCSD(T) and hybrid density functional theory MN15 calculations employing a saturated, all-electron basis set. Here, the characteristic geometry of methane adsorption on a carbon surface (“lander approach”) arises due to bonding interactions of the adsorbent π-system with the proximal H–C bonds of methane, in addition to tetrel bonding between the antibonding orbital of the distal C–H bond and the central atom of the maquette (C, B, or N). The polarization of the electron density, structural deformations, and the comprehensive energetic analysis clearly indicate a ∼3 kJ mol –1 preference for methane binding on the N-substituted maquette. The B-substituted maquette showed a comparable or lower binding energy than the unsubstituted, pure C model, depending on the level of theory employed. The calculated thermodynamic results indicate a strategy for incorporating electron-enriched substitutions (e.g., N) into carbon materials as a way to increase methane storage capacity over electron-deficient (e.g., B) modifications. The thermochemical analysis was revised for establishing a conceptual agreement between the experimental isosteric heat of adsorption and the binding enthalpies from statistical thermodynamics principles.

03 NATURAL GAS↗

Accelerating Ab Initio Simulation via Nested Monte Carlo and Machine Learned Reference Potentials

As a corollary of the rapid advances in computing, ab initio simulation is playing an increasingly important role in modeling materials at the atomic scale. Two strategies are possible, ab initio Monte Carlo (AIMC) and molecular dynamics (AIMD) simulation. The former benefits from exact sampling from the correct thermodynamic distribution, while the latter is typically more efficient with its collective all-atom coordinate updates. In this study, using a relatively simple test model comprised of helium and argon, we show that AIMC can be brought up to, and even above, the performance levels of AIMD via a hybrid nested sampling/machine learning (ML) strategy. Here, ML provides an accurate classical reference potential (up to three-body explicit interactions) that can pilot long collective Monte Carlo moves that are accepted or rejected in toto à la nested Monte Carlo (NMC); this is in contrast to the single move nature of a naive implementation. Our proposed method only requires a small up front expense from evaluating the ab initio energies and forces of (100) random configurations for training. Importantly, our method does not totally rely on the trained, assuredly imperfect, interaction. We show that high performance and exact sampling at the desired level of theory can be realized even when the trained interaction has appreciable differences from the ab initio potential. Remarkably, at the highest levels of performance realized via our approach, a pair of statistically uncorrelated atomic configurations can be generated with (1) ab initio calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗