SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Mixed finite element models for free vibrations of thin-walled beams

Simple, mixed finite element models are developed for the free vibration analysis of curved thin-walled beams with arbitrary open cross section. The analytical formulation is based on a Vlasov's type thin-walled beam theory with the effects of flexural-torsional coupling, transverse shear deformation and rotary inertia included. The fundamental unknowns consist of seven internal forces and seven generalized displacements of the beam. The element characteristic arrays are obtained by using a perturbed Lagrangian-mixed variational principle. Only C(sup o) continuity is required for the generalized displacements. The internal forces and the Lagrange multiplier are allowed to be discontinuous at interelement boundaries. Numerical results are presented to demonstrate the high accuracy and effectiveness of the elements developed. The standard of comparison is taken to be the solutions obtained by using 2-D plate/shell models for the beams.

Noor, Ahmed K.↗

NASA NTRS1989

Thermoplastic matrix composites - Finite-element analysis of mode I and mode II failure specimens

A finite-element analysis was conducted to evaluate the stress distributions within mode I and mode II failure specimens, assuming both isotropic and orthotropic elastic material properties. The effects of anisotropy on both the magnitude and the location of the highest stress concentration at the vicinity of the crack tip are significant. The results from modeling realistic blunt crack tip geometry and resin-rich zones imply that local variations in the microstructure strongly influence the stress state near the crack tip and therefore the measured fracture properties. In addition, the features of a viscoelastic model for thermoplastic matrices are described. This model will be used in future investigations of matrix-dominated failure phenomena.

Bankert, Ray J.↗

NASA NTRS1989

Finite element computations and experimental studies of flow past an array of plates

Finite element computations and laboratory experiments were carried out to study the mutual wake interaction behind two plates placed side by side. The computations are based on the vorticity-stream function formation of the two-dimensional incompressible Navier-Stokes equations. The flow visualization studies indicate independent vortex sheddings at large gap sizes and metastable biased flow at smaller gap sizes. The computational results exhibit a similar behavior.

Higuchi, H.↗

NASA NTRS1991

An intrinsically n-dimensional generalized flux vector splitting implicit finite element Euler algorithm

A generalized flux-vector splitting implicit Galerkin finite-element algorithm for the Euler equations in curvilinear coordinates for ideal and reacting gases is derived. For an arbitrary equation of state, the curvilinear-coordinate flux vector is split in kinematic and kinetic components, and the associated jacobian matrix eigenvalues explicitly depend on the metric data. After directional semidiscretization, the terminal ordinary differential-equation system is solved via a nonlinearly stable implicit Runge-Kutta scheme in concert with an accurate tensor matrix product factorization. The results for selected two-dimensional supersonic and axisymmetric hypersonic flows validate the algorithm and verify its robustness for curvilinear-coordinate computations. The evolution towards a steady state is achieved for large Courant numbers without indication of numerical instabilities.

Iannelli, G. S.↗

NASA NTRS1992

An efficient finite element method for aircraft de-icing problems

In this paper, a finite element formulation based on an assumed states method is proposed for the solution of heat conduction problems with phase change at a fixed temperature. Attention is directed toward reduction of computer cost through the use of an efficient formulation, solver and algorithm. The procedure is applied to the analysis of an electrothermally deiced aircraft surface.

Huang, J. R.↗

NASA NTRS1992

A least-squares finite element method for incompressible Navier-Stokes problems

A least-squares finite element method, based on the velocity-pressure-vorticity formulation, is developed for solving steady incompressible Navier-Stokes problems. This method leads to a minimization problem rather than to a saddle-point problem by the classic mixed method and can thus accommodate equal-order interpolations. This method has no parameter to tune. The associated algebraic system is symmetric, and positive definite. Numerical results for the cavity flow at Reynolds number up to 10,000 and the backward-facing step flow at Reynolds number up to 900 are presented.

Jiang, Bo-Nan↗

NASA NTRS1992

Finite element approach for transient analysis of multibody systems

A three-dimensional, finite element based formulation for the transient dynamics of constrained multibody systems with trusslike configurations is presented. A convected coordinate system is used to define the rigid-body motion of individual elements in the system. Deformation of each element is defined relative to its convected coordinate system. The formulation is oriented toward joint-dominated structures. Through a series of sequential transformations, the joint degree of freedom is built into the equations of motion of the element to reduce geometric constraints. Based on the derivation, a general-purpose code has been developed. Two examples are presented to illustrate the application of the code.

Wu, Shih-Chin↗

NASA NTRS1993

Finite element analysis of time-independent superconductivity

The development of electromagnetic (EM) finite elements based upon a generalized four-potential variational principle is presented. The use of the four-potential variational principle allows for downstream coupling of EM fields with the thermal, mechanical, and quantum effects exhibited by superconducting materials. The use of variational methods to model an EM system allows for a greater range of applications than just the superconducting problem. The four-potential variational principle can be used to solve a broader range of EM problems than any of the currently available formulations. It also reduces the number of independent variables from six to four while easily dealing with conductor/insulator interfaces. This methodology was applied to a range of EM field problems. Results from all these problems predict EM quantities exceptionally well and are consistent with the expected physical behavior.

Schuler, James J.↗

NASA NTRS1992

Design of an essentially non-oscillatory reconstruction procedure in finite-element type meshes

An essentially non oscillatory reconstruction for functions defined on finite element type meshes is designed. Two related problems are studied: the interpolation of possibly unsmooth multivariate functions on arbitary meshes and the reconstruction of a function from its averages in the control volumes surrounding the nodes of the mesh. Concerning the first problem, the behavior of the highest coefficients of two polynomial interpolations of a function that may admit discontinuities of locally regular curves is studied: the Lagrange interpolation and an approximation such that the mean of the polynomial on any control volume is equal to that of the function to be approximated. This enables the best stencil for the approximation to be chosen. The choice of the smallest possible number of stencils is addressed. Concerning the reconstruction problem, two methods were studied: one based on an adaptation of the so called reconstruction via deconvolution method to irregular meshes and one that lies on the approximation on the mean as defined above. The first method is conservative up to a quadrature formula and the second one is exactly conservative. The two methods have the expected order of accuracy, but the second one is much less expensive than the first one. Some numerical examples are given which demonstrate the efficiency of the reconstruction.

Abgrall, Remi↗

NASA NTRS1991

A finite element large deflection random response of a pipe containing fluid flow

A finite element approach is developed for beam type pipes undergoing large deflections subjected to random loadings. The influence of fluid velocity on the random response is investigated. The root-mean-square (rms) deflections and frequencies for different sound spectrum level values are determined for pipes with both ends either simply supported or clamped. The required number of modes to achieve accurate rms deflections is studied. The prediction of fatigue life is then based on the maximum rms stress. This analytical investigation will help to broaden the basic understanding of the role of fluid flow within structures subjected to random excitations.

Chiang, C. K.↗

NASA NTRS2001

Dynamic Finite Element Predictions for Mars Sample Return Cellular Impact Test #4

The nonlinear finite element program MSC.Dytran was used to predict the impact pulse for (he drop test of an energy absorbing cellular structure. This pre-test simulation was performed to aid in the design of an energy absorbing concept for a highly reliable passive Earth Entry Vehicle (EEV) that will directly impact the Earth without a parachute. In addition, a goal of the simulation was to bound the acceleration pulse produced and delivered to the simulated space cargo container. EEV's are designed to return materials from asteroids, comets, or planets for laboratory analysis on Earth. The EEV concept uses an energy absorbing cellular structure designed to contain and limit the acceleration of space exploration samples during Earth impact. The spherical shaped cellular structure is composed of solid hexagonal and pentagonal foam-filled cells with hybrid graphite-epoxy/Kevlar cell walls. Space samples fit inside a smaller sphere at the enter of the EEV's cellular structure. The material models and failure criteria were varied to determine their effect on the resulting acceleration pulse. Pre-test analytical predictions using MSC.Dytran were compared with the test results obtained from impact test #4 using bungee accelerator located at the NASA Langley Research Center Impact Dynamics Research Facility. The material model used to represent the foam and the proper failure criteria for the cell walls were critical in predicting the impact loads of the cellular structure. It was determined that a FOAMI model for the foam and a 20% failure strain criteria for the cell walls gave an accurate prediction of the acceleration pulse for drop test #4.

Fasanella, Edwin L.↗

NASA NTRS2007

Better Finite-Element Analysis of Composite Shell Structures

A computer program implements a finite-element-based method of predicting the deformations of thin aerospace structures made of isotropic materials or anisotropic fiber-reinforced composite materials. The technique and corresponding software are applicable to thin shell structures in general and are particularly useful for analysis of thin beamlike members having open cross-sections (e.g. I-beams and C-channels) in which significant warping can occur.

Clarke, Gregory↗

DOE OSTI2024

Gyrokinetic electromagnetic particle simulations in triangular meshes with C1 finite elements

The triangular mesh-based gyrokinetic scheme enables comprehensive axis-to-edge studies across the entire plasma volume. Our approach employs triangular finite elements with first-derivative continuity (C1), building on previous work to facilitate gyrokinetic simulations. Additionally, we have adopted the mixed variable/pullback scheme for gyrokinetic electromagnetic particle simulations. The filter-free treatment in the poloidal cross-section with triangular meshes introduces unique features and challenges compared to previous treatments using structured meshes. Our implementation has been validated through benchmarks using ITPA-toroidicity-induced Alfvén eigenmode parameters, showing its capability in moderate to small electron skin depth regimes. Additional examinations using experimental parameters confirm its applicability to realistic plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

A space-time adaptive finite element method with exponential time integrator for the phase field model of pitting corrosion

In this paper we propose a space-time adaptive finite element method for the phase field model of pitting corrosion, which is a parabolic partial differential equation system consisting of a phase variable and a concentration variable. A major challenge in solving this phase field model is that the problem is very stiff, which makes the time step size extremely small for standard temporal discretizations. Another difficulty is that a high spatial resolution is required to capture the steep gradients within the diffused interface, which results in very large number of degrees of freedom for uniform meshes. To overcome the stiffness of this model, we combine the Rosenbrock–Euler exponential integrator with Crank–Nicolson scheme for the temporal discretization. Moreover, by exploiting the fact that the speed of the corroding interface decreases with time, we derive an adaptive time stepping formula. For the spatial approximation, we propose a simple and efficient strategy to generate adaptive meshes that reduces the computational cost significantly. Thus, the proposed method utilizes local adaptivity and mesh refinement for efficient simulation of the corrosive dissolution over long times in heterogeneous media with complex microstructures. We also present an extensive set of numerical experiments in both two and three dimensional spaces to demonstrate efficiency and robustness of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2000

A Posteriori Error Estimation for Finite Volume and Finite Element Approximations Using Broken Space Approximation

We consider a posteriori error estimates for finite volume and finite element methods on arbitrary meshes subject to prescribed error functionals. Error estimates of this type are useful in a number of computational settings: (1) quantitative prediction of the numerical solution error, (2) adaptive meshing, and (3) load balancing of work on parallel computing architectures. Our analysis recasts the class of Godunov finite volumes schemes as a particular form of discontinuous Galerkin method utilizing broken space approximation obtained via reconstruction of cell-averaged data. In this general framework, weighted residual error bounds are readily obtained using duality arguments and Galerkin orthogonality. Additional consideration is given to issues such as nonlinearity, efficiency, and the relationship to other existing methods. Numerical examples are given throughout the talk to demonstrate the sharpness of the estimates and efficiency of the techniques. Additional information is contained in the original.

Barth, Timothy J.↗

DOE OSTI2025

A method for bounding high-order finite element functions: Applications to mesh validity and bounds-preserving limiters

We introduce a novel method for bounding high-order multi-dimensional polynomials in finite element approximations. The method involves precomputing optimal piecewise-linear bounding boxes for polynomial basis functions, which can then be used to locally bound any combination of these basis functions. This approach can be applied to any element/basis type at any approximation order, can provide local (i.e., subcell) extremum bounds to a desired level of accuracy, and can be evaluated efficiently on-the-fly in simulations. Furthermore, we show that this approach generally yields more accurate bounds in comparison to traditional methods based on convex hull properties (e.g., Bernstein polynomials). Furthermore, the efficacy of this technique is shown in applications such as mesh validity checks and optimization for high-order curved meshes, where positivity of the element Jacobian determinant can be ensured throughout the entire element, and continuously bounds-preserving limiters for hyperbolic systems, which can enforce maximum principle bounds across the entire solution polynomial.

Bounding box↗

NASA NTRS1984

Augmented weak forms and element-by-element preconditioners: Efficient iterative strategies for structural finite elements. A preliminary study

A weak formulation in structural analysis that provides well conditioned matrices suitable for iterative solutions is presented. A mixed formulation ensures the proper representation of the problem and the constitutive relations are added in a penalized form. The problem is solved by a double conjugate gradient algorithm combined with an element by element approximate factorization procedure. The double conjugate gradient strategy resembles Uzawa's variable-length type algorithms the main difference is the presence of quadratic terms in the mixed variables. In the case of shear deformable beams these terms ensure that the proper finite thickness solution is obtained.

Muller, A.↗

NASA NTRS1995

Finite Element Modeling on Scalable Parallel Computers

A coupled finite element-integral equation was developed to model fields scattered from inhomogenous, three-dimensional objects of arbitrary shape. This paper outlines how to implement the software on a scalable parallel processor.

modeling three-dimensional coupled finite element ↗