Search NASA⌕ Search

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.

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

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Identification of mechanisms driving heterogeneous void growth in ductile aluminum

Void growth plays a central role in ductile fracture, yet the specific mechanisms that control this remain obscure. Classical models, such as those proposed by Rice and Tracey in 1969, are able to capture average rates of void growth, but cannot capture the heterogeneity of individual void growth. Building on recent work, the present study employs laboratory-based diffraction contrast tomography and in-situ x-ray computed tomography to investigate the effect of grain structure and other microstructural factors on void growth in an Al-2219 alloy. Crystal plasticity finite element (CP-FE) modeling is used alongside experimental data to evaluate the contributions of local mechanical states, grain orientation, grain size, and neighboring microstructural features. No strong linear relationships are found with any of the considered descriptors and void growth rate. Potential complex nonlinear relationships are explored with the use of a random forest regression model, which identifies initial void volume, void aspect ratio, local normal stress state, local shear stress state, and local equivalent plastic strain (EQPS) as features that most improve void growth rate predictions. The combination of these analyses suggests that these features should be prioritized to improve models of void growth.

Diffraction contrast tomography (DCT)↗
DOE OSTI2024

Machine-Learned Linear Structural Dynamics

The tension between accuracy and computational cost is a common thread throughout computational simulation. One such example arises in the modeling of mechanical joints. Joints are typically confined to a physically small domain and yet are computationally expensive to model with a high-resolution finite element representation. A common approach is to substitute reduced-order models that can capture important aspects of the joint response and enable the use of more computationally efficient techniques overall. Unfortunately, such reduced-order models are often difficult to use, error prone, and have a narrow range of application. In contrast, we propose a new type of reduced-order model, leveraging machine learning, that would be both user-friendly and extensible to a wide range of applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Coupled phase field damage and crystal plasticity analysis of intragranular fracture: The role of crystallographic orientation and voids

Damage evolution in engineering metal alloys at the grain scale exhibits significant microstructural heterogeneity and anisotropy. These heterogeneities create local hotspots for stress and strain localization, leading to void nucleation. Crystal orientation influences the active slip systems around voids, affecting lattice rotation and potentially forming discontinuities. At low triaxiality, voids may change shape due to lower stress, rotation, elongation, and coalescence. At high triaxiality, the correlation between crystal orientation and void growth rate becomes stronger, resembling the behavior observed in isolated single crystals. Therefore, understanding the effects of crystal orientation, heterogeneous strain, and defect evolution is crucial for single crystal fracture characterization. Here, in this work, a coupled phase-field damage (PFD) and crystal plasticity (CP) model is implemented within a finite element framework to analyze crystal deformation and failure. The CP method employs a dislocation density-based constitutive model, while intragranular failure is modeled using an anisotropic PFD method. The PFD model considers both the stored energy due to elastic stretching and the energy release due to defect formation and crack formation. A single crystal Al2219 with an intracrystalline spherical void is chosen to analyze fracture. The study finds that fracture propagation is strongly correlated with crystal orientations. This coupled CP-PFD model provides accurate failure prediction in crystalline materials by incorporating the effects of crystal orientations and existing voids. This study demonstrates how the local microstructure and defects influence plastic deformation and failure mechanisms in metal alloys.

Aluminum alloy↗
DOE OSTI2024

Multilaminate Energy Storage Films from Entropy‐Driven Self‐Assembled Supramolecular Nanocomposites

Abstract Composite materials comprising polymers and inorganic nanoparticles (NPs) are promising for energy storage applications, though challenges in controlling NP dispersion often result in performance bottlenecks. Realizing nanocomposites with controlled NP locations and distributions within polymer microdomains is highly desirable for improving energy storage capabilities but is a persistent challenge, impeding the in‐depth understanding of the structure–performance relationship. In this study, a facile entropy‐driven self‐assembly approach is employed to fabricate block copolymer‐based supramolecular nanocomposite films with highly ordered lamellar structures, which are then used in electrostatic film capacitors. The oriented interfacial barriers and well‐distributed inorganic NPs within the self‐assembled multilaminate nanocomposites effectively suppress leakage current and mitigate the risk of breakdown, showing superior dielectric strength compared to their disordered counterparts. Consequently, the lamellar nanocomposite films with optimized composition exhibit high energy efficiency (>90% at 650 MV m −1 ), along with remarkable energy density and power density. Moreover, finite element simulations and statistical modeling have provided theoretical insights into the impact of the lamellar structure on electrical conduction, electric field distribution, and electrical tree propagation. This work marks a significant advancement in the design of organic–inorganic hybrids for energy storage, establishing a well‐defined correlation between microstructure and performance.

Li, He↗
DOE OSTI2025

Control of machining-induced residual stress via tool geometry and process parameter modification

Distortion generated in machined, monolithic, thin-walled aerospace components due to residual stresses leads to significant material and economic waste in the manufacturing industry. Inherent residual stress (IRS) present in stock materials combines with machining-induced residual stress (MIRS) to influence the final machined part distortion. It is hypothesized that MIRS can be controlled, based on the part geometry, through deliberate cutting tool geometry and process parameter modifications to negate the effect of IRS on distortion, consequently resulting in distortion-free parts. A finite element (FE) orthogonal cutting model is developed to study how tool geometry and process parameters influence machining-induced residual stress (MIRS). Orthogonal cutting experiments are performed on Al 7075-T651 samples to measure cutting forces and MIRS. A cutting force dynamometer is used to measure forces during cutting and a novel digital image correlation (DIC) based hole drilling technique is employed to measure the near-surface residual stress (RS) in the cut samples. These data are subsequently utilized to validate the FE prediction model. Various cases of cutting simulations involving different depths of cut, tool tip radii, and rake angles are performed to study their effects on RS. Similar to prior literature, increasing the depth of cut, tool tip radius, or rake angle is found to promote the formation of near-surface tensile stresses. The competing effects of material plowing and temperature are shown to determine the type of RS at the end of the cut. Moreover, a window of variation of RS (up to ± 400 MPa) is estimated within the given range of conditions, allowing for the control of MIRS through tool and process modification.

Mathews, Ritin [ORNL] (ORCID:0000000301440828)↗
DOE OSTI2026

An International Round-Robin Study on Thermoelectric Module Testing and Development of Standard Power Generation Modules

An international round-robin study on thermoelectric power generation modules was conducted with nine participating laboratories. Two types of commercially available bismuth telluride modules, 30 mm × 30 mm and 40 mm × 40 mm, were used. A test protocol was followed with five temperature set points from 50°C to 150°C. Graphite sheets were used as thermal interface materials with test pressure at 100 psi (0.69 MPa). The results showed large lab-to-lab variations and the key source of uncertainty for module efficiency was identified as the heat flux measurement. In the meantime, significant uncertainty was also found in maximum electrical power (P max ) measurements. As a result of the round-robin, a “standard module” with 4 × 4 legs on a 20 mm × 20 mm platform was suggested. A skutterudite module and a half-Heusler module were produced with identical geometry and 4 mm × 4 mm × 8 mm legs. All transport properties to calculate the figure-of-merit, zT, were measured from ambient temperature to 500°C. Module performance was measured by two laboratories. Two finite-element-analysis (FEA)-based models were developed independently to simulate and predict the module performance. In conclusion, the standard modules eliminated significant test uncertainties and are aimed at assisting device design and achieving more accurate performance predictions.

Round-robin↗
DOE OSTI2024

Electrification Options for Multi-Family Water Heating in Cold Climates - Final Report

In multi-family buildings, large water storage tanks in centralized domestic hot water (DHW) systems can serve as thermal energy storage (TES) batteries to mitigate grid impact. These systems offer demand shift and efficiency benefits, significantly reducing peak power consumption, particularly in cold climates. This study evaluates the load-shifting benefits of a centralized heat pump water heater (HPWH) system equipped with a CO2 heat pump in multi-family buildings through simulation. The heat pump system and water storage tank are sized using design-day sizing. A finite-element-based stratified tank model and CO2 heat pump performance map from a commercial DHW product are used. Annual simulations are conducted to assess the benefits of the centralized DHW system for energy efficiency improvements, load shifting, and emission reductions. These simulations incorporate utility tariffs and marginal grid emission data from Los Angeles and Chicago. In Los Angeles, using a water tank as a thermal battery achieves 7.4% utility cost savings and 10.2% emission reduction. In Chicago, compared to HPWH conventional operation without preheating, TES-enabled central HPWH provides 15% utility cost savings and 13% emission reduction. The case study demonstrates that the demand reduction potential of central CO2 HPWHs is significant in cold climate regions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Modeling the formation of Sedan Crater using the FLAG and HOSS codes

Numerical modeling of explosion crater formation requires accounting for complex physical processes. Numerical validation of explosion cratering is an important step in modeling and requires experimental data for comparison. Models using discrete elements and continuum models have both benefits and drawbacks to their approaches. In this work, we consider both an arbitrary Lagrangian–Eulerian (ALE) hydrocode and a finite discrete element method (FDEM) approach to modeling the formation of the Sedan crater, the largest human-made crater in the United States. The Sedan crater formed from an underground nuclear detonation in the Nevada desert as part of Project Plowshare. Our models show that the continuum approach of the hydrocode matched well compared to early test time prior to the mound rupture and subsequent fireball venting, when most of the alluvium exhibited fluid behavior. Our FDEM approach matched the final crater dimensions well, after material had settled back into the crater, when material strength and solid mechanics play key roles. Our work shows how leveraging the benefits of multiple numerical approaches can lead to better understanding of complex physical problems, especially problems with limited experimental data. By using a continuum approach to early-time hydrodynamics and an FDEM approach to later-time solid mechanics, we can better understand the different physical regimes of explosion crater formation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling of Vertical Motor-driven Pump for Simulation of a Fault Signature \\ for Condition Monitoring

As part of the ongoing effort to transition from preventive maintenance strategies to condition-based maintenance strategies in nuclear power plants, there is significant reliance on using machine learning techniques. To develop a robust machine learning model that can diagnose all the fault modes of a vertical motor-driven pump, data capturing the unique signature of each fault mode is required. In practice, it is difficult to collect or capture data that captures all the fault modes from a single plant site. So to address this situation, a computational model of a vertical motor-driven pump is developed using the multipurpose finite element software COMSOL Multiphysics. The developed model is used to generate simulated data under normal operation and is compared with the vibration data collected using vibration sensors. Once the simulation model is verified under normal operating condition, simulated data for the fault mode for which minimal or no evidence is available in historical plant process data is developed. This simulated data is used to develop fault signatures to achieve robust predictive models. This paper presents modeling details and verification of the model that can used to generate data for fault modes that are not available at a plant site for condition monitoring purpose.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
DOE OSTI2025

Macro-micro multiscale modeling to assist the design of HPDC Al castings microstructure and alloys for EV super-large body structures (Phase 1)

Implementation of High Pressure Die Casting (HPDC) Aluminum (Al) body structures for high volume electrified vehicles (EV) to improve electric efficiency remains a key strategy within many original equipment manufacturer (OEM)s. In addition to high strength for safety requirements, superior Self-Piercing Riveting (SPR) performance is demanded for HPDC Al alloys to be compatible with high volume SPR joining. In this work, it is proposed to extend and validate an existing Contractor finite element multiscale macro-micro modeling approach to quantify the influence of the microstructure of HPDC alloys on the fracture strain/displacement under 3-point bend and clinch testing. The success of this work will allow to replace solution treatment stage with low energy consumption heat treatment (HT) processes, or to design new non heat treatable (NHT) HPDC Al alloys to eliminate HT requirements. Ultimately, this project will facilitate the application of HPDC Al alloys for super-large vehicle structures to significantly reduce vehicle weight, and thus improving energy efficiency. The purpose of this project is to extend and validate an existing finite element code, which is based on the Contractor developed macro-micro multi-scale modeling approach, to numerically simulate the three-point bending and clinch test and study the influences of material microstructural characteristics and phase properties on the rivetability. The macro-micro modeling approach begins with a sample scale model and identify the location, which is mostly prone to failure, the deformation history of the boundaries of that location calculated will be used to drive a microstructure-based sub-models where the material microstructure and microscale properties are considered. Using this approach, the wrap-bending failure for two Al alloys are correctly predicted for the first time. This will start with phase I effort of building a framework of macro-micro three point bending test and clinch test of Al10SiMgMn HPDC alloy in the as-cast and T7 heat treated conditions. Those results will then be validated with experimental test results. The phase II effort will involve the utilization of the knowledge learned in phase I to establish the quantitative correlation between the microstructure characteristics and the riveting performance, which will be further used to guide the optimization of HPDC Al alloy microstructure using heat treatment process to achieve sufficient rivetability to join large thin-wall HPDC alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Review of recent activities with MOOSE, an open-source finite element & finite volume multi-fidelity simulation framework

Modeling and simulation are an increasing part of engineering. This is undoubtedly driven by the high costs of constructing experimental facilities, but also enabled by the exponential increase in computing powers over the last decades, which allows computational models to be closer than ever to reality. One of the main drivers for the development of MOOSE is supporting advanced nuclear reactor simulations. A challenging aspect of modeling advanced nuclear reactors is the plurality of physics involved, including neutronics, thermal hydraulics and fuel performance. These physics are all coupled to some extent and are generally solved in a sequential but iterative fashion. The United States (U.S.) national laboratories have been developing MOOSE, an open source multiphysics framework since its inception at the Idaho National Laboratory (INL) in 2008. This framework enables seamless coupling of multiphysics simulations and facilitates the implementation of new physics and material governing laws. It is continuously expanded with novel numerical methods and new pre-implemented physics module. Numerous applications, developed within the Department of Energy (DOE) laboratories, academia, and industry, including outside of nuclear engineering, have been developed to study specialized physics problems. International collaborations are welcome on this open-source modeling and simulation project.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Machine learning-enabled multiscale modeling of mechanical deformation of aluminum and Al-SiC nanocomposites

A machine learning-enabled multiscale framework is developed for modeling the mechanical response of both pure metal and nanoparticle-reinforced metal matrix nanocomposites (MMNCs). Using aluminum–silicon carbide (Al-SiC) as an example MMNC, atomistic simulations reveal three distinct deformation mechanisms (i.e., defect-free, dislocation-based, and interface separation) governed by the interfaces between the Al matrix and SiC nanoparticles. As compared with single crystal Al, the lattice undergoes a more abrupt failure once the dislocation network becomes extensive and void nucleation initiates, whereas in Al-SiC, nanoparticle interfaces enable a more gradual progression of damage. These mechanisms are captured through a combined classification-regression neural network surrogate model that bridges atomic-scale insights with continuum-scale finite element analysis. Machine learning-enabled multiscale modeling of pure Al accurately predicted strain localization and confirmed by in-situ scanning electron microscopic tensile testing on perforated Al specimens. This study underscores the promise of integrating physics-informed machine learning with hierarchical modeling to capture the interface dominated phenomena and guide the design of advanced MMNCs.

Al-SiC↗
DOE OSTI2025

Feasibility of an Accelerometer-Based Structural Health Monitoring System for the LANL Blast Tube

A modeling- and simulation-based study was conducted on the feasibility of implementing an accelerometer-based SHM system on the Los Alamos National Laboratory blast tube. A blast tube experiment was modeled using the Abaqus explicit finite element solver. A custom user subroutine was written to apply test-like pressure loading to the inside surface of the blast tube. The subroutine applies analytically defined pressure loads derived from tracer output taken from a Compressible Flow Computational Fluid Dynamics Solver model of the blast tube. Five unique versions of the model were created: an undamaged reference model at 65°F was used as the baseline and compared to equivalent models at 10°F and 100°F. These three models were compared to models with small damage at the reference temperature. The two types of damage considered were a radial (circumferential) crack in the main tube body and a longitudinal crack in the supports. Acceleration outputs were extracted from accelerometer bodies included in the model and were post processed using a variety of standard SHM techniques. Different potential features signaling failure were extracted and compared using statistical methods in the time and frequency domains. A method was identified that clearly shows that differences in structural response resulting from the modeled damage can be differentiated from the structural response resulting from changing environmental conditions. However, the amount of damage applied to create observable differences in the accelerometer data was so large that simpler methods of damage detection would be more cost effective in locating damage.

42 ENGINEERING↗
DOE OSTI2025

Curvature effects in the reaction-diffusion system governing frontal polymerization

Frontal polymerization is a fast and energy-efficient curing process of thermoset polymers involving the self-sustained propagation of a sharp front that transforms a liquid monomer or gel into a polymer. The temperature, speed, and geometry of the polymerization front are closely intertwined quantities that can affect the curing process and the mechanical properties of the polymer. This work investigates the effect of the front curvature on the front propagation speed based on a reaction-diffusion model. The proposed analytical relations between front speed and curvature are verified by finite element simulations and validated through experiments. The models show that a concave front has a higher front temperature and speed than its convex counterpart. When the curvature of a convex front reaches a critical value, the reaction can no longer be sustained, and the front is quenched. The relation between the speed of a steadily propagating front and its geometry is derived, showing good agreement with experimental observations. This study also provides analytical estimates of the thickness of the polymerization front and insight into the relationship between the front curvature and the kinetics of the exothermic reaction.

42 ENGINEERING↗
DOE OSTI2025

Regularizing the linearly extrapolated BDF2 scheme for incompressible flows with time relaxation

This paper presents a highly-efficient finite element scheme for the time relaxation model (TRM). The efficiency is achieved through the second-order BDF2 time-stepping scheme with linear extrapolation (BDF2LE). The accuracy of the scheme is also greatly enhanced through the use of the divergence-free Scott-Vogeulis finite elements, and van Cittert approximate deconvolution. A complete finite element analysis is provided, which includes rigorous proofs for the stability, well-possessedness, and convergence of both velocity and pressure solutions. Furthermore, we also demonstrate that the inclusion of the linear time relaxation term preserves the long-time stability of the unregularized BDF2LE scheme. Finally, numerical experiments are presented that demonstrate the added stability and accuracy that time relaxation can provide.

97 MATHEMATICS AND COMPUTING↗