ALBERT

Description: Measurements of wing buffeting, using root strain gages, were made in the NASA Langley 0.3 m cryogenic wind tunnel to refine techniques which will be used in larger cryogenic facilities such as the United States National Transonic Facility (NTF) and the European Transonic Wind Tunnel (ETW). The questions addressed included the relative importance variations in frequency parameter and Reynolds number, the choice of model material (considering both stiffness and damping) and the effects of static aeroelastic distortion. The main series of tests was made on three half models of slender 65 deg delta wings with a sharp leading edge. The three delta wings had the same planform but widely differing bending stiffnesses and frequencies (obtained by varying both the material and the thickness of the wings). It was known that the steady flow on this configuration would be insensitive to variations in Reynolds number. On this wing at vortex breakdown the spectrum of the unsteady excitation is unusual, having a sharp peak at particular frequency parameter. Additional tests were made on one unswept half-wing of aspect ratio 1.5 with an NPL 9510 aerofoil section, known to be sensitive to variations in Reynolds number at transonic speeds. The test Mach numbers were M = 0.21 and 0.35 for the delta wings and to M = 0.30 for the unswept wing. On this wing the unsteady excitation spectrum is fairly flat (as on most wings). Hence correct representation of the frequency parameter is not particularly important.

Description: In the design of an airframe, the effect of changing the geometry on resulting computations is necessary for design optimization. The geometry is defined in terms of a series of design variables, including design variables to define the wing planform, tail, canard, pylon, and nacelle. Design optimization in this research is based on how these design variable affect the potential flow. The potential flow is computed as a function of the geometry and location of a series of panels describing the airframe, which are in turn a function of the design variables. Multipole accelerated panel methods improve the computational complexity of the problem and thus are an attractive approach. To utilize the methods in design optimization, it was necessary to define the appropriate sensitivity derivatives. The overhead incurred from finding the sensitivity derivatives in conjunction with the original computation should be small. This research developed the background for multipole-accelerated panel methods and the framework for finding sensitivity derivatives in the methods. Potential flow panel codes are commonly used for powered-lift aerodynamic predictions for three dimensional geometries. Given an airframe which has been discretized into a series of panels to define the airframe geometry, potential is computed as a function of the influence of all panels on all other panels. This is a computationally intensive problem for which efficient solutions are desired to improve the computational time and to allow greater resolution by use of more panels. One such solution is the use of hierarchical multipole methods which entail approximations of the effects of far-field terms. Hierarchical multipole methods have become prevalent in molecular dynamics and gravitational physics, and have been introduced into the fields of capacitance calculations, computational fluid dynamics, and electromagnetics. The methods utilize multipole expansions to describe the effect of bodies (i.e. particles, astrophysical bodies, panels, etc.) within a sphere on points distant from the sphere, where the influence diminishes as a function of distance. The expansions are exact with infinite series, however, for practical computations, the series are truncated and accuracy is selected based on the number of terms retained in the expansions. A hierarchical tree structure groups bodies together based on proximity to allow definition of multipole expansions for each group. The multipole expansions are then used to compute the effect of the bodies in a group on distant bodies.

Description: An evaluation of existing models for Large Eddy Simulations (LES) of incompressible turbulent flows has been completed. LES is a computation in which the large, energy-carrying structures to momentum and energy transfer is computed exactly, and only the effect of the smallest scales of turbulence is modeled. That is, the large eddies are computed and the smaller eddies are modeled. The dynamics of the largest eddies are believed to account for most of sound generation and transport properties in a turbulent flow. LES analysis is based on an observation that pressure, velocity, temperature, and other variables are the sum of their large-scale and small-scale parts. For instance, u(i) (velocity) can be written as the sum of bar-u(i) and u(i)-prime, where bar-u(i) is the large-scale and u(i)-prime is the subgrid-scale (SGS). The governing equations for large eddies in compressible flows are obtained after filtering the continuity, momentum, and energy equations, and recasting in terms of Favre averages. The filtering operation maintains only large scales. The effects of the small-scales are present in the governing equations through the SGS stress tensor tau(ij) and SGS heat flux q(i). The mathematical formulation of the Favre-averaged equations of motion for LES is complete.

Description: This project primarily dealt with integral boundary-layer solution techniques that are directly applicable to the problem of determining aerodynamic heating rates of hypersonic vehicles like X-33 in the vicinity of stagnation points, windward centerlines, and swept-wing leading edges. The analyses include effects of finite-rate gas chemistry across the boundary layer and finite-rate catalysis of atom recombination at the surface. A new approach for combining the insight afforded by integral boundary-layer analysis with comprehensive (and expensive) computational fluid dynamic (CFD) flowfield solutions of the thin-layer Navier-Stokes equations was developed. The approach extracts CFD derived quantities at the wall and at the boundary layer edge for inclusion in a post-processing boundary-layer analysis. The post-processed data base allows a designer at a workstation to ask and answer the following questions: (1) How much does the heating change if one uses a thermal protection system (TPS) with different catalytic properties than was used in the original CFD solution? (2) How does the heating change when one moves the interface of two different TPS materials with different catalytic efficiencies for the purpose of reducing vehicle weight and expense? The answer to the second question is particularly critical, because abrupt changes from low catalytic efficiency to high catalytic efficiency can lead to localized increase in heating which exceeds the usually conservative estimate provided by a fully catalytic wall assumption. A secondary issue that was addressed involves the prediction of heating levels in the vicinity of sharp corners that are transverse to or aligned with the flow. An example of the first case is heating at the edge of the COMET reentry module. An example of the second case is heating along the side edge of a deflected body flap on an SSV. The difficulty of putting grids in the vicinity of such corners with continuously varying metric coefficients causes problems in CFD predictions. A preliminary theory for prediction that says the heating at the corner is X percent of the heating N boundary-layer thicknesses inboard was developed. This will prove useful to analytically evaluate the possible benefits of rounding the edges of these configurations and defining how much rounding is sufficient.

Description: Work continued on two projects which had been started during previous years. Both projects involve calculations of the subsonic, turbulent far wake of a two-dimensional object at a Reynolds number of 1000 (based on wake momentum thickness). This flow was used as a test case for direct comparison of various turbulence models and a direct numerical simulations (DNS) of this flow were undertaken. In the turbulence model comparison studies, for any particular model tested, a unique self-similar solution was obtained far enough downstream, regardless of inlet conditions. Furthermore, different turbulence models led to different far-wake equilibrium solutions. No turbulence model could correctly predict all features of the turbulent far wake. For example, the spreading rate and turbulent shear stresses were underpredicted by all the standard models (both two-equation and full Reynolds stress models). In cases where a more correct spreading rate was achieved, it was at the expense of the turbulent kinetic energy, which was overpredicted. In general, the Algebraic Dissipation Rate Model of Gatski and Speziale, 1992, when added to any of the standard models, improved the results dramatically. Also, full Reynolds stress closure models did a much better job at predicting the shapes of both the mean and turbulence profiles, but the spreading rate was not significantly improved over that predicted by the simpler two-equation models. There are two main conclusions from these studies: First, in a comparison such as this, it is not enough to compare just one parameter, like the spreading rate. A good prediction for one parameter does not necessarily imply good predictions for all parameters in a flow. Second, since no turbulence model could correctly predict the turbulent far wake, much of the important physics of turbulent free shear flows is apparently lost by the assumptions inherent in today's methods of turbulence modeling; turbulence models must be improved. Direct simulations of this flow were begun last year in order to provide a data base through which some of the deficiencies of the existing turbulence models could be identified. Quantities such as the pressure-strain correlation, turbulent diffusion, and the dissipation rate tensor can be easily calculated from the DNS results, whereas these quantities are nearly impossible to measure experimentally. Improvements to existing turbulence models (and development of new models) require knowledge about flow quantities such as these. During this summer, diagnostics codes were written which will calculate the parameters mentioned above, along with other single-point and multi-point statistics. The DNS calculations are still in progress at the time of this writing. When these calculations are complete, the diagnostics codes will be applied so that the results can aid turbulence modelers. In addition, the results will show whether or not there exists a universal equilibrium turbulent far wake, independent of initial conditions.

Description: Over the past quarter century, the NASA Langley Research Center (LaRC) and the NASA Dryden Flight Research Center (DFRC) have played major roles in the development, demonstration, and validation of aeroservoelastic modeling, analysis, design, and testing methods. Many of their contributions resulted from their participation in wind-tunnel and flight-test programs aimed at demonstrating advanced active control concepts that interact with and/or exploit the aeroelastic characteristics of flexible structures. Other contributions are a result of their interest in identifying and solving adverse aeroservoelastic interactions that allow unique flight-test demonstrations or flight envelope clearance programs to be successfully completed. This paper provides an overview of some of the more interesting aeroservoelastic investigations conducted in the transonic dynamics tunnel (TDT) at LaRC and in flight at DFRC. Four flight-test projects are reviewed in this paper. These test projects were selected because of their contributions to the state-of-the-art in active controls technology (ACT) or because of the knowledge gained in further understanding the complex mechanisms that cause adverse aeroservoelastic interactions.

Description: Analytical investigation of dynamic stall on HAWT (horizontal-axis wind turbines) rotor loads was conducted. Dynamic stall was modeled using the Gormont approach on the MOD-2 rotor, treating the blade as a rigid body teetering about a fixed axis. Blade flapwise bending moments at station 370 were determined with and without dynamic stall for spatial variations in local wind speed due to wind shear and yaw. The predicted mean flapwise bending moments were found to be in good agreement with test results. Results obtained with and without dynamic stall showed no significant difference for the mean flapwise bending moment. The cyclic bending moments calculated with and without dynamic stall effects were substantially the same. None of the calculated cyclic loads reached the level of the cyclic loads measured on the MOD-2 using the Boeing five-minute-average technique.

Description: A coordinated effort has been underway over the past four years to elevate unstructured-grid methodology to a mature level. The goal of this endeavor is to provide a validated capability to non-expert users for performing rapid aerodynamic analysis and design of complex configurations. The Euler component of the system is well developed, and is impacting a broad spectrum of engineering needs with capabilities such as rapid grid generation and inviscid flow analysis, inverse design, interactive boundary layers, and propulsion effects. Progress is also being made in the more tenuous Navier-Stokes component of the system. A robust grid generator is under development for constructing quality thin-layer tetrahedral grids, along with a companion Navier-Stokes flow solver. This paper presents an overview of this effort, along with a perspective on the present and future status of the methodology.

Description: The National Aeronautics and Space Administration (NASA) is conducting research to improve airport capacity by reducing the separation distance between aircraft. The limiting factor in reducing separation distances and improving airport capacity is the wake vortex hazard. The ability to accurately model wake vortices and predict the outcome of a vortex encounter is critical in developing a system to safely improve airport capacity. This is the focus of the wake vortex research being done at NASA Langley Research Center (LaRC). This paper will concentrate on two topics. The first topic is the control system developed for the Boeing 737 freeflight model in support of vortex encounter tests to be conducted in the 30- by 60- foot tunnel at NASA Langley Research Center later this year. The second topic discussed is the limited degree of freedom (DOF) trajectory generation study that is being conducted to determine the relative severity of a multitude of paths through a wake vortex.

Description: A thin-layer Navier-Stokes code, CFL3D, was utilized to compute the flow over a high-lift multi-element airfoil. This study was conducted to improve the prediction of high-lift flowfields using various turbulence models and improved glidding techniques. An overset Chimera grid system is used to model the three element airfoil geometry. The effects of wind tunnel wall modeling, changes to the grid density and distribution, and embedded grids are discussed. Computed pressure and lift coefficients using Spalart-Allmaras, Baldwin-Barth, and Menter's kappa-omega - Shear Stress Transport (SST) turbulence models are compared with experimental data. The ability of CFL3D to predict the effects on lift coefficient due to changes in Reynolds number changes is also discussed.

Description: The Federal Aviation Administration (FAA) and NASA have initiated a joint study in the development of reliable means of tracking, detecting, measuring, and predicting trailing wake-vortices of commercial aircraft. Being sought is an accurate model of the wake-vortex hazard, sufficient to increase airport capacity by reducing minimum safe spacings between planes. Several means of measurement are being evaluated for application to wake-vortex detection and tracking, including Doppler RADAR (Radio Detection and Ranging) systems, 2-micron Doppler LIDAR (Light Detection And Ranging) systems, and SODAR (Sound Detection And Ranging) systems. Of specific interest there is the lidar system, which has demonstrated numerous valuable capabilities as a vortex sensor Aerosols entrained in the vortex flow make the wake velocity signature visible to the lidar, (the observable lidar signal is essentially a measurement of the line-of-sight velocity of the aerosols). Measurement of the occurrence of a wake vortex requires effective reception and monitoring of the beat signal which results from the frequency-offset between the transmitted pulse and the backscattered radiation. This paper discusses the mounting, analysis, troubleshooting, and possible use of an analog processing assembly designed for such an application.

Description: As both computers and flow-analyzing equations have increased in sophistication, Computational Fluid Dynamics (CFD) has evolved into a fixture for advanced aircraft design. While CFD codes have improved in accuracy and efficiency, their ability to encompass viscous effects is lacking in certain areas. For example, current CFD codes cannot accurately predict or correct for the increased drag due to these viscous effects at some flow conditions. However, by analyzing an airfoil's turbulent boundary layer, one can predict not only flow separation via the shape factor parameter, but also viscous drag via the momentum thickness. Various codes have been written which can calculate turbulent boundary layer parameters. The goal of my research is to develop procedures for modifying an airfoil (via its local pressure distribution) to eliminate boundary layer separation and/or to reduce viscous drag. The modifications to the local pressure distribution necessary to achieve these objectives will be determined using a direct-iterative method installed into a turbulent boundary layer analyzer. Furthermore, the modifications should preserve the basic characteristics of the original airfoil.

Description: The aim of this project is to develop and implement a sniffer that is capable of measuring the mass flow rate of air through a small area of pinholes whose diameters are on the magnitude of thousandths of an inch. The sniffer is used to scan a strip of a leading edge panel, which is being used in a hybrid laminar flow control experiment, in order to survey the variations in the amount of air that passes through the porous surface at different locations. Spanwise scans are taken at different chord locations by increasing the pressure in a control volume that is connected to the sniffer head, and recording the drop in pressure as the air is allowed to flow through the tiny holes. This information is used to obtain the mass flow through the structure. More importantly, the deviations from the mean flow rate are found and used to determine whether there are any significant variations in the flow rate from one area to the next. The preliminary results show little deviation in the spanwise direction. These results are important when dealing with the location and amount of suction that will be applied to the leading edge in the active laminar flow control experiment.

Description: The conversion of the Aerodynamic Preliminary Analysis System (APAS) software from a Silicon Graphics UNIX-based platform to a DOS-based IBM PC compatible is discussed. Relevant background information is given, followed by a discussion of the steps taken to accomplish the conversion and a discussion of the type of problems encountered during the conversion. A brief comparison of aerodynamic data obtained using APAS with data from another source is also made.

Description: A computer program that models the Lidar return signal for Wake Vortex experiments conducted by the Aerosol Research Branch was written. The specifications of the program and basic theory behind the calculations are briefly discussed. Results of the research and possible future improvements on it are also discussed.

Description: A survey to determine the characteristics of a boundary layer that forms on the wall of the Subsonic Basic Research Tunnel has been performed. Early results showed significant differences in the velocity profiles as measured spanwise across the wall. An investigation of the flow in the upstream contraction revealed the presence of a separation bubble at the beginning of the contraction which caused much of the observed unsteadiness. Vortex generators were successfully applied to the contraction inlet to alleviate the separation. A final survey of the wall boundary layer revealed variations in the displacement and momentum thicknesses to be less than +/- 5% for all but the most upper portion of the wall. The flow quality was deemed adequate to continue the planned follow-on tests to help develop the semi-span test technique.

Description: This report accounts details of two research projects for the Langley Aerospace Research Summer Scholars (LARSS) program. The first project, with the Office of Mission Assurance, involved subjectively predicting the probable success of two aeronautics programs by means of a tool called a Figure of Merit. The figure of merit bases program success on the quality and reliability of the following factors: parts, complexity of research, quality programs, hazards elimination, and single point failures elimination. The second project, for the Office of Safety and Facilities Assurance, required planning, layouts, and source seeking for an addition to the fire house. Forecasted changes in facility layout necessitate this addition which will serve as housing for the fire fighters.

Description: Two contrasting models of gas-surface interactions are studied using the Direct Simulation Monte Carlo (DSMC) method. The DSMC calculations examine differences in predictions of aerodynamic forces and heat transfer between the Maxwell and Cercignani-Lampis-Lord (CLL) models for flat plate configurations at freestream conditions corresponding to a 140 km orbit around Venus. The size of the flat plate is that of one of the solar panels on the Magellan spacecraft, and the freestream conditions are one of those experienced during aerobraking maneuvers. Results are presented for both a single flat plate and a two-plate configuration as a function of angle of attack and gas-surface accommodation coefficients. The two plate system is not representative of the Magellan geometry, but is studied to explore possible experiments that might be used to differentiate between the two gas surface interaction models.

Description: The wind tunnel investigation into the acoustic nature of boundary layer transition using miniature microphones. This research is the groundwork for entry into the National Transonic Facility (NTF) at the NASA Langley Research Center (LaRC). Due to the extreme environmental conditions of NTF testing, low temperatures and high pressures, traditional boundary layer detection methods are not available. The emphasis of this project and further studies is acoustical sampling of a typical boundary layer and environmental durability of the miniature microphones. The research was conducted with the 14 by 22 Foot Subsonic Tunnel, concurrent with another wind tunnel test. Using the resources of LaRC, a full inquiry into the feasibility of using Knowles Electronics, Inc. EM-3086 microphones to detect the surface boundary layer, under differing conditions, was completed. This report shall discuss the difficulties encountered, product performance and observations, and future research adaptability of this method.

Description: In investigating the potential of a new actuator for use in an active flow control system, several objectives had to be accomplished, the largest of which was the experimental setup. The work was conducted at the NASA Langley 20x28 Shear Flow Control Tunnel. The actuator named Thunder, is a high deflection piezo device recently developed at Langley Research Center. This research involved setting up the instrumentation, the lighting, the smoke, and the recording devices. The instrumentation was automated by means of a Power Macintosh running LabVIEW, a graphical instrumentation package developed by National Instruments. Routines were written to allow the tunnel conditions to be determined at a given instant at the push of a button. This included determination of tunnel pressures, speed, density, temperature, and viscosity. Other aspects of the experimental equipment included the set up of a CCD video camera with a video frame grabber, monitor, and VCR to capture the motion. A strobe light was used to highlight the smoke that was used to visualize the flow. Additional effort was put into creating a scale drawing of another tunnel on site and a limited literature search in the area of active flow control.

Description: The coding of high-performance fluids simulations requires significant knowledge of both numerical and computational details. The magnitude and complexity of low-level details is often enough to discourage many users of turbulence data wishing to study more important, higher-level fluid dynamical questions. These same complexities are often a practical barrier to simulation experts who develop, verify and maintain the codes which generate this data. Future fluids codes, with high resolution and complex geometries, are likely to involve far more coding complexity. The research--the design and implementation of the Tensoral computer language--aims to greatly ease the coding of today's simulation and post-processing codes and at the same time provide a general computational tool for future simulations.

Description: This two-month visit at CTR was devoted to investigating possibilities in LES modeling in the context of the 3-D vortex particle method (=vortex element method, VEM) for unbounded flows. A dedicated code was developed for that purpose. Although O(N(sup 2)) and thus slow, it offers the advantage that it can easily be modified to try out many ideas on problems involving up to N approx. 10(exp 4) particles. Energy spectrums (which require O(N(sup 2)) operations per wavenumber) are also computed. Progress was realized in the following areas: particle redistribution schemes, relaxation schemes to maintain the solenoidal condition on the particle vorticity field, simple LES models and their VEM extension, possible new avenues in LES. Model problems that involve strong interaction between vortex tubes were computed, together with diagnostics: total vorticity, linear and angular impulse, energy and energy spectrum, enstrophy. More work is needed, however, especially regarding relaxation schemes and further validation and development of LES models for VEM. Finally, what works well will eventually have to be incorporated into the fast parallel tree code.

Description: The objective is the development of a high-order and high-resolution method for the direct numerical simulation of shock turbulent-boundary-layer interaction. Details concerning the spatial discretization of the convective terms can be found in Adams and Shariff (1995). The computer code based on this method as introduced in Adams (1994) was formulated in Cartesian coordinates and thus has been limited to simple rectangular domains. For more general two-dimensional geometries, as a compression corner, an extension to generalized coordinates is necessary. To keep the requirements or limitations for grid generation low, the extended formulation should allow for non-orthogonal grids. Still, for simplicity and cost efficiency, periodicity can be assumed in one cross-flow direction. For easy vectorization, the compact-ENO coupling algorithm as used in Adams (1994) treated whole planes normal to the derivative direction with the ENO scheme whenever at least one point of this plane satisfied the detection criterion. This is apparently too restrictive for more general geometries and more complex shock patterns. Here we introduce a localized compact-ENO coupling algorithm, which is efficient as long as the overall number of grid points treated by the ENO scheme is small compared to the total number of grid points. Validation and test computations with the final code are performed to assess the efficiency and suitability of the computer code for the problems of interest. We define a set of parameters where a direct numerical simulation of a turbulent boundary layer along a compression corner with reasonably fine resolution is affordable.

Description: We are developing computational tools for the simulations of three-dimensional flows past bodies undergoing arbitrary motions. High resolution viscous vortex methods have been developed that allow for extended simulations of two-dimensional configurations such as vortex generators. Our objective is to extend this methodology to three dimensions and develop a robust computational scheme for the simulation of such flows. A fundamental issue in the use of vortex methods is the ability of employing efficiently large numbers of computational elements to resolve the large range of scales that exist in complex flows. The traditional cost of the method scales as Omicron (N(sup 2)) as the N computational elements/particles induce velocities at each other, making the method unacceptable for simulations involving more than a few tens of thousands of particles. In the last decade fast methods have been developed that have operation counts of Omicron (N log N) or Omicron (N) (referred to as BH and GR respectively) depending on the details of the algorithm. These methods are based on the observation that the effect of a cluster of particles at a certain distance may be approximated by a finite series expansion. In order to exploit this observation we need to decompose the element population spatially into clusters of particles and build a hierarchy of clusters (a tree data structure) - smaller neighboring clusters combine to form a cluster of the next size up in the hierarchy and so on. This hierarchy of clusters allows one to determine efficiently when the approximation is valid. This algorithm is an N-body solver that appears in many fields of engineering and science. Some examples of its diverse use are in astrophysics, molecular dynamics, micro-magnetics, boundary element simulations of electromagnetic problems, and computer animation. More recently these N-body solvers have been implemented and applied in simulations involving vortex methods. Koumoutsakos and Leonard (1995) implemented the GR scheme in two dimensions for vector computer architectures allowing for simulations of bluff body flows using millions of particles. Winckelmans presented three-dimensional, viscous simulations of interacting vortex rings, using vortons and an implementation of a BH scheme for parallel computer architectures. Bhatt presented a vortex filament method to perform inviscid vortex ring interactions, with an alternative implementation of a BH scheme for a Connection Machine parallel computer architecture.

Description: We are interested in the study, via direct numerical simulations, of active vortex generators. Vortex generators may be used to modify the inner part of the boundary layer or to control separation thus enhancing the performance and maneuverability of aerodynamic configurations. We consider generators that consist of a surface cavity elongated in the stream direction and partially covered with a moving lid that at rest lies flush with the boundary. Streamwise vorticity is generated and ejected due to the oscillatory motion of the lid. The present simulations complement relevant experimental investigations of active vortex generators at NASA Ames and Stanford University (Saddoughi, 1994, and Jacobson and Reynolds, 1993). Jacobson and Reynolds (1993) used a piezoelectric device in water, allowing for small amplitude high frequency oscillations. They placed the lid asymmetrically on the cavity and observed a strong outward velocity at the small gap of the cavity. Saddoughi used a larger mechanically driven device in air to investigate this flow and he observed a jet emerging from the wide gap of the configuration, contrary to the findings of Jacobson and Reynolds. Our task is to simulate the flows generated by these devices and to conduct a parametric study that would help us elucidate the physical mechanisms present in the flow. Conventional computational schemes encounter difficulties when simulating flows around complex configurations undergoing arbitrary motions. Here we present a formulation that achieves this task on a purely Lagrangian frame by extending the formulation presented by Koumoutsakos, Leonard and Pepin (1994). The viscous effects are taken into account by modifying the strength of the particles, whereas fast multipole schemes employing hundreds of thousands of particles allow for high resolution simulations. The results of the present simulations would help us assess some of the effects of three-dimensionality in experiments and investigate the role of two-dimensional vortex generation due to an oscillating lid.

Description: This is a report on the continuation of our experimental investigations (Saddoughi 1994) of 'on-demand' vortex generators. Conventional vortex generators as found on aircraft wings are mainly for suppression of separation during the off-design conditions. In cruise they perform no useful function and exert a significant drag penalty. Therefore, replacement of fixed rectangular or delta-wing generators by devices that could be activated when needed would be of interest. Also in our previous report, we described one example of an 'on-demand' device, which was developed by Jacobson & Reynolds (1995) at Stanford University, suitable for manufacture by micro-electro-mechanical technology. This device consists of a surface cavity elongated in the stream direction and covered with a lid cantilevered at the upstream end. The lid, which is a metal sheet with a sheet of piezoelectric ceramic bonded to it, lies flush with the boundary. On application of a voltage the ceramic expands or contracts; however, adequate amplitude can be obtained only by running at the cantilever resonance frequency and applying amplitude modulation: for 2.5 mm x 20 mm cantilevered lids, they obtained maximum tip displacements of the order of 100 pm. Thus fluid is expelled from the cavity through the gap around the lid on the downstroke. They used an asymmetrical gap configuration and found that periodic emerging jets on the narrow side induced periodic longitudinal vorticity into the boundary layer. Their device was used to modify the inner layer of the boundary layer for skin-friction reduction. The same method could be implemented for the replacement of the conventional vortex generators; however, to promote mixing and suppress separation we needed to deposit longitudinal vortices into the outer layer of the boundary layer, which required a larger vortex generator than the device built by Jacobson & Reynolds. Our vortex generator was built with a mechanically-driven cantilevered lid with an adjustable frequency. The device was made about ten times the size of Jacobson & Reynolds', the shape or size of the cavity and lid (28 mm x 250 mm) could be easily changed. The cavity depth, the cantilever-tip displacement, and the maximum lid frequency were 20 mm, 10 mm, and 60 Hz respectively. Our vortex generator was mounted on a turntable so that its yaw angle could be changed. Finally, tests over a range of ratios of vortex generator size to boundary-layer thickness could be carried out simply by changing the streamwise location of the device.

Description: Our understanding of aerodynamic noise has its foundations in the work of Sir James Lighthill (1952), which was the first major advance in acoustics since the pioneering work of Lord Rayleigh in the last century. The combination of Lighthill's theory of aerodynamic noise as applied to turbulent flows and the experimental growing database from the early 1950's was quickly exploited by various jet propulsion engine designers in reducing the noise of jet engines at takeoff and landing to levels marginally acceptable to communities living in the neighborhoods of airports. The success in this noise containment led to the rapid growth of fast economical subsonic civil transport aircraft worldwide throughout the 1960's and has continued to the present day. One important factor in this success story has been the improvements in the engine cycle that have led to both reductions in specific fuel consumption and noise. The second is the introduction of Noise Certification, which specifies the maximum noise levels at takeoff and landing that all aircraft must meet before they can be entered on the Civil Aircraft Register. The growing interest in the development of a new supersonic civil transport to replace 'Concorde' in the early years of the next century has led to a resurgence of interest in the more challenging problem of predicting the noise of hot supersonic jets and developing means of aircraft noise reduction at takeoff and landing to meet the standards now accepted for subsonic Noise Certification. The prediction of aircraft noise to the accuracy required to meet Noise Certification requirements has necessitated reliance upon experimental measurements and empirically derived laws based on the available experimental data bases. These laws have their foundation in the results from Lighthill's theory, but in the case of jet noise, where the noise is generated in the turbulent mixing region with the external ambient fluid, the complexity of the turbulent motion has prevented the full deployment of Lighthill's theory from being achieved. However, the growth of the supercomputer and its applications in the study of the structure of turbulent shear flows in both unbounded and wall bounded flows, which complements and in certain cases extends the work of the few dedicated experimental groups working in this field for the past forty years, provides an opportunity and challenge to accurately predict the noise from jets. Moreover a combination of numerical and laboratory experiments offers the hope that in the not too distant future the physics of noise generation and flow interaction will be better understood and it will then be possible to not only improve the accuracy of noise prediction but also to explore and optimize schemes for noise reduction. The present challenge is to provide time and space accurate numerical databases for heated subsonic and supersonic jets to provide information on the fourth-order space-time covariance of Lighthill's equivalent stress tensor, T(ij), which governs the characteristics of the farfield radiated noise and the total acoustic power. Validation with available experimental databases will establish how close Lighthill's theory is to the accurate prediction of the directivity and spectrum of jet noise and the total acoustic power, and the need, in the applications of the theory, to include the effects of flow-acoustic interaction.

Description: This note describes an attempt to use standard neural network tools to fashion a means of detecting eddy patterns in the wall region of a turbulent flow. The research was motivated by the desire to formulate a means to use only flow parameters that can be sensed on the wall to describe the passing eddy structure. If a simple formulation can be obtained, it could conceivably be utilized to control actuators embedded in the wall. Such actuators have been developed by Jacobson and Reynolds (1993a), Blackwelder and Liu (1994), Tung et al. (1995), and others. These actuators have the common characteristics that they are small and are typically flush with the wall when not deployed. When they are activated, it is assumed that they will be able to interact constructively with the turbulent eddies near their location to either decrease the wall shear stress, enhance or reduce the mixing, etc. At present, there is only a nascent understanding of the interaction dynamics between the actuators and the eddies in the flow. Nevertheless, for such interaction to succeed, methods to couple the actuators to the oncoming flow must be obtained. General methods must be found that will detect the space and temporal location of the desired structure. In particular, it will be necessary to know when the eddies will arrive at the location of the actuator. This research attempted to use the shear stress measurements on the wall in the vicinity of an actuator location to predict when a particular eddy pattern would arrive and/or occur at the designated location. In this work the eddy pattern to be detected was identified by its velocity signature only.

Description: The general objective of the present work is to explore the use of Rapid Distortion Theory (RDT) in analysis of the two-point statistics of the log-layer. RDT is applicable only to unsteady flows where the non-linear turbulence-turbulence interaction can be neglected in comparison to linear turbulence-mean interactions. Here we propose to use RDT to examine the structure of the large energy-containing scales and their interaction with the mean flow in the log-region. The contents of the work are twofold: First, two-point analysis methods will be used to derive the law-of-the-wall for the special case of zero mean pressure gradient. The basic assumptions needed are one-dimensionality in the mean flow and homogeneity of the fluctuations. It will be shown that a formal solution of the two-point correlation equation can be obtained as a power series in the von Karman constant, known to be on the order of 0.4. In the second part, a detailed analysis of the two-point correlation function in the log-layer will be given. The fundamental set of equations and a functional relation for the two-point correlation function will be derived. An asymptotic expansion procedure will be used in the log-layer to match Kolmogorov's universal range and the one-point correlations to the inviscid outer region valid for large correlation distances.

Description: The objective of this effort is to carry the analysis of Lee et al. (1990) to the case of shear with rotation. We apply the RDT approximation to turbulence submitted to frame rotation for the case of a uniformly sheared flow and compare its mean statistics to results of high resolution DNS of a rotating plane channel flow. In the latter, the mean velocity profile is modified by the Coriolis force, and accordingly, different regions in the channel can be identified. The properties of the plane pure strain turbulence submitted to frame rotation are, in addition, investigated in spectral space, which shows the usefulness of the spectral RDT approach. This latter case is investigated here. Among the general class of quadratic flows, this case does not follow the same stability properties as the others since the related mean vorticity is zero.

Description: Finite-difference second-order accurate direct simulation of a turbulent pipe has been used to investigate how the turbulence production and dissipation change when a solid body rotation is applied. It is shown that when the helicity increases, the dissipation is reduced. It is asserted that to have a drag reduction the external action should be such as to disrupt the symmetry of right- and left-handed helical structures. In this study the Navier-Stokes equations in rotational form permit the turbulent energy production to be split into a part related to the energy cascade from large to small scales and into a part related to the convection by large scales. The full simulation data have shown the latter is greater than the former in the wall region and that, on the contrary, these two terms balance each other in the central region. From the pdf of the former, it has been shown how the vortical structures are changed in the wall region by the background radiation and how they are related to the changes in the energy production.

Description: This report shows that if a stochastic differential equation (Langevin equation) for velocity fluctuation vector is known, it is possible to derive the equations for scalar flux transport. Durbin and Speziale (1994) showed that the second moment of this stochastic differential equation gives an equation for the evolution of Reynolds stress tensor. Similarly, the stochastic equation will give an equation for scalar flux. Therefore, a coupling between these two is present. The basis for the present work is that there should be Langevin equations that can produce acceptable models for both the Reynolds stress tensor and the scalar flux vector. Having found this basic Langevin equation, the amount of work needed to model the second order closure problems is reduced; using the well developed models for Reynolds stress equations, it will be possible to derive corresponding models for scalar flux equation.

Description: The analysis demonstrated that the governing equations for the two-point velocity correlation tensor in the temporally evolving wake admit similarity solutions, which include the similarity solutions for the single-point moment as a special case. The resulting equations for the similarity solutions include two constants, beta and Re(sub sigma), that are ratios of three characteristic time scales of processes in the flow: a viscous time scale, a time scale characteristic of the spread rate of the flow, and a characteristic time scale of the mean strain rate. The values of these ratios depend on the initial conditions of the flow and are most likely measures of the coherent structures in the initial conditions. The occurrences of these constants in the governing equations for the similarity solutions indicates that these solutions, in general, will only be the same for two flows if these two constants are equal (and hence the coherent structures in the flows are related). The comparisons between the predictions of the similarity hypothesis and the data presented here and elsewhere indicate that the similarity solutions for the two-point correlation tensors provide a good approximation of the measures of those motions that are not significantly affected by the boundary conditions caused by the finite extent of real flows. Thus, the two-point similarity hypothesis provides a useful tool for both numerical and physical experimentalist that can be used to examine how the finite extent of real flows affect the evolution of the different scales of motion in the flow.

Description: The work performed during this year has involved further assessment and extension of the k-epsilon-v(exp 2) model, and initiation of work on scalar transport. The latter is introduced by the contribution of Y. Shabany to this volume. Flexible, computationally tractable models are needed for engineering CFD. As computational technology has progressed, the ability and need to use elaborate turbulence closure models has increased. The objective of our work is to explore and develop new analytical frameworks that might extend the applicability of the modeling techniques. In past years the development of a method for near-wall modeling was described. The method has been implemented into a CFD code and its viability has been demonstrated by various test cases. Further tests are reported herein. Non-equilibrium near-wall models are needed for some heat transfer applications. Scalar transport seems generally to be more sensitive to non-equilibrium effects than is momentum transport. For some applications turbulence anisotropy plays a role and an estimate of the full Reynolds stress tensor is needed. We have begun work on scalar transport per se, but in this brief I will only report on an extension of the k-epsilon-v(exp 2) model to predict the Reynolds stress tensor.

Description: We have proposed and implemented an extension of the structure-based model for weak deformations. It was shown that the extended model will correctly reduce to the form of standard k-e models for the case of equilibrium under weak mean strain. The realizability of the extended model is guaranteed by the method of its construction. The predictions of the proposed model were very good for rotating homogeneous shear flows and for irrotational axisymmetric contraction, but were seriously deficient in the case of plane strain and axisymmetric expansion. We have concluded that the problem behind these difficulties lies in the algebraic constitutive equation relating the Reynolds stresses to the structure parameters rather than in the slow model developed here. In its present form, this equation assumes that under irrotational strain the principal axes of the Reynolds stresses remain locked onto those of the eddy-axis tensor. This is correct in the RDT limit, but inappropriate under weaker mean strains, when the non-linear eddy-eddy interactions tend to misalign the two sets of principal axes and create some non-zero theta and gamma.

Description: The purpose of this research is to construct accurate finite difference schemes for incompressible unsteady flow simulations such as LES (large-eddy simulation) or DNS (direct numerical simulation). In this report, conservation properties of the continuity, momentum, and kinetic energy equations for incompressible flow are specified as analytical requirements for a proper set of discretized equations. Existing finite difference schemes in staggered grid systems are checked for satisfaction of the requirements. Proper higher order accurate finite difference schemes in a staggered grid system are then proposed. Plane channel flow is simulated using the proposed fourth order accurate finite difference scheme and the results compared with those of the second order accurate Harlow and Welch algorithm.

Description: There is a need for experimental measurements in complex turbulent flows that originate from very well-defined initial conditions. Testing of large-eddy simulations and other higher-order computation schemes requires inlet boundary condition data that are not normally measured. The use of fully developed upstream conditions offers a solution to this dilemma so that the upstream conditions can be adequately computed at any level of sophistication. The plane diffuser experiment by Obi et al. (1993) has received a lot of attention because it has fully-developed inlet conditions and it includes separation from a smooth wall, subsequent reattachment and redevelopment of the downstream boundary layer. The objective of this study is to provide careful qualification and detailed measurements in a recreation of the Obi experiment. The work will include extensive documentation of the flow two-dimensionality and detailed measurements required for testing of flow computations.

Description: Some of the most challenging applications of large-eddy simulation are those in complex geometries where spectral methods are of limited use. For such applications more conventional methods such as finite difference or finite element have to be used. However, it has become clear in recent years that dissipative numerical schemes which are routinely used in viscous flow simulations are not good candidates for use in LES of turbulent flows. Except in cases where the flow is extremely well resolved, it has been found that upwind schemes tend to damp out a significant portion of the small scales that can be resolved on the grid. Furthermore, it has been found that even specially designed higher-order upwind schemes that have been used successfully in the direct numerical simulation of turbulent flows produce too much dissipation when used in conjunction with large-eddy simulation. The objective of the current study is to perform a LES of incompressible flow past a circular cylinder at a Reynolds number of 3900 using a solver which employs an energy-conservative second-order central difference scheme for spatial discretization and compare the results obtained with those of Beaudan & Moin (1994) and with the experiments in order to assess the performance of the central scheme for this relatively complex geometry.

Description: To date, most large-eddy simulations (LES) have been carried out with eddy viscosity subgrid scale (SGS) models, with only a few exceptions that used the mixed model. Even though the assumptions behind Smagorinsky's model are rather stringent, it has been applied successfully to a variety of turbulent flows. This success is attributed to the ability of eddy viscosity models to drain energy from large scales, thus simulating the dissipative nature of turbulence. Most SGS models are absolutely dissipative, i.e. they remove energy from the large scales at every instant. However, SGS stresses may transfer energy back to the large scales intermittently; this reverse transfer or backscatter is especially important in geophysical flows and in transition. In a fully developed channel flow, there is reverse flow of energy from small to large scales near the walls, but eddy viscosity models are unable to account for this important feature. The dynamic localization eddy viscosity model of Ghosal et al. (1995) allows backscatter by co-evolving an auxiliary equation for the SGS energy; however, the computational cost is considerably larger than for conventional SGS models (Cabot 1994). In this report, a new non-eddy viscosity model based on local approximation of total quantities in terms of filtered ones is introduced; the scale similarity model of Bardina (1983) is a special case of this model. This procedure does not require the assumption of homogeneity, permits backscatter of energy from small to large scales, and is readily implemented in finite difference codes. The results of applying the proposed model to second order finite volume simulation of plane channel flow at high Reynolds numbers (Re(sub b) = 38,000) is described in this report. Greater emphasis is placed on the high Reynolds number flow since it provides a more rigorous test of the SGS model and its potential application. The results are compared to ones produced by the conventional and dynamic Smagorinsky models and the spectral LES of Piomelli (1993).

Description: The equations for large-eddy simulation (LES) are derived formally by applying a spatial filter to the Navier-Stokes equations. The filter width as well as the details of the filter shape are free parameters in LES, and these can be used both to control the effective resolution of the simulation and to establish the relative importance of different portions of the resolved spectrum. An analogous, but less well justified, approach to filtering is more or less universally used in conjunction with LES using finite-difference methods. In this approach, the finite support provided by the computational mesh as well as the wavenumber-dependent truncation errors associated with the finite-difference operators are assumed to define the filter operation. This approach has the advantage that it is also 'automatic' in the sense that no explicit filtering: operations need to be performed. While it is certainly convenient to avoid the explicit filtering operation, there are some practical considerations associated with finite-difference methods that favor the use of an explicit filter. Foremost among these considerations is the issue of truncation error. All finite-difference approximations have an associated truncation error that increases with increasing wavenumber. These errors can be quite severe for the smallest resolved scales, and these errors will interfere with the dynamics of the small eddies if no corrective action is taken. Years of experience at CTR with a second-order finite-difference scheme for high Reynolds number LES has repeatedly indicated that truncation errors must be minimized in order to obtain acceptable simulation results. While the potential advantages of explicit filtering are rather clear, there is a significant cost associated with its implementation. In particular, explicit filtering reduces the effective resolution of the simulation compared with that afforded by the mesh. The resolution requirements for LES are usually set by the need to capture most of the energy-containing eddies, and if explicit filtering is used, the mesh must be enlarged so that these motions are passed by the filter. Given the high cost of explicit filtering, the following interesting question arises. Since the mesh must be expanded in order to perform the explicit filter, might it be better to take advantage of the increased resolution and simply perform an unfiltered simulation on the larger mesh? The cost of the two approaches is roughly the same, but the philosophy is rather different. In the filtered simulation, resolution is sacrificed in order to minimize the various forms of numerical error. In the unfiltered simulation, the errors are left intact, but they are concentrated at very small scales that could be dynamically unimportant from a LES perspective. Very little is known about this tradeoff and the objective of this work is to study this relationship in high Reynolds number channel flow simulations using a second-order finite-difference method.

Description: The diversity of flow characteristics encountered in a flow over an airfoil near maximum lift taxes the presently available statistical turbulence models. This work describes our first attempt to apply the technique of large-eddy simulation to a flow of aeronautical interest. The challenge for this simulation comes from the high Reynolds number of the flow as well as the variety of flow regimes encountered, including a thin laminar boundary layer at the nose, transition, boundary layer growth under adverse pressure gradient, incipient separation near the trailing edge, and merging of two shear layers at the trailing edge. The flow configuration chosen is a NACA 4412 airfoil near maximum lift. The corresponding angle of attack was determined independently by Wadcock (1987) and Hastings & Williams (1984, 1987) to be close to 12 deg. The simulation matches the chord Reynolds number U(sub infinity)c/v = 1.64 x 10(exp 6) of Wadcock's experiment.

Description: Large-eddy simulation (LES) has matured to the point where application to complex flows is desirable. The extension to higher Reynolds numbers leads to an impractical number of grid points with existing structured-grid methods. Furthermore, most real world flows are rather difficult to represent geometrically with structured grids. Unstructured-grid methods offer a release from both of these constraints. However, just as it took many years for structured-grid methods to be well understood and reliable tools for LES, unstructured-grid methods must be carefully studied before we can expect them to attain their full potential. In the past two years, important building blocks have been put into place making possible a careful study of LES on unstructured grids. The first building block was an efficient mesh generator which allowed the placement of points according to smooth variation of physical length scales. This variation of length scales is in all three directions independently, which allows a large reduction in points when compared to structured-grid methods, which can only vary length scales in one direction at a time. The second building block was the development of a dynamic model appropriate for unstructured grids. The principle obstacle was the development of an unstructured-grid filtering operator. In the past year, some of the new filters developed by Jansen have been implemented into a highly parallelized finite element code based on the Galerkin/least-squares finite element method. We have chosen the NACA 4412 airfoil at maximum lift as the first simulation for a variety of reasons. First, it is a problem of significant interest since it would be the first LES of an aircraft component. Second, this flow has been the subject of three experimental studies. The third reason for considering this flow is the variety of flow features which provide an important test of the dynamic model. Only the dynamic model can be expected to perform satisfactorily in this variety of situations: from the laminar regions where it must not modify the flow at all to the turbulent boundary layers and wake where it must represent a wide variety of subgrid-scale structures. The flow configuration we have chosen is that of Wadcock (1987) at Reynolds number based on chord Re(sub c) = u(sub infinity)c/v = 1.64 x 10(exp 6), Mach number M = 0.2, and 12 deg angle of attack.

Description: The near-wall viscous and buffer regions of wall-bounded flows generally require a large expenditure of computational resources to be resolved adequately, even in large-eddy simulation (LES). Often as much as 50% of the grid points in a computational domain are devoted to these regions. The dense grids that this implies also generally require small time steps for numerical stability and/or accuracy. It is commonly assumed that the inner wall layers are near equilibrium, so that the standard logarithmic law can be applied as the boundary condition for the wall stress well away from the wall, for example, in the logarithmic region, obviating the need to expend large amounts of grid points and computational time in this region. This approach is commonly employed in LES of planetary boundary layers, and it has also been used for some simple engineering flows. In order to calculate accurately a wall-bounded flow with coarse wall resolution, one requires the wall stress as a boundary condition. The goal of this work is to determine the extent to which equilibrium and boundary layer assumptions are valid in the near-wall regions, to develop models for the inner layer based on such assumptions, and to test these modeling ideas in some relatively simple flows with different pressure gradients, such as channel flow and flow over a backward-facing step. Ultimately, models that perform adequately in these situations will be applied to more complex flow configurations, such as an airfoil.

Description: Since its first application, the dynamic procedure has been recognized as an effective means to compute rather than prescribe the unknown coefficients that appear in a subgrid-scale model for Large-Eddy Simulation (LES). The dynamic procedure is usually used to determine the nondimensional coefficient in the Smagorinsky (1963) model. In reality the procedure is quite general and it is not limited to the Smagorinsky model by any theoretical or practical constraints. The purpose of this note is to consider a generalized family of dynamic eddy viscosity models that do not necessarily rely on the local equilibrium assumption built into the Smagorinsky model. By invoking an inertial range assumption, it will be shown that the coefficients in the new models need not be nondimensional. This additional degree of freedom allows the use of models that are scaled on traditionally unknown quantities such as the dissipation rate. In certain cases, the dynamic models with dimensional coefficients are simpler to implement, and allow for a 30% reduction in the number of required filtering operations.

Description: Dynamic subgrid models have proved to be remarkably successful in predicting the behavior of turbulent flows. Part of the reasons for their success are well understood. Since they are constructed to generate an effective viscosity which is proportional to some measure of the turbulent energy at the high wavenumber end of the spectrum, their eddy viscosity vanishes as the flow becomes laminar. This alone would justify their use over simpler models. But beyond this obvious advantage, which is confined to inhomogeneous and evolving flows, the reason why they also work better in simpler homogeneous cases, and how they do it without any obvious adjustable parameter, is not clear. This lack of understanding of the internal mechanisms of a useful tool is disturbing, not only as an intellectual challenge, but because it raises the doubt of whether it will work in all cases. This note is an attempt to clarify those mechanisms. We will see why dynamic models are robust and how they can get away with even comparatively gross errors in their formulations. This will suggest that they are only particular cases of a larger family of robust models, all of which would be relatively insensitive to large simplifications in the physics of the flow. We will also construct some such models, although mostly as research tools. It will turn out, however, that the standard dynamic formulation is not only robust to errors, but also behaves as if it were substantially well formulated. The details of why this is so will still not be clear at the end of this note, specially since it will be shown that the 'a priori' testing of the stresses gives, as is usual in most subgrid models, very poor results. But it will be argued that the basic reason is that the dynamic formulation mimics the condition that the total dissipation is approximately equal to the production measured at the test filter level.

Description: All numerical simulations of turbulence (DNS or LES) involve some discretization errors. The integrity of such simulations therefore depend on our ability to quantify and control such errors. In the classical literature on analysis of errors in partial differential equations, one typically studies simple linear equations (such as the wave equation or Laplace's equation). The qualitative insight gained from studying such simple situations is then used to design numerical methods for more complex problems such as the Navier-Stokes equations. Though such an approach may seem reasonable as a first approximation, it should be recognized that strongly nonlinear problems, such as turbulence, have a feature that is absent in linear problems. This feature is the simultaneous presence of a continuum of space and time scales. Thus, in an analysis of errors in the one dimensional wave equation, one may, without loss of generality, rescale the equations so that the dependent variable is always of order unity. This is not possible in the turbulence problem since the amplitudes of the Fourier modes of the velocity field have a continuous distribution. The objective of the present research is to provide some quantitative measures of numerical errors in such situations. Though the focus of this work is LES, the methods introduced here can be just as easily applied to DNS. Errors due to discretization of the time-variable are neglected for the purpose of this analysis.

Description: The structure of pressure-pressure correlations at the interface of an incompressible steady-state turbulent flow with a rigid boundary was investigated. For the sake of completeness, the absolute value of the correlation between two random varying functions is herein defined as a number greater than or equal to zero and less than or equal to unity which is a measure of that fraction of one of the functions that 'follows' the second function (or vice versa). It was found that the soughtafter correlations can be determined by consideration of the high Re Navier-Stokes equation, but that the complexity of boundary layer turbulence, in particular the inhomogeneity perpendicular to the boundary and the anisotropy due to convective flow gradients, makes the structure of said correlations extremely difficult to assess. One of the earlier researchers in this field described the quantity under present consideration as 'a quantity which is beyond assessment.' Nonetheless, it was found that under some rather simplifying assumptions the determination of the required structure necessitates the formulation of the related structure of second order two-point correlations of turbulent velocity gradients, as well as third order two-point correlations of velocity gradients. The presence of these latter gradients is due to the nonlinearity in the turbulence ('turbulence self-interaction'). Both of these correlations are scaled, although not similarly, by factors dependent upon the magnitude of the convective flow, which can be modeled using a log law approximation. Fourth order correlations, although present, can be ignored, since they constitute 'higher order terms.' In a slightly more complex situation, it was found that convective flow gradients also have to be incorporated. At the moment, no definitive algebraic information peculiar to pressure-pressure correlations is available in the most highly idealized cases.

Description: Grid generation plays an integral part in the solution of computational fluid dynamics problems for aerodynamics applications. A major difficulty with standard structured grid generation, which produces quadrilateral (or hexahedral) elements with implicit connectivity, has been the requirement for a great deal of human intervention in developing grids around complex configurations. This has led to investigations into unstructured grids with explicit connectivities, which are primarily composed of triangular (or tetrahedral) elements, although other subdivisions of convex cells may be used. The existence of large gradients in the solution of aerodynamic problems may be exploited to reduce the computational effort by using high aspect ratio elements in high gradient regions. However, the heuristic approaches currently in use do not adequately address this need for high aspect ratio unstructured grids. High aspect ratio triangulations very often produce the large angles that are to be avoided. Point generation techniques based on contour or front generation are judged to be the most promising in terms of being able to handle complicated multiple body objects, with this technique lending itself well to adaptivity. The eventual goal encompasses several phases: first, a partitioning phase, in which the Voronoi diagram of a set of points and line segments (the input set) will be generated to partition the input domain; second, a contour generation phase in which body-conforming contours are used to subdivide the partition further as well as introduce the foundation for aspect ratio control, and; third, a Steiner triangulation phase in which points are added to the partition to enable triangulation while controlling angle bounds and aspect ratio. This provides a combination of the advancing front/contour techniques and refinement. By using a front, aspect ratio can be better controlled. By using refinement, bounds on angles can be maintained, while attempting to minimize the number of Steiner points.

Description: Asymptotic methods are used to describe the nonlinear self-interaction between pairs of oblique instability modes that eventually develops when initially linear spatially growing instability waves evolve downstream in nominally two-dimensional laminar boundary layers. The first nonlinear reaction takes place locally within a so-called 'critical layer', with the flow outside this layer consisting of a locally parallel mean flow plus a pair of oblique instability waves - which may or may not be accompanied by an associated plane wave. The amplitudes of these waves, which are completely determined by nonlinear effects within the critical layer, satisfy either a single integro-differential equation or a pair of integro-differential equations with quadratic to quartic-type nonlinearities. The physical implications of these equations are discussed.

Description: A new optical instrument, the liquid crystal point diffraction interferometer (LCPDI), is used to measure the temperature distribution across a heated chamber filled with silicone oil. Data taken using the LCPDI are compared to equivalent measurements made with a traversing thermocouple and the two data sets show excellent agreement This instrument maintains the compact, robust design of Linnik's point diffraction interferometer and adds to it phase stepping capability for quantitative interferogram analysis. The result is a compact, simple to align, environmentally insensitive interferometer capable of accurately measuring optical wavefronts with very high data density and with automated data reduction.

Description: Using a comprehensive flight test database and a parameter identification software program produced at NASA Ames Research Center, a math model of the longitudinal aerodynamics of the Harrier aircraft was formulated. The identification program employed the equation error method using multiple linear regression to estimate the nonlinear parameters. The formulated math model structure adhered closely to aerodynamic and stability/control theory, particularly with regard to compressibility and dynamic manoeuvring. Validation was accomplished by using a three degree-of-freedom nonlinear flight simulator with pilot inputs from flight test data. The simulation models agreed quite well with the measured states. It is important to note that the flight test data used for the validation of the model was not used in the model identification.

Description: A direct numerical simulation (DNS) algorithm has been developed and validated for use in the investigation of crossflow instability on supersonic swept wings, an application of potential relevance to the design of the High-Speed Civil Transport (HSCT). The algorithm is applied to the investigation of stationary crossflow instability on an infinitely long 77-degree swept wing in Mach 3.5 flow. The results of the DNS are compared with the predictions of linear parabolized stability equation (PSE) methodology. In-general, the DNS and PSE results agree closely in terms of modal growth rate, structure, and orientation angle. Although further validation is needed for large-amplitude (nonlinear) disturbances, the close agreement between independently derived methods offers preliminary validation of both DNS and PSE approaches.

Description: Two F-18 aircraft were flown, one above the other, in two formations, in order for the shock systems of the two aircraft to merge and propagate to the ground. The first formation had the canopy of the lower F-18 in the inlet shock of the upper F-18 (called inlet-canopy). The flight conditions were Mach 1.22 and an altitude of 23,500 ft. An array of five sonic boom recorders was used on the ground to record the sonic boom signatures. This paper describes the flight test technique and the ground level sonic boom signatures. The tail-canopy formation resulted in two, separated, N-wave signatures. Such signatures probably resulted from aircraft positioning error. The inlet-canopy formation yielded a single modified signature; two recorders measured an approximate flattop signature. Loudness calculations indicated that the single inlet-canopy signatures were quieter than the two, separated tail-canopy signatures. Significant loudness occurs after a sonic boom signature. Such loudness probably comes from the aircraft engines.

Description: A multiblock, discrete sensitivity analysis method is used to couple a direct optimization method and a flow analysis method. The domain is divided into smaller subdomains for which the sensitivities are obtained separately. Then, an effective sensitivity equation is solved to complete the coupling of all the sensitivity information. The flow analysis is based on the thin-layer Navier-Stokes equations solved by an implicit, upwind-biased, finite-volume method. The method of feasible directions is used for the present gradient-based optimization approach. First, a transonic airfoil is optimized to investigate the behavior of the method in highly nonlinear flows as well as the effect of different blocking strategies on the procedure. A supercritical airfoil is produced from an initially symmetric airfoil with multiblocking affecting the path but not the final shape. Secondly, a two-element airfoil is shape optimized in subsonic flow to demonstrate the present method's capability of shaping aerodynamically interfering elements simultaneously. For a very low and a very high Reynolds number cases, the shape of the main airfoil and the flap are optimized to yield improved lift-to-drag ratios.

Description: This report summarizes some NASA Lewis (i.e., government owned) computer codes capable of being used for airbreathing propulsion system studies to determine the design geometry and to predict the design/off-design performance of compressors and turbines. These are not CFD codes; velocity-diagram energy and continuity computations are performed fore and aft of the blade rows using meanline, spanline, or streamline analyses. Losses are provided by empirical methods. Both axial-flow and radial-flow configurations are included.

Description: A far-wing theory in which the validity of the detailed balance principle is maintained in each step of the derivation is presented. The role of the total density matrix including the initial correlations is analyzed rigorously. By factoring out the rapidly varying terms in the complex-time development operator in the interaction representation, better approximate expressions can be obtained. As a result, the spectral density can be expressed in terms of the line-coupling functions in which two coupled lines are arranged symmetrically and whose frequency detunings are omega - 1/2(omega(sub ji) + omega (sub j'i'). Using the approximate values omega - omega(sub ji) results in expressions that do not satisfy the detailed balance principle. However, this principle remains satisfied for the symmetrized spectral density in which not only the coupled lines are arranged symmetrically, but also the initial and final states belonging to the same lines are arranged symmetrically as well.

Description: To identify planform characteristics which have promise for a highly maneuverable vehicle, an investigation was conducted in the Langley Subsonic Basic Research Tunnel to determine the low-speed longitudinal aerodynamics of 21 planform geometries. Concepts studied included twin bodies, double wings, cutout wings, and serrated forebodies. The planform models tested were all 1/4-in.-thick flat plates with beveled edges on the lower surface to ensure uniform flow separation at angle of attack. A 1.0-in.-diameter cylindrical metric body with a hemispherical nose was used to house the six-component strain gauge balance for each configuration. Aerodynamic force and moment data were obtained across an angle-of-attack range of 0 to 70 deg with zero sideslip at a free-stream dynamic pressure of 30 psf. Surface flow visualization studies were also conducted on selected configurations using fluorescent minitufts. Results from the investigation indicate that a cutout wing planform can improve lift characteristics; however, cutout size, shape, and position and wing leading-edge sweep will all influence the effectiveness of the cutout configuration. Tests of serrated forebodies identified this concept as an extremely effective means of improving configuration lift characteristics; increases of up to 25 percent in the value of maximum lift coefficient were obtained.

Description: A general multi-block three-dimensional volume grid generator is presented which is suitable for Multi-Disciplinary Design Optimization. The code is fast, robust, highly automated, and written in ANSI C for platform independence. Algebraic techniques are used to generate and/or modify block face and volume grids to reflect geometric changes resulting from design optimization. Volume grids are generated/modified in a batch environment and controlled via an ASCII user input deck. This allows the code to be incorporated directly into the design loop. Generated volume grids are presented for a High Speed Civil Transport (HSCT) Wing/Body geometry as well a complex HSCT configuration including horizontal and vertical tails, engine nacelles and pylons, and canard surfaces.

Description: An evaluation of the effect of model inlet air temperature drift during a test run was performed to aid in the decision on the need for and/or the schedule for including heaters in the SRMAFTE. The Sverdrup acceptance test data was used to determine the drift in air temperature during runs over the entire range of delivered flow rates and pressures. The effect of this temperature drift on the model Reynolds number was also calculated. It was concluded from this study that a 2% change in absolute temperature during a test run could be adequately accounted for by the data analysis program. A handout package of these results was prepared and presented to ED35 management.

Description: The purpose of the RSRM Nozzle Slag Ejection Precursor Test is to investigate the effect that slag ejection from the RSRM nozzle has on the chamber pressure and trust of the SRB's. In past firings of the Reusable Solid Rocket Motor (RSRM) both static test and flight motors have shown small pressure perturbations occurring primarily between 65 and 80 seconds. A joint NASA/Thiokol team investigation concluded that the cause of the pressure perturbations was the periodic ingestion and ejection of molten aluminum oxide slag from the cavity around the submerged nozzle nose which tends to trap and collect individual aluminum oxide droplets from the approach flow. The conclusions of the team were supported by numerous data and observations from special tests including high speed photographic films, real time radiography, plume calorimeters, accelerometers, strain gauges, nozzle TVC system force gauges, and motor pressure and thrust data. A simplistic slag ballistics model was formulated to relate a given pressure perturbation to a required slag quantity. Also, a cold flow model using air and water was developed to provide data on the relationship between the slag flow rate and the chamber pressure increase. Both the motor and the cold flow model exhibited low frequency oscillations in conjunction with periods of slag ejection. Motor and model frequencies were related to scaling parameters. The data indicate that there is a periodicity to the slag entrainment and ejection phenomena which is possibly related to organized oscillations from instabilities in the dividing streamline shear layer which impinges on the underneath surface of the nozzle.

Description: We describe a system for interactive visualization and tracking of a 3-D unsteady computational fluid dynamics (CFD) simulation on a parallel computer. CM/AVS, a distributed, parallel implementation of a visualization environment (AVS) runs on the CM-5 parallel supercomputer. A CFD solver is run as a CM/AVS module on the CM-5. Data communication between the solver, other parallel visualization modules, and a graphics workstation, which is running AVS, are handled by CM/AVS. Partitioning of the visualization task, between CM-5 and the workstation, can be done interactively in the visual programming environment provided by AVS. Flow solver parameters can also be altered by programmable interactive widgets. This system partially removes the requirement of storing large solution files at frequent time steps, a characteristic of the traditional 'simulate (yields) store (yields) visualize' post-processing approach.

Description: An investigation was conducted in the Langley 16-Foot Transonic Tunnel to determine the multiaxis thrust-vectoring characteristics of the F-18 High-Alpha Research Vehicle (HARV). A wingtip supported, partially metric, 0.10-scale jet-effects model of an F-18 prototype aircraft was modified with hardware to simulate the thrust-vectoring control system of the HARV. Testing was conducted at free-stream Mach numbers ranging from 0.30 to 0.70, at angles of attack from O' to 70', and at nozzle pressure ratios from 1.0 to approximately 5.0. Results indicate that the thrust-vectoring control system of the HARV can successfully generate multiaxis thrust-vectoring forces and moments. During vectoring, resultant thrust vector angles were always less than the corresponding geometric vane deflection angle and were accompanied by large thrust losses. Significant external flow effects that were dependent on Mach number and angle of attack were noted during vectoring operation. Comparisons of the aerodynamic and propulsive control capabilities of the HARV configuration indicate that substantial gains in controllability are provided by the multiaxis thrust-vectoring control system.

Description: A small scale ground effect test rig was used to study the ground plane flow field generated by a STOVL aircraft in hover. The objective of the research was to support NASA-Ames Research Center planning for the Large Scale Powered Model (LSPM) test for the ARPA-sponsored ASTOVL program. Specifically, small scale oil flow visualization studies were conducted to make a relative assessment of the aerodynamic interference of a proposed strut configuration and a wall configuration on the ground plane stagnation line. A simplified flat plate model representative of a generic jet-powered STOVL aircraft was used to simulate the LSPM. Cold air jets were used to simulate both the lift fan and the twin rear engines. Nozzle Pressure Ratios were used that closely represented those used on the LSPM tests. The flow visualization data clearly identified a shift in the stagnation line location for both the strut and the wall configuration. Considering the experimental uncertainty, it was concluded that either the strut configuration o r the wall configuration caused only a minor aerodynamic interference.

Description: An all-at-once reduced Hessian Successive Quadratic Programming (SQP) scheme has been shown to be efficient for solving aerodynamic design optimization problems with a moderate number of design variables. This paper extends this scheme to allow solution refining. In particular, we introduce a reduced Hessian refining technique that is critical for making a smooth transition of the Hessian information from coarse grids to fine grids. Test results on a nozzle design using quasi-one-dimensional Euler equations show that through solution refining the efficiency and the robustness of the all-at-once reduced Hessian SQP scheme are significantly improved.

Description: This paper introduces a computational scheme for solving a class of aerodynamic design problems that can be posed as nonlinear equality constrained optimizations. The scheme treats the flow and design variables as independent variables, and solves the constrained optimization problem via reduced Hessian successive quadratic programming. It updates the design and flow variables simultaneously at each iteration and allows flow variables to be infeasible before convergence. The solution of an adjoint flow equation is never needed. In addition, a range space basis is chosen so that in a certain sense the 'cross term' ignored in reduced Hessian SQP methods is minimized. Numerical results for a nozzle design using the quasi-one-dimensional Euler equations show that this scheme is computationally efficient and robust. The computational cost of a typical nozzle design is only a fraction more than that of the corresponding analysis flow calculation. Superlinear convergence is also observed, which agrees with the theoretical properties of this scheme. All optimal solutions are obtained by starting far away from the final solution.

Description: A simple parameterization scheme for a complex turbulent flow using nondimensional parameters coming from the Reynolds stress equations is given. Definitions and brief descriptions of the physical significance of several nondimensional parameters that are used to characterize turbulence from the viewpoint of single-point turbulence closures are given. These nondimensional parameters reflect measures of (1) the spectral band width of the turbulence; (2) deviations from the ideal Kolmogorov behavior; (3) the relative magnitude, orientation, and temporal duration of the deformation to which the turbulence is subjected; (4) one and two-point measures of the large and small scale anisotropy of the turbulence; and (5) inhomogeneity. This is an attempt to create a more systematic methodology for the diagnosis and classification of turbulent flows as well as in the development, validation, and application of turbulence model strategies. The parameters serve also to indicate the adequacy of various assumptions made in single-point turbulence models and in suggesting the appropriate turbulence strategy for a particular complex flow. The compilation will be of interest to experimentalists and to those involved in either computing turbulent flows or whose interests lies in verifying the adequacy of the phenomenological beliefs used in turbulence closures.

Description: A fully-explicit, self-consistent algebraic expression for the Reynolds stress, which is the exact solution to the Reynolds stress transport equation in the 'weak equilibrium' limit for two-dimensional mean flows for all linear and some quasi-linear pressure-strain models, is derived. Current explicit algebraic Reynolds stress models derived by employing the 'weak equilibrium' assumption treat the production-to-dissipation (P/epsilon) ratio implicitly, resulting in an effective viscosity that can be singular away from the equilibrium limit. In the present paper, the set of simultaneous algebraic Reynolds stress equations are solved in the full non-linear form and the eddy viscosity is found to be non-singular. Preliminary tests indicate that the model performs adequately, even for three dimensional mean flow cases. Due to the explicit and non-singular nature of the effective viscosity, this model should mitigate many of the difficulties encountered in computing complex turbulent flows with the algebraic Reynolds stress models.

Description: A design method has been developed by which an airfoil with a substantial amount of natural laminar flow can be designed, while maintaining other aerodynamic and geometric constraints. After obtaining the initial airfoil's pressure distribution at the design lift coefficient using an Euler solver coupled with an integml turbulent boundary layer method, the calculations from a laminar boundary layer solver are used by a stability analysis code to obtain estimates of the transition location (using N-Factors) for the starting airfoil. A new design method then calculates a target pressure distribution that will increase the larninar flow toward the desired amounl An airfoil design method is then iteratively used to design an airfoil that possesses that target pressure distribution. The new airfoil's boundary layer stability characteristics are determined, and this iterative process continues until an airfoil is designed that meets the laminar flow requirement and as many of the other constraints as possible.

Description: The invention is a method and apparatus for monitoring the presence, concentration, and the movement of fluids. It is based on utilizing electromagnetic measurements of the complex permittivity of the fluids for detecting and monitoring the fluid. More particularly the apparatus uses one or more microwave probes which are placed at the locations where the measurements are to be made. A radio frequency signal is transmitted to the probe and the reflected signal is phase and amplitude detected at a rapid rate for the purpose of identifying the fluids, based on their dielectric constant at the probe. The apparatus can be used for multiple purposes including measures of flow rates, turbulence, dispersion, fluid identification, and changes in flow conditions of multiple fluids or multiple states of a single fluid in a flowline or a holding container. The apparatus includes a probe consisting of two electrical conductors separated by an insulator. A radio frequency signal is communicated to the probe and is reflected back from the portion of the probe exposed to the fluid. The radio frequency signal also provides a reference signal. An oscillator generates a second signal which combined with each of the reference signal and the reflected signal to produce signals of lower frequencies to facilitate filtering and amplifying those signals. The two signals are then mixed in a detector to produce an output signal that is representative of the phase and amplitude change caused by the reflection of the signal at the probe exposed to the fluid. The detector may be a dual phase detector that provides two such output signals that are in phase quadrature. A phase shifter may be provided for selectively changing the phase of the reference signal to improve the sensitivity of at least one of the output signals for more accurate readings and/or for calibration purposes. The two outputs that are in quadrature with respect to each other may be simultaneously monitored to account for drift errors. The output signals are digitized and provided to a computer at a sample rate which may be very high. The computer is operable to identify the fluid based on its complex permittivity as may be useful for identifying the flow rates, determining the fluid mixture ratio, detecting impurities in the fluid, and so forth. Novelty is believed to reside in the use of the real part of complex permittivity to measure small difference in permittivity of the fluid.

Description: In this paper, a fast Poisson solver for unsteady, incompressible Navier-Stokes equations with finite difference methods on the non-uniform, half-staggered grid is presented. To achieve this, new algorithms for diagonalizing a semi-definite pair are developed. Our fast solver can also be extended to the three dimensional case. The motivation and related issues in using this second kind of staggered grid are also discussed. Numerical testing has indicated the effectiveness of this algorithm.

Description: The absorption of vacuum ultraviolet light by atomic oxygen has been measured in the Electric Arc-driven Shock Tube (EAST) Facility at NASA-Ames Research Center. This investigation demonstrates the instrumentation required to determine atomic oxygen concentrations from absorption measurements in impulse facilities. A shock wave dissociates molecular oxygen, producing a high temperature sample of atomic oxygen in the shock tube. A probe beam is generated with a Raman-shifted ArF excimer laser. By suitable tuning of the laser, absorption is measured over a range of wavelengths in the region of the atomic line at 130.49 nm. The line shape function is determined from measurements at atomic oxygen densities of 3 x 10(exp 17) and 9 x 10(exp 17) cm(exp -3). The broadening coefficient for resonance interactions is deduced from this data, and this value is in accord with available theoretical models.

Description: Design for prevention of aeroelastic instability (that is, the critical speeds leading to aeroelastic instability lie outside the operating range) is an integral part of the wing design process. Availability of the sensitivity derivatives of the various critical speeds with respect to shape parameters of the wing could be very useful to a designer in the initial design phase, when several design changes are made and the shape of the final configuration is not yet frozen. These derivatives are also indispensable for a gradient-based optimization with aeroelastic constraints. In this study, flutter characteristic of a typical section in subsonic compressible flow is examined using a state-space unsteady aerodynamic representation. The sensitivity of the flutter speed of the typical section with respect to its mass and stiffness parameters, namely, mass ratio, static unbalance, radius of gyration, bending frequency, and torsional frequency is calculated analytically. A strip theory formulation is newly developed to represent the unsteady aerodynamic forces on a wing. This is coupled with an equivalent plate structural model and solved as an eigenvalue problem to determine the critical speed of the wing. Flutter analysis of the wing is also carried out using a lifting-surface subsonic kernel function aerodynamic theory (FAST) and an equivalent plate structural model. Finite element modeling of the wing is done using NASTRAN so that wing structures made of spars and ribs and top and bottom wing skins could be analyzed. The free vibration modes of the wing obtained from NASTRAN are input into FAST to compute the flutter speed. An equivalent plate model which incorporates first-order shear deformation theory is then examined so it can be used to model thick wings, where shear deformations are important. The sensitivity of natural frequencies to changes in shape parameters is obtained using ADIFOR. A simple optimization effort is made towards obtaining a minimum weight design of the wing, subject to flutter constraints, lift requirement constraints for level flight and side constraints on the planform parameters of the wing using the IMSL subroutine NCONG, which uses successive quadratic programming.

Description: The work in this report was conducted at NASA Ames Research Center during the period from August 1993 to January 1995 deals with the direct numerical simulation of transitional and turbulent flow at low Mach numbers using high-order-accurate finite-difference techniques. A computation of transition to turbulence of the spatially-evolving boundary layer on a heated flat plate in the presence of relatively high freestream turbulence was performed. The geometry and flow conditions were chosen to match earlier experiments. The development of the momentum and thermal boundary layers was documented. Velocity and temperature profiles, as well as distributions of skin friction, surface heat transfer rate, Reynolds shear stress, and turbulent heat flux were shown to compare well with experiment. The numerical method used here can be applied to complex geometries in a straightforward manner.

Description: An investigation has been conducted in the Langley 7- by 10-Foot High Speed Wind Tunnel to determine the longitudinal and lateral directional aerodynamic characteristics of a series of personnel launch system concepts. This series of configurations evolved during an effort to improve the subsonic characteristics of a proposed lifting entry vehicle (designated the HL-20). The primary purpose of the overall investigation was to provide a vehicle concept which was inherently stable and trimable from entry to landing while examining methods of improving subsonic aerodynamic performance.

Description: The process by which a laminar boundary layer internalizes the external disturbances in the form of instability waves is known as boundary-layer receptivity. The objective of the present research was to determine the effect of acoustic excitation on boundary-layer receptivity for a flat plate with distributed variable-amplitude surface roughness through measurements with a hot-wire probe. Tollmien-Schlichting mode shapes due to surface roughness receptivity have also been determined, analyzed, and shown to be in agreement with theory and other experimental work. It has been shown that there is a linear relationship between the surface roughness and receptivity for certain roughness configurations with constant roughness wavelength. In addition, strong non-linear receptivity effects exist for certain surface roughness configurations over a band where the surface roughness and T-S wavelength are matched. The results from the present experiment follow the trends predicted by theory and other experimental work for linear receptivity. In addition, the results show the existence of non-linear receptivity effects for certain combinations of surface roughness elements.

Description: A three-dimensional multiblock Navier-Stokes code, PAB3D, which was developed for propulsion integration and general aerodynamic analysis, has been used extensively by NASA Langley and other organizations to perform both internal (exhaust) and external flow analysis of complex aircraft configurations. This code was designed to solve the simplified Reynolds Averaged Navier-Stokes equations. A two-equation k-epsilon turbulence model has been used with considerable success, especially for attached flows. Accurate predicting of transonic shock wave location and pressure recovery in separated flow regions has been more difficult. Two algebraic Reynolds stress models (ASM) have been recently implemented in the code that greatly improved the code's ability to predict these difficult flow conditions. Good agreement with Direct Numerical Simulation (DNS) for a subsonic flat plate was achieved with ASM's developed by Shih, Zhu, and Lumley and Gatski and Speziale. Good predictions were also achieved at subsonic and transonic Mach numbers for shock location and trailing edge boattail pressure recovery on a single-engine afterbody/nozzle model.

Description: One source contributing to wind tunnel background noise is microphone self-noise. An experiment was conducted to investigate the flow-induced acoustic oscillations of Bruel & Kjaer (B&K) in-flow microphones. The results strongly suggest the B&K microphone cavity behaves more like an open cavity. Their cavity acoustic oscillations are likely caused by strong interactions between the cavity shear layer and the cavity trailing edge. But the results also suggest that cavity shear layer oscillations could be coupled with cavity acoustic resonance to generate tones. Detailed flow velocity measurements over the cavity screen have shown inflection points in the mean velocity profiles and high disturbance and spectral intensities in the vicinity of the cavity trailing edge. These results are the evidence for strong interactions between cavity shear layer oscillations and the cavity trailing edge. They also suggest that beside acoustic signals, the microphone inside the cavity has likely recorded hydrodynamic pressure oscillations, too. The results also suggest that the forebody shape does not have a direct effect on cavity oscillations. For the FITE (Flow Induced Tone Eliminator) microphone, it is probably the forebody length and the resulting boundary layer turbulence that have made it work. Turbulence might have thickened the boundary layer at the separation point, weakened the shear layer vortices, or lifted them to miss impinging on the cavity trailing edge. In addition, the study shows that the cavity screen can modulate the oscillation frequency but not the cavity acoustic oscillation mechanisms.

Description: Three-dimensional supersonic viscous laminar flows over symmetric corners are considered in this paper. The characteristic features of such configurations are discussed and an historical survey on the past research work is presented. A new contribution based on a numerical technique that solves the parabolized form of the Navier-Stokes equations is presented. Such a method makes it possible to obtain very detailed descriptions of the flowfield with relatively modest CPU time and memory storage requirements. The numerical approach is based on a space-marching technique, uses a finite volume discretization and an upwind flux-difference splitting scheme (developed for the steady flow equations) for the evaluation of the inviscid fluxes. Second order accuracy is reached following the guidelines of the ENO schemes. Different free-stream conditions and geometrical configurations are considered. Primary and secondary streamwise vortical structures embedded in the boundary layer and originated by the interaction of the latter with shock waves are detected and studied. Computed results are compared with experimental data taken from literature.

Description: The computational fluid dynamics code, PARC3D, is tested to see if its use of non-physical artificial dissipation affects the accuracy of its results. This is accomplished by simulating a shock-laminar boundary layer interaction and several hypersonic flight conditions of the Pegasus(TM) launch vehicle using full artificial dissipation, low artificial dissipation, and the Engquist filter. Before the filter is applied to the PARC3D code, it is validated in one-dimensional and two-dimensional form in a MacCormack scheme against the Riemann and convergent duct problem. For this explicit scheme, the filter shows great improvements in accuracy and computational time as opposed to the nonfiltered solutions. However, for the implicit PARC3D code it is found that the best estimate of the Pegasus experimental heat fluxes and surface pressures is the simulation utilizing low artificial dissipation and no filter. The filter does improve accuracy over the artificially dissipative case but at a computational expense greater than that achieved by the low artificial dissipation case which has no computational time penalty and shows better results. For the shock-boundary layer simulation, the filter does well in terms of accuracy for a strong impingement shock but not as well for weaker shock strengths. Furthermore, for the latter problem the filter reduces the required computational time to convergence by 18.7 percent.

Description: An optimization procedure is developed for the simultaneous improvement of the aerodynamic and sonic boom characteristics of high speed aircraft. From a sonic boom perspective, it is desirable to minimize the first peak in the overpressure signal at a specified distance away from the aircraft. From aerodynamic point of view, the aerodynamic drag coefficient ratio must be minimized while maintaining the lift coefficient at desired level. The optimization procedure is applied to wing-body configurations related to high speed aircraft. The objectives of this current research are: (1) development of a multiobjective optimization procedure for aerospace vehicles with the integration of sonic boom and aerodynamic performance criteria; and (2) development of semi-analytical approach for calculating sonic boom design sensitivities.

Description: During the Higher Harmonic Control Aeroacoustic Rotor Test, extensive measurements of the rotor aerodynamics, the far-field acoustics, the wake geometry, and the blade motion for powered, descent, flight conditions were made. These measurements have been used to validate and improve the prediction of blade-vortex interaction (BVI) noise. The improvements made to the BVI modeling after the evaluation of the test data are discussed. The effects of these improvements on the acoustic-pressure predictions are shown. These improvements include restructuring the wake, modifying the core size, incorporating the measured blade motion into the calculations, and attempting to improve the dynamic blade response. A comparison of four different implementations of the Ffowcs Williams and Hawkings equation is presented. A common set of aerodynamic input has been used for this comparison.

Description: Measurements from transitional, heated boundary layers along a concave-curved test wall are presented and discussed. A boundary layer subject to low free-stream turbulence intensity (FSTI), which contains stationary streamwise (Gortler) vortices, is documented. The low FSTI measurements are followed by measurements in boundary layers subject to high (initially 8%) free-stream turbulence intensity and moderate to strong streamwise acceleration. Conditions were chosen to simulate those present on the downstream half of the pressure side of a gas turbine airfoil. Mean flow characteristics as well as turbulence statistics, including the turbulent shear stress, turbulent heat flux, and turbulent Prandtl number, are documented. A technique called "octant analysis" is introduced and applied to several cases from the literature as well as to data from the present study. Spectral analysis was applied to describe the effects of turbulence scales of different sizes during transition. To the authors'knowledge, this is the first detailed documentation of boundary layer transition under such high free-stream turbulence conditions.

Description: Preliminary measurements have been made of the flow over the tip of an unswept wing flap. To achieve an acceptable Reynolds number based on flap chord, the flap chord was chosen equal to the chord of the main airfoil (c = 19 in. approx. 0.48 m). The model was mounted in a 30 in. x 30 in. wind tunnel running at up to 100 ft/sec. (30 m/s): severe wind-tunnel interference was accepted, and any computations would be done using the tunnel walls as the boundaries of the computational domain. Maximum Reynolds number based on flap chord and tunnel speed was about 1.O x lO(exp 6). The grant ended before a full set of measurements could be made, but the work done so far yields a useful picture of the flow. The vortex originates at about mid-chord on the flap and rises rapidly above the chord line. It has a concentrated core, with total pressure lower than the ambient static pressure, and there is no evidence of large-scale wandering. A simple method of model construction, giving light weight and excellent surface finish, was developed.

Description: A method using a Green's function is developed for computing transient temperatures in a semitransparent layer by using the two-flux method coupled with the transient energy equation. Each boundary of the layer is exposed to a hot or cold radiative environment, and is heated or cooled by convection. The layer refractive index is larger than one, and the effect of internal reflections is included with the boundaries assumed diffuse. The analysis accounts for internal emission, absorption, heat conduction, and isotropic scattering. Spectrally dependent radiative properties are included, and transient results are given to illustrate two-band spectral behavior with optically thin and thick bands. Transient results using the present Green's function method are verified for a gray layer by comparison with a finite difference solution of the exact radiative transfer equations; excellent agreement is obtained. The present method requires only moderate computing times and incorporates isotropic scattering without additional complexity. Typical temperature distributions are given to illustrate application of the method by examining the effect of strong radiative heating on one side of a layer with convective cooling on the other side, and the interaction of strong convective heating with radiative cooling from the layer interior.

Description: An extensive quantity of airload measurements was obtained for a pressure-instrumented model of the BO-105 main rotor for a large number of higher-harmonic control (HHC) settings at Duits-Nederlandse Wind Tunnel (DNW). The wake geometry, vortex strength, and vortex core size were also measured through a laser light sheet technique and LDV. These results are used to verify the BVI airload prediction methodologies developed by AFDD, DLR, NASA Langley, and ONERA. The comparisons show that an accurate prediction of the blade motion and the wake geometry is the most important aspect of the BVI airload predictions.

Description: This is a guide for the use of the pressure disk rotor model that has been placed in the incompressible Navier-Stokes code INS3D-UP. The pressure disk rotor model approximates a helicopter rotor or propeller in a time averaged manner and is intended to simulate the effect of a rotor in forward flight on the fuselage or the effect of a propeller on other aerodynamic components. The model uses a modified actuator disk that allows the pressure jump across the disk to vary with radius and azimuth. The cyclic and collective blade pitch angles needed to achieve a specified thrust coefficient and zero moment about the hub are predicted. The method has been validated with experimentally measured mean induced inflow velocities as well as surface pressures on a generic fuselage. Overset grids, sometimes referred to as Chimera grids, are used to simplify the grid generation process. The pressure disk model is applied to a cylindrical grid which is embedded in the grid or grids used for the rest of the configuration. This document will outline the development of the method, and present input and results for a sample case.

Description: Future hypersonic vehicles are going to be designed largely with computational fluid dynamic methods based on appropriate physical models. The question on how much of this design process can be completed with the present state of computational aerothermodynamics is addressed. Some limitations of current models are discussed. It is shown that much more research is required before it will be possible to accurately design a hypersonic vehicle for all of its flight conditions. The quantities that must be computed accurately so that a minimum weight hypersonic vehicle can be designed are discussed. The use of computational fluid dynamics methods coupled with current thermochemical models in order to compute the quantities under specific flow conditions is considered.

Description: An experimental investigation was conducted to determine the aerodynamic characteristics of a store as it was separated from the lee side of a flat plate inclined at 15 deg to the free-stream flow at Mach 6. Two store models were tested: a cone cylinder and a roof delta. Force and moment data were obtained for both stores as they were moved in 0.5-in. increments away from the flat plate lee-side separated flow region into the free-stream flow while the store angle of attack was held constant at either 0 deg or 15 deg. The results indicate that both stores had adverse separation characteristics (i.e., negative normal force and pitching moment) at an angle of attack of 0 deg, and the cone cylinder had favorable separation characteristics (i.e., positive normal force and pitching moment) at an angle of attack of 15 deg. At an angle of attack of 15 deg, the separation characteristics of the roof delta are indeterminate at small separation distances and favorable at greater separation distances. These characteristics are the result of the local flow inclination relative to the stores as they traversed through the flat plate lee-side flow field. In addition to plotted data, force and moment data are tabulated and schlieren photographs of the stores and flat plate are presented.

Description: A high Reynolds number investigation of a commercial transport model was conducted in the National Transonic Facility (NTF) at Langley Research Center. This investigation was part of a cooperative effort to test a 0.03-scale model of a Boeing 767 airplane in the NTF over a Mach number range of 0.70 to 0.86 and a Reynolds number range of 2.38 to 40.0 x 10(exp 6) based on the mean aerodynamic chord. One of several specific objectives of the current investigation was to evaluate the level of data repeatability attainable in the NTF. Data repeatability studies were performed at a Mach number of 0.80 with Reynolds numbers of 2.38, 4.45, and 40.0 x 10(exp 6) and also at a Mach number of 0.70 with a Reynolds number of 40.0 x 10(exp 6). Many test procedures and data corrections are addressed in this report, but the data presented do not include corrections for wall interference, model support interference, or model aeroelastic effects. Application of corrections for these three effects would not affect the results of this study because the corrections are systematic in nature and are more appropriately classified as sources of bias error. The repeatability of the longitudinal stability-axis force and moment data has been accessed. Coefficients of lift, drag, and pitching moment are shown to repeat well within the pretest goals of plus or minus 0.005, plus or minus 0.0001, and plus or minus 0.001, respectively, at a 95-percent confidence level over both short- and near-term periods.

Description: An experimental investigation was conducted to determine the effect of diverter wedge half-angle and nacelle lip height on the drag characteristics of an assembly consisting of a nacelle fore cowl from a typical high-speed civil transport (HSCT) and a diverter mounted on a flat plate. Data were obtained for diverter wedge half-angles of 4.0 deg, 6.0 deg, and 8.0 deg and ratios of the nacelle lip height above a flat plate to the boundary-layer thickness (h(sub n)/delta) of approximately 0.87 to 2.45. Limited drag data were also obtained on a complete nacelle/diverter configuration that included fore and aft cowls. Although the nacelle/diverter drag data were not corrected for base pressures or internal flow drag, the data are useful for comparing the relative drag of the configuration tested. The tests were conducted in the Langley Unitary Plan Wind Tunnel at Mach numbers of 1.50, 1.80, 2.10, and 2.40 and Reynolds numbers ranging from 2.00 x 10(exp 6) to 5.00 x 10(exp 6) per foot. The results of this investigation showed that the nacelle/diverter drag essentially increased linearly with increasing h(sub n)/delta except near 1.0 where the data showed a nonlinear behavior. This nonlinear behavior was probably caused by the interaction of the shock waves from the nacelle/diverter configuration with the flat-plate boundary layer. At the lowest h(sub n)/delta tested, the diverter wedge half-angle had virtually no effect on the nacelle/diverter drag. However, as h(sub n)/delta increased, the nacelle/diverter drag increased as diverter wedge half-angle increased.

Description: Water droplet trajectories within the NASA Lewis Research Center's Icing Research Tunnel (IRT) were studied through computer analysis. Of interest was the influence of the wind tunnel contraction and wind tunnel model blockage on the water droplet trajectories. The computer analysis was carried out with a program package consisting of a three-dimensional potential panel code and a three-dimensional droplet trajectory code. The wind tunnel contraction was found to influence the droplet size distribution and liquid water content distribution across the test section from that at the inlet. The wind tunnel walls were found to have negligible influence upon the impingement of water droplets upon a wing model.

Description: This document outlines the tests performed to make aerodynamic force and torque measurements on the SOFIA wind tunnel model telescope. These tests were performed during the SOFIA 2 wind tunnel test in the 14 ft wind tunnel during the months of June through August 1994. The test was designed to measure the dynamic cross elevation moment acting on the SOFIA model telescope due to aerodynamic loading. The measurements were taken with the telescope mounted in an open cavity in the tail section of the SOFIA model 747. The purpose of the test was to obtain an estimate of the full scale aerodynamic disturbance spectrum, by scaling up the wind tunnel results (taking into account differences in sail area, air density, cavity dimension, etc.). An estimate of the full scale cross elevation moment spectrum was needed to help determine the impact this disturbance would have on the telescope positioning system requirements. A model of the telescope structure, made of a light weight composite material, was mounted in the open cavity of the SOFIA wind tunnel model. This model was mounted via a force balance to the cavity bulkhead. Despite efforts to use a 'stiff' balance, and a lightweight model, the balance/telescope system had a very low resonant frequency (37 Hz) compared to the desired measurement bandwidth (1000 Hz). Due to this mechanical resonance of the balance/telescope system, the balance alone could not provide an accurate measure of applied aerodynamic force at the high frequencies desired. A method of measurement was developed that incorporated accelerometers in addition to the balance signal, to calculate the aerodynamic force.

Description: A computational fluid dynamics (CFD) analysis has been performed on the aft slot region of the Titan 4 Solid Rocket Motor Upgrade (SRMU). This analysis was performed in conjunction with MSFC structural modeling of the propellant grain to determine if the flow field induced stresses would adversely alter the propellant geometry to the extent of causing motor failure. The results of the coupled CFD/stress analysis have shown that there is a continual increase of flow field resistance at the aft slot due to the aft segment propellant grain being progressively moved radially toward the centerline of the motor port. This 'bootstrapping' effect between grain radial movement and internal flow resistance is conducive to causing a rapid motor failure.

Description: An all-at-once reduced Hessian Successive Quadratic Programming (SQP) scheme has been shown to be efficient for solving aerodynamic design optimization problems with a moderate number of design variables. This paper extends this scheme to allow solution refining. In particular, we introduce a reduced Hessian refining technique that is critical for making a smooth transition of the Hessian information from coarse grids to fine grids. Test results on a nozzle design using quasi-one-dimensional Euler equations show that through solution refining the efficiency and the robustness of the all-at-once reduced Hessian SQP scheme are significantly improved.

Description: Wind tunnel tests were made with a scale model of the HL-20 in the Langley Unitary Plan Wind Tunnel. Pitch control was investigated by deflecting the elevon surfaces on the outboard fins and body flaps on the fuselage. Yaw control tests were made with the all movable center fin deflected 5 deg. Almost full negative body flap deflection (-30 deg) was required to trim the HL-20 (moment reference center at 0.54-percent body length from nose) to positive values of life in the Mach number range from 1.6 to 2.5. Elevons were twice as effective as body flaps as a longitudinal trim device. The elevons were effective as a roll control, but because of tip-fin dihedral angle, produced about as much adverse yawing moment as rolling moment. The body flaps were less effective in producing rolling moment, but produced little adverse yawing moment. The yaw effectiveness of the all movable center fin was essentially constant over the angle-of-attack range at each Mach number. The value of yawing moment, however, was small. Center-fin deflection produced almost no rolling moments. The model was directionally unstable over most of the Mach number range with tip-fin dihedral angles less than the baseline value of 50 deg.

Description: Three-dimensional transient flow profiles of spin-up in a fully liquid filled cylinder from rest with gravity acceleration at various direction are numerically simulated and studied. Particular interests are concentrated on the development of temporary reverse flow zones and Ekman layer right after the impulsive start of spin-up from rest, and decay before the flow reaching to the solid rotation. Relationship of these flow developments and differences in the Reynolds numbers of the flow and its size selection of grid points concerning the numerical instabilities of flow computations are also discussed. In addition to the gravitational acceleration along the axial direction of the cylindrical container, a series of complicated flow profiles accompanied by three-dimensional transient flows with oblique gravitational acceleration has been studies.

Description: The work reported here pertains only to the first year of research for a three year proposal period. As a prelude to this two dimensional interface element, the one dimensional element was tested and errors were discovered in the code for built-up structures and curved interfaces. These errors were corrected and the benchmark Boeing composite crown panel was analyzed successfully. A study of various splines led to the conclusion that cubic B-splines best suit this interface element application. A least squares approach combined with cubic B-splines was constructed to make a smooth function from the noisy data obtained with random error in the coordinate data points of the Boeing crown panel analysis. Preliminary investigations for the formulation of discontinuous 2-D shell and 3-D solid elements were conducted.

Description: Test flights were conducted to evaluate the capability of Differential Global Positioning System (DGPS) to provide the accuracy and integrity required for International Civil Aviation Organization (ICAO) Category (CAT) III precision approach and landings. These test flights were part of a Federal Aviation Administration (FAA) program to evaluate the technical feasibility of using DGPS based technology for CAT III precision approach and landing applications. An IAI Westwind 1124 aircraft (N24RH) was equipped with DGPS receiving equipment and additional computing capability provided by E-Systems. The test flights were conducted at NASA Ames Research Center's Crows Landing Flight Facility, Crows Landing, California. The flight test evaluation was based on completing 100 approaches and landings. The navigation sensor error accuracy requirements were based on ICAO requirements for the Microwave Landing System (MLS). All of the approaches and landings were evaluated against ground truth reference data provided by a laser tracker. Analysis of these approaches and landings shows that the E-Systems DGPS system met the navigation sensor error requirements for a successful approach and landing 98 out of 100 approaches and landings, based on the requirements specified in the FAA CAT III Level 2 Flight Test Plan. In addition, the E-Systems DGPS system met the integrity requirements for a successful approach and landing or stationary trial for all 100 approaches and landings and all ten stationary trials, based on the requirements specified in the FAA CAT III Level 2 Flight Test Plan.

Description: This monograph provides an extensive list of formulas for airfoil polynomials. These polynomials provide convenient expansion functions for the description of the downwash and pressure distributions of linear theory for airfoils in both steady and unsteady subsonic flow.

Description: Low-speed wind-tunnel tests were conducted in the Langley 12-Foot Low-Speed Tunnel on a model of the Boeing Multirole Fighter (BMRF) aircraft. This single-seat, single-engine configuration was intended to be an F-16 replacement that would incorporate many of the design goals and advanced technologies of the F-22. Its mission requirements included supersonic cruise without afterburner, reduced observability, and the ability to attack both air-to-air and air-to-ground targets. So that it would be effective in all phases of air combat, the ability to maneuver at angles of attack up to and beyond maximum lift was also desired. Traditional aerodynamic yaw controls, such as rudders, are typically ineffective at these higher angles of attack because they are usually located in the wake from the wings and fuselage. For this reason, this study focused on investigating forebody-mounted controls that produces yawing moments by modifying the strong vortex flowfield being shed from the forebody at high angles of attack. Two forebody strakes were tested that varied in planform and chordwise location. Various patterns of porosity in the forebody skin were also tested that differed in their radial coverage and chordwise location. The tests were performed at a dynamic pressure of 4 lb/ft(exp 2) over an angle-of-attack range of -4 deg to 72 deg and a sideslip range of -10 deg to 10 deg. Static force data, static pressures on the surface of the forebody, and videotapes of flow-visualization using laser-illuminated smoke were obtained.