SEARCH · Engineering Papers

Results for “Finite”

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 703 records · Page 39

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Composite grid and finite-volume LU implicit scheme for turbine flow analysis

A composite grid was generated in an attempt to improve grid quality for a typical turbine blade with large camber in terms of mesh control, smoothness, and orthogonality. This composite grid consists of the C grid (or O grid) in the immediate vicinity of the blade and the H grid in the upstream region and in the middle of the blade passage between the C grids. It provides a good boundary layer resolution around the leading edge region for viscous calculation, has orthogonality at the blade surface and slope continuity at the C-H (or O-H) interface, and has flexibility in controlling the mesh distribution in the upstream region without using excessive grid points. This composite grid eliminates the undesirable qualities of a single grid when generated for a typical turbine geometry. A finite-volume lower-upper (LU) implicit schemes can be used in solving for the turbine flows on the composite grid. This grid has a special grid node that is connected to more than four neighboring nodes in two dimensions and to more than six nodes in three dimensions. But the finite-volume approach poses no problem at the special point because each interior cell has only four neighboring cells in two dimensions and only six cells in three dimensions. The finite-volume LU implicit scheme was demonstrated to be robust and efficient for both external and internal flows in a broad flow regime.

Choo, Yung K.↗

NASA NTRS1987

Computation of viscous blast wave solutions with an upwind finite volume method

A fully conservative, viscous, implicit, upwind, finite-volume scheme for the thin-layer Navier-Stokes equations is described with application to blast wave flow fields. In this scheme, shocks are captured without the oscillations typical of central differencing techniques and wave speeds are accurately predicted. The finite volume philosophy ensures conservation and since boundary conditions are also treated conservatively, accurate reflections of waves from surfaces are assured. Viscous terms in the governing equations are treated in a manner consistent with the finite volume philosophy, resulting in very accurate prediction of boundary layer quantities. Numerical results are presented for four viscous problems: a steady boundary layer, a shock-induced boundary layer, a blast wave/cylinder interaction and a blast wave/supersonic missile interaction. Comparisons of the results with an established boundary layer code, similarity solution, and experimental data show excellent agreement.

Molvik, Gregory A.↗

NASA NTRS1987

Quality assessment and control of finite element solutions

Status and some recent developments in the techniques for assessing the reliability of finite element solutions are summarized. Discussion focuses on a number of aspects including: the major types of errors in the finite element solutions; techniques used for a posteriori error estimation and the reliability of these estimators; the feedback and adaptive strategies for improving the finite element solutions; and postprocessing approaches used for improving the accuracy of stresses and other important engineering data. Also, future directions for research needed to make error estimation and adaptive movement practical are identified.

Noor, Ahmed K.↗

NASA NTRS1985

Finite state modelling of unsteady aerodynamics and its application to a rotor dynamic problem

The paper presents a method for formulating finite state unsteady aerodynamic models in the time domain from frequency domain unsteady aerodynamics. The method is based on recognizing that the lift deficiency function represents an aerodynamic transfer function and utilizes the Bode plot technique, used in control systems engineering, to construct approximation to the lift deficiency function. Indicial response functions for both fixed wing and rotary wing applications are obtained, using these finite state unsteady aerodynamic models. It is shown that the rotary wing indicial response function is oscillatory and thus it is fundamentally different when compared to the fixed wing indicial response function which is nonoscillatory. Certain aspects of the finite state aerodynamic model are demonstrated by applying it to the flapping dynamics of an articulated helicopter rotor blade. The influence of unsteady aerodynamics on the damping characteristics of the rotor is examined. The same problem is also treated by using a different unsteady aerodynamic model, namely dynamic inflow. Based on a comparison of the results obtained with these two unsteady aerodynamic models, useful conclusions are drawn regarding some fundamental features of these theories.

Friedmann, P. P.↗

NASA NTRS1987

Strain-energy-release rate analysis of the end-notched flexure specimen using the finite-element method

Two-dimensional finite-element analysis of the end-notched flexure specimen was performed using 8-node isoparametric, parabolic elements to evaluate compliance and mode II strain energy release rates, G sub II. The G sub II values were computed using two different techniques: the virtural crack-closure technique (VCCT) and the rate of change of compliance with crack length (compliance derivative method). The analysis was performed for various crack-length-to-semi-span (a/L) ratios ranging from 0.2 to 0.9. Three material systems representing a wide range of material properties were analyzed. The compliance and strain energy release rates of the specimen calculated with the present finite-element analysis agree very well with beam theory equations including transverse shear. The G sub II values calculated using the compliance derivative method compared extremely well with those calculated using the VCCT. The G sub II values obtained by the compliance derivative method using the top or bottom beam deflections agreed closely with each other. The strain energy release rates from a plane-stress analysis were higher than the plane-strain values by only a small percentage, indicating that either assumption may be used in the analysis. The G sub II values for one material system calculated from the finite-element analysis agreed with one solution in the literature and disagreed with the other solution in the literature.

Salpekar, S. A.↗

NASA NTRS1988

A velocity-pressure integrated, mixed interpolation, Galerkin finite element method for high Reynolds number laminar flows

A velocity-pressure integrated, mixed interpolation, Galerkin finite element method for the Navier-Stokes equations is presented. In the method, the velocity variables were interpolated using complete quadratic shape functions and the pressure was interpolated using linear shape functions. For the two dimensional case, the pressure is defined on a triangular element which is contained inside the complete biquadratic element for velocity variables; and for the three dimensional case, the pressure is defined on a tetrahedral element which is again contained inside the complete tri-quadratic element. Thus the pressure is discontinuous across the element boundaries. Example problems considered include: a cavity flow for Reynolds number of 400 through 10,000; a laminar backward facing step flow; and a laminar flow in a square duct of strong curvature. The computational results compared favorable with those of the finite difference methods as well as experimental data available. A finite elememt computer program for incompressible, laminar flows is presented.

Kim, Sang-Wook↗

NASA NTRS1988

A preliminary investigation of finite-element modeling for composite rotor blades

The results from an initial phase of an in-house study aimed at improving the dynamic and aerodynamic characteristics of composite rotor blades through the use of elastic couplings are presented. Large degree of freedom shell finite element models of an extension twist coupled composite tube were developed and analyzed using MSC/NASTRAN. An analysis employing a simplified beam finite element representation of the specimen with the equivalent engineering stiffness was additionally performed. Results from the shell finite element normal modes and frequency analysis were compared to those obtained experimentally, showing an agreement within 13 percent. There was appreciable degradation in the frequency prediction for the torsional mode, which is elastically coupled. This was due to the absence of off-diagonal coupling terms in the formulation of the equivalent engineering stiffness. Parametric studies of frequency variation due to small changes in ply orientation angle and ply thickness were also performed. Results showed linear frequency variations less than 2 percent per 1 degree variation in the ply orientation angle, and 1 percent per 0.0001 inch variation in the ply thickness.

Lake, Renee C.↗

NASA NTRS1987

On a finite element CFD algorithm for compressible, viscous and turbulent aerodynamic flows

This paper develops and analyses individual construction aspects of an efficient and accurate finite element algorithm for prediction of viscous and turbulent flow fields of impact in aerodynamics. The theoretical construction employs a Taylor weak statement (TWS) for coincident embedding of stability mechanisms within a classic Galerkin finite element formulation of semidiscrete approximation error orthogonalization. A wide variety of the stabilizing mechanisms of independently derived CFD algorithms are contained within the TWS theory. An implicit construction that meets the requirement of efficient convergence to steady state is developed. The theoretical asymptotic error estimates of the TWS finite element algorithm for supersonic and viscous boundary layer flows are verified. Application to a three-dimensional turbulent flow is cited.

Baker, A. J.↗

NASA NTRS1988

Numerical simulation of hypersonic inlet flows with equilibrium or finite rate chemistry

An efficient numerical program incorporated with comprehensive high temperature gas property models has been developed to simulate hypersonic inlet flows. The computer program employs an implicit lower-upper time marching scheme to solve the two-dimensional Navier-Stokes equations with variable thermodynamic and transport properties. Both finite-rate and local-equilibrium approaches are adopted in the chemical reaction model for dissociation and ionization of the inlet air. In the finite rate approach, eleven species equations coupled with fluid dynamic equations are solved simultaneously. In the local-equilibrium approach, instead of solving species equations, an efficient chemical equilibrium package has been developed and incorporated into the flow code to obtain chemical compositions directly. Gas properties for the reaction products species are calculated by methods of statistical mechanics and fit to a polynomial form for C(p). In the present study, since the chemical reaction time is comparable to the flow residence time, the local-equilibrium model underpredicts the temperature in the shock layer. Significant differences of predicted chemical compositions in shock layer between finite rate and local-equilibrium approaches have been observed.

Yu, Sheng-Tao↗

NASA NTRS1988

Adaptive p-version based finite element formulations for thermal modeling/analysis of structural configurations

Adaptive p-version based hierarchical finite element formulations in conjunction with a posteriori error estimation concepts are described with emphasis on applicability for thermal modeling/analysis of structural configurations. The basic concepts and formulations of hierarchical p-versiion finite element for thermal analysis are first described. A posteriori error estimation features are utililzed to steer the process of adaptive refinement. Several configurations comprised of one-dimensional structures are evaluated to validate the applicability of the proposed formulations and to demonstrate the potential of the p-version adaptive formulations for thermal modeling/analysis. The methodology offers potential and promises to be an attractive alternative to conventional finite element thermal modeling/analysis approaches.

Tamma, Kumar K.↗

NASA NTRS1988

The effects of finite rate chemical processes on high enthalpy nozzle performance - A comparison between SPARK and SEAGULL

The generalized chemistry version of the computer code SPARK is extended to include two higher-order numerical schemes, yielding fourth-order spatial accuracy for the inviscid terms. The new and old formulations are used to study the influences of finite rate chemical processes on nozzle performance. A determination is made of the computationally optimum reaction scheme for use in high-enthalpy nozzles. Finite rate calculations are compared with the frozen and equilibrium limits to assess the validity of each formulation. In addition, the finite rate SPARK results are compared with the constant ratio of specific heats (gamma) SEAGULL code, to determine its accuracy in variable gamma flow situations. Finally, the higher-order SPARK code is used to calculate nozzle flows having species stratification. Flame quenching occurs at low nozzle pressures, while for high pressures, significant burning continues in the nozzle.

Carpenter, M. H.↗

NASA NTRS1987

Solving finite element equations on concurrent computers

This paper discusses the development of a concurrent algorithm for the solution of systems of equations arising in finite element applications. The approach is based on a hybrid of direct elimination method and preconditioned conjugate iteration. Two different preconditioners are used; diagonal scaling and a concurrent implementation of incomplete LU factorization. First, an automatic procedure is used to partition the finite element mesh into sub-structures. The particular mesh partition is chosen to minimize an estimate of the cost for evaluating the solution using this algorithm on a concurrent computer. These procedures are implemented in a finite element program on the JPL/CalTech MARK III hypercube computer. An overview of the structure of this program is presented. The performance of the solution method is demonstrated with the aid of a number of numerical test runs, and its advantages for concurrent implementations are discussed. Efficiency and speed-up factors over sequential machines for the numerical examples are highlighted.

Nour-Omid, B.↗

NASA NTRS1987

Optimal mapping of irregular finite element domains to parallel processors

Mapping the solution domain of n-finite elements into N-subdomains that may be processed in parallel by N-processors is an optimal one if the subdomain decomposition results in a well-balanced workload distribution among the processors. The problem is discussed in the context of irregular finite element domains as an important aspect of the efficient utilization of the capabilities of emerging multiprocessor computers. Finding the optimal mapping is an intractable combinatorial optimization problem, for which a satisfactory approximate solution is obtained here by analogy to a method used in statistical mechanics for simulating the annealing process in solids. The simulated annealing analogy and algorithm are described, and numerical results are given for mapping an irregular two-dimensional finite element domain containing a singularity onto the Hypercube computer.

Flower, J.↗

NASA NTRS1987

Parallel processing in finite element structural analysis

A brief review is made of the fundamental concepts and basic issues of parallel processing. Discussion focuses on parallel numerical algorithms, performance evaluation of machines and algorithms, and parallelism in finite element computations. A computational strategy is proposed for maximizing the degree of parallelism at different levels of the finite element analysis process including: 1) formulation level (through the use of mixed finite element models); 2) analysis level (through additive decomposition of the different arrays in the governing equations into the contributions to a symmetrized response plus correction terms); 3) numerical algorithm level (through the use of operator splitting techniques and application of iterative processes); and 4) implementation level (through the effective combination of vectorization, multitasking and microtasking, whenever available).

Noor, Ahmed K.↗

NASA NTRS1988

Some effects of finite spatial resolution on skin friction measurements in turbulent boundary layers

The effects of finite spatial resolution often cause serious errors in measurements in turbulent boundary layers, with particularly large effects for measurements of fluctuating skin friction and velocities within the sublayer. However, classical analyses of finite spatial resolution effects have generally not accounted for the substantial inhomogeneity and anisotropy of near-wall turbulence. The present study has made use of results from recent computational simulations of wall-bounded turbulent flows to examine spatial resolution effects for measurements made at a wall using both single-sensor probes and those employing two sensing volumes in a V shape. Results are presented to show the effects of finite spatial resolution on a variety of quantitites deduced from the skin friction field.

Westphal, Russell V.↗

NASA NTRS1988

A preliminary investigation of finite-element modeling for composite rotor blades

Lake, Renee C.↗

NASA NTRS1985

Finite element analysis of thermal barrier coatings

The near-term objective is to develop an understanding of the states of stresses and strains in a Zirconia-yttria thermal barrier coating (TBC) experiencing a given temperature drop. Results so obtained are expected to facilitate experimental work. In order to gain realistic insights into the distribution of stresses and strains in a complex TBC, the finite element approach was selected to model a cylindrical TBC specimen. Experimental evidence reported in the literature indicated the presence of rough interface between the ceramic coat and bond coat. Oxidation of the bond coat at ceramic-bond interface was observed, as was a small amount of cracking in the ceramics near the ceramic-bond interface. To account for these complex features, a plane-strain finite element computer program known as TBCOC was developed, taking advantage of a generic computer code known as MARC. This generic code was made available through the use of a supercomputer (Cray 1). The TBCOC model contains 1316 nodal points and 2140 finite elements. It is capable of a uniform isothermal loading. Results of a sample computer run are presented.

Chang, G. C.↗

NASA NTRS1989

Global functions in global-local finite-element analysis of localized stresses in prismatic structures

An important consideration in the global local finite-element method (GLFEM) is the availability of global functions for the given problem. The role and mathematical requirements of these global functions in a GLFEM analysis of localized stress states in prismatic structures are discussed. A method is described for determining these global functions. Underlying this method are theorems due to Toupin and Knowles on strain energy decay rates, which are related to a quantitative expression of Saint-Venant's principle. It is mentioned that a mathematically complete set of global functions can be generated, so that any arbitrary interface condition between the finite element and global subregions can be represented. Convergence to the true behavior can be achieved with increasing global functions and finite-element degrees of freedom. Specific attention is devoted to mathematically two-dimensional and three-dimensional prismatic structures. Comments are offered on the GLFEM analysis of NASA flat panel with a discontinuous stiffener. Methods for determining global functions for other effects are also indicated, such as steady-state dynamics and bodies under initial stress.

Dong, Stanley B.↗