Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element problems”

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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Stabilized bases for high-order, interpolation semi-Lagrangian, element-based tracer transport

In a computational fluid model of the atmosphere, the advective transport of trace species, or tracers, can be computationally expensive. For efficiency, models often use semi-Lagrangian advection methods. High-order interpolation semi-Lagrangian (ISL) methods, in particular, can be extremely efficient, if the problem of property preservation specific to them can be addressed. Atmosphere models often use geometrically and logically nonuniform grids for efficiency and, as a result, element-based discretizations. Such grids and discretizations make stability a particular problem for ISL methods. Generally, high-order, element-based ISL methods that use the natural polynomial interpolant associated with a nodal finite-element discretization are unstable. Here, we derive new bases having order of accuracy up to nine, with positive nodal weights, that stabilize the element-based ISL method. We use these bases to construct the linear advection operator in the property-preserving Interpolation Semi-Lagrangian Element-based Transport (Islet) method. Then we discuss key software implementation details. Finally, we show performance results for the Energy Exascale Earth System Model's atmosphere dynamical core, comparing the original and new transport methods. These simulations used up to 27,600 Graphical Processing Units (GPU) on the Oak Ridge Leadership Computing Facility's Summit supercomputer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

MAPS: the MFEM Anisotropic Plasma Solver

Simulating magnetically confined fusion plasmas presents a uniquely challenging problem due to the nonlinear anisotropic heat conduction. We introduce the MAPS (MFEM Anisotropic Plasma Solver) tool, which uses a high-order finite element method to compute transport solutions on unstructured meshes. We show results for a set of three 2-D verification tests, two of which demonstrate the expected convergence properties for various mesh resolutions and polynomial degrees. We then discuss the convergence rate for the third test.

Barnett, Rhea [ORNL] (ORCID:0000000317527979)↗
DOE OSTI2016

Mimetic Finite Differences and Algebraic Multilevel Methods for Reservoir Performance Predictions: CRADA Final Report

The need to better predict and control reservoir performance is driving the development of increasingly sophisticated reservoir models that use unstructured general polyhedral grids to capture highly detailed stratigraphy. However, both discretization methods and scalable iterative solution algorithms for these situations are an active area of research. In this work we compared the performance of a new Mimetic Finite Difference (MFD) method with two new Mixed Finite Element methods, and LANL and the CRADA Participant investigated the performance of Los Alamos National Laboratory’s algebraic multigrid code on realistic flow problems. Furthermore, uncertainty quantification requires rapid evaluations of reservoir models. To do this, fine-scale reservoir models need to be coarsened and yet important fine-scale features and processes that have a significant impact on flows need to be preserved. To this end, we investigated the use of multi-resolution methods, such as the Multilevel Multiscale Mimetic (M3) method, to coarsen reservoir models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Implementation and Verification of the Sesame Equation of State Database into Abaqus

Towards ability to perform advanced engineering calculations of shock driven processes requiring an accurate prediction of temperature inside the shock, the Sesame equation of state database has been implemented into Abaqus explicit, a finite element analysis code commonly used in the analysis community. The interface algorithm and supporting details are discussed, and results from a simple benchmark problem comparing Abaqus predictions to other previously validated calculations and published data are presented. Features of the Sesame-Abaqus interface are also discussed and a sample input file is provided for the interested user.

42 ENGINEERING↗
NASA NTRS1975

Numerical solution of quadratic matrix equations for free vibration analysis of structures

This paper is concerned with the efficient and accurate solution of the eigenvalue problem represented by quadratic matrix equations. Such matrix forms are obtained in connection with the free vibration analysis of structures, discretized by finite 'dynamic' elements, resulting in frequency-dependent stiffness and inertia matrices. The paper presents a new numerical solution procedure of the quadratic matrix equations, based on a combined Sturm sequence and inverse iteration technique enabling economical and accurate determination of a few required eigenvalues and associated vectors. An alternative procedure based on a simultaneous iteration procedure is also described when only the first few modes are the usual requirement. The employment of finite dynamic elements in conjunction with the presently developed eigenvalue routines results in a most significant economy in the dynamic analysis of structures.

Gupta, K. K.↗
NASA NTRS1985

Fundamental studies of structure borne noise for advanced turboprop applications

The transmission of sound generated by wing-mounted, advanced turboprop engines into the cabin interior via structural paths is considered. The structural model employed is a beam representation of the wing box carried into the fuselage via a representative frame type of carry through structure. The structure for the cabin cavity is a stiffened shell of rectangular or cylindrical geometry. The structure is modelled using a finite element formulation and the acoustic cavity is modelled using an analytical representation appropriate for the geometry. The structural and acoustic models are coupled by the use of hard wall cavity modes for the interior and vacuum structural modes for the shell. The coupling is accomplished using a combination of analytical and finite element models. The advantage is the substantial reduction in dimensionality achieved by modelling the interior analytically. The mathematical model for the interior noise problem is demonstrated with a simple plate/cavity system which has all of the features of the fuselage interior noise problem.

Eversman, W.↗
NASA NTRS1987

Concepts for radically increasing the numerical convergence rate of the Euler equations

Integral equation and finite difference methods have been developed for solving transonic flow problems using linearized forms of the transonic small disturbance and Euler equations. A key element is the use of a strained coordinate system in which the shock remains fixed. Additional criteria are developed to determine the free parameters in the coordinate straining; these free parameters are functions of the shock location. An integral equation analysis showed that the shock is located by ensuring that no expansion shocks exist in the solution. The expansion shock appears as oscillations in the solution near the sonic line, and the correct shock location is determined by removing these oscillations. A second objective was to study the ability of the Euler equation to model separated flow.

Nixon, David↗
NASA NTRS1988

General Rotorcraft Aeromechanical Stability Program (GRASP) version 1.03: User's manual

The Rotorcraft Dynamics Division, Aeroflightdynamics Directorate, U.S. Army Aviation Research and Technology Activity has developed the General Rotorcraft Aeromechanical Stability Program (GRASP) to perform calculations that will assess the stability of rotorcraft in hovering flight and ground contact conditions. The program is designed to be state-of-the-art, hybrid, finite-element/multibody code that can be applied to all existing and future helicopter configurations. While GRASP was specifically designed to solve rotorcraft stability problems, its innovative structure and formulation allow for application to a wide range of structures. This manual describes the preparation of the input file required by Version 1.03 of GRASP, the procedures used to invoke GRASP on the NASA Ames Research Center CRAY X-MP 48 computer, and the interpretation of the output produced by GRASP. The parameters used by the input file are defined, and summaries of the input file and the job control language are included.

Hopkins, A. Stewart↗
NASA NTRS1987

Re-analysis procedure based on the mixed formulation

An analysis procedure is presented for large-scale structural systems, with large numbers of degrees of freedom and design variables. The procedure uses a mixed formulation with the fundamental unknowns consisting of both stress and displacement parameters. Other elements of the procedure include: (1) lumping the design variables into a single tracing parameter; (2) operator splitting or additive decomposition of different arrays in the finite element equations into the corresponding arrays of the original structure plus correction terms; and (3) application of a reduction method through the use of the finite element method and the classical Bubnov-Galerkin technique. The re-analysis procedure is applied to the linear static and free vibration problems of plate and shell structures.

Noor, Ahmed K.↗
NASA NTRS1989

Probalistic structural analysis methods for select space propulsion system structural components

A summary of the status of this five-year project which is now in its third year of research and development is presented. The goal of the project is the development of several methodologies for probabilistic structural modeling. Probabilistic structural modeling consists of stochastic models of material properties, part geometries, boundary conditions, as well as loading conditions. The current presentation focuses on one methodology - coupling of an advanced finite element structural analysis code with probabilistic modeling strategies. The essential algorithm developments for combining the finite element and probabilistic analysis methods are reported. The validity of the resulting probabilistic structural analysis method is confirmed through a series of test problems with exact results based on Monte Carlo simulations. Additionally, the applicability of the method to a Space Propulsion System (a turbine blade) is demonstrated for static stresses.

Cruse, T. A.↗
NASA NTRS1990

ZIP3D: An elastic and elastic-plastic finite-element analysis program for cracked bodies

ZIP3D is an elastic and an elastic-plastic finite element program to analyze cracks in three dimensional solids. The program may also be used to analyze uncracked bodies or multi-body problems involving contacting surfaces. For crack problems, the program has several unique features including the calculation of mixed-mode strain energy release rates using the three dimensional virtual crack closure technique, the calculation of the J integral using the equivalent domain integral method, the capability to extend the crack front under monotonic or cyclic loading, and the capability to close or open the crack surfaces during cyclic loading. The theories behind the various aspects of the program are explained briefly. Line-by-line data preparation is presented. Input data and results for an elastic analysis of a surface crack in a plate and for an elastic-plastic analysis of a single-edge-crack-tension specimen are also presented.

Shivakumar, K. N.↗
NASA NTRS1992

Automating FEA programming

In this paper we briefly describe a combined symbolic and numeric approach for solving mathematical models on parallel computers. An experimental software system, PIER, is being developed in Common Lisp to synthesize computationally intensive and domain formulation dependent phases of finite element analysis (FEA) solution methods. Quantities for domain formulation like shape functions, element stiffness matrices, etc., are automatically derived using symbolic mathematical computations. The problem specific information and derived formulae are then used to generate (parallel) numerical code for FEA solution steps. A constructive approach to specify a numerical program design is taken. The code generator compiles application oriented input specifications into (parallel) FORTRAN77 routines with the help of built-in knowledge of the particular problem, numerical solution methods and the target computer.

Sharma, Naveen↗
NASA NTRS1993

Large-scale computation of incompressible viscous flow by least-squares finite element method

The least-squares finite element method (LSFEM) based on the velocity-pressure-vorticity formulation is applied to large-scale/three-dimensional steady incompressible Navier-Stokes problems. This method can accommodate equal-order interpolations and results in symmetric, positive definite algebraic system which can be solved effectively by simple iterative methods. The first-order velocity-Bernoulli function-vorticity formulation for incompressible viscous flows is also tested. For three-dimensional cases, an additional compatibility equation, i.e., the divergence of the vorticity vector should be zero, is included to make the first-order system elliptic. The simple substitution of the Newton's method is employed to linearize the partial differential equations, the LSFEM is used to obtain discretized equations, and the system of algebraic equations is solved using the Jacobi preconditioned conjugate gradient method which avoids formation of either element or global matrices (matrix-free) to achieve high efficiency. To show the validity of this scheme for large-scale computation, we give numerical results for 2D driven cavity problem at Re = 10000 with 408 x 400 bilinear elements. The flow in a 3D cavity is calculated at Re = 100, 400, and 1,000 with 50 x 50 x 50 trilinear elements. The Taylor-Goertler-like vortices are observed for Re = 1,000.

Jiang, Bo-Nan↗
NASA NTRS2002

Genetic Optimization of a Tensegrity Structure

Marshall Space Flight Center (MSFC) is charged with developing advanced technologies for space telescopes. The next generation of space optics will be very large and lightweight. Tensegrity structures are built of compressive members (bars), and tensile members (strings). For most materials, the tensile strength of a longitudinal member is larger than its buckling strength; therefore a large stiffness to mass ratio can be achieved by increasing the use of tensile members. Tensegrities are the epitome of lightweight structures, since they take advantage of the larger tensile strength of materials. The compressive members of tensegrity structures are disjoint allowing compact storage of the structure. The structure has the potential to eliminate the requirement for assembly by man in space; it can be deployed by adjustments in its cable tension. A tensegrity structure can be more reliably modeled since none of the individual members experience bending moments. (Members that experience deformation in more than one dimension are much harder to model.) Structures that can be more precisely modeled can be more precisely controlled. Furthermore, an astoundingly wide variety of natural systems, including carbon atoms, water molecules, proteins, viruses, cells, tissues and even human and other living creatures are tensegrity structures. Through the process of evolution, nature continually improves the design of living creatures for the environment they live in. Since tensegrities are nature's structure of choice, it is conceivable that they have other benefits we are unaware of. A. Keane and S. Brown designed a satellite boom truss system with an enhanced vibration performance. They started with a standard truss system, then used a genetic algorithm to alter the design, while optimizing the vibration performance. An improvement of over 20,000% in frequency-averaged energy levels was obtained using this approach. In this report an introduction to tensegrity structures is given, along with a description of how to generate the nodal coordinates and connectivity of a multiple stage cylindrical tensegrity structure. A description of how finite elements can be used to develop a stiffness and mass matrix so that the modes of vibration can be determined from the eigenvalue problem is shown. A brief description of a micro genetic algorithm is then presented.

Jaime R. Taylor↗
NASA NTRS2002

Nonlinear Pressurization and Modal Analysis Procedure for Dynamic Modeling of Inflatable Structures

An introduction and set of guidelines for finite element dynamic modeling of nonrigidized inflatable structures is provided. A two-step approach is presented, involving 1) nonlinear static pressurization of the structure and updating of the stiffness matrix and 2) hear normal modes analysis using the updated stiffness. Advantages of this approach are that it provides physical realism in modeling of pressure stiffening, and it maintains the analytical convenience of a standard bear eigensolution once the stiffness has been modified. Demonstration of the approach is accomplished through the creation and test verification of an inflated cylinder model using a large commercial finite element code. Good frequency and mode shape comparisons are obtained with test data and previous modeling efforts, verifying the accuracy of the technique. Problems encountered in the application of the approach, as well as their solutions, are discussed in detail.

Smalley, Kurt B.↗
NASA NTRS2006

A Theoretical Investigation of Composite Overwrapped Pressure Vessel (COPV) Mechanics Applied to NASA Full Scale Tests

A theoretical investigation of the factors controlling the stress rupture life of the National Aeronautics and Space Agency's (NASA) composite overwrapped pressure vessels (COPVs) continues. Kevlar(TradeMark) fiber overwrapped tanks are of particular concern due to their long usage and the poorly understood stress rupture process in Kevlar(TradeMark) filaments. Existing long term data show that the rupture process is a function of stress, temperature and time. However, due to the presence of a load sharing liner, the manufacturing induced residual stresses and the complex mechanical response, the state of actual fiber stress in flight hardware and test articles is not clearly known. This paper is a companion to the experimental investigation reported in [1] and develops a theoretical framework necessary to design full-scale pathfinder experiments and accurately interpret the experimentally observed deformation and failure mechanisms leading up to static burst in COPVs. The fundamental mechanical response of COPVs is described using linear elasticity and thin shell theory and discussed in comparison to existing experimental observations. These comparisons reveal discrepancies between physical data and the current analytical results and suggest that the vessel's residual stress state and the spatial stress distribution as a function of pressure may be completely different from predictions based upon existing linear elastic analyses. The 3D elasticity of transversely isotropic spherical shells demonstrates that an overly compliant transverse stiffness relative to membrane stiffness can account for some of this by shifting a thin shell problem well into the realm of thick shell response. The use of calibration procedures are demonstrated as calibrated thin shell model results and finite element results are shown to be in good agreement with the experimental results. The successes reported here have lead to continuing work with full scale testing of larger NASA COPV hardware.

Greene, N.↗
NASA NTRS2009

A Theoretical Investigation of Composite Overwrapped Pressure Vessel (COPV) Mechanics Applied to NASA Full Scale Tests

A theoretical investigation of the factors controlling the stress rupture life of the National Aeronautics and Space Administration's (NASA) composite overwrapped pressure vessels (COPVs) continues. Kevlar (DuPont) fiber overwrapped tanks are of particular concern due to their long usage and the poorly understood stress rupture process in Kevlar filaments. Existing long term data show that the rupture process is a function of stress, temperature and time. However due to the presence of a load sharing liner, the manufacturing induced residual stresses and the complex mechanical response, the state of actual fiber stress in flight hardware and test articles is not clearly known. This paper is a companion to a previously reported experimental investigation and develops a theoretical framework necessary to design full-scale pathfinder experiments and accurately interpret the experimentally observed deformation and failure mechanisms leading up to static burst in COPVs. The fundamental mechanical response of COPVs is described using linear elasticity and thin shell theory and discussed in comparison to existing experimental observations. These comparisons reveal discrepancies between physical data and the current analytical results and suggest that the vessel s residual stress state and the spatial stress distribution as a function of pressure may be completely different from predictions based upon existing linear elastic analyses. The 3D elasticity of transversely isotropic spherical shells demonstrates that an overly compliant transverse stiffness relative to membrane stiffness can account for some of this by shifting a thin shell problem well into the realm of thick shell response. The use of calibration procedures are demonstrated as calibrated thin shell model results and finite element results are shown to be in good agreement with the experimental results. The successes reported here have lead to continuing work with full scale testing of larger NASA COPV hardware.

Thesken, John C.↗