SEARCH · Engineering Papers
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.
Galerkin/least-squares finite-element methods are presented for advective-diffusive equations. Galerkin/least-squares represents a conceptual simplification of streamline-upwind Petrov-Galerkin methods, and is in fact applicable to a wide variety of other problem types. A convergence analysis and error estimates are presented. Some numerical results for compressible Navier-Stokes flows are presented.
New error estimation techniques are described for two-dimensional convection-dominated flows. Examples of applications to linear convection problems are given. The methods are based on discontinuous finite-element techniques. Results show that conventional interpolation methods may lead to gross underestimates of actual error and be inadequate for use in adaptive methods.
The adaptive implicit-explicit (AIE) approach is presented for the finite-element solution of various problems in computational fluid mechanics. In the AIE approach, the elements are dynamically (adaptively) arranged into differently treated groups. The differences in treatment could be based on considerations such as the cost efficiency, the type of spatial or temporal discretization employed, the choice of field equations, etc. Several numerical tests are performed to demonstrate that this approach can achieve substantial savings in CPU time and memory.
The principles and operation of MMAX, a symbolic-computation program which automates the process of generating property matrices for structural materials, are briefly described and illustrated with sample analyses of a rubberlike material and an elastoplastic material. MMAX is written in LISP under the symbolic finite-element generator FINGER and the general symbolic manipulator MACSYMA; it first derives the formulas required by mathematical manipulation, and then translates the formulas into FORTRAN code, adapted to the particular type of machine to be used for the numerical calculations. This approach is shown to combine efficiently the advantages of symbolic and numerical computation for engineering applications.
The methodology used by the General Rotorcraft Aeromechanical Stability Program (GRASP) to model the characteristics of the flow through a helicopter rotor in hovering or axial flight is described. Since the induced flow plays a significant role in determining the aeroelastic properties of rotorcraft, the computation of the induced flow is an important aspect of the program. Because of the combined finite-element/multibody methodology used as the basis for GRASP, the implementation of induced velocity calculations presented an unusual challenge to the developers. To preserve the modelling flexibility and generality of the code, it was necessary to depart from the traditional methods of computing the induced velocity. This is accomplished by calculating the actuator disk contributions to the rotor loads in a separate element called the air mass element, and then performing the calculations of the aerodynamic forces on individual blade elements within the aeroelastic beam element.
Retrofitment of commercial aircraft with propfans could introduce undesirable aerodynamic sources of structure-borne noise that are absent for current turbojet powerplants. This paper theoretically examines the whipping action of the vortex wake from a generic propeller on the downstream rigid wing that supports it. The model addresses the high-frequency/compressible regime of most anticipated propfan implementations and produces an analytic solution for the distributed wing airload due to the periodic vortex impingement. The analysis also yields an expression for the local unsteady lift obtained from integration over an arbitrary internal patch of wing surface, for the purpose of applying a practical number of such forces at the nodes of a finite-element model for the corresponding structure (wing response results are not included in the present paper). Reported estimates of induced wing loads for a conventional-propeller example of demonstration appear to be in the reasonable range of expectation.
A novel method for computing the frequency-domain electromagnetic fields scattered from, and penetrating into, arbitrarily shaped, three-dimensional, lossy, inhomogeneous anisotropic scatters is presented. The method is based on a general volume integrodifferential formulation of the scattering problem and consists of the numerical solution of the coupled integral equations by the moment method and point matching. The numerical model of the scatterer is obtained by parametric volume elements, and the basis functions used to represent the field within each element are the same used in the finite-element method. Element integration problems due to the singular kernel of the integral equations are treated in some detail. Numerical results for both the isotropic and the anisotropic spherical scatterer are presented, including comparisons with results obtained by different numerical methods for the isotropic cases considered. The capability of the numerical code presented to deal with cases where the material parameters of the scatterer are given by singular matrices is discussed for two particular examples.
The results of the third year of work on a program which is part of the NASA Hot Section Technology program (HOST) are presented. The goals of this program are: (1) the development of unified constitutive models for rate dependent isotropic materials; and (2) the demonstration of the use of unified models in structural analyses of hot section components of gas turbine engines. The unified models selected for development and evaluation are those of Bodner-Partom and of Walker. A test procedure was developed for assisting the generation of a data base for the Bodner-Partom model using a relatively small number of specimens. This test procedure involved performing a tensile test at a temperature of interest that involves a succession of strain-rate changes. The results for B1900+Hf indicate that material constants related to hardening and thermal recovery can be obtained on the basis of such a procedure. Strain aging, thermal recovery, and unexpected material variations, however, preluded an accurate determination of the strain-rate sensitivity parameter is this exercise. The effects of casting grain size on the constitutive behavior of B1900+Hf were studied and no particular grain size effect was observed. A systematic procedure was also developed for determining the material constants in the Bodner-Partom model. Both the new test procedure and the method for determining material constants were applied to the alternate material, Mar-M247 . Test data including tensile, creep, cyclic and nonproportional biaxial (tension/torsion) loading were collected. Good correlations were obtained between the Bodner-Partom model and experiments. A literature survey was conducted to assess the effects of thermal history on the constitutive behavior of metals. Thermal history effects are expected to be present at temperature regimes where strain aging and change of microstructure are important. Possible modifications to the Bodner-Partom model to account for these effects are outlined. The use of a unified constitutive model for hot section component analyses was demonstrated by applying the Walker model and the MARC finite-element code to a B1900+Hf airfoil problem.
An upwind finite-element technique that uses cell-centered quantities and implicit and/or explicit time marching has been developed for computing hypersonic laminar viscous flows using adaptive unstructured grids in two and three dimensions. A perfect gas model as well as an equilibrium air model is implemented for solving high-speed flows. A first-order basic scheme and a higher-order flux-corrected transport (FCT) scheme have been implemented. This technique has been used to predict 'Type III and IV' shock interactions on a cylinder in two dimensions and a swept cylinder in three dimensions, with a view to determine the pressure and heating rate augmentation caused by an impinging shock on the leading edge of a cowl lip of an engine inlet. The predictions of wall pressure and heating rates compare very well with experimental data. The flow features are very distinctly captured with a sequence of adaptively-generated grids. Three-dimensional corner flow, typically encountered in engine inlets due to compression of the flow by ramps in the walls, is also modeled. This procedure is the first step in developing an integrated fluid, thermal, structural analysis capability for hypersonic flight vehicles like the National Aero-Space Plane.
The application of Galerkin/least-squares methodology to computational fluid dynamics is reviewed. Applications to advective-diffusive systems of equations and incompressible and compressible Euler and Navier-Stokes equations are described. Space-time finite-element methodology for time-dependent problems is emphasized and the status and role of so-called 'discontinuity-capturing operators' is reviewed.
Six methods for calculating shape sensitivity derivatives are compared for a two-material beam problem with a moving interface. It is found that as the finite-element mesh is refined, displacement sensitivity derivatives converge more slowly than the displacement themselves. Five of the methods agree fairly well, but the adjoint variational surface method provides substantially different results. However, the difference is found to reflect convergence from another direction to the same answer, rather than reduced accuracy.
The technology for controlling elastic deformations of flexible structures is one of the key considerations for future space initiatives. A vital area needed to achieve this objective is the development of a control design methodology applicable to future structures. The mu-synthesis technique is employed to design a high-performance vibration attenuation controller for the JPL/AFAL experimental flexible antenna structure. The results presented deal primarily with the control of first two global flexible modes using only two hub actuators and two hub sensors. Implementation of the multivariable control laws based on a finite-element model is presented. All results are from actual implementation on the JPL/AFAL flexible structure testbed.
The present conference on the dynamics and control methods of multibody robotic systems applicable to outer space applications first gives attention to such topics in multibody dynamics and control as the substructure synthesis approach to flexible multibody system control, an elastic finite-element kinematic analysis of a multirigid body, computational methods for high speed vehicles on flexible guideways, and multiarm coordination and control. Attention is then given to such issues in space manipulator dynamics and control methods as the global navigation of a free-flying space robot, dynamic singularities in the control of free-floating space manipulators, experimental results for compact space robot actuator control, and a 17-degree-of-freedom dexterous manipulator.
An elastic-plastic finite-element analysis with a critical crack-tip-opening displacement criterion was used to simulate fracture of various size compact and bend specimens made of HY-130 steel. From the calculated load-crack-extension and load-displacement curves, J-resistance (J-R) curves were determined by several methods. The simulated 3-R curves were insensitive to specimen size up to maximum load but were sensitive to specimen configuration for crack extensions greater than 10 percent of the initial uncracked ligament length.
A summary of the second phase of a design study fo a cycling spacecraft for regular and frequency personnel transport between Earth and a Mars base is presented. The first-phase study, called project CAMELOT, was a mission analysis and configuration design and was completed by the Space Systems Design class at the University of Michigan in the winter term, 1987. In the second phase, the configuration and component design has been developed by the Space Systems Design class in the winter term, 1988, beyond the initial configuration study. Analysis has led to numerous design modifications and improvements, as well as providing more detailed system definition. Major components of the second phase study include innovative detailed designs such as an interface between the rotating and non-rotating portions of the spacecraft, an elevator system that changes orientation to eliminate Coriolis forces, an electromagnetic radiation shield that significantly reduces the mass of the spacecraft, and a detailed assembly sequence and cost analysis. In addition, modifications to the propulsion and the power generation subsystems, taxi-docking facilities, and internal layout design have led to a more efficient, more reliable, and more comfortable vehicle. Finite-element stress and dynamic analysis has verified much of the system and indicates the validity of the design. Computer simulation of many moving components, including torus rotation, elevator and interface operations, and attitude control mechanisms further validate and support the various systems.
The objective of system identification is to determine reliable dynamical models of a structure by systematically using its measured excitation and response. It brings together in an integrated fashion, experimental, analytical, and computational techniques in structural dynamics. Areas of application for system identification include the following: (1) Model Evaluation--assessing assumptions (linearity and equivalent viscous damping) and techniques (finite-element modeling) used to construct theoretical models of a structure; (2) Model Improvement--updating of a theoretical model to enable more accurate response predictions for possible future loads on the structure, or for control of the structure; (3) Empirical Modelling--developing empirical relationships (nonlinear models) or empirical parameter values (modal damping) because the present state of the art does not provide theoretical results; and (4) Damage Detection and Assessment--continual or episodic updating of a structural model through vibration monitoring to detect and locate any structural damage. It can be argued that since the construction or modification of models using test data is subject to inherent uncertainties, the above problems should be properly treated within a Bayesian probabilistic framework. Such a methodology is presented which allows the precision of the estimates of the model parameters to be computed. It also leads to a guiding principle in applications. Namely, when selecting a single model from a given class of models, one should take the most probable model in the class based on the experimental data. Practical applications of this principle are given which are based on the utilization of measured seismic motions in large civil structures. Examples include the application of a computer program MODE-ID to identify modal properties directly from seismic excitation and response time histories from a nine-story steel-frame building at JPL and from a freeway overpass bridge.
The influence of pretwist, nonuniformities in mass and flexural stiffness, rotatory inertia and shear deformation on the natural frequencies of intact bones is evaluated by means of a linear elastic, finite-element model which has been programmed for solution on the digital computer. Theoretical results are compared to the results on the forced vibration of intact canine radii obtained experimentally by Thompson. Surprisingly, inclusion of fairly large pretwist angles (from -14 to 12 deg for one specimen) had little affect on the first three frequencies of transverse vibration in either the cranial or lateral directions. Inclusion of shear deformation reduced the third-mode frequency in the stiffest (lateral) direction by about six percent, otherwise shear deformation played a minor role in determining natural frequencies. Similarly. rotatory inertia had negligible influence up to the third natural frequency. The predominant influence on the first three natural frequencies of transverse vibration could be attributed to the variations in mass and flexural stiffness along the length of the test specimens. Different effective moduli of elasticity are required to yield correct absolute values for the frequencies which correspond to experimental findings. thus implying the presence of some inhomogeneities in material properties around the bone cross-section and/or along its length.
Observations from finite-element computer models, together with analytic developments based on percolation theory have suggested that subtle fluctuations of ECG morphology might serve as an indicator diminished cardiac electrical stability. With fixed-rate atrial pacing in canines, we have previously observed a pattern of alternation in T wave energy which correlated with cardiac electrical stability. We report here on a series of 20 canine experiments in which cardiac electrical stability (measured via Ventricular Fibrillation Threshold determination) was compared to a non-degenerate, multidimensional measurement of the degree of alternating activity present in the ECG complex morphology. The decrease in cardiac electrical stability brought on by both coronary artery occlusion and systemic hypothermia was consistently accompanied by subtle alternation in ECG morphology, with the absolute degree of alternating activity being significantly (negatively) correlated with cardiac electrical stability.