ALBERT

Description: The aerodynamic development of an engine inlet requires a comprehensive program of both wind tunnel testing and Computational Fluid Dynamics (CFD) simulations. To save time and resources, much "testing" is done using CFD before any design ever enters a wind tunnel. The focus of my project this summer is on CFD analysis tool development. In particular, I am working to further develop the capabilities of the Planar Inlet Design and Analysis Process (PINDAP). "PINDAP" is a collection of computational tools that allow for efficient and accurate design and analysis of the aerodynamics about and through inlets that can make use of a planar (two-dimensional or axisymmetric) geometric and flow assumption. PINDAP utilizes the WIND CFD flow solver, which is capable of simulating the turbulent, compressible flow field. My project this summer is a continuation of work that I performed for two previous summers. Two years ago, I used basic features of the PINDAP to design a Mach 5 hypersonic scramjet engine inlet and to demonstrate the feasibility of the PINDAP. The following summer, I worked to develop its geometry and grid generation capabilities to include subsonic and supersonic inlets, complete bodies and cowls, conic leading and trailing edges, as well as airfoils. These additions allowed for much more design flexibility when using the program.

Description: We have designed and built an instrument to measure and monitor the "nightglow" of the Earth's atmosphere in the near ultraviolet (NUV). In this paper we describe the design of this instrument, called NIGHTGLOW. NIGHTGLOW is designed to be flown-from a high altitude research balloon, and circumnavigate the globe. NIGHTGLOW is a NASA, University of Utah, and New Mexico State University project. A test flight took place from Palestine, Texas on July 5, 2000, lasting about 8 hours. The instrument performed well and landed safely in Stiles, Texas with little damage. The resulting measurements of the NUV nightglow are consistent with previous measurements from sounding rockets and balloons. The results will be presented and discussed.

Description: A brief review of the evolutionary progress in computational aerothermodynamics is presented. The current status of computational aerothermodynamics is then discussed, with emphasis on its capabilities and limitations for contributions to the design process of hypersonic vehicles. Some topics to be highlighted include: (1) aerodynamic coefficient predictions with emphasis on high temperature gas effects; (2) surface heating and temperature predictions for thermal protection system (TPS) design in a high temperature, thermochemical nonequilibrium environment; (3) methods for extracting and extending computational fluid dynamic (CFD) solutions for efficient utilization by all members of a multidisciplinary design team; (4) physical models; (5) validation process and error estimation; and (6) gridding and solution generation strategies. Recent experiences in the design of X-33 will be featured. Computational aerothermodynamic contributions to Mars Pathfinder, METEOR, and Stardust (Comet Sample return) will also provide context for this discussion. Some of the barriers that currently limit computational aerothermodynamics to a predominantly reactive mode in the design process will also be discussed, with the goal of providing focus for future research.

Description: Flow separation in the low-pressure turbine (LPT) is a major barrier that limits further improvements of aerodynamic designs of turbine airfoils. The separation is responsible for performance degradation, and it prevents the design of highly loaded airfoils. The separation can be delayed, reduced, or eliminated completely if flow control techniques are used. Successful flow control technology will enable breakthrough improvements in gas turbine performance and design. The focus of this research project was the development and experimental demonstration of active separation control using glow discharge plasma (GDP) actuators in flow conditions simulating the LPT. The separation delay was shown to be successful, laying the foundation for further development of the technologies to practical application in the LPT. In a fluid mechanics context, the term "flow control" means a technology by which a very small input results in a very large effect on the flow. In this project, the interest is to eliminate or delay flow separation on LPT airfoils by using an active flow control approach, in which disturbances are dynamically inserted into the flow, they interact with the flow, and they delay separation. The disturbances can be inserted using a localized, externally powered, actuating device, examples are acoustic, pneumatic, or mechanical devices that generate vibrations, flow oscillations, or pulses. A variety of flow control devices have been demonstrated in recent years in the context of the external aerodynamics of aircraft wings and airframes, where the incoming flow is quiescent or of a very low turbulence level. However, the flow conditions in the LPT are significantly different because there are high levels of disturbances in the incoming flow that are characterized by high free-stream turbulence intensity. In addition, the Reynolds number, which characterizes the viscous forces in the flow and is related to the flow speed, is very low in the LPT passages.

Description: Researchers are constantly pursuing technologies that will increase the performance of gas turbine engines. The aspirated compressor concept discussed here would allow the compression system to perform its task with about one-half of the compressor blades. To accomplish this, the researchers applied boundary layer control to the blades, casing, and hub. This method of boundary layer control consisted of removing small amounts of air from the main flow path at critical areas of the compressor. This bleed air could be used by other systems such as engine cooling or could be re-injected into lower pressure areas that require air for enhanced performance. This effort was initiated by the Massachusetts Institute of Technology (MIT) in response to a solicitation from the Defense Advanced Research Projects Agency (DARPA) who sought to advance research in flow control technology. The NASA Glenn Research Center partnered with MIT (principal investigator), Honeywell Aircraft Engines (cycle analysis, structural analysis, and mechanical design), and Pratt & Whitney (cycle analysis and aero-analysis) to conceptualize, design, analyze, build, and test the aspirated fan stage. The aero-design and aero-analysis of this fan stage were jointly executed by MIT and Glenn to minimize the amount of bleed flow needed and to maintain the highest efficiency possible (ref. 1). Mechanical design issues were complicated by the need to have a shrouded rotor with hollow blades, with rotor stress levels beyond the capabilities of titanium. The high stress issues were addressed by designing a shroud that was filament wound with a carbon fiber/epoxy matrix, resulting in an assembly that was strong enough to handle the high stresses. Both the rotor (preceding photographs) and stator (following photograph) were fabricated in two halves and then bolted together at the hub and tip, permitting the bleed passages to be machined into each half before assembly.

Description: The benchmark solution for the cascade-gust interaction problem is computed using a linearized Euler code called LINFLUX. The inherently three-dimensional code is run in the thin-annulus limit to compute the two-dimensional cascade response. The calculations are carried out in the frequency-domain and the unsteady response at each of the gust s three frequency component is computed. The results are presented on modal basis for pressure perturbations (i.e., acoustic modes) as well as velocity perturbations (i.e., convected gust modes) at each frequency.

Description: An unsteady aerodynamic code, called GUST3D (ref. 3), has been developed to solve equation (8) for flows with periodic vortical disturbances. The code uses a frequency-domain approach with second-order central differences and a pressure radiation condition in the far field. GUST3D requires as input certain mean flow quantities which are calculated separately by a potential flow solver. The solver calculates the mean ow using a Gothert's Rule approximation (ref. 3). On the airfoil surface, it uses the solution calculated by the potential code FLO36 (ref. 4). Figures 1-2 show the mean pressure along the airfoil surface for the two airfoil geometries. In Figures 3 - 8, we present the RMS pressure on the airfoil surface. Each figure shows three GUST3D solutions (calculated on grids with different far-field boundary locations). Three solutions are shown to provide some indication of the numerical uncertainty in the results. Figures 9 - 13 present the acoustic intensity. We again show three solutions per case. Note that no results are presented for the k1 = k2 = 2.0 loaded airfoil case, as an acceptable solution could not be obtained. A few comments need to be made about the results shown. First, since the last Workshop, the GUST3D code has been substantially upgraded. This includes implementing a more accurate far-field boundary condition (ref. 5) and developing improved gridding capabilities. This is the reason for any differences that may exist between the present results and results from the last Workshop. Second, the intensity results on the circle R = 4C were obtained using a Kirchoff method (ref. 6). The Kirchoff surface was the circle R = 2C. Finally, the GUST3D code is most accurate for low reduced frequencies. A new domain decomposition approach (ref. 7) has been developed to improve accuracy. Both the single domain and domain decomposition approaches were used in generating the present results.

Description: The TURBO-AE aeroelastic code has been used to verify the flutter stability of the Efficient Low Noise Fan (ELNF), which is also referred to as the trailing-edge blowing fan. The ELNF is a unique technology demonstrator being designed and fabricated at the NASA Glenn Research Center for testing in Glenn's 9-by-15-Foot Low-Speed Wind Tunnel. In the ELNF, air can be blown out of slots near the trailing edges of the fan blades to fill in the wakes downstream of the rotating blades. This filling of the wakes leads to a reduction of the rotor-stator interaction (tone) noise that results from the interaction of wakes with the downstream stators. The ELNF will demonstrate a 1.6-EPNdB1 reduction in tone noise through wake filling, without increasing the broadband noise. Furthermore, the reduced blade row interaction will decrease the possibility of forced response and enable closer spacing of blade rows, thus reducing engine length and weight. During the design of the ELNF, the rotor blades were checked for flutter stability using the detailed aeroelastic analysis capability of the three-dimensional Navier-Stokes TURBOAE code. The aeroelastic calculations were preceded by steady calculations in which the blades were not allowed to vibrate. For each rotational speed, as the back-pressure was increased, the mass flow rate decreased, and the operating point moved along the constant speed characteristic (speed-line) from choke to stall as shown on the fan map. The TURBO-AE aeroelastic analyses were performed separately for the first two vibration modes (bending and torsion) and covered the complete range of interblade phase angles or nodal diameters at which flutter can occur. The results indicated that the ELNF blades would not encounter flutter at takeoff conditions. The calculations were then repeated for a part-speed condition (70-percent rotational speed), and the results again showed no flutter in the operating region. On the fan map (shown), the predicted flutter point at part speed condition was located beyond the stall line, which means that the ELNF will not encounter flutter since it will never operate beyond the stall line. All the calculations done so far have been for the nonblowing case, and selected calculations will be repeated with air blowing from the trailing edge of the fan.

Description: Flow structure of momentum-dominated pure helium jets discharged vertically into ambient air was investigated using high-speed rainbow schlieren deflectometry (RSD) technique. Effects of the operating parameters, i.e., Reynolds number (Re) and Richardson number (Ri), on the oscillatory behavior of the flow were examined over a range of experimental conditions. To seek the individual effect of these parameters, one of them was fixed and the other was varied with certain constraints. Measurements revealed highly periodic oscillations in the laminar region as well as high regularity in transition and turbulent regions. Maximum spectral power profiles at different axial locations indicated the oscillation amplitude increasing until the breakdown of the jet in the turbulent regime. The transition from the laminar to turbulent flow was also investigated. Fast Fourier transform analysis performed in the transition regime showed that the flow oscillates at a unique frequency, which was the same in the upstream laminar flow region. Measured deflection angle data were used in Abel inversion algorithm to construct the helium concentration fields. Instantaneous helium concentration contours revealed changes in the flow structure and evolution of vortical structures during an oscillation cycle. Temporal evolution plots of helium concentration at different axial location showed repeatable oscillations at all axial and radial locations up to the turbulent regime. A cross-correlation technique, applied to find the spatial displacements of the vortical structures, provided correlation coefficient peaks between consecutive schlieren images. Results show that the vortical structure convected and accelerated only in the axial direction.

Description: The NASA-Ames ballistic range provides a unique capability for aerothermodynamic testing of configurations in hypersonic, real-gas, free-flight environments. The facility can closely simulate conditions at any point along practically any trajectory of interest experienced by a spacecraft entering an atmosphere. Sub-scale models of blunt atmospheric entry vehicles are accelerated by a two-stage light-gas gun to speeds as high as 20 times the speed of sound to fly ballistic trajectories through an 24 m long vacuum-rated test section. The test-section pressure (effective altitude), the launch velocity of the model (flight Mach number), and the test-section working gas (planetary atmosphere) are independently variable. The model travels at hypersonic speeds through a quiescent test gas, creating a strong bow-shock wave and real-gas effects that closely match conditions achieved during actual atmospheric entry. The challenge with ballistic range experiments is to obtain quantitative surface measurements from a model traveling at hypersonic speeds. The models are relatively small (less than 3.8 cm in diameter), which limits the spatial resolution possible with surface mounted sensors. Furthermore, since the model is in flight, surface-mounted sensors require some form of on-board telemetry, which must survive the massive acceleration loads experienced during launch (up to 500,000 gravities). Finally, the model and any on-board instrumentation will be destroyed at the terminal wall of the range. For these reasons, optical measurement techniques are the most practical means of acquiring data. High-speed thermal imaging has been employed in the Ames ballistic range to measure global surface temperature distributions and to visualize the onset of transition to turbulent-flow on the forward regions of hypersonic blunt bodies. Both visible wavelength and infrared high-speed cameras are in use. The visible wavelength cameras are intensified CCD imagers capable of integration times as short as 2 ns. The infrared camera uses an Indium Antimonide (InSb) sensor in the 3 to 5 micron band and is capable of integration times as short as 500 ns. The projectiles are imaged nearly head-on using expendable mirrors offset slightly from the flight path. The proposed paper will discuss the application of high-speed digital imaging systems in the NASA-Ames hypersonic ballistic range, and the challenges encountered when applying these systems. Example images of the thermal radiation from the blunt nose of projectiles flying at nearly 14 times the speed of sound will be given.

Description: The Flight Dynamics Analysis Branch (FDAB) of NASNGSFC is tasked with determining the orbital lifetime of several developmental and operational satellites, which include the Hubble Space Telescope. A DOS based program developed by the FDAB many years ago, called PC Lifetime, is used to determine a satellite s lifetime and could soon be in need of a replacement. STK s Lifetime Object Tool is a possible candidate. Due to the reduced support of the PC Lifetime program, and the growing incompatibility of older programs with new operating systems, a comparative analysis was done to determine if STWLifetime could meet the stringent requirements that were laid before it. The use of highly accurate numerical propagators such as STK s High Precision Orbit Propagator ( OP) and the Goddard Trajectory Determination System (GTDS) provided a basis on which to compare STWLifetime s results. Several test cases were run, but the main four test cases would determine whether or not STWLifetime could be PC- Lifetime s replacement. These four cases include a geotransfer orbit, two circular LEOS, and a Poiar LEO. Following rigorous testmg procedures, a conclusion will be determined. STK has proved to be a versatile program on many satellite missions and the FDAB has high hopes that it can pass FDAB s requirements for orbital lifetime prediction.

Description: Effective delay of turbulent boundary layer separation could be achieved via closed-loop control. Constructing such a system requires that sensor data be processed, real-time, and fed into the controller to determine the output. Current methods for detection of turbulent boundary layer separation are lacking the capability of localized, fast and reliable identification of the boundary layer state. A method is proposed for short-time FFT processing of time series, measured by hot-film sensors, with the purpose of identifying the alternation of the balance between small and large scales as the boundary layer separates, favoring the large scales. The method has been validated by comparison to other criteria of separation detection and over a range of baseline and controlled flow conditions on a simplified high-lift system, incorporating active flow control.

Description: Results of an experimental investigation on isolated synthetic jets in crossflow from rectangular orifices of different aspect ratio and orientation are presented. Three aspect ratios, AR = 4, 8, and 16, with pitch a = 90deg., were investigated, Additionally, the AR = 8 case was pitched at 20 deg.. The yaw angle, , was varied through 0 deg., 10 deg., 45 deg. and 90 deg.. All orifices had same exit area and the data were compared with synthetic as well as steady jet from a circular orifice of same area. Hotwire measurements were performed to obtain all three components of mean velocity and turbulent stresses. Data were acquired for momentum-flux ratio up to J = 50. Distributions of time- and phase-averaged data were obtained on the cross sectional plane at x/D = 0.5, 5 and 10, as well as on the axial plane of the symmetry. Qualitative flowfield similarity between synthetic and steady jets is observed. However, high-momentum 'cap' above the low-momentum 'dome', characteristic of steady jets, does not necessarily appear in the synthetic jet. The position and shape of the high-momentum region depend on the distance from the orifice, pitch, yaw as well as momentum-flux ratio. Consequently, the location of the minimum velocity in the 'dome' measured at the plane of symmetry, y(submin), is adopted as a reference for penetration estimate and trajectory comparison. For AR = 16, the dome is the largest in area with maximum velocity deficit. However, the penetration is somewhat higher for AR = 4. Increase in yaw reduces the spatial extent of the dome and the penetration height but augments the velocity deficit. At low J the dome is connected to the boundary layer and traces of the cap of high momentum fluid are visible above it. Increase in J lifts the dome and reorganizes the high-momentum fluid around its perimeter, eventually bringing it underneath. Phase-averaged data document dynamic topological changes within the cycle. Phase-averaged streamwise velocity contours on the cross-sectional plane exhibit behavior commensurate with that seen in time-averaged data at various J.

Description: This paper presents the results of a computational fluid dynamics validation exercise using the General Aerodynamic Simulation Program (GASP) code for vehicles under consideration for the X-38 program. The scope of the exercise involved simulation and comparison with wind tunnel aeroheating and aerodynamic data in the hypersonic regime. In addition, simulations were performed and results are presented for the full-scale flight vehicle along an expected re- entry trajectory.

Description: A NASA team of engineers has been organized to design a crew return vehicle for returning International Space Station crew members from orbit. The hypersonic aerothermodynamic characteristics of the X-23/X-24A derived X-38 crew return vehicle are being evaluated in various wind tunnels in support of this effort. Aerothermodynamic data from two NASA hypersonic tunnels at Mach 6 and Mach 10 has been obtained with cast ceramic models and a thermographic phosphorus digital imaging system. General windward surface heating features are described based on experimental surface heating images and surface oil flow patterns for the nominal hypersonic aerodynamic orientation. Body flap reattachment heating levels are examined. Computational Fluid Dynamics tools have been applied at the appropriate wind tunnel conditions to make comparisons with this data.

Description: This final report summarizes the research accomplished under NASA-Ames Cooperative Agreement NCC2-5517 during the period which extended from August 1, 2002 to July 31,2004.

Description: A variety of novel control effector concepts have recently emerged that may enable new approaches to flight control. In particular, the potential exists to shift the composition of the typical aircraft control effector suite from a small number of high authority, specialized devices (rudder, aileron, elevator, flaps), toward larger numbers of smaller, less specialized, distributed device arrays. The concept envisions effector and sensor networks composed of relatively small high-bandwidth devices able to simultaneously perform a variety of control functions using feedback from disparate data sources. To investigate this concept, a remotely piloted flight vehicle has been equipped with an array of 24 trailing edge shape-change effectors and associated pressure measurements. The vehicle, called the Multifunctional Effector and Sensor Array (MESA) testbed, was recently tested in NASA Langley's 12-ft Low Speed wind tunnel to characterize its stability properties, control authorities, and distributed pressure sensitivities for use in a dynamic simulation prior to flight testing. Another objective was to implement and evaluate a scheme for actively controlling the spanwise pressure distribution using the shape-change array. This report describes the MESA testbed, design of the pressure distribution controller, and results of the wind tunnel test.

Description: Several flight-critical aerodynamic problems such as buffet, flutter, stall, and wing rock are strongly affected or caused by abrupt changes in unsteady aerodynamic loads and moments. Advanced sensing and flow diagnostic techniques have made possible simultaneous identification and tracking, in realtime, of the critical surface, viscosity-related aerodynamic phenomena under both steady and unsteady flight conditions. The wind tunnel study reported here correlates surface hot-film measurements of leading edge stagnation point and separation point, with unsteady aerodynamic loads on a NACA 0015 airfoil. Lift predicted from the correlation model matches lift obtained from pressure sensors for an airfoil undergoing harmonic pitchup and pitchdown motions. An analytical model was developed that demonstrates expected stall trends for pitchup and pitchdown motions. This report demonstrates an ability to obtain unsteady aerodynamic loads in real time, which could lead to advances in air vehicle safety, performance, ride-quality, control, and health management.

Description: The motivation for Computational Aeroelasticity (CA) and the elements of one type of the analysis or simulation process are briefly reviewed. The need for streamlining and improving the overall process to reduce elapsed time and improve overall accuracy is discussed. Further effort is needed to establish the credibility of the methodology, obtain experience, and to incorporate the experience base to simplify the method for future use. Experience with the application of a variety of Computational Aeroelasticity programs is summarized for the transonic flutter of two wings, the AGARD 445.6 wing and a typical business jet wing. There is a compelling need for a broad range of additional flutter test cases for further comparisons. Some existing data sets that may offer CA challenges are presented.

Description: This paper presents novel analytical results for eigenvalues and eigenvectors produced using discrete time aerodynamic and aeroelastic models. An unsteady, incompressible vortex lattice aerodynamic model is formulated in discrete time; the importance of several modeling parameters is examined. A detailed study is made of the behavior of the aerodynamic eigenvalues both in discrete and continuous time. The aerodynamic model is then incorporated into aeroelastic equations of motion. Eigenanalyses of the coupled equations produce stability results and modal characteristics which are valid for critical and non-critical velocities. Insight into the modeling and physics associated with aeroelastic system behavior is gained by examining both the eigenvalues and the eigenvectors. Potential pitfalls in discrete time model construction and analysis are examined.

Description: A Momentum Source Model has been revised and implemented in the flow solver OVERFLOW-D. In this approach, the fan forces are evaluated from two-dimensional airfoil tables as a function of local Mach number and angle-of-attack and applied as source terms in the discretized Navier-Stokes equations. The model revisions include a new model for forces in the tip region and axial distribution of the source terms. The model revisions improve the results significantly. The Momentum Source Model agrees well with a discrete blade model for all computed collective pitch angles. The two models agree well with experimental data for thrust vs. torque. The Momentum Source Model is a good complement to Discrete Blade Models for ducted fan computations. The lower computational and labor costs make parametric studies, optimization studies and interactional aerodynamics studies feasible for cases beyond what is practical with a Discrete Blade Model today.

Description: The effect of nozzle trailing edge thickness on broadband acoustic radiation and the production of tones is investigated for coannular nozzles. Experiments were performed for a core nozzle trailing edge thickness between 0.38 mm and 3.17 mm. The on-set of discrete tones was found to be predominantly affected by the velocity ratio, the ratio of the fan velocity to the core velocity, although some dependency on trailing edge thickness was also noted. For a core nozzle trailing edge thickness greater than or equal to 0.89 mm, tones were produced for velocity ratios between 0.91 and 1.61. For a constant nozzle trailing edge thickness, the frequency varied almost linearly with the core velocity. The Strouhal number based on the core velocity changed with nozzle trailing edge thickness and varied between 0.16 and 0.2 for the core nozzles used in the experiments. Increases in broadband noise with increasing trailing edge thickness were observed for tone producing and non-tone producing conditions. A variable thickness trailing edge (crenellated) nozzle resulted in no tonal production and a reduction of the broadband trailing edge noise relative to that of the corresponding constant thickness trailing edge.

Description: Aerodynamic tests in support of the Columbia accident investigation were conducted in two hypersonic wind tunnels at the NASA Langley Research Center, the 20-Inch Mach 6 Air Tunnel and the 20-Inch CF4 Tunnel. The primary purpose of these tests was to measure the forces and moments generated by a variety of outer mold line alterations (damage scenarios) using 0.0075-scale models of the Space Shuttle Orbiter. Simultaneously acquired global heat transfer mappings were obtained for a majority of the configurations tested. Test parametrics included angles of attack from 38 to 42 deg, unit Reynolds numbers from 0.3 x 10(exp 6) to 3.0 x 10(exp 6) per foot, and normal shock density ratios of 5 (Mach 6 air) and 12 (CF4). The damage scenarios evaluated included asymmetric boundary layer transition, gouges in the windward surface thermal protection system tiles, wing leading edge damage (partially and fully missing reinforced carbon-carbon (RCC) panels), deformation of the wing windward surface, and main landing gear and/or door deployment. The measured aerodynamic increments for the damage scenarios examined were generally small in magnitude, as were the flight-derived values during most of the entry prior to loss of communication. A progressive damage scenario is presented that qualitatively matches the flight observations for the STS-107 entry.

Description: The main objective of this work is to predict the Aerothermal Performance Constraint (APC) for a class of Crew Transfer Vehicles (CTV) with shap noses and wing leading edges made out of UHTC which is a family of Ultra High Temperature Ceramics materials developed at NASA Ames. The APC is based on the theoretical temperature limit of the material which is usually encountered at the CTV nose or wing leading edge. The APC places a lower limit on the trajectory of the CTV in the altitude velocity space. The APC is used as one of the constraints in developing reentry and abort trajectories for the CTV. The trajectories are then used to generate transient thermal response of the nosecaps and wing leading edges which are represented as either a one piece of UHTC or two piece (UHTC + RCC) with perfect axial contact. The final paper will include more details about the analysis procedure and will also include results for reentry and abort design trajectories.

Description: A collection of computer programs has been developed that implements a procedure for predicting multiple-pure-tone (MPT) noise generated by fan blades of an aircraft engine (e.g., a turbofan engine). MPT noise arises when the fan is operating with supersonic relative tip Mach No. Under this flow condition, there is a strong upstream running shock. The strength and position of this shock are very sensitive to blade geometry variations. For a fan where all the blades are identical, the primary tone observed upstream of the fan will be the blade passing frequency. If there are small variations in geometry between blades, then tones below the blade passing frequency arise MPTs. Stagger angle differences as small as 0.1 can give rise to significant MPT. It is also noted that MPT noise is more pronounced when the fan is operating in an unstarted mode. Computational results using a three-dimensional flow solver to compute the complete annulus flow with non-uniform fans indicate that MPT noise can be estimated in a relatively simple way. Hence, once the effect of a typical geometry variation of one blade in an otherwise uniform blade row is known, the effect of all the blades being different can be quickly computed via superposition. Two computer programs that were developed as part of this work are used in conjunction with a user s computational fluid dynamics (CFD) code to predict MPT spectra for a fan with a specified set of geometric variations: (1) The first program ROTBLD reads the users CFD solution files for a single blade passage via an API (Application Program Interface). There are options to replicate and perturb the geometry with typical variations stagger, camber, thickness, and pitch. The multi-passage CFD solution files are then written in the user s file format using the API. (2) The second program SUPERPOSE requires two input files: the first is the circumferential upstream pressure distribution extracted from the CFD solution on the multi-passage mesh, the second file defines the geometry variations of each blade in a complete fan. Superposition is used to predict the spectra resulting from the geometric variations.

Description: Boundary Layer separation causes the airfoil to stall and therefore imposes dramatic performance degradation on the airfoil. In recent years, flow separation control has been one of the active research areas in the field of aerodynamics due to its promising performance improvements on the lifting device. These active flow separation control techniques include steady and unsteady air injection as well as suction on the airfoil surface etc. This paper will be focusing on the steady and unsteady air injection on the airfoil. Although wind tunnel experiments revealed that the performance improvements on the airfoil using injection techniques, the details of how the key variables such as air injection slot geometry and air injection angle etc impact the effectiveness of flow separation control via air injection has not been studied. A parametric study of both steady and unsteady air injection active flow control will be the main objective for this summer. For steady injection, the key variables include the slot geometry, orientation, spacing, air injection velocity as well as the injection angle. For unsteady injection, the injection frequency will also be investigated. Key metrics such as lift coefficient, drag coefficient, total pressure loss and total injection mass will be used to measure the effectiveness of the control technique. A design of experiments using the Box-Behnken Design is set up in order to determine how each of the variables affects each of the key metrics. Design of experiment is used so that the number of experimental runs will be at minimum and still be able to predict which variables are the key contributors to the responses. The experiments will then be conducted in the 1ft by 1ft wind tunnel according to the design of experiment settings. The data obtained from the experiments will be imported into JMP, statistical software, to generate sets of response surface equations which represent the statistical empirical model for each of the metrics as a function of the key variables. Next, the variables such as the slot geometry can be optimized using the build-in optimizer within JMP. Finally, a wind tunnel testing will be conducted using the optimized slot geometry and other key variables to verify the empirical statistical model. The long term goal for this effort is to assess the impacts of active flow control using air injection at system level as one of the task plan included in the NASAs URETI program with Georgia Institute of Technology.

Description: A performance evaluation of the Hybrid Predictor corrector Aerocapture Scheme (HYPAS) guidance algorithm for aerocapture at Neptune is presented in this paper for a Mission to Neptune and the Neptune moon Triton'. This mission has several challenges not experienced in previous aerocapture guidance assessments. These challengers are a very high Neptune arrival speed, atmospheric exit into a high energy orbit about Neptune, and a very high ballistic coefficient that results in a low altitude acceleration capability when combined with the aeroshell LD. The evaluation includes a definition of the entry corridor, a comparison to the theoretical optimum performance, and guidance responses to variations in atmospheric density, aerodynamic coefficients and flight path angle for various vehicle configurations (ballistic numbers). The benefits of utilizing angle-of-attack modulation in addition to bank angle modulation to improve flight performance is also discussed. The results show that despite large sensitivities in apoapsis targeting, the algorithm performs within the allocated AV budget for the Neptune mission bank angle only modulation. The addition of angle-of-attack modulation with as little as 5 degrees of amplitude significantly improves the scatter in final orbit apoapsis. Although the angle-of-attack modulation complicates the vehicle design, the performance enhancement reduces aerocapture risk and reduces the propellant consumption needed to reach the high energy target orbit for a conventional propulsion system.

Description: NASA's Strategic Plan for the Aerospace Technology Enterprise includes ambitious objectives focused on affordable air travel, reduced emissions, and expanded aviation-system capacity. NASA Dryden Flight Research Center, in cooperation with NASA Ames Research Center, the Boeing Company, and the University of California, Los Angeles, has embarked on an autonomous-formation-flight project that promises to make significant strides towards these goals. For millions of years, birds have taken advantage of the aerodynamic benefit of flying in formation. The traditional "V" formation flown by many species of birds (including gulls, pelicans, and geese) enables each of the trailing birds to fly in the upwash flow field that exists just outboard of the bird immediately ahead in the formation. The result for each trailing bird is a decrease in induced drag and thus a reduction in the energy needed to maintain a given speed. Hence, for migratory birds, formation flight extends the range of the system of birds over the range of birds flying solo. The Autonomous Formation Flight (AFF) Project is seeking to extend this symbiotic relationship to aircraft.

Description: This final report will document the accomplishments of the work of this project. 1) The incremental-iterative (II) form of the reverse-mode (adjoint) method for computing first-order (FO) aerodynamic sensitivity derivatives (SDs) has been successfully implemented and tested in a 2D CFD code (called ANSERS) using the reverse-mode capability of ADIFOR 3.0. These preceding results compared very well with similar SDS computed via a black-box (BB) application of the reverse-mode capability of ADIFOR 3.0, and also with similar SDs calculated via the method of finite differences. 2) Second-order (SO) SDs have been implemented in the 2D ASNWERS code using the very efficient strategy that was originally proposed (but not previously tested) of Reference 3, Appendix A. Furthermore, these SO SOs have been validated for accuracy and computational efficiency. 3) Studies were conducted in Quasi-1D and 2D concerning the smoothness (or lack of smoothness) of the FO and SO SD's for flows with shock waves. The phenomenon is documented in the publications of this study (listed subsequently), however, the specific numerical mechanism which is responsible for this unsmoothness phenomenon was not discovered. 4) The FO and SO derivatives for Quasi-1D and 2D flows were applied to predict aerodynamic design uncertainties, and were also applied in robust design optimization studies.

Description: A study to determine a limiting distance to span ratio for the extrapolation of near-field pressure signatures is described and discussed. This study was to be done in two wind-tunnel facilities with two wind-tunnel models. At this time, only the first half had been completed, so the scope of this report is limited to the design of the models, and to an analysis of the first set of measured pressure signatures. The results from this analysis showed that the pressure signatures measured at separation distances of 2 to 5 span lengths did not show the desired low-boom shapes. However, there were indications that the pressure signature shapes were becoming 'flat-topped'. This trend toward a 'flat-top' pressure signatures shape was seen to be a gradual one at the distance ratios employed in this first series of wind-tunnel tests.

Description: Closed-form direct and inverse kinematics of a new three degree-of-freedom (DOF) parallel manipulator with inextensible limbs and base-mounted actuators are presented. The manipulator has higher resolution and precision than the existing three DOF mechanisms with extensible limbs. Since all of the manipulator actuators are base-mounted; higher payload capacity, smaller actuator sizes, and lower power dissipation can be obtained. The manipulator is suitable for alignment applications where only tip, tilt, and piston motions are significant. The direct kinematics of the manipulator is reduced to solving an eighth-degree polynomial in the square of tangent of half-angle between one of the limbs and the base plane. Hence, there are at most 16 assembly configurations for the manipulator. In addition, it is shown that the 16 solutions are eight pairs of reflected configurations with respect to the base plane. Numerical examples for the direct and inverse kinematics of the manipulator are also presented.

Description: A detailed boundary layer investigation of an NACA 0012 airfoil was made by applying the WIND code at the NASA Glenn Research Center. Structured C-grids of selected configurations were fitted around the airfoil. Grids with dense spacing to capture the boundary layer provided for obtaining accurate results in this region. Concurrent with the boundary layer investigation, contour plots of overall performance of the airfoil were also obtained. A Mach number of 0.80, a Reynolds number of 4,000,000, and a zero angle of attack were set to correspond with numerical and experimental investigations reported in the literature, for comparison with the overall performance. The results were favorable for both the boundary layer investigations and the performance. As WIND evolves, it is intended to run WIND using unstructured grids to compare such results with the current structured grid results.

Description: Langley Stability and Transition Analysis Code (LASTRAC) is a general-purpose, physics-based transition prediction code released by NASA for laminar flow control studies and transition research. This paper describes the LASTRAC extension to general three-dimensional (3D) boundary layers such as finite swept wings, cones, or bodies at an angle of attack. The stability problem is formulated by using a body-fitted nonorthogonal curvilinear coordinate system constructed on the body surface. The nonorthogonal coordinate system offers a variety of marching paths and spanwise waveforms. In the extreme case of an infinite swept wing boundary layer, marching with a nonorthogonal coordinate produces identical solutions to those obtained with an orthogonal coordinate system using the earlier release of LASTRAC. Several methods to formulate the 3D parabolized stability equations (PSE) are discussed. A surface-marching procedure akin to that for 3D boundary layer equations may be used to solve the 3D parabolized disturbance equations. On the other hand, the local line-marching PSE method, formulated as an easy extension from its 2D counterpart and capable of handling the spanwise mean flow and disturbance variation, offers an alternative. A linear stability theory or parabolized stability equations based N-factor analysis carried out along the streamline direction with a fixed wavelength and downstream-varying spanwise direction constitutes an efficient engineering approach to study instability wave evolution in a 3D boundary layer. The surface-marching PSE method enables a consistent treatment of the disturbance evolution along both streamwise and spanwise directions but requires more stringent initial conditions. Both PSE methods and the traditional LST approach are implemented in the LASTRAC.3d code. Several test cases for tapered or finite swept wings and cones at an angle of attack are discussed.

Description: A pressure-sensitive paint (PSP) technique was applied in a wind tunnel experiment in the NASA Langley Research Center 8-Foot Transonic Pressure Tunnel to quantify the vortex-induced surface static pressures on a slender, faceted missile model at subsonic and transonic speeds. Satisfactory global calibrations of the PSP were obtained at =0.70, 0.90, and 1.20, angles of attack from 10 degrees to 20 degrees, and angles of sideslip of 0 and 2.5 degrees using an in-situ method featuring the simultaneous acquisition of electronically-scanned pressures (ESP) at 57 discrete locations on the model. Both techniques clearly revealed the significant influence on the surface pressure distributions of the vortices shed from the sharp, chine-like leading edges. The mean error in the PSP measurements relative to the ESP data was approximately 0.6 percent at M infinity =0.70 and 2.6 percent at M infinity =0.90 and 1.20. The vortex surface pressure signatures obtained from the PSP and ESP techniques were correlated with the off-surface vortex cross-flow structures obtained using a laser vapor screen (LVS) flow visualization technique. The on-surface and off-surface techniques were complementary, since each provided details of the vortex-dominated flow that were not clear or apparent in the other.

Description: A pressure-sensitive paint (PSP) technique was applied in a wind tunnel experiment in the NASA Langley Research Center 8-Foot Transonic Pressure Tunnel to quantify the vortex-induced surface static pressures on a slender, faceted missile model at subsonic and transonic speeds. Global PSP calibrations were obtained using an in-situ method featuring the simultaneous electronically-scanned pressures (ESP) measurements. Both techniques revealed the significant influence leading-edge vortices on the surface pressure distributions. The mean error in the PSP measurements relative to the ESP data was approximately 0.6 percent at M(sub infinity)=0.70 and 2.6 percent at M(sub infinity)=0.90 and 1.20. The vortex surface pressure signatures obtained from the PSP and ESP techniques were correlated with the off-surface vortex cross-flow structures obtained using a laser vapor screen (LVS) flow visualization technique. The on-surface and off-surface techniques were complementary, since each provided details of the vortex-dominated flow that were not clear or apparent in the other.

Description: Time scales in turbulent and transitional flow provide a link between experimental data and modeling, both in terms of physical content and for quantitative assessment. The problem of interest here is the definition of time scales in an unsteady flow. Using representative samples of data from GEAE low pressure turbine experiment in low speed research turbine facility with wake-induced transition, we document several methods to extract dominant frequencies, and compare the results. We show that conventional methods of time scale evaluation (based on autocorrelation functions and on Fourier spectra) and wavelet-based methods provide similar information when applied to stationary signals. We also show the greater flexibility of the wavelet-based methods when dealing with intermittent or strongly modulated data, as are encountered in transitioning boundary layers and in flows with unsteady forcing associated with wake passing. We define phase-averaged dominant frequencies that characterize the turbulence associated with freestream conditions and with the passing wakes downstream of a rotor. The relevance of these results for modeling is discussed in the paper.

Description: A computational investigation was performed to assess the effectiveness of Detached-Eddy Simulation (DES) as a tool for predicting icing effects. The AVUS code, a public domain flow solver, was employed to compute solutions for an iced wing configuration using DES and steady Reynolds Averaged Navier-Stokes (RANS) equation methodologies. The configuration was an extruded GLC305/944-ice shape section with a rectangular planform. The model was mounted between two walls so no tip effects were considered. The numerical results were validated by comparison with experimental data for the same configuration. The time-averaged DES computations showed some improvement in lift and drag results near stall when compared to steady RANS results. However, comparisons of the flow field details did not show the level of agreement suggested by the integrated quantities. Based on our results, we believe that DES may prove useful in a limited sense to provide analysis of iced wing configurations when there is significant flow separation, e.g., near stall, where steady RANS computations are demonstrably ineffective. However, more validation is needed to determine what role DES can play as part of an overall icing effects prediction strategy. We conclude the report with an assessment of existing computational tools for application to the iced wing problem and a discussion of issues that merit further study.

Description: This report presents results of the work completed on the preliminary design of Fan 3 of NASA s 22-inch Fan Low Noise Research project. Fan 3 was intended to build on the experience gained from Fans 1 and 2 by demonstrating noise reduction technology that surpasses 1992 levels by 6 dB. The work was performed as part of NASA s Advanced Subsonic Technology (AST) program. Work on this task was conducted in the areas of CFD code validation, acoustic prediction and validation, rotor parametric studies, and fan exit guide vane (FEGV) studies up to the time when a NASA decision was made to cancel the design, fabrication and testing phases of the work. The scope of the program changed accordingly to concentrate on two subtasks: (1) Rig data analysis and CFD code validation and (2) Fan and FEGV optimization studies. The results of the CFD code validation work showed that this tool predicts 3D flowfield features well from the blade trailing edge to about a chord downstream. The CFD tool loses accuracy as the distance from the trailing edge increases beyond a blade chord. The comparisons of noise predictions to rig test data showed that both the tone noise tool and the broadband noise tool demonstrated reasonable agreement with the data to the degree that these tools can reliably be used for design work. The section on rig airflow and inlet separation analysis describes the method used to determine total fan airflow, shows the good agreement of predicted boundary layer profiles to measured profiles, and shows separation angles of attack ranging from 29.5 to 27deg for the range of airflows tested. The results of the rotor parametric studies were significant in leading to the decision not to pursue a new rotor design for Fan 3 and resulted in recommendations to concentrate efforts on FEGV stator designs. The ensuing parametric study on FEGV designs showed the potential for 8 to 10 EPNdB noise reduction relative to the baseline.

Description: Present simulation of turbulent flows involving shock wave/boundary layer interaction invariably overestimates heat flux by almost a factor of two. One possible reason for such a performance is a result of the fact that the turbulence models employed make use of Morkovin's hypothesis. This hypothesis is valid for non-hypersonic Mach numbers and moderate rates of heat transfer. At hypersonic Mach numbers, high rates of heat transfer exist in regions where shock wave/boundary layer interactions are important. As a result, one should not expect traditional turbulence models to yield accurate results. The goal of this investigation is to explore the role of a variable Prandtl number formulation in predicting heat flux in flows dominated by strong shock wave/boundary layer interactions. The intended applications involve external flows in the absence of combustion such as those encountered in supersonic inlets. This can be achieved by adding equations for the temperature variance and its dissipation rate. Such equations can be derived from the exact Navier-Stokes equations. Traditionally, modeled equations are based on the low speed energy equation where the pressure gradient term and the term responsible for energy dissipation are ignored. It is clear that such assumptions are not valid for hypersonic flows. The approach used here is based on the procedure used in deriving the k-zeta model, in which the exact equations that governed k, the variance of velocity, and zeta, the variance of vorticity, were derived and modeled. For the variable turbulent Prandtl number, the exact equations that govern the temperature variance and its dissipation rate are derived and modeled term by term. The resulting set of equations are free of damping and wall functions and are coordinate-system independent. Moreover, modeled correlations are tensorially consistent and invariant under Galilean transformation. The final set of equations will be given in the paper.

Description: A mode choice model that generates on-demand air travel forecasts at a set of GA airports based on changes in economic characteristics, vehicle performance characteristics such as speed and cost, and demographic trends has been integrated with a model to generate itinerate aircraft operations by airplane category at a set of 3227 airports. Numerous intermediate outputs can be generated, such as the number of additional trips diverted from automobiles and schedule air by the improved performance and cost of on-demand air vehicles. The total number of transported passenger miles that are diverted is also available. From these results the number of new aircraft to service the increased demand can be calculated. Output from the models discussed is in the format to generate the origin and destination traffic flow between the 3227 airports based on solutions to a gravity model.

Description: This paper describes the avian wing geometry (Seagull, Merganser, Teal and Owl) extracted from non-contact surface measurements using a three-dimensional laser scanner. The geometric quantities, including the camber line and thickness distribution of airfoil, wing planform, chord distribution, and twist distribution, are given in convenient analytical expressions. Thus, the avian wing surfaces can be generated and the wing kinematics can be simulated. The aerodynamic characteristics of avian airfoils in steady inviscid flows are briefly discussed. The avian wing kinematics is recovered from videos of three level-flying birds (Crane, Seagull and Goose) based on a two-jointed arm model. A flapping seagull wing in the 3D physical space is re-constructed from the extracted wing geometry and kinematics.

Description: Several flight-critical aerodynamic problems such as buffet, flutter, stall, and wing rock are strongly affected or caused by abrupt changes in unsteady aerodynamic loads and moments. Advanced sensing and flow diagnostic techniques have made possible simultaneous identification and tracking, in real-time, of the critical surface, viscosity-related aerodynamic phenomena under both steady and unsteady flight conditions. The wind tunnel study reported here correlates surface hot-film measurements of leading edge stagnation point and separation point, with unsteady aerodynamic loads on a NACA 0015 airfoil. Lift predicted from the correlation model matches lift obtained from pressure sensors for an airfoil undergoing harmonic pitchup and pitchdown motions. An analytical model was developed that demonstrates expected stall trends for pitchup and pitchdown motions. This report demonstrates an ability to obtain unsteady aerodynamic loads in real-time, which could lead to advances in air vehicle safety, performance, ride-quality, control, and health management.

Description: A genetic algorithm approach suitable for solving multi-objective optimization problems is described and evaluated using a series of aerodynamic shape optimization problems. Several new features including two variations of a binning selection algorithm and a gene-space transformation procedure are included. The genetic algorithm is suitable for finding pareto optimal solutions in search spaces that are defined by any number of genes and that contain any number of local extrema. A new masking array capability is included allowing any gene or gene subset to be eliminated as decision variables from the design space. This allows determination of the effect of a single gene or gene subset on the pareto optimal solution. Results indicate that the genetic algorithm optimization approach is flexible in application and reliable. The binning selection algorithms generally provide pareto front quality enhancements and moderate convergence efficiency improvements for most of the problems solved.

Description: The aerodynamic drag of an aerial refueling assembly was calculated during the Automated Aerial Refueling project at the NASA Dryden Flight Research Center. An F/A-18A airplane was specially instrumented to obtain accurate fuel flow measurements and to determine engine thrust. A standard Navy air refueling store with a retractable refueling hose and paradrogue was mounted to the centerline pylon of the F/A-18A airplane. As the paradrogue assembly was deployed and stowed, changes in the calculated thrust of the airplane occurred and were equated to changes in vehicle drag. These drag changes were attributable to the drag of the paradrogue assembly. The drag of the paradrogue assembly was determined to range from 200 to 450 lbf at airspeeds from 170 to 250 KIAS. Analysis of the drag data resulted in a single drag coefficient of 0.0056 for the paradrogue assembly that adequately matched the calculated drag for all flight conditions. The drag relief provided to the tanker airplane when a receiver airplane engaged the paradrogue is also documented from 35 to 270 lbf at the various flight conditions tested. The results support the development of accurate aerodynamic models to be used in refueling simulations and control laws for fully autonomous refueling.

Description: Development of improved modeling methods to provide increased fidelity of flight predictions for aircraft motions during flight in flow regimes with large nonlinearities requires improvements in test techniques for measuring and characterizing wind tunnel data. This paper presents a method for providing a measure of data integrity for static and forced oscillation test techniques. Data integrity is particularly important when attempting to accurately model and predict flight of today s high performance aircraft which are operating in expanded flight envelopes, often maneuvering at high angular rates at high angles-of-attack, even above maximum lift. Current aerodynamic models are inadequate in predicting flight characteristics in the expanded envelope, such as rapid aircraft departures and other unusual motions. Present wind tunnel test methods do not factor changes of flow physics into data acquisition schemes, so in many cases data are obtained over more iterations than required, or insufficient data may be obtained to determine a valid estimate with statistical significance. Additionally, forced oscillation test techniques, one of the primary tools used to develop dynamic models, do not currently provide estimates of the uncertainty of the results during an oscillation cycle. A method to optimize the required number of forced oscillation cycles based on decay of uncertainty gradients and balance tolerances is also presented.

Description: An overview of basic aeroheating design issues for Reusable Launch Vehicles (RLV), which addresses the application of hypersonic ground-based testing, and computational fluid dynamic (CFD) and engineering codes, is presented. Challenges inherent to the prediction of aeroheating environments required for the successful design of the RLV Thermal Protection System (TPS) are discussed in conjunction with the importance of employing appropriate experimental/computational tools. The impact of the information garnered by using these tools in the resulting analyses, ultimately enhancing the RLV TPS design is illustrated. A wide range of topics is presented in this overview; e.g. the impact of flow physics issues such as boundary-layer transition, including effects of distributed and discrete roughness, shockshock interactions, and flow separation/reattachment. Also, the benefit of integrating experimental and computational studies to gain an improved understanding of flow phenomena is illustrated. From computational studies, the effect of low-density conditions and of uncertainties in material surface properties on the computed heating rates are highlighted as well as the significant role of CFD in improving the Outer Mold Line (OML) definition to reduce aeroheating while maintaining aerodynamic performance. Appropriate selection of the TPS design trajectories and trajectory shaping to mitigate aeroheating levels and loads are discussed. Lastly, an illustration of an aeroheating design process is presented whereby data from hypersonic wind-tunnel tests are integrated with predictions from CFD codes and engineering methods to provide heating environments along an entry trajectory as required for TPS design.

Description: This paper is a description of the analysis of blockage corrections for bodies of revolution for the slotted-wall configuration of the National Transonic Facility (NTF) at the NASA Langley Research Center (LaRC). A wall correction method based on the measured wall signature is used. Test data from three different-sized blockage bodies and four wall ventilation settings were analyzed at various Mach numbers and unit Reynolds numbers. The results indicate that with the proper selection of the boundary condition parameters, the wall correction method can predict blockage corrections consistent with the wall measurements for Mach numbers as high as 0.95.

Description: A multigrid solver is defined as having textbook multigrid efficiency (TME) if the solutions to the governing system of equations are attained in a computational work which is a small (less than 10) multiple of the operation count in evaluating the discrete residuals. TME in solving the incompressible inviscid fluid equations is demonstrated for leading- edge stagnation flows. The contributions of this paper include (1) a special formulation of the boundary conditions near stagnation allowing convergence of the Newton iterations on coarse grids, (2) the boundary relaxation technique to facilitate relaxation and residual restriction near the boundaries, (3) a modified relaxation scheme to prevent initial error amplification, and (4) new general analysis techniques for multigrid solvers. Convergence of algebraic errors below the level of discretization errors is attained by a full multigrid (FMG) solver with one full approximation scheme (F.4S) cycle per grid. Asymptotic convergence rates of the F.4S cycles for the full system of flow equations are very fast, approaching those for scalar elliptic equations.

Description: A structured-grid quality measure is proposed, combining three traditional measurements: intersection angles, stretching, and curvature. Quality assesses whether the grid generated provides the best possible tradeoffs in grid stretching and skewness that enable accurate flow predictions, whereas the grid density is assumed to be a constraint imposed by the available computational resources and the desired resolution of the flow field. The usefulness of this quality measure is assessed by comparing heat transfer predictions from grid convergence studies for grids of varying quality in the range of [0.6-0.8] on an 8'half-angle sphere-cone, at laminar, perfect gas, Mach 10 wind tunnel conditions.

Description: Development of a general aerodynamic model that is adequate for predicting the forces and moments in the nonlinear and unsteady portions of the flight envelope has not been accomplished to a satisfactory degree. Predicting aerodynamic response during arbitrary motion of an aircraft over the complete flight envelope requires further development of the mathematical model and the associated methods for ground-based testing in order to allow identification of the model. In this study, a general nonlinear unsteady aerodynamic model is presented, followed by a summary of a linear modeling methodology that includes test and identification methods, and then a progressive series of steps suggesting a roadmap to develop a general nonlinear methodology that defines modeling, testing, and identification methods. Initial steps of the general methodology were applied to static and oscillatory test data to identify rolling-moment coefficient. Static measurements uncovered complicated dependencies of the aerodynamic coefficient on angle of attack and sideslip in the stall region making it difficult to find a simple analytical expression for the measurement data. In order to assess the effect of sideslip on the damping and unsteady terms, oscillatory tests in roll were conducted at different values of an initial offset in sideslip. Candidate runs for analyses were selected where higher order harmonics were required for the model and where in-phase and out-of-phase components varied with frequency. From these results it was found that only data in the angle-of-attack range of 35 degrees to 37.5 degrees met these requirements. From the limited results it was observed that the identified models fit the data well and both the damping-in-roll and the unsteady term gain are decreasing with increasing sideslip and motion amplitude. Limited similarity between parameter values in the nonlinear model and the linear model suggest that identifiability of parameters in both terms may be a problem. However, the proposed methodology can still be used with careful experiment design and carefully selected values of angle of attack, sideslip, amplitude, and frequency of the oscillatory data.

Description: The ability of the Laura flow solver to predict local heating augmentation factors for shallow cavities is assessed. This assessment is part of a larger e ort within the Space Shuttle return-to-flight program to develop technologies to support on-orbit tile repair decisions. The comparison is made against global phosphor thermography images taken in the Langley Aerothermodynamic Laboratory 20-Inch Mach 6 Air Tunnel. The cavities are rectangular in shape, with lengths L/H of 14 20 and depths H/ of 1.1 5.2. The fully laminar results, for Re = 300, show good agreement between the data sets. For Re = 503, the wind tunnel data indicates boundary layer transition with turbulent flow both within and downstream of the cavity. The turbulent flow structures are significantly di erent from the laminar predictions, with order of magnitude increases in the heating augmentations. Because of the di erent flow structures, no simple bump factor can be used to correct the laminar calculations to account for the turbulent heating levels. A fine gradation in wind tunnel cases will be required to clearly delineate the laminar-to-turbulent transition point, and hence the limits of applicability of the laminar numerical approach.

Description: A simple optical technique for the structural characterization of lightweight space columns is presented. The technique is useful for determining the coefficient of thermal expansion during cool down as well as the induced strain during tension and compression testing. The technique is based upon object-to-image plane scaling and does not require any photogrammetric calibrations or computations. Examples of the measurement of the coefficient of thermal expansion are presented for several lightweight space columns. Examples of strain measured during tension and compression testing are presented along with comparisons to results obtained with Linear Variable Differential Transformer (LVDT) position transducers.

Description: The capabilities of a relatively new direct simulation Monte Carlo (DSMC) code are examined for the problem of hypersonic laminar shock/shock and shock/boundary layer interactions, where boundary layer separation is an important feature of the flow. Flow about two model configurations is considered, where both configurations (a biconic and a hollow cylinder-flare) have recent published experimental measurements. The computations are made by using the DS2V code of Bird, a general two-dimensional/axisymmetric time accurate code that incorporates many of the advances in DSMC over the past decade. The current focus is on flows produced in ground-based facilities at Mach 12 and 16 test conditions with nitrogen as the test gas and the test models at zero incidence. Results presented highlight the sensitivity of the calculations to grid resolutions, sensitivity to physical modeling parameters, and comparison with experimental measurements. Information is provided concerning the flow structure and surface results for the extent of separation, heating, pressure, and skin friction.

Description: A series of active flow control experiments were recently conducted on a simplified high-lift system. The purpose of the experiments was to explore the prospects of eliminating all but simply hinged leading and trailing edge flaps, while controlling separation on the supercritical airfoil using multiple periodic excitation slots. Excitation was provided by three. independently controlled, self-contained, piezoelectric actuators. Low frequency excitation was generated through amplitude modulation of the high frequency carrier wave, the actuators' resonant frequencies. It was demonstrated, for the first time, that pulsed modulated signal from two neighboring slots interact favorably to increase lift. Phase sensitivity at the low frequency was measured, even though the excitation was synthesized from the high-frequency carrier wave. The measurements were performed at low Reynolds numbers and included mean and unsteady surface pressures, surface hot-films, wake pressures and particle image velocimetry. A modest (6%) increase in maximum lift (compared to the optimal baseline) was obtained due t o the activation of two of the three actuators.

Description: The capabilities of a relatively new direct simulation Monte Carlo (DSMC) code are examined for the problem of hypersonic laminar shock/shock and shock/boundary layer interactions, where boundary layer separation is an important feature of the flow. Flow about two model configurations is considered, where both configurations (a biconic and a hollow cylinder-flare) have recent published experimental measurements. The computations are made by using the DS2V code of Bird, a general two-dimensional/axisymmetric time accurate code that incorporates many of the advances in DSMC over the past decade. The current focus is on flows produced in ground-based facilities at Mach 12 and 16 test conditions with nitrogen as the test gas and the test models at zero incidence. Results presented highlight the sensitivity of the calculations to grid resolution, sensitivity to physical modeling parameters, and comparison with experimental measurements. Information is provided concerning the flow structure and surface results for the extent of separation, heating, pressure, and skin friction.

Description: Differential Evolution (DE) is a simple and robust evolutionary strategy that has been provEn effective in determining the global optimum for several difficult optimization problems. Although DE offers several advantages over traditional optimization approaches, its use in applications such as aerodynamic shape optimization where the objective function evaluations are computationally expensive is limited by the large number of function evaluations often required. In this paper various approaches for improving the efficiency of DE are reviewed and discussed. Several approaches that have proven effective for other evolutionary algorithms are modified and implemented in a DE-based aerodynamic shape optimization method that uses a Navier-Stokes solver for the objective function evaluations. Parallelization techniques on distributed computers are used to reduce turnaround times. Results are presented for standard test optimization problems and for the inverse design of a turbine airfoil. The efficiency improvements achieved by the different approaches are evaluated and compared.

Description: The Videogrammetric Model Deformation (VMD) technique, developed at NASA Langley Research Center, was recently used to measure displacements and local surface angle changes on a static aeroelastic wind-tunnel model. The results were assessed for consistency, accuracy and usefulness. Vertical displacement measurements and surface angular deflections (derived from vertical displacements) taken at no-wind/no-load conditions were analyzed. For accuracy assessment, angular measurements were compared to those from a highly accurate accelerometer. Shewhart's Variables Control Charts were used in the assessment of consistency and uncertainty. Some bad data points were discovered, and it is shown that the measurement results at certain targets were more consistent than at other targets. Physical explanations for this lack of consistency have not been determined. However, overall the measurements were sufficiently accurate to be very useful in monitoring wind-tunnel model aeroelastic deformation and determining flexible stability and control derivatives. After a structural model component failed during a highly loaded condition, analysis of VMD data clearly indicated progressive structural deterioration as the wind-tunnel condition where failure occurred was approached. As a result, subsequent testing successfully incorporated near- real-time monitoring of VMD data in order to ensure structural integrity. The potential for higher levels of consistency and accuracy through the use of statistical quality control practices are discussed and recommended for future applications.

Description: A full-engine simulation of the three-dimensional flow in the GE90 94B high-bypass ratio turbofan engine has been achieved. It would take less than 11 hr of wall clock time if starting from scratch through the exploitation of parallel processing. The simulation of the compressor components, the cooled high-pressure turbine, and the low-pressure turbine was performed using the APNASA turbomachinery flow code. The combustor flow and chemistry were simulated using the National Combustor Code (NCC). The engine simulation matches the engine thermodynamic cycle for a sea-level takeoff condition. The simulation is started at the inlet of the fan and progresses downstream. Comparisons with the cycle point are presented. A detailed look at the blockage in the turbomachinery is presented as one measure to assess and view the solution and the multistage interaction effects.

Description: Aerodynamic design is an iterative process involving geometry manipulation and complex computational analysis subject to physical constraints and aerodynamic objectives. A design cycle consists of first establishing the performance of a baseline design, which is usually created with low-fidelity engineering tools, and then progressively optimizing the design to maximize its performance. Optimization techniques have evolved from relying exclusively on designer intuition and insight in traditional trial and error methods, to sophisticated local and global search methods. Recent attempts at automating the search through a large design space with formal optimization methods include both database driven and direct evaluation schemes. Databases are being used in conjunction with surrogate and neural network models as a basis on which to run optimization algorithms. Optimization algorithms are also being driven by the direct evaluation of objectives and constraints using high-fidelity simulations. Surrogate methods use data points obtained from simulations, and possibly gradients evaluated at the data points, to create mathematical approximations of a database. Neural network models work in a similar fashion, using a number of high-fidelity database calculations as training iterations to create a database model. Optimal designs are obtained by coupling an optimization algorithm to the database model. Evaluation of the current best design then gives either a new local optima and/or increases the fidelity of the approximation model for the next iteration. Surrogate methods have also been developed that iterate on the selection of data points to decrease the uncertainty of the approximation model prior to searching for an optimal design. The database approximation models for each of these cases, however, become computationally expensive with increase in dimensionality. Thus the method of using optimization algorithms to search a database model becomes problematic as the number of design variables is increased.

Description: The achievement of flight with manned, powered, heavier-than-air aircraft in 1903 marked the beginning of a new era in the means of transportation. A special advantage for aircraft was in speed. However, when an aircraft penetrates the air at very high speeds, the disturbed air is compressed and there are changes in the density, pressure and temperature of the air. These compressibility effects change the aerodynamic characteristics of an aircraft and introduce problems in drag, stability and control. Many aircraft designed in the post-World War II era were plagued with the effects of compressibility. Accordingly, the study of the aerodynamic behavior of aircraft, spacecraft and missiles at high-speed became a major part of the research activity of the NACA/NASA. The intent of the research was to determine the causes and provide some solutions for the aerodynamic problems resulting from the effects of compressibility. The purpose of this paper is to review some of the high-speed aerodynamic research work conducted at the Langley Research Center from the viewpoint of the author who has been active in much of the effort.

Description: The procedure for investigating the effect of deviation from flush mounting of pressure transducers on the exterior of Tu-144LL in flight is described. Experimental data in the mach-number range 0.58 - 2.0 are presented for distortion of the measured wall-pressure fluctuation spectra of the turbulent boundary layer by recessed and protruding transducers. The results of flight experiments are compared with data of wind tunnel experiments. The distortion of measured turbulent boundary layer wall pressure fluctuations caused by transducer-surface deviation from the surrounding surface as a function of dimensionless parameters is predicted and presented on the basis of dimensional analysis.

Description: Blade section normal force and pitching moment were investigated for six rotors operating at transition and high speeds: H-34 in flight and wind tunnel, SA 330 (research Puma), SA 349/2, UH-60A full-scale and BO-105 model (HART-I). The measured data from flight and wind tunnel tests were compared with calculations obtained using the comprehensive analysis CAMRAD II. The calculations were made using two free wake models: rolled-up and multiple-trailer with consolidation models. At transition speed, there is fair to good agreement for the blade section normal force between the test data and analysis for the H-34, research Puma, and SA 349/2 with the rolled-up wake. The calculated airloads differ significantly from the measurements for the UH-60A and BO-105. Better correlation is obtained for the UH-60A and BO-105 by using the multiple-trailer with consolidation wake model. In the high speed condition, the analysis shows generally good agreement with the research Puma flight data in both magnitude and phase. However, poor agreement is obtained for the other rotors examined. The analysis shows that the aerodynamic tip design (chord length and quarter chord location) of the Puma has an important influence on the phase correlation.

Description: This paper deals with the computational analysis of buoyancy effects in the near field of an isothermal helium jet injected into quiescent ambient air environment. The transport equations of helium mass fraction coupled with the conservation equations of mixture mass and momentum were solved using a staggered grid finite volume method. Laminar, axisymmetric, unsteady flow conditions were considered for the analysis. An orthogonal system with non-uniform grids was used to capture the instability phenomena. Computations were performed for Earth gravity and during transition from Earth to different gravitational levels. The flow physics was described by simultaneous visualizations of velocity and concentration fields at Earth and microgravity conditions. Computed results were validated by comparing with experimental data substantiating that buoyancy induced global flow oscillations present in Earth gravity are absent in microgravity. The dependence of oscillation frequency and amplitude on gravitational forcing was presented to further quantify the buoyancy effects.

Description: Rotor active control using higher harmonic blade pitch has been proposed as a means to reduce both rotor radiated noise and airframe vibration and to enhance rotor performance. The higher harmonic input, however, can affect rotor thrust and cyclic flapping - the basic trim characteristics of the rotor. Some of the trim changes can negate the active control benefits. For example, wind tunnel test results of a full scale BO-105 rotor with individual-blade control indicate some rotor performance improvements, accompanied with changes in rotor trim, using two-per-rev blade pitch input. The observed performance benefits could therefore be a simple manifestation of the trim change rather than an efficient redistribution of the rotor airloads. More recently, the flight test of the BO-105 helicopter equip,ped with individual-blade-control actuators also reported trim changes whenever the two-per-rev blade pitch for noise reduction was activated. The pilot had to adjust the trim control to maintain the aircraft under a constant flight path. These two cases highlight the, importance of trim considerations in the application of active control to rotorcraft.

Description: LASTRAC is a general-purposed, physics-based transition prediction code released by NASA for Laminar Flow Control studies and transition research. The design and development of the LASTRAC code is aimed at providing an engineering tool that is easy to use and yet capable of dealing with a broad range of transition related issues. It was written from scratch based on the state-of-the-art numerical methods for stability analysis and modern software technologies. At low fidelity, it allows users to perform linear stability analysis and N-factor transition correlation for a broad range of flow regimes and configurations by using either the linear stability theory or linear parabolized stability equations method. At high fidelity, users may use nonlinear PSE to track finite-amplitude disturbances until the skin friction rise. This document describes the governing equations, numerical methods, code development, detailed description of input/output parameters, and case studies for the current release of LASTRAC.

Description: Numerical simulations are used to examine the roughness-induced transient growth in a laminar boundary-layer flow. Based on the spectral element method, these simulations model the stationary disturbance field associated with a nonsmooth roughness geometry, such as the spanwise periodic array of circular disks used by White and co-workers during a series of wind tunnel experiments at Case Western Reserve University. Besides capturing the major trends from the recent measurements by White and Ergin, the simulations provide additional information concerning the relative accuracy of the experimental findings derived from two separate wall-finding procedures. The paper also explores the dependence of transient growth on geometric characteristics of the roughness distribution, including the height and planform shape of the roughness element and the ratio of roughness due to spacing between an adjacent pair of elements. Results are used for a preliminary assessment of the differences between recently reported theoretical results of Tumin and Reshotko and the measurements by White and Ergin.

Description: Inviscid aerodynamic characteristics of the Space Shuttle Orbiter were computed in support of the Columbia Accident Investigation. The unstructured grid software FELISA was used and computations were done using freestream conditions corresponding to those in the NASA Langley 20-Inch Mach 6 CF4 tunnel test section. The angle of attack was held constant at 40 degrees. The baseline (undamaged) configuration and a large number of damaged configurations of the Orbiter were studied. Most of the computations were done on a half model. However, one set of computations was done using the full-model to study the effect of sideslip. The differences in the aerodynamic coefficients for the damaged and the baseline configurations were computed. Simultaneously with the computation reported here, tests were being done on a scale model of the Orbiter in the 20-Inch Mach 6 CF4 tunnel to measure the deltas . The present computations complemented the CF4 tunnel test, and provided aerodynamic coefficients of the Orbiter as well as its components. Further, they also provided details of the flow field.

Description: An experimental influence coefficient technique was used to obtain unsteady aerodynamic influence coefficients and, consequently, unsteady pressures for a cascade of symmetric airfoils oscillating in pitch about mid-chord. Stagger angles of 0 deg and 10 deg were investigated for a cascade with a gap-to-chord ratio of 0.417 operating at an axial Mach number of 1.9, resulting in a supersonic leading-edge locus. Reduced frequencies ranged from 0.056 to 0.2. The influence coefficients obtained determine the unsteady pressures for any interblade phase angle. The unsteady pressures were compared with those predicted by several algorithms for interblade phase angles of 0 deg and 180 deg.

Description: A fast multigrid solver for the steady, incompressible Navier-Stokes equations is presented. The multigrid solver is based upon a factorizable discrete scheme for the velocity-pressure form of the Navier-Stokes equations. This scheme correctly distinguishes between the advection-diffusion and elliptic parts of the operator, allowing efficient smoothers to be constructed. To evaluate the multigrid algorithm, solutions are computed for flow over a flat plate, parabola, and a Karman-Trefftz airfoil. Both nonlifting and lifting airfoil flows are considered, with a Reynolds number range of 200 to 800. Convergence and accuracy of the algorithm are discussed. Using Gauss-Seidel line relaxation in alternating directions, multigrid convergence behavior approaching that of O(N) methods is achieved. The computational efficiency of the numerical scheme is compared with that of Runge-Kutta and implicit upwind based multigrid methods.

Description: Software was developed to construct approximating NURBS curves for iced airfoil geometries. Users specify a tolerance that determines the extent to which the approximating curve follows the rough ice. The user can therefore smooth the ice geometry in a controlled manner, thereby enabling the generation of grids suitable for numerical aerodynamic simulations. Ultimately, this ability to smooth the ice geometry will permit studies of the effects of smoothing upon the aerodynamics of iced airfoils. The software was applied to several different types of iced airfoil data collected in the Icing Research Tunnel at NASA Glenn Research Center, and in all cases was found to efficiently generate suitable approximating NURBS curves. This method is an improvement over the current "control point formulation" of Smaggice (v.1.2). In this report, we present the relevant theory of approximating NURBS curves and discuss typical results of the software.

Description: A control surface for an air vehicle (e.g., an aircraft, rocket, or missile) is useful for flight control at both subsonic and supersonic speeds. The control surface defines the outboardmost tip of a flight structure (e.g., a wing, tail or other stabilizer) of the air vehicle. Hence, the control surface is referred to as a `tiperon`. The tiperon has an approximately L-shaped configuration, and can be rotated relative to a fixed portion of the flight structure about a control axis. The respective surface areas of the tiperon sections forward and aft of the control axis are proportioned to place the subsonic center of pressure aft of the control axis to enhance aircraft control, and preferably also forward of the centroid of tiperon surface area. Also, the control surface sections forward and aft of the control axis are preferably mass-balanced, or at least nearly so, to enhance aircraft control at supersonic speeds. Either of the tiperon sections forward and aft of the control axis can be tapered to reduce the dependence of the moment exerted by air flow about the control axis, upon the tiperon's angle-of-attack. The tiperon also has enough surface area to control the air vehicle, even at low airspeeds. The invention is also directed to air vehicles incorporating one or more such control surfaces.

Description: The vortex flowfield over an advanced twin-tailed fighter configuration was measured in a low-speed wind tunnel at two angles of attack. The primary test data consisted of 3-component velocity surveys obtained using a Laser Doppler Velocimeter. Laser light sheet and surface flow visualization were also obtained to provide insight into the flowfield structure. Time-averaged velocities and the root mean square of the velocity fluctuations were obtained at two cross-sections above the model. At 15 degrees angle of attack, the vortices generated by the wing leading edge extension (LEX) were unburst over the model and passed outboard of the vertical tail. At 25 degrees angle of attack, the vortices burst in the vicinity of the wing-LEX intersection and impact directly on the vertical tails. The RMS levels of the velocity fluctuations reach values of approximately 30% in the region of the vertical tails.

Description: The efficient prediction of helicopter rotor performance, vibratory loads, and aeroelastic properties still relies heavily on the use of comprehensive analysis codes by the rotorcraft industry. These comprehensive codes utilize look-up tables to provide two-dimensional aerodynamic characteristics. Typically these tables are comprised of a combination of wind tunnel data, empirical data and numerical analyses. The potential to rely more heavily on numerical computations based on Computational Fluid Dynamics (CFD) simulations has become more of a reality with the advent of faster computers and more sophisticated physical models. The ability of five different CFD codes applied independently to predict the lift, drag and pitching moments of rotor airfoils is examined for the SC1095 airfoil, which is utilized in the UH-60A main rotor. Extensive comparisons with the results of ten wind tunnel tests are performed. These CFD computations are found to be as good as experimental data in predicting many of the aerodynamic performance characteristics. Four turbulence models were examined (Baldwin-Lomax, Spalart-Allmaras, Menter SST, and k-omega).

Description: Experiments were conducted in earth gravity and micro gravity to acquire quantitative data on near field flow structure of helium jets injected into air. Microgravity conditions were simulated in the 2.2 s drop tower at NASA John H. Glenn Research Center. The jet flow was observed by quantitative rainbow schlieren deflectometry, a non-intrusive line of sight measurement technique suited for the microgravity environment. The flow structure was characterized by distribution of helium mole fraction obtained from color schlieren images taken at 60 Hz. Results show that the jet in microgravity was up to 70% wider than that in Earth gravity. Experiments reveal that the global flow oscillations observed in Earth are absent in microgravity. Quantitative deatails are provided of the evolution as the experiment undergoes changes in gravity in the drop tower.

Description: A pressure-sensitive paint (PSP) technique was applied in a wind tunnel experiment in the NASA Langley Research Center 8-Foot Transonic Pressure Tunnel to quantify the effect of wing fillets on the global vortex-induced surface static pressure field about a sharp leading-edge 76o/40o double delta wing, or strake-wing, model at subsonic and transonic speeds. Global calibrations of the PSP were obtained at M = 0.50, 0.70, 0.85, 0.95, and 1.20, a Reynolds number per unit length of 2.0 million, and angles of attack from 10 degrees to 20 degrees using an in-situ method featuring the simultaneous acquisition of electronically-scanned pressures (ESP) at discrete locations on the model. The mean error in the PSP measurements relative to the ESP data was approximately 2 percent or less at M = 0.50 to 0.85 but increased to several percent at M =0.95 and 1.20. The PSP pressure distributions and pseudo-colored planform view pressure maps clearly revealed the vortex-induced pressure signatures at all Mach numbers and angles of attack. Small fillets having a parabolic or diamond planform situated at the strake-wing intersection were designed to manipulate the vortical flows by, respectively, removing the leading-edge discontinuity or introducing additional discontinuities. The fillets caused global changes in the vortex-dominated surface pressure field that were effectively captured in the PSP measurements. The vortex surface pressure signatures were compared to available off-surface vortex cross-flow structures obtained using a laser vapor screen (LVS) flow visualization technique. The fillet effects on the PSP pressure distributions and the observed leading-edge vortex flow characteristics were consistent with the trends in the measured lift, drag, and pitching moment coefficients.

Description: NASA has conducted a series of Computational Aeroacoustics (CAA) Workshops on Benchmark Problems to develop a set of realistic CAA problems that can be used for code validation. In the Third (1999) and Fourth (2003) Workshops, the single airfoil gust response problem, with real geometry effects, was included as one of the benchmark problems. Respondents were asked to calculate the airfoil RMS pressure and far-field acoustic intensity for different airfoil geometries and a wide range of gust frequencies. This paper presents the validated that have been obtained to the benchmark problem, and in addition, compares them with classical flat plate results. It is seen that airfoil geometry has a strong effect on the airfoil unsteady pressure, and a significant effect on the far-field acoustic intensity. Those parts of the benchmark problem that have not yet been adequately solved are identified and presented as a challenge to the CAA research community.

Description: A proposed new class of aircraft - the Advanced Theater Transport (ATT) will combine strategic range and high payload with 'Super-STOL' (short take-off and landing) capability. It is also proposed to modify a YC-15 into a technology demonstrator with a 20-deg tilt wing; four, eight-bladed propellers; cross-shafted gearboxes and V-22 engines. These constitute a unique combination of design features that potentially affect performance, loads and whirl-mode stability (whirl flutter). NASA Ames Research Center is working with Boeing and Hamilton Sundstrand on technology challenges presented by the concept; the purpose of NASA involvement is to establish requirements for the demonstrator and for early design guidance, with emphasis on whirl flutter. CAMRAD II is being used to study the effects of various design features on whirl flutter, with special attention to areas where such features differ from existing aircraft, notably tiltrotors. Although the stability margins appear to be more than adequate, the concept requires significantly different analytical methods, principally including far more blade modes, than typically used for tiltrotors.

Description: The paper summarizes the state of the art in aeronautical drag reduction across the speed range for the conventional drag components of viscous drag, drag due to lift and wave drag. It also describes several emerging drag-reduction approaches that are either active or reactive/interactive and the drag reduction potentially available from synergistic combinations of advanced configuration aerodynamics, viscous drag-reduction approaches, revolutionary structural concepts and propulsion integration.

Description: In an effort to discover the causes for disagreement between previous two-dimensional (2-D) computations and nominally 2-D experiment for flow over the three-element McDonnell Douglas 30P-30N airfoil configuration at high lift, a combined experimental/CFD investigation is described. The experiment explores several different side-wall boundary layer control venting patterns, documents venting mass flow rates, and looks at corner surface flow patterns. The experimental angle of attack at maximum lift is found to be sensitive to the side-wall venting pattern: a particular pattern increases the angle of attack at maximum lift by at least 2 deg. A significant amount of spanwise pressure variation is present at angles of attack near maximum lift. A CFD study using three-dimensional (3-D) structured-grid computations, which includes the modeling of side-wall venting, is employed to investigate 3-D effects on the flow. Side-wall suction strength is found to affect the angle at which maximum lift is predicted. Maximum lift in the CFD is shown to be limited by the growth of an off-body corner flow vortex and consequent increase in spanwise pressure variation and decrease in circulation. The 3-D computations with and without wall venting predict similar trends to experiment at low angles of attack, but either stall too early or else overpredict lift levels near maximum lift by as much as 5%. Unstructured-grid computations demonstrate that mounting brackets lower the lift levels near maximum lift conditions.

Description: Circulation control wings are a type of pneumatic high-lift device that have been extensively researched as to their aerodynamic benefits. However, there has been little research into the possible airframe noise reduction benefits of a circulation control wing. The key element of noise is the jet noise associated with the jet sheet emitted from the blowing slot. This jet sheet is essentially a high aspect-ratio rectangular jet. A recent study on high aspect-ratio jet noise was performed on a nozzle with aspect-ratios ranging from 100 to 3,000. In addition to the acoustic data, fluid dynamic measurements were made as well. This paper uses the results of these two studies and attempts to develop a prediction scheme for high aspect-ratio jet noise

Description: The aerodynamic characteristics of a Circulation Control Wing (CCW) airfoil have been numerically investigated, and comparisons with experimental data have been made. The configuration chosen was a supercritical airfoil with a 30 degree dual-radius CCW flap. Steady and pulsed jet calculations were performed. It was found that the use of steady jets, even at very small mass flow rates, yielded a lift coefficient that is comparable or superior to conventional high-lift systems. The attached flow over the flap also gave rise to lower drag coefficients, and high L/D ratios. Pulsed jets with a 50% duty cycle were also studied. It was found that they were effective in generating lift at lower reduced mass flow rates compared to a steady jet, provided the pulse frequency was sufficiently high. This benefit was attributable to the fact that the momentum coefficient of the pulsed jet, during the portions of the cycle when the jet was on, was typically twice as much as that of a steady jet.

Description: This portion of the report documents the results of an experimental program, which focused on pulsed blowing from the trailing edge of a CCW. The main objective of this study was to assess whether pulsed blowing resulted in more, less, or the same amount of radiated noise to the farfield. Results show that a reduction in far-field noise of up to 5 dB is measured when pulse flow is compared to steady flow for an equivalent lift configuration. This reduction is in the spectral region associated with the trailing edge jet noise. This result is due to the unique advantage that pulsed flow has over steady flow. For a range of frequencies, more lift is experienced with the same mass flow as the steady case. Thus, for an equivalent lift and slot height, the pulsed system can operate at lower jet velocities, and hence lower jet noise. At low frequencies (below 1 kHz), the pulsed flow configuration generated more noise in the farfield. This is most likely due to the pulsing mechanism itself. Since the high pressure air feeding the pulsing mechanism was first passed through a high performance muffler, it is likely that this increase in not due to upstream valve noise. Most likely, the impulsive component of the air that periodically fills the plenum causes a broadband source that reaches the farfield. Although the benefit of a pulse trailing edge jet is evident from a mass flow usage and jet noise perspective, attention should be paid towards the design of a viable pulsing system. Future research program in this area should concentrate on the development of a "quiet" pulsing device.

Description: Circulation control wings are a type of pneumatic high-lift device that have been extensively researched as to their aerodynamic benefits. However, there has been little research into the possible airframe noise reduction benefits of a circulation control wing. The key element of noise is the jet noise associated with the jet sheet emitted from the blowing slot. High aspect-ratio jet acoustic results (aspect-ratios from 100 to 3,000) from a related study showed that the jet noise of this type of jet was proportional to the slot height to the 3/2 power and slot width to the 1/2 power. Fluid dynamic experiments were performed in the present study on the high aspect-ratio nozzle to gain understanding of the flow characteristics in an effort to relate the acoustic results to flow parameters. Single hot-wire experiments indicated that the jet exhaust from the high aspect-ratio nozzle was similar to a 2-d turbulent jet. Two-wire space-correlation measurements were performed to attempt to find a relationship between the slot height of the jet and the length-scale of the flow noise generating turbulence structure. The turbulent eddy convection velocity was also calculated, and was found to vary with the local centerline velocity, and also as a function of the frequency of the eddy.

Description: This Appendix documents the salient results from an effort to mitigate the so-called flap-edge noise generated at the split between a flap edge that is deployed and the undeployed flap. Utilizing a Coanda surface installed at the flap edge, steady blowing was used in an attempt to diminish the vortex strength resulting from the uneven lift distribution. The strength of this lifting vortex was augmented by steady blowing over the deployed flap. The test article for this study was the same 2D airfoil used in the steady blowing program reported earlier (also used in pulsed blowing tests, see Appendix G), however its trailing edge geometry was modified. An exact duplicate of the airfoil shape was made out of fiberglass with no flap, and in the clean configuration. It was attached to the existing airfoil to make an airfoil that has half of its flap deployed and half un-deployed. Figure 1 shows a schematic of the planform showing the two areas where steady blowing was introduced. The flap-edge blowing or the auxiliary blowing was in the direction normal to the freestream velocity vector. Slot heights for the blowing chambers were on the order of 0.0 14 inches.

Description: Six direct numerical simulations of turbulent time-evolving strained plane wakes have been examined to investigate the response of a wake to successive irrotational plane strains of opposite sign. The orientation of the applied strain field has been selected so that the flow is the time-developing analogue of a spatially developing wake evolving in the presence of either a favourable or an adverse streamwise pressure gradient. The magnitude of the applied strain rate a is constant in time t until the total strain e(sup at) reaches about four. At this point, a new simulation is begun with the sign of the applied strain being reversed (the original simulation is continued as well). When the total strain is reduced back to its original value of one, yet another simulation is begun with the sign of the strain being reversed again back to its original sign. This process is done for both initially "favourable" and initially "adverse" strains, providing simulations for each of these strain types from three different initial conditions. The evolution of the wake mean velocity deficit and width is found to be very similar for all the adversely strained cases, with both measures rapidly achieving exponential growth at the rate associated with the cross-stream expansive strain e(sup at). In the "favourably" strained cases, the wake widths approach a constant and the velocity deficits ultimately decay rapidly as e(sup -2at). Although all three of these cases do exhibit the same asymptotic exponential behaviour, the time required to achieve this is longer for the cases that have been previously adversely strained (by at approx. equals 1). These simulations confirm the generality of the conclusions drawn in Rogers (2002) regarding the response of plane wakes to strain. The evolution of strained wakes is not consistent with the predictions of classical self-similar analysis; a more general equilibrium similarity solution is required to describe the results. At least for the cases considered here, the wake Reynolds number and the ratio of the turbulent kinetic energy to the square of the wake mean velocity deficit are determined nearly entirely by the total strain. For these measures the order in which the strains are applied does not matter and the changes brought about by the strain are nearly reversible. The wake mean velocity deficit and width, on the other hand, differ by about a factor of three when the total strain returns to one, depending on whether the wake was first "favourably" or "adversely" strained. The strain history is important for predicting the evolution of these quantities.

Description: Probabilistic CFD design is needed because we are asked to do more with less. To cost effectively accomplish the design task, we need to formally quantify the effect of uncertainties (variables) in the design. Probabilistic design is one effective method to formally quantify the effect of uncertainties. Our objective is to establish a revolutionary new early design process, by developing non-deterministic physics-based probabilistic design tools, which will include all the life cycle processes. This work was concerned with the usefulness of parametric optimization method coupled with a Navier-Stokes analysis code for the aero-thermodynamic design of turbomachinery combustor liner. The interconnection between the CFD code and NESSUS codes will facilitate the coupling between the thermal profiles and structural design. We have developed new concepts for reducing the computational cost of unsteady, three-dimensional, compressible aerodynamic analyses for multistage turbomachinery flows. The flow was modeled by the three-dimensional Favre-Reynolds-averaged Navier-Stokes equations using the k-E turbulence closure, which was integrated using an implicit third-order upwind solver. The methodology developed in this work is expected to lead to the design optimization of turbomachinery blades.

Description: Circulation Control Wing (CCW) technology is a very effective way of achieving very high lift coefficients needed by aircraft during take-off and landing. This technology can also be used to directly control the flow field over the wing. Compared to a conventional high-lift system, a Circulation Control Wing (CCW) can generate the required values of lift coefficient C(sub L,max) during take-off/landing with fewer or no moving parts and much less complexity. Earlier designs of CCW configurations used airfoils with a large radius rounded trailing edge to maximize the lift benefit. However, these designs also produced very high drag. These high drag levels associated with the blunt, large radius trailing edge can be prohibitive under cruise conditions when Circulation Control is no longer necessary. To overcome this difficulty, an advanced CCW section, i.e., a circulation hinged flap was developed to replace the original rounded trailing edge CC airfoil. This concept developed by Englar is shown. The upper surface of the CCW flap is a large-radius arc surface, but the lower surface of the flap is flat. The flap could be deflected from 0 degrees to 90 degrees. When an aircraft takes-off or lands, the flap is deflected as in a conventional high lift system. Then this large radius on the upper surface produces a large jet turning angle, leading to high lift. When the aircraft is in cruise, the flap is retracted and a conventional sharp trailing edge shape results, greatly reducing the drag. This kind of flap does have some moving elements that increase the weight and complexity over an earlier CCW design. But overall, the hinged flap design still maintains most of the Circulation Control high lift advantages, while greatly reducing the drag in cruising condition associated with the rounded trailing edge CCW design. In the present work, an unsteady three-dimensional Navier-Stokes analysis procedure has been developed and applied to this advanced CCW configuration. The solver can be used in both a 2-D and a 3-D mode, and can thus model airfoils as well as finite wings. The jet slot location, slot height, and the flap angle can all be varied easily and individually in the grid generator and the flow solver. Steady jets, pulsed jets, the leading edge and trailing edge blowing can all be studied with this solver.

Description: Circulation control technology uses tangential blowing around a rounded trailing edge or a leading edge to change the force and moment characteristics of an aerodynamic body. This technology has been applied to circular cylinders, wings, helicopter rotors, and even to automobiles for improved aerodynamic performance. Only limited research has been conducted on the acoustics of this technology. Since wing flaps contribute to the environmental noise of an aircraft, an alternate blown high lift system without complex mechanical flaps could prove beneficial in reducing the noise of an approaching aircraft. Thus, in this study, a direct comparison of the acoustic characteristics of high lift systems employing a circulation control wing configuration and a conventional wing flapped configuration has been made. These results indicate that acoustically, a circulation control wing high lift system could be considerably more acceptable than a wing with conventional mechanical flaps.

Description: The U.S. National Aeronautics and Space Administration (NASA) Sounding Rockets and Balloon Programs conduct a total of 50 to 60 missions per year in support of the NASA scientific community. These missions support investigations sponsored by NASA's Offices of Space Science, Life and Microgravity Sciences & Applications, and Earth Science. The Goddard Space Flight Center has management and implementation responsibility for these programs. The NASA Sounding Rockets Program provides the science community with payload development support, environmental testing, launch vehicles, and launch operations from fixed and mobile launch ranges. Sounding rockets continue to provide a cost-effective way to make in situ observations from 50 to 1500 km in the near-earth environment and to uniquely cover the altitude regime between 50 km and 130 km above the Earth's surface. New technology efforts include GPS payload event triggering, tailored trajectories, new vehicle configuration development to expand current capabilities, and the feasibility assessment of an ultra high altitude sounding rocket vehicle. The NASA Balloon Program continues to make advancements and developments in its capabilities for support of the scientific ballooning community. The Long Duration Balloon (LDB) is capable of providing flight durations in excess of two weeks and has had many successful flights since its development. The NASA Balloon Program is currently engaged in the development of the Ultra Long Duration Balloon (ULDB), which will be capable of providing flight times up to 100-days. Additional development efforts are focusing on ultra high altitude balloons, station keeping techniques and planetary balloon technologies.

Description: An analysis of the optimal non-modal growth of perturbations in a boundary layer in the presence of a streamwise pressure gradient is presented. The analysis is based on PSE equations for an incompressible fluid. Examples with Falkner-Scan profiles indicate that a favorable pressure gradient decreases the non-modal growth, while an unfavorable pressure gradient leads to an increase of the amplification. It is suggested that the transient growth mechanism be utilized to choose optimal parameters of tripping elements on a low-pressure turbine (LPT) airfoil. As an example, a boundary layer flow with a streamwise pressure gradient corresponding to the pressure distribution over a LPT airfoil is considered. It is shown that there is an optimal spacing of the tripping elements and that the transient growth effect depends on the starting point.

Description: Computational Fluid Dynamics (CFD) analyses of axisymmetric circular-arc boattail nozzles operating off-design at transonic Mach numbers have been completed. These computations span the very difficult transonic flight regime with shock-induced separations and strong adverse pressure gradients. External afterbody and internal nozzle pressure distributions computed with the Wind code are compared with experimental data. A range of turbulence models were examined, including the Explicit Algebraic Stress model. Computations have been completed at freestream Mach numbers of 0.9 and 1.2, and nozzle pressure ratios (NPR) of 4 and 6. Calculations completed with variable time-stepping (steady-state) did not converge to a true steady-state solution. Calculations obtained using constant timestepping (timeaccurate) indicate less variations in flow properties compared with steady-state solutions. This failure to converge to a steady-state solution was the result of using variable time-stepping with large-scale separations present in the flow. Nevertheless, time-averaged boattail surface pressure coefficient and internal nozzle pressures show reasonable agreement with experimental data. The SST turbulence model demonstrates the best overall agreement with experimental data.

Description: Model-scale aeroacoustic tests of large civil transports point to the leading-edge slat as a dominant high-lift noise source in the low- to mid-frequencies during aircraft approach and landing. Using generic multi-element high-lift models, complementary experimental and numerical tests were carefully planned and executed at NASA in order to isolate slat noise sources and the underlying noise generation mechanisms. In this paper, a brief overview of the supporting computational effort undertaken at NASA Langley Research Center, is provided. Both tonal and broadband aspects of slat noise are discussed. Recent gains in predicting a slat s far-field acoustic noise, current shortcomings of numerical simulations, and other remaining open issues, are presented. Finally, an example of the ever-expanding role of computational simulations in noise reduction studies also is given.

Description: Trailing Ballute Aerocapture offers the potential to obtain orbit insertion around a planetary body at a fraction of the mass of traditional methods. This allows for lower costs for launch, faster flight times and additional mass available for science payloads. The technique involves an inflated ballute (balloon-parachute) that provides aerodynamic drag area for use in the atmosphere of a planetary body to provide for orbit insertion in a relatively benign heating environment. To account for atmospheric, navigation and other uncertainties, the ballute is oversized and detached once the desired velocity change (Delta V) has been achieved. Analysis and trades have been performed for the purpose of assessing the feasibility of the technique including aerophysics, material assessments, inflation system and deployment sequence and dynamics, configuration trades, ballute separation and trajectory analysis. Outlined is the technology development required for advancing the technique to a level that would allow it to be viable for use in space exploration missions.

Description: The wake behind a two-bladed model rotor in light climb was measured using particle image velocimetry, with particular emphasis on the development of the trailing vortex during the first revolution of the rotor. The distribution of vorticity was distinguished from the slightly elliptical swirl pattern. Peculiar dynamics within the void region may explain why the peak vorticity appeared to shift away from the center as the vortex aged, suggesting the onset of instability. The swirl and axial velocities (which reached 44 and 12 percent of the rotor-tip speed, respectively) were found to be asymmetric relative to the vortex center. In particular, the axial flow was composed of two concentrated zones moving in opposite directions. The radial distribution of the circulation rapidly increased in magnitude until reaching a point just beyond the core radius, after which the rate of growth decreased significantly. The core-radius circulation increased slightly with wake age, but the large-radius circulation appeared to remain relatively constant. The radial distributions of swirl velocity and vorticity exhibit self-similar behaviors, especially within the core. The diameter of the vortex core was initially about 10 percent of the rotor-blade chord, but more than doubled its size after one revolution of the rotor. According to vortex models that approximate the measured data, the core-radius circulation was about 79 percent of the large-radius circulation, and the large-radius circulation was about 67 percent of the maximum bound circulation on the rotor blade. On average, about 53 percent of the maximum bound circulation resides within the vortex core during the first revolution of the rotor.

Description: A generic spinning missile with dithering canards is used to demonstrate the utility of an overset structured grid approach for simulating the aerodynamics of rolling airframe missile systems. The approach is used to generate a modest aerodynamic database for the generic missile. The database is populated with solutions to the Euler and Navier-Stokes equations. It is used to evaluate grid resolution requirements for accurate prediction of instantaneous missile loads and the relative aerodynamic significance of angle-of-attack, canard pitching sequence, viscous effects, and roll-rate effects. A novel analytical method for inter- and extrapolation of database results is also given.

Description: The objective for this paper is to present the development of an optimization capability for the Cartesian inviscid-flow analysis package of Aftosmis et al. We evaluate and characterize the following modules within the new optimization framework: (1) A component-based geometry parameterization approach using a CAD solid representation and the CAPRI interface. (2) The use of Cartesian methods in the development Optimization techniques using a genetic algorithm. The discussion and investigations focus on several real world problems of the optimization process. We examine the architectural issues associated with the deployment of a CAD-based design approach in a heterogeneous parallel computing environment that contains both CAD workstations and dedicated compute nodes. In addition, we study the influence of noise on the performance of optimization techniques, and the overall efficiency of the optimization process for aerodynamic design of complex three-dimensional configurations. of automated optimization tools. rithm and a gradient-based algorithm.

Description: The gust response of a 2 D cascade is studied by solving the full nonlinear Euler equations employing higher order accurate spatial differencing and time stepping techniques. The solutions exhibit the exponential decay of the two circumferential mode orders of the cutoff blade passing frequency (BPF) tone and propagation of one circumferential mode order at 2BPF, as would be expected for the flow configuration considered. Two frequency excitations indicate that the interaction between the frequencies and the self interaction contribute to the amplitude of the propagating mode.

Description: A new high Reynolds number test capability for boundary layer ingesting inlets has been developed for the NASA Langley Research Center 0.3-Meter Transonic Cryogenic Tunnel. Using this new capability, an experimental investigation of four S-duct inlet configurations with large amounts of boundary layer ingestion (nominal boundary layer thickness of about 40% of inlet height) was conducted at realistic operating conditions (high subsonic Mach numbers and full-scale Reynolds numbers). The objectives of this investigation were to 1) develop a new high Reynolds number, boundary-layer ingesting inlet test capability, 2) evaluate the performance of several boundary layer ingesting S-duct inlets, 3) provide a database for CFD tool validation, and 4) provide a baseline inlet for future inlet flow-control studies. Tests were conducted at Mach numbers from 0.25 to 0.83, Reynolds numbers (based on duct exit diameter) from 5.1 million to a fullscale value of 13.9 million, and inlet mass-flow ratios from 0.39 to 1.58 depending on Mach number. Results of this investigation indicate that inlet pressure recovery generally decreased and inlet distortion generally increased with increasing Mach number. Except at low Mach numbers, increasing inlet mass-flow increased pressure recovery and increased distortion. Increasing the amount of boundary layer ingestion (by decreasing inlet throat height and increasing inlet throat width) or ingesting a boundary layer with a distorted profile decreased pressure recovery and increased distortion. Finally, increasing Reynolds number had almost no effect on inlet distortion but increased inlet recovery by about one-half percent at a Mach number near cruise.

Description: A thin pressure sensitive paint (PSP) coating can slightly modify the overall shape of a wind-tunnel model and produce surface roughness or smoothness that does not exist on the unpainted model. These undesirable changes in model geometry may alter flow over the model, and affect the pressure distribution and aerodynamic forces and moments on the model. This study quantifies the effects of PSP on three models in low-speed, transonic and supersonic flow regimes. At a 95% confidence level, the PSP effects on the integrated forces are insignificant for a slender arrow-wing-fuselage model and delta wing model with two different paints at Mach 0.2, 1.8, and 2.16 relative to the total balance accuracy limit. The data displayed a repeatability of 2.5 drag counts, while the balance accuracy limit was about 5.5 drag counts. At transonic speeds, the paint has a localized effect at high angles of attack and has a resolvable effect on the normal force, which is significant relative to the balance accuracy limit. For low speeds, the PSP coating has a localized effect on the pressure tap measurements, which leads to an appreciable decrease in the pressure tap reading. Moreover, the force and moment measurements had a poor precision, which precluded the ability to measure the PSP effect for this particular test.

Description: This paper reviews past work and presents new data to evaluate how changes in similarity parameters affect ice shapes and how closely scale values of the parameters should match reference values. Experimental ice shapes presented are from tests by various researchers in the NASA Glenn Icing Research Tunnel. The parameters reviewed are the modified inertia parameter (which determines the stagnation collection efficiency), accumulation parameter, freezing fraction, Reynolds number, and Weber number. It was demonstrated that a good match of scale and reference ice shapes could sometimes be achieved even when values of the modified inertia parameter did not match precisely. Consequently, there can be some flexibility in setting scale droplet size, which is the test condition determined from the modified inertia parameter. A recommended guideline is that the modified inertia parameter be chosen so that the scale stagnation collection efficiency is within 10 percent of the reference value. The scale accumulation parameter and freezing fraction should also be within 10 percent of their reference values. The Weber number based on droplet size and water properties appears to be a more important scaling parameter than one based on model size and air properties. Scale values of both the Reynolds and Weber numbers need to be in the range of 60 to 160 percent of the corresponding reference values. The effects of variations in other similarity parameters have yet to be established.