Search NASA⌕ Search

SEARCH · Search NASA

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Assessment of a One-Dimensional Finite Element Charring Ablation Material Response Model for Phenolic-Impregnated Carbon Ablator

In this study, mathematical formulations to model charring ablation problems were numerically implemented using finite element analysis (FEA) with ABAQUS, which account for the material decomposition and progressive surface removal in the heat conduction and the surface energy balance equations. FEA was performed for a one-dimensional model to predict the temperature and ablation histories of a phenolic-impregnated carbon ablator sample (i.e., a common heat shield material for hypersonic vehicles and spacecraft) subjected to oxy-acetylene torch flame (i.e., 0.8 SLPM acetylene gas to 2.7 SLPM oxygen gas). The recovery enthalpy and convective heat transfer coefficient for the ablation model were calculated based on gas compositions and two assumed surface conditions (i.e., equilibrium and frozen). Simulations using the calculated recovery enthalpy and convective heat transfer coefficient resulted in a recession rate of 6.38 times (equilibrium) and 14.08 times (frozen) higher than the experimental data, despite fair agreement of the surface temperature. In addition, the effect of the heat transfer coefficient was investigated through a steady-state ablation analysis. The results of the analysis indicate that there is not one single value for the heat transfer coefficient that would allow the prediction to match both measured recession rate and surface temperature. Possible reasons for such an inconsistency are provided and discussed.

Thermal Protection System↗
NASA NTRS1984

Surface cracks in a plate of finite width under tension or bending

The problem of a finite plate containing collinear surface cracks is considered and solved by using the line spring model with plane elasticity and Reissner's plate theory. The main focus is on the effect of interaction between two cracks or between cracks and stress-free plate boundaries on the stress intensity factors in an effort to provide extensive numerical results which may be useful in applications. Some sample results are obtained and are compared with the existing finite element results. Then the problem is solved for a single (internal) crack, two collinear cracks, and two corner cracks for wide range of relative dimensions. Particularly in corner cracks, the agreement with the finite element solution is surprisingly very good. The results are obtained for semi-elliptic and rectangular crack profiles which may, in practice, correspond to two limiting cases of the actual profile of a subcritically growing surface crack.

Erdogan, F.↗
DOE OSTI2025

Additional considerations in analytical solution for time-dependent heat conduction in a three-dimensional multilayer sphere

This work presents an analytical method to solve the heat conduction equation in three dimensions for problems consisting of multilayer concentric spheres. The method can be used to treat time-varying heat conduction problems where the heat source that drives the transient is time-invariant. Equally applicable to all Poisson-type problems with concentric spherical geometry, the method consists of representing the solution as a summation of weighted eigenfunctions. The weights for each eigenfunction are computed algebraically. Previous work has already established the core constituents of the methodology. The current work augments the existing methods by including consideration of nonzero interface resistance between layers and explicit discussion on the boundary condition homogenization required to treat inhomogeneous problems. Also, two demonstration problems are presented. One demonstration problem is based on the method of manufactured solutions and therefore allows for comparison with exact expressions for the solution temperature distribution. The second, more complex, demonstration problem relies on the finite element method for comparisons. The expected convergence behavior is observed for both demonstration problems.

97 - MATHEMATICS AND COMPUTING↗
NASA NTRS1991

Configuration optimization of space structures

The objective is to develop a computer aid for the conceptual/initial design of aerospace structures, allowing configurations and shape to be apriori design variables. The topics are presented in viewgraph form and include the following: Kikuchi's homogenization method; a classical shape design problem; homogenization method steps; a 3D mechanical component design example; forming a homogenized finite element; a 2D optimization problem; treatment of volume inequality constraint; algorithms for the volume inequality constraint; object function derivatives--taking advantage of design locality; stiffness variations; variations of potential; and schematics of the optimization problem.

Felippa, Carlos↗
NASA NTRS1990

A finite element solver for 3-D compressible viscous flows

Computation of the flow field inside a space shuttle main engine (SSME) requires the application of state of the art computational fluid dynamic (CFD) technology. Several computer codes are under development to solve 3-D flow through the hot gas manifold. Some algorithms were designed to solve the unsteady compressible Navier-Stokes equations, either by implicit or explicit factorization methods, using several hundred or thousands of time steps to reach a steady state solution. A new iterative algorithm is being developed for the solution of the implicit finite element equations without assembling global matrices. It is an efficient iteration scheme based on a modified nonlinear Gauss-Seidel iteration with symmetric sweeps. The algorithm is analyzed for a model equation and is shown to be unconditionally stable. Results from a series of test problems are presented. The finite element code was tested for couette flow, which is flow under a pressure gradient between two parallel plates in relative motion. Another problem that was solved is viscous laminar flow over a flat plate. The general 3-D finite element code was used to compute the flow in an axisymmetric turnaround duct at low Mach numbers.

Reddy, K. C.↗
NASA NTRS2010

Fully-Implicit Navier-Stokes (FIN-S)

FIN-S is a SUPG finite element code for flow problems under active development at NASA Lyndon B. Johnson Space Center and within PECOS: a) The code is built on top of the libMesh parallel, adaptive finite element library. b) The initial implementation of the code targeted supersonic/hypersonic laminar calorically perfect gas flows & conjugate heat transfer. c) Initial extension to thermochemical nonequilibrium about 9 months ago. d) The technologies in FIN-S have been enhanced through a strongly collaborative research effort with Sandia National Labs.

Kirk, Benjamin S.↗
NASA NTRS1993

Structural/control system optimization with variable actuator masses

A method is presented to integrate the design space for structural/control system optimization problems in the case of linear state feedback control. Nonstructural lumped masses and control system design variables as well as structural sizing variables are all treated equally as independent design variables in the optimization process. Structural and control design variable linking schemes are used in order to avoid a prohibitively large increase in the total number of independent design variables. When actuator masses are treated as nonstructural lumped mass design variables, special consideration is given to the relation between the transient peak responses and the required actuator masses which is formulated as a behavior constraint form. The original nonlinear mathematical programming problem based on a finite element formulation and linear state feedback is replaced by a sequence of explicit approximate problems exploiting various approximation concepts such as design variable linkings, temporary constraint deletion and first order Taylor series expansion of nonlinear behavior constraints in terms of intermediate design variables. Examples which involve a variety of dynamic behavior constraints (including constraints on closed-loop eigenvalues, peak transient displacements, peak actuator forces, and relations between the peak responses and the actuator masses) are effectively solved by using the method presented.

Jin, Ik M.↗
NASA NTRS1986

Cyclic creep analysis from elastic finite-element solutions

A uniaxial approach was developed for calculating cyclic creep and stress relaxation at the critical location of a structure subjected to cyclic thermomechanical loading. This approach was incorporated into a simplified analytical procedure for predicting the stress-strain history at a crack initiation site for life prediction purposes. An elastic finite-element solution for the problem was used as input for the simplified procedure. The creep analysis includes a self-adaptive time incrementing scheme. Cumulative creep is the sum of the initial creep, the recovery from the stress relaxation and the incremental creep. The simplified analysis was exercised for four cases involving a benchmark notched plate problem. Comparisons were made with elastic-plastic-creep solutions for these cases using the MARC nonlinear finite-element computer code.

Kaufman, A.↗
NASA NTRS1982

Application of the Finite Element Method to Rotary Wing Aeroelasticity

A finite element method for the spatial discretization of the dynamic equations of equilibrium governing rotary-wing aeroelastic problems is presented. Formulation of the finite element equations is based on weighted Galerkin residuals. This Galerkin finite element method reduces algebraic manipulative labor significantly, when compared to the application of the global Galerkin method in similar problems. The coupled flap-lag aeroelastic stability boundaries of hingeless helicopter rotor blades in hover are calculated. The linearized dynamic equations are reduced to the standard eigenvalue problem from which the aeroelastic stability boundaries are obtained. The convergence properties of the Galerkin finite element method are studied numerically by refining the discretization process. Results indicate that four or five elements suffice to capture the dynamics of the blade with the same accuracy as the global Galerkin method.

Straub, F. K.↗
NASA NTRS1984

Finite element modelling of buried structures

In many structures the final stress states are dependent on the sequence of construction or the stress states at various stages of construction are of interest. Such problems can be analyzed using finite element programs that have the capability of adding (birthing) elements to simulate the progress of construction. However, the usual procedure of assembling elements may lead to numerical instabilities or stress states that are unrealistic. Both problems are demonstrated in the analysis of a structure using the program ADINA. A technique which combines application of a preload with element birthing to overcome these problems is described and illustrated.

Playdon, D. K.↗
NASA NTRS1996

Towards and FVE-FAC Method for Determining Thermocapillary Effects on Weld Pool Shape

Several practical materials processes, e.g., welding, float-zone purification, and Czochralski crystal growth, involve a pool of molten metal with a free surface, with strong temperature gradients along the surface. In some cases, the resulting thermocapillary flow is vigorous enough to convect heat toward the edges of the pool, increasing the driving force in a sort of positive feedback. In this work we examine this mechanism and its effect on the solid-liquid interface through a model problem: a half space of pure substance with concentrated axisymmetric surface heating, where surface tension is strong enough to keep the liquid free surface flat. The numerical method proposed for this problem utilizes a finite volume element (FVE) discretization in cylindrical coordinates. Because of the axisymmetric nature of the model problem, the control volumes used are torroidal prisms, formed by taking a polygonal cross-section in the (r, z) plane and sweeping it completely around the z-axis. Conservation of energy (in the solid), and conservation of energy, momentum, and mass (in the liquid) are enforced globally by integrating these quantities and enforcing conservation over each control volume. Judicious application of the Divergence Theorem and Stokes' Theorem, combined with a Crank-Nicolson time-stepping scheme leads to an implicit algebraic system to be solved at each time step. It is known that near the boundary of the pool, that is, near the solid-liquid interface, the full conduction-convection solution will require extremely fine length scales to resolve the physical behavior of the system. Furthermore, this boundary moves as a function of time. Accordingly, we develop the foundation of an adaptive refinement scheme based on the principles of Fast Adaptive Composite Grid methods (FAC). Implementation of the method and numerical results will appear in a later report.

Canright, David↗
NASA NTRS1990

An error analysis of least-squares finite element method of velocity-pressure-vorticity formulation for Stokes problem

A theoretical proof of the optimal rate of convergence for the least-squares method is developed for the Stokes problem based on the velocity-pressure-vorticity formula. The 2D Stokes problem is analyzed to define the product space and its inner product, and the a priori estimates are derived to give the finite-element approximation. The least-squares method is found to converge at the optimal rate for equal-order interpolation.

Chang, Ching L.↗
NASA NTRS1985

Development of non-linear finite element computer code

Recent work has shown that the use of separable symmetric functions of the principal stretches can adequately describe the response of certain propellant materials and, further, that a data reduction scheme gives a convenient way of obtaining the values of the functions from experimental data. Based on representation of the energy, a computational scheme was developed that allows finite element analysis of boundary value problems of arbitrary shape and loading. The computational procedure was implemental in a three-dimensional finite element code, TEXLESP-S, which is documented herein.

Becker, E. B.↗
NASA NTRS1991

Solution of geometrically nonlinear statics problems by the p-version of the finite element method

This project is concerned with the possibility of using computers for the simulation of structural systems with the same degree of reliability as full scale physical experiments. Reliable numerical simulation will make it possible to reduce the costs of engineering and improve the quality of engineering decisions based on computed information. An error of idealization is an error between the actual physical quantities on which engineering decisions are based (e.g., maximum principal stress, first natural frequency, etc.) and the same data corresponding to the exact solution of the mathematical model. An error of discretization is an error between the quantities of interest corresponding to the exact and approximate solutions of a mathematical model. A high degree of reliability can be achieved in numerical simulation only if both the errors of idealization and errors of discretization can be shown to be small.

Szabo, Barna A.↗
NASA NTRS1992

Finite element method for optimal guidance of an advanced launch vehicle

A temporal finite element based on a mixed form of Hamilton's weak principle is summarized for optimal control problems. The resulting weak Hamiltonian finite element method is extended to allow for discontinuities in the states and/or discontinuities in the system equations. An extension of the formulation to allow for control inequality constraints is also presented. The formulation does not require element quadrature, and it produces a sparse system of nonlinear algebraic equations. To evaluate its feasibility for real-time guidance applications, this approach is applied to the trajectory optimization of a four-state, two-stage model with inequality constraints for an advanced launch vehicle. Numerical results for this model are presented and compared to results from a multiple-shooting code. The results show the accuracy and computational efficiency of the finite element method.

Hodges, Dewey H.↗
NASA NTRS1994

BUCKY instruction manual, version 3.3

The computer program BUCKY is a p-version finite element package for the solution of structural problems. The current version of BUCKY solves the 2-D plane stress, 3-D plane stress plasticity, 3-D axisymmetric, Mindlin and Kirchoff plate bending, and buckling problems. The p-version of the finite element method is a highly accurate version of the traditional finite element method. Example cases are presented to show the accuracy and application of BUCKY.

Smith, James P.↗
NASA NTRS1982

On nonlinear finite element analysis in single-, multi- and parallel-processors

Numerical solution of nonlinear equilibrium problems of structures by means of Newton-Raphson type iterations is reviewed. Each step of the iteration is shown to correspond to the solution of a linear problem, therefore the feasibility of the finite element method for nonlinear analysis is established. Organization and flow of data for various types of digital computers, such as single-processor/single-level memory, single-processor/two-level-memory, vector-processor/two-level-memory, and parallel-processors, with and without sub-structuring (i.e. partitioning) are given. The effect of the relative costs of computation, memory and data transfer on substructuring is shown. The idea of assigning comparable size substructures to parallel processors is exploited. Under Cholesky type factorization schemes, the efficiency of parallel processing is shown to decrease due to the occasional shared data, just as that due to the shared facilities.

Utku, S.↗