ALBERT

Description: NASA is directing research to develop technology for a high-speed civil transport. Supersonic laminar flow control has been identified as a program element, since it offers significant drag-reduction benefits and is one of the more promising technologies for producing an economically viable aircraft design. NASA is using two prototype F-16XL aircraft to research supersonic laminar flow control. The F-16XL planform is similar to design planforms of high-speed civil transports. The planform makes the aircraft ideally suited for developing technology pertinent to high-speed transports. The supersonic laminar flow control research programs for both aircraft are described. Some general results of the ship-1 program demonstrate that significant laminar flow was obtained using laminar flow control on a highly swept wing at supersonic speeds.

Description: Issues and advances in current hypersonic flow research perceived to be of interest in theoretical fluid/gas dynamics are reviewed. Particular attention is given to the hypersonic aircraft as waverider, computational methods and theoretical development in the study of viscous interaction, and boundary-layer instability and transition studies. In the present framework the study of viscous hypersonic flow faces transition problems of two kinds which represent the two major areas of current research: the turbulence transition in the high Re range and the transition to the free-molecule limit.

Description: The temporal development of a 2D viscous incompressible flow generated by a circular cylinder started impulsively into steady rotatory and rectilinear motion is studied by integration of a velocity/vorticity formulation of the governing equations, using an explicit finite-difference/pseudo-spectral technique and an implementation of the Biot-Savart law. Results are presented for a Reynolds number of 200 (based on the cylinder diameter 2a and the magnitude U of the rectilinear velocity) for several values of the angular/rectilinear speed ratio alpha = omega(a)/U (where omega is the angular speed) up to 3.25. Several aspects of the kinematics and dynamics of the flow not considered earlier are discussed. For higher values of alpha, the results indicate that for Re = 200, vortex shedding does indeed occur for alpha = 3.25. However, consecutive vortices shed by the body can be shed from the same side and be of the same sense, in contrast to the nonrotating case, in which mirror-image vortices of opposite sense are shed alternately on opposite sides of the body. The implications of the results are discussed in relation to the possibility of suppressing vortex shedding by open or closed-loop control of the rotation rate.

Description: A total variation diminishing (TVD) scheme has been developed and incorporated into an existing time-accurate high-resolution Navier-Stokes code. The accuracy and the robustness of the resulting solution procedure have been assessed by performing many calculations in four different areas: shock tube flows, regular shock reflection, supersonic boundary layer, and shock boundary layer interactions. These numerical results compare well with corresponding exact solutions or experimental data.

Description: The linear stability of the trailing line vortex model of Batchelor (1964) is studied using a spectral collocation and matrix eigenvalue method. The entire unstable region in the swirl/axial wavenumber parameter space is mapped out for various azimuthal wavenumbers for both the inviscid and viscous stability problem. The results of the study provide a direct numerical validation of the large-azimuthal-wavenumber asymptotic analysis of Leibovich and Stewartson (1983). It is shown that accurate results are obtained up to azimuthal wavenumbers of 10,000 and greater, and the agreement with the asymptotic theory is excellent.

Description: In spite of many attempts at modeling natural transition, it has not been possible to predict the streamwise intensities. A procedure is developed which incorporates some results of linear stability theory into one-equation and stress model formulations. The stresses resulting from fluctuations in the transitional region have turbulent, laminar (nonturbulent) and large eddy components. Comparison with Schubauer and Klebanoff's experiments have shown that the nonturbulent and large eddy components have a large influence on the streamwise intensities and little influence on the shear stress. Finally, predictions of the one-equation model were as good as those obtained by the stress model.

Description: A two-dimensional, nonhydrostatic, elastic numerical model has been used to study the generation of gravity waves for a stably stratified shear flow over an obstacle. When a low-level wind shear is included in the simulation, we find that the predictions for noticeable upstream effects based on Froude number for a uniform flow are no longer accurate. Upstream effects are encountered in the form of upstream propagating columnar disturbances and internal bores away from the obstacle. The limited parameter space studies conducted in this study suggest that the ratio of the shear depth to the obstacle heigh (d/H), the obstacle aspect ratio (H/L), and the Froude number (U/NH) are instrumental in determining the strength and the existence of these upstream disturbances. Thus, the present theoretical and empirical understanding of the importance of the Froude number for determining the nature of upstream effects should be modified substantially to include additional nondimensional parameters when shear is present.

Description: The interaction between a swept shock wave and a laminar boundary layer was investigated experimentally in high-enthalpy hypersonic flow. The effect of high-temperature, real gas physics on the interaction was examined by conducting tests in air and helium. Heat transfer measurements were made on the surface of a flat plate and a shock-generating fin using thin-film resistance sensors for fin incidence angles of 0, 5, and 10 deg at Mach numbers of 6.9 in air and 7.2 in helium. The experiments were conducted in the NASA HYPULSE expansion tube, an impulse-type facility capable of generating high-enthalpy, high-velocity flow with freestream levels of dissociated species that are particularly low. The measurements indicate that the swept shock wave creates high local heat transfer levels in the interaction region, with the highest heating found in the strongest interaction. The maximum measured heating rates in the interaction are order of magnitude greater than laminar flat plate boundary layer heating levels at the same location.

Description: An experimental study has been conducted to examine the flow field of the 3D crossing shock wave/turbulent boundary layer interaction. A symmetric pair of 9-deg fins were used to generate the crossing shocks. The incoming boundary layer was developed on the tunnel sidewall and thus was relatively thick, 0.49 arcsec, and suited for pitot probe surveys. The test conditions were a nominal Mach number of 3 and unit Reynolds number of 1.2 x 10 exp 7/ft. The measurements obtained included surface oil flow visualizations, surface static pressures, and boundary layer pitot pressure profiles. The results showed that downstream of the crossing shock intersection, the stagnation pressure losses were significant and the stagnation pressure profiles were highly nonuniform. Despite the severe shock disturbances, the law of the wall and the law of the wake were found to give relatively good agreement with the experimental data.

Description: A 3D CFD algorithm is used to study the effect of thermal and chemical nonequilibrium on slender and blunt body aerothermodynamics. Both perfect gas and reacting gas air models are used to compute the flow over a generic transatmospheric vehicle and a proposed lunar transfer vehicle. The reacting air is characterized by a translational-rotational temperature and a vibrational-electron-electronic temperature and includes eight chemical species. The effects of chemical reaction, vibrational excitation, and ionization on lift-to-drag ratio and trim angle are investigated. Results for the NASA Ames All-body Configuration show a significant difference in center of gravity location for a reacting gas flight case when compared to a perfect gas wind tunnel case at the same Mach number, Reynolds number, and angle of attack. For the same center of gravity location, the wind tunnel model trims at lower angle of attack than the full-scale flight case. Nonionized and ionized results for a proposed lunar transfer vehicle compare well to computational results obtained from a previously validated reacting gas algorithm. Under the conditions investigated, effects of weak ionization on the heat transfer and aerodynamic coefficients were minimal.

Description: The Direct Simulation Monte Carlo (DSMC) method is applied to a radiating, hypersonic, axisymmetric flow over a blunt body in the near continuum regime. The ability of the method to predict the flowfield radiation and the radiative heating is investigated for flow over the Project Fire II configuration at 11.36 kilometers per second at an altitude of 76.42 kilometers. Two methods that differ in the manner in which they treat ionization and estimate electronic excitation are employed. The calculated results are presented and compared with both experimental data and solutions where radiation effects were not included. Differences in the results are discussed. Both methods ignore self absorption and, as a result, overpredict measured radiative heating.

Description: An implicit finite element based algorithm for the compressible Navier-Stokes equations is outlined, and the solution of the resulting equation by a line relaxation on general meshes of triangles or tetrahedra is described. The problem of generating and adapting unstructured meshes for viscous flows is reexamined, and an approach for both 2D and 3D simulations is proposed. An efficient approach appears to be the use of an implicit/explicit procedure, with the implicit treatment being restricted to those regions of the mesh where viscous effects are known to be dominant. Numerical examples demonstrating the computational performance of the proposed techniques are given.

Description: An overview is presented of the most compelling technological and economic arguments for NASA's agressive coordination of an SST-development program that would enlist all available U.S. aerospace industry resources. Attention is given to the minimization of upper atmosphere pollution through the use of low-NO(x) emission combustors and the reduction of sonic boom through wing/fuselage optimization. It is projected that a successful SST program would boost U.S. civil aircraft market share to nearly 80 percent; this represents the creation of 140,000 new jobs.

Description: In a hypersonic boundary layer over a wall of variable curvature, the region most susceptible to Goertler vortices is the temperature adjustment layer sitting at the edge of the boundary layer. This temperature adjustment layer is also the most dangerous site for Reyleigh instability. We investigate how the existence of large amplitude Goertler vortices affects the growth rate of Rayleigh instability. The effects of wall cooling and gas dissociation on this instability are also studied. We find that all these mechanisms increase the growth rate of Rayleigh instability and are therefore destabilizing.

Description: The flowfield structure of a range of symmetric crossing-shock wave/turbulent boundary-layer interactions of varying strength is presented. The test geometry, consisting of a symmetric pair of opposing sharp fins at angle of attack, alpha, mounted to a flat plate, is studied experimentally for a range of alpha from 7 to 15 degrees at Mach numbers of 3 and 4. Results reveal that the basic flowfield shock structure remains similar in nature over the range of interaction strengths examined, with the only changes being in the scale and location of the various features present. The separated flow regions are classified as being either completely or partially separated, the completely separated case being the one in which the entire incoming boundary layer separates from the plate surface. For the current experiments, all but the weakest of the interactions exhibited complete boundary layer separation. Finally, the effects of model geometry are analyzed by comparing data for shock generators of varying lengths, with the results showing no evidence of upstream influence due to the shock generator trailing edges.

Description: The results of a joint experimental and computational study on the flowfield over a periodically pitched NACA0012 airfoil, and the resultant lift variation, are reported in this paper. The lift variation over a cycle of oscillation, and hence the lift hysteresis loop, is estimated from the velocity distribution in the wake measured or computed for successive phases of the cycle. Experimentally, the estimated lift hysteresis loops are compared with available data from the literature as well as with limited force balance measurements. Computationally, the estimated lift variations are compared with the corresponding variation obtained from the surface pressure distribution. Four analytical formulations for the lift estimation from wake surveys are considered and relative successes of the four are discussed.

Description: A practical solution algorithm for steady 3D Euler flows is presented. This algorithm employs coupling of a surface triangulator, an automatic tetrahedral mesh generator, an unstructured grid flow solver, and an error estimation procedure. The performance of the method is illustrated using a shock interaction problem in high Mach number flow over a swept circular cylinder.

Description: The Baldwin-Lomax model is used in many CFD codes because it is quick and easy to implement. In this paper, we discuss implementing the Baldwin-Lomax turbulence model for both steady and unsteady compressible flows. In addition, these flows may be either separated or attached. In order to apply this turbulence model to flows which may be subjected to these conditions, certain modifications should be made to the original Baldwin-Lomax model. We discuss these modifications and determine whether the Baldwin-Lomax model is a viable turbulence model that produces reasonably accurate results for high speed flows that can be found in engine inlets.

Description: The use is considered of a multigrid method with central differencing to solve the Navier-Stokes equations for high speed flows. The time dependent form of the equations is integrated with a Runge-Kutta scheme accelerated by local time stepping and variable coefficient implicit residual smoothing. Of particular importance are the details of the numerical dissipation formulation, especially the switch between the second and fourth difference terms. Solutions are given for 2-D laminar flow over a circular cylinder and a 15 deg compression ramp.

Description: A unified finite element algorithm is developed which is applicable to a wide range of problems of fluid mechanics without recourse to artificial, empirically determined factors. In its explicit form, the algorithm is similar to the Taylor-Galerkin scheme and is easily adopted to standard codes. The scheme proposed here possesses sufficient natural balancing diffusion and thus reduces and sometimes eliminates the need for special 'shock capturing' diffusion. The efficiency of the algorithm is demonstrated using several examples ranging from incompressible through transonic regions to supersonic flows.

Description: Examination of the literature shows that the comparison between experiment and computation for highly separated unswept compression ramp flows is generally poor, irrespective of the turbulence model used. In general, the upstream influence is not correct, the wall pressure rise through separation is too steep, and the pressures under the separated shear layer are too high. In the current study, the objective is to determine if these discrepancies might be attributed more to other factors such as flowfield unsteadiness or three-dimensionality, rather than to inadequate turbulence modeling. To examine this possibility, multichannel wall pressure fluctuations were measured under the unsteady separation shock wave in a 28-deg unswept compression ramp flow at Mach 5. The results show that the large scale, low frequency separation shock unsteadiness controls the distribution of time-averaged surface properties and that neglect of the unsteadiness is probably the primary cause of the discrepancy between experiment and computation.

Description: Adaptive wall research at the University of Southampton has been directed towards the development of testing techniques for use in nonporous test sections where two flexible walls are profiled in single curvature. This paper highlights the recent advances that have been made in the testing of 2D airfoils through the speed of sound and the testing of 3D models at high subsonic speeds. Techniques have been developed to accommodate the variety of flow regimes encountered in near sonic airfoil tests. The experimental evidence to date suggests that the new techniques coupled with established procedures allow airfoil data, free from top and bottom wall interference, to be gathered from adaptive flexible walled test sections throughout the entire subsonic, transonic and supersonic speed ranges. Techniques applicable to the testing of 3D models have evolved primarily from experience gained by testing sidewall mounted half-wings. Emphasis has been placed upon models with planforms similar to those of current transport wings. Techniques for high subsonic speeds have now been developed to the point where the residual levels of interference are low.

Description: Computational and experimental studies are conducted to investigate the influence of a trailing edge jet on flow separation and subsequent vortex formation over steady and accelerated airfoils at high angles of attack. A computer code, employing the stream function-vorticity approach, is developed and utilized to conduct numerical experiments on the flow problem. To verify and economize such efforts, an experimental system is developed and incorporated into a subsonic wind tunnel where streamline and vortex flow visualization experiments are conducted. The study demonstrates the role of the trailing edge jet in controlling flow separation and subsequent vortex development for steady and accelerating flow at angles past the static stall angle of attack. The results suggest that the concept of the trailing edge jet may be utilized to control the characteristics of unsteady separated flows over lifting surfaces. This control possibility seems to be quite effective and could have a significant role in controlling unsteady separated flows.

Description: A new flux splitting scheme is proposed. The scheme is remarkably simple and yet its accuracy rivals and in some cases surpasses that of Roe's solver in the Euler and Navier-Stokes solutions performed in this study. The scheme is robust and converges as fast as the Roe splitting. An approximately defined cell-face advection Mach number is proposed using values from the two straddling cells via associated characteristic speeds. This interface Mach number is then used to determine the upwind extrapolation for the convective quantities. Accordingly, the name of the scheme is coined as Advection Upstream Splitting Method (AUSM). A new pressure splitting is introduced which is shown to behave successfully, yielding much smoother results than other existing pressure splittings. Of particular interest is the supersonic blunt body problem in which the Roe scheme gives anomalous solutions. The AUSM produces correct solutions without difficulty for a wide range of flow conditions as well as grids.

Description: This is the second of two papers on the interaction between a longitudinal vortex pair, produced by a delta-wing at angle of attack, and a turbulent boundary layer developing on a flat plate. In the first paper only the outer parts of the vortices entered the boundary layer whereas in this paper the vortices merge with it. In the resultant interaction, the boundary layer between the vortices is kept thin by lateral divergence and a three-dimensional separation line is formed outboard of each vortex. Turbulent, momentum-deficient fluid containing longitudinal vorticity is entrained from the boundary layer along these lines and wrapped around the vortices. As a consequence, the turbulent region of the vortices increases in size and the circulation slowly decreases. It is shown that the flow near the separation line and in the vortices is complicated, and this interaction is expected to be more difficult to calculate than the first. Detailed mean flow and turbulence measurements are reported.

Description: A detailed comparison is made between Navier-Stokes and DSMC calculations for flows near the continuum limit to assess the accuracy of the continuum equations in this regime. Meaningful comparisons require the use of similar physical models. This necessitates the inclusion of a separate rotational energy equation and use of slip boundary conditions. Inclusion of slip boundary conditions resulted in improved agreement between surface properties. Moreover, good agreement was obtained for the various temperatures in the nonequilibrium portion of the flow field that does not contain the shock region. Departures are noted in the shock region and in regions where thermal diffusion effects are important.

Description: Aerodynamic heating rates are calculated from time-dependent temperature measurements in the vicinity of shock-wave boundary-layer interactions due to conical compression ramps on an axisymmetric body. The data were acquired at the Ohio State University Aeronautical and Astronautical Research Laboratory and at the Air Force Flight Dynamics Laboratory at Mach numbers of 6 and 10. The model is a cylindrical body with a 10 deg conical nose. Conical ramps with half-angles of 10, 20, 25, 30, and 35 deg serve as shock-wave generators. Flowfield surveys are made in the vicinity of the ramp vertices, separation points, and reattachment points. Experimental results quantify temperature response and the resulting heat transfer rates as a function of ramp angle, Reynolds number and freestream Mach number. The temperature responses within the flowfield appear to be steady-state for all angles and all Reynolds numbers, and hence, the heat transfer rates appear to be steady-state.

Description: A combined coarse-grid correction/upwind relaxation strategy to provide rapid convergence for 3D high-speed viscous flowfields is discussed an evaluated. The construction and analysis of a simple two-grid acceleration procedure based on 'hyperbolic' multigrid concepts is presented. Numerical simulations of a 2D compression-corner flowfield, a 3D crossing shock/turbulent boundary layer interaction, and a 3D scramjet inlet flowfield are presented to illustrate the benefits of the approach. Results indicate that the procedure generally converges two or more times faster than the baseline algorithm.

Description: A computational study has been performed of sharp fin-induced swept shock wave/turbulent boundary layer interactions at low hypersonic Mach numbers. The objective was to determine if results obtained using a conical Navier-Stokes code, particularly the peak heating and pressure, are adequate for engineering predictions. The advantage of the conical approach is that the problem becomes two-dimensional and requires much less computational effort than a fully three-dimensional calculation. In this code, the standard Baldwin-Lomax model is used for turbulent closure and its performance is studied in some detail. To assess the approach interactive flowfields generated by unswept sharp fins at two angles of attack at each of three hyprsonic freestream Mach numbers (5, 6, 11) have been calculated and the results compared with experimental wall pressure and heat transfer data. Although the conical Navier-Stokes/Baldwin-Lomax approach is reasonably successful at Mach numbers up to 5, the performance deteriorates as the Mach number is increased. Nevertheless, the approach could be a valuable tool in preliminary parametric design studies.

Description: A numerical study was conducted to analyze the performance of different turbulence models when applied to the hypersonic NASA P8 inlet. Computational results from the PARC2D code, which solves the full two-dimensional Reynolds-averaged Navier-Stokes equation, were compared with experimental data. The zero-equation models considered for the study were the Baldwin-Lomax model, the Thomas model, and a combination of the Baldwin-Lomax and Thomas models; the two-equation models considered were the Chien model, the Speziale model (both low Reynolds number), and the Launder and Spalding model (high Reynolds number). The Thomas model performed best among the zero-equation models, and predicted good pressure distributions. The Chien and Speziale models compared wery well with the experimental data, and performed better than the Thomas model near the walls.

Description: An investigation into the numerical simulation of unsteady flows is undertaken using a two-stage Runge-Kutta scheme coupled with the dynamic solution-adaptive grid algorithm developed by the authors. The inviscid fluxes are described by a modified Advective Upwind Split Method to eliminate the need for artificial dissipation. A well-documented numerical example containing moving discontinuities is presented that demonstrates the ability of the coupled grid/solver scheme to accurately capture unsteady flowfield phenomena. Applications are to a typical inlet diffuser configuration at Mach 3.0 with excessive back pressure inducing inlet unstart.

Description: The hypersonic flow over a cavity is investigated. The time-dependent compressible Navier-Stokes equations are numerically solved. An implicit algorithm, with a subiteration procedure to recover time accuracy, is used to perform the time-accurate computations. The objective of the study is to investigate the effects of Reynolds number and cavity dimensions. The comparsion of the computations with available experimental data, in terms of time mean static pressure, heat transfer, and Mach number, show good agreement. In the computations large vortex structures, which adversely affect the cavity flow characteristics, are observed at the rear of the cavity. A self-sustained oscillatory motion occurs within the cavity over a range of Reynolds number and cavity dimensions. The frequency spectra of the oscillations show good agreement with a modified semiempirical relation.

Description: This video discusses how NASA uses large helium-filled balloons to take payloads up 25 miles to the edge of space to gather data. Balloons provide a cost effective approach to reach these heights.

Description: High-angle-of-attack aerodynamic studies have been conducted on both the F18 High Alpha Research Vehicle (HARV) and the X-29A aircraft. Data obtained include on- and off-surface flow visualization and static pressure measurements on the forebody. Comparisons of similar results are made between the two aircraft where possible. The forebody shapes of the two aircraft are different and the X-29A forebody flow is affected by the addition of nose strakes and a flight test noseboom. The forebody flow field of the F-18 HARV is fairly symmetric at zero sideslip and has distinct, well-defined vortices. The X-29A forebody vortices are more diffuse and are sometimes asymmetric at zero sideslip. These asymmetries correlate with observed zero-sideslip aircraft yawing moments.

Description: As a building block in the development of smart lift-enhancement devices, a new concept for flow control using active vortex generators (AVGs) is presented. Ramp, wedge, and doublet wedge (Wheeler) VG configurations are investigated. The AVGs are designed to conform to the surface of the wing section at low alpha. As the section approaches the stall, they are deployed and accordingly, alpha(stall) and C(lmax) are increased. A qualitative analysis of the flow around the various VG configurations was conducted in a low speed wind tunnel at 1.6 ft/s and a Reynolds number of approximately 3400. The results demonstrate that ramp VGs produce vortices that have the longest distance at breakdown. The VGs were also applied to a 25-in. span, 8-in. chord NACA 4415 wing section. Optimization studies were conducted on the spanwise spacing, chordwise position, and size of statically deployed VGs. The test results demonstrate a 14-percent increase in C(lmax) while increasing alpha (stall) by up to 3.

Description: A holographic interferometer has been designed, constructed, and evaluated in an experimental study of the supersonic flow over a rearward facing step. The nominal Mach number at the corner was 2.05 +/- 0.04 and the Reynolds number per inch was 11.9 x 10 exp 6. The holographic interferometric measurements were supplemented by classical measurements of surface pressure, oil flow, and schlieren visualization. The effects of step height and step width were examined. A method to determine the reattachment point from the interferograms was examined and found to be in good agreement with the other measurement techniques. The reattachment point moved closer to the step as the step height was decreased, but its location did not change with varying step width. In addition to providing surface data for the flow over a rearward facing step, this study provides quantitative off-surface density data and Mach number data throughout the flow, obtained from the holographic interferometry measurements, which are suited for code validation.

Description: Uniform high order spectral methods to solve multi-dimensional Euler equations for gas dynamics are discussed. Uniform high order spectral approximations with spectral accuracy in smooth regions of solutions are constructed by introducing the idea of the Essentially Non-Oscillatory (ENO) polynomial interpolations into the spectral methods. The authors present numerical results for the inviscid Burgers' equation, and for the one-dimensional Euler equations including the interactions between a shock wave and density disturbance, Sod's and Lax's shock tube problems, and the blast wave problem. The interaction between a Mach 3 two-dimensional shock wave and a rotating vortex is simulated.

Description: Reynolds number and cowl position effects on the internal shock structure and the resulting performance of a generic three-dimensional sidewall compression scramjet inlet with a leading edge sweep of 45 degrees at Mach 10 have been examined both computationally and experimentally. Prior to the experiment, a three-dimensional Navier-Stokes code was adapted to perform preliminary parametric studies leading to the design of the present configuration. Following this design phase, the code was then utilized as an analysis tool to provide a better understanding of the flow field and the experimental static pressure data for the final experimental configuration. The wind tunnel model possessed 240 static pressure orifices distributed on the forebody plane, sidewalls, and cowl and was tested in the NASA Langley 31 Inch Mach 10 Tunnel.

Description: A time domain numerical scheme is developed to solve for the unsteady flow about a flat plate airfoil due to imposed upstream, small amplitude, transverse velocity perturbations. The governing equation for the resulting unsteady potential is a homogeneous, constant coefficient, convective wave equation. Accurate solution of the problem requires the development of approximate boundary conditions which correctly model the physics of the unsteady flow in the far field. A uniformly valid far field boundary condition is developed, and numerical results are presented using this condition. The stability of the scheme is discussed, and the stability restriction for the scheme is established as a function of the Mach number. Finally, comparisons are made with the frequency domain calculation by Scott and Atassi, and the relative strengths and weaknesses of each approach are assessed.

Description: A method for generating and adaptively refining a highly stretched unstructured mesh suitable for the computation of high-Reynolds-number viscous flows about arbitrary two-dimensional geometries was developed. The method is based on the Delaunay triangulation of a predetermined set of points and employs a local mapping in order to achieve the high stretching rates required in the boundary-layer and wake regions. The initial mesh-point distribution is determined in a geometry-adaptive manner which clusters points in regions of high curvature and sharp corners. Adaptive mesh refinement is achieved by adding new points in regions of large flow gradients, and locally retriangulating; thus, obviating the need for global mesh regeneration. Initial and adapted meshes about complex multi-element airfoil geometries are shown and compressible flow solutions are computed on these meshes.

Description: The propulsive effects of waves in ducts, especially at high Mach numbers, are investigated, focusing on drag and thrust and on the conversion of heat into waves which produce thrust. It is shown that essentially all of the work done by an expanding fluid passing through a duct at high Mach number is delivered in the form of waves, and that duct surface angles exist that are optimum for the production of thrust from a wave. The effects of wave phenomena on drag and thrust are considered by extending the concept of a Busemann biplane into that of a 'Busemann scramjet, taking 'off-design' performance into account. An idealized model of a streamtube with heat addition is developed, and flow mechanisms involved in generating thrust by the expansion of this streamtube in an exhaust nozzle are examined.

Description: A 3-D axisymmetric hypersonic engine inlet was investigated using PARC2D, an ideal gas Computational Fluid Dynamics code. The code was used to predict the results of tests conducted in the Rensselaer Polytechnic Institute Hypersonic Shock Tunnel which measured surface and pitot pressures, and shock positions (through Schlieren photography) at freestream Mach numbers of 10, 13, and 15. A strong viscous/shock interaction was observed in both the experiment and the CFD results, due to the model's parabolic compression ramp. Good agreement was found between the experimental results and the CFD solution both for surface pressures and shock positions. Agreement between pitot pressures was less reliable.

Description: A computational fluid dynamics algorithm is developed for the study of high-pressure axisymmetric hypersonic nozzle flows. The effects of intermolecular forces and vibrational nonequilibrium are included in the analysis. The numerical simulation of gases with an arbitrary equation of state is discussed. Simulations for a high pressure nozzle (p(0) = 138 MPa) demonstrate that both intermolecular forces and vibrational nonequilibrium have a significant affect on the flow. These nonideal effects tend to increase the Mach number at the nozzle exit plane. Thus, they must be included in the design and analysis of high pressure hypersonic nozzles.

Description: Missions to Mars require the successful development of aerobraking technology, and therefore a blunt cone representative of aerobrake shapes is investigated. Ballistic tests of the Pioneer Venus configuration are conducted in carbon dioxide and air at Mach numbers from 7 to 20 and Reynolds numbers from 0.1 x 10 exp 5 to 4 x 10 exp 6. Experimental results show that for defined conditions aerodynamic research can be conducted in air rather than carbon dioxide, providing savings in time and money. In addition, the results offer a prediction of flight aerodynamics during entry into the Martian atmosphere. Also discussed is a comparison of results from two data-reduction techniques showing that a five-degree-of-freedom routine employing weighted least-squares with differential corrections analyzes ballistic data more accurately.

Description: This paper discusses the design and testing of candidate Advanced Launch System (ALS) Propulsion and Avionics (P/A) Module configurations. The P/A Module is a key element of future launch systems because it is essential to the recovery and reuse of high-value propulsion and avionics hardware. The ALS approach involves landing of first stage (booster) and/or second stage (core) P/A modules near the launch site to minimize logistics and refurbishment cost. The key issue addressed herein is the aerodynamic design of the P/A module, including the stability characteristics and the lift-to-drag (L/D) performance required to achieve the necessary landing guidance accuracy. The reference P/A module configuration was found to be statically stable for the desired flight regime, to provide adequate L/D for targeting, and to have effective modulation of the L/D performance using a body flap. The hypersonic aerodynamic trends for nose corner radius, boattail angle and body flap deflections were consistent with pretest predictions. However, the levels for the L/D and axial force for hypersonic Mach numbers were overpredicted by impact theories.

Description: The Direct Simulation Monte Carlo method of Bird is used to investigate the characteristics of low density hypersonic flowfields for typical aerobrakes during Martian atmospheric entry. The method allows for both thermal and chemical nonequilibrium. Results are presented for a sixty-degree spherically blunt cone for various nose radii and altitudes.

Description: The implementation of two different types of turbulence models for a flow solver using the Chimera overset grid method is examined. Various turbulence model characteristics, such as length scale determination and transition modeling, are found to have a significant impact on the computed pressure distribution for a multielement airfoil case. No inherent problem is found with using either algebraic or one-equation turbulence models with an overset grid scheme, but simulation of turbulence for multiple-body or complex geometry flows is very difficult regardless of the gridding method. For complex geometry flowfields, modification of the Baldwin-Lomax turbulence model is necessary to select the appropriate length scale in wall-bounded regions. The overset grid approach presents no obstacle to use of a one- or two-equation turbulence model. Both Baldwin-Lomax and Baldwin-Barth models have problems providing accurate eddy viscosity levels for complex multiple-body flowfields such as those involving the Space Shuttle.

Description: A highly efficient implicit method for the computation of steady, two-dimensional compressible Navier-Stokes flowfields is presented. The discretization of the governing equations is hybrid in nature, with flux-vector splitting utilized in the streamwise direction and central differences with flux-limited artificial dissipation used for the transverse fluxes. Line Jacobi relaxation is used to provide a suitable initial guess for a new nonlinear iteration strategy based on line Gauss-Seidel sweeps. The applicability of quasi-Newton methods as convergence accelerators for this and other line relaxation algorithms is discussed, and efficient implementations of such techniques are presented. Convergence histories and comparisons with experimental data are presented for supersonic flow over a flat plate and for several high-speed compression corner interactions. Results indicate a marked improvement in computational efficiency over more conventional upwind relaxation strategies, particularly for flowfields containing large pockets of streamwise subsonic flow.

Description: A computational and experimental study of generic 3D sidewall compression inlets is conducted to examine the effects of fore and aft leading edge sweep on the internal shock structure. Inlets with leading edge sweeps of +30 deg and -30 deg with sidewall compression angles of 6 deg were tested in the NASA Langley Mach 4 air tunnel at a geometric contraction ratio of 1.87. The principal difference in performance was determined to be in the mass capture. Spillage was identified as having two components: a pressure induced component and a sweep induced component. It was found that while the direction of the leading edge sweep had a large influence on the spillage, the pressure effects were more important.

Description: A computational investigation has been conducted to study the effect and mechanisms of the passive control of a supersonic flow over a rectangular two-dimensional cavity. The passive control was included through the use of a porous surface over a vent chamber in the floor of the cavity. The passive control effectively suppressed the low-frequency pressure oscillations for the open type cavity, (length-to-depth ratio = 6.0). The mechanism for the suppression was observed to be the stabilization of the motion of the free shear layer. For the closed type cavity flow, (length-to-depth ratio = 17.5), the passive control modified the flowfield to nearly that of an open type cavity flow; further the cavity drag was reduced by a factor of four. The computational results of both cases showed good agreement with the available experimental data and the predictions of a semiempirical formula. This study demonstrates that the passive control concept can be used to improve the aerodynamic characteristics of open and closed cavity flowfields.

Description: A combined experimental/computational study has been performed of sharp fin induced shock wave/turbulent boundary layer interactions at Mach 5. The current paper focuses on the experiments and analysis of the results. The experimental data include mean surface heat transfer, mean surface pressure distributions and surface flow visualization for fin angles of attack of 6, 8, 10, 12, 14 and 16-degrees at Mach 5 under a moderately cooled wall condition. Comparisons between the results and correlations developed earlier show that Scuderi's correlation for the upstream influence angle (recast in a conical form) is superior to other such correlations in predicting the current results, that normal Mach number based correlations for peak pressure heat transfer are adequate and that the initial heat transfer peak can be predicted using pressure-interaction theory.

Description: The Direct Simulation Monte Carlo method is applied to a rarefied, weakly ionized, hypersonic flow over a blunt axisymmetric body. An ionization model based on the concept of ambipolar diffusion is used and a model for the sheath is presented. The effects of the new modeling techniques are investigated for flow over the Project Fire II configuration at 11.37 km/s at an altitude of 84.6 km. The calculated results are presented and compared with both experimental data and solutions where ionization effects were not included. In general, the calculated results overpredict the experimental values by about 15-20 percent.

Description: The efforts toward realistic vortex modeling for rotary wings which began under the guidance of professor A. A. Nikolsky of Princeton University in 1955-1956 are discussed. Attention is given to Nikolsky's flow-visualization studies and major theoretical considerations for vortex modeling. More recent efforts by other researchers have led to models of increasing complexity. The neglect of compressibility and viscous effects in the classical approach is noted to be a major limiting factor in full-scale rotor applications of the classical vortex theory; it has nevertheless been valuable for the delineation of problem areas and the guiding of both experimental and theoretical investigations.

Description: The development status of tactical missile dynamic behavior prediction is evaluated, with a view to the formulation of guidelines for high angle-of-attack tests where highly nonlinear aerodynamics are encountered and strong cross-coupling arises between lateral and longitudinal degrees of freedom. Attention is given to the Ericsson (1979) inviscid computer code, which is entirely general within the geometric constraints of axisymmetry and low aspect ratio lifting surfaces. A serious problem in the prediction of full-scale high-alpha dynamics is the almost total dependence on subscale experimental data, whose use is severely restricted by Reynolds number scaling as well as by model support and wind tunnel wall interference.

Description: The papers presented in this volume provide an overview of current theoretical and experimental research in the field of hypersonic waveriders. In particular, attention is given to efficient waveriders from known axisymmetric flow fields, hypersonic waverider design from given shock waves, limitations of waveriders, and aerodynamic stability theory of hypersonic waveriders. The discussion also covers momentum analysis of waverider flow fields, tethered aerothermodynamic research for hypersonic waveriders, simulation of hypersonic waveriders, and an idealized tip-to-tail waverider model.

Description: Forward swept sidewall compression inlets have been tested in the Mach 4 Blowdown Facility at the NASA Langley Research Center to study the effects of bodyside compression surfaces on inlet performance in the presence of an incoming turbulent boundary layer. The measurements include mass flow capture and mean surface pressure distributions obtained during simulated combustion pressure increases downstream of the inlet. The kerosene-lampblack surface tracer technique has been used to obtain patterns of the local wall shear stress direction. Inlet performance is evaluated using starting and unstarting characteristics, mass capture, mean surface pressure distributions and permissible back pressure limits. The results indicate that inlet performance can be improved with selected bodyside compression surfaces placed between the inlet sidewalls.

Description: Unsteady flowfields around oscillating Boeing VR7 airfoil with and without a leading-edge slat were numerically investigated by a novel zonal method using a conformal mapping technique. Numerical aero-dynamic hysteresis loops show that the leading-edge slat prevents the airfoil dynamic stall at reduced frequency of 0.15, Reynolds number of 1 million, and the oscillation range of 5 deg to 25 deg.

Description: Experiments on supersonic jet control via point disturbances inside the nozzle were conducted using a pressure-matched converging-diverging nozzle with an exit diameter of 3.11 cm and an exit Mach number of 2.1. The disturbance generators are situated on the supersonic nozzle wall between the throat and the jet exit plane; the resulting fluid-mechanical disturbances are observed at the intersection of the disturbance cone and the jet lip. It is demonstrated that a single disturbance generator leads to a polygonal jet appearance, while two symmetric disturbances lead to a symmetric rectangular jet shape and two nonsymmetric disturbances lead to a star-like appearance. Significant cusps are formed in all cases.

Description: Approximate inviscid and boundary layer techniques for aerodynamic heating calculations are discussed. An inviscid flowfield solution is needed to provide surface pressures and boundary-layer edge properties. Modified Newtonian pressures coupled with an approximate shock shape will suffice for relatively simple shapes like sphere-cones with cone half-angles between 15 and 45 deg. More accurate approximate methods have been developed which make use of modified Maslen techniques. Slender and large angle sphere-cones and more complex shapes generally require an Euler code, like HALIS, to provide that information. The boundary-layer solution is reduced significantly by using the axisymmetric analog and approximate heating relations developed by Zoby, et al. (1981). Analysis is presented for the calculation of inviscid surface streamlines and metrics. Entropy-layer swallowing effects require coupling the inviscid and boundary-layer solutions.

Description: A combined experimental and numerical investigation of strut/endwall interactions within an annular duct having a supersonic core flow has been conducted. Four diamond-shaped struts with a 7 deg half angle were positioned circumferentially equidistant within an annular duct having a gap height of 0.7 strut chords, and an inner-to-outer wall radius ratio of 0.7. Turbulent boundary layers exist on both inner and outer walls of the duct, but have not merged. The core flow upstream of the struts is uniform at a nominal Mach number of 3.0 and a Reynolds number of 3 x 10 exp 5 based on the strut chord length. Experimental results, which include Pitot pressure distributions within the flow field, static pressure distributions on the inner and outer walls of the duct, and oil flow visualization on the centerbody and strut, are presented and compared with CFD predictions. Secondary flows associated with the interactions are examined including the trajectories of the horseshoe vortices formed at the leading and trailing edges of the strut and the trajectories of the vortices formed in the corner of the strut/endwall intersection.

Description: The hypersonic flow over a cavity is investigated. The time-dependent compressible Navier-Stokes equations, in terms of mass averaged variables, are numerically solved. An implicit algorithm, with a subiteration procedure to recover time-accuracy, is used to perform the time-accurate computations. The objective of the study is to investigate the effects of Reynolds number and cavity dimensions. The comparison of the computations with available experimental data, in terms of time mean static pressure, heat transfer, and Mach number show good agreement. In the computations large vortex structures, which adversely affect the cavity flow characteristics, are observed at the rear of the cavity. A self-sustained oscillatory motion occurs within the cavity over a range of Reynolds number and cavity dimensions. The frequency spectra of the oscillations show good agreement with a modified semi-empirical relation.

Description: The boundary layer on the wall of the Hypervelocity Tunnel 9 was investigated with pitot pressure and total temperature measurements. Experimental results are presented for standard and supercooled Mach 14 runs. The boundary layer data at supercooled conditions are compared to numerical predictions made with a Navier-Stokes algorithm including vibrational nonequilibrium and intermolecular force effects. For standard tunnel conditions, the numerical solutions agree well with experimental data. For the supercooled cases, the numerical code predicts the total temperature but overpredicts the pitot pressure.

Description: This paper describes a practical method for the design of an infinite cascade in incompressible flow. The method is based on conformal mapping, and as a result it allows for multi-point design. The cascade blade to be determined is divided into a desired number of segments. Over each segment, the velocity distribution is prescribed together with an inlet or outlet flow angle at which this velocity distribution is to be achieved. In this way multi-point design requirements can be met. It is necessary to satisfy several conditions that arise to guarantee compatibility with the inlet and outlet flow as well as closure of the cascade blade. Satisfaction of these conditions does not necessarily result in a cascade with all the desired characteristics. For example, the cascade blades may be bulbous or crossed. Through Newton iteration, however, the desired characteristics may be prescribed by allowing for the adjustment of the design parameters that define the mathematical problem through conformal mapping. Several examples will be illustrated to demonstrate the capability of the method. It will be shown that the method is limited to the design of cascades with solidities of up to one.

Description: The present study represents an extension of an earlier wind tunnel experiment performed with the P&W 17-in. Advanced Ducted Propeller (ADP) Simulator operating at Mach 0.2. In order to study the effects of a rotating propeller on the inlet flow, data were obtained in the UTRC 10- by 15-Foot Large Subsonic Wind Tunnel with the same hardware and instrumentation, but with the propellar removed. These new tests were performed over a range of flow rates which duplicated flow rates in the powered simulator program. The flow through the inlet was provided by a remotely located vacuum source. A comparison of the results of this flow-through study with the previous data from the powered simulator indicated that in the conventional inlet the propeller produced an increase in the separation angle of attack between 4.0 deg at a specific flow of 22.4 lb/sec-sq ft to 2.7 deg at a higher specific flow of 33.8 lb/sec-sq ft. A similar effect on separation angle of attack was obtained by using stationary blockage rather than a propeller.

Description: Thicknesses of residual ice are presented to provide information on surface contamination and associated roughness during deicing events. Data was obtained from low power ice protection systems tests conducted in the Icing Research Tunnel at NASA Lewis Research Center (LeRC) with nine different deicing systems. Results show that roughness associated with residual ice is not characterized by uniformly distributed roughness. Results also show that deicing systems require a critical mass of ice to generate a sufficient expelling force to remove the ice.

Description: Spatial adaptation procedures for the accurate and efficient solution of steady and unsteady inviscid flow problems are described. The adaptation procedures were developed and implemented within a three-dimensional, unstructured-grid, upwind-type Euler code. These procedures involve mesh enrichment and mesh coarsening to either add points in high gradient regions of the flow or remove points where they are not needed, respectively, to produce solutions of high spatial accuracy at minimal computational cost. The paper gives a detailed description of the enrichment and coarsening procedures and presents comparisons with experimental data for an ONERA M6 wing and an exact solution for a shock-tube problem to provide an assessment of the accuracy and efficiency of the capability. Steady and unsteady results, obtained using spatial adaptation procedures, are shown to be of high spatial accuracy, primarily in that discontinuities such as shock waves are captured very sharply.

Description: Studies of the effects of variable blowing in turbulent hypersonic boundary layers are presented. Numerical calculations of the skin friction and surface heat transfer rates are compared to the experimental measurements of Holden (1990) for a slender cone at zero angle of attack in steady flows at Mach numbers of 11 and 13. An analysis of the transpiration feed system of the cone model was performed and showed that the blowing rate could be variable along the cone surface. This effect is confirmed by internal pressure measurements which were taken inside the cone model. The blowing rates are recalibrated using the internal gauge readings and used as the wall boundary condition for a compressible turbulent boundary layer calculation using the low Reynolds number k-epsilon model of Chien (1982). At low blowing rates, the boundary layer calculations indicate that a situation where both the effects of suction and blowing are present within the same flow. The results show excellent qualitative prediction of the experimental data.

Description: This experimental study examines the effects of modifications to the incoming turbulent boundary layer on the highly separated shock wave/boundary layer interaction generated by an unswept compression corner. Particular focus is placed on the motion of the unsteady separation shock wave. The flowfield was generated by a 28 deg ramp in a Mach 5 flow with a freestream Reynolds number of 50 x 10 exp 6/m. The incoming turbulent boundary layer transitioned naturally and developed under near-adiabatic wall conditions. Modification of the flow entering the interaction was effected through either a single plate boundary layer manipulator (BLM) or riblets. The BLM reduced the length of separation by 35-45 percent and reduced the streamwise extent of the separation shock motion by 36-74 percent. Examination of the flowfield downstream of the BLM showed this result to be due to the inviscid preturning of the flow by the BLM, and not by changes to the boundary layer dynamics. The riblets had no measurable effect on the compression corner interaction.

Description: Continued advances in the various elements that comprise the field of computational fluid dynamics (CFD) are promoting a radically different approach to the aerodynamic design and analysis of aerospace vehicles and systems. The elements of CFD generally include numerical algorithm development, transition and turbulence modeling, surface modeling, and grid generation, scientific visualization and validation methodologies. This paper discusses the research progress and prospects for the future in each of these elements within NASA's CFD and Experimental Validation Program. The applicability of computational methods for the purposes of understanding complex flow phenomena, exploring aerodynamic concepts, and providing vehicle-design input is also addressed.

Description: The performance of a funnel-shaped hypersonic vehicle forebody is shown at its design Mach number of 15 and at two lower off-design Mach numbers. Inlet flow quality is presented by contour plots of flow angularity and pressure. Performance is indicated by mass capture ratio, total pressure recovery, kinetic energy efficiency and drag.

Description: Precise, smooth aerodynamic models are required for implementing adaptive, nonlinear control strategies. Accurate representations of aerodynamic coefficients can be generated for the complete flight envelope by combining computational neural network models with an Estimation-Before-Modeling paradigm for on-line training information. A novel method of incorporating first-partial-derivative information is employed to estimate the weights in individual feedforward neural networks for each aerodynamic coefficient. The method is demonstrated by generating a model of the normal force coefficient of a twin-jet transport aircraft from simulated flight data, and promising results are obtained.

Description: Aerothermodynamic, aerodynamic, and atmospheric science data acquired between 55 and 150 km has been limited by the lack of vehicles or platforms capable of sustained operation at these altitudes. Tethered satellites, which have been under study for this purpose by NASA, the Italian Space Agency (ASI), and others for more than a decade, are expected to become a reality by mid-1991. This approach, in which an instrumented platform is maintained at a desired altitude by a tether attached to a host vehicle orbiting at higher altitudes, will provide the first opportunity to obtain steady state data over an extended period encompassing one or more orbital revolutions. This paper describes the objectives and measurement methods for the first of the facility-class satellites, the TSS-2, which is proposed for a 1995 deployment, and gives the status of the experiment definition. Monte Carlo modeling of the flow fields at 130 km around the baseline 1.6 m diameter sphere is discussed and illustrative results of the modeling given.

Description: A computational fluid dynamics code which utilizes both structured and unstructured grids was developed. The objective of this study was to develop and demonstrate the ability of such a code to achieve solutions about complex geometries in two dimensions. An unstructured grid generator and flow solver were incorporated into the PARC2D structured flow solver. This new unstructured grid generator capability allows for easier generation and manipulation of complex grids. Several examples of the grid generation capabilities are provided. The coupling of different grid topologies and the manipulation of individual grids is shown. Also, grids for realistic geometries, a NACA 0012 airfoil and a wing/nacelle installation, were created. The flow over a NACA 0012 airfoil was used as a test case for the flow solver. Eight separate cases were run. They were both the inviscid and viscous solutions for two freestream Mach numbers and airfoil angle of attacks of 0 to 3.86 degrees. The Mach numbers chosen were for a subsonic case, Mach 0.6, and a case where supersonic regions and a shock wave exists, Mach 0.8. These test case conditions were selected to match experimentally obtained data for code comparison. The results show that the code accurately predicts the flow field for all cases.

Description: An approach is presented for the generation of two-dimensional, structured, dynamic grids. The grid motion may be due to the motion of the boundaries of the computational domain or to the adaptation of the grid to the transient, physical solution. A time-dependent grid is computed through the time integration of the grid speeds which are computed from a system of grid speed equations. The grid speed equations are derived from the time-differentiation of the grid equations so as to ensure that the dynamic grid maintains the desired qualities of the static grid. The grid equations are the Euler-Lagrange equations derived from a variational statement for the grid. The dynamic grid method is demonstrated for a model problem involving boundary motion, an inviscid flow in a converging-diverging nozzle during startup, and a viscous flow over a flat plate with an impinging shock wave. It is shown that the approach is more accurate for transient flows than an approach in which the grid speeds are computed using a finite difference with respect to time of the grid. However, the approach requires significantly more computational effort.

Description: A Cartesian, cell-based scheme for solving the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, polygonal 'cut' cells are created. The geometry of the cut cells is computed using polygon-clipping algorithms. The grid is stored in a binary-tree data structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded, with a limited linear reconstruction of the primitive variables used to provide input states to an approximate Riemann solver for computing the fluxes between neighboring cells. A multi-stage time-stepping scheme is used to reach a steady-state solution. Validation of the Euler solver with benchmark numerical and exact solutions is presented. An assessment of the accuracy of the approach is made by uniform and adaptive grid refinements for a steady, transonic, exact solution to the Euler equations. The error of the approach is directly compared to a structured solver formulation. A non smooth flow is also assessed for grid convergence, comparing uniform and adaptively refined results. Several formulations of the viscous terms are assessed analytically, both for accuracy and positivity. The two best formulations are used to compute adaptively refined solutions of the Navier-Stokes equations. These solutions are compared to each other, to experimental results and/or theory for a series of low and moderate Reynolds numbers flow fields. The most suitable viscous discretization is demonstrated for geometrically-complicated internal flows. For flows at high Reynolds numbers, both an altered grid-generation procedure and a different formulation of the viscous terms are shown to be necessary. A hybrid Cartesian/body-fitted grid generation approach is demonstrated. In addition, a grid-generation procedure based on body-aligned cell cutting coupled with a viscous stensil-construction procedure based on quadratic programming is presented.

Description: The cost and time to certify or qualify a rotorcraft for flight in forecast icing has been a major impediment to the development of ice protection systems for helicopter rotors. Development and flight test programs for those aircraft that have achieved certification or qualification for flight in icing conditions have taken many years, and the costs have been very high. NASA, Sikorsky, and others have been conducting research into alternative means for providing information for the development of ice protection systems, and subsequent flight testing to substantiate the air-worthiness of a rotor ice protection system. Model rotor icing tests conducted in 1989 and 1993 have provided a data base for correlation of codes, and for the validation of wind tunnel icing test techniques. This paper summarizes this research, showing test and correlation trends as functions of cloud liquid water content, rotor lift, flight speed, and ambient temperature. Molds were made of several of the ice formations on the rotor blades. These molds were used to form simulated ice on the rotor blades, and the blades were then tested in a wind tunnel to determine flight performance characteristics. These simulated-ice rotor performance tests are discussed in the paper. The levels of correlation achieved and the role of these tools (codes and wind tunnel tests) in flight test planning, testing, and extension of flight data to the limits of the icing envelope are discussed. The potential application of simulated ice, the NASA LEWICE computer, the Sikorsky Generalized Rotor Performance aerodynamic computer code, and NASA Icing Research Tunnel rotor tests in a rotorcraft certification or qualification program are also discussed. The correlation of these computer codes with tunnel test data is presented, and a procedure or process to use these methods as part of a certification or qualification program is introduced.

Description: This paper examines a supersonic multi jet interaction problem that we believe is likely to be important for mixing enhancement and noise reduction in supersonic mixer-ejector nozzles. We demonstrate that it is possible to synchronize the screech instability of four rectangular jets by precisely adjusting the inter jet spacing. Our experimental data agrees with a theory that assumes that the phase-locking of adjacent jets occurs through a coupling at the jet lip. Although the synchronization does not change the frequency of the screech tone, its amplitude is augmented by 10 dB. The synchronized multi jets exhibit higher spreading than the unsynchronized jets, with the single jet spreading the least. We compare the nearfield noise of the four jets with synchronized screech to the noise of the sum of four jets operated individually. Our noise measurements reveal that the more rapid mixing of the synchronized multi jets causes the peak jet noise source to move up stream and to radiate noise at larger angles to the flow direction. Based on our results, we believe that screech synchronization is advantageous for noise reduction internal to a mixer-ejector nozzle, since the noise can now be suppressed by a shorter acoustically lined ejector.

Description: A full Navier-Stokes analysis was performed to evaluate the performance of the subsonic diffuser of a NASA Lewis Research Center 70/30 mixed-compression bifurcated supersonic inlet for high speed civil transport application. The PARC3D code was used in the present study. The computations were also performed when approximately 2.5 percent of the engine mass flow was allowed to bypass through the engine bypass doors. The computational results were compared with the available experimental data which consisted of detailed Mach number and total pressure distribution along the entire length of the subsonic diffuser. The total pressure recovery, flow distortion, and crossflow velocity at the engine face were also calculated. The computed surface ramp and cowl pressure distributions were compared with experiments. Overall, the computational results compared well with experimental data. The present CFD analysis demonstrated that the bypass flow improves the total pressure recovery and lessens flow distortions at the engine face.

Description: Although the plane wake is marked by the formation of strong spanwise vortices, the initially two-dimensional Karman-like vortices soon develop a three-dimensional structure in the form of secondary streamwise vortices. So far, this streamwise vortex structure has been studied mostly through flow visualization and at relatively low Reynolds numbers. The primary objective of the present program was to investigate the origin and evolution of the three-dimensional structure of straight and curved plane wakes at relatively high Reynolds numbers (Re(sub b) = 28,000) through detailed measurements of the mean and turbulent properties at several streamwise locations. The experiments were conducted in three phases. In the first phase, the development of a straight plane wake was investigated. In the second phase, the effects of imposed streamwise curvature on the wake development were examined. The streamwise curvature was of constant radius and very mild in terms of the curvature ratio (b/(square root of R) is less than 2 percent). In both the first and second phases, the role of initial conditions was examined in wakes generated from both untripped (laminar) and tripped (turbulent) initial boundany layers. In the third phase, the effects of injecting streamwise vorticity and the effects of increased Reynolds number on the tripped wake structure and development were investigated.

Description: A two-dimensional symmetric wedge configuration representative of a single high-speed intake in steady flow was investigated. The analysis presented here is intended as an engineering approach for estimating certain features of the internal shock system. The primary interest here is the prediction of the size and location of the almost-normal shock wave that develops when the leading-edge shocks intersect at angles above a certain critical value that is less than the wedge detachment angle. The almost-normal shock wave is frequently referred to as the 'Mach stem', Parametric studies enabled the sensitivity of the Mach stem height to various flowfield parameters to be examined, thus indicating how accurately these parameters must be measured in a given experiment. Results of these predictions were compared with those of a steady-flow experiment performed at nominal freestream Mach numbers from 2.8 to 5. The predicted stem heights were consistently lower than the mean experimental values, attributable both to experimental uncertainties and to certain simplifying assumptions used in the analysis. Modification of these assumptions to better represent the test environment improved the analytical results.

Description: An engineering method has been developed that couples an approximate three dimensional inviscid technique with the axisymmetric analog and a set of approximate convective heating equations. The displacement effect on the boundary layer on the outer inviscid flow is calculated and included as a boundary condition in the inviscid technique. This accounts for the viscous interaction present at lower Reynolds numbers. The method is applied to blunted axisymmetric and three dimensional elliptic cones at angle of attack for the laminar hypersonic flow of a perfect gas. The method is applied to turbulent and equilibrium-air conditions. The present technique predicts surface heating rates, pressures, and shock shapes that compare favorably with experimental (ground-test and flight) data and numerical solutions of the Navier-Stokes and viscous shock-layer equations. In addition, the inclusion of viscous interaction significantly improves results obtained at lower Reynolds numbers. The new technique represents a major improvement over current engineering aerothermal methods with only a modest increase in computational effort.

Description: A device for controlling attachment line contamination on an airfoil is presented. A fence is installed on the leading edge of the airfoil in the freestream direction perpendicular to the airfoil, outboard of the fuselage boundary layer. The inboard side of the fence arrests the spanwise movement of the turbulent boundary layer while the laminar boundary layer on the outboard side of the fence eliminates any further turbulent contamination of the attachment line.

Description: A general method is described for automatically discretizing, into unstructured assemblies of tetrahedra, the three-dimensional solution domains of complex shape which are of interest in practical computational aerodynamics. An algorithm for the solution of the compressible Euler equations which can be implemented on such general unstructured tetrahedral grids is described. This is an explicit cell-vertex scheme which follows a general Taylor-Galerkin philosophy. The approach is employed to compute a transonic inviscid flow over a standard wing and the results are shown to compare favorably with experimental observations. As a more practical demonstration, the method is then applied to the analysis of inviscid flow over a complete modern fighter configuration. The effect of using mesh adaptivity is illustrated when the method is applied to the solution of high speed flow in an engine inlet.

Description: In order to have confidence in a numerical method, the verification of its reproduction of known benchmark analytic solutions for simple model problems is of great importance. Attention is presently given to a novel benchmarking procedure for numerical models of high speed, reactive 2D flows. The procedure is illustrated by comparing asymptotic and numerical solutions for oblique detonations in which an attached oblique shock is followed by an exothermic reaction with a thick reaction zone.

Description: This paper presents an experimental database selected and compiled from aerothermal measurements obtained on basic model configurations on which fundamental flow phenomena could be most easily examined. The experimental studies were conducted in hypersonic flows in 48-inch, 96-inch, and 6-foot shock tunnels. A special computer program was constructed to provide easy access to the measurements in the database as well as the means to plot the measurements and compare them with imported data. The database contains tabulations of model configurations, freestream conditions, and measurements of heat transfer, pressure, and skin friction for each of the studies selected for inclusion. The first segment contains measurements in laminar flow emphasizing shock-wave boundary-layer interaction. In the second segment, measurements in transitional flows over flat plates and cones are given. The third segment comprises measurements in regions of shock-wave/turbulent-boundary-layer interactions. Studies of the effects of surface roughness of nosetips and conical afterbodies are presented in the fourth segment of the database. Detailed measurements in regions of shock/shock boundary layer interaction are contained in the fifth segment. Measurements in regions of wall jet and transpiration cooling are presented in the final two segments.