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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Evaluation of Stimulation by Shear Slip in Fractured Rock Using a 3D Coupled Thermo-Poro-Mechanical FEM

This paper presents a 3D finite element model for simulating shear slip on fractures in rock in response to coupled thermo-poro-mechanical processes. The modeling is part of a broader effort to investigate the role of pore pressure and cooling by cold water injection on shear slip and permeability increase in granitic rock using laboratory shearing experiments under triaxial conditions. In particular, 3D thermo-poro-mechanical finite element modeling and analysis of injection experiment in fractured rock has been carried out to analyze the role of pore pressure, temperature, fracture deformation and their interactions. 4-noded tetrahedral elements are employed for intact rock deformation and transport processes within the matrix. To represent the mechanical response of a fracture surface, zero thickness contact interface element is developed based on recently developed element partitioning algorithm and procedures for applying hydraulic pressure on the fracture surfaces. A Mohr-Coulomb type criterion is adapted to capture the slip on the fracture and to quantify its reversible and irreversible deformation. The numerical model has been calibrated using results from well controlled, advanced laboratory experiments. Excellent agreement between modeling and experimental observations is achieved. Simulation results illustrate that pore pressure increase and rock matrix/fracture surface cooling cause the fracture system to deform and slip. Fracture slip is promoted due to its normal stress reduction associated with cooling effect of cold fluid injection. The numerical model provides a physics-based understanding of the role of coupled processes on shear stimulation phenomenon and the resulting permeability enhancement.

58 GEOSCIENCES↗
DOE OSTI2019

Residual Stress Asymmetry in Thin Sheets of Double-Sided Shot Peened Aluminum

The local mechanical properties and residual stresses of a double-sided shot-peened aluminum alloy AA7050-T7451 sheet were measured by nanoindentation and simulated by the finite element method (FEM). The variation of the pileup area around the indentation due to residual stresses, as well as changes in hardness as a function of cold working, was assessed. Numerical simulations were conducted to investigate the residual stress profile. When the indentation contact area was corrected to account for pileup induced from residual stresses, the hardness increased about 10% over the range of depths at which compressive residual stresses exist, indicating some strain hardening did occur. The maximum compressive residual stress for the second shot-peened surface is on the order of 10% higher than the first surface; this is ascribed to the sequencing during the peening process, in which the first side peening influenced the stress throughout the part thickness for these thin plate structures. The ability to couple FEM simulations and rapid experimental assessments of residual stresses and hardness allows for future prediction of part thickness effects in double-sided peening processes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of sintering temperature on adhesion of spray-on piezoelectric transducers

Conventionally sol-gel spray-on transducers require a high-temperature (> 700 ◦C) sintering process; however, this process can affect the microstructure of the substrate material. For mechanical elbows and valves utilized for fluid transport in the energy sector, the components are designed to have a specific microstructure, and deviations from these specifications can create weak points in the system. For this reason it is important to investigate how the temperature of the deposition process affects the substrate. This paper investigates the effect of high-temperature and low-temperature (< 150 ◦C) processing conditions on the surface composition of the substrate. Furthermore, the resultant transducers from high- and low-temperature fabrication processes are compared to determine if a low-temperature processing method is feasible. For these studies a sol-gel spray-on process is employed to deposit piezoelectric ceramics onto a stainless-steel 316L substrate. Energy-dispersive X-ray spectroscopy is utilized to determine the composition of the substrate surface before and after transducer deposition. Results indicate that the high-temperature processing conditions may alter the surface composition of the metal due to a diffusion of the metal into the ceramic, which results in a metal surface that is bonded to the ceramic. Furthermore, it is shown that low-temperature processing of spray-on transducers is a viable method for transducer fabrication where the resultant transducers meet the industry minimum requirement of 30 dB signalto-noise ratio. In parallel simulation calculations, finite-element method (FEM) studies were performed to model the adhesive strength of the low-temperature processed transducer to the substrate surface. Comparisons between the simulations and experiments suggest that the bond strength is much greater than the commercial gel bonds and closer to hardened epoxy glue bonds. These results indicate that spray-on transducers fabricated under lowtemperature processing conditions are a viable solution for leave-in-place monitoring of structures.

M. Sinding, Kyle↗
DOE OSTI2022

Modeling of RMS Current in CSI Filter Capacitor and Minimum Conduction Loss Operation of CSI-Fed PMSM Drives for Traction Applications

This paper presents a new minimum conduction loss (MCL) torque control algorithm for current source inverter- (CSI-) fed permanent-magnet synchronous machine drives that improves the overall machine drive efficiency by minimizing the combined conduction losses in the inverter and machine at each operating point. First, analytical models of conduction losses in CSI-fed motor drives are presented, and a closed-form expression for optimal d-axis stator current to achieve MCL operation is derived. An expression for the rms current in both wye- and delta-connected CSI output capacitors is derived. A detailed simulation model that emulates the operation of a 100 kW SiC CSI-fed integrated machine drive (IMD) has been developed based on experimental results and finite element analysis. This simulation model is used to evaluate the proposed MCL control algorithm applied to a CSI-fed IMD system for a battery-electric vehicle traction drive, and the predicted total drive system loss with MCL control is compared to predicted losses with maximum-torque-per-ampere (MTPA) and minimum dc-link current (MDCC) control. Results show that MCL control can achieve loss reductions compared to the other two control algorithms over a wide range of operating conditions, with significant loss reductions >20% in the medium-speed regime.

42 ENGINEERING↗
NASA NTRS2001

Evaluation of Test/Analysis Correlation Methods for Crash Applications

A project has been initiated to improve crash test and analysis correlation. The work in this paper concentrated on the test and simulation results for a fuselage section. Two drop tests of the section were conducted. The first test was designed to excite the linear structural response for comparison with finite element modal analysis results. The second test was designed to provide data for correlation with crash simulations. An MSC.Dytran model was developed to generate nonlinear transient dynamic results. Following minor modifications, the same model was executed in MSC.Nastran to generate modal analysis results. The results presented in this paper concentrate on evaluation of correlation methodologies for crash test data and finite element simulation results.

Lyle, Karen H.↗
DOE OSTI2023

Delamination imaging in composites using cross-correlation method by non-contact air-coupled Lamb waves

Abstract This paper presents a cross-correlation damage detection technique using damage scattered waves of all directions cross-correlate with the incident waves for delamination imaging in carbon fiber reinforced polymer (CFRP) composite plate. The presented imaging can not only detect the delamination but also precisely quantify the damage with locations and sizes. First, a CFRP composite sample with interlaminar delamination by inserting the Teflon layer is designed and manufactured in house using the hot-press machine. And a three-dimensional model of composite sample is simulated using the finite element method. Next, the cross-correlation imaging algorithm is introduced in detail, and the cross-correlation imaging proof-of-concept study is conducted with the simulated Lamb waves in the composite sample. Finally, a fully non-contact air-coupled transducer and scanning laser Doppler vibrometer system with a single-mode Lamb wave method is established to actuate and sense the interrogating Lamb waves in the structure. The imaging method is experimentally implemented for the one delamination and two delaminations imaging and quantification.

Instruments & Instrumentation↗
DOE OSTI2019

Simulation of discrete cracks driven by nearly incompressible fluid via 2D combined finite-discrete element method

Here, we propose a novel hydraulic solver in order to simulate key mechanisms that control fluid-driven cracks in the framework of the combined finite-discrete element method (FDEM). The main innovative aspect of the present work is the independence of the fluid's critical time step size with the fracture opening. This advantage is extremely important because it means that very fine meshes can be used around areas of interest, such as boreholes, without penalizing the computational cost as fractures propagate (ie, open) and the fluid flows through them. This is a great advantage over other recently introduced approaches that exhibit a dependency of the time step in the form of Δt crit ∝ (l/a) 2 where l is the element size and a is the fracture aperture. This paper presents a series of benchmark cases for the proposed solver. The rationale adopted by the authors was to benchmark and validate the implementation of the hydraulic solver in an incremental fashion, starting from the simplest cases and building in complexity. The calculations shown in this work clearly demonstrate that the proposed approach is able to reproduce analytical results for fluid flow through a single crack. The results presented in this paper also demonstrate that the new approach is robust enough to deal with complex fracture patterns and complex geometries; the obtained fluid-driven fracture patterns in the vicinity of a borehole certainly stand to the scrutiny of human visual perception.

58 GEOSCIENCES↗
NASA NTRS2012

Determining Wheel-Soil Interaction Loads Using a Meshfree Finite Element Approach Assisting Future Missions with Rover Wheel Design

A wheel experiencing sinkage and slippage events poses a high risk to rover missions as evidenced by recent mobility challenges on the Mars Exploration Rover (MER) project. Because several factors contribute to wheel sinkage and slippage conditions such as soil composition, large deformation soil behavior, wheel geometry, nonlinear contact forces, terrain irregularity, etc., there are significant benefits to modeling these events to a sufficient degree of complexity. For the purposes of modeling wheel sinkage and slippage at an engineering scale, meshfree finite element approaches enable simulations that capture sufficient detail of wheel-soil interaction while remaining computationally feasible. This study demonstrates some of the large deformation modeling capability of meshfree methods and the realistic solutions obtained by accounting for the soil material properties. A benchmark wheel-soil interaction problem is developed and analyzed using a specific class of meshfree methods called Reproducing Kernel Particle Method (RKPM). The benchmark problem is also analyzed using a commercially available finite element approach with Lagrangian meshing for comparison. RKPM results are comparable to classical pressure-sinkage terramechanics relationships proposed by Bekker-Wong. Pending experimental calibration by future work, the meshfree modeling technique will be a viable simulation tool for trade studies assisting rover wheel design.

Contreras, Michael T.↗
DOE OSTI2019

Study of 3D Pores and Its Relationship with Crack Initiation Factors of Aluminum Alloy Die Castings

Observing different types of pores in aluminum alloy from A356 high-pressure die casting by X-ray computerized tomography and scanning electron microscopy, we found that equivalent spherical diameters of pores obeyed the 3-parameter lognormal distribution rather than the 2-parameter one. Statistically, the sphericity of most round gas pores with smooth inner surfaces consisting of oxide films was > 0.58, while the sphericity of shrinkage pores with unfractured dendrites was < 0.40. Sphericity range for the gas-shrinkage pores, composed of gas pores and shrinkage pores with long tails, was 0.40-0.58. Finite element analysis with simulation of the Von Mises stress distribution around the actual 3D pores among crack generation was performed. Tensile test simulation showed that cracks often originated from the shrinkage pores and the gas-shrinkage pores, rather than from the gas pores. This method also can be used to predict fracture strength of samples with pores. For our samples, results obtained using this method showed that simulated fracture strength was slightly higher than experimentally determined strength.

36 MATERIALS SCIENCE↗
NASA NTRS2009

Nonlinear Reduced-Order Simulation Using Stress-Free and Pre-Stressed Modal Bases

A study is undertaken to determine the advantages and disadvantages associated with application of stress-free and pre-stressed modal bases in a reduced-order finite-element-based nonlinear simulation. A planar beam is chosen as an application example and its response due to combined thermal and random pressure loadings is examined. Combinations of two random pressure levels and two thermal conditions are investigated. The latter consists of an ambient temperature condition and an elevated temperature condition in the post-buckled regime. It is found that stress-free normal modes establish a broadly applicable modal basis yielding accurate results for all the loading regimes considered. In contrast, the range of applicability for a thermally pre-stressed modal basis is found to be limited. The behavior is explained by scrutinizing the coupling found in the linear stiffness and the effect this coupling has on the structural response characteristics under the range of loading conditions considered.

Przekop, Adam↗
DOE OSTI2015

FEM design and simulation of a short, 10 MV, S-band Linac with Monte Carlo dose simulations

Purpose: Current commercial 10 MV Linac waveguides are 1.5 m. The authors’ current 6 MV linear accelerator–magnetic resonance imager (Linac–MR) system fits in typical radiotherapy vaults. To allow 10 MV treatments with the Linac–MR and still fit within typical vaults, the authors design a 10 MV Linac with an accelerator waveguide of the same length (27.5 cm) as current 6 MV Linacs. Methods: The first design stage is to design a cavity such that a specific experimental measurement for breakdown is applicable to the cavity. This is accomplished through the use of finite element method (FEM) simulations to match published shunt impedance, Q factor, and ratio of peak to mean-axial electric field strength from an electric breakdown study. A full waveguide is then designed and tuned in FEM simulations based on this cavity design. Electron trajectories are computed through the resulting radio frequency fields, and the waveguide geometry is modified by shifting the first coupling cavity in order to optimize the electron beam properties until the energy spread and mean energy closely match values published for an emulated 10 MV Linac. Finally, Monte Carlo dose simulations are used to compare the resulting photon beam depth dose profile and penumbra with that produced by the emulated 10 MV Linac. Results: The shunt impedance, Q factor, and ratio of peak to mean-axial electric field strength are all matched to within 0.1%. A first coupling cavity shift of 1.45 mm produces an energy spectrum width of 0.347 MeV, very close to the published value for the emulated 10 MV of 0.315 MeV, and a mean energy of 10.53 MeV, nearly identical to the published 10.5 MeV for the emulated 10 MV Linac. The depth dose profile produced by their new Linac is within 1% of that produced by the emulated 10 MV spectrum for all depths greater than 1.5 cm. The penumbra produced is 11% narrower, as measured from 80% to 20% of the central axis dose. Conclusions: The authors have successfully designed and simulated an S-band waveguide of length of 27.5 cm capable of producing a 10 MV photon beam. This waveguide operates well within the breakdown threshold determined for the cavity geometry used. The designed Linac produces depth dose profiles similar to those of the emulated 10 MV Linac (waveguide-length of 1.5 m) but yields a narrower penumbra.

07 ISOTOPE AND RADIATION SOURCES↗
NASA NTRS2000

Efficient Use of Distributed Systems for Scientific Applications

Distributed computing has been regarded as the future of high performance computing. Nationwide high speed networks such as vBNS are becoming widely available to interconnect high-speed computers, virtual environments, scientific instruments and large data sets. One of the major issues to be addressed with distributed systems is the development of computational tools that facilitate the efficient execution of parallel applications on such systems. These tools must exploit the heterogeneous resources (networks and compute nodes) in distributed systems. This paper presents a tool, called PART, which addresses this issue for mesh partitioning. PART takes advantage of the following heterogeneous system features: (1) processor speed; (2) number of processors; (3) local network performance; and (4) wide area network performance. Further, different finite element applications under consideration may have different computational complexities, different communication patterns, and different element types, which also must be taken into consideration when partitioning. PART uses parallel simulated annealing to partition the domain, taking into consideration network and processor heterogeneity. The results of using PART for an explicit finite element application executing on two IBM SPs (located at Argonne National Laboratory and the San Diego Supercomputer Center) indicate an increase in efficiency by up to 36% as compared to METIS, a widely used mesh partitioning tool. The input to METIS was modified to take into consideration heterogeneous processor performance; METIS does not take into consideration heterogeneous networks. The execution times for these applications were reduced by up to 30% as compared to METIS. These results are given in Figure 1 for four irregular meshes with number of elements ranging from 30,269 elements for the Barth5 mesh to 11,451 elements for the Barth4 mesh. Future work with PART entails using the tool with an integrated application requiring distributed systems. In particular this application, illustrated in the document entails an integration of finite element and fluid dynamic simulations to address the cooling of turbine blades of a gas turbine engine design. It is not uncommon to encounter high-temperature, film-cooled turbine airfoils with 1,000,000s of degrees of freedom. This results because of the complexity of the various components of the airfoils, requiring fine-grain meshing for accuracy. Additional information is contained in the original.

Taylor, Valerie↗
DOE OSTI2019

Potting Adhesion Investigation Using Double Cantilever Beam Tests

The adhesive properties of an elastomeric potting material were investigated experimentally using double cantilever beam (DCB) specimens. In these tests, the applied load and opening displacement were measured, with the goal of quantitatively characterizing the adhesive strength of the thin layer of potting material. The DCB tests were also simulated using a finite element model that incorporated an adhesive traction-separation model to represent the compliance of the potting layer and evolving (weakening) strength of the potting-to-metal adhesive interfaces. The test and simulation results were used to determine a range of traction-separation model parameters, expressed in the form of correlated probability distributions, that constitute a good calibration of the model over the range of variability exhibited in the test series.

36 MATERIALS SCIENCE↗
DOE OSTI2018

The effect of cell-size dispersity on the mechanical properties of closed-cell aluminum foam

Highlights: • The effect of cell-size dispersity on the Young’s modulus of foam was investigated through experiments and simulations • Crushing experiments were conducted on 27 foam specimens in order measure the relative Young’s modulus • Three types of numerical models, including monodispersed, bidispersed and polydispersed models were studied • The numerical convergence on the model size for different models was also studied • The difference among numerical results obtained from the three models was insignificant - Abstract: The mechanical properties of open-celled and closed-celled metallic foam are strongly influenced by the cellular microstructure and material properties of the base material from which the foam is made. The present study investigated the microstructure and mechanical properties of ALPORAS closed-cell Al foam. In particular, the effect of cell-size dispersity on the Young's modulus of foam was investigated through crushing experiments and numerical simulations. Crushing experiments were first conducted on two mutually perpendicular directions. Numerical models of varying cell-size dispersity and complexity were subsequently developed, including the monodispersed model that is based on the Kelvin structure, the bidispersed model, and the polydispersed model that is based on the random Laguerre tessellation algorithm. The numerical models were discretized using shell elements, and the simulations were conducted using the ABAQUS finite element method software. Relative Young's moduli with various relative densities were predicted, and the numerical convergence on the model size was studied accordingly. Numerical results obtained from all three numerical models developed were found to overestimate the experimental results by a factor of approximately 3. The difference between the numerical and experimental results can be attributed to the difference in detailed microstructural characteristics. Furthermore, the difference among numerical results obtained from the three models was insignificant.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Surrogate Model Integration with MOOSE XFEM for Creep Crack Growth

Ferritic-martensitic steels are key structural materials for advanced reactors but experience time-dependent deformation and damage under prolonged high temperature and irradiation, leading to creep-driven crack initiation and growth. High-fidelity models—crystal plasticity with irradiation mechanisms, phase-field for microstructural evolution, and continuum-damage viscoplasticity—capture the underlying physics but are too computationally intensive for broad design-space exploration and uncertainty quantification. This milestone advances a scalable alternative by integrating a microstructure-sensitive surrogate creep model into the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework and extending it to fracture via the extended finite element method (XFEM). The surrogate model, developed with collaborators at Sandia and Los Alamos National Laboratories, maps relevant microstructural descriptors to the viscoplastic response of HT9. We embed this surrogate within a coupled deformation-damage workflow in MOOSE/XFEM to simulate creep-driven crack initiation and propagation. Implementation enhancements include updates to the material interface, a plastic correction phase involving microstructure evolution, and fracture criteria to ensure numerical robustness and compatibility with the surrogate structure. Demonstrations on canonical creep benchmarks spanning uniaxial and multiaxial states show that the surrogate reproduces key trends of high-fidelity models while substantially reducing computational cost. The resulting capability bridges physics fidelity and performance, providing a practical path to a predictive, microstructure-aware assessment of creep and fracture in reactor materials.

36 - MATERIALS SCIENCE↗
DOE OSTI2018

Modeling Resonances of Dielectric Spheres in the Terahertz Range

We have studied the resonance properties of dielectric spheres with subwavelength dimensions in the optical and terahertz frequency range. Distribution of the electric field vector in a sphere irradiated by a plane polarized wave has been numerically simulated using the finite element method. The obtained results clearly demonstrate that it is in principle possible to use resonances of various orders in dielectric spheres for the amplification of current in a conducting antenna with dimensions several times smaller than the radiation wavelength. In the case of resonances within short spectral intervals, it is possible to provide a current density gain up to 25–30 dB as compared to that in the absence of a dielectric sphere surrounding the antenna.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Evaluation of Structural Lattices for a Davis Gun Earth Penetrator Impact Experiment

The advanced materials team investigated the use of additively manufactured metallic lattice structures for mitigating impact response in a Davis gun earth penetrator impact experiment. High-fidelity finite element models were developed and validated with quasistatic experiments. These models were then used to simulate the response of such lattices when subjected to the acceleration loads expected in the Davis gun experiment. Results reveal how the impact mitigation performance of lattices can change drastically at a certain relative density. Based on these observations, an experiment deck was designed to probe the response of lattices with different relative densities during the Davis gun phase 2 shots. The expected performance of these lattices is predicted before testing based on simulation results. The results of the Davis gun phase 2 shots are expected to provide data which will be used to assess the predictive capability of the finite element simulations in such a complex impact environment.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simulated effect of defect volume and location on very high cycle fatigue of laser beam powder bed fused AlSi10Mg

This study quantifies the interaction between volumetric defect location and size on the very high cycle fatigue (VHCF) of laser beam powder bed fused (LB-PBF) AlSi10Mg. Crystal plasticity finite element method (CPFEM) simulations were used to investigate the effects of defect location and size on the driving force for crack initiation. The CPFEM model was calibrated against uniaxial and cyclic experimental data of LB-PBF AlSi10Mg. Defect characteristics were informed by experimental data from the specimens produced in various geometries to create realistic representative volume elements (RVEs) with equivalent volume fractions of defects. By embedding defects of varying sizes and locations within the RVEs, fatigue indicator parameters (FIPs) were calculated to analyze the impact of defects’ characteristics on fatigue performance. Different combinations of defect volume and locations were generated for various microstructure instantiations, providing insight into extreme value fatigue responses. Larger defect volumes located on free surfaces consistently generated the highest FIPs, suggesting defect size and boundary proximity intensify stress concentration effects. RVEs with multiple smaller defects produced lower FIPs than those with single large critical defects. These findings underscore the critical role of defect characteristics on fatigue life, providing a foundation for future predictive modeling in fatigue-sensitive AM applications.

AlSi10Mg↗