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.

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Universal Design Curves for Elastomeric Direct Ink-Writing Architectures

Direct ink writing of silicone elastomers enables printing with precise control of porosity and mechan- ical properties of ordered cellular solids, suitable for shock absorption and stress mitigation applications. With the ability to manipulate structure and feedstock stiffness, the design space becomes challenging to parse to obtain a solution producing a desired mechanical response. Here we derive an analytical design approach for a specific architecture. Results from finite element simulations and quasi-static mechanical tests of two different parallel strand architectures were analyzed to understand the structure-property relationships under compressive strain. Combining effective stiffness-density scaling with least squares optimization of the stress responses yielded general response curves parameterized by resin modulus and strand spacing. An analytical expression of these curves serves as a reduced order model, which, when optimized, provides a rapid design capability for filament-based 3D printed structures. As a demonstration, the optimal design of a face-centered tetragonal architecture is computed that satisfies prescribed minimum and maximum load constraints.

36 MATERIALS SCIENCE↗
DOE OSTI2018

NDT System for Metal Components in Concrete Structures (Final Report)

Steel reinforcements embedded within concrete can corrode and crack, leading to weakening of the structures, and existing NDT (non-destructive testing) techniques do not provide the full picture. Condition assessment is critical in evaluating these structures in order to ensure safe, long-term operation of the nuclear power plants. Improved NDT techniques are needed as the nation’s generating stations age and require action based on increasingly critical inspection data. The power industry uses eddy current test (ECT) during in-service inspections to detect and monitor defects such as cracks and corrosion in metals, particularly metal pipes. However, the detection of defects in the reinforcing steel embedded in concrete is far more challenging for numerous reasons. Primary among them is the concrete layer’s thickness, which prevents the ability to have a close proximity between sensor and metallic sample. One solution to this problem is to create a higher sensitivity technique which can inspect at depth, with good spatial resolution and excellent directional sensitivity, to isolate the metal buried deep in the concrete and characterize defects in the metal. The goal of the phase II research program was to study and develop a new reinforced concrete NDT technology based on anisotropic magnetoresistive (AMR) sensors designed specifically for ECT of steel reinforcements in concrete. The first technical objective was to demonstrate the capability of inspecting steel reinforcement structures, including metal rebar and tensioning cable at or below 6” of depth embedded in concrete. Another objective was to be able to confidently identify multiple layers of reinforcements at varying depths. In addition, the technique needed to be able to provide condition assessment of the metal, providing data on the structure and positioning of the steel members, material loss, and detection of the presence of breaks. An important part of the research and development was aimed at achieving superior performance compared to existing technology. The ability to detect corrosion breaks in unbonded multi-wire and multi-strand tendons was another important objective. The final objective that led to a successful program was to build a complete prototype portable scanning system capable of demonstrating the probe performance in the field. All of these objectives were met, which provided data for the system development work going forward. RMD worked with the engineering firm Simpson Gumpertz & Heger to develop a detailed problem definition for concrete structures and the necessary sensor specifications. In addition, we developed finite element simulation models for the eddy current based NDT system through the analysis of the electromagnetic (EM) properties of sample materials and the orientation of the transmitter, sensor array (receivers), metal, and intervening materials like concrete and air. Based on the modeling and specifications, we designed and fabricated prototype sensor modules. We then acquired and fabricated test samples to characterize and optimize the sensors. Finally, we built a complete prototype for testing on actual concrete structures. This new NDT technology has high promise for greatly improving capabilities in the inspection of steel reinforcements within concrete structures, such as rebar, liners, anchors, and prestressing cables. This method could be valuable in the maintenance of nuclear power plants, and many other structures critical to our nation’s infrastructure.

42 ENGINEERING↗
DOE OSTI2018

Thermomechanical Design Optimization of Roll-to-Roll Manufactured Hybrid Metal Polymer Heat Exchanger

Metal based heat exchangers for ultra-low temperature (<100C) waste heat recovery (WHR) face high maintenance and material costs that surpass profits from recovery[1]. Low-cost polymer pipes suffer from poor thermal conductivity (~0.2 W/mK), and thus a low overall heat transfer coefficient. Thermal conductivity of polymers can be enhanced through systematic molecular alignment by stretching[2], chemical vapor deposition of extended polymers [3], surface grating [4], etc. However, these approaches suffer from scalability issues for WHR heat exchangers in terms of cost [3], and thermomechanical considerations [5]. In this work, we propose novel hybrid metal-polymer strips that are strategically designed to enhance thermal conductivity by cladding copper around polymer strips. The polymer strips with copper-clad edges are wound helically to form the pipes for a heat exchanger. Through finite element simulations, we optimize the design for enhanced thermal performance of the polymer-metal strips. We predict up to 20 % enhancement in overall heat transfer from a base polymer with a thermal conductivity of 0.2 W/mK. We show that the thermomechanical strains induced in the pipes can be reduced to < 0.2 % by strategically stacking and welding multiple polymer-metal strips. We verified the improved thermal performance of the roll-to-roll manufactured strips through a thermal conductivity measurement setup. The optimized design of hybrid metal-polymer pipes provides a scalable and cost-effective means for harvesting waste heat from ultra-low temperature sources. Thekdi, Arvind, and Sachin U. Nimbalkar. Oak Ridge National Lab, 2015 Shen, Sheng, et al. Nature nanotechnology 5.4 (2010): 251. Xu, Yanfei, et al. Science advances 4.3 (2018): 3031. Smith, M. K., Singh, V., Kalaitzidou, K., & Cola, B. A. (2015). ACS nano, 9(2), 1080-1088. Aymonier, C., Bortzmeyer, D., Thomann, R., & Mülhaupt, R. (2003). Chemistry of Materials, 15(25), 4874-4878.

C Rajagopal, Manjunath↗
DOE OSTI2019

Modeling adhesive debonding using the shim layer modeling approach

There are multiple possible methods for modeling adhesive debonding in finite element simulations. However, many of these techniques are incapable of meeting all of the requirements of a thermomechanical simulation involving assemblies with thin layers of elastomeric potting materials. In this paper, these requirements are first identified, then one particular method, termed the “shim layer” modeling approach, is described. This modeling approach is capable of meeting all of the modeling requirements.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An Initial Assessment of the Design Margins of Different ASME Section III, Division 5 Design Rules

This report develops a method for assessing the design margin of the ASME Section III, Division 5 Class A design rules, focusing on a definition of margin as the ratio between the actual, expected component life and the ASME design life. Full inelastic finite element simulations with a complete damage model capturing the available creep, creep-fatigue, fatigue, and tension test failure data produce the expected component lives, which can then be compared to corresponding ASME design calculations. The focus of this particular work is on Alloy 617 and the ASME creep-fatigue design rules, but the approach is general and could be applied to other design limits and to other materials. A margin assessment of a representative Alloy 617 component shows the ASME creep-fatigue design rules applied to this component have a margin of 10 on design life. The report describes several challenges in completing a more comprehensive margin assessment of the ASME rules, in particular challenges in developing probabilistic assessment methods for high temperature structural components. The report details these challenge and discusses the immediate possibilities of the deterministic margin assessment method developed here.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Fracture Toughness of Microstructural Gradients

Traditional singularity-based fracture mechanics theories rely on their ability to infer the crack tip driving force (local field) by surveying macroscopic physical magnitudes far from the crack tip (far field). This key capability allows engineers to employ nominal forces or displacements to estimate the potential for stable or unstable crack growth. In the case of heterogeneous or anisotropic materials, traditional fracture approaches are not fully theoretically sound and applications rely on extrapolating methodologies with ad-hoc corrections. This Express Laboratory Directed Research and Development (ELDRD) program employed mesoscale-sensitive finite element simulations to assess the impact of grain size and texture on the crack tip behavior. A dislocation-based crystal plasticity model conveys grain size effects by computing the constraint on dislocation cell structures. We assessed the effects of microstructural variability on multiple displacement-based measurements of the fracture driving forces for crack opening (Mode I) and sliding (Mode II). We also consider multiple microstructural realizations of single phase metals undergoing ductile failure. The results show that grain size and texture affect the applied fracture driving force and can induce a significant Mode II deformation under force and displacement control, which is completely neglected in homogeneous models. A large variability in driving forces upon identical far field applied conditions is attributed to a buffering effects of the microstructure. Furthermore, crack mouth opening displacement is almost insensitive to microstructure, which suggests that experimental measurements using such a magnitude (e.g., plastic hinge model) may underestimate local crack tip driving force variability.

42 ENGINEERING↗
DOE OSTI2019

Mechanical Response of Additively Manufactured Stainless Steel 304L Across a Wide Range of Loading Conditions

The mechanical response of additively manufactured (AM) stainless steel 304L has been investigated across a broad range of loading conditions, covering 11 decades of strain rate, and compared with the behaviors of traditional ingot-derived (wrought) material. In general, the AM material exhibits a greater strength and reduced ductility compared with the baseline wrought form. These differences are consistently found from quasi-static and high strain rate tests. A detailed investigation of the microstructure, the defect structure, the phase, and the composition of both forms reveals differences that may contribute to the differing mechanical behaviors. Compared with the baseline wrought material, dense AM stainless steel 304L has a more complex grain structure with substantial sub-structure, a fine dispersion of ferrite, increased dislocation density, oxide dispersions and larger amounts of nitrogen. In-situ neutron diffraction studies conducted during quasi-static loading suggest that the increased strength of AM material is due to its initially greater dislocation density. The flow strength of both forms is correlated with dislocation density through a square root dependence akin to a Taylor-like relationship. Neutron diffraction measurements of lattice strains also correlate with a crystal plasticity finite element simulations of the tensile test. Other simulations predict a significant degree of elastic and plastic anisotropy due to crystallographic texture. Hopkinson tests at higher strain rates $\dot{ε}$ = 500 and 2500 s -1 ) also show a greater strength for AM stainless steel 304L; although, the differences compared with wrought are reduced at higher strain rates. Gas gun impact tests, including reverse ballistic, forward ballistic and spall tests, consistently reveal a larger dynamic strength in the AM material. The Hugoniot Elastic Limit (HEL) of AM SS 304L exceeds that of wrought material although considerable variability is observed with the AM material. Forward ballistic testing demonstrates spall strengths of AM material (3.27 -- 3.91 GPa) that exceed that of the wrought material (2.63 -- 2.88 GPa). The Hugoniot equation-of-state for AM samples matches archived data for this metal alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Microstructural Model for Creep-Fatigue Interaction in Grade 91 Steel

This report describes the completion of a crystal plasticity finite element simulation framework for creep-fatigue deformation and damage in Grade 91 steel. In addition to a description of the framework and a set of verification calculations, the report details an initial set of simulations applying the framework to examine low cycle fatigue and creep-fatigue interaction in Grade 91 at 600° C. Grade 91 is a likely material for future advanced microreactors and creep-fatigue will be a limiting structural damage mechanism, particularly for mobile reactors or load following systems. The framework described here could provide a better understanding of creep-fatigue damage mechanisms and directly examine creep, fatigue, and creep-fatigue interaction for Grade 91 at realistic component service conditions. This in turn could lead to improved design methods and component life estimation approaches, leading to more efficient and economical component designs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Experimental Study of Mechanical Response of Polyurethane Rubber at Varying Strain Rates for W80-4 System

This report describes an experimental study to determine the mechanical behavior of the polyurethane rubber material that was used in the W80-4 MCC Shock/Breach Phase 1 test series. Compression experiments were conducted on cylindrical specimens over a wide range of loading rates to characterize the material over the range of strain rates that were experienced in the shock/breach testing. Additionally, specimen diameter was varied to determine the effect of confinement on the material response and was found to be significant. This data is used to populate a material model to enable accurate analyses and finite element simulations of the shock/breach test series.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Revisiting Multi-Material Composite Structures with Homogenized Composite Properties

Composite structures inherently develop residual stresses during their curing process. Driven predominately by mismatched thermal strains between differing materials or ply orientations, but also affected by curing process phenomena like polymer shrinkage, these residual stresses can lead to failure within composite structures. There are several methods varying in complexity that can be used to model the development of residual stresses, all of which are capable of capturing sufficient detail to understand the residual stress state at the ply level. However, explicitly modeling all plies of a layup in a composite structure can be prohibitively expensive based on the number of plies, structure size, and required element size. The computational cost can be reduced through the homogenization of the composite layup without losing much fidelity of the overall response of the structure. The homogenization process reduces the many plies of a laminate to a single lamina that reduces complexity and increases the mesh size where a single element can span multiple plies. This report focuses on verification and validation efforts for a homogenization process using a suite of finite element simulations rather than an analytic solution derived from classical laminate theory. Initial verification using representative element volumes indicated there was minimal error in the homogenization process; however, this compounded to a small, but acceptable error in strip and split ring experimental composite structures. The error does under predict the residual stress state in the strip and split ring and should be accounted for when simulating composite structures with homogenized properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Improved Axisymmetric and High Temperature Material Structural Modeling in MOOSE and NEML

This report describes improvements made to the solid mechanics formulation in the MOOSE open source finite element simulation environment and the open source Nuclear Engineering Material model Library (NEML) for mechanical constitutive models. The focus of these improvements is to improve the usability and performance of simulations involving one or both pieces of software. Specifically, this work completes a new system for solid mechanics simulations in the MOOSE ecosystem providing exact linearizations and optimal (quadratic) convergence, for a variety of coordinate systems and material types, including large deformation simulations. This work then provides users a framework to build highly efficient mechanical simulations of structures or materials or to couple in additional MOOSE physics modules to build complex, scalable multiphysics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

White Paper: Scalable Digital Twin Capabilities for Aging and Surveillance of Engineered Systems

This white paper presents a multi-year initiative to develop practical, secure, and scalable digital twin capabilities for engineered systems in aging and surveillance contexts—an approach pioneered at the National Nuclear Security Administration (NNSA) Lawrence Livermore National Laboratory (LLNL) that maps directly onto the needs and ambitions of the Navy for ship- and fleet-level digital twins. LLNL’s work in building part- and process-level digital twins for advanced manufacturing, with a vision to scale up to entire factory floors and, ultimately, enterprise-wide digital twins, offers an adaptable pathway for the Navy as it seeks to modernize lifecycle management, readiness, and predictive maintenance across ships and fleets. For our application, we integrate physics-based modeling with automated data ingestion, processing, and AI-driven calibration, creating hybrid models that are both interpretable and data responsive. We modernized legacy workflows, established centralized data infrastructure, automated experimental pipelines, and demonstrated end-to-end coupling of accelerated aging data with finite element simulations via optimization and surrogate modeling. The result is a generalizable framework that supports part-level digital twins today and lays the groundwork for future system-level twins suitable for Navy applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dynamic Multipoles in the RCS

The Rapid Cycling Synchrotron (RCS) of the Electron--Ion Collider (EIC) will accelerate polarized electrons up to 18 GeV at a repetition rate of 1 Hz. Dynamic multipoles driven by eddy currents in the vacuum chamber can affect beam dynamics and polarization preservation during the energy ramp. Building on the analytical treatment of S.Y. Lee for a ramped normal-conducting dipole magnet, we use conformal transformations to derive compact expressions for dynamic multipoles in higher-order ramped magnets with circular beam pipes. We further derive analytical frequency-domain transfer functions for both self-response and cross-response multipoles, providing compact dynamical models for their evolution during arbitrary ramp waveforms. The resulting analytical predictions are benchmarked against finite-element simulations, including the effects of practical pole geometry and beam-pipe misalignment. Finally, the expected dynamic multipoles in the RCS magnets are estimated and their implications for beam dynamics and polarization preservation are discussed.

43 PARTICLE ACCELERATORS↗
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↗