Search NASA⌕ Search

SEARCH · Search NASA

Results for “FEM 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

On Residual Stresses and Reference Temperatures in Thermomechanical Simulations of Photovoltaic Modules Using the Finite Element Method

Thermomechanical simulation of photovoltaic (PV) modules using the finite element method (FEM) is a useful tool to evaluate module design features in terms of structural integrity, reliability, and durability. One of the main challenges in the numerical modeling of a PV module is the incorporation of residual stresses induced by the manufacturing process. Modeling assumptions and abstractions are necessary to limit the model complexity and reduce the computational time. However, oversimplifications and incorrect assumptions can lead to erroneous numerical results. Unfortunately, much simulation work still neglects process-induced stresses. This can lead to incorrect predictions of the stress-strain history and erroneous conclusions during the design process. In this work, we review current modeling practices for incorporating process-induced stresses, and contrast numerical models that consider residual stresses with those that neglect them. We find that the simulation objective and available material properties dictate which process steps need to be modeled, and explore in depth the modeling of residual stresses induced by the lamination process. We demonstrate that a simplified cooldown procedure at the beginning of the simulation can increase the model accuracy and discuss appropriate choices for starting and reference temperatures in the finite element model.

14 SOLAR ENERGY↗
DOE OSTI2023

Results of the Micromorphic Upscaling Workflow for the PSAAP III Year 3 Report

Predicting the mechanical response and failure of heterogeneous materials has proven difficult. Multiscale numerical methods based in higher order continuum theories attempt to bridge the gap between microscale and macroscale structural behavior. Micromorphic continuum theories have shown promise. The Tardigrade software package is an implementation of Eringen’s micromorphic continuum theory with capabilities to support multiscale material modeling workflows. These include homogenization through the Micromorphic Filter, calibration of micromorphic material models, and macroscale simulation in Tardigrade-MOOSE. This work discusses micromorphic upscaling efforts of the University of Colorado Boulder PSAAP III multidisciplinary simulation center (MSC). Verification studies are presented that compare the accuracy of the upscaling workflow with analytical solutions for a trivial stress state and homogeneous material using direct numerical simulations (DNS) conducted in the Ratel finite element method (FEM) and GEOS material point method (MPM) codes. These verification studies consider upscaling using the Micromorphic Filter for a “single filter domain”. Finally, DNS of a heterogeneous composite material is upscaled using a “multiple filter domain” method.

36 MATERIALS SCIENCE↗
NASA NTRS2004

FEM and Multiphysics Applications at NASA/GSFC

FEM software available to the Mechanical Systems Analysis and Simulation Branch at Goddard Space Flight Center (GSFC) include: 1) MSC/Nastran; 2) Abaqus; 3) Ansys/Multiphysics; 4) COSMOS/M; 5) 'Home-grown' programs; 6) Pre/post processors such as Patran and FEMAP. This viewgraph presentation provides additional information on MSC/Nastran and Ansys/Multiphysics, and includes screen shots of analyzed equipment, including the Wilkinson Microwave Anistropy Probe, a micro-mirror, a MEMS tunable filter, and a micro-shutter array. The presentation also includes information on the verification of results.

Loughlin, James↗
DOE OSTI2023

COUPLING SMOOTHED PARTICLE HYDRODYNAMICS WITH FINITE ELEMENT METHOD TO SIMULATE RESIDUAL STRESSES FROM FRICTION STIR PROCESSING

Friction stir processing (FSP) is a solid-state material processing technique that locally modifies the microstructure but also induces undesirable residual stresses. A robust numerical model for the FSP can help in mitigating these residual stresses. Heat source models within a finite element method (FEM) framework suffer from inaccuracies. In contrast, smoothed particle hydrodynamics (SPH) model that explicitly captures the material flow near the tool and the associated heat generation are accurate. However, the computational expense of SPH simulations can be prohibitive. In this work, we propose a coupled SPH-FEM framework. SPH is used to model the heat generation accurately near the tool and which is then inserted into to the FEM model as a heat source. To verify this proposed coupling approach, a test case is set up with typical FSP conditions and it is modeled in both SPH and SPH-FEM. The temperatures profiles were compared after the simulations have reached steady-state temperatures. The similarity of the temperature profiles from SPH-FEA and SPH validated the proposed coupling approach. This proposed approach achieves the accuracy of the SPH method while potentially retaining the low computational expense of FEM.

smoothed particle hydrodynamics, Finite Element Me↗
DOE OSTI2026

Final Technical Report - Center for Simulation of Fusion Relevant RF Actuators

We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗
DOE OSTI2021

Powder spreading, densification, and part deformation in binder jetting additive manufacturing

Binder jetting additive manufacturing (AM) can print complex structures in economical and scalable manner. Binder jetting AM comprises of deposition and weak binding of particles, known as green part, at room temperature and subsequent binder removal and sintering densification at high temperatures. However, during the densification (i.e., sintering), the part significantly deforms due to volume shrinkage. The deformation during sintering is difficult to predict, which prevents the widespread application of this technology. In this research, powder spreading simulation using discrete element method (DEM) was performed first to capture local variations in powder bed configuration. Second, finite element method (FEM) with a phenomenological continuum constitutive model was used to predict part shrinkage during the sintering process. DEM simulation showed variations in packing density, particle segregation, formation of uneven powder bed surface, and shift in particle size distribution (PSD). The sintering simulation modeled part deformation with a reasonable accuracy of 3% for solid-state sintering and intermediate liquid phase sintering. A demonstration case with non-uniform initial packing density showed that inhomogeneous green part density and PSD should be accounted for prediction of part deformation in binder jetting AM.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Finite Element Modeling of Plastic Deformation During Spin Forming of Aluminum 6061-O

Spin and flow forming are metal deformation techniques in which a disk or tube of material is radially thinned and axially lengthened over a rotating mandrel. A finite-element continuum model (FEM) based on commercial software (DEFORM®) was developed to determine the stress/strain distribution and damage accumulation in a flow-formed, near-net-shape part. Computational validation and experimental verification of the model was leveraged to assess the effectiveness of computer simulations. Quantitative comparisons between experimental and computational results show promise for applying a FEM to the flow forming process. Experimental failure locations correlated well with computed damage gradients, and there was good agreement between the measured and predicted roller forces during forming.

Elizabeth Urig↗
NASA NTRS1990

Numerical modeling of flows in simulated brush seal configurations

The two-dimensional flow field in a linear brush seal configuration is modeled and the flow field determined numerically using an FEM model consisting of a matrix of 6 x 6 uniformly spaced circular bristles with a uniform flow field upstream of the simulated brush element. Results from the numerical model are compared to earlier work for a circular seal of similar configuration. Temporal disturbance of the upstream flow field is studied to elucidate system dynamics; however, only the local velocity field is investigated. The effect of tangential flow at the upstream boundary is also studied. Tangential flow produces changes in the pressure applied to the bristle but has little effect on the axial flow rates.

Mullen, R. L.↗
NASA NTRS1985

Probabilistic finite elements for transient analysis in nonlinear continua

The probabilistic finite element method (PFEM), which is a combination of finite element methods and second-moment analysis, is formulated for linear and nonlinear continua with inhomogeneous random fields. Analogous to the discretization of the displacement field in finite element methods, the random field is also discretized. The formulation is simplified by transforming the correlated variables to a set of uncorrelated variables through an eigenvalue orthogonalization. Furthermore, it is shown that a reduced set of the uncorrelated variables is sufficient for the second-moment analysis. Based on the linear formulation of the PFEM, the method is then extended to transient analysis in nonlinear continua. The accuracy and efficiency of the method is demonstrated by application to a one-dimensional, elastic/plastic wave propagation problem. The moments calculated compare favorably with those obtained by Monte Carlo simulation. Also, the procedure is amenable to implementation in deterministic FEM based computer programs.

Liu, W. K.↗
NASA NTRS2016

Suited and Unsuited Hybrid III Impact Testing and Finite Element Model Characterization

NASA spacecraft design requirements for occupant protection are a combination of the Brinkley Dynamic Response Criteria and injury assessment reference values (IARV) extracted from anthropomorphic test devices (ATD). For the ATD IARVs, the requirements specify the use of the 5th percentile female Hybrid III and the 95th percentile male Hybrid III. Each of these ATDs is required to be fitted with an articulating pelvis (also known as the aerospace pelvis) and a straight spine. The articulating pelvis is necessary for the ATD to fit into spacecraft seats, while the straight spine is required as injury metrics for vertical accelerations are better defined for this configuration. Sled testing of the Hybrid III 5th Percentile Female Anthropomorphic Test Device (ATD) was performed at Wright-Patterson Air Force Base (WAPFB). Two 5th Percentile ATDs were tested, the Air Force Research Lab (AFRL) and NASA owned Hybrid III ATDs with aerospace pelvises. Testing was also conducted with a NASA-owned 95th Percentile Male Hybrid III with aerospace pelvis at WPAFB. Testing was performed using an Orion seat prototype provided by Johnson Space Center (JSC). A 5-point harness comprised of 2 inch webbing was also provided by JSC. For suited runs, a small and extra-large Advanced Crew Escape System (ACES) suit and helmet were also provided by JSC. Impact vectors were combined frontal/spinal and rear/lateral. Some pure spinal and rear axis testing was also performed for model validation. Peak accelerations ranged between 15 and 20-g. This range was targeted because the ATD responses fell close to the IARV defined in the Human-Systems Integration Requirements (HSIR) document. Rise times varied between 70 and 110 ms to assess differences in ATD responses and model correlation for different impact energies. The purpose of the test series was to evaluate the Hybrid III ATD models in Orion-specific landing orientations both with and without a spacesuit. The results of these tests were used by the NASA Engineering and Safety Center (NESC) to validate the finite element model (FEM) of the Hybrid III 5th Percentile Female ATD. Physical test data was compared to analytical predictions from simulations, and modelling uncertainty factors have been determined for each injury metric. Additionally, the test data has been used to further improve the FEM, particularly in the areas of the ATD preload, harness, and suit and helmet effects.

Lawrence, C.↗
DOE OSTI2020

Improving Computational Efficiency of Mechanical Finite Element Method Simulations for PV Modules: Preprint

In this work we have elucidated the trade off between structural mechanics FEM model accuracy and computation time by employing lower fidelity viscoelastic models. Results indicate that computation time can easily be cut in half while only expecting a potential maximum error of 10 % by considering lower fidelity models. A novel approach to produce the Prony Series fit for viscoelastic characterization has also been presented. Employing this approach, we were able to achieve a further 10 % reduction in computation time without further sacrificing simulation accuracy.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2023

Micro-beam bending of FCC bicrystals: A comparison between defect dynamics simulations and experiments

To understand the role of the grain boundary (GB) in plasticity at small scale, a concurrently coupled mesoscale plasticity model was developed to simulate micro-bending of bicrystalline micron-sized beams. By coupling dislocation dynamics (DD) with a finite element model (FEM), a novel defect dynamics model provides the means to investigate intricate interactions between dislocations and GBs under various loading conditions. Our simulations of micro-bending agree well with corresponding micro-bending experiments, and they show that mechanical response of bicrystals could have not only hardening but also softening depending on the characters of the GB. In addition, changing the location of the GB in the microbeams results in different mechanical responses; GBs located at the neutral plane show softening compared to single crystals, while inclined GBs located halfway along the length of the beam show little effect. In conclusion, simulation results could provide a clear picture on detailed dislocation-GB interactions, and quantitative resolved shear stress analysis supplemented by dislocation density distribution is used to analyze the mechanical response of bicrystalline samples.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Numerical Simulation of Irradiation Induced Swelling for STAR 4.1 Blanket

Here, this article presents the preliminary results of FEM implementation of irradiation-induced swelling, hardening, and creep effect using ANSYS user programmable features (UPFs) applied to material properties. The first liquid metal breeding blanket model for STAR 4.1 tokamak as a part of the virtual prototyping system is under development at Princeton Plasma Physics Laboratory (PPPL). Intense neutron irradiation produces significant changes in the physical and mechanical properties of Fe-(8%–9%)Cr-based reduced activation ferritic martensitic (RAFM) steels. An ANSYS model of the breeding blanket was built to involve these irradiation-induced material property changes with neutron fluence distribution mapped from the MCNP model. Simulation results show that significant structure deformation forms from nuclear swelling, stress, and plastic strain arising from DPA gradient through wall thickness. The material degradation effect is not negligible. This modeling feature can either assist in analyzing the structural behaviors with the influence of nuclear swelling or provide guidance to design the structure to withstand irradiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗
NASA NTRS1997

Simulating and Synthesizing Substructures Using Neural Network and Genetic Algorithms

The feasibility of simulating and synthesizing substructures by computational neural network models is illustrated by investigating a statically indeterminate beam, using both a 1-D and a 2-D plane stress modelling. The beam can be decomposed into two cantilevers with free-end loads. By training neural networks to simulate the cantilever responses to different loads, the original beam problem can be solved as a match-up between two subsystems under compatible interface conditions. The genetic algorithms are successfully used to solve the match-up problem. Simulated results are found in good agreement with the analytical or FEM solutions.

Liu, Youhua↗
DOE OSTI2021

MRP and NEAMS Meeting - Current Modeling to Support Experiments

Heat-pipe cooled micro-reactor is currently being actively investigated in the NS&T department at INL. Given the high operation temperature and large temperature gradients in the system, the induced thermal stresses need to be well analyzed and kept within the safety margin. Two commercial FEM software (ABAQUS and STAR-CCM+) have been proposed to simulate and compare for both the temperature and stress distributions for the single heat pipe experiments, with the middle single heat pipe and 6 heaters surrounded. This presentation file is aimed to provide the latest updates for the researchers who is under NEAMS program and involved with the DOE micro-reactor program

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Sensor Development for Liquid Water Detection in Dry Storage Casks: FY 2023 Status

Modeling efforts were undertaken in fiscal year (FY) 2023 to evaluate the feasibility and capability for sensing of liquid water inside canister-based dry cask storage systems (DCSSs). The focus was on the development of full-scale finite element models (FEMs) of ultrasound propagation through DCSS canister components. For these initial investigations, simulation of the baseplate component was targeted, envisioning water collection at the bottom of a vertically oriented canister. The effects of internal canister components (i.e., fuel basket, fuel assemblies) on ultrasound propagation in the baseplate component were evaluated in these efforts. The environment inside a DCSS confinement is intended to be inert and free of water to prevent potential corrosion of used fuel cladding or other internal hardware. Spent fuel assemblies are dried, after storage in water pools, to make sure water has been removed from assembly cavities. However, there is some uncertainty about the amount of residual water potentially left behind in a DCSS after drying processes, because water can become trapped in cavities or other small crevices in the surfaces formed by the fuel cladding, fuel assemblies, and other internal hardware components. Considering the complex space- and time-dependent temperature profiles in DCSSs, water may be in a liquid or gas phase depending on its location in the cask and how long the cask has been in storage. Evacuating most water and oxidizing agents contained within a canister is recommended by NUREG-1536 (NRC 2010), which covers DCSSs. As summarized by Salazar et al. (2020), existing guidance typically relies on achievable vacuum pumping pressures sustained over a hold period as a signal of dryness and water removal. However, time series data about water removal from full-scale commercial drying procedures are lacking (Hanson and Alsaed 2019). A review of drying specifications from several vendors led to the conclusion that if the specifications are followed correctly, the residual moisture in DCSSs should present an insignificant risk of cladding degradation (Knoll and Gilbert 1987). A more recent analysis concluded that much larger quantities of residual water could remain in DCSSs, but the amount would still not be expected to lead to significant corrosion of fuel cladding or other internal components (Jung et al. 2013). Industry drying procedures are mostly prescriptive in nature, and operational issues arising during the process could result in incomplete drying. A summary of operational issues and potential negative effects of residual water is provided by Salazar et al. (2019). Experimental efforts are ongoing to validate the extent of water removal in a DCSS based on drying procedures used at nuclear power plants through well-designed investigations of drying process efficacy and water retention (Durbin et al. 2021; Pulido et al. 2022a, Pulido et al. 2022b). This has been approached by simulating limited portions of a DCSS internal volume and fuel assemblies. So far, these efforts have focused on the effects of potential water trapping in the dashpot region of control rod guide tubes but are expected to be expanded to include the effects of other internal hardware features and failed fuel rod cladding. A method for detection and measurement of liquid water (Meyer et al. 2022; Meyer et al. 2021), in tandem with the drying information gained through experimental investigations, provides comprehensive bases for understanding the internal conditions of DCSSs and provides technical information to support licensing decisions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗