Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite Element”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Finite Element Modeling of Spent Nuclear Fuel in Dry Storage under Seismic Excitation

Spent nuclear fuel (SNF) is currently stored at independent spent fuel storage installations (ISFSIs) at nuclear power plant sites across the US. The US Department of Energy (DOE) Spent Fuel and Waste Science and Technology (SFWST) program is currently working to understand the response of SNF during a seismic event. This investigation is part of a long-term effort by DOE to quantify the external mechanical loads on SNF. The SFWST program has recently completed campaigns evaluating SNF behavior during normal conditions of transport and in 30 cm package drop scenarios; both campaigns included a combination of testing and modeling. Likewise, the investigation into seismic response of SNF will include full-size shake table tests of a vertical storage overpack containing instrumented surrogate fuel assemblies, as well as pre- and post-test finite element modeling. This work focuses on the development of finite element models simulating the planned test configuration. The seismic response is calculated at two scales: at the package level, and at the fuel level. The package level model evaluates the behavior of the storage overpack and its contents in response to the ground motions of a representative earthquake. The motion of the package is input to the fuel level model, which consists of a detailed finite element model of a 17x17 or 16x16 pressurized water reactor (PWR) fuel assembly or of a 10x10 boiling water reactor (BWR) fuel assembly. The stress and strain on individual components such as fuel rods and guide tubes are extracted from this model. This work details the SNF response to various input ground motions, considering representative earthquakes in the Central and Eastern US under various site conditions. Also discussed is the development of a “reduced scope” finite element model of the SNF, which consists of a single fuel rod while retaining the dynamic response of the entire fuel assembly. This reduced scope model addresses the significant computational expense of modeling the SNF under long-duration input ground motions and facilitates simulating additional parameter variations of the model beyond the as-tested configuration. The modeling results provide important insight that will help to inform future testing and modeling efforts by the SFWST program to close the knowledge gap related to external loads on SNF.

Kadooka, Kevin↗
NASA NTRS1992

Progress on hybrid finite element methods for scattering by bodies of revolution

Progress on the development and implementation of hybrid finite element methods for scattering by bodies of revolution are described. It was found that earlier finite element-boundary integral formulations suffered from convergence difficulties when applied to large and thin bodies of revolution. An alternative implementation is described where the finite element method is terminated with an absorbing termination boundary. In addition, an alternative finite element-boundary integral implementation is discussed for improving the convergence of the original code.

Collins, Jeffery D.↗
NASA NTRS1985

Implementation of design sensitivity analysis with existing finite element codes

A numerical method is presented to implement structural design sensitivity analysis theory, using the versatility and convenience of existing finite element structural analysis programs. Design variables such as thickness and cross-sectional areas of components of individual members and built-up structures are considered. Structural performance functionals considered include displacement and stress. The method is also applicable for eigenvalue problem design sensitivity analysis. It is shown that calculations can be carried out outside existing finite element codes, using postprocessing data only. Thus, design sensitivity analysis software does not have to be embedded in an existing finite element code. Feasability of the method is shown through analysis of several problems, including a built-up structure. Accurate design sensitivity results are obtained without the uncertainty of numerical accuracy associated with selection of finite difference perturbations.

Choi, Kyung K.↗
NASA NTRS1986

A mesh re-zoning technique for finite element simulations of metal forming processes

Based on some fundamental properties of finite element approximations, a mesh re-zoning scheme is proposed for finite element simulations of metal forming problems. It is demonstrated that this technique is indispensable in analyzing many difficult forming processes, especially when there exist corners or very irregular shapes on the boundaries. The algorithm is tested by a backward extrusion process and direct extrusion through a square die.

Cheng, J.-C.↗
NASA NTRS1987

Workshop on the Integration of Finite Element Modeling with Geometric Modeling

The workshop on the Integration of Finite Element Modeling with Geometric Modeling was held on 12 May 1987. It was held to discuss the geometric modeling requirements of the finite element modeling process and to better understand the technical aspects of the integration of these two areas. The 11 papers are presented except for one for which only the abstract is given.

Wozny, Michael J.↗
DOE OSTI2021

Modeling of High-Temperature Corrosion of Zirconium Alloys Using the eXtended Finite Element Method (X-FEM)

Oxidation modeling in modern nuclear fuel performance codes is currently limited by the lack of coupling with mechanics, thus preventing proper description of how high-temperature oxidation impacts mechanical properties. This is mostly due to the fact that the finite difference formalism adopted in corrosion models is incompatible with the direct coupling with mechanics in the finite element modeling employed in modern nuclear fuel performance codes. In this study, a physically based zirconium alloy corrosion model called the Coupled-Current Charge Compensation (C4) model, which was initially developed for operating temperature conditions, has been updated to include high-temperature corrosion in order to provide additional critical information (e.g., oxygen concentration profile) under loss-of-coolant accident (LOCA) conditions—information lacking in existing empirical models. The C4 model was implemented in the MOOSE finite-element framework developed at Idaho National Laboratory, enabling it to be used in the BISON nuclear fuel performance code based on the MOOSE framework. To precisely track the different interfaces at a relatively low computational cost, the eXtended Finite Element Method (X-FEM) was applied in MOOSE. The model’s results were compared to those of existing empirical models as well as metallographic analysis of high-temperature oxidized Zircaloy-4 coupons. Oxygen diffusivities in the a and ß phases resulting from this comparison closely agree with those found in the literature. The C4 model implemented with X-FEM in MOOSE now has the capability to accurately predict oxide, oxygen-stabilized a, and prior ß phase layer growth kinetics under isothermal exposure at high temperature (1000–1500°C). Furthermore, in contrast with the empirical models, the C4 model accounts for the finite thickness of the fuel cladding. It can predict the oxygen concentration profile evolution through the whole cladding, enabling evaluation of the remaining ductile thickness—a crucial variable for modeling the mechanical behavior of the fuel cladding under LOCA. Furthermore, this implementation allows direct coupling with mechanics, at a low computing cost, using finite-element-based nuclear fuel performance codes such as BISON.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Finite-element analyses and fracture simulation in thin-sheet aluminum alloy

A two-dimensional, elastic-plastic finite-element analysis was used with a critical crack-tip-opening angle (CTOA) fracture criterion to model stable crack growth in thin-sheet 2024-T3 aluminum alloy under monotonic loading after precracking at different cyclic stress levels. Tests were conducted on three types of specimens: middle-crack, three-hole-crack and blunt-notch tensile specimens. An experiment technique was developed to measure CTOA during crack growth initiation and stable tearing using a high-resolution video camera and recorder. Crack front shapes were also measured during initiation and stable tearing using a fatigue marker-load technique. Three-dimensional elastic-plastic finite-element analyses of these crack shapes for stationary cracks were conducted to study the crack-front opening displacements. Predicted load against crack extension on middle-crack tension specimens agreed well with test results even for large-scale plastic deformations. The analyses were able to predict the effects of specimen size and precracking stress history on stable tearing. Predicted load against load-line displacements agreed well with test results up to maximum load bu the analyses tended to overpredict displacements as crack grew beyond the maximum load under displacement-controlled conditions. During the initiation phase, the measured CTOA values were high but decreased and remained nearly constant after a small amount of stable tearing. The constant value of CTOA agree well with the calculated value from the finite-element analysis. The larger CTOA values measured at the sheet surface during the initiation phase may be associated with the crack tunneling observed in the tests. Three-dimensional analyses for nonstraight crack fronts predicted much higher displacements near the free surface than in the interior.

Newman, J. C., Jr.↗
NASA NTRS1989

On finite element implementation and computational techniques for constitutive modeling of high temperature composites

The research work performed during the past year on finite element implementation and computational techniques pertaining to high temperature composites is outlined. In the present research, two main issues are addressed: efficient geometric modeling of composite structures and expedient numerical integration techniques dealing with constitutive rate equations. In the first issue, mixed finite elements for modeling laminated plates and shells were examined in terms of numerical accuracy, locking property and computational efficiency. Element applications include (currently available) linearly elastic analysis and future extension to material nonlinearity for damage predictions and large deformations. On the material level, various integration methods to integrate nonlinear constitutive rate equations for finite element implementation were studied. These include explicit, implicit and automatic subincrementing schemes. In all cases, examples are included to illustrate the numerical characteristics of various methods that were considered.

Saleeb, A. F.↗
NASA NTRS1989

The MHOST finite element program: 3-D inelastic analysis methods for hot section components. Volume 2: User's manual

The user options available for running the MHOST finite element analysis package is described. MHOST is a solid and structural analysis program based on the mixed finite element technology, and is specifically designed for 3-D inelastic analysis. A family of 2- and 3-D continuum elements along with beam and shell structural elements can be utilized, many options are available in the constitutive equation library, the solution algorithms and the analysis capabilities. The outline of solution algorithms is discussed along with the data input and output, analysis options including the user subroutines and the definition of the finite elements implemented in the program package.

Nakazawa, Shohei↗
DOE OSTI2021

Efficient quadrature rules for finite element discretizations of nonlocal equations

In this paper we design efficient quadrature rules for finite element discretizations of nonlocal diffusion problems with compactly supported kernel functions. Two of the main challenges in nonlocal modeling and simulations are the prohibitive computational cost and the nontrivial implementation of discretization schemes, especially in three-dimensional settings. In this work we circumvent both challenges by introducing a parametrized mollifying function that improves the regularity of the integrand, utilizing an adaptive integration technique, and exploiting parallelization. We first showthat the “mollified” solution converges to the exact one as the mollifying parameter vanishes, then we illustrate the consistency and accuracy of the proposed method on several two- and three-dimensional test cases. Furthermore, we demonstrate the good scaling properties of the parallel implementation of the adaptive algorithm and we compare the proposed method with recently developed techniques for efficient finite element assembly.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Domain decomposition and multiscale mortar mixed finite element methods for linear elasticity with weak stress symmetry

Two non-overlapping domain decomposition methods are presented for the mixed finite element formulation of linear elasticity with weakly enforced stress symmetry. The methods utilize either displacement or normal stress Lagrange multiplier to impose interface continuity of normal stress or displacement, respectively. By eliminating the interior subdomain variables, the global problem is reduced to an interface problem, which is then solved by an iterative procedure. The condition number of the resulting algebraic interface problem is analyzed for both methods. A multiscale mortar mixed finite element method for the problem of interest on non-matching multiblock grids is also studied. It uses a coarse scale mortar finite element space on the non-matching interfaces to approximate the trace of the displacement and impose weakly the continuity of normal stress. A priori error analysis is performed. It is shown that, with appropriate choice of the mortar space, optimal convergence on the fine scale is obtained for the stress, displacement, and rotation, as well as some superconvergence for the displacement. Computational results are presented in confirmation of the theory of all proposed methods.

Mathematics↗
NASA NTRS1988

Element-by-element and implicit-explicit finite element formulations for computational fluid dynamics

Preconditioner algorithms to reduce the computational effort in FEM analyses of large-scale fluid-dynamics problems are presented. A general model problem is constructed on the basis of the convection-diffusion equation and the two-dimensional vorticity/stream-function formulation of the Navier-Stokes equations; this problem is then analyzed using element-by-element, implicit-explicit, and adaptive implicit-explicit approximation schemes. Numerical results for the two-dimensional advection and rigid-body rotation of a cosine hill, flow past a circular cylinder, and driven cavity flow are presented in extensive graphs and shown to be in good agreement with those obtained using implicit methods.

Tezduyar, T. E.↗
DOE OSTI2022

A unified framework of stabilized finite element method for solving the Boltzmann transport equation

This paper presents a unified framework of stabilized finite element method for solving the Boltzmann transport equation. Unlike the traditional Petrov-Galerkin finite element method which modifies the test function to construct the stabilization term, we derive the stabilization methods from the standard Galerkin weak form with Sub-grid scale model. The basic idea of this method is to decompose the unknowns into its numerical solution and residual, with an approximation for the residual and embeds it in the Galerkin weak form to yield a stabilized variational formula. Different approximations of the residual lead to different stabilization methods, all the frequently used stabilized methods, including the Streamline Upwinding Petrov-Galerkin (SUPG) method, Galerkin/Least-Square (GLS) method, and Algebraic Sub-Grid Scale (ASGS) method can be obtained from this framework. The similarities and differences of the different approximations are compared in this paper. The numerical results show that the behaviors of the different methods area similar with the same stabilization parameter, and all these stabilized techniques can obtain a correct and stable solution. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1989

Finite element modeling of frictionally restrained composite interfaces

The use of special interface finite elements to model frictional restraint in composite interfaces is described. These elements simulate Coulomb friction at the interface, and are incorporated into a standard finite element analysis of a two-dimensional isolated fiber pullout test. Various interfacial characteristics, such as the distribution of stresses at the interface, the extent of slip and delamination, load diffusion from fiber to matrix, and the amount of fiber extraction or depression are studied for different friction coefficients. The results are compared to those obtained analytically using a singular integral equation approach, and those obtained by assuming a constant interface shear strength. The usefulness of these elements in micromechanical modeling of fiber-reinforced composite materials is highlighted.

Ballarini, Roberto↗
NASA NTRS1990

Finite element model reduction application to parametric studies and optimization of rotorcraft structures

As a result of this work, a reduction procedure has been developed which can be applied to large finite element model of airframe type structures. This procedure, which is tailored to be used with MSC/NASTRAN finite element code, is applied to the full airframe dynamic finite element model of AH-64A Attack Helicopter. The applicability of the resulting reduced model to parametric and optimization studies is examined. Through application of the design sensitivity analysis, the viability and efficiency of this reduction technique has been demonstrated in a vibration reduction study.

Hashemi-Kia, M.↗
DOE OSTI2024

Finite elements for Matérn-type random fields: Uncertainty in computational mechanics and design optimization

This work highlights an approach for incorporating realistic uncertainties into scientific computing workflows based on finite elements, focusing on prevalent applications in computational mechanics and design optimization. We leverage Matérn-type Gaussian random fields (GRFs) generated using the SPDE method to model aleatoric uncertainties, including environmental influences, variating material properties, and geometric ambiguities. Our focus lies on delivering practical GRF realizations that accurately capture imperfections and variations and understanding how they impact the predictions of computational models as well as the shape and topology of optimized designs. Here we describe a numerical algorithm based on solving a generalized SPDE to sample GRFs on arbitrary meshed domains. The algorithm leverages established techniques and integrates seamlessly with the open-source finite element library MFEM and associated scientific computing workflows, like those found in industrial and national laboratory settings. Our solver scales efficiently for large-scale problems and supports various domain types, including surfaces and embedded manifolds. We showcase its versatility through biomechanics and topology optimization applications, emphasizing the potential to influence these domains. The flexibility and efficiency of SPDE-based GRF generation empowers us to run large-scale optimization problems on 2D and 3D domains, including finding optimized designs on embedded surfaces, and to generate design features and topologies beyond the reach of conventional techniques. Moreover, these capabilities allow us to model and quantify geometric uncertainties on reconstructed submanifolds, such as the interpolated surfaces of cerebral aneurysms provided by postprocessing CT scans. In addition to offering benefits in these specific domains, the proposed techniques transcend specific applications and generalize to arbitrary forward and backward problems in uncertainty quantification involving finite elements.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

Probabilistic Finite Element: Variational Theory

The goal of this research is to provide techniques which are cost-effective and enable the engineer to evaluate the effect of uncertainties in complex finite element models. Embedding the probabilistic aspects in a variational formulation is a natural approach. In addition, a variational approach to probabilistic finite elements enables it to be incorporated within standard finite element methodologies. Therefore, once the procedures are developed, they can easily be adapted to existing general purpose programs. Furthermore, the variational basis for these methods enables them to be adapted to a wide variety of structural elements and to provide a consistent basis for incorporating probabilistic features in many aspects of the structural problem. Tasks concluded include the theoretical development of probabilistic variational equations for structural dynamics, the development of efficient numerical algorithms for probabilistic sensitivity displacement and stress analysis, and integration of methodologies into a pilot computer code.

Belytschko, T.↗
NASA NTRS1995

An evaluation of a coupled microstructural approach for the analysis of functionally graded composites via the finite-element method

A comparison is presented between the predictions of the finite-element analysis and a recently developed higher-order theory for functionally graded materials subjected to a thorough-thickness temperature gradient. In contrast to existing micromechanical theories that utilize classical (i.e., uncoupled) homogenization schemes to calculate micro-level and macro-level stress and displacement fields in materials with uniform or nonuniform fiber spacing (i.e., functionally graded materials), the new theory explicitly couples the microstructural details with the macrostructure of the composite. Previous thermo-elastic analysis has demonstrated that such coupling is necessary when: the temperature gradient is large with respect to the dimension of the reinforcement; the characteristic dimension of the reinforcement is large relative to the global dimensions of the composite and the number of reinforcing fibers or inclusions is small. In these circumstances, the standard micromechanical analyses based on the concept of the representative volume element used to determine average composite properties produce questionable results. The comparison between the predictions of the finite-element method and the higher-order theory presented herein establish the theory's accuracy in predicting thermal and stress fields within composites with a finite number of fibers in the thickness direction subjected to a thorough-thickness thermal gradient.

Pindera, Marek-Jerzy↗