SEARCH · Engineering Papers

Results for “finite element expansion”

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 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Modal element method for potential flow in non-uniform ducts: Combining closed form analysis with CFD

An analytical procedure is presented, called the modal element method, that combines numerical grid based algorithms with eigenfunction expansions developed by separation of variables. A modal element method is presented for solving potential flow in a channel with two-dimensional cylindrical like obstacles. The infinite computational region is divided into three subdomains; the bounded finite element domain, which is characterized by the cylindrical obstacle and the surrounding unbounded uniform channel entrance and exit domains. The velocity potential is represented approximately in the grid based domain by a finite element solution and is represented analytically by an eigenfunction expansion in the uniform semi-infinite entrance and exit domains. The calculated flow fields are in excellent agreement with exact analytical solutions. By eliminating the grid surrounding the obstacle, the modal element method reduces the numerical grid size, employs a more precise far field boundary condition, as well as giving theoretical insight to the interaction of the obstacle with the mean flow. Although the analysis focuses on a specific geometry, the formulation is general and can be applied to a variety of problems as seen by a comparison to companion theories in aeroacoustics and electromagnetics.

Baumeister, Kenneth J.↗

DOE OSTI2022

A coupling approach for linear elasticity problems with spatially non-coincident discretized interfaces

Here we present a new method for coupled linear elasticity problems whose finite element discretization may lead to spatially non-coincident discretized interfaces. Our approach combines the classical Dirichlet–Neumann coupling formulation with a new set of discretized interface conditions obtained through Taylor series expansions. We show that these conditions ensure linear consistency of the coupled finite element solution. We then formulate an iterative solution method for the coupled discrete system and apply the new coupling approach to two representative settings for which we also provide several numerical illustrations. The first setting is a mesh-tying problem in which both coupled structures have the same Lamé parameters whereas the second setting is an interface problem for which the Lamé parameters in the two coupled structures are different.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1990

A finite element conjugate gradient FFT method for scattering

An extension of a two dimensional formulation is presented for a three dimensional body of revolution. With the introduction of a Fourier expansion of the vector electric and magnetic fields, a coupled two dimensional system is generated and solved via the finite element method. An exact boundary condition is employed to terminate the mesh and the fast fourier transformation (FFT) is used to evaluate the boundary integrals for low O(n) memory demand when an iterative solution algorithm is used. By virtue of the finite element method, the algorithm is applicable to structures of arbitrary material composition. Several improvements to the two dimensional algorithm are also described. These include: (1) modifications for terminating the mesh at circular boundaries without distorting the convolutionality of the boundary integrals; (2) the development of nonproprietary mesh generation routines for two dimensional applications; (3) the development of preprocessors for interfacing SDRC IDEAS with the main algorithm; and (4) the development of post-processing algorithms based on the public domain package GRAFIC to generate two and three dimensional gray level and color field maps.

Collins, Jeffery D.↗

NASA NTRS1984

A spectral element method for fluid dynamics - Laminar flow in a channel expansion

A spectral element method that combines the generality of the finite element method with the accuracy of spectral techniques is proposed for the numerical solution of the incompressible Navier-Stokes equations. In the spectral element discretization, the computational domain is broken into a series of elements, and the velocity in each element is represented as a high-order Lagrangian interpolant through Chebyshev collocation points. The hyperbolic piece of the governing equations is then treated with an explicit collocation scheme, while the pressure and viscous contributions are treated implicitly with a projection operator derived from a variational principle. The implementation of the technique is demonstrated on a one-dimensional inflow-outflow advection-diffusion equation, and the method is then applied to laminar two-dimensional (separated) flow in a channel expansion. Comparisons are made with experiment and previous numerical work.

Patera, A. T.↗

NASA NTRS1990

Large scale nonlinear numerical optimal control for finite element models of flexible structures

This paper discusses the development of large scale numerical optimal control algorithms for nonlinear systems and their application to finite element models of structures. This work is based on our expansion of the optimal control algorithm (DDP) in the following steps: improvement of convergence for initial policies in non-convex regions, development of a numerically accurate penalty function method approach for constrained DDP problems, and parallel processing on supercomputers. The expanded constrained DDP algorithm was applied to the control of a four-bay, two dimensional truss with 12 soft members, which generates geometric nonlinearities. Using an explicit finite element model to describe the structural system requires 32 state variables and 10,000 time steps. Our numerical results indicate that for constrained or unconstrained structural problems with nonlinear dynamics, the results obtained by our expanded constrained DDP are significantly better than those obtained using linear-quadratic feedback control.

Shoemaker, Christine A.↗

NASA NTRS1981

Static shape determination and control for a large space antenna

An integral operator approach is used to derive solutions to static shape determination and control problems associated with large space structures. Problem assumptions include a linear self-adjoint system model, observations and control forces at discrete points, and quadratic performance criteria for the comparison of estimates or control forces. Results are illustrated by simulations with a finite element model of a large space antenna. Modal expansions for terms in the solution algorithms are presented, using modes from the static or associated dynamic model. These expansions provide approximate solutions in the event that a closed form analytical solution to the system boundary value problem is not available.

Weeks, C.↗

NASA NTRS1993

An augmented Young-Laplace model of an evaporating meniscus in a micro-channel with high heat flux

High flux evaporations from a steady meniscus formed in a 2 micron channel is modeled using the augmented Young-Laplace equation. The heat flux is found to be a function of the long range van der Waals dispersion force which represents interfacial conditions between heptane and various substrates. Heat fluxes of (1.3-1.6) x 10(exp 6) W/m(exp 2) based on the width of the channel are obtained for heptane completely wetting the substrate at 100 C. Small channels are used to obtain these large fluxes. Even though the real contact angle is 0 deg, the apparent contact angle is found to vary between 24.8 deg and 25.6 deg. The apparent contact angle, which represents viscous losses near the contact line, has a large effect on the heat flow rate because of its effect on capillary suction and the area of the meniscus. The interfacial heat flux is modeled using kinetic theory for the evaporation rate. The superheated state depends on the temperature and the pressure of the liquid phase. The liquid pressure differs from the pressure of the vapor phase due to capillarity and long range van der Waals dispersion forces which are relevant in the ultra think film formed at the leading edge of the meniscus. Important pressure gradients in the thin film cause a substantial apparent contact angle for a complete wetting system. The temperature of the liquid is related to the evaporation rate and to the substrate temperature through the steady heat conduction equation. Conduction in the liquid phase is calculated using finite element analysis except in the vicinity of the thin film. A lubrication theory solution for the thin film is combined with the finite element analysis by the method of matched asymptotic expansions.

Wayner, P. C., Jr.↗

DOE OSTI2023

A structural model of the long-term degradation of the concrete biological shield

The concrete biological shield (CBS) of light water reactors is exposed to high neutron radiation dose in the long term, which may lead to the degradation of the concrete’s mechanical properties. Given the important shielding role of the CBS, it is necessary to investigate the irradiation effects at the structural scale and provide estimates of the damage extent from the wall’s inner surface to study potential license renewals. For this purpose, we developed a mechanical model accounting for radiation-induced expansion, creep, and damage in concrete using the Grizzly finite element code, informed by ex-core neutron flux calculations using the VERA tool. The model was applied to a 3D CBS structure represented by the CBS wall, a steel liner, reinforcement bars, and a concrete base mat and evaluated damage at 40, 60, and 80 years of operation. The VERA model predicted a maximum fluence of approximately 2 x 10 19 ncm -2 at 80 years of operation. The results showed that damage is highest at the inner surface of the CBS wall and gradually decreases with depth. It extends beyond the rebar after 60 years and reaches a depth of approximately 12 cm at 80 years.

42 ENGINEERING↗

NASA NTRS1979

Rail shear test method

The results of both an experimental test program and a finite element analysis of selected graphite/polyimide rail shear test specimens are discussed. The two dimensional finite element analysis includes both mechanical and thermal loading (differential expansion) of the specimens and their elastic rails. Parameters in this analysis of unidirectional and symmetric, balanced angle-ply laminates include ply layup angles, the effect of flexible rails, the method of load introduction to the specimen and the effect of uniform heating of the specimen and rails. Two types of tensile rail shear fixtures were investigated experimentally: a uniform thickness, bolted-rail shear fixture loaded diagonally across the specimen test section; and a tapered thickness, bonded-rail shear fixture loaded axially along the center-line of the specimen test section. Test results include room-temperature and 589K strain data taken from the center of the specimen test section during loading.

Garcia, R.↗

NASA NTRS1983

Characteristics of thermally-induced transverse cracks in graphite epoxy composite laminates

The characteristics of thermally induced transverse cracks in T300/5208 graphite-epoxy cross-ply and quasi-isotropic laminates were investigated both experimentally and analytically. The formation of transverse cracks and the subsequent crack spacing present during cool down to -250 F (116K) and thermal cycling between 250 and -250 F (116 and 394K) was investigated. The state of stress in the vicinity of a transverse crack and the influence of transverse cracking on the laminate coefficient of thermal expansion (CTE) was predicted using a generalized plane strain finite element analysis and a modified shear lag analysis. A majority of the cross-ply laminates experienced transverse cracking during the initial cool down to -250 F whereas the quasi-isotropic laminates remained uncracked. The in situ transverse strength of the 90 degree layers was more than 1.9 times greater than the transverse strength of the unidirectional 90 degree material for all laminates investigated.

Adams, D. S.↗

NASA NTRS1992

Modal element method for scattering of sound by absorbing bodies

The modal element method for acoustic scattering from 2-D body is presented. The body may be acoustically soft (absorbing) or hard (reflecting). The infinite computational region is divided into two subdomains - the bounded finite element domain, which is characterized by complicated geometry and/or variable material properties, and the surrounding unbounded homogeneous domain. The acoustic pressure field is represented approximately in the finite element domain by a finite element solution, and is represented analytically by an eigenfunction expansion in the homogeneous domain. The two solutions are coupled by the continuity of pressure and velocity across the interface between the two subdomains. Also, for hard bodies, a compact modal ring grid system is introduced for which computing requirements are drastically reduced. Analysis for 2-D scattering from solid and coated (acoustically treated) bodies is presented, and several simple numerical examples are discussed. In addition, criteria are presented for determining the number of modes to accurately resolve the scattered pressure field from a solid cylinder as a function of the frequency of the incoming wave and the radius of the cylinder.

Baumeister, Kenneth J.↗

NASA NTRS1993

Modal element method for scattering and absorbing of sound by two-dimensional bodies

The modal element method for acoustic scattering from a 2-D body is presented. The body may be acoustically soft (absorbing) or hard (reflecting). The infinite computational region is divided into two subdomains - the bounded finite element domain, which is characterized by complicated geometry and/or variable material properties, and the surrounding unbounded homogeneous domain. The acoustic pressure field is represented approximately in the finite element domain by a finite element solution, and is represented analytically by an eigenfunction expansion in the homogeneous domain. The two solutions are coupled by the continuity of pressure and velocity across the interface between the two subdomains. Also, for hard bodies, a compact modal ring grid system is introduced for which computing requirements are drastically reduced. Analysis for 2-D scattering from solid and coated (acoustically treated) bodies is presented, and several simple numerical examples are discussed. In addition, criteria are presented for determining the number of modes to accurately resolve the scattered pressure field from a solid cylinder as a function of the frequency of the incoming wave and the radius of the cylinder.

Baumeister, K. J.↗

NASA NTRS2009

Analytical Predictions of Thermal Stress in the Stardust PICA Heatshield Under Reentry Flight Conditions

We performed finite element analyses on a model of the Phenolic Impregnated Carbon Ablator (PICA) heatshield from the Stardust sample return capsule (SRC) to predict the thermal stresses in the PICA material during reentry. The heatshield on the Stardust SRC was a 0.83 m sphere cone, fabricated from a single piece of 5.82 cm-thick PICA. The heatshield performed successfully during Earth reentry of the SRC in January 2006. Material response analyses of the full, axisymmetric PICA heatshield were run using the Two-Dimensional Implicit Ablation, Pyrolysis, and Thermal Response Program (TITAN). Peak surface temperatures were predicted to be 3385K, while the temperature at the PICA backface remained at the estimated initial cold-soak temperature of 278K. Surface recession and temperature distribution results from TITAN, at several points in the reentry trajectory, were mapped onto an axisymmetric finite element model of the heatshield. We used the finite element model to predict the thermal stresses in the PICA from differential thermal expansion. The predicted peak compressive stress in the PICA heatshield was 1.38 MPa. Although this level of stress exceeded the chosen design limit for compressive stresses in PICA tiles for the design of the Orion crew exploration vehicle heatshield, the Stardust heatshield exhibited no obvious mechanical failures from thermal stress. The analyses of the Stardust heatshield were used to assess and adjust the level of conservatism in the finite element analyses in support of the Orion heatshield design.

Squire, Thomas↗

NASA NTRS1981

On the relationship between engineering properties and delamination of composite laminates

The influence of the coefficient of mutal influence, Poisson's ratio and coefficients of thermal and moisture expansion on delamination is studied. Engineering theories are compared to finite element and experimental results. It is shown that the mismatch in coefficients of mutual influence can have a strong influence on delamination with fiber angles in the 10-15 degree range being critical for adjacent layer combinations. The mismatch in coefficient of mutual influence is reduced by a factor of two and the interlaminar shear stress is reduced significantly when the + or - adjacent layers are interspersed between 0 and 90 degree layers. It is shown how the results can be used for design of composite laminates.

Herakovich, C. T.↗

NASA NTRS1993

Nonlinear response of graphite-epoxy wide columns subject to eccentric load

It is shown that including the effects of bend-twist coupling and material nonlinearity on the large displacement flexural response of graphite-epoxy wide columns is necessary to adequately model their response. Results from static tests on 0(30)T and (30/0/-30 sub 5)s laminated specimens under eccentric axial compression are compared to a nonlinear analysis. The analysis is based on a Kantrovich expansion of the displacements in the width coordinate and a finite element representation in the axial coordinate. An incremental, updated Lagrangian formulation and element corotational coordinates are employed. Nonlinear elastic material behavior was measured in a separate set of material characterization tests, and these data are used to compute tangent stillnesses in the updated configuration.

Griffin, O. H., Jr.↗

NASA NTRS1979

Micromechanics of intraply hybrid composites: Elastic and thermal properties

Composite micromechanics are used to derive equations for predicting the elastic and thermal properties of unidirectional intraply hybrid composites. The results predicted using these equations are compared with those predicted using approximate equations based on the rule of mixtures, linear laminate theory, finite element analysis and limited experimental data. The comparisons for three different intraply hybrids indicate that all four methods predict approximately the same elastic properties and are in good agreement with measured data. The micromechanics equations and linear laminate theory predict about the same values for thermal expansion coefficients. The micromechanics equations predict through-the-thickness properties which are in good agreement with the finite element results.

Chamis, C. C.↗

NASA NTRS1993

Advanced reinforcement systems for intermetallic applications

A 2-D axisymmetric model was employed to determine the magnitude of the radial, axial, and hoop stresses caused by the thermal expansion difference between fiber and matrix and which result from the fabrication temperature cycle. Finite element analysis was conducted for single fiber model systems based on SCS-6/Ti3Al+Nb and Al2O3/NiAl. The stress distribution due to the imposition of a graded intermediate layer for each system was determined and included variables of layer thickness and gradation in interlayer chemistry in order to vary the expansion gradient between fiber and matrix. Thermal cycling tests were conducted on sputter coated SCS-6 fibers selectively coated with Ti3Al+Nb, with and without an intermediate layer. Cracking of the Ti3Al+Nb layers was prevented by an interlayer based on Ti-TiN-Ti. The interlayer thickness appeared critical to its efficiency. Similarly, for the case of Al2O3/NiAl, an intermediate layer consisting of a Ni bond coat on the sapphire fiber followed by a graded Al2O3-NiAl layer did not crack when given a thermal excursion to 1100 C and then cooled to room temperature. Acoustic emission tests on single fiber specimens were unsuccessful in detecting load drops associated with the successive fracture of the fiber. For the SCS-6/Ti3Al system this was the result of several factors which included the matrix/fiber ratio and poor bonding of the matrix and fiber. In the case of the Al2O3/NiAl system brittle failure of the NiAl matrix precluded fiber breakdown during tensile loading.

Merrick, Howard F.↗

NASA NTRS1986

Supercomputer implementation of finite element algorithms for high speed compressible flows

Prediction of compressible flow phenomena using the finite element method is of recent origin and considerable interest. Two shock capturing finite element formulations for high speed compressible flows are described. A Taylor-Galerkin formulation uses a Taylor series expansion in time coupled with a Galerkin weighted residual statement. The Taylor-Galerkin algorithms use explicit artificial dissipation, and the performance of three dissipation models are compared. A Petrov-Galerkin algorithm has as its basis the concepts of streamline upwinding. Vectorization strategies are developed to implement the finite element formulations on the NASA Langley VPS-32. The vectorization scheme results in finite element programs that use vectors of length of the order of the number of nodes or elements. The use of the vectorization procedure speeds up processing rates by over two orders of magnitude. The Taylor-Galerkin and Petrov-Galerkin algorithms are evaluated for 2D inviscid flows on criteria such as solution accuracy, shock resolution, computational speed and storage requirements. The convergence rates for both algorithms are enhanced by local time-stepping schemes. Extension of the vectorization procedure for predicting 2D viscous and 3D inviscid flows are demonstrated. Conclusions are drawn regarding the applicability of the finite element procedures for realistic problems that require hundreds of thousands of nodes.

Thornton, E. A.↗