ALBERT

Description: This paper presents turbulence modeling from NASA's perspective. The topics include: 1) Hierarchy of Solution Methods; 2) Turbulence Modeling Focus; 3) Linear Eddy Viscosity Models; and 4) Nonlinear Eddy Viscosity Algebraic Stress Models.

Description: This paper presents the purpose and expectations of turbulence modeling, in viewgraph form.

Description: The value of the use of the Reynolds-averaged Navier-Stokes methodology for active flow control applications is assessed. An experimental flow control database exists for a NACA0015 airfoil modified at the leading edge to implement a fluidic actuator; hence, this configuration is used. Computational results are documented for the baseline wing configuration (no control) with the experimental results and assumes two-dimensional flow. The baseline wing configuration has discontinuities at the leading edge, trailing edge, and aft of midchord on the upper surface. A limited number of active flow control applications have been tested in the laboratory and in flight. These applications include dynamic stall control using a deformable leading edge, separation control for takeoff and landing flight conditions using piezoelectric devices, pulsed vortex generators, zero-net-mass oscillations, and thrust vectoring with zero-net-mass piezoelectric-driven oscillatory actuation. As yet, there is no definitive comparison with experimental data that indicates current computational capabilities can quantitatively predict the large aerodynamic performance gains achieved with active flow control in the laboratory. However, one study using the Reynolds-averaged Navier-Stokes (RANS) methodology has shown good quantitative agreement with experimental results for an isolated zero-net-mass actuator. In addition, some recent studies have used RANS to demonstrate qualitative performance gains compared with the experimental data for separation control on an airfoil. Those quantitative comparisons for both baseline and flow control cases indicated that computational results were in poor quantitative agreement with the experiments. The current research thrust will investigate the potential use of an unstructured grid RANS approach to predict aerodynamic performance for active flow control applications building on the early studies. First the computational results must quantitatively match experiments for the no-control case before proceeding to the time-dependent flow control case. This paper documents the baseline (no-control) case using an unswept airfoil configuration. The next section describes the configurations used for the computations and the experimentals. The computational approach is then described followed by results and concluding remarks.

Description: Two numerical procedures, one based on artificial compressibility method and the other pressure projection method, are outlined for obtaining time-accurate solutions of the incompressible Navier-Stokes equations. The performance of the two method are compared by obtaining unsteady solutions for the evolution of twin vortices behind a at plate. Calculated results are compared with experimental and other numerical results. For an un- steady ow which requires small physical time step, pressure projection method was found to be computationally efficient since it does not require any subiterations procedure. It was observed that the artificial compressibility method requires a fast convergence scheme at each physical time step in order to satisfy incompressibility condition. This was obtained by using a GMRES-ILU(0) solver in our computations. When a line-relaxation scheme was used, the time accuracy was degraded and time-accurate computations became very expensive.

Description: A parallelized version of the Flowfield Dependent Variation (FDV) Method is developed to analyze a problem of current research interest, the flowfield resulting from a triple shock/boundary layer interaction. Such flowfields are often encountered in the inlets of high speed air-breathing vehicles including the NASA Hyper-X research vehicle. In order to resolve the complex shock structure and to provide adequate resolution for boundary layer computations of the convective heat transfer from surfaces inside the inlet, models containing over 500,000 nodes are needed. Efficient parallelization of the computation is essential to achieving results in a timely manner. Results from a parallelization scheme, based upon multi-threading, as implemented on multiple processor supercomputers and workstations is presented.

Description: The capability to vent in zero gravity without resettling is a technology need that involves practically all uses of sub-critical cryogenics in space, and would extend cryogenic orbital transfer vehicle capabilities. However, the lack of definition regarding liquid/ullage orientation coupled with the somewhat random nature of the thermal stratification and resulting pressure rise rates, lead to significant technical challenges. Typically a zero gravity vent concept, termed a thermodynamic vent system (TVS), consists of a tank mixer to destratify the propellant, combined with a Joule-Thomson (J-T) valve to extract thermal energy from the propellant. Marshall Space Flight Center's (MSFC's) Multipurpose Hydrogen Test Bed (MHTB) was used to test both spray-bar and axial jet TVS concepts. The axial jet system consists of a recirculation pump heat exchanger unit. The spray-bar system consists of a recirculation pump, a parallel flow concentric tube heat exchanger, and a spray-bar positioned close to the longitudinal axis of the tank. The operation of both concepts is similar. In the mixing mode, the recirculation pump withdraws liquid from the tank and sprays it into the tank liquid, ullage, and exposed tank surfaces. When energy extraction is required, a small portion of the recirculated liquid is passed sequentially through the J-T expansion valve, the heat exchanger, and is vented overboard. The vented vapor cools the circulated bulk fluid, thereby removing thermal energy and reducing tank pressure. The pump operates alone, cycling on and off, to destratify the tank liquid and ullage until the liquid vapor pressure reaches the lower set point. At that point, the J-T valve begins to cycle on and off with the pump. Thus, for short duration missions, only the mixer may operate, thus minimizing or even eliminating boil-off losses.

Description: Here we present the results of a Large Eddy Simulation of a non-buoyant jet issuing from a circular orifice in a wall, and developing in neutral surroundings. The effects of the subgrid scales on the large eddies have been modeled with the dynamic large eddy simulation model applied to the fully 3D domain in spherical coordinates. The simulation captures the unsteady motions of the large-scales within the jet as well as the laminar motions in the entrainment region surrounding the jet. The computed time-averaged statistics (mean velocity, concentration, and turbulence parameters) compare well with laboratory data without invoking an empirical entrainment coefficient as employed by line integral models. The use of the large eddy simulation technique allows examination of unsteady and inhomogeneous features such as the evolution of eddies and the details of the entrainment process.

Description: Modern high-work turbines can be compact, transonic, supersonic, counter-rotating, or use a dense drive gas. The vast majority of modern rocket turbine designs fall into these Categories. These turbines usually have large temperature variations across a given stage, and are characterized by large amounts of flow unsteadiness. The flow unsteadiness can have a major impact on the turbine performance and durability. For example, the Space Transportation Main Engine (STME) fuel turbine, a high work, transonic design, was found to have an unsteady inter-row shock which reduced efficiency by 2 points and increased dynamic loading by 24 percent. The Revolutionary Reusable Technology Turbopump (RRTT), which uses full flow oxygen for its drive gas, was found to shed vortices with such energy as to raise serious blade durability concerns. In both cases, the sources of the problems were uncovered (before turbopump testing) with the application of validated, unsteady computational fluid dynamics (CFD) to the designs. In the case of the RRTT and the Alternate Turbopump Development (ATD) turbines, the unsteady CFD codes have been used not just to identify problems, but to guide designs which mitigate problems due to unsteadiness. Using unsteady flow analyses as a part of the design process has led to turbine designs with higher performance (which affects temperature and mass flow rate) and fewer dynamics problems. One of the many assumptions made during the design and analysis of supersonic turbine stages is that the values of the specific heats are constant. In some analyses the value is based on an average of the expected upstream and downstream temperatures. In stages where the temperature can vary by 300 to 500 K, however, the assumption of constant fluid properties may lead to erroneous performance and durability predictions. In this study the suitability of assuming constant specific heats has been investigated by performing three-dimensional unsteady Navier-Stokes simulations for a supersonic turbine stage.

Description: B-spline functions are bases for piecewise polynomials that possess attractive properties for complex flow simulations : they have compact support, provide a straightforward handling of boundary conditions and grid nonuniformities, and yield numerical schemes with high resolving power, where the order of accuracy is a mere input parameter. This paper reviews the progress made on the development and application of B-spline numerical methods to computational fluid dynamics problems. Basic B-spline approximation properties is investigated, and their relationship with conventional numerical methods is reviewed. Some fundamental developments towards efficient complex geometry spline methods are covered, such as local interpolation methods, fast solution algorithms on cartesian grid, non-conformal block-structured discretization, formulation of spline bases of higher continuity over triangulation, and treatment of pressure oscillations in Navier-Stokes equations. Application of some of these techniques to the computation of viscous incompressible flows is presented.

Description: The importance of reducing jet noise in both commercial and military aircraft applications has made jet acoustics a significant area of research. A technique for jet noise prediction commonly employed in practice is the MGB approach, based on the Lighthill acoustic analogy. This technique requires as aerodynamic input mean flow quantities and turbulence quantities like the kinetic energy and the dissipation. The purpose of the present paper is to assess existing capabilities for predicting these aerodynamic inputs. Two modern Navier-Stokes flow solvers, coupled with several modern turbulence models, are evaluated by comparison with experiment for their ability to predict mean flow properties in a supersonic jet plume. Potential weaknesses are identified for further investigation. Another comparison with similar intent is discussed by Barber et al. The ultimate goal of this research is to develop a reliable flow solver applicable to the low-noise, propulsion-efficient, nozzle exhaust systems being developed in NASA focused programs. These programs address a broad range of complex nozzle geometries operating in high temperature, compressible, flows. Seiner et al. previously discussed the jet configuration examined here. This convergent-divergent nozzle with an exit diameter of 3.6 inches was designed for an exhaust Mach number of 2.0 and a total temperature of 1680 F. The acoustic and aerodynamic data reported by Seiner et al. covered a range of jet total temperatures from 104 F to 2200 F at the fully-expanded nozzle pressure ratio. The aerodynamic data included centerline mean velocity and total temperature profiles. Computations were performed independently with two computational fluid dynamics (CFD) codes, ISAAC and PAB3D. Turbulence models employed include the k-epsilon model, the Gatski-Speziale algebraic-stress model and the Girimaji model, with and without the Sarkar compressibility correction. Centerline values of mean velocity and mean temperature are compared with experimental data.

Description: Doppler global velocimetry uses the absorption characteristics of iodine vapor to provide instantaneous three-component measurements of flow velocity within a plane defined by a laser light sheet. Although the technology is straightforward, its utilization as a flow diagnostics tool requires hardening of the optical system and careful attention to detail during data acquisition and processing if routine use in wind tunnel applications is to be achieved. A development program that reaches these goals is presented. Theoretical and experimental investigations were conducted on each technology element to determine methods that increase measurement accuracy and repeatability. Enhancements resulting from these investigations included methods to ensure iodine vapor calibration stability, single frequency operation of the laser and image alignment to sub-pixel accuracies. Methods were also developed to improve system calibration, and eliminate spatial variations of optical frequency in the laser output, spatial variations in optical transmissivity and perspective and optical distortions in the data images. Each of these enhancements is described and experimental examples given to illustrate the improved measurement performance obtained by the enhancement. The culmination of this investigation was the measured velocity profile of a rotating wheel resulting in a 1.75% error in the mean with a standard deviation of 0.5 m/s. Comparing measurements of a jet flow with corresponding Pitot measurements validated the use of these methods for flow field applications.

Description: The performance of an explicit algebraic stress model (EASM) is assessed in predicting the turbulent flow and forced heat transfer in both straight and wavy ducts, with rectangular, trapezoidal and triangular cross-sections, under fully developed conditions. A comparison of secondary flow patterns. including velocity vectors and velocity and temperature contours, are shown in order to study the effect of waviness on flow dynamics, and comparisons between the hydraulic parameters. Fanning friction factor and Nusselt number, are also presented. In all cases. isothermal conditions are imposed on the duct walls, and the turbulent heat fluxes are modeled using gradient-diffusion type models. The formulation is valid for Reynolds numbers up to 10(exp 5) and this minimizes the need for wall functions that have been used with mixed success in previous studies of complex duct flows. In addition, the present formulation imposes minimal demand on the number of grid points without any convergence or stability problems. Criteria in terms of heat transfer and friction factor needed to choose the optimal wavy duct cross-section for industrial applications among the ones considered are discussed.

Description: A compliant, thermal interface material is tested to evaluate its thermal behavior at elevated temperatures, in vacuum conditions, and under varying levels of compression. Preliminary results indicate that the thermal performance of this polymer fiber-based, felt-like material is sufficient to meet thermal extraction requirements for the Quench Module Insert, a Bridgman furnace for microgravity material science investigation. This paper discusses testing and modeling approaches employed, gives of a status of characterization activities and provides preliminary test results.

Description: Experimental studies of supersonic impinging jet flows suggest that they are greatly influenced by the flow-acoustic interactions through a feedback mechanism. The self-sustained oscillations of the jet column observed in these flows result in high velocities in the ambient medium induced by the large-scale coherent vortical structures in the jet shear layers. As a consequence, the suck down force on the surface from which the jet is issuing can reach as high as 60% of the primary jet thrust. In addition, the overall sound pressure levels (OASPL) increase significantly relative to a free jet. To alleviate these undesirable flow and acoustic characteristics, a novel control technique using supersonic microjets is demonstrated. Sixteen supersonic microjets are placed around the circumference of the main jet at the nozzle exit to disrupt the feedback mechanism. As a result, significant lift loss recovery (approximately 50%) and reduced near field OASPL (approximately 7 dB) are observed.

Description: Tension in Microgravity Modeled Microgravity poses many challenges to the designer of spacecraft tanks. Chief among these are the lack of phase separation and the need to supply vapor-free liquid or liquidfree vapor to the spacecraft processes that require fluid. One of the principal problems of phase separation is the creation of liquid jets. A jet can be created by liquid filling, settling of the fluid to one end of the tank, or even closing a valve to stop the liquid flow. Anyone who has seen a fountain knows that jets occur in normal gravity also. However, in normal gravity, the gravity controls and restricts the jet flow. In microgravity, with gravity largely absent, jets must be contained by surface tension forces. Recent NASA experiments in microgravity (Tank Pressure Control Experiment, TPCE, and Vented Tank Pressure Experiment, VTRE) resulted in a wealth of data about jet behavior in microgravity. VTRE was surprising in that, although it contained a complex geometry of baffles and vanes, the limit on liquid inflow was the emergence of a liquid jet from the top of the vane structure. Clearly understanding the restraint of liquid jets by surface tension is key to managing fluids in low gravity. To model this phenomenon, we need a numerical method that can track the fluid motion and the surface tension forces. The fluid motion is modeled with the Navier-Stokes equation formulated for low-speed incompressible flows. The quantities of velocity and pressure are placed on a staggered grid, with velocity being tracked at cell faces and pressure at cell centers. The free surface is tracked via the introduction of a color function that tracks liquid as 1/2 and gas as -1/2. A phase model developed by Jacqmin is used. This model converts the discrete surface tension force into a barrier function that peaks at the free surface and decays rapidly. Previous attempts at this formulation have been criticized for smearing the interface. However, by sharpening the phase function, double gridding the fluid function, and using a higher order solution for the fluid function, interface smearing is avoided. These equations can be rewritten as two coupled Poisson equations that also include the velocity. The method of solution is as follows: first, the phase equations are solved from this solution, a velocity field is generated, then a successive overrelaxation scheme is used to solve for a pressure field consistent with the velocity solution. After the code was implemented in axisymmetric form and verified by several test cases, the drop tower runs of Aydelott were modeled. The model handed the free-surface deformation quite nicely, even to the point of modeling geyser growth in the regime where the free surface was no longer restrained. A representative run is shown.

Description: Of vital interest to aerodynamic researchers is precise knowledge of the flow velocity profile next to the surface. This information is needed for turbulence model development and the calculation of viscous shear force. Though many instruments can determine the flow velocity profile near the surface, none of them can make measurements closer than approximately 0.01 in. from the surface. The thermocouple boundary-layer rake can measure much closer to the surface than conventional instruments can, such as a total pressure boundary layer rake, hot wire, or hot film. By embedding the sensors (thermocouples) in the region where the velocity is equivalent to the velocity ahead of a constant thickness strut, the boundary-layer flow profile can be obtained. The present device fabricated at the NASA Glenn Research Center microsystem clean room has a heater made of platinum and thermocouples made of platinum and gold. Equal numbers of thermocouples are placed both upstream and downstream of the heater, so that the voltage generated by each pair at the same distance from the surface is indicative of the difference in temperature between the upstream and downstream thermocouple locations. This voltage differential is a function of the flow velocity, and like the conventional total pressure rake, it can provide the velocity profile. In order to measure flow extremely close to the surface, the strut is made of fused quartz with extremely low heat conductivity. A large size thermocouple boundary layer rake is shown in the following photo. The latest medium size sensors already provide smooth velocity profiles well into the boundary layer, as close as 0.0025 in. from the surface. This is about 4 times closer to the surface than the previously used total pressure rakes. This device also has the advantage of providing the flow profile of separated flow and also it is possible to measure simultaneous turbulence levels within the boundary layer.

Description: Aircraft compressors can suffer debilitating consequences as a result of rotating stall and surge events caused by inlet distortions. This is particularly true of aircraft during takeoff, when the compressor is operating at peak performance close to the surge line. Significant research has been conducted by the NASA Glenn Research Center in the area of compressor stability enhancement through active and passive control methods. Most recently, an experiment was conducted at the Wright Patterson Air Force Base Research Laboratory on a two-stage fan with inlet guide vanes and inlet distortion. In this joint Small Business Innovation Research effort between Scientific Systems and Glenn, control of rotating stall was demonstrated in a multistage transonic fan. This twostage fan with inlet guide vanes was tested under clean and distorted inlet conditions. The compressor was also configured with a circumferential distortion screen capable of 180 of distortion and with 14 high-velocity injectors upstream of the first rotor. Twelve of these injectors could oscillate up to frequencies of 450 Hz. The additional two injectors were located next to each other and were used in concert with each other as a single, on/off, high-authority actuator. In a first test of injection in this multistage environment, 12 of the valves were opened 50 percent of their full stroke to assess steady injection through the compressor. This baseline injection is shown in the compressor characteristic of the following figure, and stall margin improvements are tracked from this baseline condition. The compressor was then tested with clean inlet conditions using 12 injectors and active control. Pressure disturbances were tracked before rotating stall, and a constant gain control scheme reduced the stalling mass flow by 10.8 percent over the baseline. With the distortion screen present in the inlet, a pole-zero cancellation control scheme was used to achieve a 6.4-percent decrease in stalling mass flow. These improvements also are shown in the figure. In a final experiment, actively controlled, high-frequency injection from the 12 valves was used in conjunction with the high-authority actuators. In this test, the stalling mass flow of the compressor was reduced by 27 percent as indicated in the graph. These results were obtained by injecting less than 2 percent of the total compressor throughflow into the rotor tip region via 14 injection ports. These results mark the first successful demonstration of actively controlled air injection as a stall-control strategy for multistage compressors operating at speeds typical of an actual gas turbine engine. A goal of continuing research is to determine the combination of air-injection parameters and control strategies that are most effective in providing stall control for both clean and distorted inlet flow conditions for multistage environments. Other goals include the demonstration of stall control at many locations along the core compressor and development and application of active stall control strategies that will be integral flightworthy components of onboard engine hardware.

Description: A contact line is defined at the intersection of a solid surface with the interface between two immiscible fluids. When one fluid displaces another immiscible fluid along a solid surface, the process is called dynamic wetting and a "moving" contact line (one whose position relative to the solid changes in time) often appears. The physics of dynamic wetting controls such natural and industrial processes as spraying of paints and insecticides, dishwashing, film formation and rupture in the eye and in the alveoli, application of coatings, printing, drying and imbibition of fibrous materials, oil recovery from porous rocks, and microfluidics.

Description: The simultaneous flow of gas and liquid through a fixed bed of particles occurs in many unit operations of interest to the designers of space-based as well as terrestrial equipment. Examples include separation columns, gas-liquid reactors, humidification, drying, extraction, and leaching. These operations are critical to a wide variety of industries such as petroleum, pharmaceutical, mining, biological, and chemical. NASA recognizes that similar operations will need to be performed in space and on planetary bodies such as Mars if we are to achieve our goals of human exploration and the development of space. The goal of this research is to understand how to apply our current understanding of two-phase fluid flow through fixed-bed reactors to zero- or partial-gravity environments. Previous experiments by NASA have shown that reactors designed to work on Earth do not necessarily function in a similar manner in space. Two experiments, the Water Processor Assembly and the Volatile Removal Assembly have encountered difficulties in predicting and controlling the distribution of the phases (a crucial element in the operation of this type of reactor) as well as the overall pressure drop.

Description: The dynamics of miscible displacements in a cylindrical tube are being investigated experimentally and numerically, with a view to understand the complex processes that occur, for example, in enhanced oil recovery, hydrology, and filtration. We have observed complex shapes of the interface between two liquids that mix with each other when the less viscous liquid is displaced by the more viscous one in a tube. A less viscous fluid that displaces a more viscous fluid is known to propagate in the form of a "finger," and a flight experiment proposed by Maxworthy et al. to investigate the miscible-interface dynamics is currently being developed by NASA. From the current theory of miscible displacements, which was developed for a porous medium satisfying Darcy's law, it can be shown that in the absence of gravity the interface between the fluids is destabilized and thus susceptible to fingering only when a more viscous fluid is displaced by a less viscous one. Therefore, if the interface is initially flat and the more viscous fluid displaces the less viscous fluid, the interface ought to be stable and remain flat. However, numerical simulations by Chen and Meiburg for such displacement in a cylindrical tube show that the interface is unstable and a finger of the more viscous fluid is indeed formed. Preliminary experiments performed at the NASA Glenn Research Center show that not only can fingering occur when the more viscous fluid displaces a less viscous one in a cylindrical tube, but also that under certain conditions the advancing finger achieves a sinuous or snakelike shape. These experiments were performed using silicone oils in a vertical pipette of small diameter. In the initial configuration, the more viscous fluid rested on top of the less viscous one, and the interface was nominally flat. A dye was added to the upper liquid for ease of observation of the interface between the fluids. The flow was initiated by draining the lower fluid from the bottom of the pipette, at speeds less than 0.1 mm/sec.

Description: Significant research has been underway for several years in NASA Glenn Research Center's nozzle branch to develop advanced computational methods for simulating turbulent flows in exhaust nozzles. The primary efforts of this research have concentrated on improving our ability to calculate the turbulent mixing layers that dominate flows both in the exhaust systems of modern-day aircraft and in those of hypersonic vehicles under development. As part of these efforts, a hybrid numerical method was recently developed to simulate such turbulent mixing layers. The method developed here is intended for configurations in which a dominant structural feature provides an unsteady mechanism to drive the turbulent development in the mixing layer. Interest in Large Eddy Simulation (LES) methods have increased in recent years, but applying an LES method to calculate the wide range of turbulent scales from small eddies in the wall-bounded regions to large eddies in the mixing region is not yet possible with current computers. As a result, the hybrid method developed here uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall-bounded regions entering a mixing section and uses a LES procedure to calculate the mixing-dominated regions. A numerical technique was developed to enable the use of the hybrid RANS-LES method on stretched, non-Cartesian grids. With this technique, closure for the RANS equations is obtained by using the Cebeci-Smith algebraic turbulence model in conjunction with the wall-function approach of Ota and Goldberg. The LES equations are closed using the Smagorinsky subgrid scale model. Although the function of the Cebeci-Smith model to replace all of the turbulent stresses is quite different from that of the Smagorinsky subgrid model, which only replaces the small subgrid turbulent stresses, both are eddy viscosity models and both are derived at least in part from mixing-length theory. The similar formulation of these two models enables the RANS and LES equations to be solved with a single solution scheme and computational grid. The hybrid RANS-LES method has been applied to a benchmark compressible mixing layer experiment in which two isolated supersonic streams, separated by a splitter plate, provide the flows to a constant-area mixing section. Although the configuration is largely two dimensional in nature, three-dimensional calculations were found to be necessary to enable disturbances to develop in three spatial directions and to transition to turbulence. The flow in the initial part of the mixing section consists of a periodic vortex shedding downstream of the splitter plate trailing edge. This organized vortex shedding then rapidly transitions to a turbulent structure, which is very similar to the flow development observed in the experiments. Although the qualitative nature of the large-scale turbulent development in the entire mixing section is captured well by the LES part of the current hybrid method, further efforts are planned to directly calculate a greater portion of the turbulence spectrum and to limit the subgrid scale modeling to only the very small scales. This will be accomplished by the use of higher accuracy solution schemes and more powerful computers, measured both in speed and memory capabilities.

Description: Modern turbofan engines employ a highly loaded fan stage with transonic or low-supersonic velocities in the blade-tip region. The fan blades are often prone to flutter at off-design conditions. Flutter is a highly undesirable and dangerous self-excited mode of blade oscillations that can result in high-cycle fatigue blade failure. The origins of blade flutter are not fully understood yet. The latest view is that the blade oscillations are triggered by high-frequency changes in the extent of the partially separated area on the airfoil suction side. There is a lack of experimental data describing the separated flow characteristics of modern airfoils for transonic fans.

Description: Loss of lubrication in rotorcraft drive systems is a demanding requirement placed on drive system manufacturers. The drive system must operate for at least 30 minutes once the primary lubrication system has failed. This test is a military requirement that must be passed prior to certification of the aircraft. As new aircraft engines, operating at higher speeds, are fielded, the requirements for the drive system become increasingly more difficult. Also, the drive system must be lightweight, which minimizes the opportunity to use the gear bodies to absorb the tremendous amount of heating that takes place. In many cases, the amount of heat generated because of the high speed and load requires an emergency lubrication system that negatively impacts the aircraft's weight, complexity, and cost. A single mesh spur gear test rig is being used at the NASA Glenn Research Center to investigate possible emergency lubrication system improvements that will minimize the impact of having these systems onboard rotorcraft. A technique currently being investigated uses a vapor/mist system to lubricate the contacting surfaces after the primary lubrication system has been shut down. A number of tests were conducted in which the vapor/mist used the same lubricant as the primary system, but at a greatly reduced flow rate. Each test was initiated with the primary lubrication system operational and at steady-state conditions for a given speed and load. Then the primary lubrication system was shut down, and the vapor/mist lubrication system was initiated. An example of the tests conducted is shown in the figures. These preliminary tests have uncovered a mechanism that provides a lubricious, carbonaceous solid on the surface that actually reduces the surface temperature of the meshing gear teeth during operation. Surface analysis of the carbonaceous solid revealed it was graphitic. This mechanism is the synthetic lubricant "coking" on the active profile of the gears, which reduces the friction between the contacting gear surfaces. The level of load affects the onset of this mechanism: the higher the load, the sooner coking takes place. Future work will investigate several other factors that could improve the already promising results that have been attained.

Description: For presentation at Northwestern University, May 14-28, 2001, the talk will present two types of phenomena, both recognizable to students of nonlinearity, that are exhibited by axisymmetric vortex rings in numerical and laboratory experiments. (1) The first type of phenomenon is reminiscent of inelastic solitons.

Description: Particle Image Velocimetry (PIV) is used to study the slow settling motions of spheres in suspensions ranging from dilute to highly concentrated, 0.05 less than phi less than 0.50. During sedimentation, velocity fluctuations are found to be organized into regions of characteristic size L approximately 11 a phi(exp-1/3) and buoyant mass DELTA m given by the rms density fluctuations in a region of size L(sup 3) assuming random statistics with excluded volume effects. A simple model accurately predicts the magnitudes of the velocity fluctuations DELTA V. Finally, we find a new universal relation for particle diffusion during sedimentation that can be written in a Stokes-Einstein type form as D approximately (DELTA mg L)/(6 pi eta L), where the effective temperature DELTA mg L is the gravitational potential energy of density fluctuations.

Description: The NASA Reaction Control System (RCS) Plume Model (RPM) is an exhaust plume flow field and impingement heating code that has been updated and applied to components of the International Space Station (ISS). The objective of this study was to use this code to determine if plume environments from either Orbiter PRCS jets or ISS reboost and Attitude Control System (ACS) jets cause thermal issues on ISS component surfaces. This impingement analysis becomes increasingly important as the ISS is being assembled with its first permanent crew scheduled to arrive by the end of fall 2000. By early summer 2001 , the ISS will have a number of major components installed such as the Unity (Node 1), Destiny (Lab Module), Zarya (Functional Cargo Block), and Zvezda (Service Module) along with the P6 solar arrays and radiators and the Z-1 truss. Plume heating to these components has been analyzed with the RPM code as well as additional components for missions beyond Flight 6A such as the Propulsion Module (PM), Mobile Servicing System, Space Station Remote Manipulator System, Node 2, and the Cupola. For the past several years NASA/JSC has been developing the methodology to predict plume heating on ISS components. The RPM code is a modified source flow code with capabilities for scarfed nozzles and intersecting plumes that was developed for the 44 Orbiter RCS jets. This code has been validated by comparison with Shuttle Plume Impingement Flight Experiment (SPIFEX) heat flux and pressure data and with CFD and Method of Characteristics solutions. Previous analyses of plume heating predictions to the ISS using RPM have been reported, but did not consider thermal analysis for the components nor jet-firing histories as the Orbiter approaches the ISS docking ports. The RPM code has since been modified to analyze surface temperatures with a lumped mass approach and also uses jet-firing histories to produce pulsed heating rates. In addition, RPM was modified to include plume heating from ISS jets to ISS components where the jet coordinates are specified, together with the engine cant angle. These latter studies have been focused on the PM with plumes from its reboost and ACS jets impinging on various ISS components and also focused on the Japanese H2 Transfer Vehicle (HTV) with the plumes from its reboost engines impinging on the Cupola window. This paper will present plume heating and surface temperature results on a number of ISS components with and without jet-firing histories, evaluate post-flight data, and describe any potential thermal issues

Description: The NPARC Alliance is a cooperative effort between the United States Air Force's Arnold Engineering Development Center, NASA's Glenn Research Center and The Boeing Company. The mission of the Alliance is to develop, validate, and support an integrated, general purpose computational flow simulator for the U.S. aerospace community. The Alliance provides a state of the art simulation system that includes geometry manipulation, flow solution, and post-processing capabilities. The system is centered around the WIND flow solver. This presentation provides an overview of the Alliance and the flowfield simulation system. Several example computations are provided.

Description: Recent experiments have shown that low-density gas jets injected into a high-density gas undergo an instability mode leading to highly periodic oscillations in the flow field. The transition from laminar to turbulent flow in these jets is abrupt, without a gradual change in scales. Although this type of instability at high Richardson numbers has been attributed to buoyancy, direct physical evidence was not acquired through experiments. In this study, several experiments were conducted in Earth gravity and microgravity to acquire qualitative data on near field flow structure of helium jets injected into air. Microgravity conditions were simulated in the 2.2-second drop tower at NASA Glenn Research Center. The operating parameters of this study included the tube inside diameter, the jet Reynolds number, and the jet Richardson number. Tubes with inside diameters of 19.05 mm and 31.75 mm were used in the experiments conducted in the drop tower. The jet flow was analyzed using quantitative rainbow schlieren deflectometry, a non-intrusive line of sight measurement technique for the whole field. The flow structure was characterized by distributions of angular deflection and the resulting helium mole fraction obtained from color schlieren images taken at 60 Hz. Three sets of experimental data with respect to three schlieren fields of view were acquired for each tube. Results show that the jet in microgravity was up to 70 percent wider than that in Earth gravity. The global jet flow oscillations observed in Earth gravity were absent in microgravity, providing direct experimental evidence that the flow instability in the low-density jet was buoyancy-induced. This study provides quantitative details of temporal flow evolution as the experiments undergo change in gravity in the drop tower.

Description: The objective of this study was to determine the effects of buoyancy on the absolute instability of low-density gas jets injected into high-density gas mediums. Most of the existing analyses of low-density gas jets injected into a high-density ambient have been carried out neglecting effects of gravity. In order to investigate the influence of gravity on the near-injector development of the flow, a linear temporal stability analysis and a spatio-temporal stability analysis of a low-density round jet injected into a high-density ambient gas were performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The temporal growth rates and the phase velocity of the disturbances were obtained. It was found that the presence of variable density within the shear layer resulted in an increase in the temporal amplification rate of the disturbances and an increase in the range of unstable frequencies, accompanied by a reduction in the phase velocities of the disturbances. Also, the temporal growth rates of the disturbances were increased as the Froude number was reduced (i.e. gravitational effects increased), indicating the destabilizing role played by gravity. The spatio-temporal stability analysis was performed to determine the nature of the absolute instability of the jet. The roles of the density ratio, Froude number, Schmidt number, and the lateral shift between the density and velocity profiles on the jet s absolute instability were determined. Comparisons of the results with previous experimental studies show good agreement when the effects of these variables are combined together. Thus, the combination of these variables determines how absolutely unstable the jet will be. Experiments were carried out to observe the qualitative differences between a round low-density gas jet injected into a high-density gas (helium jet injected into air) and a round constant density jet (air jet injected into air). Flow visualizations and velocity measurements in the near-injector region of the helium jet show more mixing and spreading of the helium jet than the air jet. The vortex structures develop and contribute to the jet spreading causing the helium jet to oscillate.

Description: The dynamics of an ensemble of linear disturbances with a known probability distribution associated with the initial mode amplitudes are studied in boundary-layer flows through an analysis of the transport equations for the mean disturbance kinetic energy and disturbance energy dissipation rate. Effects of adverse and favorable pressure-gradients on the disturbance dynamics are also included in the analysis. Unlike the fully turbulent regime where nonlinear phase scrambling of the fluctuations affects the flow field in proximity to the wall, the laminar regime fluctuations studied here are influenced across the boundary layer by the solid boundary. In addition to the low Reynolds number, early stage transition regime, the dynamics of these disturbance fields can be related in some respects to the near-wall dynamics of the fully turbulent regime.

Description: Microgravity technologies often require aqueous phases to spread over nonwetting hydrophobic solid surfaces. Surfactants facilitate the wetting of water on hydrophobic surfaces by adsorbing on the water/air and hydrophobic solid/water interfaces and lowering the surface tensions of these interfaces. The tension reductions decrease the contact angle, which increases the equilibrium wetted area. Hydrocarbon surfactants; (i.e., amphiphiles with a hydrophobic moiety consisting of an extended chain of (aliphatic) methylene -CH2- groups attached to a large polar group to give aqueous solubility) are capable of reducing the contact angles on surfaces which are not very hydrophobic, but do not reduce significantly the contact angles of the very hydrophobic surfaces such as parafilm, polyethylene or self assembled monolayers. Trisiloxane surfactants (amphiphiles with a hydrophobe consisting of methyl groups linked to a trisiloxane backbone in the form of a disk ((CH3)3-Si-O-Si-O-Si(CH3)3) and an extended ethoxylate (-(OCH2CH2)a-) polar group in the form of a chain with four or eight units) can significantly reduce the contact angle of water on a very hydrophobic surface and cause rapid and complete (or nearly complete) spreading (termed superspreading). The overall goal of the research described in this proposal is to establish and verify a theory for how trisiloxanes cause superspreading, and then use this knowledge as a guide to developing more general hydrocarbon based surfactant systems which superspread. We propose that the trisiloxane surfactants superspread because their structure allows them to strongly lower the high hydrophobic solid/aqueous tension when they adsorb to the solid surface. When the siloxane adsorbs, the hydrophobic disk parts of the molecule adsorb onto the surface removing the surface water. Since the cross-sectional area of the disk is larger than that of the extended ethoxylate chain, the disks can form a space-filling mat on the surface which removes a significant amount of the surface water. In this presentation, we report the results of measurements of the molecular packing and rates of kinetic exchange of the trisiloxane surfactants at the air/water interface in order to confirm our picture of trisiloxane packing, and provide additional insight into the superspreading process. We used the pendant bubble technique as a Langmuir trough to measure the trisiloxane equation of state which relates the tension to the surface concentration. From these measurements we obtain accurate values for the maximum packing density. We find that trisiloxanes with 4 and 8 ethoxylate groups have the same maximum packing concentration, indicating that the maximum packing is controlled by the cross section of the head group. For trisiloxanes with larger than eight ethoxylates, the maximum packing increases with ethoxylate number, indicating that the disposition of the ethoxylate chain (i.e., its effective size) is controlling. This supports our picture of superspreading: The superspreading ability of trisiloxanes decreases considerably for trisiloxanes with larger than eight ethoxylates; the packing measurements indicate that for the higher ethoxylate number trisiloxanes, the compact nonpolar head groups are pushed apart by the ethoxylate chain. They leave spaces of surface water on adsorption and do not lower the solid tension as much as their lower chain analogues. Finally the report measurements of the dynamic tension reduction accompanying the adsorption of trisiloxanes onto an initially clean interface using the pendant bubble technique, and we obtain from these relaxations, the equation of state and a mass transfer model, the rate constants for kinetic exchange. We find that the rate constants for desorption of trisiloxanes are generally much slower than for analogous aliphatic polyethoxylate surfactants with identical ethoxylate chain lengths. When an aqueous drop of a superspreader solution is placed on a hydrophobic surface and begins to spread, the lower desorption rates allows the tension at the drop center to remain reduced relative to the tension of the expanding periphery, thus strengthening Marangoni forces which can assist the spreading. Marangoni forces can be especially significant in the case of superspreaders because their maximum equilibrium reduction in tension is several dynes/cm lower than for aliphatic surfactants.

Description: Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow that propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillary to direct the movement of bubbles in space is the fact that surfactant impurities, which are unavoidably present in the continuous phase, can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. (This can also be the case if kinetic desorption of surfactant at the back end of the bubble is much slower than convection.) For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like and bubbles translate as solid particles). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature-induced tension gradient, and can decrease to near zero the bubble's thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have be retarded by the adsorption onto the bubble surface of a surfactant impurity. Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble surface much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition, the rapid kinetic exchange and high bulk concentration maintain a saturated surface with uniform surface concentrations. This prevents retarding surface tension gradients and keeps the thermocapillary velocity high. In our reports over the first 2 years, we presented numerical simulations of the bubble motion and surfactant transport which verified theoretically the concept of remobilization, and the development of an apparatus to track and measure the velocity of rising bubbles in a glycerol/water surfactant solution. This year, we detail experimental observations of remobilization. Two polyethylene oxide surfactants were studied, C12E6 (CH3(CH2)11(OCH2)6OH) and C10E8 (CH3(CH2)4(OCH2CH2)8OH). Measurements of the kinetic exchange for these surfactants show that the one with the longer hydrophobe chain C12E6 has a lower rate of kinetic exchange. In addition, this surfactant is much less soluble in the glycerol/water mixture because of the shorter ethoxylate chain. As a result, we found that C12E6 had only a very limited ability to remobilize rising bubbles because of the limited kinetic exchange and reduced solubility. However, C10E8, with its higher solubility and more rapid exchange was found to dramatically remobilize rising bubbles. We also compared our theoretical calculations to the experimental measurements of velocity for both the non-remobilizing and remobilizing surfactants and found excellent agreement. We further observed that for C10E8 at high concentrations, which exceeded the critical micelle concentrations, additional remobilization was measured. In this case the rapid exchange of monomer between micelle and surfactant provides an additional mechanism for maintaining a uniform surface concentrations.

Description: Modern Computational Fluid Dynamics (CFD) techniques were used to compute aerodynamic forces and moments of the Space Shuttle Orbiter in specific portions of contingency abort trajectory space. The trajectory space covers a Mach number range of 3.5-15, an angle-of-attack range of 20-60 degrees, an altitude range of 100-190 kft, and several different settings of the control surfaces (elevons, body flap, and speed brake). While approximately 40 cases have been computed, only a sampling of the results is presented here. The computed results, in general, are in good agreement with the Orbiter Operational Aerodynamic Data Book (OADB) data (i.e., within the uncertainty bands) for almost all the cases. However, in a limited number of high angle-of-attack cases (at Mach 15), there are significant differences between the computed results, especially the vehicle pitching moment, and the OADB data. A preliminary analysis of the data from the CFD simulations at Mach 15 shows that these differences can be attributed to real-gas/Mach number effects.

Description: The flow of metal during Friction Stir Welding is clarified using a faying surface tracer and a nib frozen in place during welding. It is shown that material is transported in two distinct streams or currents. One stream is a wiping of material from the advancing front side of the nib onto a plug of material that rotates and advances with the nib. The material undergoes a helical motion within the plug that both rotates and advances with the plug and descends in the wash of the threads on the nib and rises on the outer part of the plug. After one or more rotations, this material is sloughed off the plug in the wake of the tool primarily on the advancing side. The second stream of material is an entrainment of material from the retreating side of the nib that fills in between the sloughed off pieces from the advancing side. These two processes produce material with different mechanical properties and the strength of a weld should depend on the relative importance of the processes.

Description: The surface turbulent fluxes of momentum, latent heat, and sensible heat over global oceans are essential to weather, climate and ocean problems. Wind stress is the major forcing for driving the oceanic circulation, while Evaporation is a key component of hydrological cycle and surface heat budget. We have produced a 7.5-year (July 1987-December 1994) dataset of daily, individual monthly-mean and climatological (1988-94) monthly-mean surface turbulent fluxes over the global oceans from measurements of the Special Sensor Microwave/Imager (SSM/I) on board the US Defense Meteorological Satellite Program F8, F10, and F11 satellites. It has a spatial resolution of 2.0x2.5 latitude-longitude. Daily turbulent fluxes are derived from daily data of SSM/I surface winds and specific humidity, National Centers for Environmental Prediction (NCEP) sea surface temperatures, and European Centre for Medium-Range Weather Forecasts (ECMWF) air-sea temperature differences, using a stability-dependent bulk scheme. The retrieved instantaneous surface air humidity (with a 25-km resolution) IS found to be generally accurate as compared to the collocated radiosonde observations over global oceans. The surface wind speed and specific humidity (latent heat flux) derived from the F10 SSM/I are found to be -encrally smaller (larger) than those retrieved from the F11 SSM/I. The F11 SSM/I appears to have slightly better retrieval accuracy for surface wind speed and humidity as compared to the F10 SSM/I. This difference may be due to the orbital drift of the F10 satellite. The daily wind stresses and latent heat fluxes retrieved from F10 and F11 SSM/Is show useful accuracy as verified against the research quality in si -neasurerrients (IMET buoy, RV Moana Wave, and RV Wecoma) in the western Pacific warm pool during the TOGA COARE Intensive observing period (November 1992-February 1993). The 1988-94 seasonal-mean turbulent fluxes and input variables derived from FS and F11 SSM/Is show reasonable patterns related to seasonal variations of atmospheric general circulation. This dataset of SSM/I-derived turbulent fluxes is useful for climate studies, forcing of ocean models, and validation of coupled ocean-atmosphere global models and can be accessed through the NASA/GSFC Distributed Active Archive Center.

Description: Two loop heat pipes (LHPs) are to be used for tight thermal control of the Geoscience Laser Altimeter System (GLAS) instrument, planned for flight in late 2001. The LHPs are charged with Propylene as a working fluid. One LHP will be used to transport 110 W from a laser a radiator, the other will transport 190 W from electronic boxes to a separate radiator. The application includes a large amount of thermal mass in each LHP system and low initial startup powers. This along with some non-ideal flight design compromises, such as a less than ideal charge level for this design concept with a symmetrical secondary wick, lead to inadequate performance of the flight LHPs during the flight thermal vacuum test in October of 2000. This presentation focuses on identifying; the sources of the flight test difficulties by modifying the charge and test setup of the successfully tested Development Model Loop Heat Pipe (DM LHP). While very similar to the flight design, the DM L14P did have several significant difference in design and method of testing. These differences were evaluated for affect on performance by conforming the DM LHP to look more like the flight units. The major difference that was evaluated was the relative fill level of the working fluid within the concentrically design LHP compensation chamber. Other differences were also assessed through performance testing including starter heater size and "hot biasing" of major interior components. Performance was assessed with respect to startup, low power operation, conductance, and control heater power. The results of the testing showed that performance improves as initial charge increases, and when the starter heater is made smaller. The "hot biasing" of the major components did not appear to have a detrimental effect. As a result of test results of the DM LHP, modifications are being made to the flight units to increase the fluid charge and increase the watt-density of the starter heater.

Description: An axisymmetric traveling magnetic field induces a meridional base flow in a cylindrical zone of an electrically conducting liquid. This remotely induced flow can be conveniently controlled, in magnitude and direction, and can have benefits for crystal growth applications. In particular, it can be used to offset natural convection. For long vertical cylinders, non-uniform and localized in the propagating direction, magnetic fields are required for this purpose. Here we investigate a particular form of this field, namely that induced by a set of a few electric current coils. An order of magnitude reduction of buoyancy convection is theoretically demonstrated for a vertical Bridgman crystal growth configuration.

Description: Measurement of three-dimensional (3-D) three-component velocity fields is of great importance in both ground and space experiments for understanding materials processing and fluid physics. The experiments in these fields most likely inhibit the application of conventional planar probes for observing 3-D phenomena. Here, we present the investigation results of stereoscopic tracking velocimetry (STV) for measuring 3-D velocity fields, which include diagnostic technology development, experimental velocity measurement, and comparison with analytical and numerical computation. STV is advantageous in system simplicity for building compact hardware and in software efficiency for continual near-real-time monitoring. It has great freedom in illuminating and observing volumetric fields from arbitrary directions. STV is based on stereoscopic observation of particles-Seeded in a flow by CCD sensors. In the approach, part of the individual particle images that provide data points is likely to be lost or cause errors when their images overlap and crisscross each other especially under a high particle density. In order to maximize the valid recovery of data points, neural networks are implemented for these two important processes. For the step of particle overlap decomposition, the back propagation neural network is utilized because of its ability in pattern recognition with pertinent particle image feature parameters. For the step of particle tracking, the Hopfield neural network is employed to find appropriate particle tracks based on global optimization. Our investigation indicates that the neural networks are very efficient and useful for stereoscopically tracking particles. As an initial assessment of the diagnostic technology performance, laminar water jets with and without pulsation are measured. The jet tip velocity profiles are in good agreement with analytical predictions. Finally, for testing in material processing applications, a simple directional solidification apparatus is built for experimenting with a metal analog of succinonitrile. Its 3-D velocity field at the liquid phase is then measured to be compared with those from numerical computation. Our theoretical, numerical, and experimental investigations have proven STV to be a viable candidate for reliably measuring 3-D flow velocities. With current activities are focused on further improving the processing efficiency, overall accuracy, and automation, the eventual efforts of broad experimental applications and concurrent numerical modeling validation will be vital to many areas in fluid flow and materials processing.

Description: A fast multigrid solver for the steady incompressible Navier-Stokes equations is presented. Unlike time-marching schemes, this approach uses relaxation of the steady equations. Application of this method results in a discretization that correctly distinguishes between the advection and elliptic parts of the operator, allowing efficient smoothers to be constructed. Numerical solutions are shown for flow over a flat plate and a Karman-Trefftz airfoil. Using collective Gauss-Seidel line relaxation in both the vertical and horizontal directions, multigrid convergence behavior approaching that of O(N) methods is achieved. The computational efficiency of the numerical scheme is compared with that of a Runge-Kutta based multigrid method.

Description: This paper describes the application of a new Direct Simulation Monte Carlo (DSMC) code, the Molecular Beam Simulator (MBS), which is designed to analyze laboratory scale molecular beam-surface (and crossed-beam) experiments. The MBS is primarily intended to model experiments associated with spacecraft contamination effects, but it can also be used to simulate a variety of surface chemistry and reactive flow measurements. The MBS code is fully three-dimensional, includes a wide-range of chemical processes, and can model one or multiple pulsed (non-steady) sources. As an example application of the MBS code, a fast, pulsed, oxygen atom-surface experiment which examines the chemistry behind erosion of graphite by oxygen atoms is analyzed. Unsteady DSMC simulations show that experimental observations of excited molecular states after the pulse has hit the surface are consistent with two distinct chemical mechanisms: a direct one where the excited molecules are formed on the surface, and a two-step mechanism where ground state molecules formed on the surface are collisionally excited after they leave the surface by trailing oxygen atoms in the pulse. Further DSMC calculations suggest experiments which can distinguish between these mechanisms.

Description: Texture advection is an effective tool for animating and investigating 2D flows. In this paper, we discuss how this technique can be extended to 3D flows. In particular, we examine the use of 3D and 4D textures on 3D synthetic and computational fluid dynamics flow fields.

Description: The advent of advanced computer architectures and parallel computing have led to a revolutionary change in the design process for turbomachinery components. Two- and three-dimensional steady-state computational flow procedures are now routinely used in the early stages of design. Unsteady flow analyses, however, are just beginning to be incorporated into design systems. This paper outlines the transition of a three-dimensional unsteady viscous flow analysis from the research environment into the design environment. The test case used to demonstrate the analysis is the full turbine system (high-pressure turbine, inter-turbine duct and low-pressure turbine) from an advanced turboprop engine.

Description: The effects of multiple mode periodic excitation on the evolution of a circular turbulent jet were studied experimentally. A short, wide-angle diffuser was attached to the jet exit. Streamwise and cross-stream excitations were introduced at the junction between the jet exit and the diffuser inlet on opposing sides of the jet. The introduction of high amplitude, periodic excitation in the streamwise direction enhances the mixing and promotes attachment of the jet shear-layer to the diffuser wall. Cross-stream excitation applied over a fraction of the jet circumference can deflect the jet away from the excitation slot. The two modes of excitation were combined using identical frequencies and varying the relative phase between the two actuators in search of an optimal response. It is shown that, for low and moderate periodic momentum input levels, the jet deflection angles depend strongly on the relative phase between the two actuators. Optimum performance is achieved when the phase difference is pi +/- pi/6. The lower effectiveness of the equal phase excitation is attributed to the generation of an azimuthally symmetric mode that does not produce the required non-axisymmetric vectoring. For high excitation levels, identical phase becomes more effective, while phase sensitivity decreases. An important finding was that with proper phase tuning, two unsteady actuators can be combined to obtain a non-linear response greater than the superposition of the individual effects.

Description: The effects of multiple mode periodic excitation on the evolution of a circular turbulent jet were studied experimentally. A short, wide-angle diffuser was attached to the jet exit. Streamwise and cross-stream excitations were introduced at the junction between the jet exit and the diffuser inlet on opposing sides of the jet. The introduction of high amplitude, periodic excitation in the streamwise direction enhances the mixing and promotes attachment of the jet shear-layer to the diffuser wall. Cross-stream excitation applied over a fraction of the jet circumference can deflect the jet away from the excitation slot. The two modes of excitation were combined using identical frequencies and varying the relative phase between the two actuators in search of an optimal response. It is shown that, for low and moderate periodic momentum input levels, the jet deflection angles depend strongly on the relative phase between the two actuators. Optimum performance is achieved when the phase difference is pi +/- pi/6. The lower effectiveness of the equal phase excitation is attributed to the generation of an azimuthally symmetric mode that does not produce the required non-axisymmetric vectoring. For high excitation levels, identical phase becomes more effective, while phase sensitivity decreases. An important finding was that with proper phase tuning, two unsteady actuators can be combined to obtain a non-linear response greater than the superposition of the individual effects.

Description: A microchannel separator, with 2.7 millimeters as the smallest dimension, was tested, and a pore throat structure captured and removed liquid from a gas-liquid stream. The microchannel device was tested over a of gas and liquid flow rates ranging from 0.0005 up to 0. 14 volume fraction of liquid. Four liquids were tested with air. The biggest factor affecting the throughput is the capacity of liquid flow through the pore throat, which is dictated by permeability, liquid viscosity, flow area, pore throat thickness, and pressure difference across the pore throat. Typically, complete separation of gas and liquid fractions was lost when the liquid flow rate reached about 40 to 60% of the pore throat capacity. However, this could occur over a range of 10 to 90% utilization of pore throat capacity. Breakthrough occurs in the microchannel phase separator at conditions similar to the annular to plug flow transition of two-phase microgravity pipe flow implying that operating in the proper flow regime is crucial. Analysis indicates that the Bond number did not affect performance, supporting the premise that hydrodynamic, interfacial, and capillary forces are more important than gravity. However, the relative importance of gravity is better discerned through testing under reduced gravity conditions.

Description: We have studied thermal expansion of the surface layers of the hexagonally reconstructed Au (001) surface using a classical Molecular Dynamics (MD) simulation technique with an Embedded Atomic Method (EAM) type many-body potential. We find that the top-most hexagonal layer contracts as temperature increases, whereas the second layer expands or contracts depending on the system size. The magnitude of expansion coefficient of the top layer is much larger than that of the other layers. The calculated thermal expansion coefficients of the top-most layer are about -4.93 x 10(exp -5)Angstroms/Kelvin for the (262 x 227)Angstrom cluster and -3.05 x 10(exp -5)Angstroms/Kelvin for (101 x 87)Angstrom cluster. The Fast Fourier Transform (FFT) image of the atomic density shows that there exists a rotated domain of the top-most hexagonal cluster with rotation angle close to 1 degree at temperature T less than 1000Kelvin. As the temperature increases this domain undergoes a surface orientational phase transition. These predictions are in good agreement with previous phenomenological theories and experimental studies.

Description: The efficiency and accuracy of several time integration schemes are investigated for the unsteady Navier-Stokes equations. This study focuses on the efficiency of higher-order Runge-Kutta schemes in comparison with the popular Backward Differencing Formulations. For this comparison an unsteady two-dimensional laminar flow problem is chosen, i.e., flow around a circular cylinder at Re = 1200. It is concluded that for realistic error tolerances (smaller than 10(exp -1)) fourth-and fifth-order Runge-Kutta schemes are the most efficient. For reasons of robustness and computer storage, the fourth-order Runge-Kutta method is recommended. The efficiency of the fourth-order Runge-Kutta scheme exceeds that of second-order Backward Difference Formula by a factor of 2.5 at engineering error tolerance levels (10(exp -1) to 10(exp -2)). Efficiency gains are more dramatic at smaller tolerances.

Description: A new boundary condition is presented.for simulating the flow over passively porous surfaces. The model builds on the prior work of R.H. Bush to eliminate the need for constructing grid within an underlying plenum, thereby simplifying the numerical modeling of passively porous flow control systems and reducing computation cost. Code experts.for two structured-grid.flow solvers, TLNS3D and CFL3D. and one unstructured solver, USM3Dns, collaborated with an experimental porosity expert to develop the model and implement it into their respective codes. Results presented,for the three codes on a slender forebody with circumferential porosity and a wing with leading-edge porosity demonstrate a good agreement with experimental data and a remarkable ability to predict the aggregate aerodynamic effects of surface porosity with a simple boundary condition.

Description: The partial pressures of Te2 in equilibrium with Ga(1-x)Te(x) samples were measured by optical absorption technique from 450 to 1100 C for compositions, x, between 0.333 and 0.612. To establish the relationship between the partial pressure of Te, and the measured optical absorbance, the calibration runs of a pure Te sample were also conducted to determine the Beer's Law constants. The partial pressures of Te2 in equilibrium with the GaTe(s) and Ga2Te3(s)compounds, or the so-called three-phase curves, were established. These partial pressure data imply the existence of the Ga3Te4(s) compound. From the partial pressures of Te2 over the Ga-Te melts, partial molar enthalpy and entropy of mixing for Te were derived and they agree reasonable well with the published data. The activities of Te in the Ga-Te melts were also derived from the measured partial pressures of Te2. These data agree well with most of the previous results. The possible reason for the high activity of Te measured for x less than 0.60 is discussed.

Description: In this report, we present some results of problems investigated during joint research between the Hampton University Fluid Mechanics and Acoustics Laboratory (HU/FM&AL), NASA GRC, and the LaRC Hyper-X Program. This work is supported by joint research between the NASA GRC and the Institute of Mechanics at Moscow State University (IM/MSU) in Russia under a CRDF grant. The main areas of current scientific interest of the HU/FM&AL include an investigation of the proposed and patented advanced methods for aircraft engine thrust and noise benefits. These methods are based on nontraditional 3D corrugated and composite nozzle, inlet, propeller and screw designs such as a Bluebell and Telescope nozzle, Mobius-shaped screw, etc. This is the main subject of our other projects, of which one is presented at the current conference. Here we analyze additional methods for exhaust jet noise reduction without essential thrust loss and even with thrust augmentation. Such additional approaches are: (1) to add some solid, fluid, or gas mass at discrete locations to the main supersonic gas stream to minimize the negative influence of strong shock waves formed in propulsion systems. This mass addition may be accompanied by heat addition to the main stream as a result of the fuel combustion or by cooling of this stream as a result of the liquid mass evaporation and boiling; (2) Use of porous or permeable nozzles and additional shells at the nozzle exit for preliminary cooling of the hot jet exhaust and pressure compensation for off-design conditions (so-called continuous ejector with small mass flow rate); and (3) to propose and analyze new effective methods of fuel injection into the flow stream in air-breathing engines. The research is focused on a wide regime of problems in the propulsion field as well as in experimental testing and theoretical and numerical simulation analyses for advanced aircraft and rocket engines. The FM&AL Team uses analytical methods, numerical simulations, and experimental tests at the Hampton University campus, NASA, and IM/MSU. The main outcomes during this reporting period are: (1) Publications: The AIAA Paper #01-1893 has been accepted for the AIAA/NAL-NASDA-ISAS 10th International Space Planes and Hypersonic Systems and Technologies Conference, 24-27 April 2001, Kyoto, Japan. The AIAA Paper #01-3204 has been accepted for presentation at the 37th AIAA/ASME/SAE/ASEE Joint Propulsion Conference being held on 08-11 July, in Salt Lake City, UT; (2) Grants and proposals: The HU/FM&AL was awarded the NASA grant NAG3-2495 in October 2000 and the CRDF award was granted to the NASA GRC-HU FM&AL and IM/MSU (Russia) in July 2000. A solicited proposal was submitted for the NASA NRA-01GRC-02 competition and two unsolicited proposals to NASA are in preparation; (3) Theory and numerical simulations: Analytical theory, numerical simulation, comparison of theoretical with experimental results, and modification of theoretical approaches, models, grids, etc. Such investigations have been conducted for three main problems: (a) Combustion efficiency optimization in the half-duct combustor system; (b) Drag reduction effects for blunt bodies with solid needles; and (c) Solid particle injection from the butt-end against a supersonic flow. The NASA CFL3D, HU/FM&AL, and IM/MSU GODUNOV codes were used and modified for solution of these problems. The codes are based on full Euler and Navier-Stokes solvers with and without nonequilibrium oxygen-nitrogen and air-hydrogen chemical reactions in laminar and turbulent gas flow regimes; (4) Experimental tests: In the IM/MSU supersonic wind tunnel, experimental tests were conducted with different number, location, and geometric parameters of solid needles mounted at the front of the butt-end in supersonic flow. Optimal parameters were determined that provide minimal butt-end drag in the stationary flow regime. Experimental and numerical simulation results are in good agreement; (5) Student Research Activity: Involvement of one graduate and two undergraduate students as research assistants in the current project.

Description: A systematic investigation of the structure of turbulent jets before their interaction with shock or expansion waves was undertaken during the last year. In particular compressibility and density effects in circular jets issuing in still air were investigated experimentally. Jets with nitrogen, helium, and krypton gases at 0.3, 0.6, and 0.9 Mach numbers were investigated in detail. Particle Image Velocimetry technique was developed, tested, and used to obtain qualitative information of the two-dimensional velocity field on a plane inside the flow field, which was illuminated by a laser sheet. The motion of particles was recorded by a CCD camera, which was appropriately triggered to capture two images within a fraction of a microsecond. Statistical averaging of the data at each location reduced the large amount of acquired data. It was found that the spreading rate of the jets was reduced with increased Mach numbers or increased density ratio. It was also found that decay rates of centerline Mach numbers are higher in gases with reduced density ratio. Mach number fluctuations appear to decrease with increasing Mach number of the flow. It has been proposed that the reason for this behavior is the reduction of vortex stretching activities with increased Mach number.

Description: A prerequisite to study phenomena in the winter stratospheric polar vortex is the separation of measurements inside and outside the dynamical barrier of the vortex edge. We describe a technique to accurately determine the inner edge of the vortex boundary region from measurements of potential temperature and a trace gas, such as N2O, and apply it to in situ aircraft and balloon measurements from the SOLVE/THESE02000 Arctic campaign. The method may be used to refine the Nash algorithm, which, due to the inherently coarser resolution of potential vorticity on which it is dependent, may misidentify the inner edge by up to 466 km, and omit the identification of small, extra-vortex filaments within the vortex.

Description: The flow of metal during Friction Stir Welding is clarified using a faying surface tracer and a nib frozen in place during welding. It is shown that material is transported by two processes. The first is a wiping of material from the advancing front side of the nib onto a plug of material that rotates and advances with the nib. The material undergoes a helical motion within the plug that both rotates and advances with the plug and descends in the wash of the threads on the nib and rises on the outer part of the plug. After one or more rotations, this material is sloughed off the plug in its wake, primarily on the advancing side. The second process is an entrainment of material from the front retreating side of the nib that fills in between the sloughed off pieces from the advancing side. These two processes produce material with different mechanical properties and the strength of a weld should depend on the relative importance of the processes.

Description: This article details an effort to improve the understanding and prediction of turbulent flow inside a droplet of molten metal levitated in an electromagnetic field. It is shown that the flow field in a test case, a nickel droplet levitated under microgravity conditions, is in the transitional regime between laminar and turbulent flow. Past research efforts have used laminar, enhanced viscosity, and k-epsilon turbulence models to describe the flow. The method highlighted in our study is the renormalization group (RNG) algorithm. We show that an accurate description of the turbulent eddy viscosity is critical in order to obtain realistic velocity fields, and that the turbulent eddy viscosity cannot be uniform in levitated droplets. The RNG method does not impose isotropic length or time scales on the flow field, thus allowing such nonuniform features to be captured. A number of other materials processing applications exhibit similarly complex flow characteristics, such as highly recirculating, transitional, and free surface flows, for which this modeling approach may prove useful.

Description: A numerical scheme for direct numerical simulation of shock-turbulence interactions of high speed compressible flows would ideally not be significantly more expensive than the standard fourth or sixth-order compact or non-compact central differencing scheme. It should be possible to resolve all scales down to scales of order of the Kolmogorov scales of turbulence accurately and efficiently, while at the same time being able to capture steep gradients occurring at much smaller scales efficiently. The goal of this lecture is to review the progress and new development of the low dissipative high order shock-capturing schemes proposed by Yee et al. Comparison on the efficiency and accuracy of this class of schemes with spectral and the fifth-order WENO (weighted essentially nonoscillatory) scheme will be presented. A new approach to dynamically sense the appropriate amount of numerical dissipation to be added at each grid point using non-orthogonal wavelets will be discussed.

Description: A comprehensive package of scalable overset grid CFD software is reviewed. The software facilitates accurate simulation of complete aircraft aerodynamics, including viscous effects, unsteadiness, and relative motion between component parts. The software significantly lowers the manpower and computer costs normally associated with such efforts. The software is discussed in terms of current capabilities and planned future enhancements.

Description: We have analyzed Doppler velocity images from the MDI instrument on SOHO to determine the latitudinal transport of angular momentum by the cellular photospheric flows. Doppler velocity images from 60-days in May to July of 1996 were processed to remove the p-mode oscillations, the convective blue shift, the axisymmetric flows, and any instrumental artifacts. The remaining cellular flows were examined for evidence of latitudinal angular momentum transport. Small cells show no evidence of any such transport. Cells the size of supergranules (30,000 km in diameter) show strong evidence for a poleward transport of angular momentum. This would be expected if supergranules are influenced by the Coriolis force, and if the cells are elongated in an east-west direction. We find good evidence for just such an east-west elongation of the supergranules. This elongation may be the result of differential rotation shearing the cellular structures. Data simulations of this effect support the conclusion that elongated supergranules transport angular momentum from the equator toward the poles, Cells somewhat larger than supergranules do not show evidence for this poleward transport. Further analysis of the data is planned to determine if the direction of angular momentum transport reverses for even larger cellular structures. The Sun's rapidly rotating equator must be maintained by such transport somewhere within the convection zone.

Description: Hysteresis and wavenumber vacillation are studied numerically in a weakly stratified quasigeostrophic model. In general, the amplitude of the most unstable wave increases, as the flow becomes more unstable. When the wave becomes saturated, the next longer wave will grow at the expanse of the most unstable wave and becomes the dominant wave. However, once the longwave state is established, it may remain in that regime even as the instability is decreased beyond the threshold where it first developed, thus constituting a hysteresis loop. In a highly unstable case, the flow may not show a preference for any single wave. Instead, the dominant wave aperiodically varies among several long waves. This phenomenon is known as wavenumber vacillation. Hysteresis is further examined in terms of eddy heat flux. It is shown that total eddy heat flux increases as the flow becomes more unstable, but displays a sharp drop when transition to a longer wave occurs. However, in a longwave state, the heat flux always decreases with decreasing instability even pass the threshold when wave transition first occurs.

Description: A comparative experimental study of air density fluctuations in the unheated plumes of a circular, 4-tabbed-circular, chevron-circular and 10-lobed rectangular nozzles was performed at a fixed Mach number of 0.95 using a recently developed Rayleigh scattering based technique. Subsequently, the flow density fluctuations are cross-correlated with the far field sound pressure fluctuations to determine sources for acoustics emission. The nearly identical noise spectra from the baseline circular and the chevron nozzles are found to be in agreement with the similarity in spreading, turbulence fluctuations, and flow-sound correlations measured in the plumes. The lobed nozzle produced the least low frequency noise, in agreement with the weakest overall density fluctuations and flow-sound correlation. The tabbed nozzle took an intermediate position in the hierarchy of noise generation, intensity of turbulent fluctuation and flow-sound correlation. Some of the features in the 4-tabbed nozzle are found to be explainable in terms of splitting of the jet in a central large core and 4 side jetlets.

Description: The thermal-fluid coupling problems are very important to aerospace and engineering applications. Instead of analyzing heat transfer and fluid flow separately, this study merged two well-accepted engineering solution methods, SINDA for thermal analysis and FDNS for fluid flow simulation, into a unified multi-disciplinary thermal fluid prediction method. A fully conservative patched grid interface algorithm for arbitrary two-dimensional and three-dimensional geometry has been developed. The state-of-the-art parallel computing concept was used to couple SINDA and FDNS for the communication of boundary conditions through PVM (Parallel Virtual Machine) libraries. Therefore, the thermal analysis performed by SINDA and the fluid flow calculated by FDNS are fully coupled to obtain steady state or transient solutions. The natural convection between two thick-walled eccentric tubes was calculated and the predicted results match the experiment data perfectly. A 3-D rocket engine model and a real 3-D SSME geometry were used to test the current model, and the reasonable temperature field was obtained.

Description: A general purpose, one dimensional fluid flow code is currently being interfaced with the thermal analysis program SINDA/G. The flow code, GFSSP, is capable of analyzing steady state and transient flow in a complex network. The flow code is capable of modeling several physical phenomena including compressibility effects, phase changes, body forces (such as gravity and centrifugal) and mixture thermodynamics for multiple species. The addition of GFSSP to SINDA/G provides a significant improvement in convective heat transfer modeling for SINDA/G. The interface development is conducted in multiple phases. This paper describes the first phase of the interface which allows for steady and quasisteady (unsteady solid, steady fluid) conjugate heat transfer modeling.

Description: The purpose of this review is to detail the current theoretical understanding of the origin of piezoelectric and ferroelectric phenomena in polymers; to present the state-of-the-art in piezoelectric polymers and emerging material systems that exhibit promising properties; and to discuss key characterization methods, fundamental modeling approaches, and applications of piezoelectric polymers. Piezoelectric polymers have been known to exist for more than forty years, but in recent years they have gained notoriety as a valuable class of smart materials.

Description: In this paper, we report a version of the Space-Time Conservation Element and Solution Element (CE/SE) Method in which the 2D and 3D unsteady Euler equations are simulated using structured or unstructured quadrilateral and hexahedral meshes, respectively. In the present method, mesh values of flow variables and their spatial derivatives are treated as independent unknowns to be solved for. At each mesh point, the value of a flow variable is obtained by imposing a flux conservation condition. On the other hand, the spatial derivatives are evaluated using a finite-difference/weighted-average procedure. Note that the present extension retains many key advantages of the original CE/SE method which uses triangular and tetrahedral meshes, respectively, for its 2D and 3D applications. These advantages include efficient parallel computing ease of implementing non-reflecting boundary conditions, high-fidelity resolution of shocks and waves, and a genuinely multidimensional formulation without using a dimensional-splitting approach. In particular, because Riemann solvers, the cornerstones of the Godunov-type upwind schemes, are not needed to capture shocks, the computational logic of the present method is considerably simpler. To demonstrate the capability of the present method, numerical results are presented for several benchmark problems including oblique shock reflection, supersonic flow over a wedge, and a 3D detonation flow.

Description: Considerable progress over the past thirty years has been made in the development of large-scale computational fluid dynamics (CFD) solvers for the Euler and Navier-Stokes equations. Computations are used routinely to design the cruise shapes of transport aircraft through complex-geometry simulations involving the solution of 25-100 million equations; in this arena the number of wind-tunnel tests for a new design has been substantially reduced. However, simulations of the entire flight envelope of the vehicle, including maximum lift, buffet onset, flutter, and control effectiveness have not been as successful in eliminating the reliance on wind-tunnel testing. These simulations involve unsteady flows with more separation and stronger shock waves than at cruise. The main reasons limiting further inroads of CFD into the design process are: (1) the reliability of turbulence models; and (2) the time and expense of the numerical simulation. Because of the prohibitive resolution requirements of direct simulations at high Reynolds numbers, transition and turbulence modeling is expected to remain an issue for the near term. The focus of this paper addresses the latter problem by attempting to attain optimal efficiencies in solving the governing equations. Typically current CFD codes based on the use of multigrid acceleration techniques and multistage Runge-Kutta time-stepping schemes are able to converge lift and drag values for cruise configurations within approximately 1000 residual evaluations. An optimally convergent method is defined as having textbook multigrid efficiency (TME), meaning the solutions to the governing system of equations are attained in a computational work which is a small (less than 10) multiple of the operation count in the discretized system of equations (residual equations). In this paper, a distributed relaxation approach to achieving TME for Reynolds-averaged Navier-Stokes (RNAS) equations are discussed along with the foundations that form the basis of this approach. Because the governing equations are a set of coupled nonlinear conservation equations with discontinuities (shocks, slip lines, etc.) and singularities (flow- or grid-induced), the difficulties are many. This paper summarizes recent progress towards the attainment of TME in basic CFD simulations.

Description: For the last several years, Glenn-HT, a three-dimensional (3D) Computational Fluid Dynamics (CFD) computer code for the analysis of gas turbine flow and convective heat transfer has been evolving at the NASA Glenn Research Center. The code is unique in the ability to give a highly detailed representation of the flow field very close to solid surfaces in order to get accurate representation of fluid heat transfer and viscous shear stresses. The code has been validated and used extensively for both internal cooling passage flow and for hot gas path flows, including detailed film cooling calculations and complex tip clearance gap flow 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, but it can be useful for detailed design analysis. In this paper, the code is described and examples of its validation and use for complex flow calculations are presented, emphasizing the applicability to turbomachinery for space launch vehicle propulsion systems.

Description: The paper presents an order-of-magnitude analysis of the physical processes occurring during the pulsations of a vapor bubble subject to a sound field and shows several numerical examples relating to vapor bubbles in water with and without a translational velocity relative to the liquid. Finally, the growth and collapse of a bubble in a small tube under the action of a heat pulse is considered and it is pointed out that, in suitable conditions, a potentially useful pumping effect without mechanical moving parts can be achieved.

Description: A 40-percent scale model of the inlet to a rocket-based combined-cycle (RBCC) engine was tested in the NASA Glenn Research Center 1- by 1-Foot Supersonic Wind Tunnel (SWT). The full-scale RBCC engine is scheduled for test in the Hypersonic Tunnel Facility (HTF) at NASA Glenn's Plum Brook Station at Mach 5 and 6. This engine will incorporate the configuration of this inlet model which achieved the best performance during the present experiment. The inlet test was conducted at Mach numbers of 4.0, 5.0, 5.5, and 6.0. The fixed-geometry inlet consists of an 8 deg.. forebody compression plate, boundary layer diverter, and two compressive struts located within 2 parallel sidewalls. These struts extend through the inlet, dividing the flowpath into three channels. Test parameters investigated included strut geometry, boundary layer ingestion, and Reynolds number (Re). Inlet axial pressure distributions and cross-sectional Pitot-pressure surveys at the base of the struts were measured at varying back-pressures. Inlet performance and starting data are presented. The inlet chosen for the RBCC engine self-started at all Mach numbers from 4 to 6. Pitot-pressure contours showed large flow nonuniformity on the body-side of the inlet. The inlet provided adequate pressure recovery and flow quality for the RBCC cycle even with the flow separation.

Description: An investigation of flow induced vibration due to coupling between the fluid flow and the propellants lines (LOX and RP-1) was performed. Various flow rate conditions were studied to check whether flow induced vibration was possible due to vortex shedding in both valves and pipe lines. Resonance test was conducted for all segments of the LOX-feedline for the preburner under test. In addition, critical values of frequency and velocity are calculated using a mass damping model. A simple chart characterizing the relation between frequency and velocity is developed for each component; i.e. propellant lines, valves and flow meters. It was found that flow induced vibration occurs for various segments with flow rates of 113 lb/s, 275 lb/s and 40 lb/s. Even more interesting using critical conditions for buckling, it was found that the valve or pipe may collapse for a flow rate of 275 lb/s and valve height of 10% of pipe diameter. Furthermore, two models for the acoustic pressure acting on the segments particularly for the valve are proposed.

Description: The Tenth Thermal arid Fluids Analysis Workshop (TFAWS 99) was held at the Bevill Center, University of Alabama in Huntsville, Huntsville, Alabama, September 13-17, 1999. The theme for the hands-on training workshop and conference was "Tools and Techniques Contributing to Engineering Excellence". Forty-seven technical papers were presented in four sessions. The sessions were: (1) Thermal Spacecraft/Payloads, (2) Thermal Propulsion/Vehicles, (3) Interdisciplinary Paper, and (4) Fluids Paper. Forty papers were published in these proceedings. The remaining seven papers were not available in electronic format at the time of publication. In addition to the technical papers, there were (a) nine hands-on classes on thermal and flow analyses software, (b) twelve short courses, (c) thirteen product overview lectures, and (d) three keynote lectures. The workshop resulted in participation of 171 persons representing NASA Centers, Government agencies, aerospace industries, academia, software providers, and private corporations.

Description: Computational fluid dynamics (CFD) is used to simulate the wake flow and afterbody heating of the Project Fire II ballistic reentry to Earth at 11.4 km/sec. Laminar results are obtained over a portion of the trajectory between the initial heat pulse and peak afterbody heating. Although non-catalytic forebody convective heating results are in excellent agreement with previous computations, initial predictions of afterbody heating were about a factor of two below the experimental values. Further analysis suggests that significant catalysis may be occurring on the afterbody heat shield. Computations including finite-rate catalysis on the afterbody surface are in good agreement with the data over the early portion of the trajectory, but are conservative near the peak afterbody heating point, especially on the rear portion of the conical frustum. Further analysis of the flight data from Fire II shows that peak afterbody heating occurs before peak forebody heating, a result that contradicts computations and flight data from other entry vehicles. This result suggests that another mechanism, possibly pyrolysis, may be occurring during the later portion of the trajectory, resulting in less total heat transfer than the current predictions.

Description: Just after formation, optical fibers are wetted stably with acrylate at capillary numbers routinely exceeding 1000. It is hypothesized that this is possible because of dissolution of air into the liquid coating. A lubrication/boundary integral analysis that includes gas diffusion and solubility is developed. It is applied using conservatively estimated solubility and diffusivity coefficients and solutions are found that are consistent with industry practice and with the hypothesis. The results also agree with the claim of Deneka, Kar & Mensah (1988) that the use of high solubility gases to bathe a wetting line allows significantly greater wetting speeds. The solutions indicate a maximum speed of wetting which increases with gas solubility and with reduction in wetting-channel diameter.

Description: Experiments were conducted in Earth gravity and microgravity to acquire quantitative data on near field flow structure of helium jets injected into air. Microgravity conditions were simulated in the 2.2-second drop tower at NASA Glenn Research Center. The jet flow was observed by quantitative rainbow schlieren deflectometry, a non-intrusive line of site measurement technique for the whole field. The flow structure was characterized by distributions of angular deflection and helium mole percentage obtained from color schlieren images taken at 60 Hz. Results show that the jet flow was significantly influenced by the gravity. The jet in microgravity was up to 70 percent wider than that in Earth gravity. The jet flow oscillations observed in Earth gravity were absent in microgravity, providing direct experimental evidence that the flow instability in the low density jet was buoyancy induced. The paper provides quantitative details of temporal flow evolution as the experiment undergoes a change in gravity in the drop tower.

Description: The long-term goal of this project is to study the nucleation of quantized vortices in helium II by investigating the behavior of rotating droplets of helium II in a reduced gravity environment. The objective of this ground-based research grant was to develop new experimental techniques to aid in accomplishing that goal. The development of an electrostatic levitator for superfluid helium, described below, and the successful suspension of charged superfluid drops in modest electric fields was the primary focus of this work. Other key technologies of general low temperature use were developed and are also discussed.

Description: Reactive sputtering was used to grow thin films of Sc(x)Ga(1-x)N with scandium concentrations of 20%-70% on quartz substrates at temperatures of 300-675 K. X-ray diffraction (XRD) of the films showed either weak or no structure, suggesting the films are amorphous or microcrystalline. Optical absorption spectra were taken of each sample and the optical band gap was determined. The band gap varied linearly with increasing Ga concentration between 2.0 and 3.5 eV. Ellipsometry was used to confirm the band gap measurements and provide optical constants in the range 250-1200 nm. ScN and GaN have different crystal structures (rocksalt and wurzite, respectively), and thus may form a heterogeneous mixture as opposed to an alloy. Since the XRD data were inconclusive, bilayers of ScN/GaN were grown and optical absorption spectra taken. A fundamental difference in the spectra between the bilayer films and alloy films was seen, suggesting the films are alloys, not heterogeneous mixtures.

Description: This paper presents the results of a numerical study that examines the surface heating discrepancies observed between computed and measured values along a sharp cone. With Mach numbers of an order of 10 and the freestream length Reynolds number of an order of 10 000, the present computations have been made with the direct simulation Monte Carlo (DSMC) method by using the G2 code of Bird. The flow conditions are those specified for two experiments conducted in the Veridian 48-inch Hypersonic Shock Tunnel. Axisymmetric simulations are made since the test model was assumed to be at zero incidence. Details of the current calculations are presented, along with comparisons between the experimental data, for surface heating and pressure distributions. Results of the comparisons show major differences in measured and calculated results for heating distributions, with differences in excess of 25 percent for the two cases examined.

Description: Research in the FY01 focused on the analysis and development of enhanced algorithms for unsteady aerodynamics and chemically reacting flowfields. The research was performed in support of NASA Ames' efforts to improve the capabilities of the in-house computational fluid dynamics code, OVERFLOW. Specifically, the research was focused on the four areas: (1) investigation of stagnation region effects; (2) unsteady preconditioning dual-time procedures; (3) dissipation formulation for combustion; and (4) time-stepping methods for combustion.

Description: For several decades the NASA Glenn Research Center has been providing a file of thermodynamic data for use in several computer programs. These data are in the form of least-squares coefficients that have been calculated from tabular thermodynamic data by means of the NASA Properties and Coefficients (PAC) program. The source thermodynamic data are obtained from the literature or from standard compilations. Most gas-phase thermodynamic functions are calculated by the authors from molecular constant data using ideal gas partition functions. The Coefficients and Properties (CAP) program described in this report permits the generation of tabulated thermodynamic functions from the NASA least-squares coefficients. CAP provides considerable flexibility in the output format, the number of temperatures to be tabulated, and the energy units of the calculated properties. This report provides a detailed description of input preparation, examples of input and output for several species, and a listing of all species in the current NASA Glenn thermodynamic data file.

Description: This report summarizes two years of work on investigating algorithms for automatically detecting fringe patterns in images acquired using oil-drop interferometry for the determination of skin friction. Several different analysis methods were tested, and a combination of a windowed Fourier transform followed by a correlation was found to be most effective. The implementation of this method is discussed and details of the process are described. The results indicate that this method shows promise for automating the fringe detection process, but further testing is required.

Description: NASA Glenn Research Center has recently applied the WIND National Code flow solver to an accelerating elbow with a 90 deg. bend to reveal aspects of secondary flow. This elbow was designed by NACA in the early 1950's such that flow separation would be avoided. Experimental testing was also done at that time. The current three dimensional CFD investigation shows that separation has indeed been avoided. Using its three-dimensional capability, this investigation provides various viewpoints in several planes that display the inception, development, and final location of a passage vortex. Its shape first becomes discernible as a vortex near the exit of the bend. This rendition of the exit passage vortex compares well with that found in the experiments. The viewpoints show that the passage vortex settles on the suction surface at the exit about one-third of the distance between the plane wall and midspan. Furthermore, it projects into the mainstream to about one-third of the channel width. Of several turbulence models used in this investigation, the Spalart Alimaras, Baldwin Lomax, and SST (Shear Stress Transport) models were by far the most successful in matching the experiments.

Description: A statistical ensemble of large eddy simulations (LES) is run simultaneously for the same flow. The information provided by the different large scale velocity fields is used to propose an ensemble averaged version of the dynamic model. This produces local model parameters that only depend on the statistical properties of the flow. An important property of the ensemble averaged dynamic procedure is that it does not require any spatial averaging and can thus be used in fully inhomogeneous flows. Also, the ensemble of LES's provides statistics of the large scale velocity that can be used for building new models for the subgrid-scale stress tensor. The ensemble averaged dynamic procedure has been implemented with various models for three flows: decaying isotropic turbulence, forced isotropic turbulence, and the time developing plane wake. It is found that the results are almost independent of the number of LES's in the statistical ensemble provided that the ensemble contains at least 16 realizations.

Description: We report on work carried out to develop a 3-D coupled Finite Volume/BEM-based temperature forward/flux back (TFFB) coupling algorithm to solve the conjugate heat transfer (CHT) which arises naturally in analysis of systems exposed to a convective environment. Here, heat conduction within a structure is coupled to heat transfer to the external fluid which is convecting heat into or out of the solid structure. There are two basic approaches to solving coupled fluid structural systems. The first is a direct coupling where the solution of the different fields is solved simultaneously in one large set of equations. The second approach is a loose coupling strategy where each set of field equations is solved to provide boundary conditions for the other. The equations are solved in turn until an iterated convergence criterion is met at the fluid-solid interface. The loose coupling strategy is particularly attractive when coupling auxiliary field equations to computational fluid dynamics codes. We adopt the latter method in which the BEM is used to solve heat conduction inside a structure which is exposed to a convective field which in turn is resolved by solving the NASA Glenn compressible Navier-Stokes finite volume code Glenn-HT. The BEM code features constant and bi-linear discontinuous elements and an ILU-preconditioned GMRES iterative solver for the resulting non-symmetric algebraic set arising in the conduction solution. Interface of flux and temperature is enforced at the solid/fluid interface, and a radial-basis function scheme is used to interpolated information between the CFD and BEM surface grids. Additionally, relaxation is implemented in passing the fluxes from the conduction solution to the fluid solution. Results from a simple test example are reported.

Description: This document presents the full-scale analyses of the CFD RSRM. The RSRM model was developed with a 20 second burn time. The following are presented as part of the full-scale analyses: (1) RSRM embedded inclusion analysis; (2) RSRM igniter nozzle design analysis; (3) Nozzle Joint 4 erosion anomaly; (4) RSRM full motor port slag accumulation analysis; (5) RSRM motor analysis of two-phase flow in the aft segment/submerged nozzle region; (6) Completion of 3-D Analysis of the hot air nozzle manifold; (7) Bates Motor distributed combustion test case; and (8) Three Dimensional Polysulfide Bump Analysis.

Description: This report outlines the major results of the grant "Fundamental Processes of Atomization in Fluid-Fluid Flows." These include: 1) the demonstration that atomization in liquid/liquid shear flow is driven by a viscous shear instability that triggers the formation of a long thin sheet; 2) discovery of a new mode of interfacial instability for oscillatory two-layer systems whereby a mode that originates within the less viscous liquid phase causes interfacial deformation as the oscillation proceeds; 3) the demonstration that rivulet formation from gravity front occurs because the local front shape specified by gravity and surface tension changes from a nose to a wedge geometry, thus triggering a large increase in viscous resistance; and 4) extension of the studies on nonlinear wave evolution on falling films and in stratified flow, particularly the evolution towards large-amplitude solitary waves that tend to generate drops.

Description: This document was originally published in 1997 as a International Space University Master of Space Studies student's project. While the specifics may be a little dated, the results gleaned are still relevant and valid. Johnson Space Center is now abuzz with preparations for the International Space Station. The work on the various systems for the U.S. Habitation Module will begin in 1998. As a prelude, the Flight Crew Support Division perceived the need to take a closer, more critical took at planning the Whole Body Cleansing function for ISS. This report is an endeavour to retrieve all data available on whole body cleansing mechanisms used, past and present, by the Russians and the Americans, analyze it and create design concepts for products/product systems for zero-g whole body cleansing on ISS Alpha, for typical duration of about 90 days. This report takes a close look at the Skylab collapsible shower, the Mir shower/sauna, the full body cleansing methods currently in use onboard the Space Shuttle and Mir and at the Whole Body Shower designed and tested for Space Station Freedom. It attempts to "listen" carefully to what the Mir astronauts (Norm Thagard, Shannon Lucid and John Blaha) have to say about their personal hygiene experiences during their recent stays on Mir. The findings in the report call for a change in paradigm. What is good for Earth conditions is not necessarily good for Zero-g! It concludes that a shower is not a good idea for the ISS. The final concept that is proposed reflects very strongly what the Mir astronauts would like to have and to use onboard a station like the ISS, The report concludes with directions of how to take the "idea" further and realize it in the form of a product system for Whole Body Cleansing onboard the ISS.

Description: THUNDER (THin UNimorph DrivER) actuators are pre-stressed piezoelectric devices developed at NASA Langley Research Center (LaRC) that exhibit enhanced strain capabilities. As a result, they are of interest in a variety of aerospace applications. Characterization of their performance as a function of electric field, temperature and frequency is needed in order to optimize their operation. Towards that end, a number of THUNDER devices were obtained from FACE International Co. with a stainless steel substrate varying in thickness from 1 mil to 20 mils. The various devices were evaluated to determine low-field and high-field displacement its well as the polarization hysteresis loops. The thermal stability of these drivers was evaluated by two different methods. First, the samples were thermally cycled under electric field by systematically increasing the maximum temperature from 25 C to 200 C while the displacement was being measured. Second, the samples were isothermally aged at 0 C, 50 C, 100 C. and 150 C in air, and the isothermal decay of the displacement was measured at room temperature as a function of time.

Description: Experimental and computational studies have been performed to investigate the detailed distribution of convective heat transfer coefficients on the first-stage blade tip surface for a geometry typical of large power generation turbines (greater than 100 MW) In a previous work the numerical heat transfer results for a sharp edge blade tip and a radiused blade tip were presented. More recently several other tip treatments have been considered for which the tip heat transfer has been measured and documented. This paper is concerned with the numerical prediction of the tip surface heat transfer for radiused blade tip equipped with mean-camberline strip (or "squealer" as it is often called). The heat transfer results are compared with the experimental results and discussed. The effectiveness of the mean-camberline strip in reducing the tip leakage and the tip heat transfer as compared to a radiused edge tip and sharp edge tip was studied. The calculations show that the sharp edge tip works best (among the cases considered) in reducing the tip leakage flow and the tip heat transfer.

Description: Viewgraphs detailing time-dependent simulations of Turbopump flows are presented. The topics include: 1) Major Drivers of Current Work; 2) Solution Methods; and 3) Unsteady Turbopump Flow.

Description: The objectives of this viewgraph presentation are to predict the following: (1) dual-engine base-heating at 57% PL at sea level, and (2) dual-engine base-heating during PPO at three ascent abort trajectories. A systematically anchored computational fluid dynamics and heat transfer three-dimensional transfer simulation is being used to study the effect of reduced power levels on base-heating environments during sea level testing and during PPO. Preliminary results show the following: (1) convective heating is higher for the 57% PL than for 100% PL on most of the pillows and flex seals during sea level testing; and (2) convective heating on pillows and flex seals on the 'off' engine side is higher than that on the 'on' engine side.

Description: This paper presents viewgraphs of time-accurate simulations of unsteady incompressible flow with moving boundaries.

Description: Temporally and spatially-resolved, two-component measurements of velocity in a supersonic hydrogen-air combustor are reported. The combustor had a single unswept ramp fuel injector and operated with an inlet Mach number of 2 and a flow total temperature approaching 1200 K. The experiment simulated the mixing and combustion processes of a dual-mode scramjet operating at a flight Mach number near 5. The velocity measurements were obtained by seeding the fuel with alumina particles and performing Particle Image Velocimetry on the mixing and combustion wake of the ramp injector. To assess the effects of combustion on the fuel air-mixing process, the distribution of time-averaged velocity and relative turbulence intensity was determined for the cases of fuel-air mixing and fuel-air reacting. Relative to the mixing case, the near field core velocity of the reacting fuel jet had a slower streamwise decay. In the far field, downstream of 4 to 6 ramp heights from the ramp base, the heat release of combustion resulted in decreased flow velocity and increased turbulence levels. The reacting measurements were also compared with a computational fluid dynamics solution of the flow field. Numerically predicted velocity magnitudes were higher than that measured and the jet penetration was lower.

Description: An object-oriented event-driven immersive virtual environment (VE) is described for the visualization of computational fluid dynamics (CFD) results. The VE incorporates the following types of primitive software objects: interface objects, support objects, geometric entities, and finite elements. The fluid domain is discretized using either a multi-block structured grid or an unstructured finite element mesh. The VE allows natural 'fly-through' visualization of the model, the CFD grid, and the model's surroundings. In order to help visualize the flow and its effects on the model, the VE incorporates the following objects: stream objects (lines, surface-restricted lines. ribbons. and volumes); colored surfaces; elevation surfaces; surface arrows; global and local iso-surfaces; vortex cores; and separation/attachment surfaces and lines. Most of these objects can be used for dynamically probing the flow. Particles and arrow animations can be displayed on top of stream objects. Primitive response quantities as well as derived quantities can be used. A recursive tree search algorithm is used for real-time point and value search in the CFD grid.

Description: The objective of the current effort is to provide a computational framework for design and analysis of the entire fuel supply system of a liquid rocket engine, including high-fidelity unsteady turbopump flow analysis. This capability is needed to support the design of pump sub-systems for advanced space transportation vehicles that are likely to involve liquid propulsion systems. To date, computational tools for design/analysis of turbopump flows are based on relatively lower fidelity methods. An unsteady, three-dimensional viscous flow analysis tool involving stationary and rotational components for the entire turbopump assembly has not been available for real-world engineering applications. The present effort provides developers with information such as transient flow phenomena at start up, and non-uniform inflows, and will eventually impact on system vibration and structures. In the proposed paper, the progress toward the capability of complete simulation of the turbo-pump for a liquid rocket engine is reported. The Space Shuttle Main Engine (SSME) turbo-pump is used as a test case for evaluation of the hybrid MPI/Open-MP and MLP versions of the INS3D code. CAD to solution auto-scripting capability is being developed for turbopump applications. The relative motion of the grid systems for the rotor-stator interaction was obtained using overset grid techniques. Unsteady computations for the SSME turbo-pump, which contains 114 zones with 34.5 million grid points, are carried out on Origin 3000 systems at NASA Ames Research Center. Results from these time-accurate simulations with moving boundary capability will be presented along with the performance of parallel versions of the code.

Description: This viewgraph presentation provides information on mathematical modelling of the SSME (space shuttle main engine). The unsteady SSME-rig1 start-up procedure from the pump at rest has been initiated by using 34.3 million grid points. The computational model for the SSME-rig1 has been completed. Moving boundary capability is obtained by using DCF module in OVERFLOW-D. MPI (Message Passing Interface)/OpenMP hybrid parallel code has been benchmarked.

Description: The objective of this viewgraph presentation is to enhance incompressible flow simulation capability for developing aerospace vehicle components, especially unsteady flow phenomena associated with high speed turbo pumps. Unsteady Space Shuttle Main Engine (SSME)-rig1 1 1/2 rotations are completed for the 34.3 million grid points model. The moving boundary capability is obtained by using the DCF module. MLP shared memory parallelism has been implemented and benchmarked in INS3D. The scripting capability from CAD geometry to solution is developed. Data compression is applied to reduce data size in post processing and fluid/structure coupling is initiated.

Description: This paper presents a report of a multi year study of the U.S. Army LONGBOW APACHE (AH-64D) aircraft. The goals of this study were to provide the Apache Project Managers Office (PMO) with a broad spectrum of calibrated comprehensive and CFD models of the AH-64D aircraft. The goal of this paper is to present an overview of the comprehensive model which has been developed. The CAMRAD II computer code was chosen to complete this task. The paper first discusses issues that must be addressed when modeling the Apache using CAMRAD. The work required the acquisition of a data base for the aircraft and the development and application of a multidisciplinary computer model. Sample results from various parts of the model are presented. Conclusions with regard to the strengths and weaknesses of simulations based on this model are discussed.

Description: The turbine gas path is a very complex flow field. This is due to a variety of flow and heat transfer phenomena encountered in turbine passages. This manuscript provides an overview of the current work in this field at the NASA Glenn Research Center. Also, based on the author's preference, more emphasis is on the computational work. There is much more experimental work in progress at GRC than that reported here. While much has been achieved, more needs to be done in terms of validating the predictions against experimental data. More experimental data, especially on film cooled and rough turbine blades, are required for code validation. Also, the combined film cooling and internal cooling flow computation for a real blade is yet to be performed. While most computational work to date has assumed steady state conditions, the flow is clearly unsteady due to the presence of wakes. All this points to a long road ahead. However, we are well on course.

Description: This paper presents the results of a numerical study of hypersonic interacting flows at flow conditions that include those for which experiments have been conducted in the Calspan-University of Buffalo Research Center (CUBRC) Large Energy National Shock (LENS) tunnel and the ONERA R5Ch low-density wind tunnel. The computations are made with the direct simulation Monte Carlo (DSMC) method of Bird. The focus is on Mach 9.3 to 11.4 flows about flared axisymmetric configurations, both hollow cylinder flares and double cones. The results presented highlight the sensitivity of the calculations to grid resolution, provide results concerning the conditions for incipient separation, and provide information concerning the flow structure and surface results for the extent of separation, heating, pressure, and skin friction.

Description: The flow over the zero-pressure-gradient So-Mellor convex curved wall is simulated using the Navier-Stokes equations. An inviscid effective outer wall shape, undocumented in the experiment, is obtained by using an adjoint optimization method with the desired pressure distribution on the inner wall as the cost function. Using this wall shape with a Navier-Stokes method, the abilities of various turbulence models to simulate the effects of curvature without the complicating factor of streamwise pressure gradient can be evaluated. The one-equation Spalart-Allmaras turbulence model overpredicts eddy viscosity, and its boundary layer profiles are too full. A curvature-corrected version of this model improves results, which are sensitive to the choice of a particular constant. An explicit algebraic stress model does a reasonable job predicting this flow field. However, results can be slightly improved by modifying the assumption on anisotropy equilibrium in the model's derivation. The resulting curvature-corrected explicit algebraic stress model possesses no heuristic functions or additional constants. It lowers slightly the computed skin friction coefficient and the turbulent stress levels for this case (in better agreement with experiment), but the effect on computed velocity profiles is very small.

Description: This report describes the research work undertaken at the Massachusetts Institute of Technology, under NASA Research Grant NAG4-157. The aim of this research is to identify effective algorithms and methodologies for the efficient and routine solution of flow simulations about complete vehicle configurations. For over ten years we have received support from NASA to develop unstructured mesh methods for Computational Fluid Dynamics. As a result of this effort a methodology based on the use of unstructured adapted meshes of tetrahedra and finite volume flow solvers has been developed. A number of gridding algorithms, flow solvers, and adaptive strategies have been proposed. The most successful algorithms developed from the basis of the unstructured mesh system FELISA. The FELISA system has been extensively for the analysis of transonic and hypersonic flows about complete vehicle configurations. The system is highly automatic and allows for the routine aerodynamic analysis of complex configurations starting from CAD data. The code has been parallelized and utilizes efficient solution algorithms. For hypersonic flows, a version of the code which incorporates real gas effects, has been produced. The FELISA system is also a component of the STARS aeroservoelastic system developed at NASA Dryden. One of the latest developments before the start of this grant was to extend the system to include viscous effects. This required the development of viscous generators, capable of generating the anisotropic grids required to represent boundary layers, and viscous flow solvers. We show some sample hypersonic viscous computations using the developed viscous generators and solvers. Although this initial results were encouraging it became apparent that in order to develop a fully functional capability for viscous flows, several advances in solution accuracy, robustness and efficiency were required. In this grant we set out to investigate some novel methodologies that could lead to the required improvements. In particular we focused on two fronts: (1) finite element methods and (2) iterative algebraic multigrid solution techniques.

Description: An investigation of the control of turbulent boundary layer flow over flexible and rigid surfaces downstream of a concave-convex geometry has been made. The concave-convex curvature induces centrifugal forces and a pressure gradient on the growth of the turbulent boundary layer. The favorable gradient is not sufficient to overcome the unfavorable; thus, the net effect is a destabilizing, of the flow into Gortler instabilities. This study shows that control of the turbulent boundary layer and structural loading can be successfully achieved by using localized surface heating because the subsequent cooling and geometrical shaping downstream over a favorable pressure gradient is effective in laminarization of the turbulence. Wires embedded in a thermally insulated substrate provide surface heating. The laminarized velocity profile adjusts to a lower Reynolds number, and the structure responds to a lower loading. In the laminarization, the turbulent energy is dissipated by molecular transport by both viscous and conductivity mechanisms. Laminarization reduces spanwise vorticity because of the longitudinal cooling gradient of the sublayer profile. The results demonstrate that the curvature-induced mean pressure gradient enhances the receptivity of the flow to localized surface heating, a potentially viable mechanism to laminarize turbulent boundary layer flow; thus, the flow reduces the response of the flexible structure and the resultant sound radiation.