ALBERT

Description: The development of shock-capturing finite difference methods for hyperbolic conservation laws has been a rapidly growing area for the last decade. Many of the fundamental concepts, state-of-the-art developments and applications to fluid dynamics problems can only be found in meeting proceedings, scientific journals and internal reports. This paper attempts to give a unified and generalized formulation of a class of high-resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock waves, perfect gases, equilibrium real gases and nonequilibrium flow computations. These numerical methods are formulated for the purpose of ease and efficient implementation into a practical computer code. The various constructions of high-resolution shock-capturing methods fall nicely into the present framework and a computer code can be implemented with the various methods as separate modules. Included is a systematic overview of the basic design principle of the various related numerical methods. Special emphasis will be on the construction of the basic nonlinear, spatially second and third-order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows will be discussed. Some perbolic conservation laws to problems containing stiff source terms and terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas-dynamics problems. The use of the Lax-Friedrichs numerical flux to obtain high-resolution shock-capturing schemes is generalized. This method can be extended to nonlinear systems of equations without the use of Riemann solvers or flux-vector splitting approaches and thus provides a large savings for multidimensional, equilibrium real gases and nonequilibrium flow computations.

Description: One of the major achievements in engineering science has been the development of computer algorithms for solving nonlinear differential equations such as the Navier-Stokes equations. These algorithms are now used in the practical engineering design of devices such as cars and airplanes as well as theoretical studies of complex phenomena such as fluid turbulence. In past years, limited computer resources have motivated the development of efficient numerical methods in computational fluid dynamics (CFD) utilizing structured meshes. These meshes are comprised of systematic arrays of quadrilateral or hexahedral cells. The use of structured meshes greatly simplifies the implementation of CFD algorithms on conventional computers. Structured meshes also permit the use of highly efficient solution techniques such as alternating direction implicit (ADI) iteration schemes or multigrid. Following the dramatic improvement in computing speed in recent years, emphasis has shifted towards the design of algorithms capable of treating complex geometries. The automatic generation of structured grids about complex geometries is problematic. Unstructured grids offer one promising alternative technique for treating these general geometries. Unstructured meshes have irregular connectivity and usually contain triangles and/or quadrilaterals in two dimensions and tetrahedra and/or hexahedra in three dimensions. The generation and use of unstructured grids poses new challenges in computational fluid dynamics. This is true for both grid generation as well as for the design of algorithms for flow solution. The purpose of these notes is to present recent developments in the unstructured grid generation and flow solution technology.

Description: One of the most important elements in laser velocimetry, yet the most neglected, is the small particle embedded in the flow field that scatters the light necessary to make velocity measurements. The characteristics of this small particle are often ignored in the effort to obtain data. This seems strange since it is the primary cause of measurement error. If the particle is too large, it will not follow the flow resulting in an inaccurate representation of the fluid velocity. If the particle is too small, it will not scatter sufficient light to provide the signal-to-noise necessary to minimize measurement uncertainty in the signal processing electronics. This lecture will attempt to remove the confusion in choosing a seeding method by assessing many of the techniques currently used. It will outline their characteristics and typical limitations imposed by various applications. The lecture will then focus on the ramifications of these methods on measurement accuracy.

Description: As a tool for large-scale computations in fluid dynamics, spectral methods were prophesized in 1944, born in 1954, virtually buried in the mid-1960's, resurrected in 1969, evangalized in the 1970's, and catholicized in the 1980's. The use of spectral methods for meteorological problems was proposed by Blinova in 1944 and the first numerical computations were conducted by Silberman (1954). By the early 1960's computers had achieved sufficient power to permit calculations with hundreds of degrees of freedom. For problems of this size the traditional way of computing the nonlinear terms in spectral methods was expensive compared with finite-difference methods. Consequently, spectral methods fell out of favor. The expense of computing nonlinear terms remained a severe drawback until Orszag (1969) and Eliasen, Machenauer, and Rasmussen (1970) developed the transform methods that still form the backbone of many large-scale spectral computations. The original proselytes of spectral methods were meteorologists involved in global weather modeling and fluid dynamicists investigating isotropic turbulence. The converts who were inspired by the successes of these pioneers remained, for the most part, confined to these and closely related fields throughout the 1970's. During that decade spectral methods appeared to be well-suited only for problems governed by ordinary diSerential eqllations or by partial differential equations with periodic boundary conditions. And, of course, the solution itself needed to be smooth. Some of the obstacles to wider application of spectral methods were: (1) poor resolution of discontinuous solutions; (2) inefficient implementation of implicit methods; and (3) drastic geometric constraints. All of these barriers have undergone some erosion during the 1980's, particularly the latter two. As a result, the applicability and appeal of spectral methods for computational fluid dynamics has broadened considerably. The motivation for the use of spectral methods in numerical calculations stems from the attractive approximation properties of orthogonal polynomial expansions.

Description: The probes described herein, in various configurations, permit the measurement of the volume fraction of two or more fluids flowing through a pipe. Each probe measures the instantaneous relative dielectric constant of the fluid in immediate proximity. As long as separation of the relative dielectric constants of each fluid is possible, several or even many fluids can be measured in the same flow steam. By using multiple probes, the velocity of each fluid can generally be determined as well as the distribution of each constituent in the pipe. The values are determined by statistical computation. There are many potential applications for probes of this type in industry and government. Possible NASA applications include measurements of helium/hydrazine flow during rocket tests at White Sands, liquid/gas flow in hydrogen or oxygen lines in Orbiter engines, and liquid/gaseous Freon flow in zero gravity tests with the KS135 aircraft at JSC. Much interest has been shown recently by the oil industry. In this industry, a good method is needed to measure the fractions of oil, water, and natural gas flowing in a pipeline and the velocity of each. This particular problem involves an extension of what has been developed to date and our plans and program to solve this problem will be discussed herein.

Description: Numerical simulations of incompressible turbulent flow through rectangular ducts with one side rib-roughened were performed to determine pressure drop. The "PHOENICS " software package was used for the computations, which required provision of a wall function for transverse rib-roughened surfaces. The present study was conducted in the range of 10(exp 5) less than or equal to Reynolds number less than or equal to 10(exp 7), 0.01 less than or equal to rib height to hydraulic diameter ratio less than or equal to 0.04, 10 less than or equal to pitch to rib height ratio less than or equal to 40. Using the numerical results, friction factor charts for various aspect ratios were generated. The numerical results agreed well with experimental data that was obtained for 10(exp 5) less than Reynolds less than 2 x 10(exp 5). In addition, a scheme for predicting friction factor using existing correlations for smooth and rough walls was developed.

Description: A NASA report detailing a wind tunnel investigation of a variable camber and twist could effectively reduce drag, thus improving performance. The resulting VooDoo fin is made of composite materials, has a rigid internal spar and a flexible polymer exterior coating. It is computer-designed and exceptionally durable.

Description: The purpose of this work was the evaluation of the use of electron-bean fluorescence for flow measurements during hypersonic flight. Both analytical and numerical models were developed in this investigation to evaluate quantitatively flow field imaging concepts based upon the electron beam fluorescence technique for use in flight research and wind tunnel applications. Specific models were developed for: (1) fluorescence excitation/emission for nitrogen, (2) rotational fluorescence spectrum for nitrogen, (3) single and multiple scattering of electrons in a variable density medium, (4) spatial and spectral distribution of fluorescence, (5) measurement of rotational temperature and density, (6) optical filter design for fluorescence imaging, and (7) temperature accuracy and signal acquisition time requirements. Application of these models to a typical hypersonic wind tunnel flow is presented. In particular, the capability of simulating the fluorescence resulting from electron impact ionization in a variable density nitrogen or air flow provides the capability to evaluate the design of imaging instruments for flow field mapping. The result of this analysis is a recommendation that quantitative measurements of hypersonic flow fields using electron-bean fluorescence is a tractable method with electron beam energies of 100 keV. With lower electron energies, electron scattering increases with significant beam divergence which makes quantitative imaging difficult. The potential application of the analytical and numerical models developed in this work is in the design of a flow field imaging instrument for use in hypersonic wind tunnels or onboard a flight research vehicle.

Description: Any aircraft preliminary design study requires a structural model of the proposed configuration. The model must be capable of estimating the structural weight of a given configuration, and of predicting the deflections which will result from foreseen flight and ground loads. The present work develops such a model for the proposed Oblique All Wing airplane. The model is based on preliminary structural work done by Jack Williams and Peter Rudolph at Mdng, and is encoded in a FORTRAN program. As a stand-alone application, the program can calculate the weight CG location, and several types of structural deflections; used in conjunction with an aerodynamics model, the program can be used for mission analysis or sizing studies.

Description: The purpose of this investigation is to provide a comprehensive data base for the validation of numerical simulations. The objective of the present paper is to provide a tabulation of the experimental data. The data were obtained in the two-dimensional, transonic flowfield surrounding a supercritical airfoil. A variety of flows were studied in which the boundary layer at the trailing edge of the model was either attached or separated. Unsteady flows were avoided by controlling the Mach number and angle of attack. Surface pressures were measured on both the model and wind tunnel walls, and the flowfield surrounding the model was documented using a laser Doppler velocimeter (LDV). Although wall interference could not be completely eliminated, its effect was minimized by employing the following techniques. Sidewall boundary layers were reduced by aspiration, and upper and lower walls were contoured to accommodate the flow around the model and the boundary-layer growth on the tunnel walls. A data base with minimal interference from a tunnel with solid walls provides an ideal basis for evaluating the development of codes for the transonic speed range because the codes can include the wall boundary conditions more precisely than interference connections can be made to the data sets.

Description: Hybrid grids, composed of structured and unstructured grids, combines the best features of both. The chimera method is a major stepstone toward a hybrid grid from which the present approach is evolved. The chimera grid composes a set of overlapped structured grids which are independently generated and body-fitted, yielding a high quality grid readily accessible for efficient solution schemes. The chimera method has been shown to be efficient to generate a grid about complex geometries and has been demonstrated to deliver accurate aerodynamic prediction of complex flows. While its geometrical flexibility is attractive, interpolation of data in the overlapped regions - which in today's practice in 3D is done in a nonconservative fashion, is not. In the present paper we propose a hybrid grid scheme that maximizes the advantages of the chimera scheme and adapts the strengths of the unstructured grid while at the same time keeps its weaknesses minimal. Like the chimera method, we first divide up the physical domain by a set of structured body-fitted grids which are separately generated and overlaid throughout a complex configuration. To eliminate any pure data manipulation which does not necessarily follow governing equations, we use non-structured grids only to directly replace the region of the arbitrarily overlapped grids. This new adaptation to the chimera thinking is coined the DRAGON grid. The nonstructured grid region sandwiched between the structured grids is limited in size, resulting in only a small increase in memory and computational effort. The DRAGON method has three important advantages: (1) preserving strengths of the chimera grid; (2) eliminating difficulties sometimes encountered in the chimera scheme, such as the orphan points and bad quality of interpolation stencils; and (3) making grid communication in a fully conservative and consistent manner insofar as the governing equations are concerned. To demonstrate its use, the governing equations are discretized using the newly proposed flux scheme, AUSM+, which will be briefly described herein. Numerical tests on representative 2D inviscid flows are given for demonstration. Finally, extension to 3D is underway, only paced by the availability of the 3D unstructured grid generator.

Description: A three-dimensional computational fluid dynamics code, RPLUS3D, which was developed for the reactive propulsive flows of ramjets and scramjets, was validated for glancing shock wave-boundary layer interactions. Both laminar and turbulent flows were studied. A supersonic flow over a wedge mounted on a flat plate was numerically simulated. For the laminar case, the static pressure distribution, velocity vectors, and particle traces on the flat plate were obtained. For turbulent flow, both the Baldwin-Lomax and Chien two-equation turbulent models were used. The static pressure distributions, pitot pressure, and yaw angle profiles were computed. In addition, the velocity vectors and particle traces on the flat plate were also obtained from the computed solution. Overall, the computed results for both laminar and turbulent cases compared very well with the experimentally obtained data.

Description: An experimental investigation of the aerodynamic characteristics of thin, moderately swept fighter wings has been conducted to evaluate the effect of camber and twist on the effectiveness of leading- and trailing-edge flaps at supersonic speeds in the Langley Unitary Plan Wind Tunnel. The study geometry consisted of a generic fuselage with camber typical of advanced fighter designs without inlets, canopy, or vertical tail. The model was tested with two wing configurations an uncambered (flat) wing and a cambered and twisted wing. Each wing had an identical clipped delta planform with an inboard leading edge swept back 65 deg and an outboard leading edge swept back 50 deg. The trailing edge was swept forward 25 deg. The leading-edge flaps were deflected 4 deg to 15 deg, and the trailing-edge flaps were deflected from -30 deg to 10 deg. Longitudinal force and moment data were obtained at Mach numbers of 1.60, 1.80, 2.00, and 2.16 for an angle-of-attack range 4 deg to 20 deg at a Reynolds number of 2.16 x 10(exp 6) per foot and for an angle-of-attack range 4 deg to 20 deg at a Reynolds number of 2.0 x 10(exp 6) per foot. Vapor screen, tuft, and oil flow visualization data are also included.

Description: The NASA Langley 8-Foot Transonic Pressure Tunnel is a continuous-flow, variable-pressure wind tunnel with control capability to independently vary Mach number, stagnation pressure, stagnation temperature, and humidity. The top and bottom walls of the test section are axially slotted to permit continuous variation of the test section Mach number from 0.2 to 1.2, the slot-width contour provides a gradient-free test section 50 in. long for Mach numbers equal to or greater than 1.0 and 100 in. long for Mach numbers less than 1.0. The stagnation pressure may be varied from 0.25 to 2.0 atm. The tunnel test section has been recalibrated to determine the relationship between the free-stream Mach number and the test chamber reference Mach number. The hardware was the same as that of an earlier calibration in 1972 but the pressure measurement instrumentation available for the recalibration was about an order of magnitude more precise. The principal result of the recalibration was a slightly different schedule of reentry flap settings for Mach numbers from 0.80 to 1.05 than that determined during the 1972 calibration. Detailed tunnel contraction geometry, test section geometry, and limited test section wall boundary layer data are presented.

Description: A numerical investigation is carried out to determine the magnitude of wake radiation for a proposed Venus composition probe. One of the scientific goals of the mission is to determine the atmospheric composition of Venus by examining the intensity of scattered sunlight through the wake of the vehicle during planetary entry. In the wake of the vehicle, excited particles generated in the bow shock and boundary layers absorb and emit radiation. Thus, the purpose of this study is to determine if the radiation sensor will be able to sense the incoming solar radiative flux relative to the radiative flux generated in the wake. During portions of the entry trajectory the incident surface heat flux will be high enough to produce significant ablation. Ablation products such as CN are known to be strong radiators. Also, the ablation will be driven by strong radiation emanating from the bow shock. Thus, radiation and ablation will be coupled into the Navier-Stokes flow solutions.

Description: The ability to control the extent of laminar flow on swept wings at supersonic speeds may be a critical element in developing the enabling technology for a High Speed Civil Transport (HSCT). Laminar boundary layers are less resistive to forward flight than their turbulent counterparts, thus the farther downstream that transition from laminar to turbulent flow in the wing boundary layer is extended can be of significant economic impact. Due to the complex processes involved experimental studies of boundary layer stability and transition are needed, and these are performed in "quiet" wind tunnels capable of simulating the low-disturbance environment of free flight. At Ames, a wind tunnel has been built to operate at flow conditions which match those of the HSCT laminar flow flight demonstration 'aircraft, the F-16XL, i.e. at a Mach number of 1.6 and a Reynolds number range of 1 to 3 million per foot. This will allow detailed studies of the attachment line and crossflow on the leading edge area of the highly swept wing. Also, use of suction as a means of control of transition due to crossflow and attachment line instabilities can be studied. Topics covered include: test operating conditions required; design requirements to efficiently make use of the existing infrastructure; development of an injector drive system using a small pilot facility; plenum chamber design; use of computational tools for tunnel and model design; and early operational results.

Description: NASA Ames Research Center is pursuing the development of SOFIA, the Stratospheric Observatory For Infrared Astronomy. SOFIA will consist of a 2.5 meter telescope mounted aft of the wing of a Boeing 747 aircraft. Since a large portion of the infrared spectrum is not visible at ground level due to absorption by water vapor in the atmosphere below 40,000 feet, it is highly desirable to make observations above this altitude. SOFIA will provide the opportunity for astronomers to conduct high-altitude research for extended periods of time. Current study is focused on wind tunnel testing for the open cavity. If not controlled, air would create resonance and damage the telescope. For this reason, SOFIA will design a boundary layer control device to achieve laminar flow over the cavity. This also provides a clearer flow for seeing, thus improving resolution on infrared sources. Other effects being tested in the wind tunnel are aerodynamic torque loads on the telescope, and flutter loads on the tail.

Description: Experimental results for a two-dimensional separated turbulent boundary layer behind a backward facing step for five different Reynolds numbers are reported. Results are presented in the form of tables, graphs and a floppy disk for an easy access of the data. Reynolds number based on the step height was varied by changing the reference velocity upstream of the step, U(sub o), and the step height, h. Hot-wire measurement techniques were used to measure three Reynolds stresses and four triple-velocity correlations. In addition, surface pressure and skin friction coefficients were measured. All hot-wire measurements were acquired in a measuring domain which excluded recirculating flow region due to the directional insensitivity of hot-wires. The downstream extent of the domain from the step was 51 h for the largest and I 14h for the smallest step height. This significant downstream length permitted extensive study of the flow recovery. Prediction of perturbed flows and their recovery is particularly attractive for popular turbulence models since variations of turbulence length and time scales and flow interactions in different regions are generally inadequately predicted. The data indicate that the flow in the free shear layer region behaves like the plane mixing layer up to about 2/3 of the mean reattachment length when the flow interaction with the wall commences the flow recovery to that of an ordinary turbulent boundary layer structure. These changes of the flow do not occur abruptly with the change of boundary conditions. A reattachment region represents a transitional region where the flow undergoes the most dramatic adjustments to the new boundary conditions. Large eddies, created in the upstream free-shear layer region, are being torn, recirculated, reentrained back into the main stream interacting with the incoming flow structure. It is foreseeable that it is quite difficult to describe the physics of this region in a rational and quantitative manner other than statistical. Downstream of the reattachment point the flow recovers at different rates near the wall, in the newly developing internal boundary layer, and in the outer part of the flow. It appears that Reynolds stresses do not fully recover up to the longest recovery length of 114 h.

Description: Tail buffet studies were conducted on a full-scale, production, F/A-18 fighter aircraft in the 80- by 120-Foot Wind Tunnel of the National Full-Scale Aerodynamic Complex at NASA Ames Research Center at Moffett Field, California. Tail buffet data were acquired over an angle-of-attack range of +20 deg to +40 deg, a side-slip range of -16 deg to + 16 deg, and at wind speeds up to 100 knots. The maximum speed corresponds to a Reynolds number of l2.3 x l0(exp 6) based on mean aerodynamic chord and a Mach number of 0. 15. The port, vertical tail fin was instrumented with ninety-six surface-pressure transducers, arranged in six by eight arrays, on each side of the fin. ne aircraft was also equipped with a removable Leading-Edge Extension (LEX) fence whose purpose is to reduce tail-buffet loads. Current analysis methods for the unsteady aerodynamic pressures and loads are described. Only results for the zero side-slip condition are to be presented, both with and without the LEX fence. Results of the time-averaged, power-spectral analysis are presented for the tail fin bending moments which are derived from the integrated pressure field. Local wave velocities on the tail surfaces are calculated from pressure correlations. It was found that the LEX fence significantly reduces the magnitude of the root-mean-square pressures and bending moments. Scaling and repeatability issues are addressed by comparing the present full scale results for pressures at the 60%-span and 45%-chord location with previous full-scale F/A-18 tail-buffet test in the 80- by 120- Foot Wind Tunnel, and with several small-scale tests. The comparisons show that the tail buffet frequency scales very well with tail chord and free-stream velocity, and that there is good agreement with the previous full-scale test. Root-mean-square pressures and power spectra do not scale as well as the frequency results. Addition of a LEX fence caused tail-buffet loads to be reduced at all model scales.

Description: The recent resurgence of interest in utilizing laminar flow on aircraft surfaces for reduction in skin friction drag has generated a considerable amount of research in natural laminar flow (NLF) and hybrid laminar flow control (HLFC) on transonic aircraft wings. This research has focused primarily on airfoil design and understanding transition behavior with little concern for the surface imperfections and manufacturing variations inherent to most production aircraft. In order for laminar flow to find wide-spread use on production aircraft, techniques for constructing the wings must be found such that the large surface imperfections present in the leading edge region of current aircraft do not occur. Toward this end, a modification to existing leading edge construction techniques was devised such that the resulting surface did not contain large gaps and steps as are common on current production aircraft of this class. A lowspeed experiment was first conducted on a simulation of the surface that would result from this construction technique. Preston tube measurements of the boundary layer downstream of the simulated joint and flow visualization using sublimation chemicals validated the literature on the effects of steps on a laminar boundary layer. These results also indicated that the construction technique was indeed compatible with laminar flow. In order to fully validate the compatibility of this construction technique with laminar flow, thus proving that it is possible to build wings that are smooth enough to be used on business jets and light transports in a manner compatible with laminar flow, a flight experiment is being conducted. In this experiment Mach number and Reynolds number will be matched in a real flight environment. The experiment is being conducted using the NASA Dryden F-104 Flight Test Fixture (FTF). The FTF is a low aspect ratio ventral fin mounted beneath an F-104G research aircraft. A new nose shape was designed and constructed for this experiment. This nose shape provides an accelerating pressure gradient in the leading edge region. By flying the aircraft at appropriate Mach numbers and altitudes, this nose shape simulates the leading edge region of a laminar flow wing for a business jet or light transport. Manufactured into the nose shape is a spanwise slot located approximately four inches downstream of the leading edge. The slot, which is an inch wide and one-eighth of an inch deep allows the simulation of surface imperfections, such as gaps and steps at skin joints, which will occur on aircraft using this new construction technique. By placing strips of aluminum of various sizes and shapes in the slot, the effect on the boundary layer of different sizes and shapes of steps and gaps will be examined. It is planned to use five different configurations, differing primarily in the size and number of gaps. Downstream of the slot, the state of the boundary layer is determined using hot film gages and Stanton gages. Agreement between these two very different techniques of measuring boundary layer properties is considered important to being able to state with confidence the effects on the boundary layer of the simulated manufacturing imperfections. To date, the aircraft has not flown. First flights of the aircraft are on schedule to begin October 4, 1993. Low-speed, preliminary experiments at matching Reynolds numbers have been completed.

Description: Available redundancy among aircraft control surfaces allows for effective wing camber modifications. As shown in the past, this fact can be used to improve aircraft performance. To date, however, algorithm developments for in-flight camber optimization have been limited. This paper presents a perturbational approach for cruise optimization through in-flight camber adaptation. The method uses, as a performance index, an indirect measurement of the instantaneous net thrust. As such, the actual performance improvement comes from the integrated effects of airframe and engine. The algorithm, whose design and robustness properties are discussed, is demonstrated on the NASA Dryden B-720 flight simulator.

Description: A highly-instrumented UH-60A aircraft was tested at NASA-Ames Research Center from August 1993 to February 1994 obtaining an extensive data base for level flight, maneuvers, acoustics (both with respect to ground microphone arrays and inflight microphones), and flight dynamics. A majority of the data obtained are now in an electronic data base, however, only a small fraction of the data have been examined. The proposed paper will examine the issue of hovering steadiness in more detail. In particular, a single set of data obtained during ground acoustic testing may provide considerable insight as the wind speeds were measured at a hover height of 250 feet and the aircraft was positioned in 15 deg. steps in heading from 0 to 180 deg. Also, hover housekeeping data were obtained for many of the 31 flights and these will also allow a characterization of the unsteadiness. The variation in section lift will be examined in terms of the induced flow angle variation and this will be related to possible physical explanations.

Description: Compressibility plays a significant role in the development of separation on airfoils experiencing unsteady motion, even at moderately compressible free-stream flow velocities. This effect can result in completely changed stall characteristics compared to those observed at incompressible speed, and can dramatically affect techniques used to control separation. There has been a significant effort in recent years directed toward better understanding; of this process, and its impact on possible techniques for control of separation in this complex environment. A review of existing research in this area will be presented, with emphasis on the physical mechanisms that play such an important role in the development of separation on airfoils. The increasing impact of compressibility on the stall process will be discussed as a function of free-stream Mach number, and an analysis of the changing flow physics will be presented. Examples of the effect of compressibility on dynamic stall will be selected from both recent and historical efforts by members of the aerospace community, as well as from the ongoing research program of the present authors. This will include a presentation of a sample of high speed filming of compressible dynamic stall which has recently been created using real-time interferometry.

Description: This paper will describe the Airbreathing Hypersonic Research Program at NASA Ames Research Center. A main theme will be the "From Computation Through Flight" research effort. General research areas covered will include systems analysis, aerodynamics and aerothermodynamics, propulsion, materials, and flight research. Illustrative results from each discipline will be presented. The synergism between computational and experimental research will be demonstrated by examples. All examples given will have been published in the open literature.

Description: A developed method has been applied to calculate accurately the viscous flow about airfoils normal to the free-stream flow. This method has special application to the analysis of tilt rotor aircraft in the evaluation of download. In particular, the flow about an XV-15 airfoil with and without deflected leading and trailing edge flaps at -90 degrees incidence is evaluated. The multi-element aspect of the method provides for the evaluation of slotted flap configurations which may lead to decreased drag. The method solves for turbulent flow at flight Reynolds numbers. The flow about the XV-15 airfoil with and without flap deflections has been calculated and compared with experimental data at a Reynolds number of one million. The comparison between the calculated and measured pressure distributions are very good, thereby, verifying the method. The aerodynamic evaluation of multielement airfoils will be conducted to determine airfoil/flap configurations for reduced airfoil drag. Comparisons between the calculated lift, drag and pitching moment on the airfoil and the airfoil surface pressure will also be presented.

Description: Steady and unsteady viscous, three-dimensional flowfields are calculated using a thin layer approximation of Navier-Stokes equations in conjunction with Chimera overset grids. The finite-difference numerical scheme uses structured grids and a pentadiagonal flow solver called "OVERFLOW". The configuration of Boeing 747-200 has been chosen as one of configurations to be used as a platform for the SOFIA (Stratospheric Observatory For Infrared Astronomy). Initially, the steady flowfield of the full aircraft is calculated for the clean configuration (without a cavity to house telescope). This solution is then used to start the unsteady flowfield of a configuration containing cavity housing the observation telescope and its peripheral units. Analysis of unsteady flowfield in the cavity and its influence on the tail empennage, as well as the noise due to turbulence and optical quality of the flow are the main focus of this study. For the configuration considered here, the telescope housing cavity is located slightly downstream of the portwing. The entire flow-field is carefully constructed using 45 overset grids and consists of nearly 4 million grid points. All the computations axe done at one freestream flow condition of M(sub infinity) = 0.85, alpha = 2.5deg, and a Reynolds of Re = 1.85x10deg

Description: The proposed paper presents flow visualization performed during experiments conducted on a full-scale F/A-18 aircraft in the 80- by 120-Foot Wind-Tunnel at NASA Ames Research Center. This investigation used both surface and off-surface flow visualization techniques to examine the flow field on the forebody, canopy, leading edge extensions (LEXs), and wings. The various techniques used to visualize the flow field were fluorescent tufts, flow cones treated with reflective material, smoke in combination with a laser light sheet, and a video imaging system. The flow visualization experiments were conducted over an angle of attack range from 20deg to 45deg and over a sideslip range from -10deg to 10deg. The results show regions of attached and separated flow on the forebody, canopy, and wings. Additionally, the vortical flow is clearly visible over the leading-edge extensions, canopy, and wings.

Description: It is stated that the aerodynamic forces on the vehicle being aerocaptured are controlled by "altering the angle of attack" and thereby controlling the lift coefficient. Furthermore, the resulting variation of drag coefficient with angle of attack was ignored. The purpose of this Comment is to point out that an aerodynamic control method that is much more effective than the pitch modulation has been studied and utilized during entries for many years. During aerocapture, it is desirable to have a large range of lift coefficients available, while keeping the vehicle's ballistic coefficients constant. This is accomplished by modulating the vehicle's bank angle, i.e., by rolling the vehicle about its velocity vector. By this method, the angle of attack can be held constant (at the trim angle, if desired), and the C(sub D) and the ballistic coefficient remain constant. Furthermore, the vertical component of the normal force vector (essentially the lift) can be varied over its entire range, from maximum positive to maximum negative values. Reaction controls, rather than aerodynamic ones, are usually utilized to change the bank angle of the vehicle, thus requiring the use of fuel. However, the fuel expenditure that is required to change the bank angle is far less than the amount that would have to be used to continuously hold the vehicle at pitch angles that differ significantly from its trim angle of attack. Also, it has been shown that bank angle modulation to vary the lift can enlarge the entry corridor by increasing the entry angle for the undershoot boundary, where both the heating rate and deceleration reach a maximum. Finally, the crew's deceleration tolerance can be increased somewhat when the bank angle is varied, as opposed to the pitch angle. For bank modulation, the deceleration force vector can be kept at a constant angle with respect to the occupants whose tolerance to g loads is highest when the force is applied in a direction normal to the upper torso. The advantages of bank angle variation to modulate the lift vector were recognized long ago, and this method of control was used successfully on the Apollo command module during lunar return' and, more recently, for the Space Shuttle Orbiter.

Description: Study of sonic and supersonic jet plumes are relevant to understanding such phenomenon as jet-noise, plume signatures, and rocket base-heating and radiation. Jet plumes are simple to simulate and yet, have complex flow structures such as Mach disks, triple points, shear-layers, barrel shocks, shock-shear-layer interaction, etc. Experimental and computational simulation of sonic and supersonic jet plumes have been performed for under- and over-expanded, axisymmetric plume conditions. The computational simulation compare very well with the experimental observations of schlieren pictures. Experimental data such as temperature measurements with hot-wire probes are yet to be measured and will be compared with computed values. Extensive analysis of the computational simulations presents a clear picture of how the complex flow structure develops and the conditions under which self-similar flow structures evolve. From the computations, the plume structure can be further classified into many sub-groups. In the proposed paper, detail results from the experimental and computational simulations for single, axisymmetric, under- and over-expanded, sonic and supersonic plumes will be compared and the fluid dynamic aspects of flow structures will be discussed.

Description: Three direct numerical simulations of time-evolving turbulent plane wakes with velocity deficit Reynolds numbers of about 2,000 have been simulated using a spectral numerical method with up to 600 x 260 x 160 modes. The initial conditions for the simulations are generated from direct numerical simulations of a turbulent boundary layer (momentum thickness Reynolds number of 670), and varying amounts of additional two- dimensional, forcing. In order to preserve the self-similar flow evolution, the forcing is implemented by multiplying all the two-dimensional modes in the initial condition by a constant factor. In the "natural" case no additional forcing is used; in the "forced" and "heavily forced" cases this factor is 5 and 20, respectively. The wake spreading rate Is increased by factors of 1.7 and 7.1 for the two forced cases. The Reynolds stresses are also increased by a similar or even larger factor. These results indicate that the plane wake is much more sensitive to initial forcing than the plane mixing layer. As in the plane mixing layer, two-dimensional forcing promotes more organized large-scale vortical flow structures and these structures axe sometimes separated by "braid regions" containing streamwise "rib" vortices, unlike in the unforced wake.

Description: Large-eddy simulation of the incompressible Navier-Stokes equations has been used to examine the long-time development of initially isotropic turbulence subjected to solid-body rotation. The simulations were carried out using a pseudo-spectral method with 128 x 128 x 512 collocation points in a computational domain that is four times larger along the rotation axis than in the other directions; subgrid-scale motions were parameterized using a spectral eddy viscosity model modified for system rotation. Simulation results show that the correlation length along the rotation am's of velocities orthogonal to the rotation vector exhibits rapid growth while the integral length-scale of velocities aligned with the rotation axis is relatively unaffected by rotation. Examination of the energy spectrum of two-dimensional, two-component motions indicates the presence of an inverse cascade of energy. System rotation also causes an alignment of vorticity along the rotation axis with relatively stronger cyclonic vorticity than anticyclonic. The onset of anisotropic effects are well characterized by Rossby numbers defined in terms of both macroscopic and microscopic quantities.

Description: This paper will review the advances made recently in the Navier-Stokes CFD methods to simulate aerodynamics and aeroacoustics of helicopter rotors and rotor-body flows. Although a complete flowfield simulation of full helicopter is currently not feasible with these methods, impressive gains have been made in analyzing individual components of this complex problem in a very detailed manner. The use of the state-of-the-art numerical algorithms in solution methods, in conjunction with powerful supercomputers, like the Cray-2, have enabled noticeable progress to be made in modeling viscous-inviscid interactions, blade-vortex interactions, tip-vortex: simulation and wake effects, as well as high speed impulsive noise in hover and forward flight for isolated rotor blades. This paper will critically evaluate the presently available Euler and Navier-Stokes methods, both finite-difference and finite volume methods using structured and unstructured grids for helicopter applications for accuracy, suitability, and computational efficiency. The review will also include the recent progress made using overset grids to model rotor-body flows. All the material for this review will be drawn from the published material shown below.

Description: In recent years significant advances have been made for parallel computers in both hardware and software. Now parallel computers have become viable tools in computational mechanics. Many application codes developed on conventional computers have been modified to benefit from parallel computers. Significant speedups in some areas have been achieved by parallel computations. For single-discipline use of both fluid dynamics and structural dynamics, computations have been made on wing-body configurations using parallel computers. However, only a limited amount of work has been completed in combining these two disciplines for multidisciplinary applications. The prime reason is the increased level of complication associated with a multidisciplinary approach. In this work, procedures to compute aeroelasticity on parallel computers using direct coupling of fluid and structural equations will be investigated for wing-body configurations. The parallel computer selected for computations is an Intel iPSC/860 computer which is a distributed-memory, multiple-instruction, multiple data (MIMD) computer with 128 processors. In this study, the computational efficiency issues of parallel integration of both fluid and structural equations will be investigated in detail. The fluid and structural domains will be modeled using finite-difference and finite-element approaches, respectively. Results from the parallel computer will be compared with those from the conventional computers using a single processor. This study will provide an efficient computational tool for the aeroelastic analysis of wing-body structures on MIMD type parallel computers.

Description: This paper summarizes the technical development and validation of a multiphase computational fluid dynamics (CFD) numerical method using the volume-of-fluid (VOF) model and a Lagrangian tracking model which can be employed to analyze general multiphase flow problems with free surface mechanism. The gas-liquid interface mass, momentum and energy conservations are modeled by continuum surface mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed flow regimes. The objectives of the present study are to develop and verify the fractional volume-of-fluid cell partitioning approach into a predictor-corrector algorithm and to demonstrate the effectiveness of the present innovative approach by simulating benchmark problems including the coaxial jet atomization.

Description: The performance of a workstation cluster used for the solution of the Reynolds-averaged Navier-Stokes equations is compared with a conventional vector supercomputer architecture. The application simulation of the steady flowfield about a transonic transport was computed using an implicit diagonal scheme in an overset mesh framework. Static load balancing was used, while coarse grain decomposition was achieved by solution of a grid zone per processor. Price/performance ratios are estimated for several scenarios in which such clusters may be utilized.

Description: The essentially nonoscillatory (ENO) shock-capturing scheme for the solution of hyperbolic equations is extended to solve a system of coupled conservation equations governing two-dimensional, time-dependent, compressible chemically reacting flow with full chemistry. The thermodynamic properties of the mixture are modeled accurately, and stiff kinetic terms are separated from the fluid motion by a fractional step algorithm. The methodology is used to study the concept of shock-induced mixing and combustion, a process by which the interaction of a shock wave with a jet of low-density hydrogen fuel enhances mixing through streamwise vorticity generation. Test cases with and without chemical reaction are explored here. Our results indicate that, in the temperature range examined, vorticity generation as well as the distribution of atomic species do not change significantly with the introduction of a chemical reaction and subsequent heat release. The actual diffusion of hydrogen is also relatively unaffected by the reaction process. This suggests that the fluid mechanics of this problem may be successfully decoupled from the combustion processes, and that computation of the mixing problem (without combustion chemistry) can elucidate much of the important physical features of the flow.

Description: Multicomponent gas diffusion is reviewed with particular emphasis on gas flows near solid boundaries-the so-called Kramers-Kistemaker effect. The aim is to derive an appropriate momentum boundary condition which governs many gaseous species diffusing together. The many species' generalization of the traditional single gas condition, either as slip or stick (no-slip), is not obvious, particularly for technologically important cases of lower gas pressures and very dissimilar molecular weight gases. No convincing theoretical case exists for why two gases should interact with solid boundaries equally but in opposite flow directions, such that the total gas flow exactly vanishes. ln this way, the multicomponent no-slip boundary requires careful treatment The approaches discussed here generally adopt a microscopic model for gas-solid contact. The method has the advantage that the mathematics remain tractable and hence experimentally testable. Two new proposals are put forward, the first building in some molecular collision physics, the second drawing on a detailed view of surface diffusion which does not unphysically extrapolate bulk gas properties to govern the adsorbed molecules. The outcome is a better accounting of previously anomalous experiments. Models predict novel slip conditions appearing even for the case of equal molecular weight components. These approaches become particularly significant in view of a conceptual contradiction found to arise in previous derivations of the appropriate boundary conditions. The analogous case of three gases, one of which is uniformly distributed and hence non-diffusing, presents a further refinement which gives unexpected flow reversals near solid boundaries. This case is investigated alone and for aggregating gas species near their condensation point. In addition to predicting new physics, this investigation carries practical implications for controlling vapor diffusion in the growth of crystals used in medical diagnosis (e.g. mercuric iodide) and semiconductors.

Description: The aim of the present investigation is to characterize the motion of dendrite fragments falling under the influence of gravity in a uniform liquid medium at low Reynolds number. In an earlier study, Zakhem, Weidman and de Groh (1992) reported on the settling speed of model equiaxed dendrite grains released along their axis of symmetry. In this follow-up study uniaxial model dendrite grains were released off-axis to observe and document their motion at different orientations. It was hypothesized that the dendrite models might rotate when released off-axis in which case an attempt would be made to document the ensuing unsteady motion. This latter event turned out to be in fact true: at the small but finite Reynolds numbers that existed, each uniaxial dendrite slowly rotated towards its equilibrium orientation while failing under the influence of gravity. In addition to completing the original goal, we have made use of a beads-on-a shell Stokes flow code to numerically determine the drag coefficient for capsules, i.e.. uniaxial dendrites without arms. The drag on horizontally and vertically falling capsules are reported and compared with measurements.

Description: Supersonic jet plumes were studied using a two-equation turbulence model employing corrections for compressible dissipation and pressure-dilatation. A space-marching procedure based on an upwind numerical scheme was used to solve the governing equations and turbulence transport equations. The computed results indicate that two-equation models employing corrections for compressible dissipation and pressure-dilatation yield improved agreement with the experimental data. In addition, the numerical study demonstrates that the computed results are sensitive to the effect of grid refinement and insensitive to the type of velocity profiles used at the inflow boundary for the cases considered in the present study.

Description: Results are obtained for cylindrical leading edges of proposed transatmospheric vehicles by employing a two-dimensional viscous shock-layer code for nonequilibrium gas flows. The accuracy and efficiency of the planar code is verified through detailed comparisons with other predictions. This study includes results for 6-deg half-angle bodies with nose radii ranging from 0.01 to 2.0 ft for both cylindrically blunted wedges and spherically blunted cones (included for comparison). Some results are presented as a ratio of the noncatalytic to the corresponding fully catalytic heating value to illustrate the maximum potential for a heating reduction in dissociated nonequilibrium flows. Generally, this ratio and the individual heating rates are smaller for cylindrically blunted wedges with small nose radii as compared to the spherically blunted cones (for the same nose radius). Therefore, a larger potential exists for heating reduction in cylindrically blunted as compared with the spherically blunted surfaces. However, the results presented at higher altitudes (where the slip effects become important) show that the spherically, blunted nose gives lower stagnation-point heating due to stronger merged shock-layer effects as compared with a cylindrically blunted nose.

Description: The three-dimensional Reynolds-averaged, Navier-Stokes (RANS) equations are used to numerically simulate nonsteady vortical flow about a 65 degree sweep delta wing at 30 degrees angle of attack. Two large-amplitude, high-rate, forced-roll motions and a damped free-to-roll motion are presented. The free-to-roll motion is computed by coupling the time-dependent RANS equations to the flight dynamic equation of motion. The computed results are compared with experimental forces, moments, and roll-angle time histories. The overall agreement is good. Vortex breakdown is present in each case, which causes significant time lags in the vortex breakdown motions relative to the body motions. This behavior strongly influences the dynamic forces and moments.

Description: A study is described that evaluates the accuracy of vortex-lattice methods when they are used to compute the loads induced on aircraft as they encounter lift-generated wakes. The evaluation is accomplished by use of measurements made in the 80- by 120-foot wind tunnel of the lift, rolling-moment, and downwash in the wake of three configurations of a model of a subsonic transport aircraft. The downwash measurements are used as input for a vortex-lattice code in order to compute the lift and rolling moment induced on wings that have a span of 0.186, 0.510, or 1.022 times the span of the wake-generating model. Comparison of the computed results with the measured lift and rolling moment distributions are used to determine the accuracy of the vortex-lattice code. It was found that the vortex-lattice method is very reliable as long as the span of the encountering of following wing is less than about 0.2 of the generator span. As the span of the following wing increases above 0.2, the vortex-lattice method continues to correctly predict the trends and nature of the induced loads, but it overpredicts the magnitude of the loads by increasing amounts. The increase in deviation of the computed from the measured loads with size of the following wing is attributed to the increase in distortion of the structure of the vortex wake as it approaches and passes the larger following wings.

Description: In the absence of gravity drainage, froth wicking draws excess fluid onto a bubble lattice. Capillary forces only cause fluid transport; a moving front moves stably and without fluid fingering along a constant velocity bubble-fluid contact line. This percolation of fluid crawling up the lattice shows fluid coverage on lattice borders varies linearly with available surface area (proportional to lattice perimeter) and fluid accelerates through regions or nests of high bubble density (number of bubbles/sq cm). The development of nearly two-dimensional bubble lattices in variable gravity (step function between 0.01 and 1.8 times earthly gravity) are examined experimentally and a zeroth-order model for froth wetting is presented, which captures many of the principal observations. Possible applications for bubble lattices include adhesion casting of metals and separation of biological cells, bacteria, and particles.

Description: The viscous flow field near the surface of a hovering rotor blade was studied for blade twist distributions typical of a till rotor blade and a conventional helicopter rotor blade. Three blade geometries were studied, including a tilt rotor blade twist distribution (baseline), conventional helicopter rotor blade twist distribution, and the baseline twist distribution with 2 deg of precone. The results give insight into the delayed stall phenomenon often observed for highly twisted rotors. Calculations were performed for a high thrust condition near stall using the thin-layer Navier-Stokes CFD code TURNS. Effects of built-in twist on section force coefficients, skin friction, velocities, surface pressures, and boundary layer shape factor are discussed. Although the rotor thrust coefficient was nominally the same for the cases using the two twist distributions, large differences were found in the section in-plane and normal force coefficients. These preliminary results imply that the blade outboard region, rather than the inboard region, provides the majority of the performance advantage of the baseline case over the low twist case. Skin friction, velocities near the blade, and surface pressures for the two twist distributions reveal significant differences in the blade outboard region.

Description: Previous experiments indicate that the bubble point pressure of spacecraft liquid hydrogen acquisition devices is reduced substantially when the ullage is pressurized with heated hydrogen vapor. The objective is to determine whether the two-phase thermocapillary convection arising from thermodynamic non-equilibrium along the porous surfaces of such devices could lead to this observed degradation in retention performance. We also examine why retention capability appears to be unaffected by pressurization with heated helium or direct heating through the porous structure. Computational assessments based on coupled solution of the flowfield and liquid free surface indicate that for highly wetting fluids in small pores, dynamic pressure and vapor recoil dictate surface morphology and drive meniscus deformation. With superheating, the two terms exert the same influence on curvature and promote mechanical equilibrium, but with subcooling, the pressure distribution produces a suction about the pore center-line that degrades retention. This result points to thermocapillary-induced deformation arising from condensation as the cause for retention loss. It also indicates that increasing the level of non-equilibrium by reducing accommodation coefficient restricts deformation and explains why retention failure does not occur with direct screen heating or helium pressurization.

Description: New capabilities have been added to a Navier-Stokes solver to perform steady-state simulations more efficiently. The flow solver for solving the Navier-Stokes equations is completely rewritten with a combination of the LU-SGS (Lower-Upper factored Symmetric Gauss-Seidel) implicit method and the modified HLLE (Harten-Lax-van Leer-Einfeldt) upwind scheme. A pseudo-time marching method is used for the directly coupled structural equations to improve overall convergence rates for static aeroelastic analysis. Results are demonstrated for transonic flows over rigid and flexible wings.

Description: A hybrid method for computing compressible viscous flows is presented. This method divides the computational domain into two zones. In the outer zone, the unsteady full-potential equation (FPE) is solved. In the inner zone, the Navier-Stokes equations are solved. The two zones are tightly coupled so that steady and unsteady flows may be efficiently solved. The resulting CPU times are less than 50 percent of the required for a full-blown Navier-Stokes analysis. Sample applications of the method to an unswept iced wing at 4 deg and 8 deg angle of attack are presented. Surface pressures are in good agreement with the measurements obtained by Bragg et al. at the University of Illinois.

Description: Strong interactions between flow about an aircraft wing and the wing structure can result in aeroelastic phenomena which significantly impact aircraft performance. Time-accurate methods for solving the unsteady Navier-Stokes equations have matured to the point where reliable results can be obtained with reasonable computational costs for complex non-linear flows with shock waves, vortices and separations. The ability to combine such a flow solver with a general finite element structural model is key to an aeroelastic analysis in these flows. Earlier work involved time-accurate integration of modal structural models based on plate elements. A finite element model was developed to handle three-dimensional wing boxes, and incorporated into the flow solver without the need for modal analysis. Static condensation is performed on the structural model to reduce the structural degrees of freedom for the aeroelastic analysis. Direct incorporation of the finite element wing-box structural model with the flow solver requires finding adequate methods for transferring aerodynamic pressures to the structural grid and returning deflections to the aerodynamic grid. Several schemes were explored for handling the grid-to-grid transfer of information. The complex, built-up nature of the wing-box complicated this transfer. Aeroelastic calculations for a sample wing in transonic flow comparing various simple transfer schemes are presented and discussed.

Description: A model of the Shuttle Orbiter rarefied-flow aerodynamic force coefficients has been derived from the ratio of flight acceleration measurements. The in-situ, low-frequency (less than 1Hz), low-level (approximately 1 x 10(exp -6) g) acceleration measurements are made during atmospheric re-entry. The experiment equipment designed and used for this task is the High Resolution Accelerometer Package (HiRAP), one of the sensor packages in the Orbiter Experiments Program. To date, 12 HiRAP re-entry mission data sets spanning a period of about 10 years have been processed. The HiRAP-derived aerodynamics model is described in detail. The model includes normal and axial hypersonic continuum coefficient equations as function of angle of attack, body-flap deflection, and elevon deflection. Normal and axial free molecule flow coefficient equations as a function of angle of attack are also presented, along with flight-derived rarefied-flow transition bridging formulae. Comparisons are made between the aerodynamics model, data from the latest Orbiter Operational Aerodynamic Design Data Book, applicable computer simulations, and wind-tunnel data.

Description: Generation of significant side forces and yawing moments on an F/A-18 fuselage through tangential slot blowing is analyzed using computational fluid dynamics. The effects of freestream Mach number, jet exit conditions, jet length, and jet location are studied. The effects of over- and under-blowing on force and moment production are analyzed. Non-time-accurate solutions are obtained to determine the steady-state side forces, yawing moments, and surface pressure distributions generated by tangential slot blowing. Time-accurate solutions are obtained to study the force onset time lag of tangential slot blowing. Comparison with available experimental data from full-scale wind tunnel and sub-scale wind tunnel tests are made. This computational analysis complements the experimental results and provides a detailed understanding of the effects of tangential slot blowing on the flow field about the isolated F/A-18 forebody. Additionally, it extends the slot-blowing database to transonic maneuvering Mach numbers.

Description: It is well known that radiative heat transport influences many types of buoyant flows due to its effect on the temperature and thus density field in the fluid medium. It is of interest to study gaseous flows driven solely by radiation in the absence of buoyancy, particularly because of its application to astrophysical flows that are well known from astronomical observations and numerical simulation. However, no laboratory-scale experiments of this phenomenon have ever been conducted. To study the possibility of obtaining such flows in the laboratory, an apparatus was built to produce large temperature differences (Delta T) up to 300 K in a gas confined between flat parallel plates. SF6 was used as the radiatively-active gas because its Planck absorption length is much shorter than that of any other common non-reactive gas. The NASA-Lewis 2.2 second drop tower was used to obtain reduced gravity in order to suppress buoyancy effects. To image the resulting flows, a laser shearing interferometer was employed. Initial results indicate the presence of flow that does not appear to be attributable to the residual flow resulting from buoyancy influences before the drop. For Delta T greater than 70 K, slight deformations in the interferometer fringes seen at lower Delta T became large unsteady swirls. Such behavior did not occur for radiatively-inactive gases, suggesting that a flow driven solely by radiation was obtained in SF6 and to a lesser extent in CO2 This was more pronounced at higher pressures and plate spacings, consistent with our scaling predictions.

Description: This paper is concerned with the downstream evolution of a resonant triad of initially non-interacting linear Instability waves in a boundary layer with a weak adverse pressure gradient. The triad consists of a two-dimensional fundamental mode and a pair of equal-amplitude oblique modes that form a subharmonic standing wave in the spanwise direction. The growth rates are small and there is a well-defined common critical layer for these waves. As in Goldstein & Lee (1992), the wave interaction takes place entirely within this critical layer and is initially of the parametric-resonance type. This enhances the spatial growth rate of the subharmonic but does not affect that of the fundamental. However, in contrast to Goldstein & Lee (1992), the initial subharmonic amplitude is assumed to be small enough so that the fundamental can become nonlinear within its own critical layer before it is affected by the subharmonic. The subharmonic evolution is then dominated by the parametric-resonance effects and occurs on a much shorter streamwise scale than that of the fundamental. The subharmonic amplitude continues to increase during this parametric-resonance stage - even as the growth rate of the fundamental approaches zero - and the subharmonic eventually becomes large enough to influence the fundamental which causes both waves to evolve on the same shorter streamwise scale.

Description: We consider the effect of an inert gas on the diffusive-convective physical vapor transport process. We investigate the case when the temperature gradient is stabilizing and the concentration gradient is destabilizing for a wide parametric range. When an inert gas is present, the thermal and solutal convection oppose each other. The solutal field is destabilizing while the thermal field and the advective-diffusive flux stabilize the flow field. When the pressure of the inert component is increased, the stabilizing effect of the advective-diffusive flux is decreased. The intensity of convection as well as the oscillatory transient time increases. Below, the critical Rayleigh number, the nonlinear dynamics of the flow field show an oscillatory approach to steady state. For parametric values in the neighborhood of the critical Rayleigh number, the flow field undergoes a chaotic transient which settles to a periodic state. The asymptotic state of the flow field shows that growth and amalgamation of cells yields an overturning motion which results in an asymmetric cellular structure. The low gravity environment yields the stabilizing advective-diffusive flow which results in uniform temperature and concentration gradients near the crystal interface.

Description: An experimental survey of supersonic wing tip vortices has been conducted at Mach 2.5 using small performed 2.25 chords down-stream of a semi-span rectangular wing at angle of attack of 5 and 10 degrees. The main objective of the experiments was to determine the Mach number, flow angularity and total pressure distribution in the core region of supersonic wing tip vortices. A secondary aim was to demonstrate the feasibility of using cone probes calibrated with a numerical flow solver to measure flow characteristics at supersonic speeds. Results showed that the numerically generated calibration curves can be used for 4-hole cone probes, but were not sufficiently accurate for conventional 5-hole probes due to nose bluntness effects. Combination of 4-hole cone probe measurements with independent pitot pressure measurements indicated a significant Mach number and total pressure deficit in the core regions of supersonic wing tip vortices, combined with an asymmetric 'Burger like' swirl distribution.

Description: Three direct numerical simulations of incompressible turbulent plane mixing layers have been performed. All the simulations were initialized with the same two velocity fields obtained from a direct numerical simulation of a turbulent boundary layer with a momentum thickness Reynolds number of 300 computed by Spalart (J. Fluid Mech. 187, 61, 1988). In addition to a baseline case with no additional disturbances, two simulations were begun with two-dimensional disturbances of varying strength in addition to the boundary layer turbulence. After a development stage, the baseline case and the case with weaker additional two-dimensional disturbances evolve self-similarly, reaching visual thickness Reynolds numbers of up to 20 000. This self-similar period is characterized by a lack of large-scale organized pairings, a lack of streamwise vortices in the 'braid' regions, and scalar mixing that is characterized by 'marching' Probability Density Functions (PDFs). The case begun with strong additional two-dimensional disturbances only becomes approximately self-similar, but exhibits sustained organized large-scale pairings, clearly defined braid regions with streamwise vortices that span them, and scalar PDFs that are 'nonmarching.' It is also characterized by much more intense vertical velocity fluctuations than the other two cases. The statistics and structures in several experiments involving turbulent mixing layers are in better agreement with those of the simulations that do not exhibit organized pairings.

Description: A new computational efficient 3-D compressible Reynolds-averaged implicit Navier-Stokes method with advanced two equation turbulence models for high speed flows is presented. All convective terms are modeled using an entropy satisfying higher-order Total Variation Diminishing (TVD) scheme based on implicit upwind flux-difference split approximations and arithmetic averaging procedure of primitive variables. This method combines the best features of data management and computational efficiency of space marching procedures with the generality and stability of time dependent Navier-Stokes procedures to solve flows with mixed supersonic and subsonic zones, including streamwise separated flows. Its robust stability derives from a combination of conservative implicit upwind flux-difference splitting with Roe's property U to provide accurate shock capturing capability that non-conservative schemes do not guarantee, alternating symmetric Gauss-Seidel 'method of planes' relaxation procedure coupled with a three-dimensional two-factor diagonal-dominant approximate factorization scheme, TVD flux limiters of higher-order flux differences satisfying realizability, and well-posed characteristic-based implicit boundary-point a'pproximations consistent with the local characteristics domain of dependence. The efficiency of the method is highly increased with Newton Raphson acceleration which allows convergence in essentially one forward sweep for supersonic flows. The method is verified by comparing with experiment and other Navier-Stokes methods. Here, results of adiabatic and cooled flat plate flows, compression corner flow, and 3-D hypersonic shock-wave/turbulent boundary layer interaction flows are presented. The robust 3-D method achieves a better computational efficiency of at least one order of magnitude over the CNS Navier-Stokes code. It provides cost-effective aerodynamic predictions in agreement with experiment, and the capability of predicting complex flow structures in complex geometries with good accuracy.

Description: The ability of one- and two-equation turbulence models to predict unsteady separated flows over airfoils is evaluated. An implicit, factorized, upwind-biased numerical scheme is used for the integration of the compressible, Reynolds averaged Navier-Stokes equations. The turbulent eddy viscosity is obtained from the computed mean flowfield by integration of the turbulent field equations. The two-equation turbulence models are discretized in space with an upwind-biased, second order accurate total variation diminishing scheme. One and two-equation turbulence models are first tested for a separated airfoil flow at fixed angle of incidence. The same models are then applied to compute the unsteady flowfields about airfoils undergoing oscillatory motion at low subsonic Mach numbers. Experimental cases where the flow has been tripped at the leading edge and where natural transition was allowed to occur naturally are considered. The more recently developed field-equation turbulence models capture the physics of unsteady separated flow significantly better than the standard kappa-epsilon and kappa-omega models. However, certain differences in the hysteresis effects are obtained. For an untripped high-Reynolds-number flow, it was found necessary to take into account the leading edge transitional flow region in order to capture the correct physical mechanism that leads to dynamic stall.

Description: Digital flight records from reported clear-air turbulence incidents are used to determine winds and turbulence, to determine maneuver g loads, and to analyze control problems. Many cases of severe turbulence are found downwind of mountains and thunderstorms where sharp, sudden jolts are associated with vortices in atmospheric waves. Other cases of severe turbulence are round in strong updrafts above thunderstorm buildups that may be undetected by onboard weather radar. An important finding is that there are large maneuvering loads in over half of the reported clear-air turbulence incidents. Maneuvering loads are determined through an analysis of the short-term variations in elevator deflection and aircraft pitch angle. For altitude control in mountain waves the results indicate that small pitch angle changes with proper timing are sufficient to counter variations in vertical wind. For airspeed control in strong mountain waves, however, there is neither the available thrust nor the quickness in engine response necessary to counter the large variations in winds.

Description: The fluid mechanics of air flow exiting from porous materials is investigated. Porous walls with filter ratings differing by about three orders of magnitude are studied. The flow behavior is investigated for its spatial and temporal stability. The measurements show that coalescence effects cause jet development which gives rise to a stable, spatial velocity non-uniformity termed pseudo-turbulence. Beyond a critical mean velocity, the pseudo-turbulence shows a leveling trend and the flow exhibits jets whose dimensions are independent of wall porosity and provide a unique length scale for this type of flow field (approx. 1 mm). A correlation is developed to predict the onset of fully developed jets. The impact of using porous wall materials in cold flow solid rocket motor testing is addressed.

Description: The Essentially NonOscillatory (ENO) shock-capturing scheme for the solution of hyperbolic equations is extended to solve a system of coupled conservation equations governing two-dimensional, time-dependent, compressible chemically reacting flow with full chemistry. The thermodynamic properties of the mixture are modeled accurately, and stiff kinetic terms are separated from the fluid motion by a fractional step algorithm. The methodology is used to study the concept of shock-induced mixing and combustion, a process by which the interaction of a shock wave with a jet of low-density hydrogen fuel enhances mixing through streamwise vorticity generation. Test cases with and without chemical reaction are explored here. Our results indicate that, in the temperature range examined, vorticity generation as well as the distribution of atomic species do not change significantly with the introduction of a chemical reaction and subsequent heat release. The actual diffusion of hydrogen is also relatively unaffected by the reaction process. This suggests that the fluid mechanics of this problem may be successfully decoupled from the combustion processes, and that computation of the mixing problem (without combustion chemistry) can elucidate much of the important physical features of the flow.

Description: A computer code nozzle in n-temperatures (NOZNT), which calculates one-dimensional flows of partially dissociated and ionized air in an expanding nozzle, is tested against three existing sets of experimental data taken in arcjet wind tunnels. The code accounts for the differences among various temperatures, i.e., translational-rotational temperature, vibrational temperatures of individual molecular species, and electron-electronic temperature, and the effects of impurities. The experimental data considered are (1) the spectroscopic emission data; (2) electron beam data on vibrational temperature; and (3) mass-spectrometric species concentration data. It is shown that the impurities are inconsequential for the arcjet flows, and the NOZNT code is validated by numerically reproducing the experimental data.

Description: Experimental data on pressure distribution and heat transfer on a turbine airfoil were obtained over a range of Reynolds numbers from 0.75 to 7.0 x 10(exp 6) and a range of turbulence intensities from 1.8 to about 15%. The purpose of this study was to obtain fundamental heat transfer and pressure distribution data over a wide range of high Reynolds numbers and to extend the heat transfer data base to include the range or Reynolds numbers encountered in the Space Shuttle main engine turbopump turbines. The results of this study indicated that Reynolds number and turbulence intensity have a large effect on both the transition from laminar to turbulent flow and the resulting heat transfer. For a given turbulence intensity, heat transfer for all Reynolds numbers at the leading edge can generally be correlated with the Frossling number developed for lower Reynolds numbers. For a given turbulence intensity, heat transfer for the airfoil surfaces downstream of the leading edge can be approximately correlated with a dimensionless parameter. Comparisons of the experimental results were also made with a numerical solution from a two-dimensional Navier-Stokes code.

Description: New boundary-layer equations are developed and the solutions match all the boundary-layer properties, except for the normal velocity, exactly with the corresponding inviscid properties. The numerical procedure solves tridiagonal matrices at each marching station. As part of the solution an inviscid transpiration velocity at the surface is calculated from the boundary-layer solution. This transpiration velocity could be used as a boundary condition to calculate a new inviscid solution. If the inviscid/boundary layer solutions are iterated, then the normal velocity from the boundary layer solution will match the inviscid values exactly. Solutions are calculated for shear flows over a flat plate. Results from the present method compared well with Navier-Stokes solutions for incompressible constant shear flow and sinusoidal shear flow. Iterations of the inviscid/boundary layer solutions were not necessary. Compressible flow of an exponential shear at Mach 4 was calculated for cold, adiabatic, and hot walls. It was found that inviscid shear flows have a significant effect on skin friction and heating rates. The present method should yield accurate boundary-layer solutions at Reynolds numbers lower than the traditional boundary-layer equations.

Description: Two different approaches, the direct simulation Monte Carlo (DSMC) method based on molecular gasdynamics, and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gasdynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental data, in which density and rotational temperature variations along the centerline and at various locations inside a low-density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data, both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity-temperature exponent, have substantial effects on the results of the DSMC method.

Description: The effects of convection on diffusive-convective physical vapor transport process are examined computationally. We analyze conditions ranging from typical laboratory conditions to conditions achievable only in a low gravity environment. This corresponds to thermal Rayleigh numbers Ra, ranging from 1.80 x 10 to 1.92 x 10(exp 6). Our results indicate that the effect of the sublimation and condensation fluxes at the boundaries is to increase the threshold of instability. For typical ground based conditions, time dependent oscillatory convection can occur. This results in unsteady transport, and non- uniform temperature and concentration gradients at the crystal interface. Spectral analysis of the flow field shows parametric regions exhibiting both an oscillatory approach to steady state and a chaotic transient to a periodic state. Low gravity conditions stabilize the flow field. Convective effects are effectively reduced, thus resulting in uniform temperature and concentration gradients at the interface, a desirable condition for crystal growth.

Description: The nonlinear interfacial instability of a liquid jet in a coflowing compressible airstream is studied numerically. A high-resolution scheme which has second-order accuracy in space and time is coupled with a Lagrangian marker particle algorithm to visualize the large-scale motion of the interfaces in compressible flow. A numerical algorithm based on an approximate equation of state of a compressible liquid is developed to allow this two-fluid system to be governed by the nonlinear unsteady Euler equations in conservative form. The initial growth of small disturbances given by the simulations agrees well with linear theory. The process of jet disruption in compressible flow is demonstrated to consist of the formation of liquid spikes, interweaving of the gas and liquid and stretching and detachment of the liquid main center core.

Description: The effect of forebody tangential slot blowing on the flowfield about an F/A-18 aircraft is investigated numerically using solutions of the Navier-Stokes equations. Computed solutions are obtained for an aircraft geometry which includes the fuselage, a wing with deflected leading-edge flap, empennage, and a faired-over engine inlet. The computational slot geometry corresponds to that used in full-scale wind-tunnel tests. Solutions are computed using flight test conditions and jet mass flow ratios equivalent to wind-tunnel test conditions. The effect of slot location is analyzed by computing two nontime-accurate solutions with a 16-in. slot located 3 in. and 11 in. aft of the nose of the aircraft. These computations resolve the trends observed in the full-scale wind-tunnel test data. The flow aft of the leading-edge extension vortex burst is unsteady. A time-accurate solution is obtained to investigate the flow characteristics aft of the vortex burst, including the effect of blowing on tail buffet.

Description: The mathematical formulation of Dewar container sloshing dynamics when partially filled with superfluid liquid He 11 driven by the gravity jitter acceleration associated with slew motion in a microgravity environment is studied. Numerical computation of the sloshing dynamics is based on the non-inertia container bound co-ordinate, and solving time-dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. This study discloses the capillary effect of the sloshing dynamics governing the liquid-vapour interface fluctuations driven by the gravity jitter acceleration associated with slew motion, which affects the stability of the fluid system in a microgravity environment. The peculiar behaviour of superfluid helium when reacting to the sloshing dynamics driven by the gravity jitter acceleration associated with slew motion in a microgravity environment is also investigated.

Description: Floquet eigenanalysis requires a few dominant eigenvalues of the Floquet transition matrix (FTM). Although the QR method is used almost exclusively, it is expensive for such partial eigenanalysis; the operation counts and, thereby, the approximate machine-time grow cubically with the matrix order. Accordingly, for Floquet eigenanalysis, the Arnold-Saad method, a subspace iteration method, is investigated as an alternative to the QR method. The two methods are compared for machine-time efficiency and the residual errors of the corresponding eigenpairs. The Arnolds-Saad method takes much less machine-time than the QR method with comparable computational reliability and offers promise fpr large-scale Floquet eigenanalysis.

Description: The experimental documentation of a three-dimensional shock-wave/boundary-layer interaction in a nominal Mach 3 cylinder, aligned with the free-stream flow, and 20 deg. half-angle conical flare offset 1.27 cm from the cylinder centerline. Surface oil flow, laser light sheet illumination, and schlieren were used to document the flow topology. The data includes surface-pressure and skin-friction measurements. A laser interferometric skin friction data. Included in the skin-friction data are measurements within separated regions and three-dimensional measurements in highly-swept regions. The skin-friction data will be particularly valuable in turbulence modeling and computational fluid dynamics validation.

Description: The rate of heat transfer on the tip of a turbine rotor blade and on the blade surface in the vicinity of the tip, was successfully predicted. The computations were performed with a multiblock computer code which solves the Reynolds Averaged Navier-Stokes equations using an efficient multigrid method. The case considered for the present calculations was the Space Shuttle Main Engine (SSME) high pressure fuel side turbine. The predictions of the blade tip heat transfer agreed reasonably well with the experimental measurements using the present level of grid refinement. On the tip surface, regions with high rate of heat transfer was found to exist close to the pressure side and suction side edges. Enhancement of the heat transfer was also observed on the blade surface near the tip. Further comparison of the predictions was performed with results obtained from correlations based on fully developed channel flow.

Description: A 3-D compressible Navier-Stokes solver has been developed and applied to 3-D viscous flow over clean and iced wings. This method uses a third order accurate finite volume scheme with flux difference splitting to model the inviscid fluxes, and second order accurate symmetric differences to model the viscous terms. The effects of turbulence are modeled using a Kappa-epsilon model. In the vicinity of the sold walls the kappa and epsilon values are modeled using Gorski's algebraic model. Sampling results are presented for surface pressure distributions, for untapered swept clean and iced wings made of NACA 0012 airfoil sections. The leading edge of these sections is modified using a simulated ice shape. Comparisons with experimental data are given.

Description: The design of the high lift system has a major impact on the performance of an aircraft yet our understanding of the physics of this flow is still weak. Flow features include interactions between the wakes shed from the upstream elements and the pressure gradients and boundary layers of the downstream elements. Interaction of the turbulent wake of the main element and the flap can cause (1) separation of the flap boundary layer or (2) 'bursting' of the main airfoil wake. Although the first factor is at least partially understood, even the qualitative aspects of (2) remain to be determined. In order to study these phenomena at Reynolds numbers approaching those of flight, a thick high Reynolds number wake is created using a 24 foot flat plate in the long rectangular test section of a 4 ft. by 6 ft subsonic wind tunnel. The design and construction of this test section, plate, and accompanying flap is described. Results obtained in a quarter-scale model were used for design purposes and are also described. Construction of the full scale facility is complete and preliminary results are presented.

Description: Decay of isotropic turbulence is computed using direct numerical simulations. Comparisons with experimental spectra at moderate and low Reynolds numbers (R(sub lambda) less than 70) show good agreement. At moderate to high Reynolds numbers (R(sub lambda) greater 50), the spectra are found to collapse with Kolmogorov scaling at high wave numbers. However, at low Reynolds numbers (R(sub lambda) less than 50) the shape of the spectra at the Kolmogorov length scales is Reynolds number dependent. Direct simulation data from flowfields of decaying isotropic turbulence are used to compute the terms in the equation for the dissipation rate of the turbulent kinetic energy. The development of the skewness and the net destruction of the turbulence dissipation rate in the limit of low Reynolds numbers are presented. The nonlinear terms are found to remain active at surprisingly low Reynolds numbers.

Description: The effects of convection on diffusive-convective physical vapor transport process are examined computationally. We analyze conditions ranging from typical laboratory conditions to conditions achievable only in a low gravity environment. This corresponds to thermal Rayleigh numbers Ra(sub tau) ranging from 1.80 x 10 to 1.92 x 10(exp 6). Our results indicate that the effect of the sublimation and condensation fluxes at the boundaries is to increase the threshold of instability. For typical ground based conditions, time dependent oscillatory convection can occur. This results in unsteady transport, and non-uniform temperature and concentration gradients at the crystal interface. Spectral analysis of the flow field shows parametric regions exhibiting both an oscillatory approach to steady state and a chaotic transient to a periodic state. Low gravity conditions stabilize the flow field. Convective effects are effectively reduced, thus resulting in uniform temperature and concentration gradients at the interface, a desirable condition for crystal growth.

Description: A 1.27 cm diameter two phase gas-liquid flow experiment has been developed with the NASA Lewis Research Center to study two-phase flows in microgravity. The experiment allows for the measurement of void fraction, pressure drop, film thickness and bubble and wave velocities as well as for high speed photography. Three liquids were used to study the effects of liquid viscosity and surface tension, and flow pattern maps are presented for each. The experimental results are used to develop mechanistically based models to predict void fraction, bubble velocity, pressure drop and flow pattern transitions in microgravity.

Description: An investigation of the numerical simulation with two-equation turbulence models of a three-dimensional hypersonic intersecting (SWTBL) shock-wave/turbulent boundary layer interaction flow is presented. The flows are solved with an efficient implicit upwind flux-difference split Reynolds-averaged Navier-Stokes code. Numerical results are compared with experimental data for a flow at Mach 8.28 and Reynolds number 5.3x10(exp 6) with crossing shock-waves and expansion fans generated by two lateral 15 fins located on top of a cold-wall plate. This experiment belongs to the hypersonic database for modeling validation. Simulations show the development of two primary counter-rotating cross-flow vortices and secondary turbulent structures under the main vortices and in each corner singularity inside the turbulent boundary layer. A significant loss of total pressure is produced by the complex interaction between the main vortices and the uplifted jet stream of the boundary layer. The overall agreement between computational and experimental data is generally good. The turbulence modeling corrections show improvements in the predictions of surface heat transfer distribution and an increase in the strength of the cross-flow vortices. Accurate predictions of the outflow flowfield is found to require accurate modeling of the laminar/turbulent boundary layers on the fin walls.

Description: Interest has increased recently in the thrust-producing capability of rotors at very high collective pitch angles. An early reference noted this behaviour in rotors and offered alternative models for section lift characteristics to explain it. The same phenomenon was coincidentally noted and used in a propeller code, resulting in very good correlation with static thrust data. The proposed paper will present experimental data demonstrating the pronounced persistence of thrust for propellers at increasing collective pitch angles. Comparisons with blade element/momentum theory will be made. These results are expected to point to the need to define (ultimately to explain) aerodynamic lift and drag behaviour in a rotating environment. Experimental measurements made by the U.S. Army Aeroflightdynamics Directorate at the Ames Research Center have shown that locally measured normal force coefficients along the span of a highly twisted rotor blade continue to increase at high values of collective pitch. In some cases these coefficients exceed expected values for the same type of airfoil tested under two dimensional conditions. To date no one to the authors' knowledge has defined the variation of C(n) with pitch for very high angles (to 45 deg) in a rotating environment and for a blade of reasonably high aspect ratio; however, total propeller thrust measurements support the idea that stalling does not occur in the same way as on a wing. This paper will present experimental data in the form of surface pressure distributions as well as flow visualization (microtufts) to explore the aerodynamic behavior of the rotating airfoil at high values of blade incidence. This paper also reviews experimental evidence and infers some high lift coefficient behavior from it. Comparisons between predicted thrust, utilizing modified airfoil characteristics and a blade element model, and measured thrust for both rotors and propellers that cover the extremes of collective pitch are shown and discussed.

Description: An instrument is presented for laser diagnostics measurements of droplet combustion in low gravity flight aircraft. Successful applications of pulsed laser-induced fluorescence and particle interferometric velocimetry are demonstrated; the latter is noted to be somewhat less accurate than the former due to gas thermal expansion deflections.

Description: Recent progress in both the linear and nonlinear aspects of stability theory has highlighted the importance of the receptivity problem. One of the most unclear aspects of receptivity study is the receptivity of boundary-layer flow normal to vortical disturbances. Some experimental and theoretical results permit the proposition that quasi-steady outer-flow vortical disturbances may trigger by-pass transition. In present work such interaction is investigated for vorticity normal to a leading edge. The interest in these types of vortical disturbances arise from theoretical work, where it was shown that small sinusoidal variations of upstream velocity along the spanwise direction can produce significant variations in the boundary-layer profile. In the experimental part of this work, such non-uniform flow was created and the laminar-turbulent transition in this flow was investigated. The experiment was carried out in a low-turbulence direct-flow wind tunnel T-361 at the Central Aerohydrodynamic Institute (TsAGI). The non-uniform flow was produced by laminar or turbulent wakes behind a wire placed normal to the plate upstream of the leading edge. The theoretical part of the work is devoted to studying the unstable disturbance evolution in a boundary layer with strongly non-uniform velocity profiles similar to that produced by outer-flow vorticity. Specifically, the Tollmien-Schlichting wave development in the boundary layer flow with spanwise variations of velocity is investigated.

Description: Turbulent flow of supercritical hydrogen through a uniformly heated circular tube has been investigated using numerical methods, for the range of 4 x 10(exp 5) less than Re less than 3 x 10(exp 6), 5 less than or equal to q(sub W) less than or equal to 10 MW/sq m, 30 less than or equal to T(sub in) less than or equal to 90 K, and 5 less than or equal to P(sub in) less than or equal to 15 MPa. The purpose is to validate a turbulence model and calculation method for the design of active cooling systems of hydrogen-fueled hypersonic aircraft, where the hydrogen fuel a used as coolant. The PHOENICS software package was used for the computations, which required special provision for evaluation of the thermophysical properties of the supercritical hydrogen, and a low Reynolds number form of the k-epsilon turbulence model. Pressure drop and heat transfer data were compared with experiment and existing correlations and good agreement was demonstrated. For the pressure range considered here a "thermal spike" was observed and shown to be due to the secondary peak in specific heat, rather than the primary peak.

Description: Fundamental studies of fluid mechanics and transport in partially saturated soils are presented. Solution of transient diffusion problems in support of the development of probes for the in-situ measurement of moisture content is given. Numerical and analytical methods are used to study the fundamental problem of meniscus and saturation front propagation in geometric models of porous media.

Description: A technique to obtain the sensitivity of the static aeroelastic response of a three dimensional wing model is designed and implemented. The formulation is quite general and accepts any aerodynamic and structural analysis capability. A program to combine the discipline level, or local, sensitivities into global sensitivity derivatives is developed. A variety of representations of the wing pressure field are developed and tested to determine the most accurate and efficient scheme for representing the field outside of the aerodynamic code. Chebyshev polynomials are used to globally fit the pressure field. This approach had some difficulties in representing local variations in the field, so a variety of local interpolation polynomial pressure representations are also implemented. These panel based representations use a constant pressure value, a bilinearly interpolated value. or a biquadraticallv interpolated value. The interpolation polynomial approaches do an excellent job of reducing the numerical problems of the global approach for comparable computational effort. Regardless of the pressure representation used. sensitivity and response results with excellent accuracy have been produced for large integrated quantities such as wing tip deflection and trim angle of attack. The sensitivities of such things as individual generalized displacements have been found with fair accuracy. In general, accuracy is found to be proportional to the relative size of the derivatives to the quantity itself.

Description: The present paper presents a numerical procedure to calculate supersonic and hypersonic flows using the compressible law of the wall. The turbulence models under consideration include the Launder-Reece-Rodi-Gibson Reynolds-stress model and the k-epsilon model. The models coupled with the proposed wall function technique have been tested in both separated and unseparated flows. The flows include (1) an insulated flat plate flow over a range of Mach numbers, (2) a Mach 5 flat plate flow with cold wall conditions, (3) a two dimensional supersonic compression corner flow, (4) a hypersonic flow over an axisymmetric flare, and (5) a hypersonic flow over a 2-D compression corner. Results indicate that the wall function technique gives improved predictions of skin friction and heat transfer in separated flows compared with models using wall dampers. Predictions of the extent of separation are not improved over the wall damper models except with the Reynolds-stress model for the supersonic compression corner flow case.

Description: Control law design for rotorcraft fly-by-wire systems normally attempts to decouple angular responses using fixed-gain crossfeeds. This approach can lead to poor decoupling over the frequency range of pilot inputs and increase the load on the feedback loops. In order to improve the decoupling performance, dynamic crossfeeds may be adopted. Moreover, because of the large changes that occur in rotorcraft dynamics due to small changes about the nominal design condition, especially for near-hovering flight, the crossfeed design must be 'robust'. A new low-order matching method is presented here to design robust crossfeed compensators for multi-input, multi-output (MIMO) systems. The technique identifies degrees-of-freedom that can be decoupled using crossfeeds, given an anticipated set of parameter variations for the range of flight conditions of concern. Cross-coupling is then reduced for degrees-of-freedom that can use crossfeed compensation by minimizing off-axis response magnitude average and variance. Results are presented for the analysis of pitch, roll, yaw and heave coupling of the UH-60 Black Hawk helicopter in near-hovering flight. Robust crossfeeds are designed that show significant improvement in decoupling performance and robustness over nominal, single design point, compensators. The design method and results are presented in an easily used graphical format that lends significant physical insight to the design procedure. This plant pre-compensation technique is an appropriate preliminary step to the design of robust feedback control laws for rotorcraft.

Description: The use of Computational Fluid Dynamics (CFD) in the design and analysis of high performance rocket engine pumps has increased in recent years. This increase has been aided by the activities of the Marshall Space Flight Center (MSFC) Pump Stage Technology Team (PSTT). The team's goals include assessing the accuracy and efficiency of several methodologies and then applying the appropriate methodology(s) to understand and improve the flow inside a pump. The PSTT's objectives, team membership, and past activities are discussed in Garcia1 and Garcia2. The PSTT is one of three teams that form the NASA/MSFC CFD Consortium for Applications in Propulsion Technology (McConnaughey3). The PSTT first applied CFD in the design of the baseline consortium impeller. This impeller was designed for the Space Transportation Main Engine's (STME) fuel turbopump. The STME fuel pump was designed with three impeller stages because a two-stage design was deemed to pose a high developmental risk. The PSTT used CFD to design an impeller whose performance allowed for a two-stage STME fuel pump design. The availability of this design would have lead to a reduction in parts, weight, and cost had the STME reached production. One sample of the baseline consortium impeller was manufactured and tested in a water rig. The test data showed that the impeller performance was as predicted and that a two-stage design for the STME fuel pump was possible with minimal risk. The test data also verified another CFD predicted characteristic of the design that was not desirable. The classical 'jet-wake' pattern at the impeller discharge was strengthened by two aspects of the design: by the high head coefficient necessary for the required pressure rise and by the relatively few impeller exit blades, 12, necessary to reduce manufacturing cost. This 'jet-wake pattern produces an unsteady loading on the diffuser vanes and has, in past rocket engine programs, lead to diffuser structural failure. In industrial applications, this problem is typically avoided by increasing the space between the impeller and the diffuser to allow the dissipation of this pattern and, hence, the reduction of diffuser vane unsteady loading. This approach leads to small performance losses and, more importantly in rocket engine applications, to significant increases in the pump's size and weight. This latter consideration typically makes this approach unacceptable in high performance rocket engines.

Description: Experiments are conducted in the Arizona State University Unsteady Wind Tunnel with a zero-pressure-gradient flat-plate model that has a 67:1 elliptical leading edge. Boundary-layer measurements are made of the streamwise fluctuating-velocity component in order to identify the amplified T-S waves that are forced by downstream-travelling, sound waves. Measurements are taken with circular 3-D roughness elements placed at the Branch 1 neutral stability point for the frequency under consideration, and then with the roughness element downstream of Branch 1. These roughness elements have a principal chord dimension equal to 2(lambda)(sub TS)/pi, of the T-S waves under study and are 'stacked' in order to resemble a Gaussian height distribution. Measurements taken just downstream of the roughness (with leading-edge T-S waves, surface roughness T-S waves, instrumentation sting vibrations and the Stokes wave subtracted) show the generation of 3-D-T-S waves, but not in the characteristic heart-shaped disturbance field predicted by 3-D asymptotic theory. Maximum disturbance amplitudes are found on the roughness centerline. However, some near-field characteristics predicted by numerical modelling are observed.

Description: An experiment has been performed to investigate the far-field hover acoustic characteristics of the XV-15 aircraft with advanced technology blades (ATB). An extensive, high-quality, far-field acoustics data base was obtained for a rotor tip speed range of 645-771 ft/s. A 12-microphone, 500-ft radius semicircular array combined with two aircraft headings provided acoustic data over the full 360-deg azimuth about the aircraft with a resolution of 15 deg. Altitude variations provided data from near in-plane to 45 deg below the rotor tip path plane. Acoustic directivity characteristics in the lower hemisphere are explored through pressure time histories, narrow-band spectra, and contour plots. Directivity patterns were found to vary greatly with azimuth angle, especially in the forward quadrants. Sharp positive pressure pulses typical of blade-vortex interactions were found to propagate aft of the aircraft and were most intense at 45 deg below the rotor plane. Modest overall sound pressure levels were measured near in-plane indicating that thickness noise is not a major problem for this aircraft when operating in the hover mode with ATB. Rotor tip speed reductions reduced the average overall sound pressure level (dB (0.0002 dyne/cm(exp 2)) by nearly 8 dB in-plane, and 12.6 deg below the rotor plane.

Description: The National Aeronautics and Space Administration (NASA) is conducting research with the goal of enabling safe improvements in the capacity of the nation's air transportation system. The wake-vortex upset hazard is an important factor in establishing the minimum safe spacing between aircraft during landing and take-off operations, thus impacting airport capacity. A batch simulation study was conducted to assess the sensitivity of various safe landing criteria in the development of an acceptable wake encounter boundary. A baseline six-degree-of-freedom simulation of a B737-100 airplane was modified to include a wake model and the vortex-induced forces and moments. The guidance and control input for the airplane was provided by an auto-land system. The wake strength and encounter geometry were varied. A sensitivity study was also conducted to assess the effects of encounter modeling methods and accuracy.

Description: An axisymmetric jet impinging on a flat plate was simulated using a low-Reynolds number k-epsilon model. The accuracy or the numerical calculations was validated by comparing the predicted flowfields with experimental data for various nozzle-plate distance to jet diameter ratios (H/d) and Reynolds numbers. Heat transfer predictions for various H/d showed good agreement with experimental data in the wall jet region, and the complex behavior of heat transfer for small H/d is predicted. However, predictions in the neighborhood of the stagnation point were poor for small H/d The results suggest that improved modeling in the equation for the dissipation rate of the turbulence kinetic energy is required.

Description: The large-eddy simulation of the spatial evolution of a stationary crossflow vortex packet in a three-dimensional boundary layer was performed. Although a coarse grid was used (compared to that required by a direct numerical simulation) the essential features of the disturbance evolution, such as the spanwise disturbance spreading and the vortex rollover, were captured accurately. The eddy viscosity became significant only in the late nonlinear stages of the simulation.

Description: Results of calculations obtained using the direct simulation Monte Carlo method for Mach 25 flow over a control surface are presented. The numerical simulations are for a 35-deg compression ramp at a low-density wind-tunnel test condition. Calculations obtained using both two- and three-dimensional solutions are reviewed, and a qualitative comparison is made with the oil flow pictures highlight separation and three-dimensional flow structure.

Description: Two new versions of the kappa-omega two-equation turbulence model will be presented. The new Baseline (BSL) model is designed to give results similar to those of the original kappa-omega model of Wilcox, but without its strong dependency on arbitrary freestream values. The BSL model is identical to the Wilcox model in the inner 50% of the boundary-layer but changes gradually to the standard kappa-epsilon model (in a kappa- omega formulation) towards the boundary-layer edge. The free shear layers. The second version of the model is called Shear-Stress Transport (SST) model. It is a variation of the BSL model with the additional ability to account for the transport of the principal turbulent shear stress in adverse pressure gradient boundary-layers. The model is based on Bradshaw's assumption that the principal shear-stress is proportional to the turbulent kinetic energy, which is introduced into the definition of the eddy-viscosity. Both models are tested for a large number of different flowfields. The results of the BSL model are similar to those of the original kappa-omega model, but without the undesirable freestream dependency. The predictions of the SST model are also independent of the freestream values but show better agreement with experimental data for adverse pressure gradient boundary-layer flows.

Description: An efficient approach for simultaneous aerodynamic analysis and design optimization is presented. This approach does not require the performance of many flow analyses at each design optimization step, which can be an expensive procedure. Thus, this approach brings us one step closer to meeting the challenge of incorporating computational fluid dynamic codes into gradient-based optimization techniques for aerodynamic design. An adjoint-variable method is introduced to nullify the effect of the increased number of design variables in the problem formulation. The method has been successfully tested on one-dimensional nozzle flow problems, including a sample problem with a normal shock. Implementations of the above algorithm are also presented that incorporate Newton iterations to secure a high-quality flow solution at the end of the design process. Implementations with iterative flow solvers are possible and will be required for large, multidimensional flow problems.

Description: This paper presents a procedure for computing the aeroelasticity of wing-body configurations on multiple-instruction, multiple-data (MIMD) parallel computers. In this procedure, fluids are modeled using Euler equations discretized by a finite difference method, and structures are modeled using finite element equations. The procedure is designed in such a way that each discipline can be developed and maintained independently by using a domain decomposition approach. A parallel integration scheme is used to compute aeroelastic responses by solving the coupled fluid and structural equations concurrently while keeping modularity of each discipline. The present procedure is validated by computing the aeroelastic response of a wing and comparing with experiment. Aeroelastic computations are illustrated for a High Speed Civil Transport type wing-body configuration.

Description: We present calculations and measurements on the shape of liquid metal droplets in electromagnetic levitation experiments. A normal stress balance model was developed to predict the shapes of liquid metal droplets that will be obtained in a microgravity experiment to measure the viscosity and surface tension of undercooled metals. This model was tested by calculating the droplet shapes in containerless experiments conducted to determine the surface tension of liquid metals. Inconsistencies associated with the results of a previous paper are elucidated. The computational results of the mathematical model are compared with the results of ground-based experiments for two different metals. The importance of the ratio of electromagnetic skin depth-to-droplet radius to the accuracy of the mathematical model is discussed. A planned alternate approach to modeling the shape by consideration of the entire droplet rather than only the surface is presented. As an example of an application. the influence of the shape on the splitting of the surface oscillation modes of levitated liquid metal droplets is discussed.

Description: A new numerical method is presented for the analysis of combined natural convection and radiation heat transfer with applications in many engineering situations such as materials processing, combustion and fire research. Because of the recent interest in the low gravity environment of space, attention is devoted to both 1-g and low-g applications. The two-dimensional mathematical model is represented by a set of coupled nonlinear integro-partial differential equations. Radiative exchange is formulated using the Discrete Exchange Factor method (DEF). This method considers point to point exchange and provides accurate results over a wide range of radiation parameters. Numerical results show that radiation significantly influences the flow and heat transfer in both low-g and 1-g applications. In the low-g environment, convection is weak, and radiation can easily become the dominant heat transfer mode. It is also shown that volumetric heating by radiation gives rise to an intricate cell pattern in the top heated enclosure.

Description: Wind tunnel tests have been conducted on two devices for the purpose of lift augmentation on a 60 deg delta wing at low speed. Lift, drag, pitching moment, and surface pressures were measured. Detailed flow visualization was also obtained. Both the leading edge fence and the Gurney flap are shown to increase lift. The fences and flap shift the lift curve by as much as 5 deg and 10 deg, respectively. The fences aid in trapping vortices on the upper surface, thereby increasing suction. The Gurney flap improves circulation at the trailing edge. The individual influences of both devices are roughly additive, creating high lift gain. However, the lower lift to drag ratio and the precipitation of vortex burst caused by the fences, and the nose down pitching moment created by the flap are also significant factors.

Description: In high performance boardsailing, demands on the vertical fin or "skeg" often produce "spinout" - when the skeg loses horizontal lift creating a force imbalance and causing the tail of the board to slide sideways. Richard Caldwell, RACE Technology, Inc. used NASA airfoil technology to solve this problem and formed a business based on his solution. After determining that the spinout resulted from air ventilating down the low pressure side of the underwater fin, he adapted the airfoil technology to the design of a short board skeg, which would overcome the problem and lower the drag, resulting in improved performance. He patented his RACE 145 foil section, formed his company and later returned to Langley for additional technical assistance. The company's newest product is a rigid sail that also incorporates NASA technology and has excellent performance. This company no longer exists - product is no longer in production.

Description: Interest in thin sheet flows has recently been renewed due to their potential application in space radiators. Theoretical and experimental studies of the fluid dynamics and stability of thin liquid sheet flows have been carried out in this thesis. A computer program was developed to determine the cross-sectional shape of the edge cylinder given the cross-sectional area of the edge cylinder. A stability analysis was performed on a non-planer liquid sheet. A study was conducted to determine the effects of air resistance on the sheet.