Search NASA⌕ Search

SEARCH · Search NASA

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.

At least 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

AL-CRADA-2025-01 Final Report: American-Made Wind Turbine Materials Recycling Prize Accelerate! Phase 2

This report summarizes the Phase 2 achievements of the Wind Turbine Materials Recycling Prize project led by Critical Materials Recycling (CMR) in partnership with Ames National Laboratory. The team demonstrated the technical feasibility and economic viability of recovering and reusing rare-earth Nd-Fe-B magnets from end-of-life (EOL) wind turbine generators. Through mechanical disassembly, demagnetization, and precision processing, magnets were harvested and characterized to assess structural, compositional, and magnetic properties. Results showed that recovered magnets performed on par with commercial-grade counterparts, enabling their integration into actuator and motor prototypes without redesign. System-level modeling and finite-element simulations validated operational performance, while Life Cycle Cost (LCC) analysis revealed up to 20% material cost savings and 2–3% reductions in overall motor system costs. Life Cycle Assessment (LCA) further confirmed substantial environmental benefits, including a ~94% reduction in global warming potential relative to virgin magnet production. The project’s success was enabled by a tightly integrated partnership between industry and national lab experts, illustrating a viable path for domestic magnet-to-magnet recycling. These findings support a scalable, circular solution for rare-earth material recovery, reducing dependence on imported critical minerals and advancing sustainable clean energy technologies.

Lograsso, Thomas [Ames Laboratory (AMES), Ames, IA↗
DOE OSTI2023

A perfused multi-well bioreactor platform to assess tumor organoid response to a chemotherapeutic gradient

There is an urgent need to develop new therapies for colorectal cancer that has metastasized to the liver and, more fundamentally, to develop improved preclinical platforms of colorectal cancer liver metastases (CRCLM) to screen therapies for efficacy. To this end, we developed a multi-well perfusable bioreactor capable of monitoring CRCLM patient-derived organoid response to a chemotherapeutic gradient. CRCLM patient-derived organoids were cultured in the multi-well bioreactor for 7 days and the subsequently established gradient in 5-fluorouracil (5-FU) concentration resulted in a lower IC 50 in the region near the perfusion channel versus the region far from the channel. We compared behaviour of organoids in this platform to two commonly used PDO culture models: organoids in media and organoids in a static (no perfusion) hydrogel. The bioreactor IC 50 values were significantly higher than IC 50 values for organoids cultured in media whereas only the IC 50 for organoids far from the channel were significantly different than organoids cultured in the static hydrogel condition. Using finite element simulations, we showed that the total dose delivered, calculated using area under the curve (AUC) was similar between platforms, however normalized viability was lower for the organoid in media condition than in the static gel and bioreactor. Our results highlight the utility of our multi-well bioreactor for studying organoid response to chemical gradients and demonstrate that comparing drug response across these different platforms is nontrivial.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Nanoindentation Stress Relaxation to Quantify Dislocation Velocity–Stress Exponent

This work reports a new methodology using indentation stress relaxation to characterize the dislocation velocity–stress exponent. Through the indentation stress relaxation process, the dislocation structure builds up at the rate governed by dislocation velocity, which is a function of the externally applied stress. The relationship between the dislocation velocity and stress can thus be derived from the indentation stress relaxation data of the stress as a function of time. In this study, instrumented nanoindentation stress relaxation experiments were performed on pure aluminum samples, following three different initial displacement rates of 100, 400, and 800 nm/s. Based on the scaling properties of dislocation kinetics, the data were interpreted to derive a dislocation velocity–stress exponent of 2.5 ± 0.5 for room-temperature aluminum. Crystal plasticity finite-element simulations were performed to illustrate the sensitivity of the proposed nanoindentation stress relaxation methodology to the dislocation velocity–stress exponent value.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

In situ imaging during compression of plastic bonded explosives for damage modeling

Here, the microstructure of plastic bonded explosives (PBXs) is known to influence behavior during mechanical deformation, but characterizing the microstructure can be challenging. For example, the explosive crystals and binder in formulations such as PBX 9501 do not have sufficient X-ray contrast to obtain three-dimensional data by in situ, absorption contrast imaging. To address this difficulty, we have formulated a series of PBXs using octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) crystals and low-density binder systems. The binders were hydroxyl-terminated polybutadiene (HTPB) or glycidyl azide polymer (GAP) cured with a commercial blend of acrylic monomers/oligomers. The binder density is approximately half of the HMX, allowing for excellent contrast using in situ X-ray computed tomography (CT) imaging. The samples were imaged during unaxial compression using micro-scale CT in an interrupted in situ modality. The rigidity of the binder was observed to significantly influence fracture, crystal-binder delamination, and flow. Additionally, 2D slices from the segmented 3D images were meshed for finite element simulation of the mesoscale response. At low stiffness, the binder and crystal do not delaminate and the crystals move with the material flow; at high stiffness, marked delamination is noted between the crystals and the binder, leading to very different mechanical properties. Initial model results exhibit qualitatively similar delamination.

36 MATERIALS SCIENCE↗
DOE OSTI2018

FEM Modeling of In-Plane Stress Distribution in Thick Brittle Coatings/Films on Ductile Substrates Subjected to Tensile Stress to Determine Interfacial Strength

The ceramic-metal interface is present in various material structures and devices that are vulnerable to failures, like cracking, which are typically due to their incompatible properties, e.g., thermal expansion mismatch. In failure of these multilayer systems, interfacial shear strength is a good measure of the robustness of interfaces, especially for planar films. There is a widely-used shear lag model and method by Agrawal and Raj to analyse and measure the interfacial shear strength of thin brittle film on ductile substrates. The use of this classical model for a type of polymer derived ceramic coatings (thickness ~18 μm) on steel substrate leads to high values of interfacial shear strength. Here, we present finite element simulations for such a coating system when it is subjected to in-plane tension. Results show that the in-plane stresses in the coating are non-uniform, i.e., varying across the thickness of the film. Therefore, they do not meet one of the basic assumptions of the classical model: uniform in-plane stress. Furthermore, effects of three significant parameters, film thickness, crack spacing, and Young’s modulus, on the in-plane stress distribution have also been investigated. ‘Thickness-averaged In-plane Stress’ (TIS), a new failure criterion, is proposed for estimating the interfacial shear strength, which leads to a more realistic estimation of the tensile strength and interfacial shear strength of thick brittle films/coatings on ductile substrates.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Study of the Printability, Microstructures, and Mechanical Performances of Laser Powder Bed Fusion Built Haynes 230

The nickel-based superalloy, Haynes 230 (H230), is widely used in high-temperature applications, e.g., heat exchangers, because of its excellent high-temperature mechanical properties and corrosion resistance. As of today, H230 is not yet in common use for 3D printing, i.e., metal additive manufacturing (AM), primarily because of its hot cracking tendency under fast solidification. The ability to additively fabricate components in H230 attracts many applications that require the additional advantages leveraged by adopting AM, e.g., higher design complexity and faster prototyping. In this study, we fabricated nearly fully dense H230 in a laser powder bed fusion (L-PBF) process through parameter optimization. The efforts revealed the optimal process space which could guide future fabrication of H230 in various metal powder bed fusion processes. The metallurgical analysis identified the cracking problem, which was resolved by increasing the pre-heat temperature from 80 °C to 200 °C. A finite element simulation suggested that the pre-heat temperature has limited impacts on the maximum stress experienced by each location during solidification. Additionally, the crack morphology and the microstructural features imply that solidification and liquation cracking are the more probable mechanisms. Both the room temperature tensile test and the creep tests under two conditions, (a) 760 °C and 100 MPa and (b) 816 °C and 121 MPa, confirmed that the AM H230 has properties comparable to its wrought counterpart. The fractography showed that the heat treatment (anneal at 1200 °C for 2 h, followed by water quench) balances the strength and the ductility, while the printing defects did not appreciably accelerate part failure.

14 SOLAR ENERGY↗
DOE OSTI2017

Crystal Plasticity Modeling and Experimental Validation with an Orientation Distribution Function for Ti-7Al Alloy

An orientation distribution function based model is used for micromechanical modeling of the titanium-aluminum alloys, Ti-0 wt % Al and Ti-7 wt % Al, which are in demand for many aerospace applications. This probability descriptor based modeling approach is different than crystal plasticity finite element techniques since it computes the averaged material properties using upper bound averaging. A rate-independent single-crystal plasticity model is implemented to compute the effect of macroscopic strain on the polycrystal. An optimization problem is defined for calibrating the basal, prismatic, pyramidal slip system and twin parameters using the available tension and compression experimental data. The crystal plasticity parameters of Ti-7 wt % Al are not studied extensively in literature, and therefore the optimization results for the crystal plasticity model realization produce unique data, which will be beneficial to future studies in the field. The sensitivities of the slip and twin parameters to the design objectives are also investigated to identify the most critical slip system parameters. Using the optimum design parameters, the microstructural textures, during the tension test, are predicted by the crystal plasticity finite element simulations, and compared to the available experimental texture and scanning electron microscope—digital image correlation data.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Investigation of Goubau Wave Propagation on Large Pipes for Sensing Applications

We examine the application of guided waves on a single conductor (Goubau waves) for sensing. In particular, the use of such waves to remotely interrogate surface acoustic wave (SAW) sensors mounted on large-radius conductors (pipes) is considered. Experimental results using a small-radius (0.0032 m) conductor at 435 MHz are reported. The applicability of published theory to conductors of large radius is examined. Finite element simulations are then used to study the propagation and launching of Goubau waves on steel conductors up to 0.254 m in radius. Simulations show that waves can be launched and received, although energy loss into radiating waves is a problem with current launcher designs.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Investigating the Electromechanical Sensitivity of Carbon-Nanotube-Coated Microfibers

The piezoresistance of carbon nanotube (CNT)-coated microfibers is examined using diametric compression. Diverse CNT forest morphologies were studied by changing the CNT length, diameter, and areal density via synthesis time and fiber surface treatment prior to CNT synthesis. Large-diameter (30–60 nm) and relatively low-density CNTs were synthesized on as-received glass fibers. Small-diameter (5–30 nm) and-high density CNTs were synthesized on glass fibers coated with 10 nm of alumina. The CNT length was controlled by adjusting synthesis time. Electromechanical compression was performed by measuring the electrical resistance in the axial direction during diametric compression. Gauge factors exceeding three were measured for small-diameter (<25 μm) coated fibers, corresponding to as much as 35% resistance change per micrometer of compression. The gauge factor for high-density, small-diameter CNT forests was generally greater than those for low-density, large-diameter forests. A finite element simulation shows that the piezoresistive response originates from both the contact resistance and intrinsic resistance of the forest itself. The change in contact and intrinsic resistance are balanced for relatively short CNT forests, while the response is dominated by CNT electrode contact resistance for taller CNT forests. These results are expected to guide the design of piezoresistive flow and tactile sensors.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Limiting Rotation Rate of Neutron Stars from Crust Breaking and Gravitational Waves

Abstract Neutron stars are not observed to spin faster than about half their breakup rate. This limiting rotational frequency may be related to the strength of their crusts. As a star spins up from accretion, centrifugal forces stress the crust. We perform finite-element simulations of rotating neutron stars and find that the crust fails at rotation rates about half the breakup rate. Given uncertainties in microphysics, we have not determined the crust configuration after this failure. Instead, we argue that the crust may fail in an asymmetric way and could produce a configuration with a significant ellipticity (fractional difference in moments of inertia). If the ellipticity is large, a rotating star will radiate gravitational waves that may limit further spin-up. These stars may be promising sources for LIGO/VIRGO and next-generation gravitational-wave detectors.

Morales, J. A.↗
DOE OSTI2017

The Role of Second Phase Hard Particles on Hole Stretchability of two AA6xxx Alloys

The hole stretchability of two Aluminum Alloys (AA6111 and AA6022) are studied by using a two stages integrated finite element framework where the edge geometry and edge damages from the hole piercing processes were considered in the subsequent hole expansion processes. Experimentally it has been found that AA6022 has higher hole expansion ratios than those of AA6111. This observation has been nicely captured by finite element simulations. The main cause of differences have been identified to the volume fractions of the random distributed second phase hard particles which play a critical role in determining the fracture strains of the materials.

Hu, Xiaohua↗
DOE OSTI2022

Local Thermomechanical Processing for Improving Formability of High Strength Aluminum Sheets

Limited room temperature formability hinders the wide-spread use of high strength aluminum alloys in structural body-in-white parts. Stamping or extrusion at warm temperatures or from softer tempers are the current solutions. In this work, our approach is to start with age-hardened sheets from 7xxx and 6xxx family of alloys and improve their formability using local thermomechanical processing only in the regions demanding highest ductility in the forming processes. We achieved local formability improvements with friction stir processing and introduce another process named roller-bending as a concept and showed its feasibility through finite element simulations. Initial results from FSP indicated significant deformation in the processed zones with minimal sheet distortion. FSP also resulted in dynamically recrystallized, fine grained (d < 5 µm) microstructures in the processed regions with textures significantly different from the base material. This resulted in formability improvements of > 30% in 7085-T76 and > 50% in 6111-T6 alloys as measured by minimum bend radius (r/t) ratios established with V-bending tests performed at 90° and 180° bend angles. Based on the encouraging initial results, we envision applicability of these processes in both pre- and post-stamping stages through integration of the process tools as end-effectors to robotic arms.

Nasim, Wahaz↗
DOE OSTI2022

pyvisco [SWR-22-30]

pyvisco is a Python library that supports the identification of Prony series parameters for linear viscoelastic materials described by a Generalized Maxwell model. The necessary material model parameters are identified by fitting a Prony series to the experimental measurement data. pyvisco allows for the identification of Prony series parameters from experimental data measured in either the frequency-domain (via Dynamic Mechanical Thermal Analysis) or time-domain (via relaxation measurements). The experimental data can be provided as raw measurement sets at different temperatures or as pre-processed master curves. An optional minimization routine is included to reduce the number of Prony elements. This routine is helpful in Finite Element simulations where reducing the computational complexity of the linear viscoelastic material models can shorten the simulation time. See also, https://pypi.org/project/pyvisco/

Springer, Martin↗
NASA NTRS1976

Finite-element analysis of dynamic fracture

Applications of the finite element method to the two dimensional elastodynamics of cracked structures are presented. Stress intensity factors are computed for two problems involving stationary cracks. The first serves as a vehicle for discussing lumped-mass and consistent-mass characterizations of inertia. In the second problem, the behavior of a photoelastic dynamic tear test specimen is determined for the time prior to crack propagation. Some results of a finite element simulation of rapid crack propagation in an infinite body are discussed.

Aberson, J. A.↗
NASA NTRS1980

Self-adaptive incremental Newton-Raphson algorithms

Multilevel self-adaptive Newton-Raphson type strategies are developed to improve the solution efficiency of nonlinear finite element simulations of statically loaded structures. The overall strategy involves three basic levels. The first level involves preliminary solution tunneling via primative operators. Secondly, the solution is constantly monitored via quality/convergence/nonlinearity tests. Lastly, the third level involves self-adaptive algorithmic update procedures aimed at improving the convergence characteristics of the Newton-Raphson strategy. Numerical experiments are included to illustrate the results of the procedure.

Padovan, J.↗
NASA NTRS1984

On the development of hierarchical solution strategies for nonlinear finite element formulations

This paper develops a hierarchical type solution scheme which can handle the field equations associated with nonlinear finite element simulations. The overall procedure possesses various levels of application namely degree of freedom, nodal, elemental, substructural as well as global. In particular iteration, updating, assembly and solution control occurs at the various hierarchical levels. Due to the manner of formulation, the degree of matrix inversion depends on the size of the various hierarchical partitioned groups. In this context, degree of freedom partitioning requires no inversion. To benchmark the overall scheme, the results of several numerical examples are presented.

Padovan, J.↗
NASA NTRS1985

Orbit on demand - Structural analysis finds vertical launchers weigh less

Structural considerations arising from favored design concepts for the next generation on-demand launch vehicles are explored. The two emerging concepts are a two stage fully reusable vertical take-off vehicle (V-2) and a horizontal take-off, two stage subsonic boost launch vehicle (H-2-Sub). Both designs have an 1100 n. mi. cross-range capability, with the V-2 orbiter having small wings with winglets for hypersonic trim and the H-2-Sub requiring larger, swept wings. The rockets would be cryogenic, while airbreathing initial boosters would be either turbofans, turbojets and/or ramjets. Dynamic loading is lower in the launch of a V-2. The TPS is a critical factor due to thinner leading edges than on the Shuttle and may require heat-pipe cooling. Airframe structures made of metal matrix composites have passed finite element simulations of projected loads and can now undergo proof-of-concept tests, although whisker-reinforced materials may be superior once long-whisker technology is developed.

Taylor, A. H.↗