ALBERT

Description: This paper highlights some of the results and issues associated with estimating models to evaluate control law design methods and design criteria for advanced high performance aircraft. Experimental fighter aircraft such as the NASA High Alpha Research Vehicle (HARV) have the capability to maneuver at very high angles of attack where nonlinear aerodynamics often predominate. HARV is an experimental F/A-18, configured with thrust vectoring and conformal actuated nose strakes. Identifying closed-loop models for this type of aircraft can be made difficult by nonlinearities and high-order characteristics of the system. In this paper only lateral-directional axes are considered since the lateral-directional control law was specifically designed to produce classical airplane responses normally expected with low-order, rigid-body systems. Evaluation of the control design methodology was made using low-order equivalent systems determined from flight and simulation. This allowed comparison of the closed-loop rigid-body dynamics achieved in flight with that designed in simulation. In flight, the On Board Excitation System was used to apply optimal inputs to lateral stick and pedals at five angles of attack: 5, 20, 30, 45, and 60 degrees. Data analysis and closed-loop model identification were done using frequency domain maximum likelihood. The structure of the identified models was a linear state-space model reflecting classical 4th-order airplane dynamics. Input time delays associated with the high-order controller and aircraft system were accounted for in data preprocessing. A comparison of flight estimated models with small perturbation linear design models highlighted nonlinearities in the system and indicated that the estimated closed-loop rigid-body dynamics were sensitive to input amplitudes at 20 and 30 degrees angle of attack.

Description: This chapter is concerned with three-dimensional imaging of fluid flows. Although relatively young, this field of research has already yielded an enormous range of techniques. These vary widely in cost and complexity, with the cheapest light sheet systems being within the budgets of most laboratories, and the most expensive Magnetic Resonance Imaging systems available to a select few. Taking the view that the most likely systems to be developed are those using light sheets, the authors will relate their knowledge and experience of such systems. Other systems will be described briefly and references provided. Flows are inherently three-dimensional in structure; even those generated around nominally 2-D surface geometry. It is becoming increasingly apparent to scientists and engineers that the three-dimensionalities, both large and small scale, are important in terms of overall flow structure and species, momentum, and energy transport. Furthermore, we are accustomed to seeing the world in three dimensions, so it is natural that we should wish to view, measure and interpret flows in three-dimensions. Unfortunately, 3-D images do not lend themselves to convenient presentation on the printed page, and this task is one of the challenges facing us.

Description: The objectives of this effort were to determine (if possible) the best method: 1) for forcing the boundary layer to transition, 2) for assessing trip effectiveness, 3) for quantifying trip drag, 4) for testing at Reynolds numbers per foot from 5 million to maximum available rather than I to 5 million, and 5) for boundary layer state determination.

Description: Two wind tunnel tests during 1995 in the National Transonic Facility (NTF 070 and 073) served to define Reynolds number effects on longitudinal and lateral-directional stability and control. Testing was completed at both high lift and transonic conditions. The effect of Reynolds number on the total airplane configuration, horizontal and vertical tail effectiveness, forebody chine performance, rudder control and model aeroelastics was investigated. This paper will present pertinent stability and control results from these two test entries. Note that while model aeroelastic effects are examined in this presentation, no corrections for these effects have been made to the data.

Description: Buffeting is an aeroelastic phenomenon occurring at high angles of attack that plagues high performance aircraft, especially those with twin vertical tails. Previous wind-tunnel and flight tests were conducted to characterize the buffet loads on the vertical tails by measuring surface pressures, bending moments, and accelerations. Following these tests, buffeting responses were computed using the measured buffet pressures and compared to the measured buffeting responses. The calculated results did not match the measured data because the assumed spatial correlation of the buffet pressures was not correct. A better understanding of the partial (spatial) correlation of the differential buffet pressures on the tail was necessary to improve the buffeting predictions. Several wind-tunnel investigations were conducted for this purpose. When compared, the results of these tests show that the partial correlation scales with flight conditions. One of the remaining questions is whether the wind-tunnel data is consistent with flight data. Presented herein, cross-spectra and coherence functions calculated from pressures that were measured on the High Alpha Research Vehicle indicate that the partial correlation of the buffet pressures in flight agrees with the partial correlation observed in the wind tunnel.

Description: The objective was to experimentally evaluate the longitudinal and lateral-directional stability and control characteristics of the Reference H configuration at supersonic and transonic speeds. A series of conventional and alternate control devices were also evaluated at supersonic and transonic speeds. A database on the conventional and alternate control devices was to be created for use in the HSR program.

Description: The stability and control issues in high speed aerodynamics of most significance for the development of a viable HSCT are identified, and the status of the Ref. H configuration with respect to these issues is discussed. The interdependence between aerodynamic requirements and assumptions about airplane system functions such as Envelope Protection and Integrated Flight/Propulsion Control is highlighted. The conclusions presented draw on results from the Ref. H Assessment and Alternate Control Concepts Assessment performed under Configuration Aerodynamics Subtask 5 during 1995.

Description: NASA Langley Research Center (LaRC) has field tested a 2.0 gm, 100 Hertz, pulsed coherent lidar to detect and characterize wake vortices and to measure atmospheric winds and turbulence. The quantification of aircraft wake-vortex hazards is being addressed by the Wake Vortex Lidar (WVL) Project as part of Aircraft Vortex Spacing System (AVOSS), which is under the Reduced Spacing Operations Element of the Terminal Area Productivity (TAP) Program. These hazards currently set the minimum, fixed separation distance between two aircraft and affect the number of takeoff and landing operations on a single runway under Instrument Meteorological Conditions (IMC). The AVOSS concept seeks to safely reduce aircraft separation distances, when weather conditions permit, to increase the operational capacity of major airports. The current NASA wake-vortex research efforts focus on developing and validating wake vortex encounter models, wake decay and advection models, and wake sensing technologies. These technologies will be incorporated into an automated AVOSS that can properly select safe separation distances for different weather conditions, based on the aircraft pair and predicted/measured vortex behavior. The sensor subsystem efforts focus on developing and validating wake sensing technologies. The lidar system has been field-tested to provide real-time wake vortex trajectory and strength data to AVOSS for wake prediction verification. Wake vortices, atmospheric winds, and turbulence products have been generated from processing the lidar data collected during deployments to Norfolk (ORF), John F. Kennedy (JFK), and Dallas/Fort Worth (DFW) International Airports.

Description: Several analytical and experimental studies clearly demonstrate that piezoelectric materials (piezoelectrics) can be used as actuators to actively control vibratory response, including aeroelastic response. However, two important issues in using piezoelectrics as actuators for active control are: 1) the potentially large amount of power required to operate the actuators, and 2) the complexities involved with active control (added hardware, control law design, and implementation). Active or passive damping augmentation using shunted piezoelectrics may provide a viable alternative. This approach requires only simple electrical circuitry and very little or no electrical power. The current study examines the feasibility of using shunted piezoelectrics to reduce aeroelastic response using a typical-section representation of a wing and piezoelectrics shunted with a parallel resistor and inductor. The aeroelastic analysis shows that shunted piezoelectrics can effectively reduce aeroelastic response below flutter and may provide a simple, low-power method of subcritical aeroelastic control.

Description: Subgrid analysis of a transitional temporal mixing layer with evaporating droplets has been performed using three sets of results from a Direct Numerical Simulation (DNS) database, with Reynolds numbers (based on initial vorticity thickness) as large as 600 and with droplet mass loadings as large as 0.5. In the DNS, the gas phase is computed using a Eulerian formulation, with Lagrangian droplet tracking. The Large Eddy Simulation (LES) equations corresponding to the DNS are first derived, and key assumptions in deriving them are first confirmed by computing the terms using the DNS database. Since LES of this flow requires the computation of unfiltered gas-phase variables at droplet locations from filtered gas-phase variables at the grid points, it is proposed to model these by assuming the gas-phase variables to be the sum of the filtered variables and a correction based on the filtered standard deviation; this correction is then computed from the Subgrid Scale (SGS) standard deviation. This model predicts the unfiltered variables at droplet locations considerably better than simply interpolating the filtered variables. Three methods are investigated for modeling the SGS standard deviation: the Smagorinsky approach, the Gradient model and the Scale-Similarity formulation. When the proportionality constant inherent in the SGS models is properly calculated, the Gradient and Scale-Similarity methods give results in excellent agreement with the DNS.

Description: The F/A-18 Active Aeroelastic Wing research aircraft will demonstrate technologies related to aeroservoelastic effects such as wing twist and load minimization. This program presents several challenges for control design that are often not considered for traditional aircraft. This paper presents a control design based on H(sub infinity) synthesis that simultaneously considers the multiple objectives associated with handling qualities, actuator limitations, and loads. A point design is presented to demonstrate a controller and the resulting closed-loop properties.

Description: A three-dimensional Navier Stokes code has been used to compute the heat transfer coefficient on a film-cooled, rotating turbine blade. The blade chosen is the ACE rotor with five rows containing 93 film cooling holes covering the entire span. This is the only film-cooled rotating blade over which experimental data is available for comparison. Over 2.278 million grid points are used to compute the flow over the blade including the tip clearance region. using Wilcox's k-omega model, Coakley's q-omega model, and the zero-equation Baldwin-Lomax (B-L) model. A reasonably good comparison with the experimental data is obtained on the suction surface for all the turbulence models. At the leading edge, the B-L model yields a better comparison than tile two-equation models. On the pressure surface however the comparison between the experimental data and the prediction from the k-omega model is much better than from the other two models. Overall, the k-omega model provides the best comparison with the experimental data. However, the two-equation models require at least 40% more computational resources than the B-L model.

Description: We report measurements of the depression of the superfluid transition temperature by a heat current (1 less than or = Q less than or = 100 microW/sq cm) along the lambda-line (SVP less than or = P less than or = 21.6 bar). At P = 21.6 bar, measurements were also performed in a reduced gravity (0.2g). Experimental results show that the pressure dependence of the depression and the gravity effect on the measurements are small, in qualitative agreement with theoretical predictions. Keywords: superfluid helium; Lambda transition; heat current

Description: The differences between subcritical liquid drop and supercritical fluid drop behavior are discussed. Under subcritical, evaporative high emission rate conditions, a film layer is present in the inner part of the drop surface which contributes to the unique determination of the boundary conditions; it is this film layer which contributes to the solution's convective-diffusive character. In contrast, under supercritical condition as the boundary conditions contain a degree of arbitrariness due to the absence of a surface, and the solution has then a purely diffusive character. Results from simulations of a free fluid drop under no-gravity conditions are compared to microgravity experimental data from suspended, large drop experiments at high, low and intermediary temperatures and in a range of pressures encompassing the sub-and supercritical regime. Despite the difference between the conditions of the simulations and experiments (suspension vs. free floating), the time rate of variation of the drop diameter square is remarkably well predicted in the linear curve regime. The drop diameter is determined in the simulations from the location of the maximum density gradient, and agrees well with the data. It is also shown that the classical calculation of the Lewis number gives qualitatively erroneous results at supercritical conditions, but that an effective Lewis number previously defined gives qualitatively correct estimates of the length scales for heat and mass transfer at all pressures.

Description: We report here, preliminary data from an experiment studying flow of superfluid helium through a slit orifice (of sub-micron width) very close to T(sub lambda). Critical supercurrent (I(sub c)) data is obtained from a step function drive to the diaphragm in a Helmholtz resonator cell. The superfluid density (rho(sub s)) data can be obtained from the resonant frequency of the Helmholtz oscillator, as determined by transfer function of the resonator or from the free ringing after the step function excitation. Preliminary data shows that I(sub c) is proportional to (rho(sub s))(exp 1.27) and rho(sub s)) is proportional to tau(exp 0.73), where tau is the reduced temperature. However, the magnitude of I(sub c) is much larger than expected, indicating a possible parallel flow path. Further investigations are in progress. Keywords: superfluid; hydrodynamics; critical exponent

Description: The requirements for increased speed and productivity for tiltrotors has spawned several investigations associated with proprotor aeroelastic stability augmentation and aerodynamic performance enhancements. Included among these investigations is a focus on passive aeroelastic tailoring concepts which exploit the anisotropic capabilities of fiber composite materials. Researchers at Langley Research Center and Bell Helicopter have devoted considerable effort to assess the potential for using these materials to obtain aeroelastic responses which are beneficial to the important stability and performance considerations of tiltrotors. Both experimental and analytical studies have been completed to examine aeroelastic tailoring concepts for the tiltrotor, applied either to the wing or to the rotor blades. This paper reviews some of the results obtained in these aeroelastic tailoring investigations and discusses the relative merits associated with these approaches.

Description: The F/A-18 Active Aeroelastic Wing research aircraft will demonstrate technologies related to aeroservoelastic effects such as wing twist and load minimization. This program presents several challenges for control design that are often not considered for traditional aircraft. This paper presents a control design based on H-infinity synthesis that simultaneously considers the multiple objectives associated with handling qualities, actuator limitations, and loads. A point design is presented to demonstrate a controller and the resulting closed-loop properties.

Description: This paper describes and discusses the textbook, Fundamentals of Computational Fluid Dynamics by Lomax, Pulliam, and Zingg, which is intended for a graduate level first course in computational fluid dynamics. This textbook emphasizes fundamental concepts in developing, analyzing, and understanding numerical methods for the partial differential equations governing the physics of fluid flow. Its underlying philosophy is that the theory of linear algebra and the attendant eigenanalysis of linear systems provides a mathematical framework to describe and unify most numerical methods in common use in the field of fluid dynamics. Two linear model equations, the linear convection and diffusion equations, are used to illustrate concepts throughout. Emphasis is on the semi-discrete approach, in which the governing partial differential equations (PDE's) are reduced to systems of ordinary differential equations (ODE's) through a discretization of the spatial derivatives. The ordinary differential equations are then reduced to ordinary difference equations (O(Delta)E's) using a time-marching method. This methodology, using the progression from PDE through ODE's to O(Delta)E's, together with the use of the eigensystems of tridiagonal matrices and the theory of O(Delta)E's, gives the book its distinctiveness and provides a sound basis for a deep understanding of fundamental concepts in computational fluid dynamics.

Description: Buoyant flows play an important role in various technological and environmental issues. For example, dispersal of pollutants, smoke, or volcano exhaust in the atmosphere, vertical motion of air, formation of clouds and other weather systems, and flows in cooling towers and fires are all determined primarily by buoyancy effects. The buoyancy force in such flows can originate from either a heat source or due to different densities between a fluid and its surroundings. Whatever the cause, the flow can be understood by studying the effects of the tight coupling between the thermal and the velocity fields since density differences can be characterized as temperature differences.

Description: Our objective is to predict droplet size distributions created by fuel injector nozzles in Jet turbines. These results will be used to determine the initial conditions for numerical simulations of the combustion process in gas turbine combustors. To predict the droplet size distribution, we are currently constructing a numerical model to understand the instability and breakup of thin conical liquid sheets. This geometry serves as a simplified model of the liquid jet emerging from a real nozzle. The physics of this process is difficult to study experimentally as the time and length scales are very short. From existing photographic data, it does seem clear that three-dimensional effects such as the formation of streamwise ligaments and the pulling back of the sheet at its edges under the action of surface tension are important.

Description: A detailed investigation of the flow physics occurring on the suction side of a simulated Low Pressure Turbine (LPT) blade was performed. A contoured upper wall was designed to simulate the y pressure distribution of an actual LPT blade onto a flat plate. The experiments were carried out at Reynolds numbers of 100,000 and 250,000 with three levels of freestream turbulence. The main emphasis in this paper is placed on flow field surveys performed at a y Reynolds number of 100,000 with levels of freestream turbulence ranging from 0.8% to 3%. Smoke-wire flow visualization data was used to confirm that the boundary layer was separated and formed a bubble. The transition process over the separated flow region is observed to be similar to a laminar free shear layer flow with the formation of a large coherent eddy structure. For each condition, the locations defining the separation bubble were determined by careful examination of pressure and mean velocity profile data. Transition onset location and length determined from intermittency profiles decrease as freestream turbulence levels increase. Additionally, the length and height of the laminar separation bubbles were observed to be inversely proportional to the levels of freestream turbulence.

Description: A distinct boundary between turbulent and non-turbulent regions in a fluid of otherwise constant properties is found in many laboratory and engineering turbulent flows, including jets, mixing layers, boundary layers and wakes. Generally, the flow has mean shear in at least one direction within t he turbulent zone, but the non-turbulent zones have no shear (adjacent laminar shear is a different case, e.g. transition in a boundary layer). There may be purely passive differences between the turbulent and non-turbulent zones, e.g. small variations in temperature or scalar concentration, for which turbulent mixing is an important issue. The boundary has several major characteristics of interest for the present study. Firstly, the boundary advances into the non-turbulent fluid, or in other words, nonturbulent fluid is entrained. Secondly, the change in turbulence properties across the boundary is remarkably abrupt; strong turbulent motions come close to the nonturbulent fluid, promoting entrainment. Thirdly, the boundary is irregular with a continually changing convoluted shape, which produces statistical intermittency. Its shape is contorted at all scales of the turbulent motion.

Description: Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This report demonstrates that a good compromise is possible with the Advected Grid Explicit (AGE) method. AGE has proven to be both efficient and accurate for simulating turbulent free-shear flows, including planar mixing layers and planar jets. Its efficiency results from its localized fully explicit finite difference formulation (Bisset 1998a,b) that is very straightforward to compute, outweighing the need for a fairly small timestep. Also, most of the successful simulations were slightly under-resolved, and therefore they were, in effect, large-eddy simulations (LES) without a sub-grid-scale (SGS) model, rather than direct numerical simulations (DNS). The principle is that the role of the smallest scales of turbulent motion (when the Reynolds number is not too low) is to dissipate turbulent energy, and therefore they do not have to be simulated when the numerical method is inherently dissipative at its resolution limits. Such simulations are termed 'auto-LES' (LES with automatic SGS modeling) in this report.

Description: Wavelets present a method for signal processing that may be useful for analyzing responses of dynamical systems. This paper describes several wavelet-based tools that have been developed to improve the efficiency of flight flutter testing. One of the tools uses correlation filtering to identify properties of several modes throughout a flight test for envelope expansion. Another tool uses features in time-frequency representations of responses to characterize nonlinearities in the system dynamics. A third tool uses modulus and phase information from a wavelet transform to estimate modal parameters that can be used to update a linear model and reduce conservatism in robust stability margins.

Description: The benchmark active controls technology and wind tunnel test program at NASA Langley Research Center was started with the objective to investigate the nonlinear, unsteady aerodynamics and active flutter suppression of wings in transonic flow. The paper will present the flutter suppression control law design process, numerical nonlinear simulation and wind tunnel test results for the NACA 0012 benchmark active control wing model. The flutter suppression control law design processes using (1) classical, (2) linear quadratic Gaussian (LQG), and (3) minimax techniques are described. A unified general formulation and solution for the LQG and minimax approaches, based on the steady state differential game theory is presented. Design considerations for improving the control law robustness and digital implementation are outlined. It was shown that simple control laws when properly designed based on physical principles, can suppress flutter with limited control power even in the presence of transonic shocks and flow separation. In wind tunnel tests in air and heavy gas medium, the closed-loop flutter dynamic pressure was increased to the tunnel upper limit of 200 psf The control law robustness and performance predictions were verified in highly nonlinear flow conditions, gain and phase perturbations, and spoiler deployment. A non-design plunge instability condition was also successfully suppressed.

Description: A common problem occurs when refrigerant fluids wick inside the insulation of thermocouple wires through a compressor's casing feedthrough and then leak into the adjacent disconnect box outside the casing. Leaking fluids create an unfavorable situation inside the disconnect box and may contaminate the fluids. To address this problem, NASA Lewis Research Center s Manufacturing Engineering Division developed a customized hermetic feedthrough for a bank of Worthington compressors. In these compressors, bearing temperatures are measured by internal thermocouples embedded in bearings located inside the compressor casings. The thermocouple wires need to be routed outside the casing and read at another location. These wires are short and are terminated to a disconnect strip inside the casing. The bearings operate at about 170 F, but because the casing is filled with R12 refrigerant oil, the casing has a maximum temperature of about 100 F. The operating conditions of these compressors permit the use of an epoxy that is compatible with the R12 fluid. The desired finished product is a stainless steel tube that has been filled solid with epoxy after thermocouple wires bonded and sealed by epoxy have been inserted through its length. Shrink tubing extends from both ends of the tube. The process that was developed to isolate the thermocouple wires from the R12 fluid follows. For this application, use an 8-in.-long piece of 0.500-in. 304 stainless steel tube with six pairs of 24-gauge stranded, PTFE-insulated (polytetrafluoroethylene) type "T" thermocouple wires for each feedthrough. Use shrink tubing to strain relief the insulated wires at their exit from the stainless steel tube. Cut the wire to length and identify the location of the stainless steel tube sleeve with masking tape. Then, remove the outer insulation from a 2-in. section of wire that will be inside the tube, and carefully strip to bare wire a 1-in. section in the middle of the section with the outer insulation removed. For an effective seal, the epoxy must penetrate between the strands when stranded conductors are used. Make the seal with epoxy bond on the bare wire. The bare wire must be encapsulated with a thin layer of the epoxy that leaves only a very low profile. These encapsulated wires must cure before the assembly can be continued. Then, inspect the cured wires for complete encapsulation before going to the next step. Insert the wires in the stainless steel tube and orient them so that the epoxied stripped sections are staggered within the tube; then, apply shrink tubing to one end of the cleaned wires, positioning it inside the edge of the tube. The small gaps between the wires on the other end will be used to inject the epoxy into the tube. Let the epoxy cure inside the tube, free of any voids. Then, continue to fill the tube until the entire 8-in. length is nearly filled, allowing room for the other strain-relieving shrink tubing. Since this first design, the process has been adjusted to fit many needs and situations. Customized feedthroughs have been assembled from various wire types, wire gauges, and/or stainless steel tube passages. The fittings selected to mount these feedthroughs allow their use in other areas, such as pressure or vacuum systems.

Description: High-Speed Civil Transport (HSCT) engine combustor liners are subjected to complex thermal environments and have to endure these for thousands of hours with assured reliability. In the past, several deterministic analyses have been performed, including detailed heat transfer analyses to obtain thermal profiles and deterministic stress analyses to identify critical locations of high stresses. Actual rig tests also have been performed for segments by simulating these loading situations as closely as possible. However, it is well known that many uncertainties exist in loading (primarily thermal loads due to heat transfer), boundary conditions (end fixity unknowns), and material properties (moduli, thermal-expansion coefficients, and conductivities). The present in-house effort at the NASA Lewis Research Center is directed toward accounting for these in a formal way to assess the performance of liner components under complex and uncertain loading conditions as well as subject to other geometry- and material-related uncertainties.

Description: This paper presents the results of a numerical study of shock/shock interactions that include both the Edney type IV and type III interactions, with emphasis on the type IV interactions. Computations are made using the direct simulation Monte Carlo (DSMC) method of Bird for Mach 10 air flow, as produced in the ONERA R5Ch low-density wind tunnel. The simulations include the flow about a shock generator which creates a relatively weak oblique shock that impinges on a much stronger cylinder bow shock. The sensitivity and characteristics of the interactions are examined by varying the horizontal distance separating the shock generator leading edge and cylinder. Results of the simulation for one separation distance are compared with wind tunnel measurements. Comparisons are made for surface heating and pressure and for flow-field values of density and rotational temperatures, as obtained with the Dual-line Coherent Anti-Stokes Scattering (DL-CARS) technique. The comparisons between experiment and calculation yield a consistent description of the shock interaction features and a consistent description of the surface heating and pressure distributions, with the exception of the peak values-the computed values being greater than the measured values.

Description: Problems concerning laminar-turbulent transition are addressed by solving a series of initial value problems. Solutions to the temporal, initial-value problem .with an inhomogeneous forcing term imposed upon the flow are sought. It is shown that: (1) A transient disturbance lying located outside of the boundary layer can lead to the growth of an unstable Tollmein-Schlicting wave; (2) A resonance with the continuous spectrum may provide a mechanism for bypass transition; and (3) The continuum modes of a disturbance feed directly into the Tollmein-Schlicting wave downstream through non-parallel effects.

Description: At the NASA Lewis Research Center, a comprehensive assessment was made of the predictive capability of the average passage flow model as applied to multistage axial-flow compressors. This model, which describes the time-averaged flow field within a typical passage of a blade row embedded in a multistage configuration, is being widely used throughout U.S. aircraft industry as an integral part of their design systems. Rotor flow-angle deviation. In this work, detailed data taken within a four and one-half stage large low-speed compressor were used to assess the weaknesses and strengths of the predictive capabilities of the average passage flow model. The low-speed compressor blading is of modern design and employs stator end-bends. Measurements were made with slow- and high response instrumentation. The high-response measurements revealed the velocity components of both the rotor and stator wakes. From the measured wake profiles, we found that the flow exiting the rotors deviated from the rotor exit metal angle to a lesser degree than was predicted by the average passage flow model. This was found to be due to blade boundary layer transition, which recently has been shown to exist on multistage axial compressor rotor and stator blades, but was not accounted for in the average passage model. Consequently, a model that mimics the effects of blade boundary layer transition, Shih k-epsilon model, was incorporated into the average passage model. Simulations that incorporated this transition model showed a dramatic improvement in agreement with data. The altered model thus improved predictive capability for multistage axial-flow compressors, and this was verified by detailed experimental measurement.

Description: A thorough understanding of dynamic interactions between inlets and compressors is extremely important to the design and development of propulsion control systems, particularly for supersonic aircraft such as the High-Speed Civil Transport (HSCT). Computational fluid dynamics (CFD) codes are routinely used to analyze individual propulsion components. By coupling the appropriate CFD component codes, it is possible to investigate inlet-compressor interactions. The objectives of this work were to gain a better understanding of inlet-compressor interaction physics, formulate a more realistic compressor-face boundary condition for time-accurate CFD simulations of inlets, and to take a first step toward the CFD simulation of an entire engine by coupling multidimensional component codes. This work was conducted at the NASA Lewis Research Center by a team of civil servants and support service contractors as part of the High Performance Computing and Communications Program (HPCCP).

Description: Electronic time-average holograms are convenient for comparing the measured vibration modes of fan blades with those calculated by finite-element models. At the NASA Lewis Research Center, neural networks recently were trained to perform what had been a simple visual comparison of the predictions of the design models with the measurements. Finite-element models were used to train neural networks to recognize damage and strain information encoded in subtle changes in the time-average patterns of cantilevers. But the design-grade finite element models were unable to train the neural networks to detect damage in complex blade shapes. The design-model-generated patterns simply did not agree well enough with the measured patterns. Instead, hybrid-training records, with measured time-average patterns as the input and model-generated strain information as the output, were used to effect successful training.

Description: NASA Lewis Research Center's Advanced Controls and Dynamics Technology Branch is investigating active control strategies to mitigate or eliminate the combustion instabilities prevalent in lean-burning, low-emission combustors. These instabilities result from coupling between the heat-release mechanisms of the burning process and the acoustic flow field of the combustor. Control design and implementation require a simulation capability that is both fast and accurate. It must capture the essential physics of the system, yet be as simple as possible. A quasi-one-dimensional, computational fluid dynamics (CFD) based simulation has been developed which may meet these requirements. The Euler equations of mass, momentum, and energy have been used, along with a single reactive species transport equation to simulate coupled thermoacoustic oscillations. A very simple numerical integration scheme was chosen to reduce computing time. Robust boundary condition procedures were incorporated to simulate various flow conditions (e.g., valves, open ends, and choked inflow) as well as to accommodate flow reversals that may arise during large flow-field oscillations. The accompanying figure shows a sample simulation result. A combustor with an open inlet, a choked outlet, and a large constriction approximately two thirds of the way down the length is shown. The middle plot shows normalized, time-averaged distributions of the relevant flow quantities, and the bottom plot illustrates the acoustic mode shape of the resulting thermoacoustic oscillation. For this simulation, the limit cycle peak-to-peak pressure fluctuations were 13 percent of the mean. The simulation used 100 numerical cells. The total normalized simulation time was 50 units (approximately 15 oscillations), which took 26 sec on a Sun Ultra2.

Description: Pulse Phase Thermography (PPT) has been reported as a novel powerful technique of the thermal NDE. It employs application of the Discrete Fourier Transform (DFT) to thermal images obtained following flash heating of the front surface of a specimen. The computed phasegrams are excellent for defect visualization in a wide range of materials. This is in part due to their low sensitivity to uneven heating. This work is an attempt to analyze advantages and limitations of PPT. Results of application of the DFT to simulated temperature decays are presented. The temperature evolution on a surface has been simulated based on an analytical solution of the one-dimensional heat diffusion problem. A more sophisticated study has been done for different sizes of defects using numerical solution of the three-dimensional mathematical model. Capabilities of PPT for in-depth scanning and for monitoring of the material loss are discussed. The recommendations for the practical application of the PPT are presented. Experimental results obtained following these recommendations are reported.

Description: Inertial and gravitational forces determine a particle's motion in a turbulent flow field. Gravity plays the dominant role in this motion by pulling the particles through adjacent regions of fluid turbulence. To better understand and model how a particle's inertia effects its displacement, one must examine the dispersion in a turbulent flow in the absence of gravity. In this paper, we present the particle experiments planned for NASA's KC-135 Reduced-Gravity Aircraft, which generates microgravity conditions for about 20 seconds. We also predict the particle behavior using simulation and ground-based experiments. We will release particles with Stokes numbers of 0.1, 1, and 10 into an enclosed tank of near-isotropic, stationary, and homogenous turbulence. These particle Stoke numbers cover a broad range of flow regimes of interest. Two opposed grids oscillating back and forth generate the turbulent field in the tank with a range of turbulence scales that covers about three orders of magnitude and with turbulence intensities of about ten times the mean velocity. The motion of the particles will be tracked using a stereo image velocimetry technique.

Description: Experimental results are presented for drops and bubbles levitated in a liquid host, with particular attention given to the effect of shape oscillations and capillary waves on the local flow fields. Some preliminary results are also presented on the use of streaming flows for the control of evaporation rate and rotation of electrostatically levitated droplets in 1 g. The results demonstrate the potential for the technological application of acoustic methods to active control of forced convection in microgravity.

Description: We have developed a new instrument that can measure fast transient birefringence and polymer chain orientation angle in complex fluids. The instrument uses a dual-crystal transverse electro-optic modulator with the second crystal's modulation voltage applied 180 deg out of phase from that of the first crystal. In this manner, the second crystal compensates for the intrinsic static birefringence of the first crystal, and it doubles the modulation depth. By incorporating a transverse electro-optic modulator with two lithium-niobate (LiNbO3) crystals oriented orthogonal to each other with a custom-designed optical system, we have produced a very small robust instrument capable of fast transient retardation measurements. By measuring the sample thickness or optical path length through the sample, we can calculate the transient birefringence. This system can also measure dichroism. We have compared the calibration results and retardation and orientation angle measurements of this instrument with those of a photoelastic modulator (PEM) based system using a quarter wave plate and a high-precision 1/16-wave plate to simulate a birefringent sample. Transient birefringence measurements on the order of 10(exp -9) can be measured using either modulator.

Description: The bulk behavior of dispersed, fluidized, or undispersed stationary granular systems cannot be fully understood in terms of adhesive/cohesive properties without understanding the role of electrostatic forces acting at the level of the grains themselves. When grains adhere to a surface, or come in contact with one another in a stationary bulk mass, it is difficult to measure the forces acting on the grains, and the forces themselves that induced the cohesion and adhesion are changed. Even if a single grain were to be scrutinized in the laboratory, it might be difficult, perhaps impossible, to define the distribution and character of surface charging and the three-dimensional relationship that charges (electrons, holes) have to one another. The hypothesis that we propose to test in microgravity (for dielectric materials) is that adhesion and cohesion of granular matter are mediated primarily by dipole forces that do not require the presence of a net charge; in fact, nominally electrically neutral materials should express adhesive and cohesive behavior when the neutrality results from a balance of positive and negative charge carriers. Moreover, the use of net charge alone as a measure of the electrical nature of grain-to-grain relationships within a granular mass may be misleading. We believe that the dipole forces arise from the presence of randomly-distributed positive and negative fixed charge carriers on grains that give rise to a resultant dipole moment. These dipole forces have long-range attraction. Random charges are created whenever there is triboelectrical activity of a granular mass, that is, whenever the grains experience contact/separation sequences or friction.

Description: The present work is aimed at the experimental studies and numerical modeling of the bubble growth mechanisms of a single bubble attached to a heating surface and of a bubble sliding along an inclined heated plate. Single artificial cavity of 10 microns in diameter was made on the polished Silicon wafer which was electrically heated at the back side in order to control the surface nucleation superheat. Experiments with a sliding bubble were conducted at different inclination angles of the downward facing heated surface for the purpose of studying the effect of magnitude of components of gravity acting parallel to and normal to the heat transfer surface. Information on the bubble shape and size, the bubble induced liquid velocities as well as the surface temperature were obtained using the high speed imaging and hydrogen bubble techniques. Analytical/numerical models were developed to describe the heat transfer through the micro-macro layer underneath and around a bubble formed at a nucleation site. In the micro layer model the capillary and disjoining pressures were included. Evolution of the bubble-liquid interface along with induced liquid motion was modeled. As a follow-up to the studies at normal gravity, experiments are being conducted in the KC-135 aircraft to understand the bubble growth/detachment under low gravity conditions. Experiments have been defined to be performed under long duration of microgravity conditions in the space shuttle. The experiment in the space shuttle will provide bubble growth and detachment data at microgravity and will lead to validation of the nucleate boiling heat transfer model developed from the preceding studies conducted at normal and low gravity (KC-135) conditions.

Description: In this work, we study steady and oscillatory thermocapillary and natural convective flows generated by a bubble on a heated solid surface. The interaction between gas and vapor bubbles with the surrounding fluid is of interest for both space and ground-based processing. A combined numerical-experimental approach is adopted here. The temperature field is visualized using Mach-Zehnder and/or Wollaston Prism Interferometry and the flow field is observed by a laser sheet flow visualization technique. A finite element numerical model is developed which solves the transient two-dimensional continuity, momentum, and energy equations and includes the effects of temperature-dependent surface tension and bubble surface deformation. Below the critical Marangoni number, the steady state low-g and 1-g temperature and velocity fields predicted by the finite element model are in excellent agreement with both the visualization experiments in our laboratory and recently published experimental results in the literature. Above the critical Marangoni number, the model predicts an oscillatory flow which is also closely confirmed by experiments. It is shown that the dynamics of the oscillatory flow are directly controlled by the thermal and hydrodynamic interactions brought about by combined natural and thermocapillary convection. Therefore, as numerical simulations show, there are considerable differences between the 1-g and low-g temperature and flow fields at both low and high Marangoni numbers. This has serious implications for both materials processing and fluid management in space.

Description: The WIND code is a general-purpose, structured, multizone, compressible flow solver that can be used to analyze steady or unsteady flow for a wide range of geometric configurations and over a wide range of flow conditions. WIND is the latest product of the NPARC Alliance, a formal partnership between the NASA Lewis Research Center and the Air Force Arnold Engineering Development Center (AEDC). WIND Version 1.0 was released in February 1998, and Version 2.0 will be released in February 1999. The WIND code represents a merger of the capabilities of three existing computational fluid dynamics codes--NPARC (the original NPARC Alliance flow solver), NXAIR (an Air Force code used primarily for unsteady store separation problems), and NASTD (the primary flow solver at McDonnell Douglas, now part of Boeing).

Description: About half of all scientific and commercial spacecraft spin during some part of their mission. Although spinning has many benefits (increasing stability, controlling the location of liquid propellants, and distributing solar heat loads), it also creates problems because a precession (or wobble) motion is unavoidable. For modern spacecraft, by far the largest source of wobble is energy dissipation in the liquid of partially filled tanks. The liquid s energy dissipation cannot, however, be quantified adequately by any ground test. Current analytical models are also inadequate because fundamental data on fluid motion in low gravity are needed to validate them. Consequently, spacecraft attitude-control systems are designed and operated very conservatively. Nonetheless, spacecraft often still perform poorly in orbit, and some have been lost because of a rapid increase of the wobble rate. The Liquid Motion Experiment (LME) was designed to provide spacecraft designers accurate data on the wobble dynamics of spacecraft that contain large quantities of mobile liquids. LME, which was flown on the space shuttle mission STS-84, was built under contract to the NASA Lewis Research Center by the Southwest Research Institute of San Antonio, Texas. Major accomplishments for 1998 include reduction of the flight data and publication of the experimental results. LME was essentially a spin table that created a realistic nutation motion of scale-model tanks containing liquid. Two spherical and two cylindrical transparent tanks were tested simultaneously, and three sets of such tanks were employed to vary liquid viscosity, fill level, and propellant management device (PMD) design. All the tanks were approximately 4.5 in. in diameter. The primary test measurements were the radial and tangential torques exerted on the tanks by the liquid. These torques could not be measured on the ground because of the masking effects of gravity.

Description: For the last several years, LeRC-HT, a three-dimensional computational fluid dynamics (CFD) computer code for analyzing gas turbine flow and convective heat transfer, has been evolving at the NASA Lewis Research Center. The code is unique in its ability to give a highly detailed representation of the flow field very close to solid surfaces. This is necessary for an accurate representation of fluid heat transfer and viscous shear stresses. The code has been used extensively for both internal cooling passage flows and hot gas path flows--including detailed film cooling calculations, complex tip-clearance gap flows, and heat transfer. In its current form, this code has a multiblock grid capability and has been validated for a number of turbine configurations. The code has been developed and used primarily as a research tool (at least 35 technical papers have been published relative to the code and its application), but it should be useful for detailed design analysis. We now plan to make this code available to selected users for further evaluation.

Description: The overall performance of a centrifugal compressor depends on the performance of the impeller and diffuser as well as on the interactions occurring between these components. Accurate measurements of the flow fields in each component are needed to develop computational models that can be used in compressor design codes. These measurements must be made simultaneously over an area that covers both components so that researchers can understand the interactions occurring between the two components. Optical measurement techniques are being used at the NASA Lewis Research Center to measure the velocity fields present in both the impeller and diffuser of a 4:1 pressure ratio centrifugal compressor operating at several conditions ranging from design flow to surge. Laser Doppler Velocimetry (LDV) was used to measure the intrablade flows present in the impeller, and the results were compared with analyses obtained from two three-dimensional viscous codes. The development of a region of low throughflow velocity fluid within this high-speed impeller was examined and compared with a similar region first observed in a large low-speed centrifugal impeller at Lewis. Particle Image Velocimetry (PIV) is a relatively new technique that has been applied to measuring the diffuser flow fields. PIV can collect data rapidly in the diffuser while avoiding the light-reflection problems that are often encountered when LDV is used. The Particle Image Velocimeter employs a sheet of pulsed laser light that is introduced into the diffuser in a quasi-radial direction through an optical probe inserted near the diffuser discharge. The light sheet is positioned such that its centerline is parallel to the hub and shroud surfaces and such that it is parallel to the diffuser vane, thereby avoiding reflections from the solid surfaces. Seed particles small enough to follow the diffuser flow are introduced into the compressor at an upstream location. A high-speed charge-coupled discharge (CCD) camera is synchronized to the laser pulse rate; this allows it to capture images of seed particle position that are separated by a small increment in time. A crosscorrelation of a particle's position in two consecutive images provides an estimate of flow velocity and direction. Multiple image pairs obtained in rapid succession at a particular flow condition provide enough measurements for statistical significance. PIV provides simultaneous velocity measurements over the entire plane that is illuminated by the light sheet instead of at a single point, as is the case when LDV is used. PIV has a further advantage in that the laser light pulse can be triggered by an external source such as a high-response pressure transducer. This feature will allow PIV to synchronize flow imaging to physical phenomena such as rotating stall or stall precursor waves. We hope that this technique can be used to obtain images of the flow field during and just prior to stall.

Description: Recently, refractive secondary solar concentrator systems were developed for solar thermal power and propulsion (ref. 1). Single-crystal oxides-such as yttria-stabilized zirconia (Y2O3-ZrO2), yttrium aluminum garnet (Y3Al5O12, or YAG), magnesium oxide (MgO), and sapphire (Al2O3)-are candidate refractive secondary concentrator materials. However, the refractive concentrator system will experience high-temperature thermal cycling in the solar thermal engine during the sun/shade transition of a space mission. The thermal mechanical reliability of these components in severe thermal environments is of great concern. Simulated mission tests are important for evaluating these candidate oxide materials under a variety of transient and steady-state heat flux conditions. In this research at the NASA Lewis Research Center, a controlled heat flux test approach was developed for investigating the thermal mechanical stability of the candidate oxide. This approach used a 3.0-kW continuous-wave (wavelength, 10.6 mm) carbon dioxide (CO2) laser (ref. 2). The CO2 laser is especially well-suited for single-crystal thermal shock tests because it can directly deliver well-characterized heat energy to the oxide surfaces. Since the oxides are opaque at the 10.6-mm wavelength of the laser beam, the light energy is absorbed at the surfaces rather than transmitting into the crystals, and thus generates the required temperature gradients within the specimens. The following figure is a schematic diagram of the test rig.

Description: The number of optical techniques that may potentially be used during a given wind tunnel test is continually growing. These include parameter sensitive paints that are sensitive to temperature or pressure, several different types of off-body and on-body flow visualization techniques, optical angle-of-attack (AoA), optical measurement of model deformation, optical techniques for determining density or velocity, and spectroscopic techniques for determining various flow field parameters. Often in the past the various optical techniques were developed independently of each other, with little or no consideration for other techniques that might also be used during a given test. Recently two optical techniques have been increasingly requested for production measurements in NASA wind tunnels. These are the video photogrammetric (or videogrammetric) technique for measuring model deformation known as the video model deformation (VMD) technique, and the parameter sensitive paints for making global pressure and temperature measurements. Considerations for, and initial attempts at, simultaneous measurements with the pressure sensitive paint (PSP) and the videogrammetric techniques have been implemented. Temperature sensitive paint (TSP) has been found to be useful for boundary-layer transition detection since turbulent boundary layers convect heat at higher rates than laminar boundary layers of comparable thickness. Transition is marked by a characteristic surface temperature change wherever there is a difference between model and flow temperatures. Recently, additional capabilities have been implemented in the target-tracking videogrammetric measurement system. These capabilities have permitted practical simultaneous measurements using parameter sensitive paint and video model deformation measurements that led to the first successful unified test with TSP for transition detection in a large production wind tunnel.

Description: Microscopic spheres suspended in liquid become highly ordered under the proper conditions. Such collections of particles, called colloidal suspensions or colloids, are the subject of a series of ongoing microgravity experiments at the NASA Lewis Research Center. By studying the way these colloidal suspensions order themselves, scientists can better understand how atoms of a liquid become ordered to form a solid. In addition, highly ordered colloids have special properties that may make them useful in future hightech applications. Work is underway at Lewis to develop an optical microscope to view these colloidal suspensions sphere by sphere in microgravity.

Description: Under the Advanced Subsonic Technology (AST) Program, an aeroelastic analysis code (TURBO-AE) based on Navier-Stokes equations is currently under development at NASA Lewis Research Center s Machine Dynamics Branch. For a blade row, aeroelastic instability can occur in any of the possible interblade phase angles (IBPA s). Analyzing small IBPA s is very computationally expensive because a large number of blade passages must be simulated. To reduce the computational cost of these analyses, we used time shifted, or phase-lagged, boundary conditions in the TURBO-AE code. These conditions can be used to reduce the computational domain to a single blade passage by requiring the boundary conditions across the passage to be lagged depending on the IBPA being analyzed. The time-shifted boundary conditions currently implemented are based on the direct-store method. This method requires large amounts of data to be stored over a period of the oscillation cycle. On CRAY computers this is not a major problem because solid-state devices can be used for fast input and output to read and write the data onto a disk instead of storing it in core memory.

Description: A complete "geometry to drag-polar" analysis capability for three-dimensional high-lift configurations is described. The approach is based on the use of unstructured meshes in order to enable rapid turnaround for complicated geometries which arise in high-lift con gurations. Special attention is devoted to creating a capability for enabling analyses on highly resolved grids. Unstructured meshes of several million vertices are initially generated on a work-station, and subsequently refined on a supercomputer. The flow is solved on these refined meshes on large parallel computers using an unstructured agglomeration multigrid algorithm. Good prediction of lift and drag throughout the range of incidences is demonstrated on a transport take-off configuration using up to 24.7 million grid points. The feasibility of using this approach in a production environment on existing parallel machines is demonstrated, as well as the scalability of the solver on machines using up to 1450 processors.

Description: Development of HPF versions of NPB and ARC3D showed that HPF has potential to be a high level language for parallelization of CFD applications. The use of HPF requires an intimate knowledge of the applications and a detailed analysis of data affinity, data movement and data granularity. Since HPF hides data movement from the user even with this knowledge it is easy to overlook pieces of the code causing low performance of the application. In order to simplify and accelerate the task of developing HPF versions of existing CFD applications we have designed and partially implemented ADAPT (Automatic Data Distribution and Placement Tool). The ADAPT analyzes a CFD application working on a single structured grid and generates HPF TEMPLATE, (RE)DISTRIBUTION, ALIGNMENT and INDEPENDENT directives. The directives can be generated on the nest level, subroutine level, application level or inter application level. ADAPT is designed to annotate existing CFD FORTRAN application performing computations on single or multiple grids. On each grid the application can considered as a sequence of operators each applied to a set of variables defined in a particular grid domain. The operators can be classified as implicit, having data dependences, and explicit, without data dependences. In order to parallelize an explicit operator it is sufficient to create a template for the domain of the operator, align arrays used in the operator with the template, distribute the template, and declare the loops over the distributed dimensions as INDEPENDENT. In order to parallelize an implicit operator, the distribution of the operator's domain should be consistent with the operator's dependences. Any dependence between sections distributed on different processors would preclude parallelization if compiler does not have an ability to pipeline computations. If a data distribution is "orthogonal" to the dependences of an implicit operator then the loop which implements the operator can be declared as INDEPENDENT.

Description: In this work the gas volume fraction and the root-mean-squared fluid velocity are measured in buoyancy driven shear flows of bubble suspensions in a tall, inclined, rectangular channel. The experiments are performed under conditions where We 〈〈 1a nd Re 〉〉 1, for which comparisons are made with kinetic theory and numerical simulations. Here Re = gamma(a(exp 2)/nu is the Reynolds number and We = rho(gamma(exp 2))a(exp 3)/sigma is the Weber number; gamma is the shear rate, a is the bubble radius, nu is the kinematic viscosity of the liquid, rho is the density of the liquid, and sigma is the surface tension of the gas/liquid interface. Kang et al. calculated the bubble phase pressure and velocity variance of sheared bubble suspensions under conditions where the bubbles are spherical and the liquid phase velocity field can be approximated using potential flow theory, i.e. We= 0 and Re 〉〉 1. Such conditions can be achieved in an experiment using gas bubbles, with a radius of O(0.5mm), in water. The theory requires that there be no average relative motion of the gas and liquid phases, hence the motivation for an experimental program in microgravity. The necessity of performing preliminary, Earth based experiments, however, requires performing experiments where the gas phase rises in the liquid, which significantly complicates the comparison of experiments with theory. Rather than comparing experimental results with theory for a uniform, homogeneous shear flow, experiments can be compared directly with solutions of the averaged equations of motion for bubble suspensions. This requires accounting for the significant lift force acting on the gas phase when the bubbles rise parallel to the average velocity of the sheared suspension. Shear flows can be produced in which the bubble phase pressure gradient, arising from shear induced collisions amongst the bubbles, balances a body force (centrifugal or gravitational) on the gas phase. A steady, non-uniform gas volume fraction can be measured, from which the bubble phase pressure gradient can be obtained and compared to theory and numerical simulations. The presence of bounding walls further complicates the experiments, since the detailed interactions of the bubbles with bounding walls is not well understood, especially in the presence of gravity, where the momentum and energy exchange depends on the inclination of the wall.

Description: A two-phase test loop has been designed and constructed to generate the necessary data for two-phase pressure drop and Critical Heat Flux (CHF) under reduced gravity conditions. A series of airplane trajectory tests aboard NASA KC-135 were performed and the data was used to evaluate the applicability of the earth gravity models for prediction of the reduced gravity data. Several commonly used correlations for the two-phase friction multiplier and critical heat flux were used to predict the data. It was generally concluded that the two-phase pressure drop can be predicted by the earth gravity correlations. The critical heat flux under reduced gravity conditions did not show a strong dependence on mass flow rate and the measured CHF were generally lower than the equivalent 1g conditions. The earth gravity models need to be modified for application to reduced gravities.

Description: In a wind tunnel facility, the direct measurement of forces and moments induced on the model are performed by a force measurement balance. The measurement balance is a precision-machined device that has strain gages at strategic locations to measure the strain (i.e., deformations) due to applied forces and moments. The strain gages convert the strain (and hence the applied force) to an electrical voltage that is measured by external meters. Thermal gradients can complicate the process, however. Thermal gradients on the balance cause differential expansion (or contraction) of various parts of the balance that induce a strain that is detected by the strain gages and is indistinguishable from an external applied force. The thermal gradients can result when testing is done at elevated temperatures or at cryogenic temperatures such as at the National Transonic Facility (NTF) at NASA Langley Research Center (LaRC).

Description: Boundary layer and aeroheating characteristics of several X-33 configurations have been experimentally examined in the Langley 20-Inch Mach 6 Air Tunnel. Global surface heat transfer distributions, surface streamline patterns, and shock shapes were measured on 0.013-scale models at Mach 6 in air. Parametric variations include angles-of-attack of 20-deg, 30-deg, and 40-deg; Reynolds numbers based on model length of 0.9 to 6.6 million; and body-flap deflections of 0, 10 and 20-deg. The effects of discrete and distributed roughness elements on boundary layer transition, which included trip height, size, location, and distribution, both on and off the windward centerline, were investigated. The discrete roughness results on centerline were used to provide a transition correlation for the X-33 flight vehicle that was applicable across the range of reentry angles of attack. The attachment line discrete roughness results were shown to be consistent with the centerline results, as no increased sensitivity to roughness along the attachment line was identified. The effect of bowed panels was qualitatively shown to be less effective than the discrete trips; however, the distributed nature of the bowed panels affected a larger percent of the aft-body windward surface than a single discrete trip.

Description: A parametric weight assessment of advanced metallic panel, ceramic blanket, and ceramic tile thermal protection systems (TPS) was conducted using an implicit, one-dimensional (1 -D) thermal finite element sizing code. This sizing code contained models to ac- count for coatings, fasteners, adhesives, and strain isolation pads. Atmospheric entry heating profiles for two vehicles, the Access to Space (ATS) rocket-powered single-stage-to-orbit (SSTO) vehicle and a proposed Reusable Launch Vehicle (RLV), were used to ensure that the trends were not unique to a particular trajectory. Eight TPS concepts were compared for a range of applied heat loads and substructural heat capacities to identify general trends. This study found the blanket TPS concepts have the lightest weights over the majority of their applicable ranges, and current technology ceramic tiles and metallic TPS concepts have similar weights. A proposed, state-of-the-art metallic system which uses a higher temperature alloy and efficient multilayer insulation was predicted to be significantly lighter than the ceramic tile systems and approaches blanket TPS weights for higher integrated heat loads.

Description: A new method has been developed to accelerate the convergence of explicit time-marching, laminar, Navier-Stokes codes through the combination of local preconditioning and multi-stage time marching optimization. Local preconditioning is a technique to modify the time-dependent equations so that all information moves or decays at nearly the same rate, thus relieving the stiffness for a system of equations. Multi-stage time marching can be optimized by modifying its coefficients to account for the presence of viscous terms, allowing larger time steps. We show it is possible to optimize the time marching scheme for a wide range of cell Reynolds numbers for the scalar advection-diffusion equation, and local preconditioning allows this optimization to be applied to the Navier-Stokes equations. Convergence acceleration of the new method is demonstrated through numerical experiments with circular advection and laminar boundary-layer flow over a flat plate.

Description: An experimental investigation was performed in the NASA Langley 16-Foot Transonic Tunnel to determine the aerodynamic effects of external convolutions, placed on the boattail of a nonaxisymmetric nozzle for drag reduction. Boattail angles of 15 and 22 were tested with convolutions placed at a forward location upstream of the boattail curvature, at a mid location along the curvature and at a full location that spanned the entire boattail flap. Each of the baseline nozzle afterbodies (no convolutions) had a parabolic, converging contour with a parabolically decreasing corner radius. Data were obtained at several Mach numbers from static conditions to 1.2 for a range of nozzle pressure ratios and angles of attack. An oil paint flow visualization technique was used to qualitatively assess the effect of the convolutions. Results indicate that afterbody drag reduction by convoluted contouring is convolution location, Mach number, boattail angle, and NPR dependent. The forward convolution location was the most effective contouring geometry for drag reduction on the 22 afterbody, but was only effective for M 〈 0.95. At M = 0.8, drag was reduced 20 and 36 percent at NPRs of 5.4 and 7, respectively, but drag was increased 10 percent for M = 0.95 at NPR = 7. Convoluted contouring along the 15 boattail angle afterbody was not effective at reducing drag because the flow was minimally separated from the baseline afterbody, unlike the massive separation along the 22 boattail angle baseline afterbody.

Description: The benchmark active controls technology and wind tunnel test program at NASA Langley Research Center was started with the objective to investigate the nonlinear, unsteady aerodynamics and active flutter suppression of wings in transonic flow. The paper will present the flutter suppression control law design process, numerical nonlinear simulation and wind tunnel test results for the NACA 0012 benchmark active control wing model. The flutter suppression control law design processes using (1) classical, (2) linear quadratic Gaussian (LQG), and (3) minimax techniques are described. A unified general formulation and solution for the LQG and minimax approaches, based on the steady state differential game theory is presented. Design considerations for improving the control law robustness and digital implementation are outlined. It was shown that simple control laws when properly designed based on physical principles, can suppress flutter with limited control power even in the presence of transonic shocks and flow separation. In wind tunnel tests in air and heavy gas medium, the closed-loop flutter dynamic pressure was increased to the tunnel upper limit of 200 psf. The control law robustness and performance predictions were verified in highly nonlinear flow conditions, gain and phase perturbations, and spoiler deployment. A non-design plunge instability condition was also successfully suppressed.

Description: This paper discusses a method for the identification and application of reduced-order models based on linear and nonlinear aerodynamic impulse responses. The Volterra theory of nonlinear systems and an appropriate kernel identification technique are described. Insight into the nature of kernels is provided by applying the method to the nonlinear Riccati equation in a non-aerodynamic application. The method is then applied to a nonlinear aerodynamic model of RAE 2822 supercritical airfoil undergoing plunge motions using the CFL3D Navier-Stokes flow solver with the Spalart-Allmaras turbulence model. Results demonstrate the computational efficiency of the technique.

Description: A three-dimensional large-eddy simulation model, TASS, is used to simulate the behavior of aircraft wake vortices in a real atmosphere. The purpose for this study is to validate the use of TASS for simulating the decay and transport of wake vortices. Three simulations are performed and the results are compared with the observed data from the 1994-1995 Memphis field experiments. The selected cases have an atmospheric environment of weak turbulence and stable stratification. The model simulations are initialized with appropriate meteorological conditions and a post roll-up vortex system. The behavior of wake vortices as they descend within the atmospheric boundary layer and interact with the ground is discussed.

Description: USM3D is a widely-used unstructured flow solver for simulating inviscid and viscous flows over complex geometries. The current version (version 5.0) of USM3D, however, does not have advanced turbulence models to accurately simulate complicated flows. We have implemented two modified versions of the original Jones and Launder k-epsilon two-equation turbulence model and the Girimaji algebraic Reynolds stress model in USM3D. Tests have been conducted for two flat plate boundary layer cases, a RAE2822 airfoil and an ONERA M6 wing. The results are compared with those of empirical formulae, theoretical results and the existing Spalart-Allmaras one-equation model.

Description: The spreading characteristics of jets from several asymmetric nozzles, and a set of rectangular orifices are compared, covering a jet Mach number range of 0.3-2.0. The effect of 'tabs' for a rectangular and a round nozzle is also included in the comparison. Compared to a round jet, the jets from the asymmetric nozzles spread only slightly more at subsonic conditions whereas at supersonic conditions, when 'screech' occurs, they spread much more. The dynamics of the azimuthal vortical structures of the jet, organized and intensified under the screeching condition, are thought to be responsible for the observed effect at supersonic conditions. Curiously, the jet from a 'lobed' nozzle spreads much less at supersonic condition compared to all other cases; this is due to the absence of screech with this nozzle. Screech stages inducing flapping, rather than varicose or helical, flow oscillation cause a more pronounced jet spreading. At subsonic conditions, only a slight increase in jet spreading with the asymmetric nozzles contrasts previous observations by others. The present results show that the spreading of most asymmetric jets is not much different from that of a round jet. This inference is further supported by data from the rectangular orifices. In fact, jets from the orifices with small aspect ratio (AR) exhibit virtually no increase in the spreading. A noticeable increase commences only when AR is larger than about 10. Thus, 'shear layer perimeter stretching', achieved with a larger AR for a given cross-sectional area of the orifice, by itself, proves to be a relatively inefficient mechanism for increasing jet spreading. In contrast, the presence of streamwise vortices or 'natural excitation' can cause a significant increase - effects that might explain the observations in the previous investigations. Thus far, the biggest increase in jet spreading is observed with the tabs. This is true in the subsonic regime, as well as in the supersonic regime, in spite of the fact that screech is eliminated by the tabs. The characteristic spreading of the tabbed jets is explained by the induced motion of the tab-generated streamwise vortex pairs. The tabs, however, incur thrust loss; the flow blockage and loss in thrust coefficient, vis-a-vis the spreading increase, are evaluated for various configurations.

Description: The Ninth Thermal and Fluids Analysis Workshop (TFAWS 98) was held at the Ohio Aerospace Institute in Cleveland, Ohio from August 31 to September 4, 1998. The theme for the hands-on training workshop and conference was "Integrating Computational Fluid Dynamics and Heat Transfer into the Design Process." Highlights of the workshop (in addition to the papers published herein) included an address by the NASA Chief Engineer, Dr. Daniel Mulville; a CFD short course by Dr. John D. Anderson of the University of Maryland; and a short course by Dr. Robert Cochran of Sandia National Laboratories. In addition, lectures and hands-on training were offered in the use of several cutting-edge engineering design and analysis-oriented CFD and Heat Transfer tools. The workshop resulted in international participation of over 125 persons representing aerospace and automotive industries, academia, software providers, government agencies, and private corporations. The papers published herein address issues and solutions related to the integration of computational fluid dynamics and heat transfer into the engineering design process. Although the primary focus is aerospace, the topics and ideas presented are applicable to many other areas where these and other disciplines are interdependent.

Description: This paper reviews the effects of 'screech', 'asymmetric nozzle shaping', 'tabs' and 'overexpansion' on the spreading of free jets. Corresponding thrust penalty for the tabs and overexpanded condition are also evaluated. The asymmetric shapes include rectangular ones with varying aspect ratio. Tabs investigated are triangular shaped 'delta-tabs' placed at the exit of a convergent circular nozzle. The effect of overexpansion is examined with circular convergent-divergent (C-D) nozzles. Tabs and overexpansion are found to yield the largest increase in jet spreading. Each, however, involves a performance penalty, i.e., a loss in thrust coefficient. Variation of the size of four delta-tabs show that there exists an optimum size for which the gain in jet spreading is the maximum per unit loss in thrust coefficient. With the C-D nozzles, the minimum in thrust coefficient is expected near the beginning of the overexpanded regime based on idealized flow calculations. The maximum increase in jet spreading, however, is found to occur at higher pressure ratios well into the overexpanded regime. The optimum benefit with the overexpanded flow, in terms of gain in spreading for unit penalty, is found to be comparable to the optimum tab case.

Description: The effect of dynamic rolling oscillations of delta-wing/twin-tail configuration on twin-tail buffet response is investigated. The computational model consists of a sharp-edged delta wing of aspect ratio one and swept-back flexible twin tail with taper ratio of 0.23. The configuration model is statically pitched at 30 deg. angle of attack and then forced to oscillate in roll around the symmetry axis at a constant amplitude of 4 deg. and reduced frequency of pi and 2(pi). The freestream Mach number and Reynolds number are 0.3 and 1.25 million, respectively. This multidisciplinary problem is solved using three sets of equations on a dynamic multi-block grid structure. The first set is the unsteady, full Navier-Stokes equations, the second set is the aeroelastic equations for coupled bending and torsion vibrations of the tails, and the third set is the grid-displacement equations. The configuration is investigated for inboard position of the twin tails which corresponds to a separation distance between the twin tails of 33% wing span. The computed results are compared with the results of stationary configuration, which previously have been validated using experimental data. The results conclusively showed that the rolling oscillations of the configuration have led to higher loads, higher deflections, and higher excitation peaks than those of the stationary configuration. Moreover, increasing the reduced frequency has led to higher loads and excitation peaks and lower bending and torsion deflections and acceleration.

Description: Continuous improvement of aerospace product development processes is a driving requirement across much of the aerospace community. As up to 90% of the cost of an aerospace product is committed during the first 10% of the development cycle, there is a strong emphasis on capturing, creating, and communicating better information (both requirements and performance) early in the product development process. The community has responded by pursuing the development of computer-based systems designed to enhance the decision-making capabilities of product development individuals and teams. Recently, the historical foci on sharing the geometrical representation and on configuration management are being augmented: Physics-based analysis tools for filling the design space database; Distributed computational resources to reduce response time and cost; Web-based technologies to relieve machine-dependence; and Artificial intelligence technologies to accelerate processes and reduce process variability. Activities such as the Advanced Design Technologies Testbed (ADTT) project at NASA Ames Research Center study the strengths and weaknesses of the technologies supporting each of these trends, as well as the overall impact of the combination of these trends on a product development event. Lessons learned and recommendations for future activities will be reported.

Description: The overall objective of this work is to characterize the flow to rough showerheads by deriving pressure drop versus velocity correlations to at can be then used in reactor scale simulations where the showerhead is approximated as a porous medium. At relatively low Reynolds numbers (less than 1-10 based on the hole length scale) and in the absence of slip flow, Darcy's Law, grad P = mu U/k, can be used to express the relation between the pressure drop and velocity where @mu@ is the fluid viscosity and it is the permeability that can be theoretically predicted as k= e R^2 /8, where e is the porosity. However, at sufficiently small hole diameters and decreased pressures (less than 5 Torr), the Knudsen number based on showerhead tube radius increases, and the flow may be in a transition regime. Different expressions have been proposed to account for this effect in the permeability by expressing k as a function of either pressure or Knudsen number. But at even higher Knudsen numbers, the pressure drop - velocity dependence is non-linear, and Darcy's Law no longer holds such that a permeability cannot be defined. The direct simulation Monte Carlo method is used along side conventional CFD techniques to determine the extent to which the CFD technique is appropriate and helps to derive correlations for the more rarefied cases of interest in these showerhead flows.

Description: This conference presents information to the scientific community on research results, future directions, and research opportunities in microgravity fluid physics and transport phenomena within NASA's microgravity research program. The conference theme is "The International Space Station." Plenary sessions provide an overview of the Microgravity Fluid Physics Program, the International Space Station and the opportunities ISS presents to fluid physics and transport phenomena researchers, and the process by which researchers may become involved in NASA's program, including information about the NASA Research Announcement in this area. Two plenary lectures present promising areas of research in electrohydrodynamics/electrokinetics in the movement of particles and in micro- and meso-scale effects on macroscopic fluid dynamics. Featured speakers in plenary sessions present results of recent flight experiments not heretofore presented. The conference publication consists of this book of abstracts and the full Proceedings of the 4th Microgravity Fluid Physics and Transport Phenomena Conference on CD-ROM, containing full papers presented at the conference (NASA/CP-1999-208526/SUPPL1).

Description: Following is a brief summary of a theoretical investigation of material (or constitutive) instability associated with shear induced particle migration in dense particulate suspensions or granular media. It is shown that one can obtain a fairly general linear-stability analysis, including the effects of shear-induced anisotropy in the base flow as well as Reynolds dilatancy. A criterion is presented here for simple shearing instability in the absence of inertia and dilatancy.

Description: Data assimilation methods are routinely used in oceanography. The statistics of the model and measurement errors need to be specified a priori. This study addresses the problem of estimating model and measurement error statistics from observations. We start by testing innovation based methods of adaptive error estimation with low-dimensional models in the North Pacific (5-60 deg N, 132-252 deg E) to TOPEX/POSEIDON (TIP) sea level anomaly data, acoustic tomography data from the ATOC project, and the MIT General Circulation Model (GCM). A reduced state linear model that describes large scale internal (baroclinic) error dynamics is used. The methods are shown to be sensitive to the initial guess for the error statistics and the type of observations. A new off-line approach is developed, the covariance matching approach (CMA), where covariance matrices of model-data residuals are "matched" to their theoretical expectations using familiar least squares methods. This method uses observations directly instead of the innovations sequence and is shown to be related to the MT method and the method of Fu et al. (1993). Twin experiments using the same linearized MIT GCM suggest that altimetric data are ill-suited to the estimation of internal GCM errors, but that such estimates can in theory be obtained using acoustic data. The CMA is then applied to T/P sea level anomaly data and a linearization of a global GFDL GCM which uses two vertical modes. We show that the CMA method can be used with a global model and a global data set, and that the estimates of the error statistics are robust. We show that the fraction of the GCM-T/P residual variance explained by the model error is larger than that derived in Fukumori et al.(1999) with the method of Fu et al.(1993). Most of the model error is explained by the barotropic mode. However, we find that impact of the change in the error statistics on the data assimilation estimates is very small. This is explained by the large representation error, i.e. the dominance of the mesoscale eddies in the T/P signal, which are not part of the 21 by 1" GCM. Therefore, the impact of the observations on the assimilation is very small even after the adjustment of the error statistics. This work demonstrates that simult&neous estimation of the model and measurement error statistics for data assimilation with global ocean data sets and linearized GCMs is possible. However, the error covariance estimation problem is in general highly underdetermined, much more so than the state estimation problem. In other words there exist a very large number of statistical models that can be made consistent with the available data. Therefore, methods for obtaining quantitative error estimates, powerful though they may be, cannot replace physical insight. Used in the right context, as a tool for guiding the choice of a small number of model error parameters, covariance matching can be a useful addition to the repertory of tools available to oceanographers.

Description: The flight control of X-33 poses a challenge to conventional gain-scheduled flight controllers due to its large attitude maneuvers from liftoff to orbit and reentry. In addition, a wide range of uncertainties in vehicle handling qualities and disturbances must be accommodated by the attitude control system. Nonlinear tracking and decoupling control by trajectory linearization can be viewed as the ideal gain-scheduling controller designed at every point on the flight trajectory. Therefore it provides robust stability and performance at all stages of flight without interpolation of controller gains, and eliminates costly controller redesigns due to minor airframe alteration or mission reconfiguration. A prototype trajectory linearization design for X-33 ascent flight controller was designed and tested with 3-DOF and 6-DOF simulations during the 10 weeks SFFP. It is noted that the 6-DOF results were obtained from the 3-DOF design with only a few hours of tuning, which demonstrates the inherent robustness of the design technique. It is this "plug-and-play" feature that is much needed by NASA for the development, test and routine operations of the RLVs. Plans for further research are also presented.

Description: Turbulence attenuation by greater than a factor of two has been observed in many practical gas flows carrying volume fractions as small as 0.01% of dispersed particles. Particles which cause such attenuation usually are smaller than the smallest scales of the turbulence and have time constants 5 to 10 times greater than the time scale of a typical turbulent eddy. That is, strongly attenuating particles usually have Stokes numbers in the range of 5 to 10, indicating that they do not respond to the turbulent fluctuations, but instead just fall through the flow responding only to the mean flow. There are two mechanisms by which free falling particles may attenuate turbulence. First, the unresponsive particles act as a drag on the turbulent eddies, passing energy from the turbulent eddies to the small scale wakes of the particles where it is quickly dissipated by viscosity. The second mechanism is more complicated. Particles falling under gravity convert gravitational potential energy to turbulent velocity fluctuations. If the particles are large, this mechanism increases the overall turbulence level. However, with moderate size particles, the small scale turbulence generated apparently distorts the turbulent eddies leading to more rapid dissipation. Unfortunately, this conclusion is supported only by circumstantial evidence to date. The objectives of the experiment are to use microgravity to separate the two mechanisms. A region of nearly-isotropic decaying turbulence with zero mean flow will be formed in a box in the microgravity environment. Different sets of particles with Stokes numbers in the range of 2 to 20 will be dispersed in the flow. With zero gravity and no mean fluid velocity the particles will have zero mean velocity. With the large Stokes numbers, the fluctuating velocities will also be small. Therefore, the only attenuation mechanism will be the direct action of the particles on the turbulence. Control experiments will also be done in which the particles fall through the measurement volume. Measurements will be acquired using a high resolution image velocimetry (PIV) system being developed specifically for work in particle-laden flows. The measurements will include the decay of the turbulence kinetic energy under various particle loadings. The spatial spectra of the turbulence will also be measured. In a second set of experiments, the interaction of a single eddy with a collection of nearly stationary particles will be examined. The eddy will be a vortex ring emitted by a jet pulse through an orifice. The distortion of the vortex under the influence of the particles will be examined to gain a better understanding of how fine particles can cause such large reductions in turbulence levels. This experiment could not be conducted in terrestrial gravity because the high particle velocities would overwhelm the relatively low speed motion of the vortex ring. This experimental program is just getting underway. The initial challenge is to build a closed facility containing reasonably homogeneous and isotropic turbulence with zero mean velocity. Our approach is to use a set of synthetic jets mounted on the periphery of a transparent plexiglass box to create the turbulence. A synthetic jet is a plenum chamber with an orifice open to the volume of interest. The volume of the chamber fluctuates periodically so alternately a jet is ejected from the volume or flow is drawn back in as a sink. The asymmetry of this situation results in a net transport of momentum and kinetic energy into the volume of interest. The present apparatus includes eight synthetic jets each powered independently by a six inch loudspeaker. The synthetic jets discharge through ejector tubes to increase the scale of the turbulence. Construction of the apparatus is now complete and preliminary flow visualization studies have been conducted. The PIV system is also under development. A compact dual-pulse YAG laser has been acquired as the light source and special software is under development to allow simultaneous measurements of both the particle phase and the fluid phase (marked by fine tracers).

Description: The Aeroelasticity Task is intended to provide demonstrated technology readiness to predict and improve flutter characteristics of an HSCT configuration. This requires aerodynamic codes that are applicable to the wide range of flight regimes in which the HSCT will operate, and are suitable to provide the higher fidelity required for evaluation of aeroservoelastic coupling effects. Prediction of these characteristics will result in reduced airplane weight and risk associated with a highly flexible, low-aspect ratio supersonic airplane with narrow fuselage, relatively thin wings, and heavy engines. This Task is subdivided into three subtasks. The first subtask includes the design, fabrication, and testing of wind-tunnel models suitable to provide an experimental database relevant to HSCT configurations. The second subtask includes validation of candidate unsteady aerodynamic codes, applicable in the Mach and frequency ranges of interest for the HSCT, through analysis test correlation with the test data. The third subtask includes efforts to develop and enhance these codes for application to HSCT configurations. The wind tunnel models designed and constructed during this program furnished data which were useful for the analysis test correlation work but there were shortcomings. There was initial uncertainty in the proper tunnel configuration for testing, there was a need for higher quality measured model geometry, and there was a need for better measured model displacements in the test data. One of the models exhibited changes in its dynamic characteristics during testing. Model design efforts were hampered by a need for more and earlier analysis support and better knowledge of material properties. Success of the analysis test correlation work was somewhat muted by the uncertainties in the wind tunnel model data. The planned extent of the test data was not achieved, partly due to the delays in the model design and fabrication which could not be extended due to termination of the HSR program.

Description: The local viscous-inviscid interaction field generated by a wall temperature jump on a flat plate in supersonic flow and on the windside of a Reusable Launch Vehicle in hypersonic flow is studied in detail by both a Navier-Stokes numerical code and an analytical triple-deck model. Treatment of the rapid heat transfer changes both upstream and downstream of the jump is included. Closed form relationships derived from the triple-deck theory are presented. The analytically predicted pressure and heating variations including upstream influence are found to be in generally good agreement with the Computational Fluid Dynamic (CFD) predictions. These analyses not only clarify the interactive physics involved but also are useful in preliminary design of thermal protection systems and as an insertable module to improve CFD code efficiency when applied to such small-scale interaction problems. The analyses only require conditions at the wall and boundary-layer edge which are easily extracted from a baseline, constant wall temperature, CFD solution.

Description: With the recent interest in novel control effectors there is a need to determine the stability and control derivatives of new aircraft configurations early in the design process. These derivatives are central to most control law design methods and would allow the determination of closed-loop control performance of the vehicle. Early determination of the static and dynamic behavior of an aircraft may permit significant improvement in configuration weight, cost, stealth, and performance through multidisciplinary design. The classical method of determining static stability and control derivatives - constructing and testing wind tunnel models - is expensive and requires a long lead time for the resultant data. Wind tunnel tests are also limited to the preselected control effectors of the model. To overcome these shortcomings, computational fluid dynamics (CFD) solvers are augmented via automatic differentiation, to directly calculate the stability and control derivatives. The CFD forces and moments are differentiated with respect to angle of attack, angle of sideslip, and aircraft shape parameters to form these derivatives. A subset of static stability and control derivatives of a tailless aircraft concept have been computed by two differentiated inviscid CFD codes and verified for accuracy with central finite-difference approximations and favorable comparisons to a simulation database.

Description: The Line Integral Convolution (LIC) algorithm has received a lot of attention and interest. Yet, only a few of the current LIC related algorithms deal specifically with color textures for automatic detection of flow features. This paper provides an overview of research in this area.

Description: The present invention provides an ozone generation and delivery system that lends itself to small scale applications and requires very low maintenance. The system includes an anode reservoir and a cathode phase separator each having a hydrophobic membrane to allow phase separation of produced gases from water. The system may be configured to operate passively with no moving parts or in a self-pressurizing manner with the inclusion of a pressure controlling device or valve in the gas outlet of the anode reservoir. The hydrogen gas, ozone gas and water containing ozone may be delivered under pressure.

Description: Gas bubbles driven in radial oscillations are subject to an instability of the spherical shape that is opposed by surface tension and viscosity. An exact linear formulation for the study of the phenomenon has been available for many years, but its complexity has discouraged a detailed investigation. With the recent theory of sonoluminescence of Lohse and co-workers, there has been a renewed interest in the problem and new data have become available. This paper presents a numerical method for the solution of the pertinent equations and compares the theory with these new data. The coupling of the strong nonlinearity of the bubble radial oscillations with the parametric mechanism of the surface instability results in a very complex structure for the stability boundary. Nevertheless, a good agreement between theory and data is found. A comparison with earlier approximate models is also made.

Description: A steady, two dimensional cellular convection modifies the morphological instability of a binary alloy that undergoes directional solidification. When the convection wavelength is far longer than that of the morphological cells, the behavior of the moving front is described by a slow, spatial-temporal dynamics obtained through a multiple-scale analysis. The resulting system has a "parametric-excitation" structure in space, with complex parameters characterizing the interactions between flow, solute diffusion, and rejection. The convection stabilizes two dimensional disturbances oriented with the flow, but destabilizes three dimensional disturbances in general. When the flow is weak, the morphological instability behaves incommensurably to the flow wavelength, but becomes quantized and forced to fit into the flow-box as the flow gets stronger. At large flow magnitudes the instability is localized, confined in narrow envelopes with cells traveling with the flow. In this case the solutions are discrete eigenstates in an unbounded space. Their stability boundary and asymptotics are obtained by the WKB analysis.

Description: Foams are extremely important in a variety of industrial applications. Foams are widely used in fire-fighting applications, and are especially effective in fighting flammable liquid fires. In fact the Fire Suppression System aboard the Space Shuttle utilizes cylinders of Halon foam, which, when fired, force a rapidly expanding foam into the convoluted spaces behind instrument panels. Foams are critical in the process of enhanced oil recovery, due to their surface-active and highly viscous nature. They are also used as drilling fluids in underpressurized geologic formations. They are used as transport agents, and as trapping agents. They are also used as separation agents, where ore refinement is accomplished by froth flotation of the typically lighter and hydrophobic contaminants. The goal of the proposed investigation is the determination of the mechanical and rheological properties of foams, utilizing the microgravity environment to explore foam rheology for foams which cannot exist, or only exist for a short time, in 1g.

Description: Integrated modeling of spacecraft systems is a rapidly evolving area in which multidisciplinary models are developed to design and analyze spacecraft configurations. These models are especially important in the early design stages where rapid trades between subsystems can substantially impact design decisions. Integrated modeling is one of the cornerstones of two of NASA's planned missions in the Origins Program -- the Next Generation Space Telescope (NGST) and the Space Interferometry Mission (SIM). Common modeling tools for control design and opto-mechanical analysis have recently emerged and are becoming increasingly widely used. A discipline that has been somewhat less integrated, but is nevertheless of critical concern for high precision optical instruments, is thermal analysis and design. A major factor contributing to this mild estrangement is that the modeling philosophies and objectives for structural and thermal systems typically do not coincide. Consequently the tools that are used in these discplines suffer a degree of incompatibility, each having developed along their own evolutionary path. Although standard thermal tools have worked relatively well in the past. integration with other disciplines requires revisiting modeling assumptions and solution methods. Over the past several years we have been developing a MATLAB based integrated modeling tool called IMOS (Integrated Modeling of Optical Systems) which integrates many aspects of structural, optical, control and dynamical analysis disciplines. Recent efforts have included developing a thermal modeling and analysis capability, which is the subject of this article. Currently, the IMOS thermal suite contains steady state and transient heat equation solvers, and the ability to set up the linear conduction network from an IMOS finite element model. The IMOS code generates linear conduction elements associated with plates and beams/rods of the thermal network directly from the finite element structural model. Conductances for temperature varying materials are accommodated. This capability both streamlines the process of developing the thermal model from the finite element model, and also makes the structural and thermal models compatible in the sense that each structural node is associated with a thermal node. This is particularly useful when the purpose of the analysis is to predict structural deformations due to thermal loads. The steady state solver uses a restricted step size Newton method, and the transient solver is an adaptive step size implicit method applicable to general differential algebraic systems. Temperature dependent conductances and capacitances are accommodated by the solvers. In addition to discussing the modeling and solution methods. applications where the thermal modeling is "in the loop" with sensitivity analysis, optimization and optical performance drawn from our experiences with the Space Interferometry Mission (SIM), and the Next Generation Space Telescope (NGST) are presented.

Description: A brief overview of a cooperative NASA/Boeing research effort, Strake Technology Research Application to Transport Aircraft (STRATA), intended to explore the potential of applying forebody strake technology to transport aircraft configurations for directional stability and control at low angles of attack, is presented. As an initial step in the STRATA program, an exploratory wind-tunnel investigation of the effect of fixed forebody strakes on the directional stability and control characteristics of a generic transport configuration was conducted in the NASA Langley 12-Foot Low-Speed Wind Tunnel. Results of parametric variations in strake chord and span, as well as the effect of strake incidence, are presented. The use of strakes for yaw control is also discussed. Results emphasize the importance of forebody/fuselage crossflow in influencing strake effectiveness. Strake effectiveness is also seen to be directly related to its span, but less sensitive to chord; a very short-chord strake with sufficient span can have a significant effect.

Description: In the present investigation, the results obtained during the ground test of a closed-loop control system conducted on a full-scale fighter to attenuate vertical fin buffeting response using strain actuation are presented. Two groups of actuators consisting of piezoelectric elements distributed over the structure were designed to achieve authority over the first and second modes of the vertical fin. The control laws were synthesized using the Linear Quadratic Gaussian (LQG) method for a time-invariant control system. Three different pairs of sensors including strain gauges and accelerometers at different locations were used to close the feedback loop. The results demonstrated that measurable reductions in the root-mean-square (RMS) values of the fin dynamic response identified by the strain transducer at the critical point for fatigue at the root were achieved under the most severe buffet condition. For less severe buffet conditions, reductions of up to 58% were achieved.

Description: The application of an unstructured grid methodology on a three-dimensional high-lift configuration is presented. The focus of this paper is on the grid generation aspect of an integrated effort for the development of an unstructured-grid computational fluid dynamics (CFD) capability at the NASA Langley Research Center. The meshing approach is based on tetrahedral grids generated by the advancing-front and the advancing-layers procedures. The capability of the method for solving high-lift problems is demonstrated on an aircraft model referred to as the energy efficient transport configuration. The grid generation issues, including the pros and cons of the present approach, are discussed in relation to the high-lift problems. Limited viscous flow results are presented to demonstrate the viability of the generated grids. A corresponding Navier-Stokes solution capability, along with further computations on the present grid, is presented in a companion SAE paper.

Description: The Thermal Radiation Group at Virginia Polytechnic Institute and State University has been working closely with scientists and engineers at NASA's Langley Research Center to develop accurate analytical and numerical models suitable for designing next generation thin-film thermal radiation detectors for earth radiation budget measurement applications. The current study provides an analytical model of the notional thermal radiation detector that takes into account thermal transport phenomena, such as the contact resistance between the layers of the detector, and is suitable for use in parameter estimation. It was found that the responsivity of the detector can increase significantly due to the presence of contact resistance between the layers of the detector. Also presented is the effect of doping the thermal impedance layer of the detector with conducting particles in order to electrically link the two junctions of the detector. It was found that the responsivity and the time response of the doped detector decrease significantly in this case. The corresponding decrease of the electrical resistance of the doped thermal impedance layer is not sufficient to significantly improve the electrical performance of the detector. Finally, the "roughness effect" is shown to be unable to explain the decrease in the thermal conductivity often reported for thin-film layers.

Description: A method and apparatus for cold gas reinjection in through-flow and reverse-flow wave rotors having a plurality of channels formed around a periphery thereof. A first port injects a supply of cool air into the channels. A second port allows the supply of cool air to exit the channels and flow to a combustor. A third port injects a supply of hot gas from the combustor into the channels. A fourth port allows the supply of hot gas to exit the channels and flow to a turbine. A diverting port and a reinjection port are connected to the second and third ports, respectively. The diverting port diverts a portion of the cool air exiting through the second port as reinjection air. The diverting port is fluidly connected to the reinjection port which reinjects the reinjection air back into the channels. The reinjection air evacuates the channels of the hot gas resident therein and cools the channel walls, a pair of end walls of the rotor, ducts communicating with the rotor and subsequent downstream components. In a second embodiment, the second port receives all of the cool air exiting the channels and the diverting port diverts a portion of the cool air just prior to the cool air flowing to the combustor.

Description: A fractional step method for the solution of steady and unsteady incompressible Navier-Stokes equations is outlined. The method is based on a finite volume formulation and uses the pressure in the cell center and the mass fluxes across the faces of each cell as dependent variables. Implicit treatment of convective and viscous terms in the momentum equations enables the numerical stability restrictions to be relaxed. The linearization error in the implicit solution of momentum equations is reduced by using three subiterations in order to achieve second order temporal accuracy for time-accurate calculations. In spatial discretizations of the momentum equations, a high-order (3rd and 5th) flux-difference splitting for the convective terms and a second-order central difference for the viscous terms are used. The resulting algebraic equations are solved with a line-relaxation scheme which allows the use of large time step. A four color ZEBRA scheme is employed after the line-relaxation procedure in the solution of the Poisson equation for pressure. This procedure is applied to a Couette flow problem using a distorted computational grid to show that the method minimizes grid effects. Additional benchmark cases include the unsteady laminar flow over a circular cylinder for Reynolds Numbers of 200, and a 3-D, steady, turbulent wingtip vortex wake propagation study. The solution algorithm does a very good job in resolving the vortex core when 5th-order upwind differencing and a modified production term in the Baldwin-Barth one-equation turbulence model are used with adequate grid resolution.

Description: This paper analytically examines the unsteady fluid dynamics of a vortex filament subsequent to a normal collision of the vortex with a solid body. In particular, the breakdown or reconnection phenomena, post-collision, for a vortex filament is studied. The paper does not investigate the collision dynamics process itself. The derived exact solution is based upon the laminar viscous form of the Helmholtz equations.

Description: In spite of a superficial similarity with gas bubbles, the intimate coupling between dynamical and thermal processes confers to oscillating vapor bubbles some unique characteristics. This paper examines numerically the validity of some asymptotic-theory predictions such as the existence of two resonant radii and a limit size for a given sound amplitude and frequency. It is found that a small vapor bubble in a sound field of sufficient amplitude grows quickly through resonance and continues to grow thereafter at a very slow rate, seemingly indefinitely. Resonance phenomena therefore play a role for a few cycles at most, and reaching a limit size-if one exists at all-is found to require far more than several tens of thousands of cycles. It is also found that some small bubbles may grow or collapse depending on the phase of the sound field. The model accounts in detail for the thermo-fluid-mechanic processes in the vapor. In the second part of the paper, an approximate formulation valid for bubbles small with respect to the thermal penetration length in the vapor is derived and its accuracy examined, The present findings have implications for acoustically enhanced boiling heat transfer and other special applications such as boiling in microgravity.

Description: Condition queries on distributed data ask where particular conditions are satisfied. It is possible to represent condition queries as geometric objects by plotting field data in various spaces derived from the data, and by selecting loci within these derived spaces which signify the desired conditions. Rather simple geometric partitions of derived spaces can represent complex condition queries because much complexity can be encapsulated in the derived space mapping itself A geometric view of condition queries provides a useful conceptual unification, allowing one to intuitively understand many existing vector field feature detection algorithms -- and to design new ones -- as variations on a common theme. A geometric representation of condition queries also provides a simple and coherent basis for computer implementation, reducing a wide variety of existing and potential vector field feature detection techniques to a few simple geometric operations.

Description: This paper describes an experiment in which a large-scale scientific application development for tightly-coupled parallel machines is adapted to the distributed execution environment of the Information Power Grid (IPG). A brief overview of the IPG and a description of the computational fluid dynamics (CFD) algorithm are given. The Globus metacomputing toolkit is used as the enabling device for the geographically-distributed computation. Modifications related to latency hiding and Load balancing were required for an efficient implementation of the CFD application in the IPG environment. Performance results on a pair of SGI Origin 2000 machines indicate that real scientific applications can be effectively implemented on the IPG; however, a significant amount of continued effort is required to make such an environment useful and accessible to scientists and engineers.

Description: A compact modulated birefringence-measurement system has been developed for use in microgravity fluid physics applications with non-Newtonian fluids such as polymer solutions. This instrument uses a dual-crystal transverse electro-optical modulator capable of modulation frequencies in excess of 100 MHz. The two crystals are modulated 180 deg. out of phase from each other, The theoretical framework governing the development of this instrument using the Mueller-Stokes polarization matrices is discussed. Several ground-based experiments are performed to compare this system with the theoretical results. Results from this transverse electro-optical modulator-based birefringence-measurement system agree well with the theory. The instrument is also very stable and robust, making it suitable for the extreme acceleration environment to be encountered in a NASA Black Brandt sounding rocket.

Description: The creeping flow of a dilute (0.025 wt%) monodisperse polystyrene/polystyrene Boger fluid through a 4:1:4 axisymmetric contraction/expansion is experimentally observed for a wide range of Deborah numbers. Pressure drop measurements across the orifice plate show a large extra pressure drop that increases monotonically with Deborah number above the value observed for a similar Newtonian fluid at the same flow rate. This enhancement in the dimensionless pressure drop is not associated with the onset of a flow instability, yet it is not predicted by existing steady-state or transient numerical computations with simple dumbbell models. It is conjectured that this extra pressure drop is the result of an additional dissipative contribution to the polymeric stress arising from a stress-conformation hysteresis in the strong non-homogeneous extensional flow near the contraction plane. Such a hysteresis has been independently measured and computed in recent studies of homogeneous transient uniaxial stretching of PS/PS Boger fluids. Flow visualization and velocity field measurements using digital particle image velocimetry (DPIV) show large upstream growth of the corner vortex with increasing Deborah number. At large Deborah numbers, the onset of an elastic instability is observed, first locally as small amplitude fluctuations in the pressure measurements, and then globally as an azimuthal precessing of the upstream corner vortex accompanied by periodic oscillations in the pressure drop across the orifice.

Description: Scale-similar models employ multiple filtering operations to identify the smallest resolved scales, which have been shown to be the most active in the interaction with the unresolved subgrid scales. They do not assume that the principal axes of the strain-rate tensor are aligned with those of the subgrid-scale stress (SGS) tensor, and allow the explicit calculation of the SGS energy. They can provide backscatter in a numerically stable and physically realistic manner, and predict SGS stresses in regions that are well correlated with the locations where large Reynolds stress occurs. In this paper, eddy viscosity and mixed models, which include an eddy-viscosity part as well as a scale-similar contribution, are applied to the simulation of two flows, a high Reynolds number plane channel flow, and a three-dimensional, nonequilibrium flow. The results show that simulations without models or with the Smagorinsky model are unable to predict nonequilibrium effects. Dynamic models provide an improvement of the results: the adjustment of the coefficient results in more accurate prediction of the perturbation from equilibrium. The Lagrangian-ensemble approach [Meneveau et al., J. Fluid Mech. 319, 353 (1996)] is found to be very beneficial. Models that included a scale-similar term and a dissipative one, as well as the Lagrangian ensemble averaging, gave results in the best agreement with the direct simulation and experimental data.

Description: Experiments were designed, implemented, and evaluated in the thermal and fluid sciences at the NASA Langley Research Center. This research was conducted cooperatively with NASA employees using, where necessary, equipment and facilities provided by the U.S. Government. The research fell within the scope of the University Agreement between the NASA Langley Research Center and The George Washington University for Joint Research and Education Projects dated June 7, 8, 1994, which continues the Joint Institute for the Advancement of Flight Sciences (JIAFS).

Description: In the largest solar energetic particle (SEP) events, acceleration occurs at shock waves driven out from the Sun by coronal mass ejections (CMEs). Peak particle intensities are a strong function of CME speed, although the intensities, spectra, and angular distributions of particles escaping the shock are highly modified by scattering on Alfven waves produced by the streaming particles themselves. Element abundances vary in complex ways because ions with different values of Q/A resonate with different parts of the wave spectrum, which varies with space and time. Just recently, we have begun to model these systematic variations theoretically and to explore other consequences of proton-generated waves.

Description: Computational Fluid Dynamics (CFD) analysis results are compared with benchmark quality test data from the Propulsion Engineering Research Center's (PERC) Rocket Based Combined Cycle (RBCC) experiments to verify fluid dynamic code and application procedures. RBCC engine flowpath development will rely on CFD applications to capture the multi-dimensional fluid dynamic interactions and to quantify their effect on the RBCC system performance. Therefore, the accuracy of these CFD codes must be determined through detailed comparisons with test data. The PERC experiments build upon the well-known 1968 rocket-ejector experiments of Odegaard and Stroup by employing advanced optical and laser based diagnostics to evaluate mixing and secondary combustion. The Finite Difference Navier Stokes (FDNS) code was used to model the fluid dynamics of the PERC RBCC ejector mode configuration. Analyses were performed for both Diffusion and Afterburning (DAB) and Simultaneous Mixing and Combustion (SMC) test conditions. Results from both the 2D and the 3D models are presented.

Description: In this paper we describe the numerical analysis underlying our efforts to develop an accurate and reliable code for simulating flame propagation using complex physical and chemical models. We discuss our spatial and temporal discretization schemes, which in our current implementations range in order from two to six. In space we use staggered meshes to define discrete divergence and gradient operators, allowing us to approximate complex diffusion operators while maintaining ellipticity. Our temporal discretization is based on the use of preconditioning to produce a highly efficient linearly implicit method with good stability properties. High order for time accurate simulations is obtained through the use of extrapolation or deferred correction procedures. We also discuss our techniques for computing stationary flames. The primary issue here is the automatic generation of initial approximations for the application of Newton's method. We use a novel time-stepping procedure, which allows the dynamic updating of the flame speed and forces the flame front towards a specified location. Numerical experiments are presented, primarily for the stationary flame problem. These illustrate the reliability of our techniques, and the dependence of the results on various code parameters.

Description: We present an Abel-inversion algorithm to reconstruct mean and rms refractive-index profiles from spatially resolved statistical measurements of the beam-deflection angle in time-dependent, axisymmetric flows. An oscillating gas-jet diffusion flame was investigated as a test case for applying the algorithm. Experimental data were obtained across the whole field by a rainbow schlieren apparatus. Results show that simultaneous multipoint measurements are necessary to reconstruct the rms refractive index accurately.

Description: A detailed analysis of two of the dynamic maneuvers, the pushover and elevator doublet, from the NASA/FAA Tailplane Icing Program are discussed. For this series of flight tests, artificial ice shapes were attached to the leading edge of the horizontal stabilizer of the NASA Lewis Research Center icing aircraft, a DHC-6 Twin Otter. The purpose of these tests was to learn more about ice-contaminated tailplane stall (ICTS), the known cause of 16 accidents resulting in 139 fatalities. The pushover has been employed by the FAA, JAA and Transport Canada for tailplane icing certification. This research analyzes the pushover and reports on the maneuver performance degradation due to ice shape severity and flap deflection. A repeatability analysis suggests tolerances for meeting the required targets of the maneuver. A second maneuver, the elevator doublet, is also studied.