ALBERT

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 SUBAERF2 program was developed to provide for the aerodynamic analysis and design of low speed wing flap systems. SUBAERF2 is based on a linearized theory lifting surface solution. It is particularly well suited to configurations which, because of high speed flight requirements, must employ thin wings with highly swept leading edges. The program is applicable to wings with either sharp or rounded leading edges. This program is a new and improved version of LAR-13116 and LAR-12987, which it replaces. The low speed aerodynamic analysis method used in SUBAERF2 provides estimates of wing performance which include the effects of attainable leading-edge thrust and vortex lift. This basic aerodynamic analysis method has been improved to provide for the convenient, efficient and accurate treatment of simple leading-edge and trailing-edge flap systems. The user inputs flap geometry directly. Solutions can be found for various combinations of leading and trailing edge flap deflections. The program provides for the simultaneous analysis of up to 25 pairs of leading-edge and trailing-edge flap deflection schedules. A revised attainable thrust algorithm improves accuracy at the low Mach numbers sometimes encountered in wind tunnel testing. Also added is a means of estimating the distribution of leading edge separation vortex forces. The revised program has been particularly useful in the subsonic analysis of vehicles designed for supersonic cruise. The SUBAERF2 program is written in FORTRAN V for batch execution and has been implemented on a CDC 175 computer operating under NOS 2.4 with a central memory requirement of approximately 115K (octal) of 60 bit words. This program was originally developed in 1983 and later revised in 1988.

Description: WAVDRAG calculates the supersonic zero-lift wave drag of complex aircraft configurations. The numerical model of an aircraft is used throughout the design process from concept to manufacturing. WAVDRAG incorporates extended geometric input capabilities to permit use of a more accurate mathematical model. With WAVDRAG, the engineer can define aircraft components as fusiform or nonfusiform in terms of non-intersecting contours in any direction or more traditional parallel contours. In addition, laterally asymmetric configurations can be simulated. The calculations in WAVDRAG are based on Whitcomb's area-rule computation of equivalent-bodies, with modifications for supersonic speed. Instead of using a single equivalent-body, WAVDRAG calculates a series of equivalent-bodies, one for each roll angle. The total aircraft configuration wave drag is the integrated average of the equivalent-body wave drags through the full roll range of 360 degrees. WAVDRAG currently accepts up to 30 user-defined components containing a maximum of 50 contours as geometric input. Each contour contains a maximum of 50 points. The Mach number, angle-of-attack, and coordinates of angle-of-attack rotation are also input. The program warns of any fusiform-body line segments having a slope larger than the Mach angle. WAVDRAG calculates total drag and the wave-drag coefficient of the specified aircraft configuration. WAVDRAG is written in FORTRAN 77 for batch execution and has been implemented on a CDC CYBER 170 series computer with a central memory requirement of approximately 63K (octal) of 60 bit words. This program was developed in 1983.

Description: An approximate inverse solution is presented for the nonequilibrium flow in the inviscid shock layer about a vehicle in hypersonic flight. The method is based upon a thin-shock-layer approximation and has the advantage of being applicable to both subsonic and supersonic regions of the shock layer. The relative simplicity of the method makes it ideally suited for programming on a digital computer with a significant reduction in storage capacity and computing time required by other more exact methods. Comparison of nonequilibrium solutions for an air mixture obtained by the present method is made with solutions obtained by two other methods. Additional cases are presented for entry of spherical nose cones into representative Venusian and Matrian atmospheres. A digital computer program written in FORTRAN language is presented that permits an arbitrary gas mixture to be employed in the solution. The effects of vibration, dissociation, recombination, electronic excitation, and ionization are included in the program.

Description: Small areas of high heat transfer and pressure can occur on a vehicle surface due to the influence of an impinging shock on the local flow. A method was needed to determine peak pressure and heating of these areas. This package is a system of computer programs designed to calculate two-dimensional shock interference patterns for six types of interference flows. Results also include properties of the inviscid flow field and the inviscid-viscous interaction at the surface along with peak pressure and peak heating at the impingement point. The six types of interference flow patterns considered are: 1) Type I interference patterns, occurring when two weak shocks of opposite families, BS (bow shock) and IS (impingment shock), intersect when the flow upstream of the impingement point is supersonic, or in the case of a blunt body, takes place well below the sonic point. 2) Type II interference pattern occurs when two shocks of opposite families (bow shock and impinging shock) intersect. Both shocks are weak as in type I, but are of such strength that in order to turn the flow, a Mach reflection must exist in the center of the flow field with an embedded subsonic region occurring between the intersection points (A & B) and the accompanying shear layers. Type II interference occurs on a blunt body when the impinging shock intersects the bow shock near the sonic point. 3) Type III shock interference pattern occurs when a weak impinging shock intersects a strong detached bow shock. On a blunt body the shock intersection occurs near or above the lower sonic point. 4) Type IV interference can occur when the impinging shock intersects a strong bow shock ahead of a subsonic flow region. On a blunt body this shock intersection is located between the lower sonic point and just above the body axis. The impinging shock causes a displacement of the bow shock and the formation of a supersonic jet that is embedded in the subsonic region. A jet bow shock is produced when the jet impinges on the surface, creating a small region with high stagnation heating. 5) Type V interference involves the interaction of two weak shocks of the same family. The interaction produces a shear layer, a supersonic jet, and a transmitted impinging shock. On a blunt body the shock interaction occurs near the upper sonic point. 6) Type VI interference involves the intersection of two weak shocks of the same family, which leads to an entirely supersonic flow field. This type of interference is important because it provides a means for predicting the onset of type V. Peak-heating correlations for laminar and turbulent shock-boundary-layer interactions are included in the programs for types I, II, V, and VI interference patterns. Heating correlations for laminar and turbulent reattaching shear layers obtained from separation studies are included in the program for type III interference. This program is written in FORTRAN IV for batch execution and has been implemented on a CDC 6000 Series computer. This program was developed in 1973.

Description: A computer program has been written to obtain the wave and friction drag of configurations with bodies of revolution and fins. These inviscid flow fields are superimposed and the wave drag of the configuration is obtained by integration of the surface pressures. The friction drag is obtained from the viscous flow field of the body and a flat-plate friction analysis of the fins. The numerical solution of these flow fields, superposition, and integration to obtain total drag have been programmed for high-speed digital computation. A large portion of the input required by the program is involved with the description of the configuration geometry and the specific surface positions where pressures are to be evaluated. In addition to drag forces, an output is available whereby the pressure distributions on the body and fins can be obtained.

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: This paper is on the control of nonlinear-nonstationary vibration of a frame-stringer structure resulting from high levels of excitaation from a nearby supersonic jet exhaust. The structure exhibits periodic, chaotic, or random behaviors when forced by high-intensity sound from a supersonic jet exhaust with shock loading superimposed on the broadband response. The time history of the pressure, showing the rotation and flapping of the shock structure in the jet column due to large-scale instabilities, indicates that the response is not only nonlinear but also nonstationary. The acoustic pressure radiated by the structure also contains shocks and the formation of harmonics with distance. Control of the structural response is achieved by actively forcing the structure with an actuator at the shock oscillation frequency whose amplitude is locked into a self-control cycle. Results show that the peak power level is reduced by a factor of 63, or 18 dB. As a result, new broadband components emerge with at least four harmonics. At accelerating and decelerating supersonic speeds, the exhaust from the jet induces higher transient loading on the nearby flexible structure due to the occurence of multiple shocks from the jet.

Description: Prandtl's theory is used to determine the airflow over bodies and wings adapted to supersonic flight. By making use of these results, and by incorporating in them an allowance for the probable skin friction, some estimates of expected lift-drag ratios are made for various flight speeds with the best configuration. At each speed a slender body and wings having the best angle of sweepback are considered. For the range of supersonic speeds shown an airplane of normal density and loading would be required to operate at an altitude of the order of 60,000 feet. The limiting value of 1-1/2 times the speed of sound corresponds to a flight speed of 1000 miles per hour. At this speed about 1.5 miles per gallon of fuel are expected. It is interesting to note that this value corresponds to a value of more than 15 miles per gallon when the weight is reduced to correspond to that of an ordinary automobile.

Description: In theory, the most efficient wing shape for transonic and low supersonic speeds is simply a long narrow straight subsonic wing turned at an oblique angle to the flight direction. This theory has been verified by tests at Mach numbers from .6 to 1.4 in supersonic wind tunnel and by comparative studies of transonic transport designs.

Description: Recent theoretical and experimental work in supersonic aerodynamics is reviewed with its practical application in mind. Several arrangements of supporting surfaces and bodies are discussed and in some cases comparisons of theory and experiment are made. Finally, certain phenomena connected with lift and drag in a rarefied medium are considered briefly.

Description: A method is reported for determining mathematically the combined disturbance field, and in certain cases the minimum drag, of wings at supersonic speeds. The simplest analytic example is provided by the wing of elliptic planform, which achieves its minimum drag when the lift is distributed uniformly over the surface. With a symmetrical distribution of thickness, the requirement of minimum drag for a given total volume is found to lead to profiles of constant curvature.

Description: The assumptions of the thin airfoil theory are found to provide certain necessary conditions for the minimum drag of airfoils having a given total lift, a given maximum thickness, or a given volume. The conditions are applicable to steady or unsteady motions and to subsonic or supersonic speeds without restriction on the planform. The computation of drag and the statement of the conditions for minimum drag depend on the consideration of a combined flow field, which is obtained by superimposing the disturbance velocities in forward and reversed motions. If the planform of the airfoil and its total lift are given, it is found that, for minimum drag, the lift must be distributed in such a way that the downwash in the combined field is constant over the entire planform. If the planform is given and the thickness of the airfoil is required to contain a specified volume, then the thickness must be distributed over the planform in such a way that the pressure gradient of the combined field in the direction of flight is constant at all points of the wing.

Description: The application of mathematical advances made in electricity and other branches to problems of airplane dynamics is demonstrated. The Heaviside-Bromwich methods of solution of linear differential equations are described and it is shown how these methods avoid the consideration of boundary conditions and of particular or complementary integrals. It is pointed out that if the solution of the differential equation is obtained for the case of a unit disturbance, the effect of varying disturbances may be found therefrom by Carson's theorem. A graphical solution of Carson's integral for irregular disturbances is given. The procedure of obtaining unit solutions of the equations is then taken up and the analogy between Heaviside's symbolic series solution and a physical procedure of approximation is shown. It is suggested that a fictitious impulsive disturbance be used in the treatment of initial motions.

Description: Accounts are presented of tangential, semisubmerged, and internal store carriage drag, as well as of the results to date of stores-separation investigations employing both computational and experimental methods. It is demonstrated that CFD can contribute to the data base needed for internal stores carriage design and trade studies. Attention is given to the case of an internally carried store separating from its bay at supersonic speeds, where various kinds of interference are encountered. A code is developed for the simulation of these unsteady flows through time-accurate computations; computation results are noted to closely reproduce experiment data.

Description: A development status evaluation is presented for the aerodynamics of missile configurations with noncircular cross-sections and bank-to-turn maneuvering systems, giving attention to cases with elliptical and square cross-sections, as well as bodies with variable cross-sections. The assessment of bank-to-turn missile performance notes inherent stability/control problems. A summary and index are provided for aerodynamic data on monoplanar configurations, including those which incorporate airbreathing propulsion systems.

Description: A survey is conducted of the results of investigations into the flowfields and aerodynamic forces associated with low aspect ratio wings at high angles of attack. Attention is given to criteria for the cataloging of these flowfields, the phenomenon of vortex breakdown, with varying wing incidence and Mach number, and the effects of aspect ratio and compressibility. The planforms treated are of rectangular, clipped-delta, and strake-wing combination geometries. Extensive graphic representations of performance trends with varying parameters are furnished.

Description: A central artificial-viscosity and an upwind-biased difference method are contrived to solve the Euler equations for flowfields over typical spacecrafts. The spatial discretization is based on either nodal or cell-vertex formulation in the domain extending from free stream to the end of the vehicle. The outer boundary is treated as a bow shock in the first method but is placed in the free stream in the second, which captures both bow and internal shocks using an approximate Riemann solver based on high-order extrapolation to the cell face. These methods were tested for the Shuttle and Hermes orbiters at wind-tunnel conditions and angles of attack ranging from 0 to 60 deg. The artificial-viscosity method incorporated with a shock-fitting procedure shows smeared crossflow and wing-shock positions and required 15 percent more CPU per node than the upwind method. Greater flexibility and robustness is demonstrated by the latter on a fixed grid for all cases considered.

Description: Three-dimensional Navier-Stokes simulations have been done for transonic and low supersonic flow past a nonaxisymmetric nozzle typical of those advocated for advanced fighter airplanes. The jet exhaust is included in the calculations. The investigation compares the performance of the unmodified Baldwin-Lomax turbulence model with its performance when enhanced by the Degani-Schiff and the Goldberg modifications. Solutions are presented for Mach numbers of 0.80, 0.94, and 1.20 at 0-deg angle of attack and a Reynolds number of 20 x 10 to the 6th. The numerical results, which are compared to the wind-tunnel data, show that the three turbulence models predict considerably different shock locations, separated-flow regions, and flowfields.

Description: A new kind of C-type grid is proposed, this grid is non-periodic on the wake and allows minimum skewness for cascades with high turning and large camber. Reynolds-averaged Navier-Stokes equations are solved on this type of grid using a finite volume discretization and a full multigrid method which uses Runge-Kutta stepping as the driving scheme. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. A detailed numerical study is proposed for a highly loaded transonic blade. A grid independence analysis is presented in terms of pressure distribution, exit flow angles, and loss coefficient. Comparison with experiments clearly demonstrates the capability of the proposed procedure.

Description: The flow through a highly offset subsonic diffuser with cross-sectional profiles that varied from rectangular at the duct entrance to circular at the engine face was numerically simulated. A multizonal approach combined with a two-grid topology was used to represent both the internal and external flowfields, and an implicit, approximately-factored, partially flux-split finite-difference algorithm was used to solve the three-dimensional thin-layer Navier-Stokes equations. The computed static pressures along the inlet wall and total pressures on the engine face were compared with experimental data. In addition, the overall flowfield within the duct was examined in detail. Good agreement is shown between experiment and computations, with the limiting factor being the lack of a reliable turbulence model for internal flow problems.

Description: This paper discusses NASA's Computational Aerosciences (CAS) Project of the High Performance Computing and Communications Program (HPCCP). The project is aimed at developing advanced, multidisciplinary simulation capabilities for aerospace vehicle and propulsion system design. It is also aimed at overcoming computational performance barriers by accelerating the development of parallel computer technology. The goals and approach of the CAS Project are described and the challenges to its implementation are addressed. Specific vehicle class simulations to be demonstrated and the principal mutidisciplinary modeling approaches to be emphasized are described. The computational speed and memory requirements for representative multidisciplinary applications are estimated. Finally, the state of parallel computer technology including programming issues and the results of performance measurements are explored.

Description: Previous work on the use of a vortex trapped above a wing in order to produce high lift at low angles of attack is extended here. It is first postulated that the optimum way to trap a vortex is to design the airfoil section and wing so that the flow along the vortex core is minimized. It is then shown that a vertical fence both in front of and behind the separation bubble generated by the trapped vortex is an effective way to reduce the mass flow removal and its associated drag to a negligible amount. In order to show that vertical surfaces upstream and downstream of the vortex separation bubble have an opposite effect on the source requirements for vortex trapping, conformal mapping methods are used to obtain the solutions for a variety of simple two-dimensional, inviscid, incompressible flow configurations. Trapped-vortex flowfield solutions for the flow over flat plate and Clark-Y airfoils are then used to demonstrate that the heights of the fences can be tailored to make the required mass withdrawal (and therefore, the drag due to trapping) to be vanishingly small.

Description: Viscous, heat-conducting flow with chemical and vibrational relaxation processes of the constituent gases surrounding winged spacecraft is considered in the continuum regime. The Navier-Stokes equations are appended by additional vibrational energy and species rate equations and supplemented by the equations of state and the phenomenological laws based on mixture rules or collisional cross sections. Numerical convective flux can be obtained from several forms of one-dimensional Riemann solver, with or without entropy correction. High-order accuracy is obtained from two types of reconstructive interpolation. A number of explicit and implicit numerical schemes have been implemented as a means to yield converged solutions. Both shock-fitting, finite-difference and shock-capturing, finite-volume techniques have been tested for configurations such as a sphere, double ellipsoid, blunt-edge delta wing, a European Hermes vehicle, and the U.S. Shuttle Orbiter. The shock-fitting code provides excellent results only for simple configurations, whereas the shock-capturing code leads to overall satisfying solutions for complex geometries.

Description: The effect of nonequilibrium thermochemical processes on the aerodynamic and heat transfer characteristics of a hypersonic vehicle is illustrated with examples. It is shown that the conventional method of predicting chemical reactions always predict that the flow is closer to equilibrium than it actually is, leading to incorrect predictions of aerodynamic characteristics of a vehicle. A method is presented for predicting chemical processes using a two-temperature kinetic model, which is developed on the basis of combining conservation equations for vibrational energy and for electron electronic energy, which are also derived. The validity of the two-temperature model is proven using two types of experimental data: radiation and shock shapes.

Description: The meaning of catalysis and its relation to aerodynamic heating in nonequilibrium hypersonic flows are discussed. The species equations are described and boundary conditions for them are derived for a multicomponent gas and for a binary gas. Slip effects are included for application of continuum methods to low-density flows. Measurement techniques for determining catalytic wall recombination rates are discussed. Among them are experiments carried out in arc jets as well as flow reactors. Diagnostic methods for determining the atom or molecule concentrations in the flow are included. Results are given for a number of materials of interest to the aerospace community, including glassy coatings such as the RCG coating of the Space Shuttle and for high temperature refractory metals such as coated niobium. Methods of calculating the heat flux to space vehicles in nonequilibrium flows are described. These methods are applied to the Space Shuttle, the planned Aeroassist Flight Experiment, and a hypersonic slender vehicle such as a transatmospheric vehicle.

Description: Turbulence modeling for high-speed compressible flows is described and discussed. Starting with the compressible Navier-Stokes equations, methods of statistical averaging are described by means of which the Reynolds-averaged Navier-Stokes equations are developed. Unknown averages in these equations are approximated using various closure concepts. Zero-, one-, and two-equation eddy viscosity models, algebraic stress models, and Reynolds stress transport models are discussed. Computations of supersonic and hypersonic flows obtained using several of the models are discussed and compared with experimental results. Specific examples include attached boundary-layer flows, shock-wave boundary-layer interactions, and compressible shear layers. From these examples, conclusions regarding the status of modeling and recommendations for future studies are discussed.

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 linearized flow theory, certain very interesting extremal properties of wings can be derived under rather broad conditions without the use of a complicated mathematical apparatus. The present chapter reviews certain results of this theory and indicates some rather obvious extensions to incorporate various auxiliary conditions. Several examples illustrating the relation between the geometrical features of the wing and the lift distribution for minimum drag are given.

Description: The items discussed are: (1) a recently proposed correction formula for the effect of compressibility in two dimensional subsonic flow; (2) the equivalence rule and the area rule for transonic speeds; (3) reciprocal relations in linearized wing theory; and (4) some general results connected with the problem of minimum wave resistance. The paper concludes with an example showing indentation of the fuselage to obtain favorable interference with the wing at supersonic speeds.

Description: In the wing section theory the magnitude of the circulation, and hence of the lift, is determined by the velocity that would be induced near the trailing edge of the section in a non-lifting potential flow. In three dimensional flow the problem is complicated by the presence of the wake and no simple basic solution has been found. Treatment of the problem of a wing of finite span is reported on the basis of the two dimensional theory, corrected for the effect of the wake.

Description: In theory, antisymmetric arrangements of wings and bodies can have smaller wave drag than corresponding mirror-symmetric arrangements. Thus, a long narrow oblique wing which presents the same aspect for two opposite directions of flight is potentially more efficient than corresponding (i.e., structurally equivalent) swept wing. The single continuous wing panel also adapts itself more readily to varying angles of obliquity, and hence, to varying flight speeds. Previous work on the aerodynamics and flight stability of oblique wing combinations is reviewed and a possible mode of application to transport aircraft operating at moderate supersonic speeds is suggested.

Description: It is shown that the drag of any semi-infinite airfoil section in purely subsonic inviscid flow follows precisely the Prandtl-Glauert compressibility rule. The result for the parabola has application to leading edge corrections in thin airfoil theory.

Description: Comparisons of wing-body combinations may not disclose the full effect of a loss in aerodynamic efficiency. If the thrust needs to be increased at a given altitude then more or larger engines will have to be used and the possibility of concealing them becomes less. In this process the lift drag ratio of the complete airplane may become still more unfavorable than indicated by the comparison. Primarily aerodynamic and structural considerations point toward the development of turbojet engines specifically adapted to operation in an atmosphere of one tenth normal density. In addition to the numerous other technological problems associated with operation at these high altitudes, the problems of safe descent and effective limitation to low speeds at low altitudes seem important.

Description: If the thin airfoil theory is applied to an airfoil having a rounded leading edge, a certain error will arise in the determination of the pressure distribution around the nose. It is shown that the evaluation of the drag of such a blunt nosed airfoil by the thin airfoil theory requires the addition of a leading edge force, analogous to the leading edge thrust of the lifting airfoil. The method of calculation is illustrated by application to: (1) The Joukowski airfoil in subsonic flow; and (2) the thin elliptic cone in supersonic flow. A general formula for the edge force is provided which is applicable to a variety of wing forms.

Description: Some of the recent advances in the theory of thin airfoils are presented with particular reference to extensions of the theory to three dimensional flows and to supersonic speeds. The problem discussed herein is the calculation of the small disturbance velocities u, v, and w in the external field produced by the flight velocity V of the airfoil.

Description: Three optimization-based methods for solving aerodynamic design problems are compared. The Euler equations for one-dimensional duct flow was used as a model problem, and the three methods are compared for efficiency, robustness, and implementation difficulty. The smoothness of the design problem with respect to different shock-capturing finite difference schemes, and in the presence of grid refinement, is investigated.