ALBERT

Description: The paper reviews a combined numerical and experimental activity on the Shuttle Orbiter, first performed at NASA Langley within the Orbiter Experiment (OEX) and subsequently at ESA, as part of the AGARD FDP WG 18 activities. The study at Langley was undertaken to resolve the pitch up anomaly observed during the entry of the first flight of the Shuttle Orbiter. The present paper will focus on real gas effects on aerodynamics and not on heating. The facilities used at NASA Langley were the 15-in. Mach 6, the 20-in, Mach 6, the 31-in. Mach 10 and the 20-in. Mach 6 CF4 facility. The paper focuses on the high Mach, high altitude portion of the first entry of the Shuttle where the vehicle exhibited a nose-up pitching moment relative to pre-flight prediction of (Delta C(sub m)) = 0.03. In order to study the relative contribution of compressibility, viscous interaction and real gas effects on basic body pitching moment and flap efficiency, an experimental study was undertaken to examine the effects of Mach, Reynolds and ratio of specific heats at NASA. At high Mach, a decrease of gamma occurs in the shock layer due to high temperature effects. The primary effect of this lower specific heat ratio is a decrease of the pressure on the aft windward expansion surface of the Orbiter causing the nose-up pitching moment. Testing in the heavy gas, Mach 6 CF4 tunnel, gave a good simulation of high temperature effects. The facilities used at ESA were the lm Mach 10 at ONERA Modane, the 0.7 m hot shot F4 at ONERA Le Fauga and the 0.88 m piston driven shock tube HEG at DLR Goettingen. Encouraging good force measurements were obtained in the F4 facility on the Orbiter configuration. Testing of the same model in the perfect gas Mach 10 S4 Modane facility was performed so as to have "reference" conditions. When one compares the F4 and S4 test results, the data suggests that the Orbiter "pitch up" is due to real gas effects. In addition, pressure measurements, performed on the aft portion of the windward side of the Halis configuration in HEG and F4, confirm that the pitch up is mainly attributed to a reduction of pressure due to a local decrease in gamma.

Description: This paper presents new trends in Particle Image Velocimetry and practical aspects relevant to the application of the technique to large scale wind tunnel testing. The various problems and their solutions to the operation of PIV in large scale wind tunnels are discussed. Application of the technique in mapping complex flows are also presented.

Description: A 3D CFD code, SCISEAL, has been developed and validated. Its capabilities include cylindrical seals, and it is employed on labyrinth seals, rim seals, and disc cavities. State-of-the-art numerical methods include colocated grids, high-order differencing, and turbulence models which account for wall roughness. SCISEAL computes efficient solutions for complicated flow geometries and seal-specific capabilities (rotor loads, torques, etc.).

Description: Pressure and temperature sensitive paints have been utilized for the measurement of blade surface pressure and temperature distributions in a high speed axial compressor and an Allied Signal F109 gas turbine engine. Alternate blades were painted with temperature sensitive paints and then pressure sensitive paint. This combination allows temperature distributions to be accounted for when determining the blade suction surface pressure distribution. Measurements were taken and pressure maps on the suction surface of a blade were obtained over a range of rotational speeds. Pressure maps of the suction surface show-strong shock waves at the higher speeds.

Description: Optical pressure measurements have been made on a NACA 0012 airfoil coated with Pressure Sensitive Paint (PSP) at very low flow speeds (less than 50 m/s). Angle of attack was limited to 5 deg. for most measurements. Effects of temperature gradients and mis-registration errors on PSP response have been established and minimized. By reducing measurement error caused by these effects. PSP sensitivity has been enhanced. Acceptable aerodynamic data at flow speeds down to 20 m/s have been obtained and valid pressure paint response was observed down to 10 m/s. Measurement errors (in terms of pressure and pressure coefficient) using PSP with pressure taps as a reference are provided for the range of flow speeds from 50 m/s to 10 m/s.

Description: The oil film interferometry skin-friction technique is described and applied to flows in some of the NASA Ames large wind tunnel facilities. Various schemes for applying the technique are discussed. Results are shown for tests in several wind tunnels which illustrate the oil film's ability to measure a variety of flow features such as shock waves separation, and 3D flow.

Description: Only recently have large amounts of model deformation data been acquired in NASA wind tunnels. This acquisition of model deformation data was made possible by the development of an automated video photogrammetric system to measure the changes in wing twist and bending under aerodynamic load. The measurement technique is based upon a single view photogrammetric determination of two dimensional coordinates of wing targets with a fixed third dimensional coordinate, namely the spanwise location. A major consideration in the development of the measurement system was that use of the technique must not appreciably reduce wind tunnel productivity. The measurement technique has been used successfully for a number of tests at four large production wind tunnels at NASA and a dedicated system is nearing completion for a fifth facility. These facilities are the National Transonic Facility, the Transonic Dynamics Tunnel, and the Unitary Plan Wind Tunnel at NASA Langley, and the 12-FT Pressure Tunnel at NASA Ames. A dedicated system for the Langley 16-Foot Transonic Tunnel is scheduled to be used for the first time for a test in September. The advantages, limitations, and strategy of the technique as currently used in NASA wind tunnels are presented. Model deformation data are presented which illustrate the value of these measurements. Plans for further enhancements to the technique are presented.

Description: The test capabilities of the Stability Wind Tunnel of the Virginia Polytechnic Institute and State University are described, and calibrations for curved and rolling flow techniques are given. Oscillatory snaking tests to determine pure yawing derivatives are considered. Representative aerodynamic data obtained for a current fighter configuration using the curved and rolling flow techniques are presented. The application of dynamic derivatives obtained in such tests to the analysis of airplane motions in general, and to high angle of attack flight conditions in particular, is discussed.

Description: The flow in a hydrostatic pocket is described by a mathematical model that uses the three dimensional Navier-Stokes equations written in terms of the primary variables, u, v, w, and p. Using a conservative formulation, a finite volume multi-block method is applied through a collocated, body fitted grid. The flow is simulated in a shallow pocket with a depth/length ratio of 0.02. The flow structures obtained and described by the authors in their previous two dimensional models are made visible in their three dimensional aspect for the Couette flow. It has been found that the flow regimes formed central and secondary vortical cells with three dimensional corkscrew-like structures that lead the fluid on an outward bound path in the axial direction of the pocket. The position of the central vortical cell center is at the exit region of the capillary restrictor feedline. It has also been determined that a fluid turn around zone occupies all the upstream space between the floor of the pocket and the runner, thus preventing any flow exit through the upstream port. The corresponding pressure distribution under the shaft presented as well. It was clearly established that for the Couette dominated case the pressure varies significantly in the pocket in the circumferential direction, while its variation is less pronounced axially.

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: An experimental program was conducted at NASA Langley Research Center that included development and evaluation of an operational facility for wall drag measurement of potential scramjet fuel injection or wall cooling configurations. The facility consisted of a supersonic tunnel, with one wall composed of a series of interchangeable aluminum plates attached to an air bearing suspension system. The system was equipped with load cells that measured drag forces of 115 psia (793 kPa). This flow field contained a train of weak, unsteady, reflecting shock waves that were produced in the Mach 2 nozzle flows, the effect of reflecting shocks (which are to be expected in scramjet combustors) in internal flows has not previously been documented.

Description: The wall drag test tunnel at NASA Langley Research Center was used to evaluate simulated scramjet fuel injection into a wall cavity. In this tunnel, one wall consists of interchangeable aluminum plates attached to an air bearing suspension system. The plates were equipped with load cells to measure drag forces and static taps to determine pressure distributions. The plates were exposed to a Mach 2 air stream at a total pressure of 115 psia (793 kPa). This flow field contained a train of weak unsteady, reflecting shock waves that were produced in the nozzle assembly located upstream of the test section.

Description: The work to be described was performed at the NASA Langley UPWT (4-ft supersonic), test section #2, during 21-24 May 1996. The configuration being tested was the 1.675% Ref H controls model; test conditions were Ma = 2.40, Re = 3 million/ft. This was an exploration of a new technique, and it was not intended to provide definitive comparison of measured and computed skin friction results. It is, however, hoped that the experience gained will make such a rigorous comparison possible in the future.

Description: To summarize the significant highlights in this report: (1) Data quality, determined by multiple repeat runs performed on the TCA baseline configuration, and long-term repeatability, determined by comparing baseline Reference H data from this test to a previous test, have been shown to be good. (2) The longitudinal stability of the TCA is more non-linear than for the Reference H, and while it is similar at normal lift values, the TCA has considerably more pitch-up at higher lift. (3) Longitudinal control effectiveness of the TCA is similar to the Reference H and the ratio of elevator effectiveness to horizontal tail effectiveness is approximately 0.3. 4) The directional stability of the TCA is improved relative to Reference H at higher angles-of attack. The chine is effective for improving directional stability. (5) The directional control effectiveness 'of the TCA rudder is the same as that of the Reference H rudder at low angles-of-attack, after taking factors, such as number of rudder panels deflected and vertical tail volume into account. However, rudder effectiveness was shown to be reduced at higher angles-of-attack. (6) The lateral stability was shown to be reduced relative to the Reference H, which may be beneficial at low speeds for alleviating lateral control saturation. (7) Lateral control effectiveness for the TCA was shown to be similar to the Reference H for negative trailing-edge flap deflections and was reduced by approximately 25% for positive trailing-edge flap deflections.

Description: This paper gives the results of a grid study, a turbulence model study, and a Reynolds number effect study for transonic flows over a high-speed aircraft using the thin-layer, upwind, Navier-Stokes CFL3D code. The four turbulence models evaluated are the algebraic Baldwin-Lomax model with the Degani-Schiff modifications, the one-equation Baldwin-Barth model, the one-equation Spalart-Allmaras model, and Menter's two-equation Shear-Stress-Transport (SST) model. The flow conditions, which correspond to tests performed in the NASA Langley National Transonic Facility (NTF), are a Mach number of 0.90 and a Reynolds number of 30 million based on chord for a range of angle-of-attacks (1 degree to 10 degrees). For the Reynolds number effect study, Reynolds numbers of 10 and 80 million based on chord were also evaluated. Computed forces and surface pressures compare reasonably well with the experimental data for all four of the turbulence models. The Baldwin-Lomax model with the Degani-Schiff modifications and the one-equation Baldwin-Barth model show the best agreement with experiment overall. The Reynolds number effects are evaluated using the Baldwin-Lomax with the Degani-Schiff modifications and the Baldwin-Barth turbulence models. Five angles-of-attack were evaluated for the Reynolds number effect study at three different Reynolds numbers. More work is needed to determine the ability of CFL3D to accurately predict Reynolds number effects.

Description: The NASA High Speed Research (HSR) Program is intended to establish a technology base enabling industry development of an economically viable and environmentally acceptable second generation high speed civil transport (HSCT). The HSR program consists of work directed towards several broad technology areas, one of which is aerodynamic performance. The objective of the Configuration Aerodynamics task of the Aerodynamic Performance technology area is the development of aerodynamic drag reduction, stability and control, and propulsion airframe integration technologies required to support the HSCT development process. Towards this goal, computational and empirical based aerodynamic design tools are being developed, evaluated, and validated through ground based experimental testing. In addition, methods for ground to flight scaling are being developed and refined. Successful development of validated design and scaling methodologies will result in improved economy of operation for an HSCT and reduce uncertainty in full-scale flight predictions throughout the development process.

Description: It is not unusual when comparing CFD data to experimental data to find discrepancies between the results. Sometimes forces and moments compare well, while surface pressures do not, and vice versa. It is commonplace for the researcher to believe that the flow field has been accurately simulated when these types of measurements compare well. However, being able to routinely predict boundary layer transition and separated flows are not guaranteed. In fact accurate simulation of these types of flow physics has been a challenge to the CFD community. In order to improve Navier-Stokes predictions for complex vortical flow fields, more detailed information about the flow physics is necessary. Unfortunately, the many wind-tunnel tests performed in Langley's NTF and 14x22 facilities as well as in the Ames' 12 ft. Tunnel provided little information about the detailed flow physics, and no priority was given to obtaining any CFD measurements. Using the latest experimental techniques, this information can and should be obtained for present and future use.

Description: This paper focuses on the parallel computation of aerodynamic derivatives via automatic differentiation of the Euler/Navier-Stokes solver CFL3D. The comparison with derivatives obtained by finite differences is presented and the scaling of the time required to obtain the derivatives relative to the number of processors employed for the computation is shown. Finally, the derivative computations are coupled with an optimizer and surface/volume grid deformation tools to perform an optimization to reduce the drag of a three-dimensional wing.

Description: This paper provides an overview of the NASA High-Speed Research (HSR) Program dedicated to establishing the technology foundation to support the US transport industry's decision for an environmentally acceptable, economically viable 300 passenger, 5000 n.mi., Mach 2.4 aircraft. The HSR program, begun in 1990, is supported by a team of US aerospace companies. The international economic stakes are high. The projected market for more than 500 High-Speed Civil Transport (HSCT) airplanes introduced between the years 2000 and 2015 translates to more than $200 billion in aircraft sales, and the potential of 140,000 new jobs. The paper addresses the history of supersonic commercial air transportation beginning with the Concorde and TU-144 developments in the early 1960 time period. The technology goals for the HSR program are derived from market study results, projections on environmental requirements, and technical goals for each discipline area referenced to the design and operational features of the Concorde. Progress since the inception of the program is reviewed and a summary of some of the lessons learned will be highlighted. An outline is presented of the remaining technological challenges. Emphasis in this paper will be on the traditional aeronautical technologies that lead to higher performance to ensure economic viability. Specific discussion will center around aerodynamic performance, flight deck research, materials and structures development and propulsion systems. The environmental barriers to the HSCT and that part of the HSR program that addresses those technologies are reviewed and assessed in a companion paper.

Description: Current parallel computational approaches involve distributed and shared memory paradigms. In the distributed memory paradigm, each processor has its own independent memory. Message passing typically uses a function library such as MPI or PVM. In the shared memory paradigm, such as that used on the SGI Origin 2000 machine, compiler directives are used to instruct the compiler to schedule multiple threads to perform calculations. In this paradigm, it must be assured that processors (threads) do not simultaneously access regions of memory in such away that errors would occur. This paper utilizes the latest version of the SGI MPI function library to combine the two parallelization paradigms to perform aerodynamic shape optimization of a generic wing/body.

Description: Minuscule amounts (e.g., 1 part in 10,000) of a surface-active material in a liquid can drastically affect the surface behavior of the liquid, influencing how the material flows and mixes with other liquid and solid materials. In many respects, the science of surfactants has been empirical, with trial and error dominating over the ability to predict how surfactant type and concentration influence surface behavior. A program for the modeling of surfactant behavior has been established at Yale. This program combines experimental work performed both on the ground and in space, and theoretical and numerical modeling. By levitating a drop of liquid in air, away from solid container surfaces, and by manipulating the drop with acoustic radiation forces, we have been able to establish idealized conditions for surface behavior studies. The primary experiments involve the study of the free oscillations of initially deformed drops. In STS-73, the USML-2 mission of the Space Shuttle, we performed the following measurements: 1) the oscillation of a spherical drop in its quadrupole mode; 2) the oscillation of a drop about a deformed (oblate) shape; 3) the slow static squeezing of the drop from spherical to nearly flat; and 4) the superoscillations of drops when the radiation forces maintaining the drop in a flattened state are suddenly reduced. Analytic and numerical studies have enabled us to understand the physics of these oscillations and to extract material properties such as the dynamic surface tension and the surface viscosities (shear and dilatational). The relation to ground-based studies is essential, because the knowledge and understanding gleaned from our space studies enable us to interpret ground-based data.

Description: Progress in aerodynamics over the past 50 years has been evidenced by the development of increasingly sophisticated and efficient flight vehicles throughout the flight spectrum. Advances have generally arisen in an evolutionary manner from experience gained in wind tunnel testing, flight testing, and improvements in analytical and computational capabilities. As a result of this evolutionary development, both military and commercial vehicles operate at a relatively high efficiency level. This observation plus the fact that airplanes have not changed appreciably in outward appearance over recent years has led some skeptics to conclude incorrectly that aerodynamics is a mature technology, with little to be gained from further developments in the field. It is of interest to note that progress in aerodynamics has occurred without a thorough understanding of the fundamental physics of flow, turbulence, vortex dynamics, and separated flow, for example. The present understanding of transition, turbulence, and boundary layer separation is actually very limited. However, these fundamental flow phenomena provide the key to reducing the viscous drag of aircraft. Drag reduction provides the greatest potential for increased flight efficiency from the standpoint of both saving energy and maximizing performance. Recent advances have led to innovative concepts for reducing turbulent friction drag by modifying the turbulent structure within the boundary layer. Further advances in this basic area should lead to methods for reducing skin friction drag significantly. The current challenges for military aircraft open entirely new fields of investigation for the aerodynamicist. The ability through very high speed information processing technology to totally integrate the flight and propulsion controls can permit an aircraft to fly with "complete abandon," avoiding departure, buffet, and other undesirable characteristics. To utilize these new control concepts, complex aerodynamic phenomena will have to be understood, predicted, and controlled. Current requirements for military aircraft include configuration optimization through a widened envelope from subsonic to supersonic and from low to high angles of attack. This task is further complicated by requirements for control of observables. These challenging new designs do not have the luxury of a large experimental data base from which to optimize for various parameter combinations. Consequently, there exists a strong need for better techniques, both experimental and computational, to permit design optimization in a complete sense.

Description: A methodology for the simulation of turbine cooling flows is being developed. The methodology seeks to combine numerical techniques that optimize both accuracy and computational efficiency. Key components of the methodology include the use of multiblock grid systems for modeling complex geometries, and multigrid convergence acceleration for enhancing computational efficiency in highly resolved fluid flow simulations. The use of the methodology has been demonstrated in several turbo machinery flow and heat transfer studies. Ongoing and future work involves implementing additional turbulence models, improving computational efficiency, adding AMR.

Description: The following are conclusions and recomendations from the study. Primary wake effect is linear reduction in (eta) with St. Secondary wake effect is skewing of suction/pressure side cooling. Steady computations match experimental Nu, but overpredict (eta). Unsteady computations elucidate wake/film interaction. Model may be used to estimate wake passing effect. Need boundary layer and full stage experiments. Need resolved film hole and full stage unsteady computations. Need validated turbulence models for film cooling.

Description: The following topics are discussed: (1) The Transient Liquid-Crystal Heat-Transfer Technique; (2) 2-D Film-Cooling Heat-Transfer on an AlliedSignal Vane; and (3) Effects of Tab Vortex Generators on Surface Heat Transfer. Downstream of a Jet in Crossflow.

Description: A tab placed on the leeward side of the nozzle was expected to increase jet penetration into the cross-flow. An experiment at UTRC showed insignificant effect. The primary objective of the present study was to confirm and explain the ineffectiveness. The overall approach of the study was to conduct experiments in a low-speed wind tunnel and to conduct hot-wire measurements for mean velocity and streamwise vorticity fields.

Description: Testing is currently in progress to acquire simple geometry surface heat transfer data for internal channels with trips and bleed holes which can be used in the development and validation of models. The transient liquid crystal technique is used on a simple multipass model with rectangular channels and normal ribs. Normal bleed holes are located on the floor of the model in the first channel. Each hole is attached to a flow meter, allowing various bleed flow rates to simulate external pressures on the blade.

Description: Reasonable heat transfer prediction can be achieved in complex geometries. Multi-block grid allows efficient placement of grid points, and efficient use of computer resources. Wilcox k-(omega) turbulence model predicts heat transfer well, and has good numerical behavior.

Description: The flow through a plane asymmetric diffuser is a good test case for assessing the capability of LES since it contains features such as large scale unsteady separation and strong intermittency which are difficult to capture using conventional modeling approaches. Previous attempts to simulate this flow (Kaltenbach, 1994) have significantly underpredicted the extent of separation. The objective of the present research is to understand why the previous simulations did not predict the flow separation correctly. This study focuses on mesh refinement and matching of the inlet velocity profile. In order to perform this study, the flow solver of Kaltenbach (1994) was modified to increase its accuracy and efficiency. The improved algorithm allows for better resolution at affordable CPU cost.

Description: The objective of the present research is to assess the usefulness of large-eddy simulation (LES) methodology for flows in complex geometries. Flow past a circular cylinder has been calculated using a central-difference based solver, and the results have been compared to those obtained by a solver that employs higher-order upwind biased schemes (Beaudan & Moin, 1994). This comparison allows us to assess the suitability of these schemes for LES in complex geometry flows.

Description: Accurate SGS models must be capable of correctly representing the energy transfer between GS and SGS. Recent direct assessment of the energy transfer carried out using direct numerical simulation (DNS) data for wall-bounded flows revealed that the energy exchange is not unidirectional. Although GS kinetic energy is transferred to the SGS (forward scatter (F-scatter) on average, SGS energy is also transferred to the GS. The latter energy exchange (backward scatter (B-scatter) is very significant, i.e., the local energy exchange can be backward nearly as often as forward and the local rate of B-scatter is considerably higher than the net rate of energy dissipation.

Description: An efficient scheme for the direct numerical simulation of 3D transitional and developed turbulent flow is presented. Explicit and implicit time integration schemes for the compressible Navier-Stokes equations are compared. The nonlinear system resulting from the implicit time discretization is solved with an iterative method and accelerated by the application of a multigrid technique. Since we use central spatial discretizations and no artificial dissipation is added to the equations, the smoothing method is less effective than in the more traditional use of multigrid in steady-state calculations. Therefore, a special prolongation method is needed in order to obtain an effective multigrid method. This simulation scheme was studied in detail for compressible flow over a flat plate. In the laminar regime and in the first stages of turbulent flow the implicit method provides a speed-up of a factor 2 relative to the explicit method on a relatively coarse grid. At increased resolution this speed-up is enhanced correspondingly.

Description: The role of unsteady flow processes in establishing the performance of axial flow turbomachinery was investigated. The development of the flow model, as defined by the time average flow equations associated with the passage of the blade row embedded in a multistage configuration, is presented. The mechanisms for unsteady momentum and energy transport is outlined. The modeling of the unsteady momentum and energy transport are discussed. The procedure for simulating unsteady multistage turbomachinery flows is described.

Description: A mathematical model for closing or mathematically completing the system of equations is proposed. The model describes the time average flow field through the blade passages of multistage turbomachinery. These average-passage equation systems govern a conceptual model useful in turbomachinery aerodynamic design and analysis. The closure model was developed to insure a consistency between these equations and the axisymmetric through-flow equations. The closure model was incorporated into a calculation code for use in the simulation of the flow field about a high-speed counter rotating propeller and a high-speed fan stage.

Description: The primary objectives of this study were to expand the data base showing the effects of LE radius distribution and corresponding sensitivity to Rn at subsonic and transonic conditions, and to assess the predictive capability of CFD for these effects. Several key elements led to the initiation of this project: 1) the necessity of meeting multipoint design requirements to enable a viable HSCT, 2) the demonstration that blunt supersonic leading-edges can be associated with performance gain at supersonic speeds , and 3) limited data. A test of a modified Reference H model with the TCA planform and 2 LE radius distributions was performed in the NTF, in addition to Navier-Stokes analysis for an additional 3 LE radius distributions. Results indicate that there is a tremendous potential to improve high-lift performance through the use of a blunt LE across the span given an integrated, fully optimized design, and that low Rn data alone is probably not sufficient to demonstrate the benefit.

Description: This paper presents results of a study which attempted to provide some understanding of the relationship between skin friction drag estimates produced by flat plate methods and those produced by Navier-Stokes computations. A brief introduction is followed by analysis, including a flat plate grid study, analysis of the wing flow, an analysis of the fuselage flow. Other results of interest are then presented, including turbulence model sensitivities, and brief analysis of other configurations.

Description: Efforts towards understanding boundary layer transition characteristics on a High Speed Civil Transport (HSCT)-class configuration in the National Transonic Facility (NTF) are ongoing. The majority of the High Speed Research (HSR) data base in the NTF has free transition on the wing, even at low Reynolds numbers (Rn) attainable in conventional facilities. Limited data has been obtained and is described herein showing the effects of a conventional, Braslow method based wing boundary-layer trip on drag. Comparisons are made using force data polars and surface flow visualization at selected angles-of-attack and Mach number. Minimum drag data obtained in this study suggest that boundary layer transition occurred very near the wing leading edge by a chord Rn of 30 million. Sublimating chemicals were used in the air mode of operation only at low Rn and low angles-of-attack with no flap deflections; sublimation results suggest that the forebody and outboard wing panel are the only regions with significant laminar flow. The process and issues related to the sublimating chemical technique as applied in the NTF are discussed. Beyond the existing experience, status of efforts to develop a production transition detection system applicable to both air and cryogenic nitrogen environments is presented.

Description: Model deformation measurement techniques have been investigated and developed at NASA's Langley Research Center. The current technique is based upon a single video camera photogrammetric determination of two dimensional coordinates of wing targets with a fixed (and known) third dimensional coordinate, namely the spanwise location. Variations of this technique have been used to measure wing twist and bending at a few selected spanwise locations near the wing tip on HSR models at the National Transonic Facility, the Transonic Dynamics Tunnel, and the Unitary Plan Wind Tunnel. Automated measurements have been made at both the Transonic Dynamics Tunnel and at Unitary Plan Wind Tunnel during the past year. Automated measurements were made for the first time at the NTF during the recently completed HSR Reference H Test 78 in early 1996. A major problem in automation for the NTF has been the need for high contrast targets which do not exceed the stringent surface finish requirements. The advantages and limitations (including targeting) of the technique as well as the rationale for selection of this particular technique are discussed. Wing twist examples from the HSR Reference H model are presented to illustrate the run-to-run and test-to-test repeatability of the technique in air mode at the NTF. Examples of wing twist in cryogenic nitrogen mode at the NTF are also presented.

Description: To develop full scale flight performance predictions an understanding of Reynolds number effects on HSCT-class configurations is essential. A wind tunnel database utilizing a 2.2% scale Reference H model in NASA Langley Research Centers National Transonic Facility is being developed to assess these Reynolds number effects. In developing this database temperature and aeroelastic corrections to the wind tunnel data have been identified and are being analyzed. Once final corrections have been developed and applied, then pure Reynolds number effects can be determined. In addition, final corrections will yield the data required for CFD validation at q = 0. Presented in this report are the results of seven tests involving the wing/body configuration. This includes summaries of data acquired in these tests, uncorrected Reynolds number effects, and temperature and aeroelastic corrections. The data presented herein illustrates the successes achieved to date as well as the challenges that will be faced in obtaining full scale flight performance predictions.

Description: Experience with afterbody closure effects and accompanying test techniques issues on a High Speed Civil Transport (HSCT)-class configuration is described. An experimental data base has been developed which includes force, moment, and surface pressure data for the High Speed Research (HSR) Reference H configuration with a closed afterbody at subsonic and transonic speeds, and with a cylindrical afterbody at transonic and supersonic speeds. A supporting computational study has been performed using the USM3D unstructured Euler solver for the purposes of computational fluid dynamics (CFD) method assessment and model support system interference assessment with a focus on lower blade mount effects on longitudinal data at transonic speeds. Test technique issues related to a lower blade sting mount strategy are described based on experience in the National Transonic Facility (NTF). The assessment and application of the USM3D code to the afterbody/sting interference problem is discussed. Finally, status and plans to address critical test technique issues and for continuation of the computational study are presented.

Description: The Boeing Reference H configuration was tested in the NASA Ames 9x7 Supersonic Wind Tunnel. A simulated unstarted inlet was evaluated as well as the aerodynamic performance of the configuration with and without nacelle and diverter components. These experimental results were compared with computational results from the unstructured grid Euler flow solver AIRPLANE. The comparisons between computational and experimental results were good, and demonstrated that the Euler code is capable of efficiently and accurately predicting the changes in the aerodynamic coefficients associated with inlet unstart and the effects of the nacelle and diverter components.

Description: This presentation will describe the organization and conduct of the workshops, list the topics discussed, and conclude with a more-detailed examination of a related set of issues dear to the presenters heart. Because the current HSCT configuration is expected to have (mostly) turbulent flow over the wings, and because current CFD predictions assume fully-turbulent flow, the wind tunnel testing to date has attempted to duplicate this condition at the lower Reynolds numbers attainable on the ground. This frequently requires some form of artificial boundary layer trip to induce transition near the wing's leading edge. But this innocent-sounding goal leads to a number of complications, and it is not clear that present-day testing technology is adequate to the task. An description of some of the difficulties, and work underway to address them, forms the "Results" section of this talk. Additional results of the testing workshop will be covered in presentations by other team members.

Description: This paper presents The Propulsion Airframe Integration Advisory report in viewgraph form. The approach of the advisory group is to identify and prioritize technology elements (1.0 Inlet Issues, 2.0 Nozzle Issues, 3.0 Nacelle Design, and 4.0 Airframe Integration).

Description: Preliminary human acceptability studies of sonic booms indicate that supersonic flight is unlikely to be acceptable even at noise levels significantly below 1994 low boom designs (reference 1, p. 288). Further, these low boom designs represent considerable changes to baseline configurations, and changes translate into additional effort and uncertain structural weight penalties that may provide no annoyance benefit, increasing the risk of including low boom technology. Since over land sonic boom designs were so risky (and yet the acceptability studies highlight how annoying sonic booms are), boom softening studies were undertaken to reduce the boom of baseline configurations using minor modifications that would not significantly change the designs. The goal of this work is to reduce boom levels over water. Even though Concorde over water boom has not been found to have any adverse environmental impact, boom levels for baseline HSCT designs are 50% higher in overpressure than the Concorde (due to a doubling in configuration weight with only a 50% increase in length),

Description: Officially, the Tu-144 was the first supersonic-cruise, passenger-carrying aircraft to enter commercial service. Design, construction, and testing were carried out by the Soviet Union, flight certification was by the Soviet Union, and the only regular passenger flights were scheduled and flown across the territory of the Soviet Union. Although it was not introduced to international passenger service, there were many significant engineering accomplishments achieved in the design, production, and flight of this aircraft. Development of the aircraft began with a prototype stage. Systematic testing and redesign led to a production aircraft in discrete stages that measurably improved the performance of the aircraft from the starting concept to final aircraft certification. It flew in competition with the English-French Concorde for a short time, but was withdrawn from national commercial service due to a lack of interest by airlines outside the Soviet Union. NASA became interested in the Tu- 144 aircraft when it was offered for use as a flying "testbed" in the study of operating characteristics of a supersonic-cruise commercial airplane. Since it had been in supersonic-cruise service, the Tu- 144 had operational characteris'tics similar to those anticipated in the conceptual aircraft designs being studied by the United States aircraft companies. In addition to the other operational tests being conducted on the Tu-144 aircraft, it was proposed that two sets of sonic-boom pressure signature measurements be made. The first set would be made on the ground, using techniques and devices similar to those in reference I and many other subsequent studies. A second set would be made in the air with an instrumented aircraft flying close under the Tu-144 in supersonic flight. Such in-flight measurements would require pressure gages that were capable of accurately recording the flow-field overpressures generated by the Tu- 144 at relatively close distances under the vehicle. Therefore, an analysis of the Tu-144 was made to obtain predictions of pressure signature shape and shock strengths at cruise conditions so that the range and characteristics of the required pressure gages could be determined well in advance of the tests. Cancellation of the sonic-boom signature measurement part of the tests removed the need for these pressure gages. Since CFD methods would be used to analyze the aerodynamic performance of the Tu-144 and make similar pressure signature predictions, the relatively quick and simple Whitham-theory pressure signature predictions presented in this paper could be used for comparisons. Pressure signature predictions of sonic-boom disturbances from the Tu- 144 aircraft were obtained from geometry derived from a three-view description of the production aircraft. The geometry was used to calculate aerodynamic performance characteristics at supersonic-cruise conditions. These characteristics and Whitham/Walkden sonic-boom theory were employed to obtain F-functions and flow-field pressure signature predictions at a Mach number of 2.2, at a cruise altitude of 61000 feet, and at a cruise weight of 350000 pounds.

Description: This document contains the details of the thermal analysis of the X-38 aft fin during re-entry. This analysis was performed in order to calculate temperature response of the aft fin components. This would be provided as input to a structural analysis and would also define the operating environment for the electromechanical actuator (EMA). The calculated structural temperature response would verify the performance of the thermal protection system (TPS). The geometric representation of the aft fin was derived from an I-DEAS finite element model that was used for structural analysis. The thermal mass network model was derived from the geometric representation.

Description: While turbulence models are being developed which show promise for simulating the transition region on a turbine blade or vane, it is believed that the best approach with the greatest potential for practical use is the use of models which incorporate the physics of turbulent spots present in the transition region. This type of modeling results in the prediction of transition region intermittency which when incorporated in turbulence models give a good to excellent prediction of the transition region heat transfer. Some models are presented which show how turbulent spot characteristics and behavior can be employed to predict the effect of pressure gradient and Mach number on the transition region. The models predict the spot formation rate which is needed, in addition to the transition onset location, in the Narasimha concentrated breakdown intermittency equation. A simplified approach is taken for modeling turbulent spot growth and interaction in the transition region which utilizes the turbulent spot variables governing transition length and spot generation rate. The models are expressed in terms of spot spreading angle, dimensionless spot velocity, dimensionless spot area, disturbance frequency and Mach number. The models are used in conjunction with a computer code to predict the effects of pressure gradient and Mach number on the transition region and compared with VKI experimental turbine data.

Description: NASA and FAA initiated a program in 1994 to develop methods of setting spacings for landing aircraft by incorporating information on the real-time behavior of aircraft wake vortices. The current wake separation standards were developed in the 1970's when there was relatively light airport traffic and a logical break point by which to categorize aircraft. Today's continuum of aircraft sizes and increased airport packing densities have created a need for re-evaluation of wake separation standards. The goals of this effort are to ensure that separation standards are adequate for safety and to reduce aircraft spacing for higher airport capacity. Of particular interest are the different requirements for landing under visual flight conditions and instrument flight conditions. Over the years, greater spacings have been established for instrument flight than are allowed for visual flight conditions. Preliminary studies indicate that the airline industry would save considerable money and incur fewer passenger delays if a dynamic spacing system could reduce separations at major hubs during inclement weather to the levels routinely achieved under visual flight conditions. The sensor described herein may become part of this dynamic spacing system known as the "Aircraft VOrtex Spacing System" (AVOSS) that will interface with a future air traffic control system. AVOSS will use vortex behavioral models and short-term weather prediction models in order to predict vortex behavior sufficiently into the future to allow dynamic separation standards to be generated. The wake vortex sensor will periodically provide data to validate AVOSS predictions. Feasibility of measuring wake vortices using a lidar was first demonstrated using a continuous wave (CW) system from NASA Marshall Space Flight Sensor and tested at the Volpe National Transportation Systems Center's wake vortex test site at JFK International Airport. Other applications of CW lidar for wake vortex measurement have been made more recently, including a system developed by the MIT Lincoln Laboratory. This lidar has been used for detailed measurements of wake vortex velocities in support of wake vortex model validation. The first measurements of wake vortices using a pulsed, lidar were made by Coherent Technologies, Inc. (CTI) using a 2 micron solid-state, flashlamp-pumped system operating at 5 Hz. This system was first deployed at Denver's Stapleton Airport. Pulsed lidar has been selected as the baseline technology for an operational sensor due to its longer range capability.

Description: A research is performed: to define wake non-encounter & hazard, to provide requirements for sensors, and to obtain input from the user community. This research includes: validating wake encounter simulation models, establishing a metric to quantify the upset potential of a wake encounter, applying hazard metric and simulation models to the commercial fleet for development of candidate acceptable encounter limits, and applying technology to near term problems to evaluate current status of technology. The following lessons are learned from this project: technology is not adequate to determine absolute spacing requirements; time, not distance, determines the duration of the wake hazard; Optimum standards depend on the traffic; Wing span is an important parameter for characterizing both generator and follower; and Short span "biz jets" are easily rolled.

Description: The overall approach should be to: (1) Seek simplest, sufficiently robust, integrated ground based sensor systems (wakes and weather) for AVOSS; (2) Expand all sensor performance cross-comparisons and data mergings in on-going field deployments; and (3) Achieve maximal cost effectiveness through hardware/info sharing. An effective team is in place to accomplish the above tasks.

Description: In the overview, a description of the LaRC trailer facility, lasers and transceivers, scanners, data systems and deployment are presented.

Description: Included in the overview is a discussion of the 1.5 micron laser specifications, eye safety and cost, scan rates, pulselength, range capability issues, Raman beam cleanup, receiver layout, and the real-time processor and display.

Description: Results from parametric runs using two-dimensional TASS are presented. First, a set of experiments are presented that examine the sensitivity of the aircraft initiation height for an "in ground effect" case with weak crosswind. Interaction between the ground and the wake vortex produces an oscillatory rebound whose phase and amplitude are a function of the generation height. A second set of experiments are presented which examine the influence on crosswind shear. Shear layers, such as may be found between the nocturnal stable layer and the residual layer, can act to deflect vortices upward. Further investigation reveals that the second derivative of the crosswind can differentially reduce the descent speed of each member of a vortex pair, causing tilting of the vortex pair. If sufficiently large, the second derivative of crosswind can deflect the vortex pair upwards, with the sign of the second derivative determining which of the two vortices rises to a higher altitude. Linear shear, on the other hand, caused no change in the descent speed of the vortices; thus having no effect on the orientation of the vortices. Observed and model data from an actual case are presented in support of the conclusion regarding the influence of shear on rising vortices.

Description: Crow instability can develop in most atmospheric turbulence levels, however, the ring vortices may not form in extremely strong turbulence cases due to strong dissipation of the vortices. It appears that strong turbulence tends to accelerate the occurrences of Crow instability. The wavelength of the most unstable mode is estimated to be about 5b(sub 0), which is less than the theoretical value of 8.6b(sub 0) (Crow, 1970) and may be due to limited domain size and highly nonlinear turbulent flow characteristics. Three-dimensional turbulence can decay wake vortices more rapidly. Axial velocity may be developed by vertical distortion of a vortex pair due to Crow instability or large turbulent eddy motion. More experiments with various non-dimensional turbulence levels are necessary to get useful statistics of wake vortex behavior due to turbulence. Need to investigate larger turbulence length scale effects by enlarging domain size or using grid nesting.

Description: The vortices produced by an aircraft in flight are a complex phenomena created from a 'sheet of vorticity' leaving the trailing edge of the aircraft surfaces. This sheet tends to roll-up into two counter-rotating vortices. After a few spans downstream of the aircraft, the roll-up process is complete and the vortex pair may be characterized in a simple manner for modeling purposes. Our research will focus on what happens to these post roll-up vortices in the vicinity of an airport terminal. As the aircraft wake vortices descend, they are transported by the air mass which they are embedded and are decayed by both internal and external processes. In the vicinity of the airport, these external influences are usually due to planetary boundary layer (PBL) turbulence. Using large-eddy simulation (LES), one may simulate a variety of PBL conditions. In the LES method, turbulence is generated in the PBL as a response to surface heat flux, horizontal pressure gradient, wind shear, and/or stratification, and may produce convective or unstably stratified, neutral, or stably stratified PBL's. Each of these PBL types can occur during a typical diurnal cycle of the PBL. Thus it is important to be able to characterize these conditions with the LES method. Once this turbulent environment has been generated, a vortex pair will be introduced and the interactions are observed. The objective is to be able to quantify the PBL turbulence vortex interaction and be able to draw some conclusions of vortex behavior from the various scale interactions. This research is ongoing, and we will focus on what has been accomplished to date and the future direction of this research. We will discuss the model being used, show results that validate its use in the PBL, and present a nested-grid method proposed to analyze the entire PBL and vortex pair simultaneously.

Description: The purpose of the modelling effort at NASA Langley, including goals, is outlined in this presentation. Included, is a description of the numerical model that is used for the NASA wake vortex modeling effort and the approach that is taken in order to achieve the stated goals. Also shown are: 1) a demonstration of using the model in a fog environment; 2) preliminary results from a 3-D simulation in a nonturbulent and thermally-stable environment with comparison to a comparable 2-D simulation of the same event; and 3) several validation cases from the Idaho-Falls and Memphis field studies where results from the 2-D version of the model are compared with Lidar and tower data.

Description: The objective of this study is to provide careful qualification and detailed measurements in a re-creation of the Obi experiment. The work will include extensive documentation of the flow two-dimensionality and detailed measurements required for testing of flow computations. Also important to this study is the close interaction of the experimental and computational groups to improve the utility of the data obtained and the accuracy of computation.

Description: The control of vortical flows is gaining significance in the design of aeronautical and marine structures. While passive devices have been used effectively in the past, active control strategies have the potential of allowing a leap in the performance of future configurations. The efficiency of control schemes is strongly dependent on the development of accurate flow models that can be devised using information that is available not only from numerical solutions of the governing Navier-Stokes equations but also can be measured experimentally. In that context it is desirable to construct adaptive control schemes using information that can be measured at the wall.

Description: Stability of the elliptic vortex attracted interest in the past decade. Cambon (1982), and Cambon, Teissedre and Jeandel (1985) have studied the stability of such flows with spatially uniform velocity gradient, and have provided RDT solutions for a wide range of the parameter S/Omega (where the strain rate S and the vorticity 2 Omega define the velocity gradient matrix). The range studied included those of hyperbolic streamlines (strain dominated, S/Omega is greater than 1), linear streamlines (simple shear, S/Omega = 1), and elliptical streamlines (vorticity dominated, S/Omega is less than 1). The latter class has more recently attracted interest and several studies appeared (Pierrehumbert 1986, Bayly 1986, Craik and Criminale 1986). These studies will be collectively referred to as PBCC.

Description: The objective is the development of a new theory which enables the algorithmic computation of all self-similar mean velocity profiles. The theory is based on Liegroup analysis and unifies a large set of self-similar solutions for the mean velocity of stationary parallel turbulent shear flows. The results include the logarithmic law of the wall, an algebraic law, the viscous sublayer, the linear region in the middle of a Couette flow and in the middle of a rotating channel flow, and a new exponential mean velocity profile not previously reported. Experimental results taken in the outer parts of a high Reynolds number flat-plate boundary layer, strongly support the exponential profile. From experimental as well as from DNS data of a turbulent channel flow the algebraic scaling law could be confirmed in both the center region and in the near wall region. In the case of the logarithmic law of the wall, the scaling with the wall distance arises as a result of the analysis and has not been assumed in the derivation. The crucial part of the derivation of all the different mean velocity profiles is to consider the invariance of the equation for the velocity fluctuations at the same time as the invariance of the equation for the velocity product equations. The latter is the dyad product of the velocity fluctuations with the equation for the velocity fluctuations. It has been proven that all the invariant solutions are also consistent with similarity of all velocity moment equations up to any arbitrary order.

Description: Nonpremixed combustion is the process whereby fuel and oxidizer species, which are each nonflammable in isolation, concurrently mix to burn a flammable mixture, and chemically react in the flammable mixture. In cases of practical industrial interest, the bulk of nonpremixed combustion occurs in a turbulent mixing regime where enhanced mass transfer rates flow the maximum power density to be achieved in any given thermochemical device.

Description: Direct numerical simulation of turbulent reacting flows places extreme demands on computational resources. At the present time, simulations can be performed only for greatly simplified reaction systems and for very low Reynolds numbers. Direct simulation of more realistic cases occurring at higher Reynolds number and including multiple species and numerous chemical reactions will exceed available computational resources far into the future. Because of this, there is a clear need to develop the technique of large eddy simulation for reacting flows. Unfortunately this task is complicated by the fact that combustion arises from chemical reactions that occur at the smallest scales of the flow. Capturing the large-scale behavior without resolving the small-scale details is extremely difficult in combustion problems. Thus LES modeling for turbulent combustion encounters difficulties not present in modeling momentum transport, in which the main effect of the small scales is to provide dissipation. The difficulty is more pronounced in premixed combustion, where detailed chemistry plays an essential role in determining the flame speed (or overall burning rate); in nonpremixed combustion infinite rate chemistry can be assumed, eliminating the small scale features to a first approximation.

Description: Many industrial applications in such fields as aeronautical, mechanical, thermal, and environmental engineering involve complex turbulent flows containing global separations and subsequent reattachment zones. Accurate prediction of this phenomena is very important because separations influence the whole fluid flow and may have an even bigger impact on surface heat transfer. In particular, reattaching flows are known to be responsible for large local variations of the local wall heat transfer coefficient as well as modifying the overall heat transfer. For incompressible, non-buoyant situations, the fluid mechanics have to be accurately predicted in order to have a good resolution of the temperature field.

Description: In simple flows, where the mean deformation rates are mild and the turbulence has time to come to equilibrium with the mean flow, the Reynolds stresses are determined by the applied strain rate. Hence in these flows, it is often adequate to use an eddy-viscosity representation. The modern family of kappa-epsilon models has been very useful in predicting near equilibrium turbulent flows, where the rms deformation rate S is small compared to the reciprocal time scale of the turbulence (epsilon/kappa). In modern engineering applications, turbulence models are quite often required to predict flows with very rapid deformations (large S kappa/epsilon). In these flows, the structure takes some time to respond and eddy viscosity models are inadequate. The response of turbulence to rapid deformations is given by rapid distortion theory (RDT). Under RDT the nonlinear effects due to turbulence-turbulence interactions are neglected in the governing equations, but even when linearized in this fashion, the governing equations are unclosed at the one-point level due to the non-locality of the pressure fluctuations.

Description: Direct numerical simulation (DNS) was used to study homogeneous turbulence with a mean flow having elliptic streamlines. This flow combines the effects of rotation and strain on the turbulence. There are many important reasons for studying the elliptic streamline flow. This flow contains the effects of both rotation and strain and is therefore similar to the mean flow in a vortex strained in the plane perpendicular to its axis. Such flows provide insight into the fundamental vortical interactions within turbulence. A strained vortex also occurs in airplane wakes, in which each wingtip vortex induces a strain field on the other. The strain field can affect the stability of these vortices and thereby their turbulent structure downstream. The ability to understand and predict the turbulent structure of the vortices is important to the wake hazard problem which is of major concern for the safety of commercial aircraft.

Description: A combined experimental and computational effort to test the performance of turbulence models in prediction of three-dimensional separation is under way. The current phase pf the project is the acquisition of flow field data using laser Doppler velocimetry. Accomplishments for summer 1996 include development of the instrument and signal processing system and acquisition of mean velocity profiles.

Description: The model is a replica of a modified MAX-103 kit aircraft that our Parks College of St. Louis University Student Design Group built and modified from a tail wheel to a tricycle configuration. A model was tested in the Parks College low-speed wind tunnel. I hope to initiate flight-testing upon my second return to. St. Louis. The combined data using wind tunnel, water tunnel, RC, flight-testing and analytical results will be very valuable for assessing the correlation between the different methods of analyses, since at present it is almost impossible to accurately predict flight characteristics from anything but in-situ tests. Unfortunately, political/financial reasons dictated using a generic wing rather than a specific model in the NASA-DFRC water tunnel.

Description: The mean and phase averaged pressure and wall shear stress distributions were measured on the stator wall of a 50% eccentric annular seal which was whirling in a circular orbit at the same speed as the shaft rotation. The shear stresses were measured using flush mounted hot-film probes. Four different operating conditions were considered consisting of Reynolds numbers of 12,000 and 24,000 and Taylor numbers of 3,300 and 6,600. At each of the operating conditions the axial distribution (from Z/L = -0.2 to 1.2) of the mean pressure, shear stress magnitude, and shear stress direction on the stator wall were measured. Also measured were the phase averaged pressure and shear stress. These data were combined to calculate the force distributions along the seal length. Integration of the force distributions result in the net forces and moments generated by the pressure and shear stresses. The flow field inside the seal operating at a Reynolds number of 24,000 and a Taylor number of 6,600 has been measured using a 3-D laser Doppler anemometer system. Phase averaged wall pressure and wall shear stress are presented along with phase averaged mean velocity and turbulence kinetic energy distributions located 0.16c from the stator wall where c is the seal clearance. The relationships between the velocity, turbulence, wall pressure and wall shear stress are very complex and do not follow simple bulk flow predictions.

Description: Computational programs developed for the thermal analysis of tilting and flexure-pad hybrid bearings, and the unsteady flow and transient response of a point mass rotor supported on fluid film bearings are described. The motion of a cryogenic liquid on the thin film annular region of a fluid film bearing is described by a set of mass and momentum conservation, and energy transport equations for the turbulent bulk-flow velocities and pressure, and accompanied by thermophysical state equations for evaluation of the fluid material properties. Zeroth-order equations describe the fluid flow field for a journal static equilibrium position, while first-order (linear) equations govern the fluid flow for small amplitude-journal center translational motions. Solution to the zeroth-order flow field equations provides the bearing flow rate, load capacity, drag torque and temperature rise. Solution to the first-order equations determines the rotordynamic force coefficients due to journal radial motions.

Description: A hybrid boundary element finite volume method for unsteady transonic flow computation has been developed. In this method, the unsteady Euler equations in a moving frame of reference are solved in a small embedded domain (inner domain) around the airfoil using an implicit finite volume scheme. The unsteady full-potential equation, written in the same frame of reference and in the form of the Poisson equation. is solved in the outer domain using the integral equation boundary element method to provide the boundary conditions for the inner Euler domain. The solution procedure is a time-accurate stepping procedure, where the outer boundary conditions for the inner domain are updated using the integral equation -- boundary element solution over the outer domain. The method is applied to unsteady transonic flows around the NACA0012 airfoil undergoing pitching oscillation and ramp motion. The results are compared with those of an implicit Euler equation solver, which is used throughout a large computational domain, and experimental data.

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: To summarize the significant highlights in this report: (1) Data quality, determined by multiple repeat runs performed on the TCA baseline configuration, and long-term repeatability, determined by comparing baseline Reference H data from this test to a previous test, have been shown to be good. (2) The longitudinal stability of the TCA is more non-linear than for the Reference H, and while it is similar at normal lift values, the TCA has considerably more pitch-up at higher lift. (3) Longitudinal control effectiveness of the TCA is similar to the Reference H and the ratio of elevator effectiveness to horizontal tail effectiveness is approximately 0.3. (4) The directional stability of the TCA is improved relative to Reference H at higher angles-of attack. The chine is effective for improving directional stability.

Description: The NASA-industry team has sponsored several studies in the last two years to address the installed nozzle boattail drag issues. Some early studies suggested that nozzle boattail drag could be as much as 25 to 40 percent of the subsonic cruise. As part of this study tests have been conducted at NASA-Langley to determine the uninstalled drag characteristics of a proposed nozzle. The overall objective was to determine the effects of nozzle external flap curvature and sidewall boattail variations. This test would also provide data for validating CFD predictions of nozzle boattail drag.

Description: AIRPLANE (Jameson/Baker) is a steady inviscid unstructured Euler flow solver. It has been validated on many HSR geometries. It is implemented as MESHPLANE, an unstructured mesh generator, and FLOPLANE, an iterative flow solver. The surface description from an Intergraph CAD system goes into MESHPLANE as collections of polygonal curves to generate the 3D mesh. The flow solver uses a multistage time stepping scheme with residual averaging to approach steady state, but R is not time accurate. The flow solver was ported from Cray to IBM SP2 by Wu-Sun Cheng (IBM); it could only be run on 4 CPUs at a time because of memory limitations. Meshes for the four cases had about 655,000 points in the flow field, about 3.9 million tetrahedra, about 77,500 points on the surface. The flow solver took about 23 wall seconds per iteration when using 4 CPUs. It took about eight and a half wall hours to run 1,300 iterations at a time (the queue limit is 10 hours). A revised version of FLOPLANE (Thomas) was used on up to 64 CPUs to finish up some calculations at the end. We had to turn on more communication when using more processors to eliminate noise that was contaminating the flow field; this added about 50% to the elapsed wall time per iteration when using 64 CPUs. This study involved computing lift and drag for a wing/body/nacelle configuration at Mach 0.9 and 4 degrees pitch. Four cases were considered, corresponding to four nacelle mass flow conditions.

Description: The primary objectives of this study were to expand the data base showing the effects of LE radius distribution and corresponding . sensitivity to Rn at subsonic and transonic conditions, and to assess the predictive capability of CFD for these effects. Several key elements led to the initiation of this project: 1) the necessity of meeting multipoint design requirements to enable a viable HSCT, 2) the demonstration that blunt supersonic leading-edges can be associated with performance gain at supersonic speeds , and 3) limited data. A test of a modified Reference H model with the TCA planform and 2 LE radius distributions was performed in the NTF, in addition to Navier-Stokes analysis for an additional 3 LE radius distributions. Results indicate that there is a tremendous potential to improve high-lift performance through the use of a blunt LE across the span given an integrated, fully optimized design, and that low Rn data alone is probably not sufficient to demonstrate the benefit.

Description: In cooperation with personnel from the Boeing ANP Laboratory and NASA Langley, a performance test was conducted using the Reference-H 1.675% model ("NASA Modular Model") without nacelles at the NASA Langley 16-Ft Transonic Tunnel. The main objective of the test was to determine the drag reduction achievable with leading-edge and trailing-edge flaps deflected along the outboard wing span at transonic Mach numbers (M = 0.9 to 1.2) for purpose of preliminary design and for comparison with computational predictions. The obtained drag data with flap deflections for Mach numbers of 1.07 to 1.20 are unique for the Reference H wing. Four leading-edge and two trailing-edge flap deflection angles were tested at a mean-wing chord-Reynolds number of about 5.7 million. An outboard-wing leading-edge flap deflection of 81 provides a 4.5 percent drag reduction at M = 1.2 A = 0.2), and much larger values at lower Mach numbers with larger flap deflections. The present results for the baseline (no flaps deflected) compare reasonably well with previous Boeing and NASA Ref-H tunnel tests, including high-Reynolds number NTF results. Viscous CFD simulations using the OVERFLOW thin-layer N.S. method properly predict the observed trend in drag reduction at M = 1.2 as function of leading-edge flap deflection. Modified linear theory properly predicts the flap effects on drag at subsonic conditions (Aero2S code), and properly predicts the absolute drag for the 40 and 80 leading-edge deflection at M = 1.2 (A389 code).

Description: Unstructured grid Euler computations, performed at supersonic cruise speed, are presented for a proposed high speed civil transport configuration, designated as the Technology Concept Airplane (TCA) within the High Speed Research (HSR) Program. The numerical results are obtained for the complete TCA cruise configuration which includes the wing, fuselage, empennage, diverters, and flow through nacelles at Mach 2.4 for a range of angles-of-attack and sideslip. The computed surface and off-surface flow characteristics are analyzed and the pressure coefficient contours on the wing lower surface are shown to correlate reasonably well with the available pressure sensitive paint results, particularly, for the complex shock wave structures around the nacelles. The predicted longitudinal and lateral/directional performance characteristics are shown to correlate very well with the measured data across the examined range of angles-of-attack and sideslip. The results from the present effort have been documented into a NASA Controlled-Distribution report which is being presently reviewed for publication.

Description: The objectives of the Cycle 2 Nonlinear Design Optimization Anlaytical Cross Checks are to: 1) Understand the variability in the predicted performance levels of the nonlinear designs arising from the use of different inviscid (full potential/Euler) and viscous (Navier-Stokes) analysis methods; and 2) Provide the information required to allow the performance levels of all three designs to be validated using the data from the NCV (nonlinear Cruise Validation) model test.

Description: During the last cycle of concept design and wind-tunnel testing, the goal of the low-boom- shaped HSCT concepts (the B-935, the LB-16, and the LB- 1 8) was to meet mission requirements and generate shaped, ground-level pressure signatures with nose shock strengths of 1.0 psf or less. The wind-tunnel tests of these concepts produced results that were partially successful and encouraging although not fully up to expectations. In spite of this, however, these conceptual designs were overly optimistic and not acceptable because: the wing planforms had excessive area; the wing structural aspect ratio was too high; one concept had aft-fuselage rather than under-the-wing engines; and the gross takeoff weights were unrealistically low because of engines that were early, high-tech versions of later, revised, more-realistic engines. The need for reducing the ground-level overpressure shock strengths still existed; a need to be met within more restrictive guidelines of mission performance and gross takeoff weight limitations. Therefore, it was decided that the next conceptual design cycle would focus on decreased nose shock strengths, "boom softening," in the signatures of the Boeing and the McDonnell Douglas baseline concepts rather than low-boom concepts with shaped-signature designs. Overly-optimistic results were not the only problem with these low-sonic-boom concepts. Papers given at the 1994 Sonic-Boom Workshop had demonstrated that the problem of successful nacelle integration on HSCT concepts had only been partially solved. Wind-tunnel pressure signature data, from the HSCT-11B (a.k.a. the LB-18) wind-tunnel model, showed that the Langley HSCT design and analysis method had been successful in reducing the nacelle-volume disturbances in the flow field. This was due.to the engine nacelles mounted behind the wing trailing-edge on the aft fuselage so that no nacelle-wing interference-lift flow-field disturbances were generated. While acceptable from a sonic-boom research point of view, this concept was unacceptable from several practical and structural considerations. Preliminary wind-tunnel pressure signature data from the LB-16 wind-tunnel model, which had the engine nacelles mounted under the wings (the usual location), indicated that the application of the Langley nacelle-integration method had been only partially successful in the reduction of the nacelle-volume with nacelle-wing interference-lift pressure disturbances. So, "boom softening" had to also address the task of successful integration of the engine nacelles, with the engines in the required under-the-wing location. Unless this problem was solved, low-sonic-boom and low-drag modifications to the wing planform, the airfoil shape, and the fuselage longitudinal area distribution could be nullified if the nacelle disturbances added increments to the nose-shock strengths that were removed through component tailoring. In this paper, an arrow-wing boom-softened HSC7 concept which incorporated modifications to a baseline McDonnell Douglas concept is discussed. The analysis of the concept's characteristics will include estimates of weight, center of gravity, takeoff field length, mission range, and predictions of its ground-level sonic-boom pressure signature. Additional modifications which enhanced the softened-boom performance of this concept are also described as well as estimates of the performance penalties induced by these modifications.

Description: A 1:300 scale wind-tunnel model of a conceptual High-Speed Civil Transport (HSCT) designed to generate a shaped, low-boom pressure signature on the ground was tested to obtain sonic-boom pressure signatures in the Langley Research Center Unitary Plan Wind Tunnel at a Mach number of 1.8 and a separation distance of about two body lengths or four wing-spans from the model. Two sets of engine nacelles representing two levels of engine technology were used on the model to determine the effects of increased nacelle volume. Pressure signatures were measured for (model lift)/(design lift) ratios of 0.5, 0.63, 0.75, and 1.0 so that the effect of lift on the pressure signature could be determined. The results of these tests were analyzed and used to discuss the agreement between experimental data and design expectations.

Description: The NASA High Speed Research (HSR) Program is intended to establish a technology base enabling industry development of an economically viable and environmentally acceptable second generation high speed civil transport (HSCT). The objective of the Configuration Aerodynamics task of the program is the development of aerodynamic drag reduction, stability and control, and propulsion airframe integration technologies required to support the HSCT development process. Aerodynamic design tools are being developed, evaluated, and validated through ground based experimental testing. In addition, methods for ground to flight scaling are being developed and refined.

Description: This paper reports on the model, test, and results from the Langley Supersonic Aftbody Closure wind tunnel test. This project is an experimental evaluation of the 1.5% Technology Concept Aircraft (TCA) aftbody closure model (Model 23) in the Langley Unitary Plan Wind Tunnel. The baseline TCA design is the result of a multidisciplinary, multipoint optimization process and was developed using linear design and analysis methods, supplemented with Euler and Navier-Stokes numerical methods. After a thorough design review, it was decided to use an upswept blade attached to the forebody as the mounting system. Structural concerns dictated that a wingtip support system would not be feasible. Only the aftbody part of the model is metric. The metric break was chosen to be at the fuselage station where prior aft-sting supported models had been truncated. Model 23 is thus a modified version of Model 20. The wing strongback, flap parts, and nacelles from Model 20 were used, whereas new aftbodies, a common forebody, and some new tails were fabricated. In summary, significant differences in longitudinal and direction stability and control characteristics between the ABF and ABB aftbody geometries were measured. Correcting the experimental data obtained for the TCA configuration with the flared aftbody to the representative of the baseline TCA closed aftbody will result in a significant reduction in longitudinal stability, a moderate reduction in stabilizer effectiveness and directional stability, and a moderate to significant reduction in rudder effectiveness. These reductions in the stability and control effectiveness levels of the baseline TCA closed aftbody are attributed to the reduction in carry-over area.

Description: There were two objectives for this test. First, was to assess the reasons why there is approximately 1.5 drag counts (cts) discrepancy between measured and computed drag improvement of the Non-linear Cruise Validation (NCV) over the Technology Concept Airplane (TCA) wing body (WB) configurations. The Navier-Stokes (N-S) pre-test predictions from Boeing Commercial Airplane Group (BCAG) show 4.5 drag cts of improvement for NCV over TCA at a lift coefficient (CL) of 0. I at Mach 2.4. The pre-test predictions from Boeing Phantom Works - Long Beach, BPW-LB, show 3.75 drag cts of improvement. BCAG used OVERFLOW and BPW-LB used CFL3D. The first test entry to validate the improvement was held at the NASA Langley Research Center (LARC) UPV;T, test number 1687. The experimental results showed that the drag improvement was only 2.6 cts, not accounting for laminar run and trip drag. This is approximately 1.5 cts less than predicted computationally. In addition to the low Reynolds Number (RN) test, there was a high RN test in the Boeing Supersonic Wind Tunnel (BSWT) of NCV and TCA. BSV@T test 647 showed that the drag improvement of NCV over TCA was also 2.6 cts, but this did account for laminar run and trip drag. Every effort needed to be done to assess if the improvement measured in LaRC UPWT and BSWT was correct. The second objective, once the first objective was met, was to assess the performance increment of NCV over TCA accounting for the associated laminar run and trip drag corrections in LaRC UPWT. We know that the configurations tested have laminar flow on portions of the wing and have trip drag due to the mechanisms used to force the flow to go from laminar to turbulent aft of the transition location.

Description: An improved laminar run and trip drag correction methodology for supersonic cruise performance testing was derived. This method required more careful analysis of the flow visualization images which revealed delayed transition particularly on the inboard upper surface, even for the largest trip disks. In addition, a new code was developed to estimate the laminar run correction. Once the data were corrected for laminar run, the correct approach to the analysis of the trip drag became evident. Although the data originally appeared confusing, the corrected data are consistent with previous results. Furthermore, the modified approach, which was described in this presentation, extends prior historical work by taking into account the delayed transition caused by the blunt leading edges.

Description: The turbulent, incompressible reattaching flow over a rearward-facing step has been studied by many researchers over the years. One of the principal quantities determined in these experiments has been the distance from the step to the point (or region) where the separated shear layer reattaches to the surface (x(r)). The values for x(r)/h, where h is the step height, have covered a wider range than can reasonably be attributed to experimental technique or inaccuracy. Often the reason for a largely different value of x(r)/h can be attributed to an incompletely developed turbulent layer, or a transitional or laminar boundary layer. However, for the majority of experiments where the boundary layer is believed to be fully developed and turbulent, x(r)/h still varies several step heights; generally, 5 1/2 approximately 〈 x(r)/h approximately 〈 7 1/2. This observed variation has usually been attributed to such variables as l/h (step length to height, h/delta (step height to initial boundary-layer thickness), R(e)(theta)), or the experimental technique for determining reattachment location. However, there are so many different combinations of variables in the previous experiments that it was not possible to sort out the effects of particular conditions on the location of reattachment. In the present experiment velocity profiles have been measured in and around the region of separated flow. Results show a large influence of adverse pressure gradient on the reattaching flow over a rearward-facing step that has not been reported previously. Further, the many previous experiments for fully developed, turbulent flow in parallel-walled channels have shown a range of reattachment location that has not been explained by differences in initial flow conditions. Although these initial flow conditions might contribute to the observed variation of reattachment location, it appears that the pressure gradient effect can explain most of that variation.

Description: Gradient-based optimization requires accurate derivatives of the objective function and constraints. These gradients may have previously been obtained by manual differentiation of analysis codes, symbolic manipulators, finite-difference approximations, or existing automatic differentiation (AD) tools such as ADIFOR (Automatic Differentiation in FORTRAN). Each of these methods has certain deficiencies, particularly when applied to complex, coupled analyses with many design variables. Recently, a new AD tool called ADJIFOR (Automatic Adjoint Generation in FORTRAN), based upon ADIFOR, was developed and demonstrated. Whereas ADIFOR implements forward-mode (direct) differentiation throughout an analysis program to obtain exact derivatives via the chain rule of calculus, ADJIFOR implements the reverse-mode counterpart of the chain rule to obtain exact adjoint form derivatives from FORTRAN code. Automatically-generated adjoint versions of the widely-used CFL3D computational fluid dynamics (CFD) code and an algebraic wing grid generation code were obtained with just a few hours processing time using the ADJIFOR tool. The codes were verified for accuracy and were shown to compute the exact gradient of the wing lift-to-drag ratio, with respect to any number of shape parameters, in about the time required for 7 to 20 function evaluations. The codes have now been executed on various computers with typical memory and disk space for problems with up to 129 x 65 x 33 grid points, and for hundreds to thousands of independent variables. These adjoint codes are now used in a gradient-based aerodynamic shape optimization problem for a swept, tapered wing. For each design iteration, the optimization package constructs an approximate, linear optimization problem, based upon the current objective function, constraints, and gradient values. The optimizer subroutines are called within a design loop employing the approximate linear problem until an optimum shape is found, the design loop limit is reached, or no further design improvement is possible due to active design variable bounds and/or constraints. The resulting shape parameters are then used by the grid generation code to define a new wing surface and computational grid. The lift-to-drag ratio and its gradient are computed for the new design by the automatically-generated adjoint codes. Several optimization iterations may be required to find an optimum wing shape. Results from two sample cases will be discussed. The reader should note that this work primarily represents a demonstration of use of automatically- generated adjoint code within an aerodynamic shape optimization. As such, little significance is placed upon the actual optimization results, relative to the method for obtaining the results.

Description: Propulsion performance and operability are key factors in the development of a successful aircraft. For high-speed supersonic aircraft, mixed-compression inlets offer high performance but are susceptible to an instability referred to as unstart. An unstart occurs when a disturbance originating in the atmosphere or the engine causes the shock system to be expelled from the inlet. This event can have adverse effects on control of the aircraft, which is unacceptable for a passenger plane such as the high speed civil transport (HSCT). The ability to predict the transient response of such inlets to flow perturbations is, therefore, important to the proper design of the inlet and the control measures used to prevent unstart. Computational fluid dynamics (CFD) is having an increasing role in the analysis of individual propulsion components. Isolated inlet studies are relatively easy to perform, but a major uncertainty is the boundary condition used at the inlet exit to represent the engine - the so-called compressor face boundary condition. A one-dimensional (1-D) Euler inlet simulation showed that the predicted inlet unstart tolerance to free-stream pressure perturbations can vary by as much as a factor of about six, depending on the boundary condition used. Obviously a thorough understanding of dynamic interactions between inlets and compressors/fans is required to provide the proper boundary condition. To aid in this understanding and to help evaluate possible boundary conditions, an inlet-engine experiment was conducted at the University of Cincinnati. The interaction of acoustic pulses, generated in the inlet, with the engine were investigated. Because of the availability of experimental data for validation, it was decided to simulate the experiment using CFD. The philosophy here is that the inlet-engine system is best simulated by coupling (existing) specialized CFD component-codes. The objectives of this work were to aid in a better understanding of inlet-compressor interaction physics and the formulation of a more realistic compressor-face boundary condition for time-accurate CFD simulations of inlets. Previous simulations have used 1-D Euler engine simulations in conjunction with 1-D Euler and axisymmetric Euler inlet simulations. This effort is a first step toward CFD simulation of an entire engine by coupling multidimensional component codes.

Description: In recent years the techniques of computational fluid dynamics (CFD) have been used to compute flows associated with geometrically complex configurations. However, success in terms of accuracy and reliability has been limited to cases where the effects of turbulence and transition could be modeled in a straightforward manner. Even in simple flows, the accurate computation of skin friction and heat transfer using existing turbulence models has proved to be a difficult task, one that has required extensive fine-tuning of the turbulence models used. In more complex flows (for example, in turbomachinery flows in which vortices and wakes impinge on airfoil surfaces causing periodic transitions from laminar to turbulent flow) the development of a model that accounts for all scales of turbulence and predicts the onset of transition is an extremely difficult task. Fortunately, current trends in computing suggest that it may be possible to perform direct simulations of turbulence and transition at moderate Reynolds numbers in some complex cases in the near future. This presentation will focus on direct simulations of transition and turbulence using high-order accurate finite-difference methods. The advantage of the finite-difference approach over spectral methods is that complex geometries can be treated in a straightforward manner. Additionally, finite-difference techniques are the prevailing methods in existing application codes. An application of high-order-accurate finite-difference methods to direct simulations of transition and turbulence in a spatially evolving boundary layer subjected to high levels of freestream turbulence will be presented.

Description: The SPAcecraft SIMulator (SPASIM) simulates the functions and resources of a spacecraft to quickly perform conceptual design (Phase A) trade-off and sensitivity analyses and uncover any operational bottlenecks during any part of the mission. Failure modes and operational contingencies can be evaluated allowing operational planning (what-if scenarios) and optimization for a range of mission scenarios. The payloads and subsystems are simulated, using a hierarchy of graphical models, in terms of how their functions affect resources such as propellant, power, and data. Any of the inputs and outputs of the payloads and subsystems can be plotted during the simulation or stored in a file so they can be used by other programs. Most trade-off analyses, including those that compare current versus advanced technology, can be performed by changing values in the parameter menus. However, when a component is replaced by one with a different functional architecture, its graphical model can also be modified or replaced by drawing from a component library. SPASIM has been validated using several spacecraft designs that were at least at the Critical Design Review level. The user and programmer guide, including figures, is available on line as a hypertext document. This is an easy-to-use and expandable tool which is based on MATLAB(R) and SIMULINK(R). It runs on Silicon Graphics Inc. workstations and personal computers with Windows 95(TM) or NT(TM).

Description: Turbulent flows that exhibit inhomogeneities in the streamwise direction pose a particular challenge to numerical simulation approaches due to the need to prescribe time-dependent turbulent inflow conditions. In most cases the flow downstream is more or less 'driven' by the conditions at the inlet, making it necessary to specify realistic turbulent fluctuations that are in equilibrium with the assumed mean flow. This requirement often dictates that the inflow data should satisfy the Navier-Stokes equations, which in turn implies that an independent simulation be used to generate them. Detailed simulations for the purpose of creating inflow conditions can be costly and thus certain levels of approximation are desirable. In this paper we shall focus on an approximate yet accurate method for generating inflow conditions for spatially-developing boundary layer simulations. The proposed method is essentially a simplification of the method of Spalart and Leonard (1985), who devised an ingenious transformation that allows for the calculation of spatially evolving boundary layers in conjunction with periodic boundary conditions applied in the streamwise direction. While this method is elegant and highly accurate, it is more complicated than is necessary for the purpose of generating inflow data. A few key approximations are used in this work to arrive at a 'modified Spalart method' that is very easy to implement and efficient to use. The new method is shown to yield results that compare well with the computations of Spalart (1988). When used as a means of generating inflow data, the modified Spalart method is shown to be superior to existing approaches.

Description: Accurate large eddy simulation (LES) of a wall-bounded flow generally requires a near-wall resolution comparable to that in direct numerical simulation (DNS). As much as 50% of the total grid points and computational costs are expended in the near-wall regions in a typical simulation. This limits LES to fairly low Reynolds numbers on current computers. To perform practical flow applications at realistically high Reynolds numbers, such as flow over an airfoil, it is desirable to replace very thin, near-wall regions in the LES with easily and inexpensively computed wall models to specify the near-wall boundary conditions.

Description: Reducing the infinite degrees of freedom inherent in fluid motion into a manageable number of modes to analyze fluid motion is presented. The concepts behind the center manifold technique are used. Study of the Blasius boundary layer and a precise description of stability within the flow field are discussed.

Description: A key ingredient in the simulation of flow in porous media is the accurate determination of the velocities that drive the flow. The large scale irregularities of the geology, such as faults, fractures, and layers suggest the use of irregular grids in the simulation. Work has been done in applying the finite volume element (FVE) methodology as developed by McCormick in conjunction with mixed methods which were developed by Raviart and Thomas. The resulting mixed finite volume element discretization scheme has the potential to generate more accurate solutions than standard approaches. The focus of this paper is on a multilevel algorithm for solving the discrete mixed FVE equations. The algorithm uses a standard cell centered finite difference scheme as the 'coarse' level and the more accurate mixed FVE scheme as the 'fine' level. The algorithm appears to have potential as a fast solver for large size simulations of flow in porous media.

Description: The effect of the unsteady flow field in a multistage compressor on the time-averaged performance was assessed. The energy transport by the unsteady deterministic flow field was taken into account. The magnitude of the body force resulting from the aerodynamic loading on a blade row was compared to the gradient of the stress tensor associated with the unsteady time-resolved flow field generated by the blade row. The magnitude of the work performed by these forces was compared to the divergence of the energy correlations produced by the unsteady time-resolved flow field. The stress tensor and the energy correlations are non-negligible in the end wall regions. The results suggest that the turbulence is the primary source of flow mixing in the midspan region.

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: Computations have been performed on the baseline Reference H wing/body configuration, as well as the Wing 704 configuration, an optimized wing and fuselage combination derived from Ref. H through automated optimization. The parabolized Navier-Stokes solver UPS was employed with viscous terms in two directions in an effort to understand the source and level of potential viscous/inviscid interactions. The paper briefly describes the UPS code and the grids used to obtain the solutions before the discussion of results. Results of these computations indicate that viscous/inviscid interaction can contribute increments to both the pressure- and friction-related drag. Computations were performed for wind tunnel conditions-1.675% scale models at a Reynolds number of 4 million per foot. Turbulent flow results were obtained using the Baldwin-Lomax algebraic turbulence model and were compared with laminar flow results. The laminar flow fields were used to obtain upper bounds on potential interaction effects.

Description: The NASA-industry team has sponsored several studies in the last two years to address the installed nozzle boattail drag issues. Some early studies suggested that nozzle boattail drag could be as much as 25 to 40 percent of the subsonic cruise. As part of this study tests have been conducted at NASA-Langley to determine the uninstalled drag characteristics of a proposed nozzle. The overall objective was to determine the effects of nozzle external flap curvature and sidewall boattail variations. This test would also provide data for validating CFD predictions of nozzle boattail drag.

Description: The objective of this milestone is to assess the propulsion/airframe integration characteristics of the Technology Concept Airplane and design variations through computational analysis and experimental subsonic through supersonic wind tunnel testing. The Milestone will generate a comprehensive CFD and wind tunnel data base of the baseline, and design variations. Emphasis will be placed on establishing the propulsion induced effects on the flight performance of the Technology Concept Airplane with all appropriate wind tunnel corrections.