SEARCH · Search NASA

Results for “Finite Element 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Predicting microstructurally sensitive fatigue‐crack path in WE43 magnesium using high‐fidelity numerical modeling and three‐dimensional experimental characterization

Abstract Microstructurally small fatigue‐crack growth in polycrystalline materials is highly three‐dimensional due to sensitivity to local microstructural features (e.g., grains). One requirement for modeling microstructurally sensitive crack propagation is establishing the criteria that govern crack evolution, including crack deflection. Here, a high‐fidelity finite‐element modeling framework is used to assess the performance and validity of various crack‐growth criteria, including slip‐based metrics (e.g., fatigue‐indicator parameters), as potential criteria for predicting three‐dimensional crack paths in polycrystalline materials. The modeling framework represents cracks as geometrically explicit discontinuities and involves voxel‐based remeshing, mesh‐gradation control, and a crystal‐plasticity constitutive model. The predictions are compared to experimental measurements of WE43 magnesium samples subject to fatigue loading, for which three‐dimensional grain structures and fatigue‐crack surfaces were measured post‐mortem using near‐field high‐energy x‐ray diffraction microscopy and x‐ray computed tomography. Findings from this work are expected to improve the predictive capabilities of simulations involving microstructurally small fatigue‐crack growth in polycrystalline materials.

Engineering↗

DOE OSTI2025

Numerical and experimental analysis of mechanically induced failure in electric vehicle battery modules

Mitigating thermal runaway and cell-to-cell propagation is essential for improving the safety of electric and hybrid vehicles. Enhancing digital twin capabilities to predict battery mechanical abuse is particularly critical for automotive and aerospace applications, where crashworthiness is a key concern. Understanding failure conditions and propagation in battery modules during mechanical abuse is complex due to interactions between structural deformation, heat transfer, electrochemical processes, exothermic reactions and mechanical fracture. While prior studies have focused on modeling cell-level behavior, extending these models to module or pack level is necessary for a system level understating of electric vehicle safety. This study develops coupled large deformation finite element models that simultaneously solve for electrochemistry, material failure, internal short circuit and thermal runaway propagation. The models account for mechanical and thermal interactions between lithium-ion cells and other battery components while the contact interfaces are evolving with time. Model-predicted voltage, temperature and force responses are compared with experimental data for validation. The results demonstrate that the approach captures key failure mechanisms, including thermal propagation through heat transfer, electrical propagation from short circuits in parallel-connected cells, and mechanical propagation via penetration and crack formation. These findings show that computational models are valuable tools for understanding battery module failure and providing insight that can reduce the need for extensive experimental testing.

25 ENERGY STORAGE↗

DOE OSTI2025

Integrating Crack Detection and Pipe Shape Optimization for Enhanced Sewage System Durability

Crack detection in underground reinforced concrete pipes has been essential in determining the state of stormwater infrastructure. Detection models have been implemented for detecting cracks and other defects in pipes using CCTV footage for stormwater drainage systems. In addition, Finite element models have been used to determine optimum shapes and pipe thickness for different boundary conditions such as header pipes in power plants. The concept of shape optimization emerges as a crucial factor in power plant design and operation, with the potential to maximize performance while minimizing the use of materials. Shape optimization not only enhances efficiency but also contributes to reducing the environmental footprint. This paper discusses the integration of both topics by using the cracks detected in underground pipes as boundary conditions for shape optimization of the pipes. A machine learning model has been developed which uses limited data for training and outlines the location of detected cracks. A shape optimization methodology is proposed in which ANSYS modules are used to analyze fluid flow and then optimize the shape of the pipe. The crack detection model developed has been applied to a crack detected in lab setting and machine learning model used has an accuracy of 98% using a random forest algorithm.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2024

Observation and Modeling of Dynamic Fracture Behaviors of Battery Cell Under Impact Loading Using Enhanced Representative Volume Element Concept

The burgeoning electric automobile industry has increased interest in battery safety. Battery cells experience significant mechanical stress during operation, including the impact of accidents and vibrations from driving. The potential for thermal runaway reactions in battery cells raises safety concerns. Although numerous researchers have defined the dynamic behavior of battery cells and proposed numerical models to describe it, few studies have focused on the high-strain rate mechanical impact phase correlated with the onset of fracture. In this study, we describe the dynamic behavior of pouch battery cells and propose a modeling method to study their mechanical failure under impact situations. Impact tests are conducted at various velocities and heights. To overcome numerical issues commonly encountered under rapid deformation scenarios, a new finite element model is developed based on the representative volume element model. The proposed approach efficiently simulates continuous crack propagation and brittleness behavior during impact by permitting the individual behavior of the cell components. Therefore, engineers can reliably design safer electric vehicle battery cells by measuring the properties of the cell components.

ENERGY STORAGE↗

DOE OSTI2024

Measurement of the Static Nonlinear Third-Order Elastic Moduli of Rocks: Problems and Applicability

The third-order elastic (TOE) model has been used to describe the widely observed nonlinear mechanical behaviors of earth materials. In addition to linear elastic constants ( λ , μ ), three nonlinear elastic moduli ( A , B , C ) are required for isotropic rocks. Contrary to previous research on dynamic TOE moduli, this study followed the protocol to measure strain and stress under uniaxial and hydrostatic compressive tests statically, which were later used to invert for the full set of TOE moduli for four standard rock types with differing pore structures of a porous oolitic limestone, quartz-rich sandstones, and a dense crystalline basalt. The applicability of the TOE model to characterize nonlinearity depends on the fulfillment of path-independence and small-strain assumptions. Using the measured static TOE moduli, the finite element model demonstrates that the stress in the vicinity of the wellbore is more amplified than the stress in the linear elastic case, which leads to a wider zone of rock failure around the wellbore. Due to the long-standing discrepancy between static and dynamic moduli, the rarely reported full set of static TOE moduli is necessary and will benefit future research in understanding the effect of rock nonlinearity on geophysical and geomechanical applications, such as long-term safe storage of CO 2 and generating process of geohazards.

58 GEOSCIENCES↗

DOE OSTI2026

AGR-5/6/7 Thermal Model with Non-uniform Gas Gaps

The daily as-run thermal analysis for the Advanced Gas Reactor (AGR)-5/6/7 experiment was documented in ECAR-5633 (Hawkes 2021). A finite element model was created for the entire test train, which consisted of five capsules. The fuel compacts, graphite holders, stainless-steel capsule walls, and all other major components were individually modeled, including all thermocouples (TCs). Daily heat rates for each compact and component in the model were input from daily as-run physics analyses. Daily gas compositions for each capsule were also input. The thermal conductivity of the compacts and graphite holders varied with fast neutron fluence. Gas mixture thermal conductivity was implemented using experimentally attained values from the literature. Fluence and temperature-dependent thermal conductivity was used for the graphite components and the fuel compacts. Radiation heat transfer was implemented. The model was tuned in an attempt to match the TC readings during the first cycle by adjusting the Neolube (graphite lubricant) thickness.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Deployment of neural-network-based neutron microscopic cross sections in the Griffin reactor physics application

The capability to utilize neural networks to predict macroscopic and microscopic cross section parametric spaces has been developed for the Griffin reactor physics application. The LibTorch interface enables Griffin's MOOSE-based materials to interact with LibTorch-trained models, allowing for the evaluation of complex macroscopic or microscopic cross section spaces, which are then used to evaluate the neutronic properties of the Griffin finite element model. This study benchmarks traditional ISOXML-formatted tabulation libraries against neural network-based models for 279 nuclides on 20,160 grid points for zero-dimensional and two-dimensional reactor models. Benchmark metrics include the fundamental mode eigenvalue, fission and absorption rates, and various temperature coefficients of reactivity (isothermal, fuel, and moderator). From the perspective of storage space, the complete set of LibTorch models uses 11 MB on disk, compared to the 10 GB for the ISOXML multigroup library that covers the same grid space. For the two-dimensional performance case considered in Griffin, the Torch model uses 97% less RAM than the reference ISOXML dataset while runtime increases by a factor of 3 when using the LibTorch model compared to the ISOXML dataset with multi-linear interpolation. The LibTorch model consistently yields errors within 0.01% for most analyzed quantities except for the temperature coefficients of reactivity where the maximum discrepancies are up to 0.3 $\frac{pcm}{K}$. Due to the neural network attempting to best predict quantities with no regard for a positive or negative bias for any given quantity, predictions may experience random fluctuations, resulting in both positive and negative errors. Future work will entail both depletion and coupled transient analysis to determine the predictive capabilities of Griffin with neural network-based cross sections.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

Inverse design of hypoeutectoid pearlite steel microstructures using a deep learning and genetic algorithm optimization framework

Goal-oriented microstructure design in metallic materials is a challenging task due to complex structure-property relationships. Traditional experimental and computational approaches are time-intensive and economically inefficient, limiting their applicability for large-scale design space exploration. Here, in this work, we propose an end-to-end framework that integrates deep learning models with genetic optimization to design microstructures with targeted mechanical properties. Deep learning models enable accurate forward design, while their integration with genetic optimization enables efficient inverse design within a few hours, compared to days or weeks using conventional finite element simulations. The framework combines experimental characterization and finite element modeling to analyze the influence of microstructural features on the mechanical behavior of hypoeutectoid steels. Data from both experiments and simulations are used to train the deep learning models. To demonstrate its effectiveness, we apply the framework to 0.63% carbon steel with proeutectoid ferrite and pearlite phases, commonly used in industrial applications. In this study, 2D microstructures were used for modeling, selected primarily for computational efficiency and to establish proof of concept. The framework successfully optimizes microstructures for targeted yield strength, ultimate strength, and stress concentration factors while significantly reducing computational time. Beyond hypoeutectoid steels, this scalable framework can be extended to other material systems and integrated with additive manufacturing, offering an efficient approach for accelerating microstructure design for specific engineering applications.

ConvLSTM↗

DOE OSTI2025

Simulation of Creep Deformation and Failure in Graded AM Microstructures

This report describes modeling tools and techniques developed to simulate the long-term material performance of 316H stainless steel manufactured using Laser Powder Bed Fusion (LPBF). A physics-based Crystal Plasticity Finite Element model is used to simulate creep in microstructures and to study the roles of grain morphology, porosity, and texture. We describe our modeling methodology, including an orientation-mapping technique to capture the spatially varying crystallographic orientation that results from the build conditions. Our study of microstructural features shows that AM microstructures produced by LPBF tend to creep faster in the build direction, while texture and grain boundaries strengthen the transverse directions. However, when grain-boundary porosity and the consequent cavity growth are included in the model, the transverse directions begin to creep faster. In examining texture, the results indicate that spatially varying orientation arising from the build conditions increases anisotropy in the material, making it critical to account for orientation gradients in the material to accurately model its mechanical behavior. We also describe a material-model calibration campaign in which we calibrated the constitutive model specifically for LPBF 316H stainless steel at 725℃ for both solution-annealed and as-built conditions. Finally, these tools and techniques are used to model creep in microstructures representing different regions of an LPBF material with graded microstructure, owing to intentional variation in processing conditions. The creep simulation results show good agreement with experimental data across all three microstructures, with future work planned to study rupture in the material.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Micromechanical Surrogate Machine Learning Model for Creep Deformation Modeling

Process variability during the manufacture of gas turbine engine hot section components can significantly affect the material’s resulting microstructure. In casting, for instance, geometric variation within a component (thin sections versus thick sections, radial location) influences cooling rates and the resulting grain size. The high temperature creep response is known to be sensitive to grain size owing to a diffusional creep mechanism which occurs more readily along grain boundaries. Microstructural variation correspondingly drives mechanical behavior which propagates into component scale performance uncertainty. These factors are essential when planning inspection, maintenance, and repair strategies within a reliability framework. These benefits provide opportunities to increase overall energy efficiency through refined margins. Critically, there is an opportunity to bolster existing data-driven reliability models using physics-driven process-structure-property relations. Here we present recent work establishing a framework for evaluating the probabilistic creep performance of high-temperature materials. A novel microstructure-sensitive crystal plasticity finite element model is established that captures both grain boundary and crystallographic deformation effects. The computationally expensive physics model is calibrated using a statistical approach and this high-fidelity model is subsequently used to train a computationally efficient machine learning surrogate model. The surrogate model is essential for sampling a large ensemble of simulated structure-property pair results. The ensemble data are then mined to extract salient trends to be incorporated into a microstructure-sensitive reliability model. The proposed approach represents a novel way to capture microstructure-sensitive trends from physics-based models within a modern reliability framework.

Fernandez-Zelaia, Patxi [ORNL]↗

DOE OSTI2024

Transient Chemo‐Mechanical Model of Lithium Plating Impacted by External Pressure

Abstract Lithium anodes show great promise in commercial applications, but are hindered by lithium plating and dendrite growth which cause safety concerns during long‐term cell operation. Stack pressure is experimentally observed to improve cell lifetime; however, the relationship between stress and lithium deposition has remained difficult to elucidate. In this work, a transient, 3D, finite‐element model of the evolution of a lithium anode due to stripping and plating is developed. The evolution of a microscale protrusion on the anode surface is tracked over one charge‐discharge cycle with respect to stack pressure and lithium yield strength. Lithium plastic deformation, nonconformal anode‐separator contact, and separator porosity effects are accounted for. Over the course of several hours of stripping/plating, the anode surface evolves to a similar morphology under pressure regardless of the initial conditions due to lithium plastic deformation and hardening. The rate of this evolution highly depends on the applied pressure and assumed lithium yield strength.

25 ENERGY STORAGE↗

DOE OSTI2025

Computationally efficient models for aqueous organic redox flow batteries

The rising usage of intermittent energy has garnered the need for large scale energy storage systems. Redox flow batteries (RFB) based energy storage system shows promising potential. Numerical simulations and machine learning approaches have been widely used to study RFB performance. The development of autonomous material discovery framework and digital twin of energy storage system usually needs to query cell performance through fast response models. In this study, two computationally efficient models are introduced: a physics-based analytical flow battery model (EZBattery), and a machine learning operator model (Deep Operator Network, denoted by DeepONet). Both models can provide cell performance near instantly, and prediction accuracy was systematically examined on an application of evaluating the performances of a 780 cm 2 aqueous organic redox flow battery (AORFB), using potential anolyte candidates in dihydroxyphenazine (DHP)-based family of organic materials. A validated computationally expansive 3-dimensional multi-physics finite element model by COMSOL was used as the ground truth and provided the training data set for the DeepONet. 1280 samples were generated with 10 properties to mimic the different possible anolyte candidates, and the cell performances were evaluated under 10 different combined operating conditions. The accuracy comparisons for the two computationally efficient models show that both models can provide comparable accuracy in predicting cell charging/discharging voltage curves. DeepONet can provide slightly higher overall accuracy than EZBattery with faster calculation speed, but highly relies on the training dataset. EZBattery does not need a training dataset and can provide interpretable physics-based explanations of the results, while being more flexible to adjust to adapt any different cell designs, flow battery architectures, and electrolyte materials.

Analytical model↗

DOE OSTI2025

Grain2Mesh: Mesh Generation for Grain-Scale Nonlinear Elasticity Modeling

The nonlinear hysteretic behavior of rocks under cyclic loading is a crucial area of study in geomechanics. The macroscopic response of a variety of materials has been found to be contingent upon the behavior of the micro-scale structure. This project aims to develop a functional and maintainable software package for generating a multi-phase numerical mesh and accompanying simulation files for finite element modeling used in computational mechanics solvers. Meshes generated from images often lack key preprocessing that reduces noise and prevents mesh element distortion that can increase computational cost. By incorporating user feedback throughout, grain2mesh ensures a high-fidelity mesh that can be used to model grain-scale interactions such as shearing, crack propagation, and interfacial material contrast. Scientific applications of this software include material fracturing, stress-strain analysis for natural and engineered materials, and nonlinear meso-scale analysis.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Modeling low cycle fatigue (LCF) of additively manufactured Hastelloy X using An accelerated crystal plasticity fatigue damage model

This paper presents a microstructure-based model for low cycle fatigue (LCF) behavior and life of Nickel-based alloy Hastelloy X manufactured using laser-powder bed fusion (L-PBF) additive manufacturing (AM). AM Hastelloy X, a solution-strengthened alloy, is tested at elevated temperature under fully reversed LCF conditions at different strain levels. A generalized plane strain finite element model is generated from electron backscatter diffraction (EBSD) characterization. The constitutive behavior of the material under fatigue is modeled using crystal plasticity and calibrated with both monotonic tensile and cyclic stress–strain data. The fatigue micro-crack initiation and propagation in the microstructure is modeled using a modified Chaboche fatigue damage model. An embedded boundary condition with a homogenous medium is used to apply the cyclic deformation and prevent numerically introduced over-constraints during fatigue simulation. A ‘cycle-jump’ method is used to accelerate the fatigue simulation and reduce the computational cost. The simulation results are compared to LCF experiments, showing satisfactory matches in cyclic stress behavior and number of cycles to macro-crack initiation for all applied strain ranges. In addition, the model illustrates the potential for quantifying microscale fatigue life impacting factors such as microstructure and surface roughness, which is needed to accurately quantify the reliability of AM components in service.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Seismic response of vertical dry storage casks under three-dimensional earthquake motions

Ensuring the long-term seismic safety of dry storage casks (DSCs) is becoming increasingly critical as these systems evolve from temporary to de facto permanent repositories for spent nuclear fuels. Traditional seismic soil–structure interaction (SSI) assessment methods use one-dimensional deconvolution or simplified boundary conditions to model incident waves. Although computationally appealing, simplifying assumptions may alter the seismic risk by neglecting the full complexity of three-dimensional (3D) wave propagation effects. To address this challenge, this paper introduces a novel high-fidelity computational framework that leverages the Domain Reduction Method (DRM) with perfectly matched layers (PML) to accurately transfer complex, 3D seismic wavefields from regional-scale fault-rupture simulations into local-scale finite element models of DSCs. Using broadband, physics-based ground motions from a generic M w 7.0 strike-slip event, both single-cask and multi-cask configurations were investigated under near- and far-field conditions. Emphasis is placed on capturing complex SSI, spatial variability in the ground motion, and nonlinear phenomena such as cask rocking and sliding. Numerical results demonstrate that near-field conditions, where forward directivity and fling-step effects dominate, lead to significantly higher DSC rocking and sliding. Far-field cases, by contrast, generally exhibit modest responses. Incorporating SSI tends to amplify or alter DSC response spectra and introduce response variability, which underscores the need for site-specific evaluations and robust modeling approaches to ensure the seismic integrity of DSCs in interim spent fuel storage installations.

Das, Tonmoy↗

DOE OSTI2024

Pretest Modeling A Spent Nuclear Fuel Seismic Shake Test

The U.S. Department of Energy Spent Fuel and Waste Science and Technology (SFWST) program is planning to conduct a series of full-scale shake table tests to simulate hypothetical earthquake conditions and record the response of surrogate spent nuclear fuel (SNF) assemblies in a dry canister storage system mockup. The shake table motions will represent a range of hypothetical earthquake conditions at hypothetical locations in the continental U.S. to generally define the range of mechanical loads that SNF can be expected to experience during extended dry storage periods. This paper describes the pretest predictions made with LS-DYNA models of the mockup storage systems. The test will use two dry storage system configurations, a mockup vertical concrete cask (VCC) and a mockup horizontal storage module (HSM). The test will use a production-quality canister and basket. Within the canister will be four instrumented fuel assemblies with fuel rods containing surrogate mass and 28 instrumented dummy assemblies that are intended to match the mass and outer dimension of a fuel assembly. The finite element models include models of the VCC and HSM on the shake table to calculate the system level dynamic responses and separate single fuel assembly models to calculate stress and strain in fuel assembly components. Both types of models include nonlinear behavior like rod-to-rod contact and the ability for VCC’s to rock and slide. This paper presents the expected response of the VCC, HSM, and fuel assemblies to the shake table testing that is planned to start in April of 2024. The earthquake conditions represent seismic hazards in the 2,000-to-20,000-year return period range. The test data is expected to confirm the expectation that fuel rod cladding will remain intact, fuel assembly structural components like guide tubes will remain intact, and no significant VCC sliding or tipping will occur in the range of conditions to be tested with the shake table.

Klymyshyn, Nicholas A.↗

DOE OSTI2026

Novel application of the modal strain energy technique for state-of-the-art damping predictions

Structural damping, which measures the energy dissipation of a vibrating structure, is a key modeling input for lightweight structures but is notoriously hard to predict. This work utilizes vibration-based measurements of centimeter-scale coupons and the modal strain energy approach to predict structural damping of a lightly damped structure. The approach was originally validated with panels shorter than a meter in length. This work extends the validation to a 2.75 m beam made of unidirectional and biaxial glass fiber laminates bonded by adhesive. The comparison between three-dimensional finite element model predictions and full-scale experimental measurements of damping show an average error of 5.2% for the first five modes. Additionally, the modal strain energy approach is newly applied with a one-dimensional geometrically exact beam theory model and a two-dimensional sectional analysis solver. This beam approach accurately predicts the damping behavior of the first bending modes but loses accuracy for higher order modes that are dominated by three-dimensional effects. This novel approach provides faster simulations while allowing arbitrary beam cross sections. The paper also investigates traditional and high-force dynamic mechanical analysis to measure structural damping of coupons. Both alternatives show significant errors in attempted validation against the theoretical thermoelastic damping of aluminum coupons.

17 WIND ENERGY↗

DOE OSTI2024

Simulating the effect of grain structure and porosity on creep for powder bed fusion 316H

This report details results from modeling studies focused on identifying the potential factors influencing creep performance of 316H stainless steel manufactured using Laser Powder Bed Fusion (LPBF), with a focus on experimentally-observed and predicted differences between the response of the AM material and conventionally manufactured wrought 316H. The studies presented here systematically look at the role of grain structure, porosity and texture in the anisotropic creep behavior of the material under conditions expected in high temperature advanced nuclear reactors. A physics-based Crystal Plasticity Finite Element modeling approach has allowed us to look at these factors in isolation and study their role in the material's long-term deformation behavior. An update to the precipitation model used in the constitutive framework for 316H is also presented here, with the aim of improving the accuracy of our results from this modeling effort. This work is a step in the direction of gaining mechanistic understanding of the long-term deformation behavior of LPBF 316H, helping us model its long term performance and reducing the qualification time for the material.

36 MATERIALS SCIENCE↗