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.
Towards a modified Immersed Finite Element Method for blast-on-structure simulations
Explore the source record for details and available documents.
Finite Element Basis Compression for Extreme-Scale Simulations and R-Adaptive Mesh Optimization to Enable Compression
Explore the source record for details and available documents.
Hyperbolic heat-conduction problems - Numerical simulations via explicit Lax-Wendroff-based finite element formulations
Explore the source record for details and available documents.
ZIP2DL: An Elastic-Plastic, Large-Rotation Finite-Element Stress Analysis and Crack-Growth Simulation Program
ZIP2DL is a two-dimensional, elastic-plastic finte element program for stress analysis and crack growth simulations, developed for the NASA Langley Research Center. It has many of the salient features of the ZIP2D program. For example, ZIP2DL contains five material models (linearly elastic, elastic-perfectly plastic, power-law hardening, linear hardening, and multi-linear hardening models), and it can simulate mixed-mode crack growth for prescribed crack growth paths under plane stress, plane strain and mixed state of stress conditions. Further, as an extension of ZIP2D, it also includes a number of new capabilities. The large-deformation kinematics in ZIP2DL will allow it to handle elastic problems with large strains and large rotations, and elastic-plastic problems with small strains and large rotations. Loading conditions in terms of surface traction, concentrated load, and nodal displacement can be applied with a default linear time dependence or they can be programmed according to a user-defined time dependence through a user subroutine. The restart capability of ZIP2DL will make it possible to stop the execution of the program at any time, analyze the results and/or modify execution options and resume and continue the execution of the program. This report includes three sectons: a theoretical manual section, a user manual section, and an example manual secton. In the theoretical secton, the mathematics behind the various aspects of the program are concisely outlined. In the user manual section, a line-by-line explanation of the input data is given. In the example manual secton, three types of examples are presented to demonstrate the accuracy and illustrate the use of this program.
Nonlinear Finite Element Model Updating Part II: Implementation and Simulation.
Abstract not provided.
Use of the Conformal Decomposition Finite Element Method for Coupled Electrochemical-Mechanical Simulations of Lithium-Ion Battery Electrodes.
Abstract not provided.
Combination of Discrete and Finite Element Methods for Coupled Electrochemical-Mechanical Simulations of Lithium-Ion Battery Electrodes.
Abstract not provided.
Two-Dimensional Dynamic Simulation of a Continuous Foil Bearing
In this paper, the two dimensional(radial and circumferential) transient Navier-Stokes equations are used to solve the hydrodynamic problem in conjunction with the time dependent motion of the journal, and the deformable, spring supported foil. The elastic deformation of the foil and its supports are simulated by a finite element model. The time-dependent Navier-Stokes formulation is used to solve for the interaction between the fluid lubricant, the motion of the journal and the deformable foil boundary. The steady state, the quasi-transient and the full transient dynamic simulation of the foil-fluid journal interaction are examined on a comparative basis. For the steady state simulation, the fluid lubricant pressures are evaluated for a particular journal position, by means of an iterative scheme until convergence is achieved in both the fluid pressures and the corresponding foil deformation. For the quasi-transient case, the transient motion of the journal is calculated using a numerical integration scheme for the velocity and displacement of the journal. The deformation of the foil is evaluated through numerical iteration in feedback mode with the fluid film pressure generated by the journal motion until convergence at every time step is achieved. For the full transient simulation, a parallel real-time integration scheme is used to evaluate simultaneously the new journal position and the new deformed shape of the foil at each time step. The pressure of the fluid lubricant is iterated jointly with the corresponding journal position and the deformed foil geometry until convergence is achieved. A variable time-stepping Newmark-Beta integration procedure is used to evaluate the transient dynamics at each time step of the bearing.
Performance analysis of higher mode spoof surface plasmon polariton for terahertz sensing
We investigated the spoof surface plasmon polaritons (SSPPs) on 1D grooved metal surface for terahertz sensing of refractive index of the filling analyte through a prism-coupling attenuated total reflection setup. From the dispersion relation analysis and the finite element method-based simulation, we revealed that the dispersion curve of SSPP got suppressed as the filling refractive index increased, which cause the coupling resonance frequency redshifting in the reflection spectrum. The simulated results for testing various refractive indexes demonstrated that the incident angle of terahertz radiation has a great effect on the performance of sensing. Smaller incident angle will result in a higher sensitive sensing with a narrower detection range. In the meanwhile, the higher order mode SSPP-based sensing has a higher sensitivity with a narrower detection range. The maximum sensitivity is 2.57 THz/RIU for the second-order mode sensing at 45° internal incident angle. The proposed SSPP-based method has great potential for high sensitive terahertz sensing.
Uncertainty Quantification Framework for Predicting Material Response with Large Number of Parameters: Application to Creep Prediction in Ferritic-Martensitic Steels Using Combined Crystal Plasticity and Grain Boundary Models
This paper presents an uncertainty quantification (UQ) framework for the physics-based model prediction of material response with a large number of parameters. The application problem presented in this work is that of predicting creep in Grade 91 steel at 600°C. The material response is defined with a physically based microstructural model with constitutive equations emulating several observed phenomena in Grade 91 and embodied into an explicit geometry mesoscale finite element model for prior austenite grains and grain boundaries. Creep within the grains and in grain boundaries are represented by crystal plasticity for dislocation motion and a physics-based model for cavity growth and nucleation, respectively. The creep behavior of this material is influenced by several parameters, some of which have a wide range of variation based on experimental data. UQ combined with microstructural modeling can discover the core microstructural causes of experimental variability, leading to improved materials with lower variability in critical long-term material properties. In this study, we investigate the model's uncertainty to identify material properties that may be modified during production to increase creep life and analyze different components of the crystal plasticity model for improvements. For this purpose, a quantity of interest is defined as time to minimum creep rate, which correlates well to the creep failure of the material. A deep neural network model was trained and validated to be used as a surrogate for the finite element model. Then, a variance-based sensitivity analysis is performed on the surrogate model to find the Sobol indices of the input parameters in respect to the output quantity of interest. The Sobol indices are used to reduce the dimensionality of the model. Generalized polynomial chaos expansion is used on the reduced basis models to propagate the uncertainty from the input parameters to the quantity of interest using the deep neural network surrogate model. These results are benchmarked against uncertainty propagation using Monte Carlo simulations. In conclusion, the UQ performed through the reduced basis model captures almost all the uncertainty in the model with significantly fewer simulations, making it possible to perform the UQ directly via simulations with the finite element model rather than surrogate machine-learned models.
Simulation of the Fast Reactor Fuel Assembly Duct-Bowing Reactivity Effect Using Monte Carlo Neutron Transport and Finite Element Analysis
This paper discusses a new method of simulating the fuel assembly duct-bowing reactivity coefficient for EBR-II run 138B. Quantification of the fuel assembly duct-bowing reactivity effect in liquid metal–cooled fast reactors has been a persistent problem since they were first designed and operated. Simulation of the duct-bowing reactivity effect is difficult because the level of detail required to simulate the effect has exceeded most modeling capabilities. The new method outlined in this paper utilizes the finite element analysis code ANSYS to analyze the thermal and structural components. Here, the displacement of the fuel assembly duct due to thermal expansion and mechanical interaction was calculated by ANSYS using recorded EBR-II run 138B temperature and power boundary value data. The displacement values were incorporated into to a Monte Carlo model of EBR-II run 138B and keff was calculated. Multiple Monte Carlo calculations were performed with duct displacement values corresponding to different reactor temperatures. Using the calculated keff values associated with the different duct displacement results allowed calculation of the duct-bowing reactivity coefficient. The duct-bowing reactivity coefficient was calculated to be –14.5 × 10 –4 $/°C/ ± 4.4%.
Optimization of Type 4 composite pressure vessels using genetic algorithms and simulated annealing
We propose and demonstrate two methodologies for weight minimization of Type 4 (plastic lined, fiber wrapped) compressed hydrogen pressure vessels: genetic algorithms and simulated annealing. We consider 70 MPa vessels with a safety factor of 2.25, and analyze the vessels with classical laminate theory. We propose an objective function based on Tsai-Wu criterion, composite thickness, a safety factor, and a penalization factor. Here, the optimum results are analyzed and compared by a high resolution finite element model. Computer simulations show that the proposed methodology produces more efficient designs by reducing the weight by up to 9.8% and 11.2% when compared to previously published vessel optimization research.
Stress Relaxation Cracking of Alloys at Temperatures Higher Than 540°C
Type 347H stainless steel (347H SS), used in commercial concentrating solar power (CSP) thermal energy storage to store solar-salt at a temperature of 565°C, has been reported in the literature to be susceptible to stress-relaxation cracking (SRC). The welded heat-affected zone (HAZ) and fusion zone (FZ) of 347H SS, particularly in thick sections, are known to be susceptible to failure during post-weld heat treatment (reheat cracking). SRC could also occur after months or years under an elevated-temperature service environment. Two conditions must be present for failure to occur in the HAZ and/or FZ: 1) a modified or sensitized microstructure and 2) sufficiently high tensile residual stresses present at the elevated service temperature. The overarching goal of this project is to recommend SRC mitigation protocols to avoid susceptibility to fail through SRC at temperatures relevant for CSP. Post weld heat treatment conditions and alternative alloys are investigated as potential mitigation solutions to SRC. We used Gleeble thermomechanical simulation tests and finite element (FE) models to understand the susceptibility of 347H SS to SRC as a function of temperature, stress, and microstructure. We started with a literature review of weldability issues with 347H SS and techniques to mitigate SRC in 347H SS weldments. Next, we used Gleeble experiments to determine reheat cracking susceptibility in the simulated HAZ of 347SS weldments and an alternative alloy, 316L SS (NUCL 167 SPH) with boron added. We also performed Gleeble experiments to compare the reheat cracking susceptibility of 347H cross welded with two different fillers: E347, which is used in some commercial CSP TES tanks, and E16.8.2. We validated the experimental results with FE models of the residual stresses. We found that although the 316L (NUCL 167 SPH) is less susceptible to reheat cracking, it is slightly weaker than 347H and ASME BP&V codes limit its use to a service condition of 565°C. We also found that welds using E16.8.2 as the weld filler are less susceptible to failure than those using E347, likely due to the higher creep ductility of E16.8.2, and that it may be used as an alternative filler for repair welding of 347H welds or as the primary choice of filler for newly developed weld joints. We also found that post weld heat treatment could be a viable solution for mitigating stress in E347-347H SS thick, constrained welds, like those found in CSP tanks, and propose several options for mitigating SRC in existing and future TES tanks.
A numerical homogenization method for heterogeneous, anisotropic elastic media based on multiscale theory
The development of reliable methods for upscaling fine scale models of elastic media has long been an important topic for rock physics and applied seismology. Several effective medium theories have been developed to provide elastic parameters for materials such as finely layered media or randomly oriented or aligned fractures. In such cases, the analytic solutions for upscaled properties can be used for accurate prediction of wave propagation. However, such theories cannot be applied directly to homogenize elastic media with more complex, arbitrary spatial heterogeneity. We therefore propose a numerical homogenization algorithm based on multiscale finite element methods for simulating elastic wave propagation in heterogeneous, anisotropic elastic media. Specifically, our method used multiscale basis functions obtained from a local linear elasticity problem with appropriately defined boundary conditions. Homogenized, effective medium parameters were then computed using these basis functions, and the approach applied a numerical discretization that is similar to the rotated staggered-grid finite difference scheme. Comparisons of the results from our method and from conventional, analytical approaches for finely layered media showed that the homogenization reliably estimated elastic parameters for this simple geometry. Additional tests examined anisotropic models with arbitrary spatial heterogeneity where the average size of the heterogeneities ranged from several centimeters to several meters, and the ratio between the dominant wavelength and the average size of the arbitrary heterogeneities ranged from 10 to 100. Comparisons to finite-difference simulations proved that the numerical homogenization was equally accurate for these complex cases.
Cluster Dynamics Modeling Needs for the Advanced Materials and Manufacturing Technologies Program
This milestone report aims to identify and assess the cluster dynamics (CD) modeling requirements within the Department of Energy's Office of Nuclear Energy (DOE-NE) Advanced Materials and Manufacturing Technologies (AMMT) program and to communicate these needs to the DOE-NE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The goal is to ensure NEAMS is well-informed about the CD modeling requirements to support AMMT's mission of accelerating the development, qualification, demonstration, and deployment of advanced structural materials and manufacturing for nuclear energy applications. CD modeling is an essential tool for predicting the degradation of structural materials under irradiation, which is a key component of AMMT's accelerated qualification process. The AMMT program focuses on both additively manufactured and wrought structural alloys, such as laser powder-bed fusion 316H austenitic stainless steel, alloy 709, Haynes 244, and alloy 617. These materials require a generalized CD modeling framework to facilitate rapid model development and computational simulation. A flexible, generalized CD software, similar to the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework, would enable modeling of various cluster types, including defect clusters, defect-solute clusters, and multicomponent clusters, incorporating thermodynamics and kinetics parameters. Radiation effects, microstructural feature evolution, and multi-dimensional modeling are critical considerations for the CD model. The usability of the CD code should allow for easy modification and coupling with MOOSE-based simulations. Additionally, the software should adhere to Nuclear Quality Assurance-1 standards, include a testing suite for verification and validation, and be version-controlled within a national laboratory-managed Git repository. Benchmark problems are needed to assess code predictions and performance.
Hob Strength in Gear Cutting
The strength of gear hobs in cutting worm gears is calculated by the finite-element method, on the basis of SolidWorks/Simulation software.
Radiation-MHD simulations for the development of a spark discharge channel.
The growth of a cylindrical s park discharge channel in water and Lexan is studied using a series of one - dimensional simulations with the finite - element radiation - magnetohydrodynamics code ALEGRA. Computed solutions are analyzed in order to characterize the rate of growth and dynamics of the spark c hannels during the rising - current phase of the drive pulse. The current ramp rate is varied between 0.2 and 3.0 kA/ns, and values of the mechanical coupling coefficient K p are extracted for each case. The simulations predict spark channel expansion veloc ities primarily in the range of 2000 to 3500 m/s, channel pressures primarily in the range 10 - 40 GPa, and K p values primarily between 1.1 and 1.4. When Lexan is preheated, slightly larger expansion velocities and smaller K p values are predicted , but the o verall behavior is unchanged.