Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials modelling”

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 1,369 records · Page 76

Discrete-Direct Model Calibration and Uncertainty Propagation Method Confirmed on Multi-Parameter Plasticity Model Calibrated to Sparse Random Field Data

A discrete direct (DD) model calibration and uncertainty propagation approach is explained and demonstrated on a 4-parameter Johnson-Cook (J-C) strain-rate dependent material strength model for an aluminum alloy. The methodology's performance is characterized in many trials involving four random realizations of strain-rate dependent material-test data curves per trial, drawn from a large synthetic population. The J-C model is calibrated to particular combinations of the data curves to obtain calibration parameter sets which are then propagated to “Can Crush” structural model predictions to produce samples of predicted response variability. These are processed with appropriate sparse-sample uncertainty quantification (UQ) methods to estimate various statistics of response with an appropriate level of conservatism. This is tested on 16 output quantities (von Mises stresses and equivalent plastic strains) and it is shown that important statistics of the true variabilities of the 16 quantities are bounded with a high success rate that is reasonably predictable and controllable. The DD approach has several advantages over other calibration-UQ approaches like Bayesian inference for capturing and utilizing the information obtained from typically small numbers of replicate experiments in model calibration situations—especially when sparse replicate functional data are involved like force–displacement curves from material tests. The DD methodology is straightforward and efficient for calibration and propagation problems involving aleatory and epistemic uncertainties in calibration experiments, models, and procedures.

42 ENGINEERING↗

Modeling Single-Crystal Battery Materials: From Fundamental Understanding to Performance Evaluation

The performance of rechargeable batteries is fundamentally influenced by the physicochemical properties and microstructural features of their key material components. Recent experimental advancements have highlighted the potential of single-crystal (SC) morphologies to address inherent limitations of polycrystalline (PC) electrodes and solid-state electrolytes, offering tunable charge transport kinetics and improved cell cycling performance. Here, this review examines how state-of-the-art computational modeling, from atomistic and mesoscale to continuum-level approaches, including machine learning methodologies, has been utilized to investigate the critical factors governing the electrochemical behavior of SC battery materials. We explore how predictive modeling can elucidate the processing–structure–property–performance relationships of SC cathodes, anodes, and solid-state electrolytes, with a focus on unique SC characteristics such as crystallographic anisotropy, size effects, and facet-dependent properties. Additionally, we identify limitations in commonly used modeling techniques and discuss strategies to address these challenges. By integrating high-fidelity simulations with experimental insights, this review aims to outline a clear path for the rational design and optimization of SC battery components, paving the way for accelerated advancements in energy storage technologies.

Materials science↗

An Eulerian crystal plasticity framework for modeling large anisotropic deformations in energetic materials under shocks

Here, this paper demonstrates a novel Eulerian computational framework for modeling anisotropic elastoplastic deformations of organic crystalline energetic materials (EM) under shocks. While Eulerian formulations are advantageous for handling large deformations, constitutive laws in such formulations have been limited to isotropic elastoplastic models, which may not fully capture the shock response of crystalline EM. The present Eulerian framework for high-strain rates, large deformation material dynamics of EM incorporates anisotropic isochoric elasticity via a hypo-elastic constitutive law and visco-plastic single-crystal models. The calculations are validated against atomistic calculations and experimental data and benchmarked against Lagrangian (finite element) crystal plasticity computations for shock-propagation in a monoclinic organic crystal, octahydro-1,3,5,7-tetranitro-1,3,5,7 tetrazocine (β-HMX). The Cauchy stress components and the resolved shear stresses calculated using the present Eulerian approach are shown to be in good agreement with the Lagrangian computations for different crystal orientations. The Eulerian framework is then used for computations of shock-induced inert void collapse in β-HMX to study the effects of crystal orientations on hotspot formation under different loading intensities. The computations show that the hotspot temperature distributions and the collapse profiles are sensitive to the crystal orientations at lower impact velocities (viz., 500 m/s); when the impact velocity is increased to 1000 m/s, the collapse is predominantly hydrodynamic and the role of anisotropy is modest. The present methodology will be useful to simulate energy localization in shocked porous energetic material microstructures and other situations where large deformations of single and polycrystals govern the thermomechanical response.

42 ENGINEERING↗

Engineering and fabrication cost considerations for cryogenic wind tunnel models

Design and fabrication cost drivers for cryogenic transonic wind tunnel models are defined. The major cost factors for wind tunnel models are model complexity, tolerances, surface finishes, materials, material validation, and model inspection. The cryogenic temperatures require the use of materials with relatively high fracture toughness but at the same time high strength. Some of these materials are very difficult to machine, requiring extensive machine hours which can add significantly to the manufacturing costs. Some additional engineering costs are incurred to certify the materials through mechanical tests and nondestructive evaluation techniques, which are not normally required with conventional models. When instrumentation such as accelerometers and electronically scanned pressure modules is required, temperature control of these devices needs to be incorporated into the design, which requires added effort. Additional thermal analyses and subsystem tests may be necessary, which also adds to the design costs. The largest driver to the design costs is potentially the additional static and dynamic analyses required to insure structural integrity of the model and support system.

Boykin, R. M., Jr.↗

First-principles investigation of the resistive switching energetics in monolayer MoS 2 : insights into metal diffusion and adsorption

A deeper understanding of resistive switching (RS) in 2D materials is essential for advancing neuromorphic computing. The Dissociation-Diffusion-Adsorption (DDA) model offers a useful framework for probing RS mechanisms in non-volatile memory (NVM) and in-memory computing. We have employed first-principles density functional theory (DFT) to explore dissociation, diffusion, and adsorption phenomena within the DDA model, focusing on the interactions between exemplary metal atoms (Au, Ag, Cu) and monolayer MoS 2 . Nudged elastic band (NEB) calculations evaluated diffusion barriers in pristine and sulfur-vacancy MoS 2 . Charged systems were modeled to assess the impact of applied bias on migration pathways. We also examined metal dissociation from bulk electrodes and adsorption at S vacancies. Ag/MoS 2 shows the lowest dissociation barrier (~0.034 eV), while Au and Cu exhibit similar values (~0.32 eV). These insights highlight Ag as a promising candidate for low-energy RS applications and provide guidance for optimizing switching efficiency in 2D memory devices.

Atomistic models↗

Graphite Material Behavior: NRC Graphite Behavior Model

Introduction to graphite Fabrication of graphite Unique features of nuclear grade graphite Providing the properties As-fabricated material properties of interest Thermal diff, thermal stability, mechanical strength, etc. General graphite behavior and degradation Microstructure: key to behavior Crystallographic structure Basal plane (covalent) bonding Porosity and pore microstructure Degradation Irradiation behavior Oxidation behavior Molten salt issues

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Influence of Catalysis and Oxidation on Slug Calorimeter Measurements in Arc Jets

Arc jet tests play a critical role in the characterization and certification of thermal protection materials and systems (TPS). The results from these arc jet tests feed directly into computational models of material response and aerothermodynamics to predict the performance of the TPS in flight. Thus the precise knowledge of the plasma environment to which the test material is subjected, is invaluable. As one of the environmental parameters, the heat flux is commonly measured. The measured heat flux is used to determine the plasma enthalpy through analytical or computational models. At NASA Ames Research Center (ARC), slug calorimeters of a geometrically similar body to the test article are routinely used to determine the heat flux. A slug calorimeter is a thermal capacitance-type calorimeter that uses the temperature rise in a thermally insulated slug to determine the heat transfer rate, see Figure 1(left). Current best practices for measuring the heat flux with a slug calorimeter are described in ASTM E457 - 96. Both the calorimeter body and slug are made of Oxygen Free High Conductivity Copper, and are cleaned before each run.

Nawaz, Anuscheh↗

A phase-field model of quasi-brittle fracture for pressurized cracks: Application to UO 2 high-burnup microstructure fragmentation

Here, in this paper, we present a phase-field model of quasi-brittle fracture with pressurized cracks, with dedicated applications for polycrystalline materials. The model is formulated as a minimization problem within the variational framework. The external work done by pressure on the crack surfaces is included in the objective function. Several careful modeling choices lead to a regularization-length-independent critical strength. The model is constructed to give a softening response with an underlying linear traction-separation law. The pressure-dependent softening response and the regularization of the prescribed pressure are demonstrated with a (quasi) one-dimensional numerical analysis. In a two-dimensional numerical analysis under plane strain assumptions, the critical stress corresponding to crack propagation (as predicted by our quasi-brittle fracture model) is compared with linear elastic fracture mechanics (LEFM) analytical solutions. Our model is further utilized to simulate fission-gas-induced fragmentation of the UO 2 high-burnup structure (HBS). Simulation results show that pressurized bubbles can cause crack nucleation and propagation, and that the bubble size and the surrounding external pressure affect the critical pressure corresponding to crack nucleation. Simulations of a partial HBS at different recrystallization stages show that different grain structures (due to recrystallization) also influence crack paths and fragmentation morphology.

36 MATERIALS SCIENCE↗

A model for the distribution of material generating the soft X-ray background

The observational evidence relating to the soft X-ray diffuse background is discussed, and a simple model for its source and spatial structure is presented. In this simple model with one free parameter, the observed 1/4 keV X-ray intensity originates as thermal emission from a uniform hot plasma filling a cavity in the neutral material of the Galactic disk which contains the sun. Variations in the observed X-ray intensity are due to variations in the extent of the emission volume and therefore the emission measure of the plasma. The model reproduces the observed negative correlation between X-ray intensity and H I column density and predicts reasonable values for interstellar medium parameters.

Snowden, S. L.↗

Physical processes in directed ion beam sputtering

The general operation of a discharge chamber for the production of ions is described. A model is presented for the magnetic containment of both primary and secondary or Maxwellian electrons in the discharge plasma. Cross sections were calculated for energy and momentum transfer in binary collisions between like pairs of Ar, Kr, and Xe atoms in the energy range from about 1 eV to 1000 eV. These calculations were made from available pair interaction potentials using a classical model. Experimental data from the literature were fit to a theoretical expression for the Ar resonance charge exchange cross section over the same energy range. A model was developed that describes the processes of conical texturing of a surface due to simultaneous directed ion beam etching and sputter deposition of an impurity material. This model accurately predicts both a minimum temperature for texturing to take place and the variation of cone density with temperature. It also provides the correct order of magnitude of cone separation. It was predicted from the model, and subsequently verified experimentally, that a high sputter yield material could serve as a seed for coning of a lower sputter yield substrate. Seeding geometries and seed deposition rates were studied to obtain an important input to the theoretical texturing model.

Robinson, R. S.↗

Discussion on molecular dynamics (MD) simulations of the asphalt materials

The application of asphalt materials in pavement engineering has been increasingly widespread and sophisticated over the past several decades. Variations in the properties of asphalt binder during mixing, transportation, and paving can affect the performance of asphalt pavement. However, the asphalt material is a non-homogeneous and complex organic substance, consisting of various molecules with widely various molecular weights, elemental compositions, and structures. This complexity leads to difficulties for researchers to clearly and immediately understand the properties of asphalt materials and their variations. The multi-scale research approach combines macroscopic experimental data and microscopic simulation results from a practical engineering perspective. It helps to improve the understanding of asphalt materials. The molecular dynamics (MD) simulation proposes a corresponding molecular model of asphalt material based on experimental data, and the simulation algorithm is able to derive properties similar to those of real asphalt. Here, this paper provides a comprehensive review of the current studies on MD simulation of asphalt materials, including modeling, properties, and multi-scale analysis. As a key part of the computational simulation, this paper discusses the typical asphalt binder and asphalt-aggregate interface models constructed by different groups, and also presents their differences from real samples and their feasibility based on fundamental properties. After the introduction of molecular models, the extensive work made by researchers based on molecular models is categorically reviewed and discussed. The strengths and weaknesses of MD simulation methods in the study of asphalt materials are also summarized in order to provide the reader with a more comprehensive understanding of the relevant contents and to guide subsequent research.

36 MATERIALS SCIENCE↗

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗

Computational Discovery of Intermolecular Singlet Fission Materials Using Many-Body Perturbation Theory

Intermolecular singlet fission (SF) is the conversion of a photogenerated singlet exciton into two triplet excitons residing on different molecules. SF has the potential to enhance the conversion efficiency of solar cells by harvesting two charge carriers from one high-energy photon, whose surplus energy would otherwise be lost to heat. The development of commercial SF-augmented modules is hindered by the limited selection of molecular crystals that exhibit intermolecular SF in the solid state. Computational exploration may accelerate the discovery of new SF materials. The GW approximation and Bethe–Salpeter equation (GW+BSE) within the framework of many-body perturbation theory is the current state-of-the-art method for calculating the excited-state properties of molecular crystals with periodic boundary conditions. In this Review, we discuss the usage of GW+BSE to assess candidate SF materials as well as its combination with low-cost physical or machine learned models in materials discovery workflows. We demonstrate three successful strategies for the discovery of new SF materials: (i) functionalization of known materials to tune their properties, (ii) finding potential polymorphs with improved crystal packing, and (iii) exploring new classes of materials. In addition, three new candidate SF materials are proposed here, which have not been published previously.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experiments at TA-55 have used optical pyrometry to measure temperatures of shocked plutonium on the 40mm gun

Researchers at LANL, along with collaborators from MSTS, have worked for more than a decade to develop and implement an optical pyrometry method to measure the temperature of metals shocked to the high pressures and temperatures. Because traditional shock wave measurements only provide information on the mechanical state of the material, temperature measurements are required to validate and improve equation-of-state (EOS) models for materials. The conditions sampled are complex; commenserate with those found in planetary impacts, and accressed by conventional and nuclear weapons. Optical pyrometry is currently the best available method for inferring the temperature of shocked metals, by careful measurement of the calibrated radiance emitted by a surface at finite temperature. These measurements are difficult often complicated by the short time scales associated with shock wave experiments (<1 µs) coupled with the many sources of non-thermal light (impact flash, fracture light, etc.) that pollutes the radiant light measured from the sample.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Validation and Sensitivity Analyses of Arc-Jet Performance of Woven Thermal Protection Entry Systems

A material response model is developed for the determination of erosion rates and thermal response of woven materials to enable accurate thermal protection material sizing and future missions to Venus and giant planets. The stagnation-point thermal response is further evaluated using a Monte Carlo global sensitivity analysis to determine the interactions between key parameters in the material surface energy balance, such as, convective blowing correction parameter and weave material properties including isotropic and orthotropic thermal conductivites. The time-accurate material response for the dual layer weave erosion and surface temperatures show excellent agreement with stagnation flat face arc-jet tests involving time-varying aeroheating-cooling cycles. The high-fidelity Monte Carlo sensitivity analysis technique was used to address challenges in thermal protection sizing for the Adaptive Deployable Entry and Placement Technology (ADEPT) system, a low ballistic coefficient hypersonic decelerator. The results represent new correlations for such weaves involving Carbon-air oxidation equilibrium chemistry. Through the multi-variate regression analyses and uncertainty rankings performed, only three properties were found to contribute to the transient thermal response of the gore, with emissivity contributing to nearly 95% of the total output uncertainty.

Pratibha Raghunandan↗

Validation and Sensitivity Analyses of Arc-Jet Performance of Woven Thermal Protection Entry Systems

A material response model is developed for the determination of erosion rates and thermal response of woven materials to enable accurate thermal protection material sizing and future missions to Venus and giant planets. The stagnation-point thermal response is further evaluated using a Monte Carlo global sensitivity analysis to determine the interactions between key parameters in the material surface energy balance, such as, convective blowing correction parameter and weave material properties including isotropic and orthotropic thermal conductivites. The time-accurate material response for the dual layer weave erosion and surface temperatures show excellent agreement with stagnation flat face arc-jet tests involving time-varying aeroheating-cooling cycles. The high-fidelity Monte Carlo sensitivity analysis technique was used to address challenges in thermal protection sizing for the Adaptive Deployable Entry and Placement Technology (ADEPT) system, a low ballistic coefficient hypersonic decelerator. The results represent new correlations for such weaves involving Carbon-air oxidation equilibrium chemistry. Through the multi-variate regression analyses and uncertainty rankings performed, only three properties were found to contribute to the transient thermal response of the gore, with emissivity contributing to nearly 95% of the total output uncertainty.

Pratibha Raghunandan↗

Reliability Quantification of Advanced Stirling Convertor (ASC) Components

The Advanced Stirling Convertor, is intended to provide power for an unmanned planetary spacecraft and has an operational life requirement of 17 years. Over this 17 year mission, the ASC must provide power with desired performance and efficiency and require no corrective maintenance. Reliability demonstration testing for the ASC was found to be very limited due to schedule and resource constraints. Reliability demonstration must involve the application of analysis, system and component level testing, and simulation models, taken collectively. Therefore, computer simulation with limited test data verification is a viable approach to assess the reliability of ASC components. This approach is based on physics-of-failure mechanisms and involves the relationship among the design variables based on physics, mechanics, material behavior models, interaction of different components and their respective disciplines such as structures, materials, fluid, thermal, mechanical, electrical, etc. In addition, these models are based on the available test data, which can be updated, and analysis refined as more data and information becomes available. The failure mechanisms and causes of failure are included in the analysis, especially in light of the new information, in order to develop guidelines to improve design reliability and better operating controls to reduce the probability of failure. Quantified reliability assessment based on fundamental physical behavior of components and their relationship with other components has demonstrated itself to be a superior technique to conventional reliability approaches based on utilizing failure rates derived from similar equipment or simply expert judgment.

Shah, Ashwin R.↗