Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational 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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Integrated Quantum-Classical Protocol for the Realistic Description of Solvated Multinuclear Mixed-Valence Transition-Metal Complexes and Their Solvatochromic Properties

Linear cyanide-bridged polymetallic complexes, which undergo photoinduced metal-to-metal charge transfer, represent prototypical systems for studying long-range electron-transfer reactions and understanding the role played by specific solute–solvent interactions in modulating the excited-state dynamics. Here, to tackle this problem, while achieving a statistically meaningful description of the solvent and of its relaxation, one needs a computational approach capable of handling large polynuclear transition-metal complexes, both in their ground and excited states, as well as the ability to follow their dynamics in several environments up to nanosecond time scales. Here, we present a mixed quantum classical approach, which combines large-scale molecular dynamics (MD) simulations based on an accurate quantum mechanically derived force field (QMD-FF) and self-consistent QMD polarized point charges, with IR and UV–vis spectral calculations to model the solvation dynamics and optical properties of a cyano-bridged trinuclear mixed-valence compound (trans-[(NC) 5 Fe III (μ-CN)Ru II (pyridine) 4 (μ-NC)Fe III (CN) 5 ] 4– ). We demonstrate the reliability of the QMD-FF/MD approach in sampling the solute conformational space and capturing the local solute–solvent interactions by comparing the results with higher-level quantum mechanics/molecular mechanics (QM/MM) MD reference data. The IR spectra calculated along the classical MD trajectories in different solvents correctly predict the red shift of the CN stretching band in the aprotic medium (acetonitrile) and the subtle differences measured in water and methanol, respectively. By explicitly including the solvent molecules around the cyanide ligands and calculating the thermal averaged absorption spectra using time-dependent density functional theory calculations within the Tamm–Dancoff approximation, the experimental solvatochromic shift is quantitatively reproduced going from water to methanol, while it is overestimated for acetonitrile. This discrepancy can likely be traced back to the lack of important dispersion interactions between the solvent cyano groups and the pyridine substituents in our micro solvation model. The proposed protocol is applied to the ground state in water, methanol, and acetonitrile and can be flexibly generalized to study excited-state nonequilibrium solvation dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

GRCop-84 Developed for Rocket Engines

GRCop-84 (Cu-8 at.% Cr-4 at.% Nb) was developed at the NASA Glenn Research Center for use in regeneratively cooled rocket engines. This alloy possesses an excellent combination of conductivity, thermal expansion, strength, creep resistance, ductility, and low-cycle fatigue (LCF) life. In comparison to the currently used alloy, NARloy-Z (Cu-3 wt % Ag-0.5 wt % Zr), GRCop-84's significantly better properties give it the potential for significant gains in engine performance and reliability. For GRCop-84 to be safely used in a rocket engine, it is critical to develop a detailed data base of the important thermophysical and mechanical properties. Work has focused on five major properties: thermal conductivity, thermal expansion, tensile strength, creep resistance, and LCF life. The analysis went beyond simply measuring the data and reporting averages. A detailed statistical analysis was conducted that allowed regression of the data over the entire temperature range tested and establishment of design minimums. The design values were expressed as simple mathematical formulas that are highly amenable to use in computer codes such as in finite element analysis and related computer modeling work.

Ellis, David L.↗
DOE OSTI2024

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↗
DOE OSTI2020

A practical and efficient approach for Bayesian quantum state estimation

Bayesian inference is a powerful paradigm for quantum state tomography, treating uncertainty in meaningful and informative ways. Yet the numerical challenges associated with sampling from complex probability distributions hampers Bayesian tomography in practical settings. In this article, we introduce an improved, self-contained approach for Bayesian quantum state estimation. Leveraging advances in machine learning and statistics, our formulation relies on highly efficient preconditioned Crank–Nicolson sampling and a pseudo-likelihood. We theoretically analyze the computational cost, and provide explicit examples of inference for both actual and simulated datasets, illustrating improved performance with respect to existing approaches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1961

THEORETICAL STUDY OF ATOMIC AND MOLECULAR GASES AND THEIR REACTIONS IN THE UPPER ATMOSPHERE. SUMMARY OF WORK DONE IN THE PERIOD NOV. 1, 1960 TO APR. 30, 1961

In the previous report it was stated that the major mathematical problems of computations were nearly solved and an attempt is being made to start application of "Tables of Integrals" to the computation of the different states of atoms and molecules. In this respect I like to report that the computation of a number of states of atoms and molecules has been carried out. The results are so far confined only to the equilibrium distance of molecules. We are extending the calculation to other distances so as to compute potential energy curves of different states of molecules. We are also developing Thomas Fermi Statistical model of the atom to verify our results of more accurate calculation. The results so far have been very encouraging and I hope to report some of these results in the next paper.

Gas↗
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 OSTI2025

Universal energy-speed-accuracy trade-offs in driven nonequilibrium systems

The connection between measure theoretic optimal transport and dissipative nonequilibrium dynamics provides a language for quantifying nonequilibrium control costs, leading to a collection of thermodynamic speed limits, which rely on the assumption that the target probability distribution is perfectly realized. This is almost never the case in experiments or numerical simulations, so here we address the situation in which the external controller is imperfect. We obtain a lower bound for the dissipated work in generic nonequilibrium control problems that (1) is asymptotically tight and (2) matches the thermodynamic speed limit in the case of optimal driving. Along with analytically solvable examples, we refine this imperfect driving notion to systems in which the controlled degrees of freedom are slow relative to the nonequilibrium relaxation rate, and identify independent energy contributions from fast and slow degrees of freedom. Furthermore, we develop a strategy for optimizing minimally dissipative protocols based on optimal transport flow matching, a generative machine learning technique. Furthermore, this latter approach ensures the scalability of both the theoretical and computational framework we put forth. Crucially, we demonstrate that we can compute the terms in our bound numerically using efficient algorithms from the computational optimal transport literature and that the protocols we learn saturate the bound.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Quantum Simulations of Radiation Damage in a Molecular Polyethylene Analog

Abstract An atomic‐level understanding of radiation‐induced damage in simple polymers like polyethylene is essential for determining how these chemical changes can alter the physical and mechanical properties of important technological materials such as plastics. Ensembles of quantum simulations of radiation damage in a polyethylene analog are performed using the Density Functional Tight Binding method to help bind its radiolysis and subsequent degradation as a function of radiation dose. Chemical degradation products are categorized with a graph theory approach, and occurrence rates of unsaturated carbon bond formation, crosslinking, cycle formation, chain scission reactions, and out‐gassing products are computed. Statistical correlations between product pairs show significant correlations between chain scission reactions, unsaturated carbon bond formation, and out‐gassing products, though these correlations decrease with increasing atom recoil energy. The results present relatively simple chemical descriptors as possible indications of network rearrangements in the middle range of excitation energies. Ultimately, the work provides a computational framework for determining the coupling between nonequilibrium chemistry in polymers and potential changes to macro‐scale properties that can aid in the interpretation of future radiation damage experiments on plastic materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Thermalization in the one-dimensional Salerno model lattice

The Salerno model constitutes an intriguing interpolation between the integrable Ablowitz-Ladik (AL) model and the more standard (nonintegrable) discrete nonlinear Schrödinger (DNLS) one. The competition of local on-site nonlinearity and nonlinear dispersion governs the thermalization of this model. Here, we investigate the statistical mechanics of the Salerno one-dimensional lattice model in the nonintegrable case and illustrate the thermalization in the Gibbs regime. As the parameter interpolating between the two limits (from DNLS toward AL) is varied, the region in the space of initial energy and norm densities leading to thermalization expands. The thermalization in the non-Gibbs regime heavily depends on the finite system size; we explore this feature via direct numerical computations for different parametric regimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Deterministic Calibration of MFiX-PIC, Part 1: Settling Bed

The Particle-in-cell (PIC) numerical approach for modeling granular solids in fluid flow has gained significant interest in recent years. Valued for its often shorter time-to-solution, the PIC formulation relies on modeling statistical groupings of particles called parcels in cooperation with a solids stress model to affect local solids velocity. This is in contrast to the discrete element model (DEM) where every particle in a system is modelled individually and directly coupled to local solids velocity through Newtonian mechanics. The U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) develops and maintains Multiphase Flow with Interphase eXchanges (MFiX), a collection of open-source computational fluid dynamics (CFD) solvers. Included in the MFiX suite are traditional two-fluid model (TFM) and DEM solvers, and a recently added PIC solver (NETL, 2021). In general, PIC methodologies offer an accuracy trade-off in lieu of computational speed; and therefore, it is important to assess the credibility of MFiX-PIC simulations. For this purpose, a systematic verification, validation and uncertainty quantification (VVUQ) effort was initiated at NETL to assess the new PIC solver

42 ENGINEERING↗
DOE OSTI2020

Photoproduction of Cascade Baryons Using the GlueX Detector at Jefferson Laboratory

Multi-strange baryons have played an important role in understanding the strong interaction. Despite their importance to advancing quantum chromodynamics (QCD), little is known about such hyperons due to the lack of a strange probe. Almost all knowledge of Cascade baryons today stems from kaon-nucleon interactions in bubble chamber experiments performed in the 1960s and 1970s. The 1980s saw advancement in electronics and computational resources allowing for higher statistic experiments such as the MultiParticle Spectrometer (MPS) at Brookhaven National Laboratory. From these previous experiments, the octet and decuplet ground states, Xi(1320) and Xi(1530) respectively, were well established with known properties and four excited states were well established. This research uses the GlueX experiment at Jefferson Laboratory to extend photoproduction cross section measurements for Xi(1320)- to higher beam energies to further understand the production mechanism, as well as to search for the first excited Cascade resonance observation in photoproduction through the decay Xi-* -> K-Lambda. The GlueX Spectrometer is a nearly hermetic detector used to study exclusive photoproduction reactions with photon beam energies up to 12 GeV. The first phase of GlueX running has recently been completed. Preliminary results for the total cross section and differential cross sections in terms of momentum transfer -t are shown within. Additionally, the first observation of an excited Cascade baryon, in particular Xi(1820)-, is shown.

Ernst, Ashley↗
NASA NTRS1985

Theoretical Studies of the Surface Tension of Liquid Metal System

A major goal of this project is to understand the surface tension and other thermophysical properties of liquid metals and alloys from a fundamental viewpoint. The approach is to calculate these quantities by a first principles technique which combines the statistical-mechanical theory of the liquid state with an electronic pseudopotential theory of electrons in metals. The inhomogeneity of the surface is treated using an ionic-density-functional formalism developed with the support of NASA. Of particular interest are the variation of surface tension with temperature and impurity concentration: such variations strongly influence the types of convection which make take place in a low-gravity environment. Some progress has already been achieved in computing the reduction of surface tension due to the presence of low-surface-tension impurities, and the corresponding surface segregation of such impurities. In the coming year, it is planned to concentrate on the surface properties of materials of particular interest to the MSA program: Si, Ga and GaSn alloys. An additional goal is to gain some theoretical understanding of the high temperature thermophysical properties of liquid metals, particularly high melting point materials which have not been studied extensively from a theoretical viewpoint.

Stroud, D. G.↗
NASA NTRSDate not supplied

Correlations Between Porosity, Spatter, and Process Metrics for Powder Bed Fusion Laser Beam Metallic Additive Manufacturing

Components fabricated using the powder bed fusion laser beam metallic (PBF-LB/M)additive manufacturing process are the result of a multitude of weld passes conducted sequentially. Qualifying components for aerospace applications requires a thorough understanding of the process-structure-properties relationships. Porosity defects are known to have a strong adverse effect on the mechanical properties of a component. In particular, porosity defects created by lack of fusion have high aspect ratio morphologies leading to stress concentrations that become crack initiation sites. In the present work, the occurrence of spatter induced lack of fusion porosity was studied using synchronized in-situ process monitoring, additive manufacturing model-based process metrics, and high-resolution X-ray computed tomography. The results show that lack of fusion porosity is statistically correlated with unremoved welding spatter ejecta of the PBF-LB/M process and process metrics related to the hatching strategy.

Qualification↗
DOE OSTI2022

Software Development for Data Visualization and Analysis of PN-Diodes & MOSFET Devices

The U.S. Department of Energy has a broad mission to ensure the nation’s security by addressing ongoing environmental challenges, developing novel energy production technologies, and mitigating nuclear security concerns. These efforts benefit from developing more effective semiconductor materials and devices. The research discussed in this report aims to further the development of these new electronics by improving implementations of vertical device structures, which offer theoretically better performance compared to lateral structures. Specifically, it reviews the background and development of a new and easy-to-use program for data processing and analysis from experiments run on vertical gallium nitride power devices. Current vertical structures fall short in critical performance metrics, such as On-Resistance and Breakdown Voltage, due to poor management of the device’s electric field. Junction termination extensions, or JTEs, are a field management technique that increases a device’s resilience to expected failure modes by controlling its surface electric fields. Each JTE design requires significant experimental validation, consisting of hundreds of tested devices, each outputting an enormous set of data points. The new software includes multiple methods of filtering large sets of device testing data to identify and remove flawed devices, automate device analyses, and provide statistical breakdowns for four metrics used to define the effectiveness of the JTE and the current carrying capacity of the device. These results can then be related back to fabrication techniques and device design, and in the future, be combined with processing simulations to further improve our understanding of the JTE mechanisms.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Polarimetric microwave scattering from known randomly rough surfaces

Scattering from rough surfaces is reviewed, and artificial randomly rough surfaces are described. Measurements were made from 4 to 18 GHz for the 4 linear polarization states. The target consists of a metalized randomly rough surface constructed from known surface statistical information. The surface was generated using surface generation routines on the computer, and then the data was converted for the machining of the surface on a computer controlled mill. Bistatic and monostatic measurements were acquired at incident angles from 25 to 60 deg. At the 13.5 to 18 GHz frequencies, the effects of backscatter enhancement and depolarization effects can be studied in a regime which does not satisfy Kirchoff or small perturbation assumptions. This allows scattering mechanisms to be introduced into models, which may show missing elements in the more standard formulations. At lower frequency bands such as 4 to 8.5 GHz the scattering behavior was also examined to see the behavior as assumptions of other formulations are satisfied.

Nance, C. E.↗
DOE OSTI2021

Guiding the Design of Heterogeneous Electrode Microstructures for Li-Ion Batteries: Microscopic Imaging, Predictive Modeling, and Machine Learning

Electrochemical and mechanical properties of lithium-ion battery materials are heavily dependent on their 3D microstructure characteristics. A quantitative understanding of the role played by stochastic microstructures is critical for the prediction of material properties and for guiding synthesis processes. Furthermore, tailoring microstructure morphology is also a viable way of achieving optimal electrochemical and mechanical performances of lithium-ion cells. To facilitate the establishment of microstructure-resolved modeling and design methods, a review covering spatially and temporally resolved imaging of microstructure and electrochemical phenomena, microstructure statistical characterization and stochastic reconstruction, microstructure-resolved modeling for property prediction, and machine learning for microstructure design is presented here. The perspectives on the unresolved challenges and opportunities in applying experimental data, modeling, and machine learning to improve the understanding of materials and identify paths toward enhanced performance of lithium-ion cells are presented.

25 ENERGY STORAGE↗
NASA NTRS1987

Advanced probabilistic method of development

Advanced structural reliability methods are utilized on the Probabilistic Structural Analysis Methods (PSAM) project to provide a tool for analysis and design of space propulsion system hardware. The role of the effort at the University of Arizona is to provide reliability technology support to this project. PSAM computer programs will provide a design tool for analyzing uncertainty associated with thermal and mechanical loading, material behavior, geometry, and the analysis methods used. Specifically, reliability methods are employed to perform sensitivity analyses, to establish the distribution of a critical response variable (e.g., stress, deflection), to perform reliability assessment, and ultimately to produce a design which will minimize cost and/or weight. Uncertainties in the design factors of space propulsion hardware are described by probability models constructed using statistical analysis of data. Statistical methods are employed to produce a probability model, i.e., a statistical synthesis or summary of each design variable in a format suitable for reliability analysis and ultimately, design decisions.

Wirsching, P. H.↗
NASA NTRS2009

A Biomechanical Approach to Assessing Hip Fracture Risk

Bone loss in microgravity is well documented, but it is difficult to quantify how declines in bone mineral density (BMD) contribute to an astronaut's overall risk of fracture upon return. This study uses a biomechanical approach to assessing hip fracture risk, or Factor of Risk (Phi), which is defined as the ratio of applied load to bone strength. All long-duration NASA astronauts from Expeditions 1-18 were included in this study (n=25), while crewmembers who flew twice (n=2) were treated as separate subjects. Bone strength was estimated based on an empirical relationship between areal BMD at the hip, as measured by DXA, and failure load, as determined by mechanical testing of cadaver femora. Fall load during a sideways fall was calculated from a previously developed biomechanical model, which takes into account body weight, height, gender, and soft tissue thickness overlying the lateral aspect of the hip that serves to attenuate the impact force. While no statistical analyses have been performed yet, preliminary results show that males in this population have a higher FOR than females, with a post- flight Phi of 0.87 and 0.36, respectively. FOR increases 5.1% from preflight to postflight, while only one subject crossed the fracture "threshold" of Phi = 1, for a total of 2 subjects with a postflight Phi > 1. These results suggest that men may be at greater risk for hip fracture due largely in part to their relatively thin soft tissue padding as compared to women, since soft tissue thickness has the highest correlation (R(exp 2)= .53) with FOR of all subject-specific parameters. Future work will investigate changes in FOR during recovery to see if baseline risk levels are restored upon return to 1-g activity. While dual x-ray absorptiometry (DXA) is the most commonly used clinical measure of bone health, it fails to provide compartment-specific information that is useful in assessing changes to bone quality as a result of microgravity exposure. Peripheral quantitative computed tomography (pQCT) accomplishes this by imaging transverse "slices" of the long bones. This project was a re-analysis of a 90 day bed rest study to determine if changes to cortical and trabecular compartments could be detected in the distal tibia with statistical significance using a new pQCT image analysis method. Nearly all changes in bone mineral density (BMD) and cross sectional area (CSA) measures were seen with statistical significance, with the exception of a change in cortical BMD. Total bone CSA increased by 1.1 % (p =0.01), cortical CSA decreased by - 5.6% (p<0.001) and trabecular CSA increased by 1.76% (p=0.007); the combination of which suggests bone resorption occurred at the endocortical surface in response to mechanical unloading by bed rest. Furthermore, total BMD and trabecular BMD decreased (-3.8%, p=0.001 and -2.8%, p =0.007, respectively), while decreases in cortical BMD failed to reach significance (-1.2%, p=0.07). Given that compartment-specific changes are seen with significance and are likely to influence bone strength, it is recommended that pQCT remain a standard measure used in bed rest because it provides a unique measure by which to better evaluate the efficacy of countermeasures to microgravity-induced bone loss.

Ellman, Rachel↗