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

The relation between lignin sequence and its 3D structure

Background: Lignin, the second most abundant biopolymer on earth, plays a major structural role in plants, conferring mechanical strength and regulating water conduction. Understanding the three-dimensional structure of lignin is important for fundamental reasons as well as engineering plants towards lignin valorization. Lignin lacks a specific primary sequence, making its average chemical composition the focus of most recent studies. However, it remains unclear whether the 3D structure of lignin molecules depends on their sequence. Methods: In this work, we performed all-atom molecular dynamics simulation of three S/G-lignin molecules with the same average composition but different sequence. Results: A detailed statistical analysis of the radius of gyration and relative shape anisotropy reveals that the lignin sequence has no statistically significant effect on the global three-dimensional structure. We found however, that homopolymers of C-lignin with the same molecular weight have smaller radii of gyration than S/G-lignin. We attribute this to lower hydroxyl content of C-lignin, which makes it more compact and rigid. Conclusions: The 3D structure of lignin is influenced by the overall content of monomeric units and interunit linkages and not by its precise primary sequence. General Significance: Lignin is assumed to not have a well-defined primary structure. The results presented here demonstrate there are no significant differences in the global 3D structure of lignin molecules with the same average composition but different primary sequence.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Polishing the Gold Standard: The Role of Orbital Choice in CCSD(T) Vibrational Frequency Prediction

While CCSD(T) with spin-restricted Hartree-Fock (RHF) orbitals has long been lauded for its ability to accurately describe closed-shell interactions, the performance of CCSD(T) on open-shell species is much more erratic, especially when using a spin-unrestricted HF (UHF) reference. Previous studies have shown improved treatment of open-shell systems when a non-HF set of molecular orbitals, like Brueckner or Kohn-Sham density functional theory (DFT) orbitals, is used as a reference. Inspired by the success of regularized orbital-optimized second-order Møller-Plesset perturbation theory (κ-OOMP2) orbitals as reference orbitals for MP3, we investigate the use of κ-OOMP2 orbitals and various DFT orbitals as reference orbitals for CCSD(T) calculations of the corrected ground-state harmonic vibrational frequencies of a set of 36 closed-shell (29 neutrals, 6 cations, 1 anion) and 59 open-shell diatomic species (38 neutrals, 15 cations, 6 anions). The aug-cc-pwCVTZ basis set is used for all calculations. The use of κ-OOMP2 orbitals in this context alleviates difficult cases observed for both UHF orbitals and OOMP2 orbitals. Removing two multireference systems and 12 systems with ambiguous experimental data leaves a pruned data set. Overall performance on the pruned data set highlights CCSD(T) with a B97 orbital reference (CCSD(T):B97), CCSD(T) with a κ-OOMP2 orbital reference (CCSD(T):κ-OOMP2), and CCSD(T) with a B97M-rV orbital reference (CCSD(T):B97M-rV) with RMSDs of 8.48 cm -1 , and 8.50 cm -1 , and 8.75 cm -1 respectively, outperforming CCSD(T):UHF by nearly a factor of 5. Moreover, the performance on the closed- and open-shell subsets shows these methods are able to treat open-shell and closed-shell systems with comparable accuracy and robustness. CCSD(T) with RHF orbitals is seen to improve upon UHF for the closed-shell species, while spatial symmetry breaking in a number of restricted open-shell HF (ROHF) references leads CCSD(T) with ROHF reference orbitals to exhibit the poorest statistical performance of all methods surveyed for open-shell species. The use of κ-OOMP2 orbitals has also proven useful in diagnosing multireference character that can hinder the reliability of CCSD(T).

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

Hall conductance and the statistics of flux insertions in gapped interacting lattice systems

We study charge transport for zero-temperature infinite-volume gapped lattice systems in two dimensions with short-range interactions. We show that the Hall conductance is locally computable and is the same for all systems that are in the same gapped phase. Here, we provide a rigorous version of Laughlin’s flux-insertion argument, which shows that for short-range entangled systems, the Hall conductance is an integer multiple of e 2 / h . We show that the Hall conductance determines the statistics of flux insertions. For bosonic short-range entangled systems, this implies that the Hall conductance is an even multiple of e 2 / h . Finally, we adapt a proof of quantization of the Thouless charge pump to the case of infinite-volume gapped lattice systems in one dimension.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Compression behavior of microcrystalline cellulose spheres: Single particle compression and confined bulk compression across regimes

Particle characteristics can drastically influence the process-structure-property-performance aspects of granular materials in compression. In this study, we aim to computationally simulate the mechanical processes of stress redistribution in compacts including the kinematics of particle rearrangement during densification and particle deformation leading to fragmentation. Confined compression experiments are conducted with three sets of commercial microcrystalline cellulose particles nearly spherical in shape with different mean particle size. Experimentally measured compression curves from tall powder columns are fitted with the Kenkre et al. (J. of American Chemical Society, Vol. 79, No. 12) model. This model provides a basis to derive several common two-parameter literature models and as a framework to incorporate statistical representations of critical particle behaviors. We focus on the low-stress compression data and the model comparisons typically not discussed in the literature. Furthermore, additional single particle compressions report fracture strength with particle size for comparison to the apparent particle strength extracted from bulk compression data.

42 ENGINEERING↗
DOE OSTI2020

Pass-efficient methods for compression of high-dimensional turbulent flow data

The future of high-performance computing, specifically on future Exascale computers, will presumably see memory capacity and bandwidth fail to keep pace with data generated, for instance, from massively parallel partial differential equation (PDE) systems. Current strategies proposed to address this bottleneck entail the omission of large fractions of data, as well as the incorporation of in situ compression algorithms to avoid overuse of memory. To ensure that post-processing operations are successful, this must be done in a way that a sufficiently accurate representation of the solution is stored. Moreover, in situations where the input/output system becomes a bottleneck in analysis, visualization, etc., or the execution of the PDE solver is expensive, the number of passes made over the data must be minimized. In the interest of addressing this problem, this work focuses on the utility of pass-efficient, parallelizable, low-rank, matrix decomposition methods in compressing high-dimensional simulation data from turbulent flows. Additionally, a particular emphasis is placed on using coarse representation of the data – compatible with the PDE discretization grid – to accelerate the construction of the low-rank factorization. This includes the presentation of a novel single-pass matrix decomposition algorithm for computing the so-called interpolative decomposition. The methods are described extensively and numerical experiments on two turbulent channel flow data are performed. In the first (unladen) channel flow case, compression factors exceeding 400 are achieved while maintaining accuracy with respect to first- and second-order flow statistics. In the particle-laden case, compression factors of 100 are achieved and the compressed data is used to recover particle velocities. These results show that these compression methods can enable efficient computation of various quantities of interest in both the carrier and disperse phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2023

Large Eddy Simulation of Rotationally Induced Ingress and Egress around an Axial Seal between Rotor and Stator Disks

In gas turbines, the hot gas exiting the combustor can have temperatures as high as 2000 °C, and some of this hot gas enter into the space between the stator and rotor disks (wheelspace). Since the entering hot gas could damage the disks, its ingestion must be minimized. This is carried out by rim seals and by introducing a cooler flow from the compressor (sealing flow) into the wheelspace. Ingress and egress into rim seals are driven by the stator vanes, the rotor and its rotation, and the rotor blades. This study focuses on the ingress and egress driven by the rotor and its rotation. This is carried out by performing wall-resolved large eddy simulation (LES) around an axial seal in a rotor–stator configuration without vanes and blades. Results obtained show the mechanisms by which the rotor and its rotation induce ingress, egress, and flow trajectories. Kelvin–Helmholtz instability was found to create a wavy shear layer and displacement thickness that produces alternating regions of high and low pressures around the rotor side of the seal. Vortex shedding on the backward-facing side of the seal and its impingement on the rotor side of the seal also produces alternating regions of high and low pressures. The locations of the alternating regions of high and low pressures were found to be statistically stationary and to cause ingress to start on the rotor side of the seal. Vortex shedding and recirculating flow in the seal clearance also cause ingress by entrainment. With the effects of the rotor and its rotation on ingress and egress isolated, this study enables the effects of stator vanes and rotor blades to be assessed.

42 ENGINEERING↗
DOE OSTI2025

A Standardized Analysis Process Using Digital Image Correlation to Calculate In Situ Cladding Strain from Modified Burst Tests for Fuel Performance Code Validation

Historical data collection on nuclear fuel cladding materials has focused on generating a statistically significant amount of data to assess the material and its failure behavior. Furthermore, data generated to support material model and failure criteria development were previously posttest evaluations, so a large number of tests was required to gain new understanding. A way to expedite this process is to develop techniques capable of generating large, high-fidelity data sets from a single test with lower uncertainty or quantified uncertainty. One such example of this approach is Oak Ridge National Laboratory’s use of modified burst tests (MBTs) to analyze the mechanical behavior and failure conditions of cladding during a simulated reactivity-initiated accident (RIA). Each test incorporates digital image correlation (DIC) analysis techniques that are used to assess the accumulated strain in situ, as well as eventual cladding failure. This work has been fruitful in defining strain-to-failure conditions for materials like silicon carbide (SiC) fiber–reinforced/SiC matrix composite tubes (SiC/SiC), iron-chromium-aluminum (FeCrAl) alloy tubes, and chromium-coated Zircaloy-4 tubes. However, there are numerous DIC software available, including open-source and proprietary software. The different DIC software use various algorithms to process images and calculate displacement values. Using these different software and algorithms can lead to varying results, and perhaps larger-than-expected uncertainties. In the present study, previously published MBT data encompassing a variety of test conditions were reanalyzed with two different DIC software to assess the variance in the calculated strain results. The data consisted of SiC/SiC, FeCrAl, and chromium-coated Zircaloy-4 tubes. Plots of the calculated strains during the transient revealed good agreement between the two DIC software. The average root-mean-square errors between the two software was 0.20% strain, which is slightly larger than a previously reported error value for these tests. In conclusion, this variance in results is low enough that this analysis method can be used for code validation.

Reactivity-initiated accident↗
DOE OSTI2024

Surface Pressure Fluctuations Induced by a Hypersonic Turbulent Boundary Layer on a Sharp Cone at Angle of Attack

High-fidelity simulations are performed to characterize the turbulence-induced wall pressure fluctuations on a sharp cone at a 5.5° angle-of-attack in a Mach 8 flow. Wall-resolved large-eddy simulation (LES) and wall-modeled large-eddy simulation (WMLES) results are compared to measurements at several locations on the cone body, where comparisons with high-frequency PCB sensors are good, while comparison with low-frequency kulite sensors varies depending on the location. For a given streamwise location, simulation results show significant azimuthal variation in the wall pressure fluctuation statistics. Comparisons between LES and WMLES results indicate that WMLES can accurately reproduce the low-frequency component of the autospectra quite well, while slight sensitivity is apparent in the coherence. A grid sensitivity study for WMLES further demonstrated sensitivity in the wall pressure fluctuation coherence, but not in the autospectra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Federated Learning for Efficient Condition Monitoring and Anomaly Detection in Industrial Cyber-Physical Systems

Detecting and localizing anomalies in cyber-physical systems (CPS) has become increasingly challenging as systems grow in complexity, particularly due to varying sensor reliability and node failures in distributed environments. While federated learning (FL) offers a foundation for distributed model training, existing approaches lack mechanisms to handle these CPS-specific challenges. This paper presents an enhanced FL framework that introduces three key innovations: adaptive model aggregation based on sensor reliability, dynamic node selection for resource optimization, and Weibull-based checkpointing for fault tolerance. Our framework enables reliable condition monitoring while addressing the computational and reliability challenges of industrial CPS deployments. Experiments on NASA Bearing and Hydraulic System Datasets demonstrate superior performance over state-of-the-art FL methods, achieving 99.5% AUC-ROC in anomaly detection and maintaining accuracy under node failures. Statistical validation using Mann-Whitney (U) test confirms significant improvements (p < 0.05) in both detection accuracy and computational efficiency across diverse operational scenarios.1

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

Hadronic structure in high-energy collisions

Here, Parton distribution functions (PDFs) describe the structure of hadrons as composed of quarks and gluons. They are needed to make predictions for short-distance processes in high-energy collisions and are determined by fitting to cross-section data. Definitions of the PDFs and their relations to high-energy cross sections are reviewed. The focus is on the PDFs in protons, but PDFs in nuclei are also discussed. The standard statistical treatment needed to fit the PDFs to data using the Hessian method is reviewed in some detail. Tests are discussed that critically examine whether the needed assumptions are indeed valid. Also presented are some ideas of what one can do in case tests indicate that the assumptions fail.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Boundary-induced classical generalized Gibbs ensemble with angular momentum

We investigate how confinement geometry leads to the emergence of a Generalized Gibbs Ensemble (GGE) in classical systems. Unlike the standard Gibbs ensemble, the GGE includes additional conserved quantities, such as angular momentum, that arise from boundary-induced symmetries. Using analytical arguments based on the maximum entropy principle, we show that circular boundaries preserve angular momentum and drive the system toward a chiral, non-ergodic GGE that violates time-reversal symmetry. This ensemble differs fundamentally from the Gibbs case, producing near-boundary condensation and revealing how geometry alone can alter thermal equilibration. To quantify these effects, we introduce an order parameter measuring deviations from Gibbs behavior and demonstrate that conventional Monte Carlo methods must incorporate angular momentum conservation under such conditions. Our study highlights how geometric constraints shape non-equilibrium statistical ensembles and lead to subtle departures from the Bohr-van Leeuwen theorem. These predictions are validated through detailed simulations of confined classical hard-disk gases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Identifying geological structures through microseismic cluster and burst analyses complementing active seismic interpretation

At the Decatur carbon capture and storage site (IL, USA) CO 2 has been injected from 2011–2014 and from 2017 to present near the base of the Lower Mt. Simon Sandstone saline reservoir, resulting in microseismicity. Microseismicity is mainly located in the basement and distributed in distinct spatial clusters. The lack of significant impedance contrasts within the basement makes the interpretation of active-source seismic reflection data challenging, however, recent reprocessing allowed to resolve faults above and at the top of the basement. These faults generally do not coincide with the location of microseismic events and their continuation to the general depth of the seismic events cannot be assumed. This paper shows how the interpretation of the microseismicity can complement structural interpretations of active-source seismic reflection data. In particular, we analyze clusters and bursts (abrupt increases) of microseismicity, identify unresolved, smaller-scale weaknesses and extract statistical parameters. These parameters allow comparisons with the interpreted faults, and with fracture sets intercepted by boreholes. During injection at the Decatur site, the injection pressure was kept far below fracture pressure, nevertheless, seismic events were induced and spread far beyond the expected extent of the CO 2 plume. We argue that local stress transfers related to the CO 2 injection reactivated pre-existing fractures within the critically stressed basement. Finally, we conducted a slip tendency analysis for faults interpreted from active seismic, selected cluster, bursts and nodal planes from focal mechanisms to determine if the interpreted structures are optimally oriented with respect to the stress regime. Our results suggest that the orientation of fractures close to the injection well, generally shows slight deviations from the optimal orientation for slip. This might indicate either slight local deviations of the maximum horizontal stress azimuth from the average direction used in the analysis, or the lack of optimally oriented fractures at this location.

58 GEOSCIENCES↗
DOE OSTI2022

Tea consumption and risk of bladder cancer in the Bladder Cancer Epidemiology and Nutritional Determinants (BLEND) Study: Pooled analysis of 12 international cohort studies

Tea has been shown to be associated with reduced risk of several diseases including cardiovascular diseases, stroke, metabolic syndrome, and obesity. However, the results on the relationship between tea consumption and bladder cancer are conflicting. This research aimed to assess the association between tea consumption and risk of bladder cancer using a pooled analysis of prospective cohort data. Individual data from 532,949 participants in 12 cohort studies, were pooled for analyses. Cox regression models stratified by study centre was used to estimate hazard ratios (HR) and corresponding 95% CIs. Fractional polynomial regression models were used to examine the dose–response relationship. A higher level of tea consumption was associated with lower risk of bladder cancer incidence (compared with no tea consumption: HR = 0.87, 95% C.I. = 0.77–0.98 for low consumption; HR = 0.86, 95% C.I. = 0.77–0.96 for moderate consumption; HR = 0.84, 95% C.I. = 0.75–0.95 for high consumption). When stratified by sex and smoking status, this reduced risk was statistically significant among men and current and former smokers. In addition, dose–response analyses showed a lower bladder cancer risk with increment of 100 ml of tea consumption per day (HR-increment = 0.97; 95% CI = 0.96–0.98). A similar inverse association was found among males, current and former smokers while never smokers and females showed non-significant results, suggesting potential sex-dependent effect. Higher consumption of tea is associated with reduced risk of bladder cancer with potential interaction with sex and smoking status. Further studies are needed to clarify the mechanisms for a protective effect of tea (e.g. inhibition of the survival and proliferation of cancer cells and anti-inflammatory mechanisms) and its interaction with smoking and sex.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Challenges in quantifying unparameterized spatial uncertainties in deep geologic repositories for nuclear waste

Spatially heterogeneous uncertainties are prevalent in geophysical modeling applications, such as probabilistic post-closure performance assessment (PA) of deep geologic repositories for nuclear waste. Such uncertainties are often highly influential to model outputs, so it is desirable to identify the most important mechanisms by which they influence model predictions. However, these uncertainties are often unparameterized in the sense that there is no set of parameters that can be specified to yield a particular realization of the uncertainty. Additionally, the uncertainty is not intrinsically endowed with a parameterization that captures a realization’s mechanistic influence on model outputs. Therefore, in this work we present a novel methodology to develop and assess a set of proxy variables that aim to represent this influence. We show how they can be computed, downselected, and incorporated into the construction of statistical surrogate models mapping model inputs to outputs. We present our methodology in the context of a motivating application problem in deep geologic repository PA and discuss the challenges in capturing the effects of these spatial heterogeneities in uncertainty analyses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Analysis of Deformation and Fracture Mechanisms in the Harvested Low-Dose Baffle Former Bolt via Advanced Mechanical Tests

This report details the production of tensile specimens from in-service irradiated baffle former bolts (BFBs) and presents the results of a pilot in-situ SEM/EBSD test conducted at Oak Ridge National Laboratory’s Low Activation Materials Development and Analysis (LAMDA) facility. The SEM/EDS results confirmed that the material is AISI 316L steel, not AISI 347, and that no Nb addition was detected in several analyzed locations. EBSD analysis of the microstructure revealed in-grain misorientation gradients indicative of cold work, likely caused by material processing, and no retained ferrite was identified. During straining, dislocation channel formation was the primary deformation mechanism, with twinning also observed in favorably oriented grains at local strains as small as ~1.5%. No strain-induced phase instability was observed, likely due to the fact that the Ni and Cr content reached the upper limits for the 316L steel specification. Multiple microfracture events were observed during the tensile test. Fine strain-induced pores formed at the channel-grain boundary intersection points, while coarser microcracks were associated with non-metallic inclusions. Data analysis is ongoing, and additional tests are planned in the near future to provide more statistics and details.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Identification of Life Models for Li-Ion Batteries Using Penalized Regression and Bilevel Optimization

Reduced-order physics-based life models are extremely useful for rapidly predicting battery state-of-health and for simulating battery lifetime in arbitrary aging conditions. However, identification of well-parameterized models is difficult. This is because, for maximum usefulness in predicting lifetime under a variety of conditions, aging test data exhibits many degradation mechanisms, which all need to be accurately modeled. However, because aging tests are time-consuming and expensive, especially for large-format batteries, a minimum of tests are conducted while probing many stress factors. Building a well-parameterized model is then very challenging: an under-parameterized model will neglect critical degradation modes, and an over-parameterized model will extrapolate poorly to new testing conditions. To complicate this matter, the functional form of the model for any individual degradation rate can be very difficult to identify. In this work, the statistical tools of penalized regression and bilevel optimization are used to help identify both the functional forms of and optimize the parameters of reduced-order life models, accelerating identification of robust models. Model robustness is demonstrated through traditional statistics methods of cross-validation and sensitivity analysis, uncertainty quantification through bootstrap resampling and Monte-Carlo simulation, and simulation of real-world use cases.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Quantifying precursors to void nucleation and coalescence in aluminum

Ductile rupture is a common failure mode for engineering alloys. It is generally comprised of three mechanisms: void nucleation at secondary particles, void growth, and coalescence of voids. Conventional models for ductile rupture have limited precision for predicting macroscopic failure. This is partially because they have been corroborated using ex-situ observations of these three mechanisms and calibrated by macroscopic strain or stress metrics. In this study, in-situ high-energy X-ray characterization was conducted to correlate sites of void nucleation via particle cracking and sites of void growth with grain-scale metrics. Here it is shown that particle cracking is not predicated by elevated stress metrics, which disagrees with classical models. Instead, particle cracking tended to occur for the largest, least spherical particles. This trend persisted on the aggregate-scale and within individual neighborhoods around particles. Furthermore, an Eshelby analysis illuminated a statistically significant increase in maximum principal stress within a cracked particle, compared to an uncracked particle. In-situ observations also revealed two separate examples of intragranular void growth comprised of flat, crack-like features that grew showing alignment with slip systems of the highest resolved shear stress. Post-mortem fractography revealed a secondary population of dimples on these crack-like features, implying that void sheeting may be occurring during this crack-like void growth. Furthermore, these in-situ observations imply that classical models for void growth, that assume a homogenous medium, should be extended to account for the anisotropic grain-level behavior to correctly capture distinct features of void nucleation and growth.

Al-2219↗