Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element simulation”

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.

333 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2026

Development and Experimental Validation of a Path-Dependent Spin Forming Finite Element Model

Spin forming is an advanced manufacturing process widely used in the aerospace and defense sectors to produce lightweight, high-strength cylindrical components with tight dimensional tolerances. This study explores the applicability of the path-dependent Mechanical Threshold Stress (MTS) constitutive model by simulating the evolution of geometry, machining forces, and plastic deformation during the spin forming of a 10-mm thick 6061-O aluminum cylinder. While numerical modeling of spin forming has advanced substantially over the past decade, systematic verification and experimental validation of material models remain limited, particularly in predicting through-thickness process evolution. The MTS model, incorporating a Voce hardening rule, is employed for its ability to represent cyclic loading, rapidly varying temperature fields, and strain rates characteristic of spin forming. Numerical convergence analysis indicates discretization uncertainties between 0.3% and 9.2% for key quantities of interest. Experimental validation demonstrates that the MTS model, when implemented with a verified mesh, accurately reproduces both elastic and plastic behavior of 6061-O aluminum, predicting peak roller loads within 11–18% of measurements, geometric tolerances within 3%, and plastic strain distributions within 10% of experimental values. Collectively, these results establish a validated computational framework for predictive spin-forming simulations with quantified confidence, providing a foundation for extension to other alloys, geometries, and forming conditions.

Spin forming↗
NASA NTRS2011

Thermal Cycling and Isothermal Deformation Response of Polycrystalline NiTi: Simulations vs. Experiment

A recent microstructure-based FEM model that couples crystal-based plasticity, the B2<-> MB190 phase transformation and anisotropic elasticity at the grain scale is calibrated to recent data for polycrystalline NiTi (49.9 at.% Ni). Inputs include anisotropic elastic properties, texture and differential scanning calorimetry data, as well as a subset of recent isothermal deformation and load-biased thermal cycling data. The model is assessed against additional experimental data. Several experimental trends are captured - in particular, the transformation strain during thermal cycling monotonically increases and reaches a peak with increasing bias stress. This is achieved, in part, by modifying the martensite hardening matrix proposed by Patoor et al. [Patoor E, Eberhardt A, Berveiller M. J Phys IV 1996;6:277]. Some experimental trends are underestimated - in particular, the ratcheting of macrostrain during thermal cycling. This may reflect a model limitation that transformation-plasticity coupling is captured on a coarse (grain) scale but not on a fine (martensitic plate) scale.

Phase Transformations↗
NASA NTRS1995

A Probabilistic Model of a Porous Heat Exchanger

This paper presents a probabilistic one-dimensional finite element model for heat transfer processes in porous heat exchangers. The Galerkin approach is used to develop the finite element matrices. Some of the submatrices are asymmetric due to the presence of the flow term. The Neumann expansion is used to write the temperature distribution as a series of random variables, and the expectation operator is applied to obtain the mean and deviation statistics. To demonstrate the feasibility of the formulation, a one-dimensional model of heat transfer phenomenon in superfluid flow through a porous media is considered. Results of this formulation agree well with the Monte-Carlo simulations and the analytical solutions. Although the numerical experiments are confined to parametric random variables, a formulation is presented to account for the random spatial variations.

O P Agrawal↗
NASA NTRSDate not supplied

Turbomachinery Simulation Impact on Design, Understanding, and Optimization

This presentation shows the impact of Turbomachinery Simulation from simple analytical simulation to high fidelity CFD and Finite Element Analysis on the design of turbomachinery and the understanding of flow physics that is then used to improve design approaches. The impact of Optimization is also presented. The best approach for the tool development is to work with a compressor, fan, or turbine designer or to work on the design process directly. The presentation represents the work and impact of the author over his 45-year career and provides insight for both new and experienced engineers. The presentation explores applications of distortion from a downstream fan frame, the first uses of 3D CFD for fan, compressor and turbine design, and approaches to optimization for performance and structures.

optimization↗
NASA NTRSDate not supplied

Multidisciplinary Design Optimization and Analysis of an Open Rotor Stage: Part 1

Successful design of open rotor propulsors requires effective analysis across multiple disciplines, including aerodynamics, acoustics, and structures. A viable design must not only be efficient but must also produce an acceptable level of noise and meet all static and dynamic structural requirements. For design and optimization, this is especially challenging because running high fidelity analyses is resource-intensive, and optimizing a design may require many thousands of cases to be analyzed. For this reason, the NASA team has applied design methodology that utilizes low-cost aerodynamic methods, machine learning models, and high-fidelity analyses when necessary. This includes standard two-dimensional methods such as throughflow analysis and quasi-3D blade-to-blade CFD, as well as some newly developed methods. Optimization using 3D CFD is necessary to maximize performance, and this is considered as well. All optimizations are carried out subject to structural constraints evaluated using finite element analysis. Doing this accurately requires a robust trunnion design, capable of pitching the blade between cruise and takeoff conditions while maintaining acceptable factor of safety. Hot to cold analysis must also be applied in order to correctly determine the as-manufactured shape. For acoustics, the unsteady pressures on the blade surfaces must be predicted, and this can be done either through full-annulus unsteady CFD or through a nonlinear harmonic method (NLH). NLH can run much faster, allowing some acoustic considerations to be made earlier in the design process. The design process is ongoing, and this presentation will review the current status and planned next steps. This part of the talk will focus on aerodynamic performance and be followed by a talk on structures and acoustics.

Design↗
NASA NTRSDate not supplied

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions to the finite element problem, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization↗
NASA NTRSDate not supplied

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. The transient finite-element analysis is done using TACS, and the cell-discharge battery model uses OpenMDAO and dymos. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization↗
NASA NTRSDate not supplied

High-Fidelity Modeling Approach for Coupling Thermal, Fluid, and Neutronic Analysis in an NTP Engine

Nuclear thermal propulsion (NTP), with its increased efficiency compared to traditional rocket engines, is a potentially enabling technology for space missions of higher complexity, such as delivering heavy machinery or crew to the moon or other planets. Thus, exploring modeling approaches that can effectively couple the thermal, fluid, and neutronic behaviors of NTP reactors is essential to mitigate cost and schedule challenges with hardware development. This paper details one such approach which utilizes Cardinal to unify reactor physics analysis performed using OpenMC, heat transfer analysis performed using the MOOSE framework, and thermal hydraulic analysis performed using the MOOSE Thermal Hydraulic Module. This approach is demonstrated with a proof of concept that focuses on a legacy reactor geometry from the historic NERVA program. This preliminary model found an overprediction in, when compared to the reference data, axial power peaking and the temperatures of both solids and fluids by roughly 20%. It is clear, however, that the model conserved a number of broad trends, and a few model refinements should reduce the observed error. This approach has definite merit and should enable the production of informative results for more modern NTP systems. This work is a first step in creating a transient reactor model that can be updated across a sequence of time steps based on continually fluctuating, time-dependent boundary conditions.

Nuclear Thermal Propulsion↗
NASA NTRSDate not supplied

Multiscale Modeling of Fracture Strength in Fibrous Thermal Protection System Materials

This work presents a multiscale modeling approach to predict the fracture strength of fibrous Thermal Protection System (TPS) materials. The model assumes that system failure is initiated at the joints between individual fibers. We investigated three distinct TPS compositions: amorphous silica, alumina and aluminosilicate fibers. Molecular dynamics (MD) simulations were employed to determine the fracture strength values of these fiber joints for both material systems. These fracture strength values were then integrated into simulations of 3D randomly populated fiber structures, where tensile load transfer occurs through the fiber joints. These microscale properties are upscaled through a renormalization approach [1] to predict macroscale tensile strength of 3D random fiber networks, accounting for joint-dominated failure and effective load-bearing area. The study concludes by demonstrating the resulting strength variation as a function of material composition, fiber density, and morphology. We also show validation of results by comparing them against explicit fiber finite element (FE) modeling [2] where fiber joint fracture is represented by cohesive elements.

Jaehyun Cho↗
NASA NTRSDate not supplied

The Use of Immersion Rigs for High Temperature Hydrogen Exposure Testing Within the Nuclear Thermal Rocket Element Environmental Simulator (NTREES): NTREES-5 Hydrogen Heat Exchanger (N5HHX-01)

The Nuclear Thermal Rocket Element Environmental Simulator (NTREES) facility at NASA’s Marshall Space Flight Center (MSFC) subjects materials and components typically found on nuclear thermal propulsion (NTP) systems to non-nuclear simulations of prototypic NTP reactor thermohydraulic environments. With 1.2 MW of inductive power available, the currently employed methods to simulate nuclear heating in NTREES consist of 1) inductively heating a test specimen, which typically cannot be performed on ceramic materials, or 2) heating gas that envelops the test specimen. The latter of these has provided a basis for developing so-called “immersion rigs” to be used in NTREES. Furthermore, immersion rigs are especially useful in the case of smaller, technically challenging test specimens. One such device, the “NTREES-5 Hydrogen Heat Exchanger” (N5HHX-01), designed and deployed to test small material coupons to ~0.5 g/s of 2800 K, 500 psi Hydrogen, is discussed herein. To achieve this desired delivery temperature level, the most suitable material for the N5HHX’s construction was pure tungsten, which imposed budgetary and manufacturing challenges. Also discussed in this paper is the design of the N5HHX rig’s sample container, which became the first iteration of improving NTREES small material coupon test operations by permitting rapid changeout of different material samples.

Reactor Testing↗
NASA NTRS1995

Numerical Solution of Fluid Flow and Heat Tranfer Problems with Surface Radiation

This paper presents a numerical scheme, based on the finite element method, to solve strongly coupled fluid flow and heat transfer problems. The surface radiation effect for gray, diffuse and isothermal surfaces is considered. A procedure for obtaining the view factors between the radiating surfaces is discussed. The overall solution strategy is verified by comparing the available results with those obtained using this approach. An analysis of a thermosyphon is undertaken and the effect of considering the surface radiation is clearly explained.

S Ahuja↗
NASA NTRSDate not supplied

The Use of Immersion Rigs for High Temperature Hydrogen Exposure Testing within the Nuclear Thermal Rocket Element Environmental Simulator (NTREES): Thermal Soak Rig (TSR)

The Nuclear Thermal Rocket Element Environmental Simulator (NTREES) facility was purpose constructed to perform non-nuclear evaluations of nuclear thermal propulsion (NTP) system fuel materials and structures within prototypic thermochemical environments. This system has been utilized steadily in its ability to subject test specimens to thermochemical and thermohydraulic environments simulating that of an operating nuclear rocket engine. Fission heat is simulated by induction power and experiments are conducted within a ~1000 psi pressure vessel. Hydrogen is conventionally passed through the heated fuel surrogate test specimen while pressure, temperature, and gas species data are collected at various points along the experiment. In order to test fuel and material coupon samples, a class of test apparatus named “immersion rigs” are being developed and employed to more rapidly test these smaller and more technically challenging test specimen. One example of a promising potential fuel structure, Tristructural-isotropic (TRISO) particles, presents unique challenges for testing of this type. TRISO fuel micro-particles are spheroids typically on the order of 500 – 1000 μm in diameter, and exposing a batch sample to hot hydrogen requires purpose-built special test equipment. Thusly, an immersion rig was developed and successfully demonstrated to expose ~1 g of micro-particles to hydrogen gas at temperatures and pressures relevant to NTP systems for the purpose of fuel evaluation. The rig, comprised primarily of graphite and pure tungsten, houses in its core a batch of micro-particles between pucks of porous silicon carbide (SiC). This approach permits gas flow while simultaneously retaining the particles in place. Herein is a discussion of the design, analysis, fabrication, and testing of the NTREES Thermal Soak Rig (TSR).

Space Nuclear Propulsion↗
NASA NTRSDate not supplied

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Discrete Element Method↗
NASA NTRSDate not supplied

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials (3MDCP)

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Woven↗
NASA NTRSDate not supplied

Thermomechanical Modeling of Woven Materials With Particle-Based, Explicit-Fiber Simulations

Fiber-based materials are extensively used to protect spacecraft during entry. Insulative fibers, often in a fiber network or woven, provide rigidity, strength, and control of material anisotropy and density. Woven thermal protection materials, such as ADEPT (Adaptable, Deployable Entry and Placement Technology), 3D-MAT (3-Dimensional Multifunctional Ablative Thermal Protection), and 3MDCP (3D Woven Mid-Density Carbon Phenolic), enable missions with stronger and denser materials for entry profiles with high shear and heat flux. Vulnerabilities to woven thermal protection materials include manufacturing-induced material property variation, and impact from micrometeoroids. Simulating woven materials under these conditions require models that can resolve hierarchal structures, thermomechanical behavior, and failure. To address this, we simulate weave thermal conduction and mechanical deformation. We simulate the full weave with a coarse-grained yarn model is presented. The model combines a validated, high-resolution single 3MDCP yarn model and phenolic resin model. Instead of modeling every fiber, each yarn ply with order 10, instead of order 1000, fibers. The discrete element bonded particle model (DEM-BPM) of fibers captures the thermal and mechanical behavior within and between fibers. We study the proportion of heat transfer and stress via the contact network, fiber bonds, and overall weave geometry.

bonded particle↗
NASA NTRS1995

The Benard Problem: A Comparison of Finite Difference and Spectral Collocation Eigen Value Solutions

The application of spectral methods, using a Chebyshev collocation scheme, to solve hydrodynamic stability problems is demonstrated on the Benard problem. Implementation of the Chebyshev collocation formulation is described. The performance of the spectral scheme is compared with that of a 2 nd order finite difference scheme. An exact solution to the Marangoni-Benard problem is used to evaluate the performance of both schemes. The error of the spectral scheme is at least seven orders of magnitude smaller than finite difference error for a grid resolution of N = 15 (number of points used). The performance of the spectral formulation far exceeded the performance of the finite difference formulation for this problem. The spectral scheme required only slightly more effort to set up than the 2 nd order finite difference scheme. This suggests that the spectral scheme may actually be faster to implement than higher order finite difference schemes.

J Raymond Lee Skarda↗