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 217 records · Page 12

Computing Thermodynamic Properties of Fluids Augmented by Nanoconfinement: Application to Pressurized Methane

Nanoconfined fluids exhibit remarkably different thermodynamic behavior compared to the bulk phase. These confinement effects render predictions of thermodynamic quantities of nanoconfined fluids challenging. In particular, confinement creates a spatially varying density profile near the wall that is primarily responsible for adsorption and capillary condensation behavior. Significant fluctuations in thermodynamic quantities, inherent in such nanoscale systems, coupled to strong fluid–wall interactions give rise to this near-wall density profile. Empirical models have been proposed to explain and model these effects, yet no first-principles based formulation has been developed. We present a statistical mechanics framework that embeds such a coupling to describe the effect of the fluid–wall interaction in amplifying the near-wall density behavior for compressible gases at elevated pressures such as pressurized methane in confinement. We show that the proposed theory predicts accurately the adsorbed layer thickness as obtained with small-angle neutron scattering measurements. Furthermore, the predictions of density under confinement from the proposed theory are shown to be in excellent agreement with available experimental and atomistic simulations data for a range of temperatures for nanoconfined methane. While the framework is presented for evaluating the near-wall density, owing to its rigorous foundation in statistical mechanics, the proposed theory can also be generalized for predicting phase-transition and nonequilibrium transport of nanoconfined fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lieb-Robinson Light Cone for Power-Law Interactions

The Lieb-Robinson theorem states that information propagates with a finite velocity in quantum systems on a lattice with nearest-neighbor interactions. What are the speed limits on information propagation in quantum systems with power-law interactions, which decay as 1 / r α at distance r ? Here, we present a definitive answer to this question for all exponents α > 2 d and all spatial dimensions d . Schematically, information takes time at least r min { 1 , α - 2 d } to propagate a distance r . As recent state transfer protocols saturate this bound, our work closes a decades-long hunt for optimal Lieb-Robinson bounds on quantum information dynamics with power-law interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamics of particle network in composite battery cathodes

We report that improving composite battery electrodes requires a delicate control of active materials and electrode formulation. The electrochemically active particles fulfill their role as energy exchange reservoirs through interacting with the surrounding conductive network. We formulate a network evolution model to interpret the regulation and equilibration between electrochemical activity and mechanical damage of these particles. Through statistical analysis of thousands of particles using x-ray phase contrast holotomography in a LiNi 0.8 Mn 0.1 Co 0.1 O 2 -based cathode, we found that the local network heterogeneity results in asynchronous activities in the early cycles, and subsequently the particle assemblies move toward a synchronous behavior. Our study pinpoints the chemomechanical behavior of individual particles and enables better designs of the conductive network to optimize the utility of all the particles during operation.

25 ENERGY STORAGE↗

Analysis of atmospheric flow over a surface protrusion using the turbulence kinetic energy equation

Atmospheric flow fields resulting from a semi-elliptical surface obstruction in an otherwise horizontally homogeneous statistically stationary flow are modelled with the boundary-layer/Boussinesq-approximation of the governing equation of fluid mechanics. The turbulence kinetic energy equation is used to determine the dissipative effects of turbulent shear on the mean flow. Mean-flow results are compared with those given in a previous paper where the same problem was attacked using a Prandtl mixing-length hypothesis. Iso-lines of turbulence kinetic energy and turbulence intensity are plotted in the plane of the flow. They highlight regions of high turbulence intensity in the stagnation zone and sharp gradients in intensity along the transition from adverse to favourable pressure gradient.

Frost, W.↗

Analysis of atmospheric flow over a surface protrusion using the turbulence kinetic energy equation with reference to aeronautical operating systems

Flow over surface obstructions can produce significantly large wind shears such that adverse flying conditions can occur for aeronautical systems (helicopters, STOL vehicles, etc.). Atmospheric flow fields resulting from a semi-elliptical surface obstruction in an otherwise horizontally homogeneous statistically stationary flow are modelled with the boundary-layer/Boussinesq-approximation of the governing equation of fluid mechanics. The turbulence kinetic energy equation is used to determine the dissipative effects of turbulent shear on the mean flow. Iso-lines of turbulence kinetic energy and turbulence intensity are plotted in the plane of the flow and highlight regions of high turbulence intensity in the stagnation zone and sharp gradients in intensity along the transition from adverse to favourable pressure gradient. Discussion of the effects of the disturbed wind field in CTOL and STOL aircraft flight path and obstruction clearance standards is given. The results indicate that closer inspection of these presently recommended standards as influenced by wind over irregular terrains is required.

Frost, W.↗

Physical concepts in the development of constitutive equations

Proposed viscoplastic material models include in their formulation observed material response but do not generally incorporate principles from thermodynamics, statistical mechanics, and quantum mechanics. Numerous hypotheses were made for material response based on first principles. Many of these hypotheses were tested experimentally. The proposed viscoplastic theories and the experimental basis of these hypotheses must be checked against the hypotheses. The physics of thermodynamics, statistical mechanics and quantum mechanics, and the effects of defects, are reviewed for their application to the development of constitutive laws.

Cassenti, B. N.↗

Some design considerations for ceramic components in heat engine applications

The design methodology for brittle material structures which is being developed and used at the Lewis Research Center for sizing ceramic components in heat engine applications is reviewed. Theoretical aspects of designing with structural ceramics are discussed, and a general purpose reliability program for predicting fast fracture response due to volume distributed flaws is described. Statistical treatment of brittle behavior, based on the Weibull model, is reviewed and its advantages, as well as drawbacks, are listed. A mechanistic statistical fracture theory, proposed by Batdorf to overcome the Weibull model limitations and based on Griffith fracture mechanics, is summarized. Failure probability predictions are made for rotating annular Si3N4 disks using various fracture models, and the results are compared to actual failure data. The application of these design methods to Government funded ceramics engine demonstration programs is surveyed. The uncertainty in observed component performance emphasizes the need for proof testing and improved nondestructive evaluation to guarantee adequate structural integrity.

Gyekenyesi, John P.↗

Visualizing characteristics of ocean data collected during the Shuttle Imaging Radar-B experiment

Topographic measurements of sea surface elevation collected by the Surface Contour Radar (SCR) during NASA's Shuttle Imaging Radar (SIR-B) experiment are plotted as three dimensional surface plots to observe wave height variance along the track of a P-3 aircraft. Ocean wave spectra were computed from rotating altimeter measurements acquired by the Radar Ocean Wave Spectrometer (ROWS). Fourier power spectra computed from SIR-B synthetic aperture radar (SAR) images of the ocean are compared to ROWS surface wave spectra. Fourier inversion of SAR spectra, after subtraction of spectral noise and modeling of wave height modulation, yields topography similar to direct measurements made by SCR. Visual perspectives on the SCR and SAR ocean data are compared. Threshold distinctions between surface elevation and texture modulations of SAR data are considered within the context of a dynamic statistical model of rough surface scattering. The result of these endeavors is insight as to the physical mechanism governing the imaging of ocean waves with SAR.

Tilley, David G.↗

Simulation of the Onset of the Southeast Asian Monsoon during 1997 and 1998: The Impact of Surface Processes

The onset of the southeast Asian monsoon during 1997 and 1998 was simulated by coupling a mesoscale atmospheric model (MM5) and a detailed, land surface model, PLACE (the Parameterization for Land-Atmosphere-Cloud Exchange). The rainfall results from the simulations were compared with observed satellite data from the TRMM (Tropical Rainfall Measuring Mission) TMI (TRMM Microwave Imager) and GPCP (Global Precipitation Climatology Project). The control simulation with the PLACE land surface model and variable sea surface temperature captured the basic signatures of the monsoon onset processes and associated rainfall statistics. Sensitivity tests indicated that simulations were sigmficantly improved by including the PLACE land surface model. The mechanism by which the land surface processes affect the moisture transport and the convection during the onset of the southeast Asian monsoon were analyzed. The results indicated that land surface processes played an important role in modifying the low-level wind field over two major branches of the circulation: the southwest low-level flow over the Indo-china peninsula and the northern, cold frontal intrusion from southern China. The surface sensible and latent heat fluxes modified the low-level temperature distribution and gradient, and therefore the low-level wind due to the thermal wind effect. The more realistic forcing of the sensible and latent heat fluxes from the detailed, land surface model improved the low-level wind simulation apd associated moisture transport and convection.

Wang, Yansen↗

Simulation of the Onset of the Southeast Asian Monsoon during 1997 and 1998: The Impact of Surface Processes

The onset of the southeast Asian monsoon during 1997 and 1998 was simulated by coupling a mesoscale atmospheric model (MM5) and a detailed, land surface model, PLACE (the Parameterization for Land-Atmosphere-Cloud Exchange). The rainfall results from the simulations were compared with observed satellite data from the TRMM (Tropical Rainfall Measuring Mission) TMI (TRMM Microwave Imager) and GPCP (Global Precipitation Climatology Project). The control simulation with the PLACE land surface model and variable sea surface temperature captured the basic signatures of the monsoon onset processes and associated rainfall statistics. Sensitivity tests indicated that simulations were significantly improved by including the PLACE land surface model. The mechanism by which the land surface processes affect the moisture transport and the convection during the onset of the southeast Asian monsoon were analyzed. The results indicated that land surface processes played an important role in modifying the low-level wind field over two major branches of the circulation: the southwest low-level flow over the Indo-China peninsula and the northern, cold frontal intrusion from southern China. The surface sensible and latent heat fluxes modified the low-level temperature distribution and merit, and therefore the low-level wind due to the thermal wind effect. The more realistic forcing of the sensible and latent heat fluxes from the detailed, land surface model improved the low-level wind simulation and associated moisture transport and convection.

Wang, Yansen↗

Benchmarking Microscale Ductility Measurements (Final Report of the Project DE-NE0008799)

Conventional macroscale experimentation is generally considered to be straightforward with few limitations. Conversely, micro/nanoscale experimentation presents numerous challenges in loading device design, sample preparation and handling, as well as accurate understanding of grain size and local texture effects on recorded measurements. Despite these challenges, nanopillar compression, MEMs based micro-tension, and nanoindentation approaches have been able to provide fundamental contributions to the understanding of material behavior at small lengthscales. However, the overarching shortcoming of these micro/nanoscale experimentation approaches, is the inability to directly translate measurements evaluated at the nm and µm length scales (e.g., hardness) to macroscale tensile material behavior (i.e., elastic modulus, yield strength, and ductility). The objectives of the proposed study are, 1) to establish best practices for obtaining tensile microscale ductility measurements, and 2) to validate methodologies to for comparing microscale ductility measurements to macroscale ductility measurements. In order to achieve these objectives, a multi-lengthscale, multi-temperature testing protocol and simulation framework are executed first on copper as a model material to validate the following approach, and second on reactor grade Zircaloy-2. Experiments are conducted on specimens extracted from the same test piece to ensure nominally identical grain size and texture from specimens to specimen. Motivated by the need to isolate the contribution of size-effects on obtained mechanical property measurements, specimens are manufactured with thicknesses at the micro- (1-10 µm), meso- (10-100's µm), and macroscales (sub-sized ASTM E8). In-situ full-field deformation techniques (scanning electron microscopy (SEM) grid methods and optical DIC) are incorporated into testing at each specimen length-scale to capture plasticity localization and evolution. Experimental testing for all specimens is conducted at both room temperature and elevated temperatures to probe the role of thermal activation on plastic deformation accommodation processes. Simulation efforts focus on examining the mechanical behavior of microscale specimens using a finite element approach with explicitly resolved grain morphologies, and an embedded crystal plasticity model. The cost-efficient implementation method allows for the modeling of a statistically significant number of both real (i.e., digital twin) and generated microstructures to obtain an understanding of the interrelationships between specimen microstructure and geometric variables (grain size, texture, specimen geometry, etc.) on microscale mechanical behavior.

36 MATERIALS SCIENCE↗

pyTCR: A tropical cyclone rainfall model for python

pyTCR is a climatology software package developed in the Python programming language. It integrates the capabilities of several legacy physical models and increases computational efficiency to allow rapid estimation of tropical cyclone (TC) rainfall consistent with the large-scale environment. Specifically, pyTCR implements a horizontally distributed and vertically integrated model [Zhu et al., 2013] for simulating rainfall driven by TCs. Along storm tracks, rainfall is estimated by computing the cross-boundary-layer, upward water vapor transport caused by different mechanisms including frictional convergence, vortex stretching, large-scale baroclinic effect (i.e., wind shear), topographic forcing, and radiative cooling [Lu et al., 2018]. The package provides essential functionalities for modeling and interpreting spatio-temporal TC rainfall data. pyTCR requires a limited number of model input parameters, making it a convenient and useful tool for analyzing rainfall mechanisms driven by TCs. To sample rare (most intense) rainfall events that are often of great societal interest, pyTCR adapts and leverages outputs from a statistical-dynamical TC downscaling model [Lin et al., 2023] capable of rapidly generating a large number of synthetic TCs given a certain climate. As a result, pyTCR significantly reduces computational effort and improves the efficiency in capturing extreme TC rainfall events at the tail of the distributions from limited datasets. Furthermore, the TC downscaling model is forced entirely by large-scale environmental conditions from reanalysis data or coupled General Circulation Models (GCMs), simplifying the projection of TC-induced rainfall and wind speed under future climate using pyTCR. Finally, pyTCR can be coupled with hydrological and wind models to assess risks associated with independent and compound events (e.g., storm surges and freshwater flooding).

54 ENVIRONMENTAL SCIENCES↗

A Multiphysics Multiscale Simulation Platform for Damage, Environmental Degradation, and Life Prediction of CMCs in Extreme Environments

This project successfully developed a multiphysics, multiscale computational framework to enhance the design and development of CMCs, with a focus on modeling highly nonlinear, time-dependent damage mechanisms and material degradation under extreme conditions, such as those experienced in turbine service environments. The project made significant advances in improving our understanding of progressive damage, oxidative degradation, and time-dependent inelastic deformation in CMCs, with particular attention to the role of uncertainties in predictions. Key outcomes include the integration of advanced material characterization, uncertainty quantification, and multiphysics constitutive models to predict the behavior of CMCs over their service life. A novel multiscale methodology was employed, which integrated microscale constituent behaviors with structural-scale responses, enabling the manufacturing defects in the microstructure that are prone to damage nucleation. Through the development of DL algorithms, the project advanced the prediction of damage initiation and crack propagation, taking into account the defect morphology and statistical variations across multiple scales. The framework was rigorously validated using thermomechanical experiments, which tested CMCs under various mechanical loadings at elevated temperatures, further enhancing the model's predictive capability. Overall, the research outcomes have provided a more accurate, reliable method for predicting CMC component life, significantly advancing material design, and improving component reliability in extreme environments. This work has strong implications for the optimization of turbine components and other high-performance applications where CMCs are used.

03 NATURAL GAS↗

The endpoint of partial deconfinement

We study the matrix quantum mechanics of two free hermitian N × N matrices subject to a singlet constraint in the microcanonical ensemble. This is the simplest example of a theory that at large N has a confinement/deconfinement transition. In the microcanonical ensemble, it also exhibits partial deconfinement with a Hagedorn density of states. We argue that the entropy of these configurations, based on a combinatorial counting of Young diagrams, are dominated by Young diagrams that have the VKLS shape. When the shape gets to the maximal depth allowed for a Young diagram of SU(N), namely N, we argue that the system stops exhibiting the Hagedorn behavior. The number of boxes (energy) at the transition is N 2 /4, independent of the charge of the state.

1/N Expansion↗

Research on decreasing the noise of a jet

An investigation was conducted to determine the amount of noise reduction that could be obtained by communicating the axial vibrations of the nozzle of a jet engine to the suspension mechanism. Mathematical models were developed to establish the relationship of the parameters. The pressure pulsations in the jet were recorded by measuring microphones and the data were processed by statistical dynamics methods. It was determined that the jet noise spectrum was successfully changed with various mechanical suspension impedances.

Akelis, M. E.↗

Simultaneous Development and Robust Optimization of a Microstructure Dependent Material

Recent microstructure characterization techniques combined with Symbolic Regression(SR)analysis has been proven to generate white box plasticity models well suited for incorporation into FEA software.The current work builds upon those efforts and demonstrates the applicability of Sequential Monte-Carlo (SMC) methods within SR analysis to condense model development and robust optimization into a single, co-dependent process. In this project, SMC methods provide a mechanism through which the observed microstructure features and associated variability can be incorporated into the discovery phase of model development and simultaneously recover approximate parameter distributions through SR analysis. The demonstration utilized a data set consisting of tensile test results from a limited number of sample specimens with corresponding EBSD data from which microstructure features were characterized.The maximum threshold stress model in the Visco-Plastic Self-Consistent (VPSC) code developed by Los Alamos National Laboratories was calibrated using mechanical test data.Synthetic volume elements with statistically equivalent microstructure were generated with DREAM3Dbased on the observed EBSD data. VPSC was used to simulate the corresponding tensile test response for each of the synthetic volume elements. The simulated microstructure and tensile test data was used astraining datafor SMC-SR algorithm and the resulting model was validated with data from the original empirical data set.

Karl Garbrecht↗

Learning to Predict Crystal Plasticity at the Nanoscale: Deep Residual Networks and Size Effects in Uniaxial Compression Discrete Dislocation Simulations

The density and configurational changes of crystal dislocations during plastic deformation influence the mechanical properties of materials. These influences have become clearest in nanoscale experiments, in terms of strength, hardness and work hardening size effects in small volumes. The mechanical characterization of a model crystal may be cast as an inverse problem of deducing the defect population characteristics (density, correlations) in small volumes from the mechanical behavior. In this work, we demonstrate how a deep residual network can be used to deduce the dislocation characteristics of a sample of interest using only its surface strain profiles at small deformations, and then statistically predict the mechanical response of size-affected samples at larger deformations. As a testbed of our approach, we utilize high-throughput discrete dislocation simulations for systems of widths that range from nano- to micro- meters. We show that the proposed deep learning model significantly outperforms a traditional machine learning model, as well as accurately produces statistical predictions of the size effects in samples of various widths. By visualizing the filters in convolutional layers and saliency maps, we find that the proposed model is able to learn the significant features of sample strain profiles.

36 MATERIALS SCIENCE↗

Sterilization of N95 Respirators via Gamma Radiation: Comparison of Post-sterilization Efficacy

This study evaluated gamma irradiation for sterilization and reuse of two models of N95 respirators after gamma radiation sterilization as a method to increase availability of N95 respirators during a shortage. The Sandia National Laboratories Gamma Irradiation Facility was used to irradiate two different models of N95 filtering facepiece respirators at doses ranging from 0 kGy(tissue) to 50 kGy(tissue). The following tests were used to determine the efficacy of the respirator after irradiation sterilization: Ambient Aerosol Condensation Nuclei Counter Quantitative Fit Test, tensile test, strain cycling, oscillatory dynamic mechanical analysis, microscopic image analysis of fiber layers, and electrostatic field measurements. Both of the respirator models exhibited statistically significant changes after gamma irradiation as shown by the Quantitative Fit Test, electrostatic testing and the aerosol testing. The change in electrostatic capability of the filter reduced the efficiency of challenging particles near the 200 nm size by approximately 40-50%. Both tested respirators showed statistically significant changes associated with gamma sterilization. However, our results indicate that choices in materials and manufacturing methods to achieve N95 filtration lead to different magnitudes of damage when exposed to gamma radiation at sterilization relevant doses. This damage results in lower filtration performance. While our sample size (2 different types of respirators) was small, we did observe a change in electrostatic properties on a filter layer that coincided with the failure on the Quantitative Fit Test.

61 RADIATION PROTECTION AND DOSIMETRY↗