Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element Method”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Determination of Strain Rate Sensitivity of Micro-struts Manufactured Using the Selective Laser Melting Method

Micro-lattice structures manufactured using the selective laser melting (SLM) process provides the opportunity to realize optimal cellular materials for impact energy absorption. In this paper, strain rate-dependent material properties are measured for stainless steel 316L SLM micro-lattice struts in the strain rate range of 10{sup −3} to 6000 s{sup −1}. At high strain rates, a novel version of the split Hopkinson Bar has been developed. Strain rate-dependent materials data have been used in Cowper–Symonds material model, and the scope and limit of this model in the context of SLM struts have been discussed. Strain rate material data and the Cowper–Symonds model have been applied to the finite element analysis of a micro-lattice block subjected to drop weight impact loading. The model output has been compared to experimental results, and it has been shown that the increase in crush stress due to impact loading is mainly the result of strain rate material behavior. Hence, a systematic methodology has been developed to investigate the impact energy absorption of a micro-lattice structure manufactured using additive layer manufacture (SLM). This methodology can be extended to other micro-lattice materials and configurations, and to other impact conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2016

WE-AB-202-05: Validation of Lung Stress Maps for CT-Ventilation Imaging

Purpose: To date, lung CT-ventilation imaging has been based on quantification of local breathing-induced changes in Hounsfield Units (HU) or volume. This work investigates the use of a stress map resulting from a biomechanical deformable image registration (DIR) algorithm as a metric of the ventilation function. Method: Eight lung cancer patients presenting different kinds of ventilation defects were retrospectively analyzed. Additionally, to the 4DCT acquired for radiotherapy planning, five of them had PET and three had SPECT imaging following inhalation of Ga-68 and Tc-99m, respectively. For each patient, the inhale phase of the 4DCT was registered to the exhale phase using Morfeus, a biomechanical DIR algorithm based on the determination of boundary conditions on the lung surfaces and vessel tree. To take into account the heterogeneity of the tissue stiffness in the stress map estimation, each tetrahedral element of the finite-element model was assigned a Young’s modulus ranging from 60kPa to 12MPa, as a function of the HU in the inhale CT. The node displacements and element stresses resulting from the numerical simulation were used to generate three CT-ventilation maps based on: (i) volume changes (Jacobian determinant), (ii) changes in HU, (iii) the maximum principal stress. The voxel-wise correlation between each CT-ventilation map and the PET or SPECT V image was computed in a lung mask. Results: For patients with PET, the mean (min-max) Spearman correlation coefficients r were: 0.33 (0.19–0.45), 0.36 (0.16–0.51) and 0.42 (0.21–0.59) considering the Jacobian, changes in HU and maximum principal stress, respectively. For patients with SPECT V, the mean r were: 0.12 (−0.12–0.43), 0.29 (0.22–0.45) and 0.33 (0.25–0.39). Conclusion: The maximum principal stress maps showed a stronger correlation with the ventilation images than the previously proposed Jacobian or change in HU maps. This metric thus appears promising for CT-ventilation imaging. This work was funded in part by NIH P01CA059827.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Review on failure prediction techniques of composite single lap joint

Adhesive bonding is the most appropriate joining method in construction of composite structures. The use of reliable design and prediction technique will produce better performance of bonded joints. Several papers from recent papers and journals have been reviewed and synthesized to understand the current state of the art in this area. It is done by studying the most relevant analytical solutions for composite adherends with start of reviewing the most fundamental ones involving beam/plate theory. It is then extended to review single lap joint non linearity and failure prediction and finally on the failure prediction on composite single lap joint. The review also encompasses the finite element modelling part as tool to predict the elastic response of composite single lap joint and failure prediction numerically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Multifidelity uncertainty quantification with models based on dissimilar parameters

Multifidelity uncertainty quantification (MF UQ) sampling approaches have been shown to significantly reduce the variance of statistical estimators while preserving the bias of the highest-fidelity model, provided that the low-fidelity models are well correlated. However, maintaining a high level of correlation can be challenging, especially when models depend on different input uncertain parameters, which drastically reduces the correlation. Existing MF UQ approaches do not adequately address this issue. In this work, we propose a new sampling strategy that exploits a shared space to improve the correlation among models with dissimilar parameterization. We achieve this by transforming the original coordinates onto an auxiliary manifold using the adaptive basis (AB) method (Tipireddy and Ghanem, 2014). The AB method has two main benefits: (1) it provides an effective tool to identify the low-dimensional manifold on which each model can be represented, and (2) it enables easy transformation of polynomial chaos representations from high- to low-dimensional spaces. This latter feature is used to identify a shared manifold among models without requiring additional evaluations. Here we present two algorithmic flavors of the new estimator to cover different analysis scenarios, including those with legacy and non-legacy high-fidelity (HF) data. We provide numerical results for analytical examples, a direct field acoustic test, and a finite element model of a nuclear fuel assembly. For all examples, we compare the proposed strategy against both single-fidelity and MF estimators based on the original model parameterization.

42 ENGINEERING↗
DOE OSTI2018

Finite Element Based Full-Life Cyclic Stress Analysis of 316 Grade Nuclear Reactor Stainless Steel Under Constant, Variable, and Random Fatigue Loading

Although S~N curve-based approaches are widely followed for fatigue evaluation of nuclear reactor components and other safety critical structural systems, there is a chance of large uncertainty in estimated fatigue lives. This uncertainty may be reduced by using a more mechanistic approach such as physics based three-dimensional (3D) finite element (FE) methods. In a recent paper (Barua et al., 2018, ASME J. Pressure Vessel Technol., 140(1), p. 011403), a fully mechanistic fatigue modeling approach which is based on time-dependent stress–strain evolution of material over the entire fatigue life was presented. Based on this approach, FE-based cyclic stress analysis was performed on 316 nuclear grade reactor stainless steel (SS) fatigue specimens, subjected to constant, variable, and random amplitude loading, for their entire fatigue lives. The simulated results are found to be in good agreement with experimental observation. An elastic-plastic analysis of a pressurized water reactor (PWR) surge line (SL) pipe under idealistic fatigue loading condition was performed and compared with experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dyn$\mathrm{AMO}$: Multi-agent reinforcement learning for dynamic anticipatory mesh optimization with applications to hyperbolic conservation laws

Here we introduce DynAMO, a reinforcement learning paradigm for Dynamic Anticipatory Mesh Optimization. Adaptive mesh refinement is an effective tool for optimizing computational cost and solution accuracy in numerical methods for partial differential equations. However, traditional adaptive mesh refinement approaches for time-dependent problems typically rely only on instantaneous error indicators to guide adaptivity. As a result, standard strategies often require frequent remeshing to maintain accuracy. In the DynAMO approach, multi-agent reinforcement learning is used to discover new local refinement policies that can anticipate and respond to future solution states by producing meshes that deliver more accurate solutions for longer time intervals. By applying DynAMO to discontinuous Galerkin methods for the linear advection and compressible Euler equations in two dimensions, we demonstrate that this new mesh refinement paradigm can outperform conventional threshold-based strategies while also generalizing to different mesh sizes, remeshing and simulation times, and initial conditions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Characterization of Flame Cut Heavy Steel: Modeling of Temperature History and Residual Stress Formation

Heavy steel plates are used in demanding applications that require both high strength and hardness. An important step in the production of such components is cutting the plates with a cost-effective thermal cutting method such as flame cutting. Flame cutting is performed with a controlled flame and oxygen jet, which burns the steel and forms a cutting edge. However, the thermal cutting of heavy steel plates causes several problems. A heat-affected zone (HAZ) is generated at the cut edge due to the steep temperature gradient. Consequently, volume changes, hardness variations, and microstructural changes occur in the HAZ. In addition, residual stresses are formed at the cut edge during the process. In the worst case, unsuitable flame cutting practices generate cracks at the cut edge. The flame cutting of thick steel plate was modeled using the commercial finite element software ABAQUS. The results of modeling were verified by X-ray diffraction-based residual stress measurements and microstructural analysis. The model provides several outcomes, such as obtaining more information related to the formation of residual stresses and the temperature history during the flame cutting process. In addition, an extensive series of flame cut samples was designed with the assistance of the model.

36 MATERIALS SCIENCE↗
DOE OSTI2021

On the competitive antagonism effect in combined chloride-sulfate attack: A numerical exploration

Highlights: • A chemo-physical-mechanical method based on thermodynamic modelling is implemented. • The contradicting findings from reported experiments are numerically explored. • The competitive antagonism effect is time-dependent rather than a conclusive term. • The competitive antagonism effect is quantifiable with a novel assessment criterion. It is of tremendous difficulty to draw a conclusive verdict on the complex process of mutual interference regarding the combined chloride-sulfate attack on concrete, where many have explored the topic, but contradicting results have been reported. In this paper, by taking advantages of the chemo-physical-mechanical method, based on the robust technique of thermodynamic modelling, the competitive antagonism effect in combined chloride-sulfate attack is investigated. Firstly, a range of reported tests are numerically modelled to demonstrate the effectiveness of the employed method in interpreting those opposing and diverse experimental outcomes. More importantly, a fresh conjecture incorporating a novel assessment criterion for the competitive antagonism effect is proposed and supported through a carefully designed illustrative example. Based on the numerical exploration, unique findings and implications for future engineering practices are revealed.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Predicting Metallic Defects in Casting BrA9ZhZL Alloy Tool Electrodes

A method is proposed for modeling the lost-wax casting of tool electrodes by means of fusible BrA9ZhZL copper-based alloy. The functionality of the casting system is analyzed. The heating and cooling of a ceramic mold is taken into account in visualizing the temperature fields. Sections where shrinkage pores may form in the casting are identified.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Space-Time Block Preconditioning for Incompressible Flow

Parallel-in-time methods have become increasingly popular in the simulation of time-dependent numerical PDEs, allowing for the efficient use of additional message passing interface processes when spatial parallelism saturates. Most methods treat the solution and parallelism in space and time separately. In contrast, all-at-once methods solve the full space-time system directly, largely treating time as simply another spatial dimension. All-at-once methods offer a number of benefits over separate treatment of space and time, most notably significantly increased parallelism and faster time to solution (when applicable). However, the development of fast, scalable all-at-once methods has largely been limited to time-dependent (advection-)diffusion problems. This paper introduces the concept of space-time block preconditioning for the all-at-once solution of incompressible flow. By extending well-known concepts of spatial block preconditioning to the space-time setting, we develop a block preconditioner whose application requires the solution of a space-time (advection-)diffusion equation in the velocity block, coupled with a pressure Schur complement approximation consisting of independent spatial solves at each time-step, and a space-time matrix-vector multiplication. The new method is tested on four classical models in incompressible flow. Finally, the results indicate perfect scalability in refinement of spatial and temporal mesh spacing, perfect scalability in nonlinear Picard iteration count when applied to a nonlinear Navier--Stokes problem, and minimal overhead in terms of number of preconditioner applications compared with sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Pronghorn Theory Manual

This manual is intended to provide a thorough description of the governing equations; the relations used for material, heat transfer, and fluid flow parameters; and the numerical method utilized in Pronghorn. Detailed derivations of the governing equations are given, beginning from first principles in order to be as explicit as possible about assumptions made along the way and guide any future methods extensions. Additional explanatory or tangential information that is not required for detailed understanding of the governing equations used in Pronghorn is shown in gray boxes in the text. The convention used in this document is to represent all class names (or file names) in typewriter font.

42 ENGINEERING↗
DOE OSTI2026

Polynomial range estimation as a troubled-cell indicator for high-order methods

Two troubled-cell indicators based on polynomial range estimation methods are used to flag cells that may violate positivity constraints. One method uses interval extension, and the second uses the range enclosure property of the Bernstein polynomial basis. Furthermore, both methods reduce compute time for the positivity preserver by limiting its application to a subset of cells. The Bernstein polynomial method remains effective as the problem dimensionality increases. Interval extension applied to the internal energy equation permits the use of the troubled-cell indicators for rational functions, though performance suffers compared to directly applying the indicators to polynomial functions.

42 ENGINEERING↗
DOE OSTI2021

Monitoring of corrosion effects in pipes with multi-mode acoustic signals

Current commercial long-range guided wave inspection and monitoring technologies for pipes usually require establishment of certain specific acoustic modes for wave propagation and use a large number of sensors in a certain geometric arrangement (for example, uniformly-spaced arrays of sensors on a collar, placed around a pipe). This paper, in contrast, presents a long-range corrosion detection method in pipes based on multi-mode guided acoustic waves, generated and received using simple, single-element sensors. In this technique, available guided wave modes in any interrogated zone of the structure are automatically down-selected by the medium, and changes in energy distribution over the selected modes are analyzed through signal processing to monitor material loss or addition. Here, to investigate the feasibility of the proposed method, we have performed both experimental and numerical simulation studies on 50.8-mm (2-in)-diameter schedule-40 carbon steel pipes. In these studies, a linear chirp signal is injected in compression mode, perpendicular to the surface of the pipe, generating multi-mode, multi-directional guided waves that travel within the pipe wall. Corrosion effects in the pipes are simulated numerically and experimentally. To reduce errors in the measurement technique due to changes in the local or global temperature of the measured media, a universal temperature compensation algorithm is introduced. The experimental results show good agreement with simulation results, indicating that the proposed multi-mode wave generation and detection in a pipe could be used for long-range detection and monitoring of corrosion effects, such as pipe wall losses or solid material deposits.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Contact resonances of U-shaped atomic force microscope probes

Recent approaches used to characterize the elastic or viscoelastic properties of materials with nanoscale resolution have focused on the contact resonances of atomic force microscope (CR-AFM) probes. The experiments for these CR-AFM methods involve measurement of several contact resonances from which the resonant frequency and peak width are found. The contact resonance values are then compared with the noncontact values in order for the sample properties to be evaluated. The data analysis requires vibration models associated with the probe during contact in order for the beam response to be deconvolved from the measured spectra. To date, the majority of CR-AFM research has used rectangular probes that have a relatively simple vibration response. Recently, U-shaped AFM probes have created much interest because they allow local sample heating. However, the vibration response of these probes is much more complex such that CR-AFM is still in its infancy. In this article, a simplified analytical model of U-shaped probes is evaluated for contact resonance applications relative to a more complex finite element (FE) computational model. The tip-sample contact is modeled using three orthogonal Kelvin-Voigt elements such that the resonant frequency and peak width of each mode are functions of the contact conditions. For the purely elastic case, the frequency results of the simple model are within 8% of the FE model for the lowest six modes over a wide range of contact stiffness values. Results for the viscoelastic contact problem for which the quality factor of the lowest six modes is compared show agreement to within 13%. These results suggest that this simple model can be used effectively to evaluate CR-AFM experimental results during AFM scanning such that quantitative mapping of viscoelastic properties may be possible using U-shaped probes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

FY 2016 Status Report: Documentation of All CIRFT Data including Hydride Reorientation Tests (Draft M2)

The first portion of this report provides a detailed description of fiscal year (FY) 2015 test result corrections and analysis updates based on FY 2016 updates to the Cyclic Integrated Reversible-Bending Fatigue Tester (CIRFT) program methodology, which is used to evaluate the vibration integrity of spent nuclear fuel (SNF) under normal conditions of transport (NCT). The CIRFT consists of a U-frame test setup and a real-time curvature measurement method. The three-component U-frame setup of the CIRFT has two rigid arms and linkages connecting to a universal testing machine. The curvature SNF rod bending is obtained through a three-point deflection measurement method. Three linear variable differential transformers (LVDTs) are clamped to the side connecting plates of the U-frame and used to capture deformation of the rod. The second portion of this report provides the latest CIRFT data, including data for the hydride reorientation test. The variations in fatigue life are provided in terms of moment, equivalent stress, curvature, and equivalent strain for the tested SNFs. The equivalent stress plot collapsed the data points from all of the SNF samples into a single zone. A detailed examination revealed that, at the same stress level, fatigue lives display a descending order as follows: H. B. Robinson Nuclear Power Station (HBR), LMK, and mixed uranium-plutonium oxide (MOX). Just looking at the strain, LMK fuel has a slightly longer fatigue life than HBR fuel, but the difference is subtle. The third portion of this report provides finite element analysis (FEA) dynamic deformation simulation of SNF assemblies . In a horizontal layout under NCT, the fuel assembly’s skeleton, which is formed by guide tubes and spacer grids, is the primary load bearing apparatus carrying and transferring vibration loads within an SNF assembly. These vibration loads include interaction forces between the SNF assembly and the canister basket walls. Therefore, the integrity of the guide tubes and spacer grids critically affects the vibration intensity of the fuel assembly during transport and must be considered when developing the multipurpose purpose canister (MPC) design for safe SNF transport.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

IMEX and exact sequence discretization of the multi-fluid plasma model

Multi-fluid plasma models, where an electron fluid is modeled in addition to multiple ion and neutral species as well as the full set of Maxwell's equations, are useful for representing physics beyond the scope of classic MHD. This advantage presents challenges in appropriately dealing with electron dynamics and electromagnetic behavior characterized by the plasma and cyclotron frequencies and the speed of light. For physical systems, such as those near the MHD asymptotic regime, this requirement drastically increases runtimes for explicit time integration even though resolving fast dynamics may not be critical for accuracy. Implicit time integration methods, with efficient solvers, can help to step over fast time-scales that constrain stability, but do not strongly influence accuracy. As an extension, Implicit-explicit (IMEX) schemes provide an additional mechanism to choose which dynamics are evolved using an expensive implicit solve or resolved using a fast explicit solve. In this study, in addition to IMEX methods we also consider a physics compatible exact sequence spatial discretization. Here, this combines nodal bases (H-Grad) for fluid dynamics with a set of vector bases (H-Curl and H-Div) for Maxwell's equations. This discretization allows for multi-fluid plasma modeling without violating Gauss' laws for the electric and magnetic fields. This initial study presents a discussion of the major elements of this formulation and focuses on demonstrating accuracy in the linear wave regime and in the MHD limit for both a visco-resistive and a dispersive ideal MHD problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Large Area Detection of Microstructural Defects with Multi-Mode Ultrasonic Signals

Cyclic loading or other stresses can lead to development of cracks and crack growth in mechanical structures, leading to eventual failure. While ultrasound imaging can be used for non-destructive testing of such structures, conventional ultrasound techniques are often limited by crack size, density, and areal coverage. An effective characterization of real-world, large-area structures is required at an early damage stage to prevent catastrophic failure and predict remaining life. In this study, a new nonlinear ultrasonic testing (NUT) method is proposed for large-area monitoring of practical structures with arbitrary complexity by using multiple-mode guided-wave ultrasonic signals. The proposed guided-wave NUT technique requires single-element transducers, simple electronics, and a mixed time-frequency domain signal processing. As a proof-of-concept demonstration, numerical simulations and experiments are performed on an A36 carbon steel beam assembly with previously formed microstructural defects that cause nonlinearities in ultrasonic response. The quadratic dependence of the nonlinear wave excitation on the input ultrasonic signal amplitude is shown by numerical simulations, and such a nonlinear ultrasonic response is experimentally observed in the zone with a high density of microstructural defects.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Multiscale modeling of twinning and detwinning behavior of HCP polycrystals

A multiscale scheme is presented here to capture the twinning and detwinning mechanisms during cyclic loading of HCP polycrystals. In the current framework, the crystal plasticity finite element (CPFE) method is employed at the microstructural scale. The response of twinned regions at a material point is calculated using a novel homogenization scheme via a Taylor-type model. The homogenized behavior is then used in the CPFE to obtain the response of a polycrystal. An efficient rate-independent crystal plasticity model that has been incorporated into the open source PRISMS-Plasticity software is introduced here to capture the response of the twinned region. Accordingly, the corresponding algorithmic tangent modulus is derived. Finally, this multiscale scheme is validated against data available for the cyclic response of Mg alloy ZK60A. The results show that the multiscale model can capture the key elements of the sample response observed during the experiment.

36 MATERIALS SCIENCE↗