Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical mechanics models”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Full-polarization millimeter wavelength variability of Sagittarius A * during the 2018 EHT campaign

Context. Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, provides a unique laboratory to study accretion dynamics and plasma processes near the event horizon. Aims. We investigated the variability and polarization properties of Sgr A* using ALMA observations during the 2018 Event Horizon Telescope campaign. Methods. We analyzed high-cadence full-polarization light curves from ALMA at millimeter wavelengths, performed time-series analysis, and investigated the temporal behavior during an X-ray flare observed by Chandra on 2018 April 24. The variability characteristics are compared with expectations from standard accretion flow models. Results. We find low variability in total intensity (σ/μ < 10%), but significantly higher variability in linear and circular polarization (∼30% and ∼50%, respectively). A time-series analysis reveals red-noise variability, with power spectral densities between −2 and −3 across all Stokes parameters. Polarized intensity shows stable intra-day timescales, while total intensity exhibits more variable timescales, suggesting distinct emission regions, with polarization likely arising from a coherent structure. On April 24, a statistically significant inter-band delay in polarized intensity coincides with a near-simultaneous X-ray and millimeter peak that deviates from the typical delayed flare scenario. This event also features enhanced millimeter variability and coherent polarization loop evolution. The observed simultaneity challenges standard models of transient synchrotron emission with cooling delays, favoring instead a scenario of continuous energy injection in an optically thin region. Conclusions. Our results offer new constraints on the physical mechanisms driving variability in Sgr A*, and provide key observational input for refining theoretical models of accretion and plasma behavior in the vicinity of supermassive black holes.

Galaxy: center↗
DOE OSTI2021

Odd Diffusivity of Chiral Random Motion

Diffusive transport is characterized by a diffusivity tensor which may, in general, contain both a symmetric and an antisymmetric component. Although the latter is often neglected, we derive Green-Kubo relations showing it to be a general characteristic of random motion breaking time-reversal and parity symmetries, as encountered in chiral active matter. In analogy with the odd viscosity appearing in chiral active fluids, we term this component the odd diffusivity. Furthermore, we show how odd diffusivity emerges in a chiral random walk model, and demonstrate the applicability of the Green-Kubo relations through molecular dynamics simulations of a passive tracer particle diffusing in a chiral active bath.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Bayesian Estimation of Oscillator Parameters: Toward Anomaly Detection and Cyber-Physical System Security

Cyber-physical system security presents unique challenges to conventional measurement science and technology. Anomaly detection in software-assisted physical systems, such as those employed in additive manufacturing or in DNA synthesis, is often hampered by the limited available parameter space of the underlying mechanism that is transducing the anomaly. As a result, the formulation of anomaly detection for such systems often leads to inverse or ill-posed problems, requiring statistical treatments. Here, we present Bayesian inference of unknown parameters associated with a generic actuator considered as a representative vital element of a cyber-physical system. Via a series of experimental input-output measurements, a transfer function for the actuator is obtained numerically, which serves as our model for the proposed method. Linear, nonlinear, and delayed dynamics may be assumed for the actuator response. By devising a code-based malicious signal, we study the efficacy of Bayesian inference for its potential to produce a detection, including uncertainty quantification, with a remarkably small number of input data points. Our approach should be adaptable to a variety of real-time cyber-physical anomaly detection scenarios.

47 OTHER INSTRUMENTATION↗
NASA NTRS2018

Cumulative Distribution of Ballistic Impact Failures of Common Twisted-Pair Data Cables at Orbital Speeds

Data wire cable runs are a significant presence on the exterior of the International Space Station (ISS), and continued ISS mission support requires detailed assessment of cables due to micrometeoroid and orbit debris (MMOD) impact. These data wire cables are twisted-pair cables consisting of two 22AWG stranded conductors and fillers inside a tight fitting braided copper shield. The copper shield and its contents are covered with a jacket that has a nominal outer diameter of 3.76 mm and beta-cloth tape. The ISS engineering community has identified two loss-of-function mechanisms for these cables: open circuits due to severed conductors within the cable, and short circuits due to contact between conductors or grounded components. As these data cables are low power systems, short circuits are not expected to burn away the contact, so both open and short circuits are considered permanent loss-of-function for the cable. A total of ninety-seven impact experiments have been performed into these cables to develop a statistical model for the failure of these cables to be used in reliability studies. The experimental work has yielded cumulative distribution functions for these cables for steel and aluminum components of the orbital debris environment at representative speeds and impact obliquities.

Miller, Joshua E.↗
DOE OSTI2022

Robust estimation of charge carrier diffusivity using transient photoluminescence microscopy

Transient microscopy has emerged as a powerful tool for imaging the diffusion of excitons and free charge carriers in optoelectronic materials. In many excitonic materials, extraction of diffusion coefficients can be simplified because of the linear relationship between signal intensity and local excited state population. However, in materials where transport is dominated by free charge carriers, extracting diffusivities accurately from multidimensional data is complicated by the nonlinear dependence of the measured signal on the local charge carrier density. To obtain accurate estimates of charge carrier diffusivity from transient microscopy data, statistically robust fitting algorithms coupled to efficient 3D numerical solvers that faithfully relate local carrier dynamics to raw experimental measurables are sometimes needed. Here, we provide a detailed numerical framework for modeling the spatiotemporal dynamics of free charge carriers in bulk semiconductors with significant solving speed reduction and for simulating the corresponding transient photoluminescence microscopy data. To demonstrate the utility of this approach, we apply a fitting algorithm using a Markov chain Monte Carlo sampler to experimental data on bulk CdS and methylammonium lead bromide (MAPbBr 3 ) crystals. Parameter analyses reveal that transient photoluminescence microscopy can be used to obtain robust estimates of charge carrier diffusivities in optoelectronic materials of interest, but that other experimental approaches should be used for obtaining carrier recombination constants. Additionally, simplifications can be made to the fitting model depending on the experimental conditions and material systems studied. Our open-source simulation code and fitting algorithm are made freely available to the scientific community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Direct cross-section measurement of the weak r -process Sr 88 ( α , n ) Zr 91 reaction in ν -driven winds of core-collapse supernovae

About half of the heavy elements beyond iron are known to be produced by the rapid neutron capture process, known as the r process. However, the astrophysical site producing the r process is still uncertain. Chemical abundances observed in several cosmic sites indicate that different mechanisms should be at play. For instance, the abundances around silver measured in a subset of metal-poor stars indicate the presence of a weak r process. This process may be active in neutrino-driven winds of core collapse supernovae where (α,n) reactions dominate the synthesis of Ζ ≈ 40 elements in the expelled materials. Scarcely measured, the rates of (α,n) reactions are determined from statistical Hauser-Feshbach calculations with α-optical-model potentials, which are still poorly constrained. Further, the uncertainties of the (α,n) reaction rates therefore make a significant contribution to the uncertainties of the abundances determined from stellar modeling. In this work, the 88 Sr (α,n)⁢ 91 Zr reaction which impacts the weak r-process abundances has been probed at astrophysics energy for the first time; directly measuring the total cross sections at astrophysical energies of 8.37–13.09 MeV in the center of mass (3.8–7.5 GK). Two measurements were performed at ATLAS with the electrically segmented ionization chamber MUSIC, in inverse kinematics, while following the active target technique. The cross sections of this α-induced reaction on 88 Sr, located at the shell closure N = 50, have been found to be lower than expected, by a factor of 3, despite recent statistical calculations validated by measurements on neighboring nuclei. This result encourages more experimental investigations of (α,n) reactions, at Ν = 50 and towards the neutron-rich side, to further test the predictive power and reliability of such calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Nanotomography for Quantitative 3D Particle Reconstruction

Particulates are ubiquitous across fuel cycle operations and carry critical information about particle formation, processing, and potential proliferation-related activities. Traditional analytical techniques, including micro-Raman spectroscopy and standard electron microscopy, are often limited in spatial resolution or dimensionality, particularly when used to examine metallic or submicron-scale features. Understanding particle morphology, phase distribution, and internal porosity is essential for constraining formation conditions, thermodynamic environments, and material transport behavior. In this report, we demonstrate the application of plasma focused ion beam nanotomography to reconstruct micron-scale particulates at nanoscale resolution. Using high-resolution backscattered electron imaging and Avizo software, we obtained 3D reconstructions that enabled quantitative analysis of particle morphology, phase composition, and internal voids. Representative examples include a Ta particle with a large central void and a composite particle with embedded tetrahedral crystalline structures. These reconstructions reveal structural and compositional details that are inaccessible through conventional 2D imaging. The results demonstrate that nanotomography provides both qualitative and quantitative insights into particle formation and behavior. Using nanotomography, porosity and phase distributions can be quantified to inform models of particle density, transport, and solidification conditions. Beyond technical insights, the workflow developed here establishes a transferable capability for analyzing heterogeneous particles and has potential applications in bulk materials studies via x-ray computed tomography or other volumetric imaging modalities. Ongoing efforts are focused on optimizing the workflow to process multiple particles simultaneously, increasing throughput and statistical robustness. Overall, this work illustrates the power of nanotomography as a tool for connecting particulate morphology to formation mechanisms, composition, and transport, thereby strengthening analytical capabilities for nuclear forensics, fuel cycle analysis, and related scientific investigations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1991

Hexagonal and nematic phases of chains. I - Correlation functions. II - Phase transitions

The statistical mechanics of a system of semiflexible chains, which can represent polymer liquid crystals, long-chain biomolecules, stiff wormlike micelles, or columns of discotic liquid crystals, are examined. A continuum theory is used to calculate static correlation functions in the hexagonal and nematic phases. Two correlation functions are considered: (1) the structure factor which describes fluctuations in the density; and (2) the director fluctuation spectrum, which describes fluctuations in the local optical axis. In addition, a model is developed for the phase transitions of a system of infinitely long, semiflexible chains which interact through a steric, excluded-volume repulsion. The model yields generic phase diagrams in terms of pressure or density vs. persistence length or temperature.

Selinger, Jonathan V.↗
DOE OSTI2023

Collisional-energy-cascade model for nonthermal velocity distributions of neutral atoms in plasmas

Nonthermal velocity distributions with much greater tails than a Maxwellian have been observed for radical atoms in plasmas for a long time. Historically, such velocity distributions have been modeled by a two-temperature Maxwell distribution. In this paper, I propose a model based on collisional energy cascade, which has been studied in the field of granular materials. In the collisional energy cascade, a particle ensemble undergoes energy input at the high-energy region, entropy production by elastic collisions among particles, and energy dissipation. For radical atoms, energy input may be caused by the Franck-Condon energy of molecular dissociation or charge-exchange collision with hot ions, and the input energy is eventually dissipated by collisions with the walls. Here, I show that the steady-state velocity distribution in the collisional energy cascade is approximated by the generalized Mittag-Leffler distribution, which is a one-parameter extension of the Maxwell distribution. This parameter indicates the degree of the nonthermality and is related to the relative importance of energy dissipation over entropy production. This model is compared with a direct molecular dynamics simulation for simplified gaseous systems with energy input, as well as some experimentally observed velocity distributions of light radicals in plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Signatures of a liquid–liquid transition in an ab initio deep neural network model for water

Significance Water is central across much of the physical and biological sciences and exhibits physical properties that are qualitatively distinct from those of most other liquids. Understanding the microscopic basis of water’s peculiar properties remains an active area of research. One intriguing hypothesis is that liquid water can separate into metastable high- and low-density liquid phases at low temperatures and high pressures, and the existence of this liquid–liquid transition could explain many of water’s anomalous properties. We used state-of-the-art approaches in computational quantum chemistry, statistical mechanics, and machine learning and obtained evidence consistent with a liquid–liquid transition, supporting the argument for the existence of this phenomenon in real water.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Spin polarization induced by magnetic field and the relativistic Barnett effect

First, I study the analogy between the magnetization of a material and the spin polarization of particles in a fluid. Using the relativistic version of the Barnett effect, i.e. the magnetization of a material induced by mechanical rotation, the spin polarization induced by thermal vorticity is obtained within a purely classical model, where spin is treated as an intrinsic magnetic moment and rotation is included as a non-inertial effect. I argue that since spin polarization induced by thermal vorticity can be obtained in a classical theory, it can not be dominated by quantum anomalies. Second, the spin polarization induced by magnetic field is obtained for a fluid at local thermal equilibrium using statistical quantum field theory. The obtained formula is valid beyond the weak field approximation and when contributions from the non-homogeneity of the magnetic field are small. The exact form of spin polarization is studied for a free Dirac field at global equilibrium, and, like magnetic susceptibility, it oscillates according to the de Haas - van Alphen effect. Finally, I briefly review how magnetic field contributes to the difference between the spin polarization of Λ and $\barΛ$ observed in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Debris-induced consequences on turbulence and vorticity in solar photovoltaic module-generated array wakes

Particle-laden flows in solar photovoltaic (PV) systems are inevitable, where wind-swept debris in open environments are carried by high winds and turbulence, coating panel surfaces or damaging structures. Particle deposition, or soiling, is a well-known issue for large-scale plants which rely on uninhibited solar rays for optimal production. But understanding the mechanisms leading to soiling requires a physical and fluid dynamics-centered focus, since turbulence dominates PV panel wakes and is also known to alter particle concentration and trajectories. This study presents an experimental campaign toward consequences of particle-laden flow between two model PV panels using time-resolved particle image velocimetry. The model array was subjected to varied particle volume fractions, including a tracer particle case and a water droplet case. Characterization of mean velocity, turbulence statistics, and mean kinetic energy within the single phase and, separately, particle phase flows showed modified features due to particle inertia. Images captured at a frequency of 1 kHz in the near wake of the upstream panel allow for a first experimental look at vorticity and convective velocity of vortex structures for single-phase and particle-phase flows which are crucial to debris transport and soiling in PV environments.

Energy & Fuels↗
NASA NTRS1993

A regional analysis of cloudy mean spherical albedo over the marine stratocumulus region and the tropical Atlantic Ocean

Since clouds are the largest variable in Earth's radiation budget, it is critical to determine both the spatial and temporal characteristics of their radiative properties. The relationships between cloud properties and cloud fraction are studied in order to supplement grid scale parameterizations. The satellite data used is from three hourly ISCCP (International Satellite Cloud Climatology Project) and monthly ERBE (Earth Radiation Budget Experiment) data on a 2.5 deg x 2.5 deg latitude-longitude grid. Mean cloud spherical albedo, the mean optical depth distribution, and cloud fraction are examined and compared off the coast of California and the mid-tropical Atlantic for July 1987 and 1988. Individual grid boxes and spatial averages over several grid boxes are correlated to Coakley's theory of reflection for uniform and broken layered cloud and to Kedem, et al.'s findings that rainfall volume and fractional area of rain in convective systems is linear. Kedem's hypothesis can be expressed in terms of cloud properties. That is, the total volume of liquid in a box is a linear function of cloud fraction. Results for the marine stratocumulus regime indicate that albedo is often invariant for cloud fractions of 20% to 80%. Coakley's satellite model of small and large clouds with cores (1 km) and edges (100 m) is consistent with this observation. The cores maintain high liquid water concentrations and large droplets while the edges contain low liquid water concentrations and small droplets. Large clouds are just a collection of cores. The mean optical depth (TAU) distributions support the above observation with TAU values of 3.55 to 9.38 favored across all cloud fractions. From these results, a method based upon Kedem, et al's theory is proposed to separate the cloud fraction and liquid water path (LWP) calculations in a general circulation model (GCM). In terms of spatial averaging, a linear relationship between albedo and cloud fraction is observed. For tropical locations outside the Intertropical Convergence Zone (ITCZ), results of cloud fraction and albedo spatial averaging followed that of the stratus boxes containing few overcast scenes. Both the ideas of Coakley and Kedem, et al. apply. Within the ITCZ, the grid boxes tended to have the same statistical properties as stratus boxes containing many overcast scenes. Because different dynamical forcing mechanisms are present, it is difficult to devise a method for determining subgrid scale variations. Neither of the theories proposed by Kedem, et al. or Coakley works well for the boxes with numerous overcast scenes.

Ginger, Kathryn M.↗
NASA NTRS1993

A Regional Analysis of Cloudy Mean Spherical Albedo over the Marine Stratocumulus Region and the Tropical Atlantic Ocean

Since clouds are the largest variable in Earth's radiation budget, it is critical to determine both the spatial and temporal characteristics of their radiative properties. This study examines the relationships between cloud properties and cloud fraction in order to supplement grid scale parameterizations. The satellite data used in this study is from three hourly ISCCP (International Satellite Cloud Climatology Project) and monthly ERBE (Earth Radiation Budget Experiment) data on a 2.50 x 2.50 latitude-longitude grid. Mean cloud spherical albedo, the mean optical depth distribution and cloud fraction are examined and compared off the coast of California and the mid-tropical Atlantic for July 1987 and 1988. Individual grid boxes and spatial averages over several grid boxes are correlated to Coakleys (1991) theory of reflection for uniform and broken layered cloud and to Kedem, et al.(1990) findings that rainfall volume and fractional area of rain in convective systems is linear. Kedem's hypothesis can be expressed in terms of cloud properties. That is, the total volume of liquid in a box is a linear function of cloud fraction. Results for the marine stratocumulus regime indicate that albedo is often invariant for cloud fractions of 20% to 80%. Coakley's satellite model of small and large clouds with cores (1 km) and edges (100 in) is consistent with this observation. The cores maintain high liquid water concentrations and large droplets while the edges contain low liquid water concentrations and small droplets. Large clouds are just a collection of cores. The mean optical depth (TAU) distributions support the above observation with TAU values of 3.55 to 9.38 favored across all cloud fractions. From these results, a method based upon Kedem, et al. theory is proposed to separate the cloud fraction and liquid water path (LWP) calculations in a general circulation model (GCM). In terms of spatial averaging, a linear relationship between albedo and cloud fraction is observed. For tropical locations outside the Intertropical Convergence Zone (ITCZ), results of cloud fraction and albedo spatial averaging followed that of the stratus boxes containing few overcast scenes. Both the ideas of Coakley and Kedem, et al. apply. Within the ITCZ, the grid boxes tended to have the same statistical properties as stratus boxes containing many overcast scenes. Because different dynamical forcing mechanisms are present, it is difficult to devise a method for determining subgrid scale variations. Neither of the theories proposed by Kedem, et al. or Coakley works well for the boxes with numerous overcast scenes.

Ginger, Kathryn M.↗
NASA NTRS2010

On the Phase-Averaged Spectrum of Pulsars and Shape of Their Cutoffs

All gamma ray pulsars exhibit an exponential cutoff in their spectra and for bright pulsars the statistics are sufficiently high to study the detailed shape of the cutoff. The phase averaged spectra of some pulsars exhibit a sub-exponential cutoff, not predicted by any single physical mechanism. Further studies clarified that (his gentler average cutoff is a consequence of having significant variations of the cutoff energy in the phase-resolved spectrum. In conclusion, the phase-averaged spectrum of a pulsar is not a physical quantity to test high-energy emission models.

Celik, O.↗
DOE OSTI2020

Jitter-induced Max-of-N fluence distribution at the National Ignition Facility

Max-of-N fluence is the maximum peak fluence at a given location over N number of shots and is important for calculating fluence dependence for intrinsic laser-induced optic damage. Previously, we observed the Max-of-N effect on the National Ignition Facility and developed an ad-hoc model to calculate its effect. In this work, we attempt to understand the fundamental mechanism that causes this Max-of-N effect. We conclude that the primary fundamental mechanism responsible for this effect is dominated by the combination of fluence variations and pointing jitter of the laser. This discovery both strengthens our model for predicting optics longevity and gives us insight into how to mitigate this effect.

42 ENGINEERING↗
DOE OSTI2021

Quantifying the Influence of Charge Rate and Cathode-Particle Architectures on Degradation of Li-Ion Cells Through 3D Continuum-Level Damage Models

In this article, we develop a 3D, continuum-level damage model implemented on statistically generated LiNi 0.5 Mn 0.3 Co 0.2 O (NMC 532) secondary cathode particles. The primary motivation of the particle-level model is to inform cathode-particle design through detailed exploration of the influence of secondary and primary particle sizes on the damage predicted during operation, and determine charging profiles that reduce cathode fracture. The model considers NMC 532 secondary particles containing an agglomeration of anisotropic, randomly oriented grains. These brittle, Ni-based cathodes are prone to mechanical degradation, which reduces overall battery cycle life. The model predicts that secondary-particle fracture is primarily due to non-ideal grain interactions and high-rate charge demands. The model predicts that small secondary-particles with large grains develop significantly less damage than larger secondary particles with small grains. The model predicts most of the chemo-mechanical damage accumulates in the first few cycles. The chemo-mechanical model predicts monotonically increasing capacity fade with cycling and rate. Comparing to experimental results, the model is well suited for capturing initial capacity fade mechanisms, but additional physics is required to capture long-term capacity fade effects.

25 ENERGY STORAGE↗
DOE OSTI2021

Regularized machine learning on molecular graph model explains systematic error in DFT enthalpies

Abstract A major goal of materials research is the discovery of novel and efficient heterogeneous catalysts for various chemical processes. In such studies, the candidate catalyst material is modeled using tens to thousands of chemical species and elementary reactions. Density Functional Theory (DFT) is widely used to calculate the thermochemistry of these species which might be surface species or gas-phase molecules. The use of an approximate exchange correlation functional in the DFT framework introduces an important source of error in such models. This is especially true in the calculation of gas phase molecules whose thermochemistry is calculated using the same planewave basis set as the rest of the surface mechanism. Unfortunately, the nature and magnitude of these errors is unknown for most practical molecules. Here, we investigate the error in the enthalpy of formation for 1676 gaseous species using two different DFT levels of theory and the ‘ground truth values’ obtained from the NIST database. We featurize molecules using graph theory. We use a regularized algorithm to discover a sparse model of the error and identify important molecular fragments that drive this error. The model is robust to rigorous statistical tests and is used to correct DFT thermochemistry, achieving more than an order of magnitude improvement.

36 MATERIALS SCIENCE↗