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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Chemical modelling of Delayed Ettringite Formation for assessment of affected concrete structures

Delayed Ettringite Formation (DEF) modelling is addressed in the context of structural analysis. First, a chemical model is specified to simulate the effects of the heating of concrete in terms of reactions involving sulfates, aluminates and alkalis. It can be decomposed into: dissolution of primary sulfoaluminates at high temperature, fixation of aluminates in hydrogarnet or carboaluminates, and formation of delayed ettringite at low temperature. The influences of alkalis, temperature and water saturation on these reactions are taken into account. The main result of this model is the amount of delayed ettringite. The parameters of the proposed chemical model are fitted on a large number of experiments taken from the literature. The chemical model is then coupled with a nonlinear mechanical model based on poro-mechanics modelling. An application of chemo-mechanical coupling shows the features and benefits of such modelling to help practitioners in the management of structures affected by DEF.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Parameter covariance and non-uniqueness in material model calibration using the Virtual Fields Method

Traditionally, material identification is performed using global load and displacement data from simple boundary-value problems such as uni-axial tensile and simple shear tests. More recently, however, inverse techniques such as the Virtual Fields Method (VFM) that capitalize on heterogeneous, full-field deformation data have gained popularity. In this work, we have written a VFM code in a finite-deformation framework for calibration of a viscoplastic (i.e. strain-rate dependent) material model for 304L stainless steel. Using simulated experimental data generated via finite-element analysis (FEA), we verified our VFM code and compared the identified parameters with the reference parameters input into the FEA. The identified material model parameters had surprisingly large error compared to the reference parameters, which was traced to parameter covariance and the existence of many essentially equivalent parameter sets. This parameter non-uniqueness and its implications for FEA predictions is discussed in detail. Lastly, we present two strategies to reduce parameter covariance – reduced parametrization of the material model and increased richness of the calibration data – which allow for the recovery of a unique solution.

304L stainless steel↗
DOE OSTI2013

A new smoothed particle hydrodynamics non-Newtonian model for friction stir welding: Process modeling and simulation of microstructure evolution in a magnesium alloy

We present a new smoothed particle hydrodynamics (SPH) model for friction stir welding (FSW). FSW has broad commercial application in the marine, aerospace, rail, and automotive industries. However, development of the FSW process for each new application has remained largely empirical. Few established numerical modeling techniques have been developed that can explain and predict important features of the process physics involved in FSW. This is particularly true in the areas of material flow and mixing mechanisms. In this paper, we present a novel modeling approach to simulate FSW that may have significant advantages over current finite element or finite difference based methods. Unlike traditional grid-based methods, Lagrangian particle methods such as SPH can simulate the dynamics of interfaces, large material deformations, and the material’s strain and temperature history without employing complex tracking schemes. Three-dimensional simulations of FSW on AZ31 Mg alloy are performed. The temperature history and distribution, grain size, microhardness as well as the texture evolution are presented. Numerical results are found to be in good agreement with experimental observations.

Pan, Wenxiao↗
DOE OSTI2025

Directed energy deposition of radially functionally graded nickel based superalloys IN625 and IN738: Microstructure and mechanical properties

This study investigates the microstructure and mechanical properties of radially functionally graded (RFG) materials fabricated from nickel-based superalloys IN625 and IN738 using directed energy deposition (DED). RFG cylindrical deposits with radial gradients of IN625-IN738 and IN738-IN625 were successfully fabricated and confirmed to be crack-free through tailored composition grading, scan strategy and process parameters that minimized thermal gradients and solidification velocities. Microstructural analysis revealed distinct variations across the gradient, characterized by differences in dendrite arm spacing, microsegregation, and the formation of distinct γ' precipitates and carbide/boride phases (MC, M 6 C, M 2 C, and MB 2 ) in as-deposited (AD) conditions, which evolved into M 23 C 6 and M 2 B upon subsequent heat treatment (HT). EBSD showed predominant columnar grains with a [100] texture, while microhardness profiles demonstrated consistent variation across the gradients. Mechanical testing highlighted the synergy of the graded structure, combining IN738's strength with IN625's ductility, resulting in enhanced properties such as a yield strength of 1131 MPa, ultimate tensile strength of 1382 MPa, and elongation of 12 % for HT IN625-IN738 grading. Numerical simulations employing CALPHAD thermodynamic calculations and finite element thermal modeling are used to establish phase evolution and microstructural feature correlation with thermal profiles, aiding in the optimization of process-structure-property relationship.

Additive manufacturing↗
DOE OSTI2022

Time-explicit Darwin PIC algorithm

A new approach to Darwin particle-in-cell plasma simulation is described. Using a finite-element approach, the vector potential is assured to be exactly solenoidal. This allows writing the system of multi-species particles as an action-at-a-distance Hamiltonian system. Applying recently-developed explicit symplectic methods to this gives a time-explicit algorithm with all desired properties. Elliptic systems for the electrostatic and magneto-static portions of the problem are inverted efficiently by an algebraic multi-grid algorithm. The algorithm is implemented in a two-dimensional Cartesian-geometry code and tested by application to Weibel instability and to Alfvén wave dynamics. Results show the effectiveness of this approach. Extensions to arbitrary meshes and to a partially time-implicit scheme are briefly discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Effect of stress evolution on microstructural behavior in U-Mo/Al dispersion fuel [Effect of stress on microstructural evolution in U-Mo/Al dispersion fuel]

U-Mo/Al dispersion fuel irradiated to high burnup at high power (high fission rate) exhibited microstructural changes such as deformation of the fuel particles, pore growth, and rupture of the Al matrix. The driving force for these microstructural changes was meat swelling caused by a combination of fuel particle swelling and interaction layer growth. Five miniplates with well-recorded fabrication data and irradiation conditions were selected, and their PIE data was analyzed. ABAQUS finite element analysis (FEA) was utilized to simulate the microstructural evolution of the plates. Using the simulation results shear stress, effective stress and hydrostatic stress exerted on both the fuel particles and the Al matrix were determined. The effects of fabrication and irradiation variables on stress-induced microstructural evolutions, such as pore growth in the interaction layers and Al matrix rupture, were investigated. The observed microstructural changes were consistent with the calculated stress distribution in the meat.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Failure analysis of pinch-torsion tests as a thermal runaway risk evaluation method of Li-Ion Cells

Recently a pinch-torsion test is developed for safety testing of Li-ion batteries (Ren et al., J. Power Source, 2013). It has been demonstrated that this test can generate small internal short-circuit spots in the separator in a controllable and repeatable manner. In the current research, the failure mechanism is examined by numerical simulations and comparisons to experimental observations. Finite element models are developed to evaluate the deformation of the separators under both pure pinch and pinch-torsion loading conditions. It is discovered that the addition of the torsion component significantly increased the maximum principal strain, which is believed to induce the internal short circuit. In addition, the applied load in the pinch-torsion test is significantly less than in the pure pinch test, thus dramatically improving the applicability of this method to ultra-thick batteries which otherwise require heavy load in excess of machine capability. It is further found that the separator failure is achieved in the early stage of torsion (within a few degree of rotation). Effect of coefficient of friction on the maximum principal strain is also examined.

finite element analysis↗
DOE OSTI2019

Towards an integrated experimental and computational framework for large-scale metal additive manufacturing

In this work using the Metal Big Area Additive Manufacturing (MBAAM) system, a thin steel wall was manufactured from a low carbon steel wire. The wall was then characterized comprehensively by high-throughput high-energy X-ray diffraction (HEXRD), electron backscatter diffraction (EBSD), and in-situ HEXRD tensile tests. With the predicted temperature histories from the finite element-based additive manufacturing process simulations, the correlations between processing parameters, microstructure, and properties were established. The correlation between the final microstructure with the predicted temperature history is well explained with the material's continuous cooling transformation (CCT) diagram calculated based on the composition of the low carbon steel wire. The final microstructure is dependent on the cooling rate during austenite to ferrite/bainite transformation during initial cooling and the subsequent reheating cycles. Fast cooling rate resulted in small ferrite grain size and fine bainite structure at the location closest to the base plate. Slower cooling rate at the side wall location and repeated reheating cycles to the ferrite-pearlite regions resulted in all allotriomorphic (equiaxed) ferrite with medium grain size with small amount of pearlite. With no reheating cycles, the top location has the slowest cooling rate and a large grained allotriomorphic ferrite and bainitic structures. The measured mechanical strength is then related to the microstructural feature size (grain or lath size) observed in those locations. Thus, a good correlation is found between the mechanical properties, microstructure features and the temperature history at various locations of the printed wall.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Improved performance of HgCdTe infrared detector focal plane arrays by modulating light field based on photonic crystal structure

An HgCdTe long-wavelength infrared focal plane array photodetector is proposed by modulating light distributions based on the photonic crystal. It is shown that a promising prospect of improving performance is better light harvest and dark current limitation. To optimize the photon field distributions of the HgCdTe-based photonic crystal structure, a numerical method is built by combining the finite-element modeling and the finite-difference time-domain simulation. The optical and electrical characteristics of designed HgCdTe mid-wavelength and long-wavelength photon-trapping infrared detector focal plane arrays are obtained numerically. The results indicate that the photon crystal structure, which is entirely compatible with the large infrared focal plane arrays, can significantly reduce the dark current without degrading the quantum efficiency compared to the regular mesa or planar structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Full-field quantitative visualization of shock-driven pore collapse and failure modes in PMMA

The dynamic collapse of pores under shock loading is thought to be directly related to hot spot generation and material failure, which is critical to the performance of porous energetic and structural materials. However, the shock compression response of porous materials at the local, individual pore scale is not well understood. This study examines, quantitatively, the collapse phenomenon of a single spherical void in PMMA at shock stresses ranging from 0.4 to 1.0 GPa. Using a newly developed internal digital image correlation technique in conjunction with plate impact experiments, full-field quantitative deformation measurements are conducted in the material surrounding the collapsing pore for the first time. The experimental results reveal two failure mode transitions as shock stress is increased: (i) the first in situ evidence of shear localization via adiabatic shear banding and (ii) dynamic fracture initiation at the pore surface. Numerical simulations using thermo-viscoplastic dynamic finite element analysis provide insights into the formation of adiabatic shear bands (ASBs) and stresses at which failure mode transitions occur. Further numerical and theoretical modeling indicates the dynamic fracture to occur along the weakened material inside an adiabatic shear band. Finally, analysis of the evolution of pore asymmetry and models for ASB spacing elucidate the mechanisms for the shear band initiation sites, and elastostatic theory explains the experimentally observed ASB and fracture paths based on the directions of maximum shear.

42 ENGINEERING↗
DOE OSTI2024

Behavior of strain stripe networks in barium titanate nanocrystals on crossing its ferroelectric phase transition

Nanoscale strain networks are reported in BaTi O 3 (BTO) crystals of 300 nm size using Bragg coherent diffractive imaging. BTO nanocrystals with clear facets were chosen to identify the crystallographic directions, allowing the strain field direction and periodicity to be studied in detail. Stripes of strain were observed, which were both stable and preserved in tetragonal and cubic phases at elevated temperatures, above the tetragonal-to-cubic phase transition. A finite element analysis approach was used to simulate the domain structures inside a BTO crystal and to understand the origins of the strain stripes as piezoelectric blocks. A 180 ∘ domain model gives a better qualitative match to the experimental images. Published by the American Physical Society 2024

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Errors in reconstruction of dichroic X-ray orientation tomography due to polarization rotation of the incident beam

Dichroic X-ray tomography is a technique in which the crystal orientation or magnetization of a sample is resolved in three dimensions. The best-known uses of this technique are for observation of magnetic moments via circular dichroism, using left- and right-handed circularly polarized X-ray beams. Another variant uses linear dichroism to resolve the crystal orientation. In both these techniques, it is assumed that the absorption of X-rays along a path inside a material can be computed as a line integral of a local absorption coefficient along the ray path. For linear dichroism, this assumption is inaccurate because the polarization of the beam changes along the propagation direction when the optic axis of the material is not aligned along the polarization. In this work, a finite-element Maxwell solver is used to simulate tomography and reconstructions. The propagation effect can lead to significant errors in the reconstructed orientations. These errors may be mitigated by taking data at additional angles or by operating at energies at which the dichroism is weak. An iterative approach is proposed which may allow accurate reconstruction with fewer data than would otherwise be required.

X-ray linear dichroism↗
DOE OSTI2025

Rapid Prototyping Techniques for Organic Direct-Bonded-Copper Power Modules

Organic Direct-Bonded-Copper (ODBC) is a novel packaging technology for power modules which allows higher flexibility in layout design. In this work, a set of rapid prototyping techniques are developed for ODBC modules based on a polyimide dielectric material. These techniques enable fast and low-cost fabrication of modules with 3-Dimensional (3D) layout features. An example half-bridge (HB) silicon carbide (SiC) metal oxide semiconductor field effect transistor (MOSFET) module is designed and prototyped using the proposed techniques aiming at ultra-low power loop inductance. Finite element analysis (FEA)-circuit co-simulations results and experimental results validate an approximately 0.71nH power loop inductance for the example module design.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Embedded Fiber Bragg Grating (FBG) Sensors Fabricated by Ultrasonic Additive Manufacturing for High-Frequency Dynamic Strain Measurements

This article demonstrates high-frequency dynamic strain measurements using fiber Bragg grating (FBG) sensors embedded in metal parts. Using an ultrasonic additive manufacturing (UAM) process, FBGs inscribed in polyimide coated optical fibers were embedded in aluminum parts. An electromagnetic shaker was used to exert dynamic events on the embedded FBG sensors with frequencies from 1 to 10 kHz. The high-speed interrogation of FBG sensors was accomplished using a tunable vertical-cavity surface-emitting laser (VCSEL) and a high-speed interrogation system sampling at 120 kHz. The strain response measured by the FBG sensors was compared with real-time measurements using a laser velocimeter. A finite-element analysis (FEA) was performed to simulate responses to both static strain and high-frequency dynamic strain. In conclusion, results show that strains as small as 2.5$με$ can be resolved at frequencies up to 10 kHz.

36 MATERIALS SCIENCE↗
DOE OSTI2018

High-Order Multi-Material ALE Hydrodynamics

Here, we present a new approach for multi-material arbitrary Lagrangian--Eulerian (ALE) hydrodynamics simulations based on high-order finite elements posed on high-order curvilinear meshes. The method builds on and extends our previous work in the Lagrangian and remap phases of ALE, and depends critically on a functional perspective that enables subzonal physics and material modeling. Curvilinear mesh relaxation is based on node movement, which is determined through the solution of an elliptic equation. The remap phase is posed in terms of advecting state variables between two meshes over a fictitious time interval. The resulting advection equation is solved by a discontinuous Galerkin (DG) formulation, combined with a customized Flux Corrected Transport (FCT) type algorithm. Because conservative fields are remapped, additional synchronization steps are introduced to preserve bounds with respect to primal fields. These steps include modification of the low-order FCT solutions, definition of conservative FCT fluxes based on primal field bounds, and monotone transitions between primal and conservative fields.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Modeling Isotope Separation in Electrochemical Lithium Deposition

Naturally occurring Li consists of two stable isotopes, 6 Li with an abundance of about 7.5%, and 7 Li making up the remainder with 92.5%. The development of a 6 Li enrichment technique, in terms of technical reliability and environmental safety to reach 6 Li future requirements, represents a key step in the roadmap for nuclear fusion energy supply worldwide. This paper uses finite element analysis-based models to simulate electrochemical Li isotope separation, which is an attractive method in terms of simplicity, safety, and scalability. In the model, we quantitatively analyze how different electrochemical factors including thermodynamics, charge-transfer kinetics, and diffusivities affect the separation process (separation factor), together with cell parameters, such as cell length and current density. The maximum separation factor of 1.128 could be obtained with the cell under the optimal thermodynamic, kinetic, and diffusive conditions, which is among the highest separation factors ever reported. Furthermore, these results will assist in designing the actual isotope separation setup with large separation factor and appropriate timing for sample collection.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2017

Affect v. 0.1

Affect is a module written in the Python language that provides functions for post-processing data that are the result of computer simulations—usually the result of the finite element method—for engineering and physics systems and environments.

Copps, Kevin↗
DOE OSTI2016

Shear dominated failure in the 'hat' specimen from the 2013 Sandia Fracture Challenge.

The objective of this memo is to present a brief report of the progress achieved during FY2016 on the investigation of ductile failure in the 2013 Sandia Fracture Challenge specimen. The experimental investigation was conducted with both the original steel A286 material used in the fracture challenge as well as with Al 7075-T651. The new results include further microscopy work for the steel A286 specimens, failure criterion verification for both materials and the implementation of a finite element model containing `material imperfections' to simulate the limit load in the response of the steel A286 specimens. Funding used to conduct the work presented here was provided by the ASC V&V program on validation of shear failure (Benjamin Reedlunn, PI) and from Sandia's LDRD program.

36 MATERIALS SCIENCE↗