Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element analysis”

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 163 records · Page 9

Optimized Gear Selection to Maximize Energy Savings in Electric Traction Drives for Medium and Heavy Duty Vehicles

Multi‑gear transmission systems are commonly used in electric traction drives for medium and heavy‑duty vehicles, while most passenger‑vehicle electric drivetrains rely on a single fixed ratio to reduce cost, weight, and complexity. Using multiple gear ratios can enable downsizing of the motor and inverter while still meeting performance requirements. Additionally, appropriately chosen ratios allow the motor to operate more frequently in high‑efficiency regions, improving overall energy usage and reducing operating costs over the drive cycle. This paper presents a systematic approach for selecting optimal gear ratios for electric drive systems. A neural‑network model is first developed to represent motor losses across the full torque–speed range using data generated from finite element analysis. This model enables fast, accurate evaluation of motor efficiency under varying operating conditions. A genetic‑algorithm‑based optimization framework is then applied to identify gear ratios that maximize energy cost savings over the drive cycle, with the resulting optimal ratios stored for real‑time implementation.

Gadiyar, Nishanth [ORNL] (ORCID:0000000348267524)↗

Optimal Hairpin Winding Configuration for EV Traction Motors to Enhance Active CMV Cancellation Capability of NPL.X Inverter

Active common-mode voltage (CMV) cancellation in neutral-point-less (NPL.X) inverter is an effective approach for reducing common-mode electromagnetic interference and alleviating voltage stress associated with high-speed switching in EV traction drives. However, the existing hairpin winding configuration often exhibits impedance imbalance due to slot leakage inductance and manufacturing constraints, limiting the effectiveness of active CMV cancellation. This paper investigates symmetric hairpin winding configurations for a dual-three-phase EV traction machine to improve impedance symmetry. Considering the manufacturing constraints of six-layer hairpin windings, four feasible symmetric winding configurations are analyzed. Two-dimensional finite-element analysis is performed to extract the impedance parameters of each configuration. The extracted impedance parameters are used to quantify the impedance mismatch between complementary phase pairs and evaluate its influence on active CMV cancellation. The results show that the fully interleaved winding configuration achieves perfectly balanced phase impedances with 0 % mismatch, whereas the mismatch in the other configurations ranges from approximately 2.4 % to 0.8 %. Consequently, the resulting total CMV overshoot is significantly suppressed to near zero voltage under high-voltage and high switching frequency operating conditions, demonstrating the importance of winding symmetry for maximizing the active CMV cancellation capability of the NPL.X inverter.

Lee, Kangbeen [Purdue Univ., West Lafayette, IN (U↗

SAVY-4000 Finite-Element Drop Test Analysis

PFE Auxiliary Systems conducted drop testing on SAVY-4000 containers to evaluate structural response under 12-foot drop conditions. In support of that effort, a finite-element modeling capability was developed to simulate drop response across multiple container sizes and impact orientations. The purpose of this work was to provide a consistent analysis framework that could support interpretation of testing, compare response trends across multiple configurations, and generate quantities of interest for later comparison with experimental data. More broadly, the analysis and testing were intended to assess whether the containers continued to perform their primary function after a 12-foot drop, namely maintaining structural integrity and containment of the contents. The modeling approach combined an implicit preload analysis with an explicit drop simulation so that each drop event began from a mechanically realistic assembled condition, including compression of the silicone O-ring. Separate models were developed for 2-quart, 5-quart, 12-quart, and 10-gallon containers. The results were evaluated in terms of strain-gauge response, collar-lid gap behavior, and accumulated plastic strain. In addition, parametric studies were performed on the 2-quart container to assess sensitivity to O-ring stiffness, friction, canister thickness, geometry tolerance, and mesh density. The simulations showed that predicted drop responses depended strongly on both container size and drop orientation. Gap metrics identified cases in which the predicted collar-lid opening exceeded the nominal O-ring cross-section threshold, while plastic strain metrics identified localized regions of elevated permanent deformation. Parametric studies showed that the predicted response was especially sensitive to the assumed O-ring stiffness and contact friction, while the geometry tolerance study produced smaller changes in the cases examined. The main value of this work was that it established a repeatable modeling and simulation workflow to support drop-test implementation, evaluate effects of future configuration changes, and understand modeling assumptions that most influenced predicted response. At the current stage, the results were viewed as preliminary model predictions rather than validated predictions. The next step would be to compare drop-test data to the model so that predictive values of the workflow could be refined and used with greater confidence to assess whether the containers maintained structural integrity and containment of the contents after a 12-foot drop.

42 ENGINEERING↗

Enabling the Broader Use of MOOSE for Nuclear Energy and Other Simulation

This Final Scientific and Technical Report summarizes work performed under the Phase IIA SBIR project “Enabling the Broader Use of MOOSE for Nuclear Energy and Other Simulation” (DE-SC0020906) from August 2023 through August 2025. The objective of the Phase IIA effort was to mature and harden capabilities developed during Phase II, with the goal of enabling practical interoperability between Coreform’s isogeometric analysis (IGA) technologies and the Multiphysics Object-Oriented Simulation Environment (MOOSE), while improving robustness, performance, and scalability for complex, nuclear-relevant geometries. Over the course of Phase IIA, the project established and validated an extraction-based interoperability pathway between Coreform tools and MOOSE. A combined mesh and matrix format was defined collaboratively with MOOSE developers and integrated into the solver, enabling standard MOOSE workflows to operate on data exported from Coreform’s IGA and Flex Representation Method (FRM) pipelines. Early demonstrations validated architectural compatibility using linear solid mechanics problems, while later efforts focused on benchmark testing and external use. By the end of the project period, engineers at BWXT were able to independently set up and execute a simulation using the Coreform–MOOSE workflow and provide direct feedback that informed further refinement. In parallel, substantial effort was devoted to improving the robustness of trimmed U-spline construction for complex CAD geometries. A growing test suite of nuclear-relevant models was compiled through collaboration with multiple stakeholders and used to drive extensive bug fixing and reliability improvements. These efforts resulted in improved robustness and performance, including the addition of fallback capabilities that enhance reliability when the underlying commercial CAD kernel fails. Performance-oriented work progressed later in the project, with the development and demonstration of methods to decompose complex geometries into structured subregions and updated data representations to support more efficient solver processing. Additionally, extensive enhancements to threadsafe parallel data structures and trimming operations established a foundation for scalable processing of large assemblies. Collaboration with Sandia National Laboratories on the SGM geometric modeling kernel advanced to a functioning interface test case, positioning the workflow for future kernel integration. Overall, the Phase IIA effort successfully transitioned the project from architectural proof-of-concept to externally exercised, solver-integrated capability, while clarifying remaining technical challenges related to standardization, performance optimization, and kernel integration.

42 ENGINEERING↗

Liquid Crystal Orientation and Shape Optimization for the Active Response of Liquid Crystal Elastomers

Liquid crystal elastomers (LCEs) are responsive materials that can undergo large reversible deformations upon exposure to external stimuli, such as electrical and thermal fields. Controlling the alignment of their liquid crystals mesogens to achieve desired shape changes unlocks a new design paradigm that is unavailable when using traditional materials. While experimental measurements can provide valuable insights into their behavior, computational analysis is essential to exploit their full potential. Accurate simulation is not, however, the end goal; rather, it is the means to achieve their optimal design. Such design optimization problems are best solved with algorithms that require gradients, i.e., sensitivities, of the cost and constraint functions with respect to the design parameters, to efficiently traverse the design space. In this work, a nonlinear LCE model and adjoint sensitivity analysis are implemented in a scalable and flexible finite element-based open source framework and integrated into a gradient-based design optimization tool. To display the versatility of the computational framework, LCE design problems that optimize both the material, i.e., liquid crystal orientation, and structural shape to reach a target actuated shapes or maximize energy absorption are solved. Multiple parameterizations, customized to address fabrication limitations, are investigated in both 2D and 3D. The case studies are followed by a discussion on the simulation and design optimization hurdles, as well as potential avenues for improving the robustness of similar computational frameworks for applications of interest.

42 ENGINEERING↗

X-point effects on the ideal MHD modes in tokamaks in the description of dual-poloidal-region safety factor

The flux coordinates with dual-region safety factor (q) in the poloidal direction are developed in this work. The X-point effects on the ideal MHD modes in tokamaks are then analyzed using this coordinate system. Since the X-point effects mainly affect the edge region, the modes localized at the tokamak edge are particularly examined. Two types of modes are studied. The first is related to the conventional peeling or peeling-ballooning modes. The mode existence aligned with the local magnetic field in the poloidally core region as observed experimentally is confirmed. The X points are shown to contribute to a stabilizing effect for the conventionally treated modes with the surface-averaged q and with the tokamak edge portion truncated. The other is the axisymmetric modes localized in the vicinity of X points, which can affect the cross-field-line transport near the X points. The existence of axisymmetric modes points to the possibility of applying a toroidally axisymmetric resonant magnetic perturbation (RMP) in the X-point area for mitigating the edge localized modes, which can be an alternative to the current RMP design. The dual q description also has important implications for the existing non-axisymmetric RMP concept. It helps to understand why the RMP suppression of edge localized modes is difficult to achieve in the double-null tokamak configurations and points to the possibility of further improving the current RMP concept by considering the alignment to the local q.

Fourier analysis↗

Insight into the interfacial microstructure and chemistry of hot isostatically pressed AA6061-AA6061 bonds for U-10Mo fuel cladding application

Here, in the present study, hot isostatic pressing (HIP) bonding was carried out using AA6061-AA6061 plates for U-10Mo monolithic fuel cladding application. HIP was performed at 560°C with an applied isostatic pressure of 103 MPa for 90 min. Interfacial microstructure characterizations were performed using scanning electron microscopy (SEM) combined with electron backscatter diffraction (EBSD) to quantify precipitate fractions and study recrystallization, grain growth, and orientation relationships across the interface. EBSD results confirmed that no interfacial grain growth occurred during the HIP process. Transmission electron microscopy (TEM) was also employed to quantify composition and identify phases of various precipitates and an oxide layer that formed at the interface during HIP bonding. Microstructural analyses revealed the formation of a large fraction of Mg 2 Si precipitate and a ~50 nm thick Mg-oxide layer at the interface. Further TEM analysis confirmed that the oxide layer was not continuous, and Mg 2 Al 2 O 5 oxide particles with an average size of 18 nm were present in the oxide layer. Finite element modeling analysis was performed to study the stress distribution at the interface during HIP processing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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↗

Effects of Material Selection on the Remanufacturability of a Swashplate in an Axial Piston Pump

A swashplate is a critical component in an axial piston pump for motion transformation that undergoes cyclic loading. Traditionally, swashplates are not designed to have multiple lives of operation. According to industry experts, a remanufacturability analysis of swashplates suggests that the scrap rate is high. Among the primary reasons for this high scrap rate is insufficient material on the running surface, especially for subtractive recovery. Therefore, the selection of the proper material and its amount plays a significant role in the fatigue behavior, the remanufacturability of the swashplate, as well as its cost. Analyses addressing fatigue life, remanufacturability, and cost of the swashplate were carried out in this paper through a comparative study utilizing various materials. Recommendations on the material selection that balance the failure consequences, remanufacturability, and cost are provided to guide designers in making informed decisions for alternative design options.

42 ENGINEERING↗

Facile Tensile Testing Platform for In Situ Transmission Electron Microscopy of Nanomaterials

In situ tensile testing using transmission electron microscopy (TEM) is a powerful technique to probe structure-property relationships of materials at the atomic scale. In this work, a facile tensile testing platform for in situ characterization of materials inside a transmission electron microscope is demonstrated. The platform consists of: 1) a commercially available, flexible, electron-transparent substrate (e.g., TEM grid) integrated with a conventional tensile testing holder, and 2) a finite element simulation providing quantification of specimen-applied strain. The flexible substrate (carbon support film of the TEM grid) mitigates strain concentrations usually found in free-standing films and enables in situ straining experiments to be performed on materials that cannot undergo localized thinning or focused ion beam lift-out. The finite element simulation enables direct correlation of holder displacement with sample strain, providing upper and lower bounds of expected strain across the substrate. The tensile testing platform is validated for three disparate material systems: sputtered gold-palladium, few-layer transferred tungsten disulfide, and electrodeposited lithium, by measuring lattice strain from experimentally recorded electron diffraction data. The results show good agreement between experiment and simulation, providing confidence in the ability to transfer strain from holder to sample and relate TEM crystal structural observations with material mechanical properties.

2D materials↗

Reinforcement induced microcracking during the conversion of polymer-derived ceramics

This study investigates the role of discontinuous reinforcement on the microcracking of a preceramic polymer matrix during polymer to ceramic conversion. The material system comprised a carbosilane-based, preceramic polymer reinforced with mullite particles. The preceramic resin slurry was formulated for photopolymerization on a digital light projection 3D printer. Micro X-ray computed tomography, using a synchrotron light source, monitored a printed composite part during pyrolysis. The conversion of the carbosilane matrix into Si(O)C generated microcracks in the matrix that radially extended from particles. The likelihood of microcrack formation positively correlated with particle size. The largest particles (>10 4 µm 3 volume or >60 µm side length) always abutted microcracks, while the matrix surrounding smaller particles (<5 µm) was free of microcracks. Numerical calculations showed the particles resist the shrinking matrix during its conversion. This created tensile stresses within the matrix. The driving force for microcrack growth was also analyzed with regard to the influence of matrix material parameters, particle form-factor and crack length. These results reveal the need to consider the form factor of discontinuous reinforcements, with regard to the material properties of the polymer-derived ceramic, in order to reduce or eliminate defects in the final part.

36 MATERIALS SCIENCE↗

High compressive energy absorption and shape recovery behavior of additively manufactured textile-inspired cylindrical braided metamaterials

Mechanical metamaterials (MMs) are engineered structures with unique mechanical properties that arise from their unique spatial arrangement or lattice-like structure. The most commonly designed MMs such as honeycomb and re-entrant auxetics are prone to failure at the sharp corners and weak joints due to the increased stress concentration under deformation. To mitigate this challenge, braided MM structures involving intertwining threads of nylon—forming curved unit cells—have been studied. These textile-inspired cylindrical braided metamaterials (CBMMs) with contrasting unit cells, namely diamond and regular CBMMs, were fabricated by 3D printing. The layer-by-layer deposited structure built by fused filament fabrication delivered an assembly of overlapped threads that are fused at the contact point. To understand deformation behavior of these MMs, finite element models were developed for various load scenarios including quasi-static compression, cyclic and creep loads at room temperature. Stress distribution, deformation mechanisms, and failure modes were analyzed and validated by experiments to analyze the geometries and associated performance. The diamond CBMMs showed stress softening at 30 % compressive strain, withstanding a load of ∼440 N, whereas the regular CBMMs at 50 % strain experienced ∼250 N. The diamond CBMMs delivered higher creep resistance under sustained load and better energy absorption under cyclic loading than the regular CBMMs. The latter, however, exhibited 94 % shape recovery in contrast to 88 % recovery in former prototype during their first cyclic load. In conclusion, this study helps design mechanical lightweight devices that endure significant sustained load and exhibit enhanced energy absorption and shape recovery characteristics in cyclic loading.

Creep↗

Multi-Objective design of interlocking metasurfaces using conditional diffusion models

Unit cell design remains a major challenge for interlocking metasurfaces, a promising joining technology for dissimilar materials, due to the complex, competing, multivariate design space and the need for rapid adaptation to varying performance requirements. This study explores Conditional Diffusion Models as a design optimization tool for interlocking metasurfaces. Given the complex, competing, multivariate design space for interlocking metasurfaces, unit cell design remains a major challenge for this joining technology. We trained a conditional diffusion model on 25,000 finite element analysis-simulated interlocking metasurface unit cells to generate designs with tailored thermo-mechanical properties (tensile strength, shear strength, and thermal conductivity) based on specified performance criteria. The model demonstrated a success rate of approximately 72 % in producing designs that met specified property bounds. The conditional diffusion model generated both thermally resistive and conductive designs, revealing clear trends in design characteristics: taller, dendritic structures were advantageous for tensile loads, while shorter, robust designs excelled in shear applications. Our findings indicate that the model's performance is more influenced by the breadth of the design space than by the quantity of training data, highlighting the importance of expansive design domains for generating innovative solutions. This work establishes conditional diffusion models as a highly efficient and adaptable tool for rapid interlocking metasurface unit cell design, paving the way for advancements in multi-material joining technologies, as well as highlighting the justification to leverage conditional diffusion models as design tools across complex design domains.

Conditional diffusion models↗

Solid face sheets enable lattice metamaterials to withstand high-amplitude impulsive loading without yielding

Owing to their ability to provide tunable mechanical responses, lattice materials are frequently studied to elucidate their response to static and dynamic loads. However, these roles are typically in opposition: static loads must be supported sufficiently far away from the onset of buckling or yielding, whereas dynamic loads are typically ameliorated by crushing of the lattice, which provides excellent energy-absorption due to the large plastic deformation accompanying densification. In contrast, this work considers the octet truss as an exemplar topology, in a structural role where it must simultaneously support static loads while enduring high-amplitude impulsive loads. This study focuses on the ability to withstand impulsive loads without yielding, an essential prerequisite to enduring dual loading. Computational studies using the ALE3D hydrocode were performed to examine the response of the octet truss under a short temporal width impulse shape associated with laser-driven shocks. A key finding was that covering the lattice with a solid face sheet and treating this face sheet thickness as a design variable allows the Taylor-like pulse to be attenuated prior to entering the weaker lattice, at the cost of added mass up front. Experimental validation was accomplished by laser-driven shock testing, using octet trusses printed out of Ti-5Al-5V-5Mo-3Cr. The results show that for a given quantity of mass, the attenuation is maximized when as much mass as possible is moved into the face sheet, leaving a more slender lattice structure. The effect of placing mass in the face sheet rather than lattice beams dominates the effect of relative density, to the point where a low-mass structure with most of the mass concentrated in the face sheet can outperform a high-mass structure with most of the mass in the lattice. Finally, by further understanding the propagation of short pulse width waves within under-dense structures, this study expand the domain of applicability of such structures, including lattice materials, to challenging dual-loading regimes spanning decades of strain rates.

36 MATERIALS SCIENCE↗

Adaptively remeshed multiphysical modeling of resistance forge welding with experimental validation of residual stress fields and measurement processes

Welding processes used in the production of pressure vessels impart residual stresses in the manufactured component. Computational modeling is critical to predicting these residual stress fields and understanding how they interact with notches and flaws to impact pressure vessel durability. Here, in this work, we present a finite element model for a resistance forge weld and validate it using laboratory measurements. Extensive microstructural changes, near-melt temperatures, and large localized deformations along the weld interface pose significant challenges to Lagrangian finite element modeling. The proposed modeling approach overcomes these roadblocks in order to provide a high-fidelity simulation that can predict the residual stress state in the manufactured pressure vessel; a rich microstructural constitutive model accounts for material recrystallization dynamics, a frictional-to-tied contact model is coordinated with the constitutive model to represent interfacial bonding, and adaptive remeshing is employed to alleviate severe mesh distortion. An interrupted-weld approach is applied to the simulation to facilitate comparison to displacement measures. Several techniques are employed for residual stress measurement in order to validate the finite element model: neutron diffraction, the contour method, and the slitting method. Model-measurement comparisons are supplemented with detailed simulations that reflect the configurations of the residual-stress measurement processes themselves. The model results show general agreement with experimental measurements, and we observe some similarities in the features around the weld region. Factors that contribute to model-measurement differences are identified. Finally, we conclude with some discussion of the model development and residual stress measurement strategies, including how to best leverage the efforts put forth here for other weld problems.

36 MATERIALS SCIENCE↗

Development of a thermal creep model for aluminum alloy 6061 cladding in U-10Mo monolithic fuel plates

Plate-type fuel elements consisting of a high-density, low-enriched uranium (LEU) U–10Mo-based fuel foil encapsulated in an aluminum alloy (AA) cladding are fabricated using the hot isostatic pressing (HIP) technique. During the HIP process, the fuel plate system is heated to 560 °C, then cooled to room temperature. This heat cycle significantly affects the mechanical properties of the aluminum cladding, and experimental investigations have shown that, post-HIP bonding, the mechanical properties of the aluminum cladding transition from those of AA 6061-T6 to something closer to the O temper. More specifically, the ultimate strength of the cladding decreases while its ductility increases, making it challenging to capture the changes in mechanical behavior and material properties. Understanding the residual stresses generated during the HIP process is critical for assessing the fuel plate’s integrity under various temperature, pressure, and irradiation. To simulate the HIP bonding process, the elastic, plastic, and thermal properties of the cladding are assumed to be similar to those of AA 6061-O temper. However, the primary challenge lies in the lack of available data for the creep model of the AA 6061 cladding during this transient process of HIP. The present study focuses on developing a computational model that predicts the creep behavior of the aluminum cladding in the fuel plates during the HIP process, as cladding creep significantly influences the residual stresses generated in U-10Mo fuel plates during HIP fabrication. Furthermore, as HIP bonding occurs at high temperatures that are nearing the melting point of aluminum, the present work considered a temperature-dependent Arrhenius-type creep model. In particular, a hyperbolic sine creep model is employed to estimate the creep properties of the as-fabricated aluminum cladding. In conclusion, the residual stresses predicted in the U-10Mo fuel when using the newly calibrated creep model closely align with the experimental measurements, validating the model’s accuracy.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Effect of H + Exchange and Surface Impurities on Bulk and Interfacial Electrochemistry of Garnet Solid Electrolytes

Contact loss and current constriction pose significant challenges at the Li metal interface of solid-state batteries. For garnet-structured Li 7 La 3 Zr 2 O 12 (LLZO), these effects are amplified by Li + /H + exchange and surface contamination reactions, which lead to conductivity losses and poor Li wetting. In this study, we utilize a variety of surface treatment processes across 37 cells to selectively induce proton exchange and contamination reactions in LLZO. The resulting bulk and surface chemistry is systematically characterized and correlated to changes in electrochemical properties. Additionally, we combine impedance analysis and finite element method modeling to deconvolute sources of impedance contributions at the Li metal interface. Specifically, we show that constriction impedance at the Li metal interface arises not solely from voids, but also from ionically-resistive surface contaminants. Further, these findings emphasize the connection between ionic conductivity and constriction, demonstrating that micron-scale ionically-resistive components increase constriction even with identical contact geometries. Finally, we leverage our comprehensive dataset to highlight unstable overpotential growth as a failure mechanism, additionally showing that the phase of a cell’s impedance is a sensitive indicator for the onset of interfacial instability. Overall, this study clarifies the impacts of proton exchange and surface contamination on electrochemical properties at the Li|solid electrolyte interface and elucidates insights that are generalizable to other solid-state battery systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗