Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite elements”

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

Solving sparse finite element problems on neuromorphic hardware

The finite element method (FEM) is one of the most important and ubiquitous numerical methods for solving partial differential equations (PDEs) on computers for scientific and engineering discovery. Applying the FEM to larger and more detailed scientific models has driven advances in high-performance computing for decades. Here we demonstrate that scalable spiking neuromorphic hardware can directly implement the FEM by constructing a spiking neural network that solves the large, sparse, linear systems of equations at the core of the FEM. We show that for the Poisson equation, a fundamental PDE in science and engineering, our neural circuit achieves meaningful levels of numerical accuracy and close to ideal scaling on modern, inherently parallel and energy-efficient neuromorphic hardware, specifically Intel’s Loihi 2 neuromorphic platform. We illustrate extensions to irregular mesh geometries in both two and three dimensions as well as other PDEs such as linear elasticity. Our spiking neural network is constructed from a recurrent network model of the brain’s motor cortex and, in contrast to black-box deep artificial neural network-based methods for PDEs, directly translates the well-understood and trusted mathematics of the FEM to a natively spiking neuromorphic algorithm.

Applied mathematics↗

Finite element modeling of a cantilever beam

Finite element modeling of uniform cantilever beams and computation of transverse frequencies due to combinations of bending with rotatory inertia and shear flexibility

CANTILEVER BEAM↗

New triangular and quadrilateral plate-bending finite elements

A nonconforming plate-bending finite element of triangular shape and associated quadrilateral elements are developed. The transverse displacement is approximated within the element by a quintic polynomial. The formulation takes into account the effects of transverse shear deformation. Results of the static and dynamic analysis of a square plate, with edges simply supported or clamped, are compared with exact solutions. Good accuracy is obtained in all calculations.

Narayanaswami, R.↗

The control data "GIRAFFE" system for interactive graphic finite element analysis

The Graphical Interface for Finite Elements (GIRAFFE) general purpose interactive graphics application package was described. This system may be used as a pre/post processor for structural analysis computer programs. It facilitates the operations of creating, editing, or reviewing all the structural input/output data on a graphics terminal in a time-sharing mode of operation. An application program for a simple three-dimensional plate problem was illustrated.

Park, S.↗

Investigation of rotational transonic flows through ducts using a finite element scheme

An application of the finite element method is presented in order to study the flow through a two-dimensional channel including the choked and nearly-choked flow conditions. The mathematical procedure provides a combined treatment of potential and Euler equations, where the steady Euler equations are solved through the integration of a pseudo-time system which is equivalent to a relaxation scheme; the isentropic, potential flow equations are embedded in this formulation. To analyze the transonic flow through a channel with near-choked flow conditions, the flow through a parallel channel with a 10 percent circular arc bump is calculated as a sample problem.

Akay, H. U.↗

Finite element analysis of wrinkling membranes

The finite element analysis of wrinkling membranes was investigated. The determination of stresses and deformations within large partly wrinkled membrane surfaces is a problem of significant technical interest in such areas as conceptual design and analysis of ultra lightweight spacecraft structures. A closed-form solution to an axisymmetric problem involving partial wrinkling of an inflated shallow membrane was obtained. In particular, a membrane in the shape of a sperical annulus was considered. The outer edge of the annulus was assumed to be fixed so that no displacements occur along the outer perimeter. The inner edge is assumed to be clamped to a rigid movable plug. Solutions for the complete stress, strain, and displacement fields under the assumption of inextensional material behavior are presented for the case of pure torsional loads applied to the plug, and for the case of pure axial loads applied to the plug.

Miller, R. K.↗

Generalizing the TRAPRG and TRAPAX finite elements

The NASTRAN TRAPRG and TRAPAX finite elements are very restrictive as to shape and grid point numbering. The elements must be trapezoidal with two sides parallel to the radial axis. In addition, the ordering of the grid points on the element connection card must follow strict rules. The paper describes the generalization of these elements so that these restrictions no longer apply.

Hurwitz, M. M.↗

Recent advances in error estimation and adaptive improvement of finite element calculations

Recent advances in adaptive finite elements are summarized. General concepts behind a posteriori error estimation, h-method data management, and algorithms for fluid-mechanics applications are then examined, and some results of numerical experiments with new adaptive codes are given. Numerical examples include supersonic flow over a 20-deg ramp, supersonic flow in expansion corners, the reflecting-shock problem, and the rotating-cone problem. Finally, future directions of research in the field are outlined.

Oden, J. T.↗

Finite element methods in fracture mechanics

Finite-element methodology specific to the analysis of fracture mechanics problems is reviewed. Primary emphasis is on the important algorithmic developments which have enhanced the numerical modeling of fracture processes. Methodologies to address elastostatic problems in two and three dimensions, elastodynamic problems, elastoplastic problems, special considerations for three-dimensional nonlinear problems, and the modeling of stable crack growth are reviewed. In addition, the future needs of the fracture community are discussed and open questions are identified.

Liebowitz, H.↗

Vector algorithms for geometrically nonlinear 3D finite element analysis

Algorithms for geometrically nonlinear finite element analysis are presented which exploit the vector processing capability of the VPS-32, which is closely related to the CYBER 205. By manipulating vectors (which are long lists of numbers) rather than individual numbers, very high processing speeds are obtained. Long vector lengths are obtained without extensive replication or reordering by storage of intermediate results in strategic patterns at all stages of the computations. Comparisons of execution times with those from programs using either scalar or other vector programming techniques indicate that the algorithms presented are quite efficient.

Whitcomb, John D.↗

A globally well-posed finite element algorithm for aerodynamics applications

A finite element CFD algorithm is developed for Euler and Navier-Stokes aerodynamic applications. For the linear basis, the resultant approximation is at least second-order-accurate in time and space for synergistic use of three procedures: (1) a Taylor weak statement, which provides for derivation of companion conservation law systems with embedded dispersion-error control mechanisms; (2) a stiffly stable second-order-accurate implicit Rosenbrock-Runge-Kutta temporal algorithm; and (3) a matrix tensor product factorization that permits efficient numerical linear algebra handling of the terminal large-matrix statement. Thorough analyses are presented regarding well-posed boundary conditions for inviscid and viscous flow specifications. Numerical solutions are generated and compared for critical evaluation of quasi-one- and two-dimensional Euler and Navier-Stokes benchmark test problems.

Iannelli, G. S.↗

Fatigue crack growth reliability by probabilistic finite elements

Fusion of the probabilistic finite-element method and reliability analysis for probabilistic fatigue-crack growth is presented. A comprehensive method for determining the probability of fatigue failure for mixed-mode cyclic loading is also presented. The loading is mixed-mode with randomness in the initial and final crack lengths, initial crack angle and position, material properties, crack-growth law, crack-direction law, and loading. The methodology consists of calculating the reliability index via an optimization procedure which is used to calculate the probability of fatigue failure. Performance of the methodology presented is demonstrated on a classical mode-I fatigue problem.

Besterfield, Glen H.↗

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↗

Analysis of Fluid Gauge Sensor for Zero or Microgravity Conditions using Finite Element Method

In this paper the Finite Element Method (FEM) is presented for mass/volume gauging of a fluid in a tank subjected to zero or microgravity conditions. In this approach first mutual capacitances between electrodes embedded inside the tank are measured. Assuming the medium properties the mutual capacitances are also estimated using FEM approach. Using proper non-linear optimization the assumed properties are updated by minimizing the mean square error between estimated and measured capacitances values. Numerical results are presented to validate the present approach.

Deshpande, Manohar D.↗

An Approach for Assessing Delamination Propagation Capabilities in Commercial Finite Element Codes

An approach for assessing the delamination propagation capabilities in commercial finite element codes is presented and demonstrated for one code. For this investigation, the Double Cantilever Beam (DCB) specimen and the Single Leg Bending (SLB) specimen were chosen for full three-dimensional finite element simulations. First, benchmark results were created for both specimens. Second, starting from an initially straight front, the delamination was allowed to propagate. Good agreement between the load-displacement relationship obtained from the propagation analysis results and the benchmark results could be achieved by selecting the appropriate input parameters. Selecting the appropriate input parameters, however, was not straightforward and often required an iterative procedure. Qualitatively, the delamination front computed for the DCB specimen did not take the shape of a curved front as expected. However, the analysis of the SLB specimen yielded a curved front as may be expected from the distribution of the energy release rate and the failure index across the width of the specimen. Overall, the results are encouraging but further assessment on a structural level is required.

Krueger, Ronald↗

An Approach for Assessing Delamination Propagation Capabilities in Commercial Finite Element Codes

An approach to assessing the delamination propagation capabilities in commercial finite element codes is presented and demonstrated for one code. For this investigation, the Double Cantilever Beam (DCB) specimen and the Single Leg Bending (SLB) specimen were chosen for full three-dimensional finite element simulations. First, benchmark results were created for both specimens. Second, starting from an initially straight front, the delamination was allowed to propagate. Good agreement between the load-displacement relationship obtained from the propagation analysis results and the benchmark results could be achieved by selecting the appropriate input parameters. Selecting the appropriate input parameters, however, was not straightforward and often required an iterative procedure. Qualitatively, the delamination front computed for the DCB specimen did not take the shape of a curved front as expected. However, the analysis of the SLB specimen yielded a curved front as may be expected from the distribution of the energy release rate and the failure index across the width of the specimen. Overall, the results are encouraging but further assessment on a structural level is required.

Krueger, Ronald↗