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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Evaluating proxies for the drivers of natural gas productivity using machine-learning models

We report the extensive development of unconventional reservoirs using horizontal drilling and multistage hydraulic fracturing has generated large volumes of reservoir characterization and production data. The analysis of this abundant data using statistical methods and advanced machine-learning (ML) techniques can provide data-driven insights into well performance. Most predictive modeling studies have focused on the impact that different well completion and stimulation strategies have on well production but have not fully exploited the available in situ rock property data to determine its role in reservoir productivity. We have used machine-learning techniques to rank rock mechanical properties, microseismic attributes, and stimulation parameters in the order of their significance for predicting natural gas production from an unconventional reservoir. The data for this study came from a hydraulically fractured well in the Marcellus Shale in Monongalia County, West Virginia. The data classes included measurements aggregated by well completion stage that included (1) gas production, (2) well-log-derived measurements including bulk density, elastic moduli, shear impedance, compressional impedance, brittleness, and gamma measurements, (3) microseismic attributes, (4) long-period long-duration (LPLD) event counts, (5) fracture counts, and (6) stimulation parameters that included the fluid injection volume and average pumping pressure. To identify observable proxies for the drivers of gas production, we evaluated five commonly used ML approaches including multivariate adaptive regression spline, Gaussian mixture model, random forest, gradient boosting, and neural network. We selected five variables including LPLD event count, seismogenic b-value, hydraulic diffusivity, cumulative moment, and fluid volume as the features most likely to impact gas productivity at the stage level in the study area. The data-driven selection of these parameters for their importance in determining gas production can help reservoir engineers design more effective hydraulic-fracture treatments in the Marcellus Shale and other similar unconventional reservoirs. Plain language summary: We use machine-learning methods and data-driven selection of reservoir parameters to rank and better understand their importance in determining gas production, which can help reservoir engineers design more effective hydraulic-fracture treatments in the Marcellus Shale and other similar unconventional reservoirs.

58 GEOSCIENCES↗
DOE OSTI2021

Data-Driven Learning for the Mori--Zwanzig Formalism: A Generalization of the Koopman Learning Framework

A theoretical framework which unifies the conventional Mori--Zwanzig formalism and the approximate Koopman learning of deterministic dynamical systems from noiseless observation is presented. In this framework, the Mori--Zwanzig formalism, developed in statistical mechanics to tackle the hard problem of construction of reduced-order dynamics for high-dimensional dynamical systems, can be considered as a natural generalization of the Koopman description of the dynamical system. We next show that, similar to the approximate Koopman learning methods, data-driven methods can be developed for the Mori--Zwanzig formalism with Mori's linear projection operator. We have developed two algorithms to extract the key operators, the Markov and the memory kernel, using time series of a reduced set of observables in a dynamical system. We have adopted the Lorenz `96 system as a test problem and solved for the above operators. These operators exhibit complex behaviors, which are unlikely to be captured by traditional modeling approaches in Mori--Zwanzig analysis. The nontrivial generalized fluctuation-dissipation relationship, which relates the memory kernel with the two-time correlation statistics of the orthogonal dynamics, was numerically verified as a validation of the solved operators. Here we present numerical evidence that the generalized Langevin equation, a key construct in the Mori--Zwanzig formalism, is more advantageous in predicting the evolution of the reduced set of observables than the conventional approximate Koopman operators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Fast strength model characterization using Bayesian statistics

A variety of flow stress models exist with new models constantly being developed. These models aim to approximate the strength of materials in a variety of regimes from quasistatic loading through shock scenarios. All models contain an array of parameters which need to be tuned to the material under study. Some models perform well under limited conditions, requiring adjustment of the parameters when venturing outside of those predefined ranges. Other models perform well over a wide range of conditions with a set of parameters, but may be outperformed by other models optimized on a tighter range of conditions. Recent research by Los Alamos demonstrated the ability to optimize the Johnson Cook (JC) model using a set of 3 plate-impact experiments on Aluminum. They utilized Bayesian statistics and emulation to determine optimal parameters for the model with a quantification of parameter uncertainty. We present an extension of this capability to incorporate velocimetry from plate-impact tests, stress-strain data from split Hopkinson pressure bar and quasistatic compression tests, plus profiles from Taylor cylinders in a unified fashion. Statistically robust comparisons of the performance and uncertainty of different realizations of the JC flow stress model were carried out based on calibration to several possible combinations of these three different experiment types.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Deviations from maximal entanglement for eigenstates of the Sachdev-Ye-Kitaev model

We consider mid-spectrum eigenstates of the Sachdev-Ye-Kiteav (SYK) model. We prove that for subsystems whose size is a constant fraction of the system size, the entanglement entropy deviates from the maximum entropy by at least a positive constant. This result highlights the difference between the entanglement entropy of mid-spectrum eigenstates of the SYK model and that of random states.

Disordered Systems and Neural Networks (cond-mat.d↗
DOE OSTI2025

Photon–photon chemical thermodynamics of frequency conversion processes in highly multimode systems

Abstract Frequency generation in highly multimode nonlinear optical systems is inherently a complex process, giving rise to an exceedingly convoluted landscape of evolution dynamics. While predicting and controlling the global conversion efficiencies in such nonlinear environments has long been considered impossible, here, we formally address this challenge even in scenarios involving a very large number of spatial modes. By utilizing fundamental notions from optical statistical mechanics, we develop a universal theoretical framework that effectively treats all frequency components as chemical reactants/products, capable of undergoing optical thermodynamic reactions facilitated by a variety of multi-wave mixing effects. These photon–photon reactions are governed by conservation laws that directly determine the optical temperatures and chemical potentials of the ensued chemical equilibria for each frequency species. In this context, we develop a comprehensive stoichiometric model and formally derive an expression that relates the chemical potentials to the optical stoichiometric coefficients, in a manner akin to atomic/molecular chemical reactions. This advancement unlocks new predictive capabilities that can facilitate the optimization of frequency generation in highly multimode photonic arrangements, surpassing the limitations of conventional schemes that rely exclusively on nonlinear optical dynamics. Notably, we identify a universal regime of Rayleigh–Jeans thermalization where an optical reaction at near-zero optical temperatures can promote the complete and entropically irreversible conversion of light to the fundamental mode at a target frequency. Our theoretical results are corroborated by numerical simulations in settings where second-harmonic generation, sum-frequency generation and four-wave mixing processes can manifest.

Optics↗
DOE OSTI2024

Solving the one-dimensional penetration problem for the fission channel in the statistical Hauser-Feshbach theory

Here we solve the Schrödinger equation for an arbitrary one-dimensional potential energy to calculate the transmission coefficient in the fission channel of compound nucleus reactions. We incorporate the calculated transmission coefficients into the statistical Hauser-Feshbach model calculation for neutron-induced reactions on 235,238 U and 239 Pu. The one-dimensional model reproduces the evaluated fission cross section data reasonably well considering the limited number of model parameters involved. A resonance-like structure appears in the transmission coefficient for a double-humped fission barrier shape that includes an intermediate well, which is understood to be a quantum mechanical effect in the fission channel. The calculated fission cross sections for the neutron-induced reactions on 235,238 U and 239 Pu all exhibit a similar structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A probabilistic-phase field model for the fracture of brittle materials

Abstract We develop a computational method to determine the failure probability of brittle materials under general mechanical loading conditions. The method is a combination of two parts: (1) numerical simulations of materials with multiple cracks using phase field theory, where the complete fracture process is viewed as ‘damage percolation’ along critical paths or clusters of cracks, rather than the traditional weak-link failure mechanism of Weibull, and (2) an extension of the Batdorf statistical theory of fracture to finite domains, where it is implemented within the finite element framework. The results of phase-field simulations at the ‘percolation threshold’ are used as failure data in the Batdorf theory to determine the overall probability of failure. The input to this approach is the size distribution of cracks in a pristine material. An example is shown, where alumina samples that were previously tested by Abe and coworkers (Abeet al2003J. Am. Ceram. Soc.861019–21) in four-point loading are compared to the results of our numerical simulations. The approach developed here has the advantage of being extendable to more complex thermomechanical loading.

Materials Science↗
DOE OSTI2024

Mechanisms of Low-Level Jet Formation in the U.S. Mid-Atlantic Offshore

Abstract Low-level jets (LLJs), in which the wind speed attains a local maximum at low altitudes, have been found to occur in the U.S. mid-Atlantic offshore, a region of active wind energy deployment as of 2023. In contrast to widely studied regions such as the U.S. southern Great Plains and the California coastline, the mechanisms that underlie LLJs in the U.S. mid-Atlantic are poorly understood. This work analyzes floating lidar data from buoys deployed in the New York Bight to understand the characteristics and causes of coastal LLJs in the region. Application of the Hilbert–Huang transform, a frequency analysis technique, to LLJ case studies reveals that mid-Atlantic jets frequently occur during times of adjustment in synoptic-scale motions, such as large-scale temperature and pressure gradients or frontal passages, and that they do not coincide with motions at the native inertial oscillation frequency. Subsequent analysis with theoretical models of inertial oscillation and thermal winds further reveals that these jets can form in the stationary geostrophic wind profile from horizontal temperature gradients alone—in contrast to canonical LLJs, which arise from low-level inertial motions. Here, inertial oscillation can further modulate the intensity and altitude of the wind speed maximum. Statistical evidence indicates that these oscillations arise from stable stratification and the associated frictional decoupling due to warmer air flowing over a cold sea surface during the springtime land–sea breeze. These results improve our conceptual understanding of mid-Atlantic jets and may be used to better predict low-level wind speed maxima. Significance Statement The purpose of this work is to identify and characterize the atmospheric mechanisms that result in an occasional low-level maximum in the wind speed off the U.S. mid-Atlantic coastline. Our findings show that these low-level jets form due to horizontal temperature gradients arising from fronts and synoptic systems, as well as from the land–sea breeze that forces warmer air over the cold ocean surface. This work aids predictability of such jets, improves our understanding of this coastal environment, and has implications for future deployment of offshore wind energy in this region.

17 WIND ENERGY↗
DOE OSTI2023

A data driven approach for cross-slip modelling in continuum dislocation dynamics

Cross-slip is a thermally activated process by which screw dislocation changes its glide plane to another slip plane sharing the same Burgers vector. The rate at which this process happens is determined by a Boltzmann type expression that is a function of the screw segment length and the stress acting on the dislocation. In continuum dislocation dynamics (CDD), the information regarding the length of the screw dislocation segment and local stress state on dislocations are lost due to the coarse-grained representation of the density. Here, in this work, a data driven approach to characterize the lost information by analyzing the discrete dislocation configurations is proposed to enable cross-slip modeling in the CDD framework in terms of the coarse-grained dislocation density and stress fields. The analysis showed that the screw segment length follows an exponential distribution, and the stress fluctuations, defined as the difference between the stress on the dislocations and the mean field stress in CDD, follows a Lorentzian distribution. A novel approach for cross slip implementation in CDD employing the screw segment length and stress fluctuation statistics was proposed and rigorously tested by comparing the CDD cross-slip rates with discrete dislocation dynamics (DDD) rates. This approach has been applied in conjunction with three cross-slip models used in DDD simulations differing mainly in the functional form of cross slip activation energy. It was found that different cross-slip activation energy formulations yielded different cross-slip rates, yet the effect on mechanical stress-strain response and dislocation density evolution was minimal for the [001] type loading.

42 ENGINEERING↗
DOE OSTI2023

Equibiaxial flexural strength determination of UO 2 using a ball-on-ring test

To increase nuclear fuel performance and reliability, their mechanical properties require an accurate and statistically relevant assessment. Biaxial flexural strength tests provide an alternative to bend bar techniques for assessing mechanical behavior; namely, the transverse rupture strength (TRS) of ceramic samples. Biaxial test samples require simple geometries and minimal surface preparation, reducing fabrication costs and handling hazards. This study investigated the TRS of polycrystalline UO 2 fuel forms at room temperature using a ball-on-ring test fixture. Here, pellets were fabricated from UO 2 powder using conventional powder processing and sintering techniques. The TRS and Weibull parameters were obtained through Weibull statistics on over 60 UO 2 samples tested under equibiaxial flexure. The larger sample size in this study enabled a more robust Weibull statistical analysis than alternative test methods, which may not capture the stochastic failure of typical ceramics. Furthermore, two different loading ball diameters were employed to assess the impact of contact damage on fracture strength. While Hertzian contact damage was observed with the smaller loading ball, the fracture strength remained unaffected. A fracture analysis of the tested UO 2 samples indicated a mixture of intergranular and transgranular fracture that transitions to transgranular fracture with increasing distance from the fracture origin. The characteristic strength of the combined data sets was determined to be 148 MPa, and the Weibull modulus was determined to be 9.1. The TRS values and Weibull parameters were close to values found in the literature for alternative testing techniques using samples with similar microstructure and density. The findings in this study validate the ball-on-ring method used to obtain the TRS of UO 2 with a sample geometry more representative of nuclear fuels. Additionally, experimental TRS results from this study can be implemented in modeling codes to predict fuel performance, which is critical to fuel burnup extension and advanced nuclear fuel technologies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Exotic Higgs Decays

Exotic decays of the Standard Model (SM)-like Higgs boson into beyond-the-SM particles are predicted in a wide range of well-motivated theories. The enormous samples of Higgs bosons that have been and will be produced at the Large Hadron Collider thus constitute one of the key discovery opportunities at that facility, particularly in the upcoming high-statistics, high-luminosity run. Here we review recent theoretical work on models that predict or accommodate exotic Higgs decays, survey the status of current experimental searches, and look forward to future capabilities at dedicated Higgs factories and beyond.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Single-Particle Studies Reveal a Nanoscale Mechanism for Elastic, Bright, and Repeatable ZnS:Mn Mechanoluminescence in a Low-Pressure Regime

Mechanoluminescent materials, which emit light in response to elastic deformation, are demanded for use as in situ stress sensors. ZnS doped with Mn is known to exhibit one of the lowest reported thresholds for appearance of mechanoluminescence, with repeatable light emission under contact pressure <10 MPa. Here, the physical basis for such behavior remains as yet unclear. Here, reliable microscopic detection of mechanoluminescence of single ZnS:Mn microparticles, in combination with nanoscale structural characterization, provides evidence that the mechanoluminescent properties of these particles result from interplay between a non-centrosymmetric crystal lattice and its defects, viz., dislocations and stacking faults. Statistical analysis of the distributions of mechanoluminescence energy release trajectories reveals two distinct mechanisms of excitation: one attributable to a piezo-phototronic effect and the other due to dislocation motion. At pressures below 8.1 MPa, both mechanisms contribute to mechanoluminescent output, with a dominant contribution from the piezo-phototronic mechanism. In contrast, above 8.1 MPa, dislocation motion is the primary excitation source. For the piezo-phototronic mechanism, we propose a specific model that accounts for elastic ZnS:Mn mechanoluminescence under very low pressure. The charged interfaces in stacking faults lead to the presence of filled traps, which otherwise would be empty in the absence of the built-in electric field. Upon application of external stress, local enhancement of the piezoelectric field at the stacking faults’ interfaces facilitates release of the trapped carriers and subsequent luminescence. This field enhancement explains how <10 MPa pressure produces thousands of photons.

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

Calculation of self-diffusion coefficients in supercritical carbon dioxide using mean force kinetic theory

This paper presents an application of mean force kinetic theory (MFT) to the calculation of the self-diffusivity of CO 2 in the supercritical fluid regime. Two modifications to the typical application of MFT are employed to allow its application to a system of molecular species. Therefore, the first is the assumption that the inter-particle potential of mean force can be obtained from the molecule center-of-mass pair correlation function, which in the case of CO 2 is the C–C pair correlation function. The second is a new definition of the Enskog factor that describes the effect of correlations at the surface of the collision volume. The new definition retains the physical picture that this quantity represents a local density increase, resulting from particle correlations, relative to that in the zero density homogeneous fluid limit. These calculations are facilitated by the calculation of pair correlation functions from molecular dynamics (MD) simulations using the FEPM2 molecular CO2 model. The self-diffusivity calculated from theory is in good agreement with that from MD simulations up to and slightly beyond the density at the location of the Frenkel line. The calculation is compared with and is found to perform similarly well to other commonly used models but has a greater potential for application to systems of mixed species and to systems of particles with long range interatomic potentials due to electrostatic interactions.

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

Mechanical Properties of Irradiated U-10wt%Mo Alloy Degraded by Porosity Development

A plate-type nuclear fuel consisting of a solid monolithic foil of U-10wt%Mo is under development for use in the United States’ high performance research reactors. In support of developing this fuel, the fuel has been fabricated for the first time by a commercial fuel vendor and subsequently irradiated in a test reactor. This provides an opportunity to evaluate post-irradiation mechanical properties of commercially fabricated fuel. Four-point bend testing was conducted on the irradiated U-10Mo fuel and the data produced includes bending strength and Young’s modulus. Although the material behaves in a brittle manner due to the developed porosity, a general trend of strength and modulus reduction are found as fission density increases. The data produced is evaluated using both Weibull statistics and a modulus degradation model with recommendations provided.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Indirect climate impacts of the Hunga eruption

Injection of sulfur and water vapour by the Hunga volcanic eruption significantly altered chemical composition and radiative budget of the stratosphere. Yet, whether the eruption could also affect surface climate, especially via indirect pathways, remains poorly understood. Here we investigate these effects using large ensembles of simulations with the CESM2(WACCM6) Earth system model, incorporating interactive chemistry and aerosols in both coupled ocean and atmosphere-only configurations. We find some statistically significant extratropical regional climate responses to the eruption driven by circulation changes; these are partially linked to the modulation of El Niño–Southern Oscillation, and its associated teleconnections, and to perturbations of the stratospheric polar vortex in both hemispheres. The stratospheric anomalies affect surface climate through modulating the North Atlantic Oscillation in the Northern Hemisphere (up to three boreal winters following the eruption) and the Southern Annular Mode in the Southern Hemisphere in late 2023. The latter is partly related to a concurrent reduction in Antarctic ozone, as increased stratospheric aerosols and water vapor reach the polar vortex. Our study suggests that the eruption could have had a non-negligible influence on regional surface climate, and discusses the mechanisms via which such an influence could occur. However, the results also highlight that this forcing is relatively weak compared to interannual variability, and is subject to model uncertainties in the representation of key processes. More research is thus needed before definitive statements on the role of the eruption in contributing to surface climate and weather events in the following years are made.

Geosciences↗
DOE OSTI2020

A Bayesian analysis of nuclear deformation properties with Skyrme energy functionals

In spite of numerous scientific and practical applications, there is still no comprehensive theoretical description of the nuclear fission process based solely on protons, neutrons and their interactions. Furthermore, the most advanced simulations of fission are currently carried out within nuclear density functional theory (DFT). In spite of being fully quantum-mechanical and rooted in the theory of nuclear forces, DFT still depends on a dozen or so parameters characterizing the energy functional. Calibrating these parameters on experimental data results in uncertainties that must be quantified for applications. This task is very challenging because of the high computational cost of DFT calculations for fission. In this paper, we use Gaussian processes to build emulators of DFT models in order to quantify and propagate statistical uncertainties of theoretical predictions for a range of nuclear deformations relevant to describing the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Adaptive Cognitive Mechanisms to Maintain Calibrated Trust and Reliance in Automation

Trust calibration for a human–machine team is the process by which a human adjusts their expectations of the automation’s reliability and trustworthiness; adaptive support for trust calibration is needed to engender appropriate reliance on automation. Herein, we leverage an instance-based learning ACT-R cognitive model of decisions to obtain and rely on an automated assistant for visual search in a UAV interface. This cognitive model matches well with the human predictive power statistics measuring reliance decisions; we obtain from the model an internal estimate of automation reliability that mirrors human subjective ratings. The model is able to predict the effect of various potential disruptions, such as environmental changes or particular classes of adversarial intrusions on human trust in automation. Finally, we consider the use of model predictions to improve automation transparency that account for human cognitive biases in order to optimize the bidirectional interaction between human and machine through supporting trust calibration. The implications of our findings for the design of reliable and trustworthy automation are discussed.

60 APPLIED LIFE SCIENCES↗