Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite Element 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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Deep Neural Network Assisted Distributed Strain and Temperature Fiber Sensor System for Natural Gas Pipeline Monitoring

Natural gas pipeline integrity monitoring is crucial to detect potential leaks, find structural issues, and prevent environmental damage. This article presents a system of natural gas pipeline monitoring that uses a specialized double Brillouin peak sensing fiber along with the Brillouin optical time domain analysis (BOTDAs) technique. The calibrated sensing fiber coefficients for strain and temperature are 41.8 kHz/ με and 0.9 MHz/°C for peak 1; and 47.2 kHz/ με , and 1.11 MHz/°C for peak 2, respectively. Initially, lab tests were performed by installing a short section of double Brillouin peak fiber (DBPF) on a 1-in steel pipe under pressure up to 1000 per square inch (psi) at elevated temperatures. Simultaneous distributed measurements of temperature and pressure-induced hoop strain were successfully measured. Considering the long processing speed to extract Brillouin frequency shift (BFS), we employ a novel probabilistic deep neural network (PDNN) framework for rapid BFS prediction. Additionally, using the Finite Element Method, the effects of the pipeline pressure on hoop strain were modeled and compared to the experimental hoop strain under the same set of pipeline conditions. Finally, an actual 4-in outer diameter steel natural gas pipeline was used for pilot-scale tests, where hoop strain was measured at various pressure levels. Leaks were simulated to demonstrate accurate pipeline integrity monitoring. At an internal pipe pressure of 1000 psi, hoop strain of approximately 300 με was observed, and the sensitivity was calculated as 0.28 με /psi. The results of this pilot-scale study demonstrated that the system is capable of performing distributed monitoring sufficient to detect pipeline pressure and the presence of leaks to ensure the safe operation of gas pipelines in the field.

03 NATURAL GAS↗
DOE OSTI2024

Transfer Tube Material Properties Obtained from Test Data

A series of tension and compression tests were performed on typical transfer tubes made of 304L corrosion-resistant steel. The results were compiled and converted to true stress and true plastic strain for comparison to numerical models. The parameters of a Johnson-Cook constitutive model and fracture criterion were calibrated to mimic the test data. Comparisons were made between the converted data and the empirical formulas in the Johnson-Cook model. Simulations were performed to evaluate the implementation of the formula in finite element analysis, and the parameters were adjusted to improve the correlation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Pre-metered coating flow models with Goma 7: Workflow Tutorial

Tutorials for modeling of slot-die and slide-die coating flows with Goma 7, an open source finite element code, are presented. The tutorials cover the workflow to attaining steady state solutions for these flows, and continuation strategies for navigating the operating windows. Advanced topics of coating window prediction, automated multiparameter continuation, non-Newtonian rheology, dynamic contact line modeling, and some more solution strategies are also covered.

08 HYDROGEN↗
DOE OSTI2024

Enhancing the cooling performance of thermocouples: a power-constrained topology optimization procedure

Abstract Heat pumping through thermoelectric devices has many advantages over traditional cooling. However, their current efficiency is a limiting factor in their implementation. In this paper, we approach the non-convex topology optimization of thermoelectrical elements for cooling applications through the method of moving asymptotes (MMA) to improve their cooling capabilities per watt usage. The optimization problem is defined for a given power budget, aiming for the minimum temperature with a known heat pumping need. The introduction of power as a constraint justifies the introduction of the voltage gradient across the thermocouple as a design variable to maintain the thermoelectrical device in its optimum power-to-heat extraction ratio. To better understand the convergence of this non-convex problem, we present a two-variable analytical thermoelectric optimization model. This example provides information on how to select the penalty parameters used to scale the three material coefficients involved in the problem to obtain lower objective values and better convergence using MMA. The analytical model shows the non-convexity of the problem and provides the recommendation to use penalization coefficients of the form $$p_k=p_{\sigma }>p_{\alpha }=1$$ p k = p σ > p α = 1 for the thermal conductivity, electrical conductivity, and Seebeck coefficients. We tested these penalization coefficients through optimizations of a model based on the 1MC10-031 commercial thermoelectric-cooler (TEC) using the finite element method (FEM). These penalization coefficients provided local minima without the need for volume constraints. With this procedure, we found designs that provided temperatures close to 10 degrees lower using 60% less semiconductor material volume compared to the initial design.

Gutiérrez, G. Reales↗
DOE OSTI2025

Calibrating a finite-strain phase-field model of fracture for bonded granular materials with uncertainty quantification

To study the mechanical behavior of mock high explosives, an experimental and simulation program was developed to calibrate, with quantified uncertainty, a material model of the bonded granular material Idoxuridine and nitroplasticized Estane-5703. This paper reports on the efficacy of such a framework as a generalizable methodology for calibrating material models against experimental data with uncertainty quantification. Additionally, this paper studies the effect of two manufacturing temperatures and three initial granular configurations on the unconfined compressive behavior of the resulting bonded granular materials. In each of these cases, the same calibration framework was used; in that, hundreds of high-fidelity direct numerical simulations using a new, graphics processing unit-enabled, high-performance finite element method software, Ratel, were run to calibrate a finite-strain phase-field fracture model against experimental data. It was found that manufacturing temperature influenced the elastic response of the mock high explosives, with higher temperatures yielding a stiffer response. By contrast, it was found that the initial configuration of the grains had a negligible impact on the overall behavior of the mock high explosives though it remains possible that local damage accumulation within the specimens could be altered by the initial configurations. Overall, the calibration framework was successful at creating well-calibrated models, showing its usefulness as an engineering and scientific tool.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simulation Tools for Characterizing Stress Distribution in Laser Welded Dissimilar Joints

This project focuses on developing a thermo-metallurgical-mechanical modeling method to accurately predict the microstructural evolution and residual stress in laser welding between dissimilar metals, such as HSLA steel and high carbon equivalent (CE) gear steel. The method leverages a comprehensive material database to model the temperature and rate dependent phase transformations, along with their associated effects on material properties, such as thermal expansion and flow stress, throughout the welding process. A key innovation is the incorporation of phase transformation and phase-specific properties, which enhances the accuracy of residual stress predictions. The mixture material in the fusion zone due to the dissimilar metals will also be addressed in the numerical model. This is especially critical in scenarios involving phase transformations in the fusion zone and heat-affected zone (HAZ), where the phase changes can induce substantial residual stress variations. The material database has been generated using JMatPro. The modeling approach is implemented through a custom User Material (UMAT) subroutine, executed with the commercial finite element software Abaqus.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Mechanical Characterization of the NIF Ignition Target TMPSA Bonding Flexure

The Thermo-Mechanical Package Sub-Assembly (TMPSA) provides a critical mechanical, thermal, and electrical interface between a silicon arm and a TMP aluminum can in a NIF ignition target. During assembly, sixteen silicon pads are bonded to the aluminum using a fixture that positions the components and applies a repeatable prescribed displacement through a compliant flexure. The flexure converts fixture interference into displacement and reaction force. Because bondline thickness variation must be maintained within +/-4 µm, consistent flexure behavior is important to the assembly process. With a recent string of TMPSAs exhibiting low bond strength, the flexures were inspected. Despite being manufactured to the same specifications, flexures were found to exhibit variation in measured stiffness. Additionally, the measured stiffness values did not always follow the presumed beam mechanics model. This work addresses two questions: whether the bonding fixture is working as intended, and whether the previously made stiffness measurements are accurate representations of the use case. The flexures are analyzed using Euler-Bernoulli beam theory, measured beam dimensions, finite element analysis, and tolerance stack-up calculations. The analysis shows that the fixed-guided beam mechanics model appropriately represents the flexure during TMPSA bonding, but the chisel-tip stiffness measurement method introduces a deformation to the inner ring of the flexure that is not represented during use. The measured stiffness values should therefore be interpreted as test-condition stiffness values rather than direct measurements of operational flexure stiffness. The discrepancy is therefore attributed primarily to the measurement boundary condition rather than to failure of the fixed-guided beam model. Recommendations are provided for GD&T, dimensional inspection, and a representative stiffness testing method to better control bondline variation.

42 ENGINEERING↗
DOE OSTI2025

Assessment of buffer-IPyC thermomechanical debonding behavior using new experimental strength data in BISON

TRIstructural ISOtropic (TRISO) fuel is a nuclear fuel commonly used in High Temperature Gas-cooled Reactors (HTGRs). A single sub-millimeter-diameter TRISO fuel particle consists of a spherical fuel kernel surrounded by four coating layers: a low-density pyrocarbon buffer layer, an inner pyrolytic carbon (IPyC) layer, a silicon carbide (SiC) layer, and an outer pyrolytic carbon (OPyC) layer. The kernel is commonly made of UO2 or a mixture of uranium carbide and uranium oxide (UCO). During reactor operation, the TRISO coating layers are subjected to irradiation-induced dimensional changes and the associated thermomechanical behavior of each layer. One of the observed behaviors is gap formation between the buffer and IPyC layer due to the porous buffer’s irradiation-induced shrinkage exceeding that of the IPyC layer. Not all irradiated particles will experience buffer-IPyC gap formation. The debonding may be partial, or it may be nearly total. However, from post-irradiation examination of UCO TRISO fuels irradiated as part of the Advanced Gas Reactor (AGR) Fuel Development and Qualification Program, it was concluded that partial buffer-IPyC debonding was the most common type of buffer-IPyC interaction. To predict TRISO thermomechanical performance, multi-physics models have been built that are being continually updated and refined. The BISON code is a finite element-based nuclear fuel performance code that may be used for 1D, 2D, and 3D TRISO particle simulations. This code is used to calculate fuel temperature, kernel swelling, buffer densification, thermal and irradiation creep, fracture, and fission gas production and release. One of the recent additions to the BISON code is the ability to model the process of layer debonding. This paper will focus on the simulation results of the improved BISON debonding model that will utilize updated strengths measured from irradiated AGR TRISO fuel particles. The new experimental strength data from micromechanical tests of irradiated TRISO fuel samples were exercised in the BISON simulations and compared to baseline strength data to assess their applicability in the models. This also includes updated buffer-IPyC bond strengths to simulate layer delamination. Based on current experimental observations it is noted that the buffer-IPyC separation occurs not exactly at the junction of these two layers, but more on the side of the buffer layer. That observation is also implemented in the TRISO interface debonding model. This improved modeling approach using experimental strength data to characterize buffer-IPyC debonding and its potential subsequent cracking will be presented in the paper along with comparisons to available experimental observations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

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

Instabilities and phase transitions in architected metamaterials: a gradient-enhanced continuum approach

Architected metamaterials such as foams and lattices exhibit a wide range of properties governed by microstructural instabilities and emerging phase transitions. Their macroscopic response–including energy dissipation during impact, large recoverable deformations, morphing between configurations, and auxetic behavior–remains difficult to capture with conventional continuum models, which often rely on discrete approaches that limit scalability. In this work, we propose a nonlocal continuum formulation that captures both stable and unstable responses of elastic architected metamaterials. The framework extends anisotropic hyperelasticity by introducing nonlocal variables and internal length scales reflective of microstructural features. Local polyconvex free-energy models are systematically augmented with two families of non-(poly)convex energies, enabling both metastable and bistable responses. Implementation in a finite element framework enables solution using a hybrid monolithic–staggered strategy. Simulations capture densification fronts, forward and reverse transitions, hysteresis loops, imperfection sensitivity, and globally coordinated auxetic modes. Overall, this framework provides a robust foundation for accelerated modeling of instability-driven phenomena in architected metamaterials, while enabling extensions to anisotropic, dissipative, and active systems as well as integration with data-driven and machine learning approaches.

42 ENGINEERING↗
DOE OSTI2025

Effect of toolpath in large-format additive manufacturing with bio-derived composites

There has been growing interest in integrating bio-derived composites into Large-format additive manufacturing (LFAM) feedstocks to reduce the use of petroleum-derived materials and reduce the overall carbon footprint of LFAM. However, these materials present unique challenges during manufacturing due to their variability, which can lead to unintended deformations and failures attributable to suboptimal process conditions. While numerical modelling has been extensively employed to simulate numerous manufacturing processes, its application in LFAM with bio-derived composites remains limited. This study addresses this gap by systematically developing a numerical model to simulate the LFAM process using bio-based materials, specifically wood fiber-reinforced polylactic acid (PLA/WF). Experimental investigations were conducted to characterise the thermal and mechanical properties of additively manufactured PLA/WF specimens. Numerical simulations were performed to predict temperature profiles and deformations during LFAM. The effect of varying infill patterns, internal structures, and tool paths on the temperature distribution and deformation of printed parts was explored using the developed model. This article aims to advance the utilisation of bio-derived composites in LFAM systems and provide a comprehensive understanding of the LFAM process. The findings offer valuable insights for optimising process parameters and enhancing the performance of LFAM with bio-based composites.

Large-format additive manufacturing↗
DOE OSTI2024

NSTXU Diagnostic Disruption Dynamic Loading Represented by Response Spectra

This article presents the results of transient dynamic simulations of loads due to disruption eddy currents on the NSTXU vacuum vessel. Dynamic loading at diagnostic mounting locations is expressed as response spectra derived from the time history results of the dynamic structural simulations of a variety of disruption scenarios. The disruption simulations draw on a history of the project assessments of worst case disruptions for specific components. Major efforts to assess disruption loading have included the vacuum vessel which is the major structural support for the machine, as well as the passive plates (PPs), high harmonic fast wave (HHFW) antenna, and centerstack casing. Each one of these efforts included transient electromagnetic simulations producing time-dependent eddy current Lorentz loads (and in some cases halo loads) which then were applied to time-dependent structural dynamic analyses intended to obtain the proper dynamic amplification factors. In some instances, the EM model and structural model were identical allowing direct transfer of EM forces to the structural model. In other cases, the EM and structural model were not identical and the vector potential (VP) transfer method was used. The results files from these analyses were available (or re-run) to post process in ANSYS Classic time history postprocessor. In conclusion, the ANSYS command is used to create response spectra from time history data at desired points on the vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Development of Graphite Thermal and Mechanical Modeling Capabilities in Grizzly

Nuclear-grade graphites are used extensively in the core designs of multiple types of advanced nuclear reactors. In the reactor environment, graphite is exposed over long durations to extreme conditions, including high temperatures, radiation and potentially molten salt and oxygen. Exposure to these conditions can cause several degradation mechanisms in graphite, including nonuniform volumetric strains induced by irradiation and thermal expansion, which lead to stresses that can compromise the performance of graphite components. Evaluating component integrity, predicting component performance over the reactor lifetime, and developing design standards all require robust tools for predicting fracture initiation and propagation in graphite structural components in nuclear reactors. This report documents progress in an ongoing effort to develop modeling and simulation tools in the Grizzly code for predicting the performance of graphite exposed to reactor conditions. Recent developments include a set of thermal and mechanical models that now include the IG-110, NBG-18, and H-451 graphite grades. Improvements have also been made to a nonlinear damaged plasticity model applicable to predicting damage under tension and compression to quasibrittle materials, including graphite. In addition, enhancements have been made to the extended finite element method implementation targeted at simulating graphite fracture. These include new capabilities for crack nucleation in the interior of a solid body, improved treatment of crack nucleation on free surfaces, and more robust modeling approaches for crack growth approaching free surfaces or other cracks.

36 MATERIALS SCIENCE↗
DOE OSTI2026

ACT University Collaboration Proposal: Role of manufacturing defects on material failure under dynamic loading for developing enhanced failure models and theories (Final Report)

Research on the tensile behavior of additively manufactured 316L stainless steel coupon specimens at increasing strain rates was conducted over the past year at Penn State. Dynamic loading rates into the 1000 1/s loading rates were performed on specimens with nearly 1.32” gage lengths. Stress-strain plots show ductile and plastic behavior well beyond 12% strain with necked specimens, having smaller effective gage lengths showing up to 56% ultimate strain. Additional tests performed on compact tension specimens helped with simulation work to understand the deformation behavior. Using finite elements, it was possible to determine the feasibility of a comprehensive experimental validation study towards an improved model for AM material failure – specifically, the Bai-Wierzbicki approach which accounts for lode angle and triaxiality. A non-significant number of tests are projected for eight specific failure nodes, with additional replicates to provide strain-rate capabilities to the existing formulation. One such correction factor is explored for quasi-static, notched specimens. Finally, stress intensity factor was explored using a set of compact-tension specimens which would also provide useful validation data for any simulation work.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mechanistic Multiscale Uncertainty Propagation in Support of Accelerated Fuel Qualification

Taking a nuclear fuel concept through the research, development, and qualification stages has historically taken on the order of 20 to 25 years because of extensive irradiation tests required for a variety of conditions. The concept of accelerated fuel qualification (AFQ) has been proposed to increase the innovation pace for nuclear fuels. The goal of AFQ is not to replace the traditional qualification approach but rather to reduce the total number of experiments required to ensure approval from the regulatory authority. Of the many AFQ approaches being explored, advanced modeling—and, in particular, mechanistic modeling—is in a uniquely cross-cutting position to reduce the number of required integral tests through the inclusion of separate-effects testing, while helping to extrapolate reactor performance during rare events. We make the case that propagation of uncertainty through various computational length scales helps contextualize mechanistic modeling. We will utilize UO 2 fission gas diffusion predictions from the atomistically informed cluster dynamics code Centipede to inform fuel performance rodlet simulations using the BISON finite element code as the metric for showing how multiscale mechanistic uncertainty quantification can help reduce uncertainty in fuel performance. In conclusion, by quantifying uncertainty and its reduction through multiscale modeling, the qualification process may be accelerated through the reduction of costly irradiation experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Reduced Order Model for Guided Wave Propagation on Gas Pipelines to Enable Real-Time Simulation

Reduced order model for simulation of Guided wave propagation is presented here. The utilization of reduced order models ensures efficient data generation for a variety of parameters where it takes huge computational effort to simulate, crucial for timely monitoring and decision-making. Autoencoder based reduced order models are proposed here, which are trained on simulated data from open-source finite element framework, Firedrake.

Bukka, Sandeep Reddy↗
DOE OSTI2025

SULI Intern Final Report: Computationally Investigating Hydrogen Thermo-Diffusion in Yttrium Hydride Using Multiscale Methods

The renaissance of nuclear energy has arrived, heralding an age of abundant inexpensive clean energy, and renewed space exploration. In nuclear-powered spacecraft and microreactors, safety and size are of utmost importance. Yttrium Hydride (YHx) is being researched for its utility as a neutron moderator in nuclear reactors; the hydrogen in YHx slows down neutrons, enabling a continuous nuclear reaction in the reactor. This has the benefit of allowing reactors to be more safe, compact, and efficient. The goal of this effort is to computationally predict the coefficient of temperature-dependent hydrogen diffusion within YHx, the Soret coefficient. This parameter is essential for determining the safe operating modes of YHx moderators. Zirconium Hydride (ZrHx) is used in the Training, Research, Isotopes, General Atomics (TRIGA) reactor, is the reference material for these calculations. In this work, nanoscale atomic modeling in the Vienna Ab initio Simulation Package (VASP) is combined with the mesoscale finite element phase-field module in the Multiphysics Object-Oriented Simulation Environment (MOOSE); this culminates in a new multiscale computational method to simulate Soret diffusion of hydrogen in YHx. This data is useful for predicting experimental outcomes. This workflow involves convergence testing followed by static, Nudged Elastic Band (NEB), Quasi-Harmonic Approximation (QHA), and Molecular Dynamics (MD) calculations - linked with phase field simulation. NEB simulates hydrogen migration, while QHA and MD predict temperature-dependent properties. The static calculations align with literature, and preliminary NEB and QHA calculations yield accurate results. Once the atomic calculations are complete, we will incorporate Electron Backscatter Diffraction (EBSD) images and VASP-generated parameters into the phase field module to simulate intra- and intergranular transport of hydrogen in ZrHx and YHx. Future research will extend our approach to fuel-moderator materials systems such as Uranium-Yttrium Hydride (U-YHx). This work contributes to the development of advanced nuclear energy solutions for space travel.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

Reduced Order Model for Guided Wave Propagation on Gas Pipelines to Enable Real-Time Simulation

Reduced order model for simulation of Guided wave propagation is presented here. The utilization of reduced order models ensures efficient data generation for a variety of parameters where it takes huge computational effort to simulate, crucial for timely monitoring and decision-making. Autoencoder based reduced order models are proposed here, which are trained on simulated data from open-source finite element framework, Firedrake.

Bukka, Sandeep Reddy↗