Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STATISTICAL MECHANICS”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Transient work function gating: A new photoemission regime

We present the theoretical basis for a new photoemission regime, transient work function gating (TWFG), which temporally and energetically gates photoemission and produces near-threshold photoelectrons with thermally limited emittance, percent-level quantum efficiency, and control over temporal coherence. The technique consists of actively gating the work function of a generalized photocathode using a non-ionizing long-wavelength optical field to produce an adiabatic modulation of the carrier density at their surface. We examine TWFG as a means to circumvent the long-standing trade-off between low emittance and high quantum efficiency, untethered to particle source or photocathode specifics. TWFG promises new opportunities in photoemission physics for next generation electron and accelerator-based x-ray photon sources.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Piecewise interaction picture density matrix quantum Monte Carlo

The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact N-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC (IP-DMQMC) method that overcomes the limitation of only sampling one inverse temperature point at a time, instead allowing for the sampling of a temperature range within a single calculation, thereby reducing the computational cost. At the target inverse temperature, instead of ending the simulation, we incorporate a change of picture away from the interaction picture. The resulting equations of motion have piecewise functions and use the interaction picture in the first phase of a simulation, followed by the application of the Bloch equation once the target inverse temperature is reached. We find that the performance of this method is similar to or better than the DMQMC and IP-DMQMC algorithms in a variety of molecular test systems.

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

Progress toward hydro-equivalent ignition in OMEGA direct-drive DT-layered implosions

Considerable progress has been made in deuterium-tritium-layered implosion experiments on the OMEGA Laser System, bringing the prospects for thermonuclear ignition in direct-drive configurations with megajoule-class lasers closer to reality. Doing so has required navigating the balance between improved 1D performance and multidimensional stability. Using statistical modeling based on over 350 cryogenic implosions to identify various degradation mechanisms, and combined with multidimensional simulations and experimental techniques such as target offsets to combat residual flows, core conditions have repeatably been achieved that extrapolate to the burning-plasma state when scaled to 2.15 MJ of symmetric laser illumination. Using high implosion velocities (⁠> 450 km/s) and moderately high adiabats (⁠~5⁠), these experiments produced record-high scaled Lawson parameters in direct drive equal to 89 ± 2% of that required for ignition with expected yields of up to 1.5 ± 0.2 MJ. To improve these results still further, focused physics studies are performed to improve physics understanding and identify routes to even greater performance. Recent studies include investigations into the impact of mounting features, laser imprint, reduced fuel temperatures, and greater on-target intensities through subscale experiments. This manuscript gives a summary of the cryogenic direct-drive program on the OMEGA laser, including routes taken to achieve the current best performance, the status of recent focused physics investigations, and future designs—such as target solutions to laser imprint and reducing vapor density to increase convergence—that are expected lead to the demonstration of hydro-equivalent ignition on OMEGA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Adaptive Client Selection in Federated Learning: A Network Anomaly Detection Use Case

Federated Learning (FL) has become a ubiquitous approach for training machine learning models on decentralized data, addressing the myriad privacy concerns inherent in traditional centralized methods. However, the efficiency of FL depends on effective client selection and robust privacy preservation mechanisms. Inadequate client selection may lead to suboptimal model performance, while insufficient privacy measures risk exposing sensitive data. This paper proposes a client selection framework for FL that integrates differential privacy and fault tolerance. Our adaptive approach dynamically adjusts the number of selected clients based on model performance and system constraints, ensuring privacy through calibrated noise addition. We evaluate our method on a network anomaly detection use case using the UNSW-NB15 and ROAD datasets. Results show up to a 7% increase in accuracy and a 25% reduction in training time compared to FedL2P. Moreover, we highlight the trade-offs between privacy budgets and model performance, with higher privacy budgets reducing noise and improving accuracy. Our fault tolerance mechanism, while causing a slight performance drop, enhances robustness to client failures. Statistical validation using Mann-Whitney U tests confirms the significance of these improvements (p < 0.05).

Marfo, William [University of Texas at El Paso,Dep↗
DOE OSTI2026

Pinned Clouds over Industrial Sources of Heat during TRACER

An analysis of stereo-camera data from the Tracking Aerosol Convection Interactions Experiment (TRACER) campaign in Houston, Texas, reveals the existence of pinned clouds sitting atop industrial heat sources. These are not plumes of vapor and condensate emanating from stacks: The pinned shallow cumuli have cloud bases at about the lifting condensation of near-surface air. On many mornings, the pinned clouds are the only shallow cumuli in the field of view of the stereo cameras, persisting by themselves for at least an hour. On those mornings, the lower atmosphere straddles the boundary between stable and unstable, and the waste heat from industrial facilities is able to pin moist convection overhead. When solar heating of the surface becomes sufficient later in the morning, those mornings transition to having widespread shallow moist convection. Occurring within the field of view of the stereo cameras and other TRACER instruments, these pinned clouds represent steady plume-like moist convection triggered by known heat sources in a well-characterized atmosphere, making them a rich target for further study. Significance Statement This article introduces steady-state cumulus clouds observed during the Tracking Aerosol Convection Interactions Experiment (TRACER) campaign that were pinned to industrial heat sources. These pinned clouds were present in the early morning only when conditions were primed for moist convection, but for which there was not yet enough solar heating to trigger widespread convection. They persist in a statistical steady state and therefore provide a natural setting to explore the mechanisms controlling convective dynamics, which play a central role in Earth’s atmospheric circulation and radiative balance.

Atmosphere↗
NASA NTRS1980

Clear air turbulence forecasting techniques

A method to improve clear air turbulence (CAT) forecasting by more effectively using the currently operational rawinsonde (RW) system is discussed. The method is called the Diagnostic Richardson Number Tendency (DRT) technique. The technique does not attempt to use the RW as a direct detector of the turbulent motion or even of the CAT mechanism structure but rather senses the synoptic scale centers of action which provide the energy to the CAT mechanism at the mesoscale level. The DRT algorithm is deterministic rather than statistical in nature, using the hydrodynamic equations (equations of motion) relevant to the synoptic scale. However, interpretation, by necessity, is probabilistic. What is most important with respect to its operational implementation is that this method uses the same input data as currently used by the operational National Meteorological Center prognostic models.

Keller, J. L.↗
NASA NTRS1986

Rotor noise due to atmospheric turbulence ingestion. I - Fluid mechanics

In the present analytical procedure for the prediction of helicopter rotor noise generation due to the ingestion of atmospheric turbulence, different models for turbulence fluid mechanics and the ingestion process are combined. The mean flow and turbulence statistics associated with the atmospheric boundary layer are modeled with attention to the effects of atmospheric stability length, windspeed, and altitude. The turbulence field can be modeled as isotropic, locally stationary, and homogeneous. For large mean flow contraction ratios, accurate predictions of turbulence vorticity components at the rotor face requires the incorporation of the differential drift of fluid particles on adjacent streamlines.

Simonich, J. C.↗
NASA NTRS1992

Low-earth-orbit effects on strength of glasses

The effects of a 5.8-y exposure to low-earth-orbit environment upon the mechanical properties of five commercial glasses and a low-expansion-coefficient glass-ceramic have been examined. The radiation components of the earth-orbit environment did not degrade the mechanical strength of the samples examined within the limits of experimental error. Statistical problems arising from the low frequency of micrometeorite or space debris impacts upon the samples precluded statistically valid measurement of impacted sample strengths. Upper bounds for the magnitude of the impact event damage upon the strengths for impacted samples have been determined.

Wiedlocher, David E.↗
NASA NTRS2002

Insights into Corona Formation through Statistical Analyses

Statistical analysis of an expanded database of coronae on Venus indicates that the populations of Type 1 (with fracture annuli) and 2 (without fracture annuli) corona diameters are statistically indistinguishable, and therefore we have no basis for assuming different formation mechanisms. Analysis of the topography and diameters of coronae shows that coronae that are depressions, rimmed depressions, and domes tend to be significantly smaller than those that are plateaus, rimmed plateaus, or domes with surrounding rims. This is consistent with the model of Smrekar and Stofan and inconsistent with predictions of the spreading drop model of Koch and Manga. The diameter range for domes, the initial stage of corona formation, provides a broad constraint on the buoyancy of corona-forming plumes. Coronae are only slightly more likely to be topographically raised than depressions, with Type 1 coronae most frequently occurring as rimmed depressions and Type 2 coronae most frequently occuring with flat interiors and raised rims. Most Type 1 coronae are located along chasmata systems or fracture belts, while Type 2 coronas are found predominantly as isolated features in the plains. Coronae at hotspot rises tend to be significantly larger than coronae in other settings, consistent with a hotter upper mantle at hotspot rises and their active state.

Glaze, L. S.↗
NASA NTRS2002

Insights into Corona Formation Through Statistical Analyses

Statistical analysis of an expanded database of coronae on Venus indicates that the populations of Type 1 (with fracture annuli) and 2 (without fracture annuli) corona diameters are statistically indistinguishable, and therefore we have no basis for assuming different formation mechanisms. Analysis of the topography and diameters of coronae shows that coronae that are depressions, rimmed depressions, and domes tend to be significantly smaller than those that are plateaus, rimmed plateaus, or domes with surrounding rims. This is consistent with the model of Smrekar and Stofan and inconsistent with predictions of the spreading drop model of Koch and Munga. The diameter range for domes, the initial stage of corona formation, provides a broad constraint on the buoyancy of corona-forming plumes. Coronae are only slightly more likely to be topographically raised than depressions, with Type 1 coronae most frequently occurring as rimmed depressions and Type 2 coronae most frequently occurring with flat interiors and raised rims. Most Type 1 coronae are located along chasmata systems or fracture belts, while Type 2 coronae are found predominantly as isolated features in the plains. Coronae at hot spot rises tend to be significantly lager than coronae in other settings, consistent with a hotter upper mantle at hot spot rises and their active state.

Glaze, L. S.↗
NASA NTRS2010

Fiber Breakage Model for Carbon Composite Stress Rupture Phenomenon: Theoretical Development and Applications

Stress rupture failure of Carbon Composite Overwrapped Pressure Vessels (COPVs) is of serious concern to Science Mission and Constellation programs since there are a number of COPVs on board space vehicles with stored gases under high pressure for long durations of time. It has become customary to establish the reliability of these vessels using the so called classic models. The classical models are based on Weibull statistics fitted to observed stress rupture data. These stochastic models cannot account for any additional damage due to the complex pressure-time histories characteristic of COPVs being supplied for NASA missions. In particular, it is suspected that the effects of proof test could significantly reduce the stress rupture lifetime of COPVs. The focus of this paper is to present an analytical appraisal of a model that incorporates damage due to proof test. The model examined in the current paper is based on physical mechanisms such as micromechanics based load sharing concepts coupled with creep rupture and Weibull statistics. For example, the classic model cannot accommodate for damage due to proof testing which every flight vessel undergoes. The paper compares current model to the classic model with a number of examples. In addition, several applications of the model to current ISS and Constellation program issues are also examined.

Murthy, Pappu L. N.↗
NASA NTRS2006

A Comparison of Wind Speed Data from Mechanical and Ultrasonic Anemometers

This study compared the performance of mechanical and ultrasonic anemometers at the Eastern Range (ER; Kennedy Space Center and Cape Canaveral Air Force Station on Florida's Atlantic coast) and the Western Range (WR; Vandenberg Air Force Base on California's Pacific coast). Launch Weather Officers, forecasters, and Range Safety analysts need to understand the performance of wind sensors at the ER and WR for weather warnings, watches, advisories, special ground processing operations, launch pad exposure forecasts, user Launch Commit Criteria (LCC) forecasts and evaluations, and toxic dispersion support. The current ER and WR weather tower wind instruments are being changed from the current propeller-and-vane (ER) and cup-and-vane (WR) sensors to ultrasonic sensors through the Range Standardization and Automation (RSA) program. The differences between mechanical and ultrasonic techniques have been found to cause differences in the statistics of peak wind speed in previous studies. The 45th Weather Squadron (45 WS) and the 30th Weather Squadron (30 WS) requested the Applied Meteorology Unit (AMU) to compare data between RSA and current sensors to determine if there are significant differences. Approximately 3 weeks of Legacy and RSA wind data from each range were used in the study, archived during May and June 2005. The ER data spanned the full diurnal cycle, while the WR data was confined to 1000-1600 local time. The sample of 1-minute data from numerous levels on 5 different towers on each range totaled more than 500,000 minutes of data (482,979 minutes of data after quality control). The 10 towers were instrumented at several levels, ranging from 12 ft to 492 ft above ground level. The RSA sensors were collocated at the same vertical levels as the present sensors and typically within 15 ft horizontally of each another. Data from a total of 53 RSA ultrasonic sensors, collocated with present sensors were compared. The 1-minute average wind speed/direction and the 1-second peak wind speed/direction were compared.

Short, D.↗
DOE OSTI2023

On the convergence of statistics in simulations of stationary incompressible turbulent flows

When reporting statistics from simulations of statistically stationary chaotic phenomenon, it is important to verify that the simulations are time-converged. This condition is connected with the statistical error or number of digits with which statistics can be reliably reported. In this work we consider numerical experiments of low Reynolds number incompressible homogeneous and isotropic turbulence as a model problem to investigate statistical convergence over finite simulation times. Specifically, we investigate the time integration requirements that allow meaningful reporting of the statistical error associated with finiteness of the temporal domain. We address two key questions: (1) How long should a simulation be performed in terms of large eddy time, and (2) How should the simulation time be divided among temporal windows over which a quantity of interest is estimated so that its statistical error could be reliably reported? We find that reliable reporting of statistical errors requires simulations on the order of 10 4 large eddy times, which is orders of magnitude longer than typically performed. Additionally, data post-processing should employ windows of at least ten times the large eddy time scale, with the most robust computation of statistical error of the mean requiring window sizes of an additional factor of ten. For practical simulations, we demonstrate that it is possible to estimate the statistical error within a factor of two under a less stringent condition in which a minimum of four windows with size at least ten large eddy times are used. In conclusion, our observations for homogeneous isotropic turbulence are also shown to hold in turbulent channel flow.

42 ENGINEERING↗
NASA NTRS1992

The compositional structure and the effects of exothermicity in an unpremixed planar jet flame

Results are presented of direct numerical simulations (DNS) of a randomly perturbed compressible, spatially developing planar jet under the influence of a finite rate Arrhenius chemical reaction of the type F + O yields Product + Heat. The objectives of the simulations are to assess the compositional structure of the flame and to determine the influence of reaction exothermicity by means of statistical sampling of the data generated by DNS. It is shown that even with this idealized kinetics model, the simulated results exhibit features in accordance with experimental data. These results indicate that the Damkoehler number is an important parameter in determining the statistical composition of the reacting field, and that the results are insensitive to the mechanism by which this parameter is varied. It is demonstrated that as the intensity of mixing is increased and the effect of finite rate chemistry is more pronounced, the magnitudes of the ensemble mean and rms of the product mass fraction decrease, and those of the reactants mass fraction increase. Also, at higher mixing rates the joint probability density functions of the reactants' mass fractions shift towards higher values within the composition domain indicating a lower reactedness. These trends are consistent with those observed experimentally and are very useful in portraying the statistical structure of non-equilibrium diffusion flames. The DNS generated data are also utilized to examine the applicability of the 'laminar diffusion flamelet model' in predicting the rate of the reactant conversion with finite rate chemistry. This examination indicates that the bounds of the product formation scale reasonably well with those obtained by the flamelet model. However, in the range of the relatively low values of the Damkoehler numbers considered, the scatter of the results is substantially more than that to be completely determined by the closure.

Steinberger, C. J.↗
DOE OSTI2023

Modes and Mechanisms of Pacific Decadal-Scale Variability

The modes of Pacific decadal-scale variability (PDV), traditionally defined as statistical patterns of variance, reflect to first order the ocean's integration (i.e., reddening) of atmospheric forcing that arises from both a shift and a change in strength of the climatological (time-mean) atmospheric circulation. While these patterns concisely describe PDV, they do not distinguish among the key dynamical processes driving the evolution of PDV anomalies, including atmospheric and ocean teleconnections and coupled feedbacks with similar spatial structures that operate on different timescales. In this review, we synthesize past analysis using an empirical dynamical model constructed from monthly ocean surface anomalies drawn from several reanalysis products, showing that the PDV modes of variance result from two fundamental low-frequency dynamical eigenmodes: the North Pacific–central Pacific (NP-CP) and Kuroshio–Oyashio Extension (KOE) modes. Both eigenmodes highlight how two-way tropical–extratropical teleconnection dynamics are the primary mechanisms energizing and synchronizing the basin-scale footprint of PDV. While the NP-CP mode captures interannual- to decadal-scale variability, the KOE mode is linked to the basin-scale expression of PDV on decadal to multidecadal timescales, including contributions from the South Pacific.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Statistical Homogenization of Elastic and Fracture Properties of a Sample Selective Laser Melting Material

Selective laser melting (SLM) is an additive manufacturing technique commonly used in the rapid prototyping of components. The complexity of the SLM microstructure poses a unique challenge to deriving effective mechanical properties at different length scales. Representative volume elements (RVEs) are often used to homogenize the material properties of composites. Instead of RVEs, we use statistical volume elements (SVEs) to homogenize the elastic and fracture properties of the material. This relates the inherent variation of a material’s microstructure to the variation in its mechanical properties at different observation scales. The convergence to the RVE limit is examined from two perspectives: the stability of the mean value as the SVE size increases for the mean-based approach, and the tendency of the normalized variation in homogenized properties to zero as the SVE size increases for the variation-based approach. Fracture properties tend to make the RVE limit slower than do elastic properties from both perspectives. There are also differences between vertical (normal to printing plane) and horizontal (in-plane) properties. While the elastic properties tend to make the RVE limit faster for the horizontal direction, i.e., having a smaller variation and more stable mean value, the fracture properties exhibit the opposite effect. We attributed these differences to the geometry of the melt pools.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Quantitative measurements of dislocations in metals for advancing predictive simulations

LLNL applications require scientists to predict how materials evolve under various thermomechanical conditions. While this is achieved through physics-based simulations, uncertainty in the predictions of mechanical properties remains a serious challenge that limits the predictive capabilities of models because we lack methods to compare predictions of atomic-scale defects (dislocations) with experimental measurements. High energy X-ray diffraction (HEXRD) is the most relevant technique that can provide the necessary statistical information on dislocations. However, this technique is not yet quantitative because we lack a precise understanding of the relationship between X-ray diffraction patterns and the underlying material dislocation content and arrangements. To address this need, we used our novel computational X-ray diffraction method to simulate the effect of dislocations on the diffraction patterns. We compared virtual and experimental diffraction patterns. Results allowed us to clearly establish the relationship between X-ray diffraction patterns and the underlying dislocation structures, proving that it is feasible to quantitatively measure dislocation statistics with HEXRD. This project delivered a method that can provide the missing piece to LLNL’s mechanical property simulations in advanced metals by obtaining experimentally long-needed quantitative dislocation data, which could fully enable predictive capabilities.

36 MATERIALS SCIENCE↗