Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Verification of a 2-D to 3-D global/local finite element method for symmetric laminates

A two-dimensional to three-dimensional global/local finite element analysis technique, as applied to both cross-ply and general symmetric composite laminates, is presented. In particular, the local areas of interest for the laminate are the straight free edge in proximity to a hole and curved free edge of a hole. Verification of this technique was accomplished by comparison to results from complete three-dimensional finite element analyses. Examination of the interlaminar stress fields and the amount of time for their determination indicated that the global/local method developed yielded accurate results and was more efficient.

Thompson, Danniella M.↗
NASA NTRS1990

Finite element methods of analysis for 3D inviscid compressible flows

The applicants have developed a finite element based approach for the solution of three-dimensional compressible flows. The procedure enables flow solutions to be obtained on tetrahedral discretizations of computational domains of complex form. A further development was the incorporation of a solution adaptive mesh strategy in which the adaptivity is achieved by complete remeshing of the solution domain. During the previous year, the applicants were working with the Advanced Aerodynamics Concepts Branch at NASA Ames Research Center with an implementation of the basic meshing and solution procedure. The objective of the work to be performed over this twelve month period was the transfer of the adaptive mesh technology and also the undertaking of basic research into alternative flow algorithms for the Euler equations on unstructured meshes.

Peraire, Jaime↗
DOE OSTI2017

Interior energy error estimates for the weak Galerkin finite element method

Consider the Poisson equation in a polytopal domain ΩRC d (d=2,3) as the model problem. In this paper, we study interior energy error estimates for the weak Galerkin finite element approximation to elliptic boundary value problems. In particular, we show that the interior error in the energy norm is bounded by three components: the best local approximation error, the error in negative norms, and the trace error on the element boundaries. This implies that the interior convergence rate can be polluted by solution singularities on the domain boundary, even when the solution is smooth in the interior region. Numerical results are reported to support the theoretical findings. Lastly, to the best of our knowledge, this is the first local energy error analysis that applies to general meshes consisting of polytopal elements and hanging nodes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Modelling NDE pulse-echo inspection of misorientated planar rough defects using an elastic finite element method

Pulse-echo ultrasonic NDE examination of large pressure vessel forgings is a design and construction code requirement in the power generation industry. Such inspections aim to size and characterise potential defects that may have formed during the forging process. Typically these defects have a range of orientations and surface roughnesses which can greatly affect ultrasonic wave scattering behaviour. Ultrasonic modelling techniques can provide insight into defect response and therefore aid in characterisation. However, analytical approaches to solving these scattering problems can become inaccurate, especially when applied to increasingly complex defect geometries. To overcome these limitations a elastic Finite Element (FE) method has been developed to simulate pulse-echo inspections of embedded planar defects. The FE model comprises a significantly reduced spatial domain allowing for a Monte-Carlo based approach to consider multiple realisations of defect orientation and surface roughness. The results confirm that defects aligned perpendicular to the path of beam propagation attenuate ultrasonic signals according to the level of surface roughness. However, for defects orientated away from this plane, surface roughness can increase the magnitude of the scattered component propagating back along the path of the incident beam. This study therefore highlights instances where defect roughness increases the magnitude of ultrasonic scattered signals, as opposed to attenuation which is more often assumed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1978

Application of the finite element method to rotary-wing aeroelasticity

Recent research in rotary-wing aeroelasticity has indicated that all fundamental problems in this area are inherently nonlinear. The non-linearities in this problem are due to the inclusion of finite slopes, due to moderate deflections, in the structural, inertia and aerodynamic operators associated with this aeroelastic problem. In this paper the equations of motion, which are both time and space dependent, for the aeroelastic problem are first formulated in P.D.E. form. Next the equations are linearized about a suitable equilibrium position. The spatial dependence in these equations is discretized using a local Galerkin method of weighted residuals resulting in a finite element formulation of the aeroelastic problem. As an illustration the method is applied to the coupled flap-lag problem of a helicopter rotor blade in hover. Comparison of the solutions with previously published solutions establishes the convergence properties of the method. It is concluded that this formulation is a practical tool for solving rotary-wing aeroelastic stability or response problems.

Friedmann, P.↗
NASA NTRS1978

Study on prestressed concrete reactor vessel structures. II-5: Crack analysis by three dimensional finite elements method of 1/20 multicavity type PCRV subjected to internal pressure

A three-dimensional finite elements analysis is reported of the nonlinear behavior of PCRV subjected to internal pressure by comparing calculated results with test results. As the first stage, an analysis considering the nonlinearity of cracking in concrete was attempted. As a result, it is found possible to make an analysis up to three times the design pressure (50 kg/sqcm), and calculated results agree well with test results.

Source record↗
NASA NTRS1987

Finite element methods of analysis for high speed viscous flows

Over the past three years a finite element based procedure for the solution of high speed viscous compressible flows was developed. The approach followed was to compute steady state solutions via a false transient, using an explicit time stepping scheme, and to attempt to improve the solution quality by incorporating adaptive mesh procedures. The main thrust of the work was to continue on the extension of the approach to the solution of some realistic compressible viscous flows. When flows at high Reynolds number are investigated, is soon becomes apparent that explicit techniwues have to be supplemented if they are to deal effectively with the large variations in element size and aspect ratio which characterize the computational grids necessary for adequate resolutin of the primary flow features. For this reason, the Taylor-Galerkin solution algorithm was rewritten in an explicit/implicit form. Solutions were computed for the problems of a flow past a flat plate, M-3, Re-100; shock/boundary layer interaction, M-2, Re-296000; flow over a compression corner, M-11.68, Re246000; and unifrom flow past a circular cylinder, M-6.34, Re-39770. A summary of the results is included and demonstrates the numerical performance of the scheme.

Source record↗
NASA NTRS1989

Application of the p-version of the finite-element method to global-local problems

A brief survey is given of some recent developments in finite-element analysis technology which bear upon the three main research areas under consideration in this workshop: (1) analysis methods; (2) software testing and quality assurance; and (3) parallel processing. The variational principle incorporated in a finite-element computer program, together with a particular set of input data, determines the exact solution corresponding to that input data. Most finite-element analysis computer programs are based on the principle of virtual work. In the following, researchers consider only programs based on the principle of virtual work and denote the exact displacement vector field corresponding to some specific set of input data by vector u(EX). The exact solution vector u(EX) is independent of the design of the mesh or the choice of elements. Except for very simple problems, or specially constructed test problems, vector u(EX) is not known. Researchers perform a finite-element analysis (or any other numerical analysis) because they wish to make conclusions concerning the response of a physical system to certain imposed conditions, as if vector u(EX) were known.

Szabo, Barna A.↗
NASA NTRS1993

A p-version finite element method for steady incompressible fluid flow and convective heat transfer

A new p-version finite element formulation for steady, incompressible fluid flow and convective heat transfer problems is presented. The steady-state residual equations are obtained by considering a limiting case of the least-squares formulation for the transient problem. The method circumvents the Babuska-Brezzi condition, permitting the use of equal-order interpolation for velocity and pressure, without requiring the use of arbitrary parameters. Numerical results are presented to demonstrate the accuracy and generality of the method.

Winterscheidt, Daniel L.↗
DOE OSTI2014

Enhancing Least-Squares Finite Element Methods Through a Quantity-of-Interest

Here, we introduce an approach that augments least-squares finite element formulations with user-specified quantities-of-interest. The method incorporates the quantity-of-interest into the least-squares functional and inherits the global approximation properties of the standard formulation as well as increased resolution of the quantity-of-interest. We establish theoretical properties such as optimality and enhanced convergence under a set of general assumptions. Central to the approach is that it offers an element-level estimate of the error in the quantity-of-interest. As a result, we introduce an adaptive approach that yields efficient, adaptively refined approximations. Several numerical experiments for a range of situations are presented to support the theory and highlight the effectiveness of our methodology. Notably, the results show that the new approach is effective at improving the accuracy per total computational cost.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

A one-dimensional shock capturing finite element method and multi-dimensional generalizations

Multi-dimensional generalizations of a one-dimensional finite element shock capturing scheme are proposed. A scalar model problem is used to emphasize that 'preferred directions' are important in multi-dimensional applications. Schemes are developed for the two-dimensional Euler equations. One, based upon characteristics, employs the Mach lines and streamlines as preferred directions.

Hughes, T. J. R.↗
NASA NTRS1994

Development of an adaptive hp-version finite element method for computational optimal control

In this research effort, the usefulness of hp-version finite elements and adaptive solution-refinement techniques in generating numerical solutions to optimal control problems has been investigated. Under NAG-939, a general FORTRAN code was developed which approximated solutions to optimal control problems with control constraints and state constraints. Within that methodology, to get high-order accuracy in solutions, the finite element mesh would have to be refined repeatedly through bisection of the entire mesh in a given phase. In the current research effort, the order of the shape functions in each element has been made a variable, giving more flexibility in error reduction and smoothing. Similarly, individual elements can each be subdivided into many pieces, depending on the local error indicator, while other parts of the mesh remain coarsely discretized. The problem remains to reduce and smooth the error while still keeping computational effort reasonable enough to calculate time histories in a short enough time for on-board applications.

Hodges, Dewey H.↗
NASA NTRS1996

Robust Hybrid Finite Element Methods for Antennas and Microwave Circuits

One of the primary goals in this dissertation is concerned with the development of robust hybrid finite element-boundary integral (FE-BI) techniques for modeling and design of conformal antennas of arbitrary shape. Both the finite element and integral equation methods will be first overviewed in this chapter with an emphasis on recently developed hybrid FE-BI methodologies for antennas, microwave and millimeter wave applications. The structure of the dissertation is then outlined. We conclude the chapter with discussions of certain fundamental concepts and methods in electromagnetics, which are important to this study.

Gong, J.↗
DOE OSTI2016

Charged particle tracking through electrostatic wire meshes using the finite element method

Wire meshes are used across many disciplines to accelerate and focus charged particles, however, analytical solutions are non-exact and few codes exist which simulate the exact fields around a mesh with physical sizes. A tracking code based in Matlab-Simulink using field maps generated using finite element software has been developed which tracks electrons or ions through electrostatic wire meshes. The fields around such a geometry are presented as an analytical expression using several basic assumptions, however, it is apparent that computational calculations are required to obtain realistic values of electric potential and fields, particularly when multiple wire meshes are deployed. The tracking code is flexible in that any quantitatively describable particle distribution can be used for both electrons and ions as well as other benefits such as ease of export to other programs for analysis. The code is made freely available and physical examples are highlighted where this code could be beneficial for different applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Higher-order adaptive finite-element methods for Kohn–Sham density functional theory

We present an efficient computational approach to perform real-space electronic structure calculations using an adaptive higher-order finite-element discretization of Kohn–Sham density-functional theory (DFT). To this end, we develop an a priori mesh-adaption technique to construct a close to optimal finite-element discretization of the problem. We further propose an efficient solution strategy for solving the discrete eigenvalue problem by using spectral finite-elements in conjunction with Gauss–Lobatto quadrature, and a Chebyshev acceleration technique for computing the occupied eigenspace. The proposed approach has been observed to provide a staggering 100–200-fold computational advantage over the solution of a generalized eigenvalue problem. Using the proposed solution procedure, we investigate the computational efficiency afforded by higher-order finite-element discretizations of the Kohn–Sham DFT problem. Our studies suggest that staggering computational savings—of the order of 1000-fold—relative to linear finite-elements can be realized, for both all-electron and local pseudopotential calculations, by using higher-order finite-element discretizations. On all the benchmark systems studied, we observe diminishing returns in computational savings beyond the sixth-order for accuracies commensurate with chemical accuracy, suggesting that the hexic spectral-element may be an optimal choice for the finite-element discretization of the Kohn–Sham DFT problem. A comparative study of the computational efficiency of the proposed higher-order finite-element discretizations suggests that the performance of finite-element basis is competing with the plane-wave discretization for non-periodic local pseudopotential calculations, and compares to the Gaussian basis for all-electron calculations to within an order of magnitude. Further, we demonstrate the capability of the proposed approach to compute the electronic structure of a metallic system containing 1688 atoms using modest computational resources, and good scalability of the present implementation up to 192 processors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Smoothness of Functions in Spaces of the Finite Element Method

We find necessary and sufficient conditions for the generalized smoothness of the coordinate functions obtained from the approximate relations. We show that for the coordinate functions the smoothness on their supports is equivalent to that on the boundaries of supports. We obtain conditions for the continuity of Courant type finite element approximations and conditions for the uniqueness of a linear space of such approximations.

97 MATHEMATICS AND COMPUTING↗