ALBERT

Description: The benefits of using a hypersonic waverider for spacecraft trajectory modification are presented. A waverider is a hypersonic vehicle specifically designed so that the undersurface bow shock is attached to the leading edge, which provides for the highest known lift-to-drag ratios achievable at high Mach number flight. Several viable space missions are suggested which could use such configurations for low-drag aero-assisted maneuvers in planetary atmospheres. It is shown that large changes in the spacecraft velocity vector can be accomplished with acceptably small losses in energy due to drag using a waverider aeroshell. The primary advantage of an aero-assist maneuver is suggested by comparison to a traditional gravity-assist trajectory. Some scaling laws are presented for comparing waveriders designed for different planetary atmospheres, and it is shown that the compositional differences between the terrestrial planets has a minimal impact on waverider design.

Description: Unsteady flow behavior and load characteristics of a 2D VR-7 airfoil with and without a leading-edge slat were studied in the water tunnel of the Aeroflightdynamics Directorate, NASA Ames Research Center. Both airfoils were oscillated sinusoidally between 5 and 25 deg at Re = 200,000 to obtain the unsteady lift, drag, and pitching moment data. A fluorescent dye was released from an orifice located at the leading edge of the airfoil for the purpose of visualizing the boundary layer and wake flow. The flowfield and load predictions of an incompressible Navier-Stokes code based on a velocity-vorticity formulation were compared with the test data. The test and predictions both confirm that the slatted VR-7 airfoil delays both static and dynamic stall as compared to the VR-7 airfoil alone.

Description: Analytic expressions have been derived for estimating the nonablating laminar and turbulent boundary-layer convective heating rates on inclined flat surfaces for the Martian atmosphere in thermochemical equilibrium. The equations are valid in the speed and altitude regime where aerobraking would occur at Mars. Comparisons with limited experimental measurements and calculations for CO2 (the Martian atmosphere is 95.6 percent CO2) yielded reasonably good agreement, especially for the ratios of heating rates in CO2 to those in air at the same conditions. In the aerobraking speed regime, the laminar flat surface boundary layer heating rates are 15-25 percent greater at Mars than in air. The differences between the turbulent heating rates are even more pronounced. The turbulent heating rates can be over 50 percent greater at Mars than in air at the same flight conditions.

Description: An incompressible three-dimensional laminar boundary-layer flow over a swept wing is used as a model to study both the wall-curvature and streamline-curvature effects on the stationary crossflow instability. The basic state is obtained by solving the full Navier-Stokes (N-S) equations numerically. The linear disturbance equations are cast on a fixed, body-intrinsic, curvilinear coordinate system. Those nonparallel terms which contribute mainly to the streamline-curvature effect are retained in the formulation of the disturbance equations and approximated by their local finite difference values. The resulting eigenvalue problem is solved by a Chebyshev collocation method. The present results indicate that the convex wall curvature has a stabilizing effect, whereas the streamline curvature has a destabilizing effect. A validation of these effects with an N-S solution for the linear disturbance flow is provided.

Description: The paper presents a general approach to constructing mean velocity profiles for compressible turbulent boundary layers with isothermal or adiabatic walls. The theory is based on a density-weighted transformation that allows the extension of the incompressible similarity laws of the wall to the compressible regions. The velocity profile family is compared to a range of experimental data, and excellent agreement is obtained. A self-consistent skin friction law, which satisfies the proposed velocity profile family, is derived and compared with the well-known Van Driest II theory for boundary layers in zero pressure gradient. The results are found to be at least as good as those obtained by using the Van Driest II transformation.

Description: A number of chemical-kinetic problems related to phenomena occurring behind a shock wave surrounding an object flying in the earth atmosphere are discussed, including the nonequilibrium thermochemical relaxation phenomena occurring behind a shock wave surrounding the flying object, problems related to aerobraking maneuver, the radiation phenomena for shock velocities of up to 12 km/sec, and the determination of rate coefficients for ionization reactions and associated electron-impact ionization reactions. Results of experiments are presented in form of graphs and tables, giving data on the reaction rate coefficients for air, the ionization distances, thermodynamic properties behind a shock wave, radiative heat flux calculations, Damkoehler numbers for the ablation-product layer, together with conclusions.

Description: The development of implicit upwind algorithms for the solution of the three-dimensional, time-dependent Euler equations on unstructured tetrahedral meshes is described. The implicit temporal discretization involves either a two-sweep Gauss-Seide relaxation procedure, a two-sweep Point-Jacobi relaxation procedure, or a single-sweep Point-Implicit procedure; the upwind spatial discretization is based on the flux-difference splitting of Roe. Detailed descriptions of the three implicit solution algorithms are given, and calculations for the Boeing 747 transport configuration are presented to demonstrate the algorithms. Advantages and disadvantages of the implicit algorithms are discussed. A steady-state solution for the 747 configuration, obtained at transonic flow conditions using a mesh of over 100,000 cells, required less than one hour of CPU time on a Cray-2 computer, thus demonstrating the speed and robustness of the general capability.

Description: Two-dimensional CFD analyses are presented related to the ground testing of hypersonic, air-breathing models which feature scramjet exhaust flow simulation. CFD analysis indicates that it is possible to test aftbody powered hypersonic airbreather configurations in a static, pumped-down environment to obtain aftbody aerodynamic performance data.

Description: Experimental studies have been conducted to assess Reynolds and Mach number effects on a supercritical multielement airfoil. The airfoil is representative of the stall-critical station of an advanced transport wing design. The experimental work was conducted as part of a cooperative program between the Douglas Aircraft Company and the NASA Langley Research Center to improve current knowledge of high-lift flows and to develop a validation data base with practical geometries/conditions for emerging computational methods. This article describes results obtained for both landing and takeoff multielement airfoils (four- and three-element configurations) for a variety of Mach/Reynolds number combinations up to flight conditions. Effects on maximum lift are considered for the landing configurations, and effects on both lift and drag are reported for the takeoff geometry. The present test results revealed considerable maximum lift effects on the three-element landing configuration for Reynolds number variations, and significant Mach number effects on the four-element airfoil.

Description: This paper presents an alternate approach for the generation of volumetric grids for supersonic and hypersonic flows about complex configurations. The method uses parametric two-dimensional block-face grid definitions, within the framework of GRIDGEN2D. The incorporation of a face decomposition reduces complex surfaces to simple shapes. These simple shapes are recombined to obtain the final face definition. The advantages of this method include the reduction of overall grid generation time through the use of vectorized computer code, the elimination of generating matching block faces, and simplified boundary conditions.

Description: An algorithm is described that calculates inviscid, laminar, and turbulent flows on triangular meshes with an upwind discretization. A brief description of the base solver and the multigrid implementation is given, followed by results that consist mainly of convergence rates for inviscid and viscous flows over a NACA four-digit airfoil section. The results show that multigrid does accelerate convergence when the same relaxation parameters that yield good single-grid performance are used; however, larger gains in performance can be realized by doing less work in the relaxation scheme.

Description: The interaction of a turbulent boundary layer with rarefaction and shock waves in flows past inclined steps was investigated experimentally and theoretically. The experiments were carried out in supersonic wind tunnels with test sections of 0.6 x 0.6 m and 0.2 x 0.2 m, with adiabatic conditions on the surface model. Based on detailed measurements of pressure fields, velocity, and surface friction, combined with results of optical visualization and analysis of limiting flow lines, three characteristic flow regimes are identified. These are (1) nonseparated flow at small step angles, (2) formation of a local separation zone with a free separation point at moderate step angles, and (3) formation of a large-scale separated flow with a fixed separation point at sufficiently large step angles. The quantitative gasdynamic schemes and surface pressure distributions for these flow regimes are shown for Mach 2.85.

Description: An account is given of experimental and computational results on the mixing of single, double, and opposed rows of jets characterized by an either isothermal or variable temperature mainstream in a confined subsonic crossflow; these flow configurations are typical of gas turbine combustor dilution chambers. It is established that the momentum-flux ratio is the most significant flow variable. Combinations of flow and geometry yielding optimum mixing were identified from experimental and computational results. The orifice spacing and momentum-flux relationships affected jet structure, which was significantly different between jets injected from the inner wall of a turn and those injected from the outer wall.

Description: The essential elements of a computational fluid dynamics analysis of the HL-20/personnel launch system aerothermal environment at hypersonic speeds including surface definition, grid generation, solution techniques, and visual representation of results are presented. Examples of solution technique validation through comparison with data from ground-based facilities are presented, along with results from computations at flight conditions. Computations at flight points indicate that real-gas effects have little or no effect on vehicle aerodynamics and, at these conditions, results from approximate techniques for determining surface heating are comparable with those obtained from Navier-Stokes solutions.

Description: Definition of the aerothermal environment is critical to any vehicle such as the HL-20 Personnel Launch System that operates within the hypersonic flight regime. Selection of an appropriate thermal protection system design is highly dependent on the accuracy of the heating-environment prediction. It is demonstrated that the entry environment determines the thermal protection system design for this vehicle. The methods used to predict the thermal environment for the HL-20 Personnel Launch System vehicle are described. Comparisons of the engineering solutions with computational fluid dynamic predictions, as well as wind-tunnel test results, show good agreement. The aeroheating predictions over several critical regions of the vehicle, including the stagnation areas of the nose and leading edges, windward centerline and wing surfaces, and leeward surfaces, are discussed. Results of predictions based on the engineering methods found within the MINIVER aerodynamic heating code are used in conjunction with the results of the extensive wind-tunnel tests on this configuration to define a flight thermal environment. Finally, the selection of the thermal protection system based on these predictions and current technology is described.

Description: Experiments were conducted to observe the cross-sectional structure and streamwise growth of round transverse liquid jets injected into a highly accelerated boundary layer in supersonic flow. The accompanying shock structure was also visualized. In one case, a round jet of acetone was injected into a fully turbulent Mach 2.5 boundary layer that was subsequently accelerated and partially laminarized through a sharp Prandtl-Meyer expansion corner. In the second case, a jet was injected into the laminarized Mach 3.2 boundary layer downstream of the expansion corner at the same jet-to-freestream momentum ratio. The jet and shock structure in both cases were visualized using schlieren optics. Wall-flow patterns were visualized using paints. It was found that the lateral spreading of jets injected downstream of the expansion fan was augmented close to the wall and had a cross-sectional structure significantly different from that of the jet injected upstream: the upstream jet spreads rapidly at the expansion corner in both the lateral and vertical directions.

Description: Computational results are presented to illustrate the powered aftbody effects of representing the scramjet inlet on a generic hypersonic vehicle with a fairing, to divert the external flow, as compared to an operating flow-through scramjet inlet. This study is pertinent to the ground testing of hypersonic, airbreathing models employing scramjet exhaust flow simulation in typical small-scale hypersonic wind tunnels. The comparison of aftbody effects due to inlet representation is well-suited for computational study, since small model size typically precludes the ability to ingest flow into the inlet and perform exhaust simulation at the same time. Two-dimensional analysis indicates that, although flowfield differences exist for the two types of inlet representations, little, if any, difference in surface aftbody characteristics is caused by fairing over the inlet.

Description: A detailed description of the bow-shock behavior associated with a conical-walled cavity with a flat circular base at M(infinity) = 10 is presented. An experimental test was performed on this configuration, and measurements of shock-oscillation frequency and amplitude, as well as shock shape, were recorded by a number of techniques. A laser-interferometer system was used for the first time during this test to determine bow-shock oscillation frequency in a nonintrusive manner. The primary behavior was a stable, periodically oscillating bow shock. Attention is also given to a violent bow-shock instability.

Description: Hypersonic reactive flowfield results are presented for the flow about an axisymmetric representation of a three-dimensional aeroassisted vehicle using the viscous shock-layer (VSL) and Navier-Stokes (NS) equations for low- and high-altitude entry conditions, respectively. Solutions are limited to the stagnation region. The NS results have been obtained by employing the surface recombination boundary condition for an 11-species air model, whereas the VSL calculations are carried out by using a 5-species model. Three widely used surface recombination rate models are implemented with these two flowfield calculation procedures. The low-density (or high altitude) NS solutions are obtained with surface slip equations and compare quite well with the direct simulation Monte Carlo predictions. For higher densities (or lower altitudes), the no-slip VSL solutions are considered adequate. Results from this study highlight the effect of surface catalyticity on surface heating and electron number density.

Description: Qualitative features of the inviscid instability characteristics of hypersonic boundary-layer flows over a flat plate are considered. The instability of a viscous hypersonic boundary layer which exists far downstream from the leading edge of the plate. It is shown that the vorticity mode of instability operates on a different lengthscale from that obtained using a Chapman viscosity law. The growth rate predicted by a linear viscosity law is found to overestimate the size of the growth rate. The inviscid instability of the boundary layer near the leading edge interaction zone is discussed focusing on the strong-interaction zone which occurs sufficiently close to the leading edge. The vorticity mode in this regime is found to be unstable.

Description: A closed-form analysis of flow in a two-dimensional subsonic wind tunnel that uses sidewall suction around the model to reduce sidewall boundary-layer effects is presented. The model problem that is treated involves a flat plate airfoil in a tunnel with a suction window shaped to permit an analytic solution. This solution shows that the lift coefficient depends explicitly on the porosity parameter of the suction window and implicitly on the suction pressure differential. For a given sidewall displacement thickness, the lift coefficient increases as the suction-window porosity decreases.

Description: In support of the NASA Balloon Program, small-scale balloons were flown with varying lifting gas and total system mass. Instrument packages were developed to measure and record acceleration and temperature data during these tests. Top fitting and instrument payload accelerations were measured from launch to steady state ascent and through ballast drop transients. The development of the small lightweight self-powered Stowaway Special instrument packages is discussed along with mathematical models developed to determine gas mass, drag and virtual mass coefficients.

Description: The flow field over an airfoil contains several distinct flow characteristics, e.g. laminar, transitional, turbulent boundary layer flow, flow separation, unstable free shear layers, and a wake. This diversity of flow regimes taxes the presently available Reynolds averaged turbulence models. Such models are generally tuned to predict a particular flow regime, and adjustments are necessary for the prediction of a different flow regime. Similar difficulties are likely to emerge when the large eddy simulation technique is applied with the widely used Smagorinsky model. This model has not been successful in correctly representing different turbulent flow fields with a single universal constant and has an incorrect near-wall behavior. Germano et al. (1991) and Ghosal, Lund & Moin have developed a new subgrid-scale model, the dynamic model, which is very promising in alleviating many of the persistent inadequacies of the Smagorinsky model: the model coefficient is computed dynamically as the calculation progresses rather than input a priori. The model has been remarkably successful in prediction of several turbulent and transitional flows. We plan to simulate turbulent flow over a '2D' airfoil using the large eddy simulation technique. Our primary objective is to assess the performance of the newly developed dynamic subgrid-scale model for computation of complex flows about aircraft components and to compare the results with those obtained using the Reynolds average approach and experiments. The present computation represents the first application of large eddy simulation to a flow of aeronautical interest and a key demonstration of the capabilities of the large eddy simulation technique.

Description: The ability to model turbulence near solid walls and other types of boundaries is important in predicting complex engineering flows. Most turbulence modeling has concentrated either on flows which are nearly homogeneous or isotropic, or on turbulent boundary layers. Boundary layer models usually rely very heavily on the presence of mean shear and the production of turbulence due to that mean shear. Most other turbulence models are based on the assumption of quasi-homogeneity. However, there are many situations of engineering interest which do not involve large shear rates and which are not quasi-homogeneous or isotropic. Shear-free turbulent boundary layers are the prototypical example of such flows, with practical situations being separation and reattachment, bluff body flow, high free-stream turbulence, and free surface flows. Although these situations are not as common as the variants of the flat plate turbulent boundary layer, they tend to be critical factors in complex engineering situations. The models developed are intended to extend classical quasi-homogeneous models into regions of large inhomogeneity. These models do not rely on the presence of mean shear or production, but are still applicable when those additional effects are included. Although the focus is on shear-free boundary layers as tests for these models, results for standard shearing boundary layers are also shown.

Description: The problem to be addressed in this work is that of modeling nondimensional force and moment aerodynamic coefficients over the entire subsonic envelope. The particular application discussed here is the Z force coefficient for the F-18 High Angle of Attack Research Vehicle (HARV).

Description: This project is a continuation of the work performed in the summers of 1991 and 1992, during which a 9-block structured grid for the computational domain around each of the four NTF 65-degree Delta Wing models with the sting mount were created. The objective of the project is to validate and supplement the test data on the wing models by computing the viscous flow field about the models.

Description: Wingtip vortices are generated by lifting airfoils; their salient features are compactness and relatively slow rate of decay. The principal motivation for studying the far field evolution of wingtip vortices is the need to understand and predict the extent of the vortex influence during aircraft take-off or landing. On submarines a wingtip vortex ingested into a propeller can be a source of undesirable noise. The main objectives of this research are (1) to establish theoretical understanding of the principal mechanisms that govern the later (diffusive) stages of a turbulent vortex, (2) to develop a turbulence closure model representing the basic physical mechanisms that control the vortex diffusive stage, and further (3) to investigate coupling between the near and far field evolutions; in other words, to study the effect of initial conditions on the vortex lifetime and the ultimate state.

Description: A mathematical model for a two-phase flow laminar airfoil in simulated heavy rain has been established. The set of non-linear partial differential equations has been converted into a set of finite difference equations; appropriate initial and boundary conditions are provided. The numerical results are compared with the experimental measurements. They show good agreement in quality.

Description: Several multigrid schemes are considered for the numerical computation of viscous hypersonic flows. For each scheme, the basic solution algorithm employs upwind spatial discretization with explicit multistage time stepping. Two-level versions of the various multigrid algorithms are applied to the two-dimensional advection equation, and Fourier analysis is used to determine their damping properties. The capabilities of the multigrid methods are assessed by solving two different hypersonic flow problems. Some new multigrid schemes, based on semicoarsening strategies, are shown to be quite effective in relieving the stiffness caused by the high-aspect-ratio cells required to resolve high Reynolds number flows. These schemes exhibit good convergence rates for Reynolds numbers up to 200 x 10(exp 6).

Description: The objective of this educational exercise is to have students build and evaluate simple wing structures, and in doing so, learn about materials choices and lightweight construction methods. A list of equipment and supplies and the procedure for the experiment are presented.

Description: Boundary-layer transition and relaminarization may have a critical effect on the flow development about multi-element high-lift systems of subsonic transport jets with swept wings. The purpose of the research is to study these transition phenomena in the leading-edge region of the various elements of a high-lift system. The flow phenomena studied include transition to the attachment-line flow, relaminarization, and crossflow instability, and transition. The calculations are based on pressure distributions measured in flight on the NASA Transport Systems Research Vehicle (Boeing 737-100) at a wing station where the flow approximated infinite swept wing conditions. The results indicate that significant regions of laminar flow can exist on all flap elements in flight. In future flight experiments (planned for January-February, 1994) the extent of these regions, the transition mechanisms and the effect of laminar flow on the high-lift characteristics of the multi-element system will be further explored.

Description: Aerodynamic heating in hypersonic space vehicles is an important factor to be considered in their design. Therefore the designers of such vehicles need reliable heat transfer data in this respect for a successful design. Such data is usually produced by testing the models of hypersonic surfaces in wind tunnels. Most of the hypersonic test facilities at present are conventional blow-down tunnels whose run times are of the order of several seconds. The surface temperatures on such models are obtained using standard techniques such as thin-film resistance gages, thin-skin transient calorimeter gages and coaxial thermocouple or video acquisition systems such as phosphor thermography and infrared thermography. The data are usually reduced assuming that the model behaves like a semi-infinite solid (SIS) with constant properties and that heat transfer is by one-dimensional conduction only. This simplifying assumption may be valid in cases where models are thick, run-times short, and thermal diffusivities small. In many instances, however, when these conditions are not met, the assumption may lead to significant errors in the heat transfer results. The purpose of the present paper is to investigate this aspect. Specifically, the objectives are as follows: (1) to determine the limiting conditions under which a model can be considered a semi-infinite body; (2) to estimate the extent of errors involved in the reduction of the data if the models violate the assumption; and (3) to come up with correlation factors which when multiplied by the results obtained under the SIS assumption will provide the results under the actual conditions.

Description: Existing aerodynamic design methods have generally concentrated on the optimization of airfoil or wing shapes to produce a minimum drag while satisfying some basic constraints such as lift, pitching moment, or thickness. Since the minimization of drag almost always precludes the existence of separated flow, the evaluation and validation of these design methods for their robustness and accuracy when separated flow is present has not been aggressively pursued. However, two new applications for these design tools may be expected to include separated flow and the issues of aerodynamic design with this feature must be addressed. The first application of the aerodynamic design tools is the design of airfoils or wings to provide an optimal performance over a wide range of flight conditions (multipoint design). While the definition of 'optimal performance' in the multipoint setting is currently being hashed out, it is recognized that given a wide range of flight conditions, it will not be possible to ensure a minimum drag constraint at all conditions, and in fact some amount of separated flow (presumably small) may have to be allowed at the more demanding flight conditions. Thus a multipoint design method must be tolerant of the existence of separated flow and may include some controls upon its extent. The second application is in the design of wings with extended high speed buffet boundaries of their flight envelopes. Buffet occurs on a wing when regions of flow separation have grown to the extent that their time varying pressures induce possible destructive effects upon the wing structure or adversely effect either the aircraft controllability or passenger comfort. A conservative approach to the expansion of the buffet flight boundary is to simply expand the flight envelope of nonseparated flow under the assumption that buffet will also thus be alleviated. However, having the ability to design a wing with separated flow and thus to control the location, extent and severity of the separated flow regions may allow aircraft manufacturers to gain an advantage in the early design stages of an aircraft, when configuration changes are relatively inexpensive to make. The goal of the summer research at NASA Langley Research Center (LaRC) was twofold: first, to investigate a particular airfoil design problem observed under conditions of strong shock induced flow separation on the upper surface of an airfoil at transonic conditions; and second, to suggest and investigate design methodologies for the prediction (or detection) and control of flow separation. The context of both investigations was to use an existing two dimensional Navier-Stokes flow solver and the constrained direct/iterative surface curvature (CDISC) design algorithm developed at LaRC. As a lead in to the primary task, it was necessary to gain a familiarity with both the design method and the computational analysis and to perform the FORTRAN coding needed to couple them together.

Description: Separation and reattachment of turbulent shear layers is observed in many important engineering applications, yet it is poorly understood. This has motivated many studies on understanding and predicting the processes of separation and reattachment of turbulent shear layers. Both of the situations in which separation is induced by adverse pressure gradient, or by discontinuities of geometry, have attracted attention of turbulence model developers. Formulation of turbulence closure models to describe the essential features of separated turbulent flows accurately is still a formidable task. Computations of separated flows associated with sharp-edged bluff bodies are described. For the past two decades, the backward-facing step flow, the simplest separated flow, has been a popular test case for turbulence models. Detailed studies on the performance of many turbulence models, including two equation turbulence models and Reynolds stress models, for flows over steps can be found in the papers by Thangam & Speziale and Lasher & Taulbee). These studies indicate that almost all the existing turbulence models fail to accurately predict many important features of back step flow such as reattachment length, recovery rate of the redeveloping boundary layers downstream of the reattachment point, streamlines near the reattachment point, and the skin friction coefficient. The main objectives are to calculate flows over backward and forward-facing steps using the NRSM and to make use of the newest DNS data for detailed comparison. This will give insights for possible improvements of the turbulence model.

Description: A computational study has been conducted to evaluate the performance of various turbulence models. The NASA P8 inlet, which represents cruise condition of a typical hypersonic air-breathing vehicle, was selected as a test case for the study; the PARC2D code, which solves the full two dimensional Reynolds-averaged Navier-Stokes equations, was used. Results are presented for a total of six versions of zero- and two-equation turbulence models. Zero-equation models tested are the Baldwin-Lomax model, the Thomas model, and a combination of the two. Two-equation models tested are low-Reynolds number models (the Chien model and the Speziale model) and a high-Reynolds number model (the Launder and Spalding model).

Description: The effectiveness of two types of hypersonic decelerators is examined: mechanically deployable flares and inflatable ballutes. Computational fluid dynamics (CFD) is used to predict the flowfield around a solid rocket motor (SRM) with a deployed decelerator. The computations are performed with an ideal gas solver using an effective specific heat ratio of 1.15. The results from the ideal gas solver are compared to computational results from a thermochemical nonequilibrium solver. The surface pressure coefficient, the drag, and the extend of the compression corner separation zone predicted by the ideal gas solver compare well with those predicted by the nonequilibrium solver. The ideal gas solver is computationally inexpensive and is shown to be well suited for preliminary design studies. The computed solutions are used to determine the size and shape of the decelerator that are required to achieve a drag coefficient of 5. Heat transfer rates to the SRM and the decelerators are predicted to estimate the amount of thermal protection required.

Description: The motivation for this work is the fact that in turbulent flows where compressibility effects are important, they are often poorly understood. A few examples of such flows are those associated with astrophysical phenomena and those found in combustion chambers, supersonic diffusers and nozzles, and over high-speed airfoils. For this project, we are primarily interested in compressibility effects near solid surfaces. Our main objective is an improved understanding of the fundamentals of compressible wall-bounded turbulence, which can in turn be used to cast light upon modeling concepts such as the Morkovin hypothesis and the Van Driest transformation. To this end, we have performed a direct numerical simulation (DNS) study of supersonic turbulent flow in a plane channel with constant-temperature walls. All of the relevant spatial and temporal scales are resolved so that no sub grid scale or turbulence model is necessary. The channel geometry was chosen so that finite Mach number effects can be isolated by comparing the present results to well established incompressible channel data. Here the fluid is assumed to be an ideal gas with constant specific heats, constant Prandtl number, and power-law temperature-dependent viscosity. Isothermal-wall boundary conditions are imposed so that a statistically stationary state may be obtained. The flow is driven by a uniform (in space) body force (rather than a mean pressure gradient) to preserve stream wise homogeneity, with the body force defined so that the total mass flux is constant.

Description: Recently, dense gases have been investigated for many engineering applications such as for turbomachinery and wind tunnels. Supersonic nozzle design for these gases is complicated by their nonclassical behavior in the transonic flow regime. In this paper a method of characteristics (MOC) is developed for two-dimensional (planar) and, primarily, axisymmetric flow of a van der Waals gas. Using a straight aortic line assumption, a centered expansion is used to generate an inviscid wall contour of minimum length. The van der Waals results are compared to previous perfect gas results to show the real gas effects on the flow properties and inviscid wall contours.

Description: A study of Navier-Stokes calculations of flows about multielement airfoils using a chimera grid approach is presented. The chimera approach utilizes structured, overlapped grids which allow great flexibility of grid arrangement and simplifies grid generation. Calculations are made for two-, three-, and four-element airfoils, and modeling of the effect of gap distance between elements is demonstrated for a two element case. Solutions are obtained using the thin-layer form of the Reynolds averaged Navier-Stokes equations with turbulence closure provided by the Baldwin-Lomax algebraic model or the Baldwin-Barth one equation model. The Baldwin-Barth turbulence model is shown to provide better agreement with experimental data and to dramatically improve convergence rates for some cases. Recently developed, improved farfield boundary conditions are incorporated into the solver for greater efficiency. Computed results show good comparison with experimental data which include aerodynamic forces, surface pressures, and boundary layer velocity profiles.

Description: A detailed comparison of numerical results obtained by solving the Euler equations for the inviscid flow over delta wings in reentry configurations is presented. The investigation involves a side-by-side comparison between independently developed upwind Euler solvers at VKI (M3D) and NASA Johnson Space Center (E3D) using identical grids. In both solvers the governing equations are integrated by means of time marching finite volume shock capturing methods, based on a cell centered upwind evaluation of the cell face fluxes and nonlinear limiters. High resolution schemes are obtained via MUSCL characteristic variable extrapolation to ensure total variation diminishing (TVD) properties and therefore monotonic discontinuity capturing. Further comparisons are made with several published results by other authors based on both upwind and central discretizations. The present results contribute in asserting the high resolution upwind TVD schemes as the most reliable numerical technique to handle the strong discontinuities typical of high speed flows. The performances of the two upwind solvers are satisfactory and the numerical results in good agreement. However, the important issue of reaching a grid converged solution for the present complex three-dimensional inviscid flow problems was not achieved.

Description: Several analytical aerodynamic design tools that were applied to the Pegasus air-launched space booster were evaluated using flight measurements. The study was limited to existing codes and was conducted with limited computational resources. The flight instrumentation was constrained to have minimal impact on the primary Pegasus missions. Where appropriate, the flight measurements were compared with computational data. Aerodynamic performance and trim data from the first two flights were correlated with predictions. Local measurements in the wing and wing-body interference region were correlated with analytical data. This complex flow region includes the effect of aerothermal heating magnification caused by the presence of a corner vortex and interaction of the wing leading edge shock and fuselage boundary layer. The operation of the first two missions indicates that the aerodynamic design approach for Pegasus was adequate, and data show that acceptable margins were available. Additionally, the correlations provide insight into the capabilities of these analytical tools for more complex vehicles in which design margins may be more stringent.

Description: Critical issues concerning the modeling of low density hypervelocity flows where thermochemical nonequilibrium effects are pronounced are discussed. Emphasis is on the development of validated analysis tools, and the activity in the NASA Ames Research Center's Aerothermodynamics Branch is described. Inherent in the process is a strong synergism between ground test and real gas computational fluid dynamics (CFD). Approaches to develop and/or enhance phenomenological models and incorporate them into computational flowfield simulation codes are discussed. These models were partially validated with experimental data for flows where the gas temperature is raised (compressive flows). Expanding flows, where temperatures drop, however, exhibit somewhat different behavior. Experimental data for these expanding flow conditions is sparse and reliance must be made on intuition and guidance from computational chemistry to model transport processes under these conditions. Ground based experimental studies used to provide necessary data for model development and validation are described. Included are the performance characteristics of high enthalpy flow facilities, such as shock tubes and ballistic ranges.

Description: A review of measurement techniques used to obtain aerodynamic forces and moments and surface/flow field pressures for models tested in impulse hypersonic-hypervelocity facilities and in conventional-type hypersonic wind tunnels is presented. Although force and moment measurement techniques presently used in hypersonic wind tunnels are relatively unchanged from the 1960's and 1970's, significant advances have recently been made for impulse facilities. For both hypersonic wind tunnels and impulse facilities, the state-of-the-art has advanced via refinements, improved test techniques, and advances in semiconductor technology, data acquisition systems, and computers. The introduction of electronically scanned pressure systems over a decade ago 'revolutionized' pressure measurements in hypersonic wind tunnels and a second 'revolution' is impending with the development and application of optical, two-dimensional, global pressure measurement techniques. The development and continued refinement of miniature piezoresistive transducers has provided the capability to perform detailed surface pressure measurements on relatively small, complex models in impulse facilities; these transducers also provided the capability for intrusive flow field pressure measurements with miniature survey rakes.

Description: A summary of three-dimensional simulations on the hypersonic rarefied flows in an effort to understand the highly nonequilibrium flows about space vehicles entering the Earth's atmosphere for a realistic estimation of the aerothermal loads is presented. Calculations are performed using the direct simulation Monte Carlo method with a five-species reacting gas model, which accounts for rotational and vibrational internal energies. Results are obtained for the external flows about various bodies in the transitional flow regime. For the cases considered, convective heating, flowfield structure and overall aerodynamic coefficients are presented and comparisons are made with the available experimental data. The agreement between the calculated and measured results are very good.

Description: An overview of the NATO Advanced Research Workshop is presented. It describes the process followed to obtain a group consensus on the main technical recommendations for each of the five technical sessions of the Workshop and presents the general conclusions and recommendations for future research agreed upon by the workshop participants.

Description: An experimental program has been conducted to assess performance of a transport multielement airfoil at flight Reynolds numbers. The studies were performed at chord Reynolds numbers as high as 16 million in the NASA Langley Low Turbulence Pressure Tunnel. Sidewall boundary-layer control to enforce flow two dimensionality was provided via an endplate suction system. The basic airfoil was an 11.55 percent thick supercritical airfoil representative of the stall critical station of a new-generation transport aircraft wing. The multielement airfoil was configured as a three-element airfoil with slat and flap chord ratios of 14.48 percent and 30 percent respectively. Testing focused on the development of landing configurations with high maximum lift capability and the assessment of Reynolds and Mach number effects. Also assessed were high-lift performance effects due to devices such as drooped spoilers and trailing-edge wedges. The present experimental studies revealed significant effects on high-lift airfoil performance due to Reynolds and Mach number variations and favorable lift increments at approach angles of attack due to the use of drooped spoilers or trailing-edge wedges. However, no substantial improvements in maximum lift capability were identified. A recently developed high performance single-segment flap was also tested and results indicated considerable improvements in lift and drag performance over existing airfoils. Additionally, it was found that this new flap shape at its optimum rigging was less sensitive to Reynolds number variations than previous designs.

Description: The incompressible, viscous, turbulent flow over single and multi-element airfoils is numerically simulated in an efficient manner by solving the incompressible Navier-Stokes equations. The solution algorithm employs the method of pseudo compressibility and utilizes an upwind differencing scheme for the convective fluxes, and an implicit line-relaxation scheme. The motivation for this work includes interest in studying high-lift take-off and landing configurations of various aircraft. In particular, accurate computation of lift and drag at various angles of attack up to stall is desired. Two different turbulence models are tested in computing the flow over an NACA 4412 airfoil; an accurate prediction of stall is obtained. The approach used for multi-element airfoils involves the use of multiple zones of structured grids fitted to each element. Two different approaches are compared; a patched system of grids, and an overlaid Chimera system of grids. Computational results are presented for two-element, three-element, and four-element airfoil configurations. Excellent agreement with experimental surface pressure coefficients is seen. The code converges in less than 200 iterations, requiring on the order of one minute of CPU time on a CRAY YMP per element in the airfoil configuration.

Description: The flow over multi-element airfoils with flat-plate lift-enhancing tabs was numerically investigated. Tabs ranging in height from 0.25 percent to 1.25 percent of the reference airfoil chord were studied near the trailing edge of the main-element. This two-dimensional numerical simulation employed an incompressible Navier-Stokes solver on a structured, embedded grid topology. New grid refinements were used to improve the accuracy of the solution near the overlapping grid boundaries. The effects of various tabs were studied at a constant Reynolds number on a two-element airfoil with a slotted flap. Both computed and measured results indicated that a tab in the main-element cove improved the maximum lift and lift-to-drag ratio relative to the baseline airfoil without a tab. Computed streamlines revealed that the additional turning caused by the tab may reduce the amount of separated flow on the flap. A three-element airfoil was also studied over a range of Reynolds numbers. For the optimized flap rigging, the computed and measured Reynolds number effects were similar. When the flap was moved from the optimum position, numerical results indicated that a tab may help to reoptimize the airfoil to within 1 percent of the optimum flap case.

Description: Airloads measured on a two-bladed helicopter rotor in flight during the Ames' Tip Aerodynamic and Acoustic Test are compared with calculations from a comprehensive helicopter analysis (CAMRAD/JA), and the pressures compared with calculations from a full-potential rotor code (FPR). The flight-test results cover an advance ratio range of 0.19 to 0.38. The lowest-speed case is characterized by the presence of significant blade-vortex interactions. Good correlation of peak-to-peak vortex-induced loads and the corresponding pressures is obtained. Results of the correlation for this two-bladed rotor are substantially similar to those for three- and four-bladed rotors, including the tip-vortex core size for best correlation, calculation of the peak-to-peak loads on the retreating side, and calculation of vortex iduced loads on inboard radial stations. The higher-speed cases are characterized by the presence of transonic flow on the outboard sections of the blade. Comparison of calculated and measured airloads on the advancing side is not considered appropriate because the presence of shocks makes chordwise integration of the measured data difficult. However, good correlation of the corresponding pressures is obtained.

Description: Experimental and computational tests were performed on a VR-12 airfoil to determine if the dynamic-stall behavior that normally accompanies high-angle pitch oscillations could be modified by segmenting the forward portion of the airfoil and extending it ahead of the main element. In the extended position the configuration would appear as an airfoil with a leading-edge slat, and in the retracted position it would appear as a conventional VR-12 airfoil. The calculations were obtained from a numerical code that models the vorticity transport equation for an incompressible fluid. These results were compared with test data from the water tunnel facility of the Aeroflightdynamics Directorate at Ames Research Center. Steady and unsteady flows around both airfoils were examined at angles of attack between 0 and 30 deg. The Reynolds number was fixed at 200,000 and the unsteady pitch oscillations followed a sinusoidal motion described by alpha = alpha(sub m) + 10 deg sin(omega t). The mean angle (alpha(sub m)) was varied from 10 to 20 deg and the reduced frequency from 0.05 to 0.20. The results from the experiment and the calculations show that the extended-slat VR-12 airfoil experiences a delay in both static and dynamic stall not experienced by the basic VR-12 airfoil.

Description: The quality of reentry simulation for Shuttle, HERMES, Sanger, and NASP systematically suffers from the strong non-equilibrium of rotational and vibrational temperature due to the rapid acceleration of the test gas in the nozzle. Therefore the determination of temperatures is necessary and, if possible, preferable by a non-intrusive technique. The specific interests of this review are optical techniques such as electron beam fluorescence, laser-induced fluorescence, and coherent anti-Stokes Raman scattering. The capabilities available for local measurements with temporal resolution and quantitative accuracy are discussed for velocity, temperature, density, species concentrations, and fluctuations due to turbulence. The applicability of these methods of measurement is presented and discussed for the coming topic in aerothermodynamics: experimental techniques of hot gases in high enthalpy flows.

Description: Design studies are underway for a variety of hypersonic flight vehicles. The National Aero-Space Plane will provide a reusable, single-stage-to-orbit capability for routine access to low earth orbit. Flight-capable satellites will dip into the atmosphere to maneuver to new orbits, while planetary probes will decelerate at their destination by atmospheric aerobraking. To supplement limited experimental capabilities in the hypersonic regime, computational fluid dynamics (CFD) is being used to analyze the flow about these configurations. The governing equations include fluid dynamic as well as chemical species equations, which are being solved with new, robust numerical algorithms. Examples of CFD applications to hypersonic vehicles suggest an important role this technology will play in the development of future aerospace systems. The computational resources needed to obtain solutions are large, but solution adaptive grids, convergence acceleration, and parallel processing may make run times manageable.

Description: The primary focus of this research centered on improving the vorticity-capturing properties of Eulerian flow solvers in convection-dominated flows. The work was motivated by the fact that excessive numerical diffusion observed in current Euler solvers on grids of feasible size is a major obstacle to accurately computing the flow about helicopter rotors, especially in the forward flight case. The work was originally supported by NASA Grant NAG2-421 for two years. The research then continued for four more years under an ONR contract. The research goal of reduction of numerical diffusion in Euler solvers remained the same, although the target application was now the prediction of internal vortical secondary flows. This was part of a larger ONR-sponsored research program at MIT aimed at reducing the noise of turbopumps. The research work culminated in the PhD thesis titled 'A Coupled Eulerian/Lagrangian Method for the Solution of Three-Dimensional Vortical Flows'.

Description: Hypersonic lateral and directional stability characteristics measured on a 60 deg half-angle elliptical cone, which was raked at an angle of 73 deg from the cone centerline and with an ellipsoid nose (ellipticity equal to 2.0 in the symmetry plane), are presented for angles of attack from -10 to 10 deg. The high normal-shock density ratio of a real gas was simulated by tests at a Mach number of 6 in air and CF4 (density ratio equal to 5.25 and 12.0, respectively). Tests were conducted in air at Mach 6 and 10 and in CF4 at Mach 6 to examine the effects of Mach number, Reynolds number, and normal-shock density ratio. Changes in Mach number from 6 to 10 in air or in Reynolds number by a factor of 4 at Mach 6 had a negligible effect on lateral and directional stability characteristics. Variations in normal-shock density ratio had a measurable effect on lateral and directional aerodynamic coefficients, but no significant effect on lateral and directional stability characteristics. Tests in air and CF4 indicated that the configuration was laterally and directionally stable through the test range of angle of attack.

Description: An experimental program was conducted in the Langley 16-Foot Transonic Tunnel to evaluate the performance effects of various types of pylons on a 1/17th-scale, low-wing transport model. The model wing was designed for cruise at a Mach number of 0.77 and a lift coefficient of 0.55. The pylons were tested at two wing semispan locations over a range of toe-in angles. The effects of toe-in angle were found to be minimal, but the variation in geometry had a more pronounced effect on the lift characteristics of the model. A pylon whose maximum thickness occurred at the wing trailing edge, known as a compression pylon, proved to be the best choice in terms of retaining the flow characteristics of the wing without pylons. Practical considerations such as structural viability may necessitate modification of the compression pylon concept in order to take advantage of its apparent benefits.

Description: The unsteady aerodynamic effects due to the separated flow around slender delta wings in motion were analyzed. By combining the unsteady flow field solution with the rigid body Euler equations of motion, self-induced wing rock motion is simulated. The aerodynamic model successfully captures the qualitative characteristics of wing rock observed in experiments. For the one degree of freedom in roll case, the model is used to look into the mechanisms of wing rock and to investigate the effects of various parameters, like angle of attack, yaw angle, displacement of the separation point, and wing inertia. To investigate the roll and yaw coupling for the delta wing, an additional degree of freedom is added. However, no limit cycle was observed in the two degree of freedom case. Nonetheless, the model can be used to apply various control laws to actively control wing rock using, for example, the displacement of the leading edge vortex separation point by inboard span wise blowing.

Description: Implicit and explicit spatial differencing techniques with fourth order accuracy have been developed. The implicit technique is based on the Pade compact scheme. A Dispersion Relation Preserving concept has been incorporated into both of the numerical schemes. Two dimensional Euler computation of a spatially-developing free shear flow, with and without external excitation, has been performed to demonstrate the capability of numerical schemes developed. Results are in good agreement with theory and experimental observation regarding the growth rate of fluctuating velocity, the convective velocity, and the vortex-pairing process.

Description: A new method of generating unstructured triangular/tetrahedral grids with high-aspect-ratio cells is proposed. The method is based on new grid-marching strategy referred to as 'advancing-layers' for construction of highly stretched cells in the boundary layer and the conventional advancing-front technique for generation of regular, equilateral cells in the inviscid-flow region. Unlike the existing semi-structured viscous grid generation techniques, the new procedure relies on a totally unstructured advancing-front grid strategy resulting in a substantially enhanced grid flexibility and efficiency. The method is conceptually simple but powerful, capable of producing high quality viscous grids for complex configurations with ease. A number of two-dimensional, triangular grids are presented to demonstrate the methodology. The basic elements of the method, however, have been primarily designed with three-dimensional problems in mind, making it extendible for tetrahedral, viscous grid generation.

Description: The receptivity theory of Goldstein and Ruban is extended within the nonasymptotic (quasi-parallel) framework of Zavol'skii et al to predict the roughness-induced generation of stationary and nonstationary instability waves in three-dimensional, incompressible boundary layers. The influence of acoustic-wave orientation, as well as that of different types of roughness geometries, including isolated roughness elements, periodic arrays, and two-dimensional lattices of compact roughness shapes, as well as random, but spatially homogeneous roughness distributions, is examined. The parametric study for the Falkner-Skan-Cooke family of boundary layers supports our earlier conjecture that the initial amplitudes of roughness-induced stationary vortices are likely to be significantly larger than the amplitudes of similarly induced nonstationary vortices in the presence of acoustic disturbances in the free stream. Maximum unsteady receptivity occurs when the acoustic velocity fluctuation is aligned with the wavenumber vector of the unsteady vortex mode. On the other hand, roughness arrays that are oriented somewhere close to the group velocity direction are likely to produce higher instability amplitudes. Limitations of the nonasymptotic theory are discussed, and future work is suggested.

Description: Panel flutter is a form of dynamic aeroelastic instability resulting from the interaction between motion of an aircraft structural panel and the aerodynamic loads exerted on that panel by air flowing past one of the faces. It differs from lifting surface flutter in the sense that it is not usually catastrophic, the panel's motion being limited by nonlinear membrane stresses produced by the transverse displacement. Above some critical airflow condition, the linear instability grows to a limit cycle . The present investigation studies panel flutter in an aerodynamic regime known as 'free molecule flow', wherein intermolecular collisions can be neglected and loads are caused by interactions between individual molecules and the bounding surface. After collision with the panel, molecules may be reflected specularly or reemitted in diffuse fashion. Two parameters characterize this process: the 'momentum accommodation coefficient', which is the fraction of the specularly reflected molecules; and the ratio between the panel temperature and that of the free airstream. This model is relevant to the case of hypersonic flight vehicles traveling at very high altitudes and especially for panels oriented parallel to the airstream or in the vehicle's lee. Under these conditions the aerodynamic shear stress turns out to be considerably larger than the surface pressures, and shear effects must be included in the model. This is accomplished by means of distributed longitudinal and bending loads. The former can cause the panel to buckle. In the example of a simply-supported panel, it turns out that the second mode of free vibration tends to dominate the flutter solution, which is carried out by a Galerkin analysis. Several parametric studies are presented. They include the effects of (1) temperature ratio; (2) momentum accommodation coefficient; (3) spring parameters, which are associated with how the panel is connected to adjacent structures; (4) a parameter which relates compressive end load to its value which would cause classical column buckling; (5) a parameter proportional to the pressure differential between the front and back faces; and (6) initial curvature. The research is completed by an investigation into the possibility of accounting for molecular collisions, which proves to be infeasible given the speeds of current mainframe supercomputers.

Description: An investigation was conducted in the Langley 14- by 22-Foot Subsonic Tunnel to establish a transition data base for an unmanned aerial vehicle utilizing a powered-lift ejector system and to evaluate alterations to the ejector system for improved vehicle performance. The model used in this investigation was a 20-percent-scale, blended-body, arrow-wing configuration with integrated twin rectangular ejectors. The test was conducted from hover through transition conditions with variations in angle of attack, angle of sideslip, free-stream dynamic pressure, nozzle pressure ratio, and model ground height. Force and moment data along with extensive surface pressure data were obtained. A laser velocimeter technique for measuring inlet flow velocities was demonstrated at a single flow condition, and also a low order panel method was successfully used to numerically simulate the ejector inlet flow.

Description: A large-scale, outdoor, ground-based test capability for acquiring aerodynamic data in a simulated rain environment was developed at the Langley Aircraft Landing Dynamics Facility (ALDF) to assess the effect of heavy rain on airfoil performance. The ALDF test carriage was modified to transport a 10-ft-chord NACA 64210 wing section along a 3000-ft track at full-scale aircraft approach speeds. An overhead rain simulation system was constructed along a 525-ft section of the track with the capability of producing simulated rain fields of 2, 10, 30, and 40 in/hr. The facility modifications, the aerodynamic testing and rain simulation capability, the design and calibration of the rain simulation system, and the operational procedures developed to minimize the effect of wind on the simulated rain field and aerodynamic data are described in detail. The data acquisition and reduction processes are also presented along with sample force data illustrating the environmental effects on data accuracy and repeatability for the 'rain-off' test condition.

Description: An experimental study was performed at supersonic speeds to measure wing and body spanwise pressure distributions on an axisymmetric-body delta wing model on which the wing vertical location on the body was systematically varied from low- to high-mounted positions. In addition, for two of these positions both horizontal and radial wing angular orientations relative to the body were tested, and roll angle effects were investigated for one of the positions. Seven different wing-body configurations and a body-alone configuration were studied. The test was conducted at Mach numbers from 1.70 to 2.86 at angles of attack from about -4 deg to 24 deg. Pressure orifices were located at three longitudinal stations on each wing-body model, and at each station the orifices were located completely around the body, along the lower surface of the right wing (looking upstream), and along the upper surface of the left wing. All pressure coefficient data are tabulated and selected samples are shown graphically to illustrate the effects of the test variables. The effects of angle of attack, roll angle, Mach number, longitudinal station, wing vertical location, wing angular orientation, and wing-body juncture are analyzed. The vertical location of the wing on the body had a very strong effect on the body pressures. For a given angle of attack at a roll angle of 0 deg, the pressures were virtually constant in the spanwise direction across the windward surfaces of the wing-body combination. Pressure-relieving, channeling, and vortex effects were noted in the data.

Description: A ground vortex, produced when a jet impinges on the ground in the presence of cross flow, is encountered by V/STOL aircraft hovering near the ground and is known to be hazardous to the aircraft. The objective of this research was to identify a ground-based technique by which both the mean size and fluctuation in size of the ground vortex could be reduced. A simple passive method was identified and examined in the laboratory. Specifically, one or two fine wire mesh screens (ground fences) bent in a horseshoe shape and located on the ground in front of the jet impingement point proved to be very effective. The ground fences work by decreasing the momentum of the upstream-traveling wall jet, effectively causing a higher freestream-to-jet velocity ratio (V(sub infinity)/V(sub j)) and thus, a ground vortex smaller in size and unsteadiness. At(V(sub infinity)/V(sub j)) = 0.15, the addition of a single ground fence resulted in a 70 percent reduction in mean size of the ground vortex. With two ground fences, the mean size decreased by about 85 percent. Fluctuations in size decreased nearly in proportion to the mean size, for both the single and double fence configurations. These results were consistent over a wide range of jet Reynolds number (10(exp 4) less than Re(sub jet) less than 10(exp 5)); further development and full-scale Reynolds number testing are required, however, to determine if this technique can be made practical for the case of actual VTOL aircraft.

Description: The summation-by-parts energy norm is used to establish a new class of high-order finite-difference techniques referred to here as 'cyclo-difference' techniques. These techniques are constructed cyclically from stable subelements, and require no numerical boundary conditions; when coupled with the simultaneous approximation term (SAT) boundary treatment, they are time asymptotically stable for an arbitrary hyperbolic system. These techniques are similar to spectral element techniques and are ideally suited for parallel implementation, but do not require special collocation points or orthogonal basis functions. The principal focus is on methods of sixth-order formal accuracy or less; however, these methods could be extended in principle to any arbitrary order of accuracy.

Description: The purpose of this paper is to present the Structural Dynamics Division's research accomplishments for F.Y. 1992 and research plans for F.Y. 1993. The work under each Branch (technical area) is described in terms of highlights of accomplishments during the past year and highlights of plans for the current year as they relate to 5-year plans for each technical area. This information will be useful in program coordination with other government organizations and industry in areas of mutual interest.