Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “characteristic equations”

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 595 records · Page 33

Kinetic theory model for the flow of a simple gas from a three-dimensional axisymmetric nozzle

A system of nonlinear integral equations equivalent to the Krook kinetic equations for the steady state is the mathematical basis used to develop a computer code to model the flowfields for low-thrust three-dimensional axisymmetric nozzles. The method of characteristics is used to solve numerically by an iteration process the approximated Boltzmann equation for the number density, temperature, and velocity profiles of a simple gas as it expands into a vacuum. Results predict backscatter and show the effect of the nozzle wall boundary layer on the external flowfields.

Riley, B. R.↗

Euler/Navier-Stokes Solvers Applied to Ducted Fan Configurations

Due to noise considerations, ultra high bypass ducted fans have become a more viable design. These ducted fans typically consist of a rotor stage containing a wide chord fan and a stator stage. One of the concerns for this design is the classical flutter that keeps occurring in various unducted fan blade designs. These flutter are catastrophic and are to be avoided in the flight envelope of the engine. Some numerical investigations by Williams, Cho and Dalton, have suggested that a duct around a propeller makes it more unstable. This needs to be further investigated. In order to design an engine to safely perform a set of desired tasks, accurate information of the stresses on the blade during the entire cycle of blade motion is required. This requirement in turn demands that accurate knowledge of steady and unsteady blade loading be available. Aerodynamic solvers based on unsteady three-dimensional analysis will provide accurate and fast solutions and are best suited for aeroelastic analysis. The Euler solvers capture significant physics of the flowfield and are reasonably fast. An aerodynamic solver Ref. based on Euler equations had been developed under a separate grant from NASA Lewis in the past. Under the current grant, this solver has been modified to calculate the aeroelastic characteristics of unducted and ducted rotors. Even though, the aeroelastic solver based on three-dimensional Euler equations is computationally efficient, it is still very expensive to investigate the effects of multiple stages on the aeroelastic characteristics. In order to investigate the effects of multiple stages, a two-dimensional multi stage aeroelastic solver was also developed under this task, in collaboration with Dr. T. S. R. Reddy of the University of Toledo. Both of these solvers were applied to several test cases and validated against experimental data, where available.

Keith, Theo G., Jr.↗

Micrometeoroid and Orbital Debris Threat Assessment: Mars Sample Return Earth Entry Vehicle

This report provides results of a Micrometeoroid and Orbital Debris (MMOD) risk assessment of the Mars Sample Return Earth Entry Vehicle (MSR EEV). The assessment was performed using standard risk assessment methodology illustrated in Figure 1-1. Central to the process is the Bumper risk assessment code (Figure 1-2), which calculates the critical penetration risk based on geometry, shielding configurations and flight parameters. The assessment process begins by building a finite element model (FEM) of the spacecraft, which defines the size and shape of the spacecraft as well as the locations of the various shielding configurations. This model is built using the NX I-deas software package from Siemens PLM Software. The FEM is constructed using triangular and quadrilateral elements that define the outer shell of the spacecraft. Bumper-II uses the model file to determine the geometry of the spacecraft for the analysis. The next step of the process is to identify the ballistic limit characteristics for the various shield types. These ballistic limits define the critical size particle that will penetrate a shield at a given impact angle and impact velocity. When the finite element model is built, each individual element is assigned a property identifier (PID) to act as an index for its shielding properties. Using the ballistic limit equations (BLEs) built into the Bumper-II code, the shield characteristics are defined for each and every PID in the model. The final stage of the analysis is to determine the probability of no penetration (PNP) on the spacecraft. This is done using the micrometeoroid and orbital debris environment definitions that are built into the Bumper-II code. These engineering models take into account orbit inclination, altitude, attitude and analysis date in order to predict an impacting particle flux on the spacecraft. Using the geometry and shielding characteristics previously defined for the spacecraft and combining that information with the environment model calculations, the Bumper-II code calculates a probability of no penetration for the spacecraft.

Christiansen, Eric L.↗

Investigations into a new algorithm for calculating H infinity optimal controllers

A new algorithm for calculating H (sup infinity) optimal controllers is investigated. The algorithm is significantly simpler than existing approaches and yields much simpler controllers. The design equations are first presented. Special system transformations required to apply the algorithm are then presented. The use of the algorithm with sampled-data systems is outlined in detail. Several constraints on the characteristics of the problem formulation are required for the application of the design equations. The consequences of these constraints are investigated by applying the algorithm to a simplified design for a subsystem of a large space structure ground test facility. The investigation of these constraints is continued by application of the design equations and constraints to an extremely simple tracking problem. The result of these investigations is the development of a frequency dependent weighting strategy that allows realistic control problems to be cast in a form compatible with the new algorithm. Further work is indicated in the area of developing strategies for choosing frequency-dependent weights to achieve specific design goals. The use of the freedom in problem formulation to achieve robustness/performance tradeoffs should also be investigated. It is not clear that the algorithm always leads to simpler controllers. The more restrictive formulation may dictate that frequency-dependent weighting adds to the controller order disproportionately. This effect must also be investigated.

Irwin, R. Dennis↗

Investigation of boundary layer and turbulence characteristics inside the passages of an axial flow inducer

A study of the boundary layer and turbulence characteristics inside the passages of an axial flow inducer is reported. The first part deals with the analytical and experimental investigation of the boundary layer characteristics in a four bladed flat plate inducer passage operated with no throttle. An approximate analysis for the prediction of radial and chordwise velocity profiles across the passage is carried out. The momentum integral technique is used to predict the gross properties of the boundary layer. Equations are given for the exact analysis of the turbulent boundary layer characteristics using the turbulent field method. Detailed measurement of boundary layer profiles, limiting streamline angle and skin friction stress on the rotating blade is also reported. Part two of this report deals with the prediction of the flow as well as blade static pressure measurements in a three bladed inducer with cambered blades operated at a flow coefficient of 0.065. In addition, the mean velocity and turbulence measurements carried out inside the passage using a rotating triaxial probe is reported.

Anand, A.↗

Dynamic characteristics of large repetitive framelike structures

Using a building block approach and starting with a single element, expressions for the energy of various two-dimensional frametype gridwork configurations are derived. These are then used to develop energy equivalent continua for the gridworks. Equations of motion and associated boundary conditions are obtained for the continua. Some dynamic characteristics of these continua are investigated and compared with corresponding results obtained from finite element codes and also with some available theoretical predictions.

Nayfeh, A. H.↗

Analysis and design of numerical schemes for gas dynamics 1: Artificial diffusion, upwind biasing, limiters and their effect on accuracy and multigrid convergence

The theory of non-oscillatory scalar schemes is developed in this paper in terms of the local extremum diminishing (LED) principle that maxima should not increase and minima should not decrease. This principle can be used for multi-dimensional problems on both structured and unstructured meshes, while it is equivalent to the total variation diminishing (TVD) principle for one-dimensional problems. A new formulation of symmetric limited positive (SLIP) schemes is presented, which can be generalized to produce schemes with arbitrary high order of accuracy in regions where the solution contains no extrema, and which can also be implemented on multi-dimensional unstructured meshes. Systems of equations lead to waves traveling with distinct speeds and possibly in opposite directions. Alternative treatments using characteristic splitting and scalar diffusive fluxes are examined, together with modification of the scalar diffusion through the addition of pressure differences to the momentum equations to produce full upwinding in supersonic flow. This convective upwind and split pressure (CUSP) scheme exhibits very rapid convergence in multigrid calculations of transonic flow, and provides excellent shock resolution at very high Mach numbers.

Jameson, Antony↗

Solution of partial differential equations on vector and parallel computers

The present status of numerical methods for partial differential equations on vector and parallel computers was reviewed. The relevant aspects of these computers are discussed and a brief review of their development is included, with particular attention paid to those characteristics that influence algorithm selection. Both direct and iterative methods are given for elliptic equations as well as explicit and implicit methods for initial boundary value problems. The intent is to point out attractive methods as well as areas where this class of computer architecture cannot be fully utilized because of either hardware restrictions or the lack of adequate algorithms. Application areas utilizing these computers are briefly discussed.

Ortega, J. M.↗

Solution of partial differential equations on vector and parallel computers

The present status of numerical methods for partial differential equations on vector and parallel computers was reviewed. The relevant aspects of these computers are discussed and a brief review of their development is included, with particular attention paid to those characteristics that influence algorithm selection. Both direct and iterative methods are given for elliptic equations as well as explicit and implicit methods for initial boundary value problems. The intent is to point out attractive methods as well as areas where this class of computer architecture cannot be fully utilized because of either hardware restrictions or the lack of adequate algorithms. Application areas utilizing these computers are briefly discussed.

Ortega, J. M.↗

Navier-Stokes Simulation of the Canard-Wing-Body Longitudinal Dynamic Stability Characteristics

Many modern aircraft are canard-configured for aircraft control and improved aerodynamic performance. Canards can often enhance aircraft cruise performance, maneuverability and agility. For close-coupled canard configurations, the aerodynamic interaction between the canard and wing significantly changes the flow characteristics of the wing. In unsteady flow, such changes in the flow structure and performance of wings can be quite pronounced. Accurate modeling of the unsteady aerodynamics is essential for potential CFD design and analysis of such configurations. A time-accurate numerical simulation is performed to study the unsteady aerodynamic interaction between a canard and wing with emphasis on the effects of the canard on the configuration's dynamic response characteristics. The thin-layer Reynolds-averaged Navier-Stokes Equations with various turbulence models are used in this study. Computations are made on a generic, analytically-defined, close-coupled canard-wing-body configuration which has been the subject of numerous previously published experimental studies during the 1970's to mid-80's. More recently, a series of steady-flow simulations has been performed and published by the author. In the current study, the configuration is given prescribed ramp and oscillatory motions in order to predict characteristics such as the damping-in-pitch and oscillatory longitudinal stability parameters. The current computations are made at high-subsonic and transonic Mach numbers, moderate angles-of- attack from -4 to 20 degrees, and at various pitch rates and reduced frequencies. Comparisons of pressures and integrated force quantities (e.g. lift, drag, pitching moment and selected dynamic parameters) are made with other published computational results and available experimental data. Results showing the unsteady effects of the canard on surface pressures, integrated forces, canard-wing vortex interaction and vortex breakdown will be presented.

Tu, Eugene L.↗

Experimental Verification of the Rudder-Free Stability Theory for an Airplane Model Equipped with Rudders Having Negative Floating Tendency and Negligible Friction

An investigation has been made in the Langley free-flight tunnel to obtain an experimental verification of the theoretical rudder-free stability characteristics of an airplane model equipped with conventional rudders having negative floating tendencies and negligible friction. The model used in the tests was equipped with a conventional single vertical tail having rudder area 40 percent of the vertical tail area. The model was tested both in free flight and mounted on a strut that allowed freedom only in yaw. Tests were made with three different amounts of rudder aerodynamic balance and with various values of mass, moment of inertia, and center-of-gravity location of the rudder. Most of the stability derivatives required for the theoretical calculations were determined from forced and free-oscillation tests of the particular model tested. The theoretical analysis showed that the rudder-free motions of an airplane consist largely of two oscillatory modes - a long-period oscillation somewhat similar to the normal rudder-fixed oscillation and a short-period oscillation introduced only when the rudder is set free. It was found possible in the tests to create lateral instability of the rudder-free short-period mode by large values of rudder mass parameters even though the rudder-fixed condition was highly stable. The results of the tests and calculation indicated that for most present-day airplanes having rudders of negative floating tendency, the rudder-free stability characteristics may be examined by simply considering the dynamic lateral stability using the value of the directional-stability parameter Cn(sub p) for the rudder-free condition in the conventional controls-fixed lateral-stability equations. For very large airplanes having relatively high values of the rudder mass parameters with respect to the rudder aerodynamic parameters, however, analysis of the rudder-free stability should be made with the complete equations of motion. Good agreement between calculated and measured rudder-free stability characteristics was obtained by use of the general rudder-free stability theory, in which four degrees of lateral freedom are considered. When this assumption is made that the rolling motions alone or the lateral and rolling motions may be neglected in the calculations of rudder-free stability, it is possible to predict satisfactorily the characteristics of the long-period (Dutch roll type) rudder-free oscillation for airplanes only when the effective-dihedral angle is small. With these simplifying assumptions, however, satisfactory prediction of the short-period oscillation may be obtained for any dihedral. Further simplification of the theory based on the assumption that the rudder moment of inertia might be disregarded was found to be invalid because this assumption made it impossible to calculate the characteristics of the short-period oscillations.

STABILITY- FAIRCHILD R2K-1↗

A fast efficient implicit scheme for the gasdynamic equations using a matrix reduction technique

An efficient implicit finite-difference algorithm for the gasdynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved without loss of the stability characteristics generality found in the Beam and Warming approximate factorization algorithm. Steady-state solutions to the conservative Euler equations in generalized coordinates are obtained for transonic flows and used to show that the method offers computational advantages over the conventional Beam and Warming scheme. Existing Beam and Warming codes can be retrofit with minimal effort. The theoretical extension of the matrix reduction technique to the full Navier-Stokes equations in Cartesian coordinates is presented in detail. Linear stability, using a Fourier stability analysis, is demonstrated and discussed for the one-dimensional Euler equations.

Barth, T. J.↗

Determination of compressor in-stall characteristics from engine surge transients

A technique for extracting the in-stall pumping characteristics for an axial flow compressor operating in an engine system environment is developed. The technique utilizes a Hybrid computer simulation of the compressor momentum equation into which actual transient data are used to provide all terms but the desired compressor characteristic. The compressor force characteristic as a function of corrected flow and speed result from the computation. The critical problem of data filtering is addressed. Results for a compressor operating in a turbofan engine are presented and comparison is made with the conventional compressor map. The relationship of the compressor surge characteristic with its rotating stall characteristic is explored. Initial interpretation of the measured results is presented.

Lorenzo, C. F.↗

Determination of compressor in-stall characteristics from engine surge transients

A technique for extracting the in-stall pumping characteristics for an axial flow compressor operating in an engine system environment is developed. The technique utilzes a Hybrid computer simulation of the compressor momentum equation in which actual transient data are used to provide all terms but the desired compressor characteristic. The compressor force characteristic as a function of corrected flow and speed results from the computation. The critical problem of data filtering is addressed. Results for a compressor operating in a turbofan engine are presented and comparison is made with the conventional compressor map. The relationship of the compressor surge charcteristic with its rotating stall characteristic is explored. Initial interpretation of the measured results is presented. Previously announced in STAR as N84-22566

Lorenzo, C. F.↗

Long-term radiation effects on GaAs solar cell characteristics

This report investigates preliminary design considerations which should be considered for a space experiment involving Gallium Arsenide (GaAs) solar cells. The electron radiation effects on GaAs solar cells were conducted in a laboratory environment, and a statistical analysis of the data is presented. In order to augment the limited laboratory data, a theoretical investigation of the effect of radiation on GaAs solar cells is also developed. The results of this study are empirical prediction equations which can be used to estimate the actual damage of electrical characteristics in a space environment. The experimental and theoretical studies also indicate how GaAs solar cell parameters should be designed in order to withstand the effects of electron radiation damage.

Heinbockel, J. H.↗

Numerical determination of the transmissibility characteristics of a squeeze film damped forced vibration system

Numerical solutions of the governing equations of motion of a liquid squeeze film damped forced vibration system were carried out to examine the feasibility of using a liquid squeeze film to cushion and protect large structures, such as buildings, located in areas of high seismic activity. The mathematical model used was that for a single degree of freedom squeeze film damped spring mass system. The input disturbance was simulated by curve fitting actual seismic data with an eleventh order Lagranging polynomial technique. Only the normal component of the seismic input was considered. The nonlinear, nonhomogeneous governing differential equation of motion was solved numerically to determine the transmissibility over a wide range of physical parameters using a fourth-order Runge-Kutta technique. It is determined that a liquid squeeze film used as a damping agent in a spring-mass system can significantly reduce the response amplitude for a seismic input disturbance.

Sutton, M. A.↗

Exact and approximate solutions to the oblique shock equations for real-time applications

The derivation of exact solutions for determining the characteristics of an oblique shock wave in a supersonic flow is investigated. Specifically, an explicit expression for the oblique shock angle in terms of the free stream Mach number, the centerbody deflection angle, and the ratio of the specific heats, is derived. A simpler approximate solution is obtained and compared to the exact solution. The primary objectives of obtaining these solutions is to provide a fast algorithm that can run in a real time environment.

Hartley, T. T.↗