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 631 records · Page 35

Trajectory characteristics and heating of hypervelocity projectiles having large ballistic coefficients

A simple, approximate equation describing the velocity-density relationship (or velocity-altitude) has been derived from the flight of large ballistic coefficient projectiles launched at high speeds. The calculations obtained by using the approximate equation compared well with results for numerical integrations of the exact equations of motion. The flightpath equation was used to parametrically calculate maximum body decelerations and stagnation pressures for initial velocities from 2 to 6 km/s. Expressions were derived for the stagnation-point convective heating rates and total heat loads. The stagnation-point heating was parametrically calculated for a nonablating wall and an ablating carbon surface. Although the heating rates were very high, the pulse decayed quickly. The total nose-region heat shield weight was conservatively estimated to be only about 1 percent of the body mass.

Tauber, Michael E.↗

Compact high-order schemes for the Euler equations

An implicit approximate factorization (AF) algorithm is constructed which has the following characteristics. In 2-D: the scheme is unconditionally stable, has a 3 x 3 stencil and at steady state has a fourth order spatial accuracy. The temporal evolution is time accurate either to first or second order through choice of parameter. In 3-D: the scheme has almost the same properties as in 2-D except that it is now only conditionally stable, with the stability condition (the CFL number) being dependent on the cell aspect ratios, delta y/delta x and delta z/delta x. The stencil is still compact and fourth order accuracy at steady state is maintained.

Abarbanel, Saul↗

Three-dimensional boundary layer calculation by a characteristic method

A numerical method for solving the three-dimensional boundary layer equations for bodies of arbitrary shape is presented. In laminar flows, the application domain extends from incompressible to hypersonic flows with the assumption of chemical equilibrium. For turbulent boundary layers, the application domain is limited by the validity of the mixing length model used. In order to respect the hyperbolic nature of the equations reduced to first order partial derivative terms, the momentum equations are discretized along the local streamlines using of the osculator tangent plane at each node of the body fitted coordinate system. With this original approach, it is possible to overcome the use of the generalized coordinates, and therefore, it is not necessary to impose an extra hypothesis about the regularity of the mesh in which the boundary conditions are given. By doing so, it is possible to limit, and sometimes to suppress, the pre-treatment of the data coming from an inviscid calculation. Although the proposed scheme is only semi-implicit, the method remains numerically very efficient.

Houdeville, R.↗

Modeling of aircraft unsteady aerodynamic characteristics. Part 1: Postulated models

A short theoretical study of aircraft aerodynamic model equations with unsteady effects is presented. The aerodynamic forces and moments are expressed in terms of indicial functions or internal state variables. The first representation leads to aircraft integro-differential equations of motion; the second preserves the state-space form of the model equations. The formulations of unsteady aerodynamics is applied in two examples. The first example deals with a one-degree-of-freedom harmonic motion about one of the aircraft body axes. In the second example, the equations for longitudinal short-period motion are developed. In these examples, only linear aerodynamic terms are considered. The indicial functions are postulated as simple exponentials and the internal state variables are governed by linear, time-invariant, first-order differential equations. It is shown that both approaches to the modeling of unsteady aerodynamics lead to identical models.

Klein, Vladislav↗

Wave theory of turbulence in compressible media (acoustic theory of turbulence)

The generation and the transmission of sound in turbulent flows are treated as one of the several aspects of wave propagation in turbulence. Fluid fluctuations are decomposed into orthogonal Fourier components, with five interacting modes of wave propagation: two vorticity modes, one entropy mode, and two acoustic modes. Wave interactions, governed by the inhomogeneous and nonlinear terms of the perturbed Navier-Stokes equations, are modeled by random functions which give the rates of change of wave amplitudes equal to the averaged interaction terms. The statistical framework adopted is a quantum-like formulation in terms of complex distribution functions. The spatial probability distributions are given by the squares of the absolute values of the complex characteristic functions. This formulation results in nonlinear diffusion-type transport equations for the probability densities of the five modes of wave propagation.

Kentzer, C. P.↗

Structural Dynamics, Stability, and Control of Helicopters

The dynamic synthesis of gyroscopic structures consisting of point-connected substructures is investigated. The objective is to develop a mathematical model capable of an adequate simulation of the modal characteristics of a helicopter using a minimum number of degrees of freedom. The basic approach is to regard the helicopter structure as an assemblage of flexible substructures. The variational equations for the perturbed motion about certain equilibrium solutions are derived. The discretized variational equations can be conveniently exhibited in matrix form, and a great deal of information about the system modal characteristics can be extracted from the coefficient matrices. The derivation of the variational equations requires a monumental amount of algebraic operations. To automate this task a symbolic manipulation program on a digital computer is developed.

Meirovitch, L.↗

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↗