ALBERT

Description: When a pool of flammable liquid is ignited, the flame spread rate can vary widely depending on the initial fuel temperature, pool geometry, ambient atmospheric conditions, and gravitational level. There is substantial decades-old debate in the scientific literature about the role of gravity in these phenomena. The objective of the research was to measure ignition and flame spread across liquid pools in both normal and microgravity with and without forced airflow, and to obtain detailed thermal and velocity field data for comparison to a predictive numerical model that is concurrently being developed. To that end, an experiment known as Spread Across Liquids 5 (SAL 5) was designed to be conducted in a low-gravity environment on a sounding rocket. Unfortunately, during the sounding rocket flight the fuel tray for the experiment failed to fill properly prior to ignition. The primary cause of the failure to fill properly is the geometry of the fuel tray; that is, the presence of the gaps on the side walls and the sharp edge around the thermocouple through-holes. The contamination resulting from the cleaning process is a strong secondary factor which would have prevented proper, timely filling had the primary cause not been present.

Description: The application of a shell/3D modeling technique for the simulation of skin/stringer debond in a specimen subjected to tension and three-point bending was studied. The global structure was modeled with shell elements. A local three-dimensional model, extending to about three specimen thicknesses on either side of the delamination front was used to model the details of the damaged section. Computed total strain energy release rates and mixed-mode ratios obtained from shell/3D simulations were in good agreement with results obtained from full solid models. The good correlation of the results demonstrated the effectiveness of the shell/3D modeling technique for the investigation of skin/stiffener separation due to delamination in the adherents. In addition, the application of the submodeling technique for the simulation of skin/stringer debond was also studied. Global models made of shell elements and solid elements were studied. Solid elements were used for local submodels, which extended between three and six specimen thicknesses on either side of the delamination front to model the details of the damaged section. Computed total strain energy release rates and mixed-mode ratios obtained from the simulations using the submodeling technique were not in agreement with results obtained from full solid models.

Description: This study focuses on localized ignition by external radiant flux and subsequent flame growth over thin polymeric materials (plastic and paper) in microgravity. Two transition stages were observed. The first transition stage covers the period from the onset of ignition to the formation of stabilized flame near the ignited area. This is followed by the second transition of the flame growth stage from the initial stabilized flame to sustained fire growth away from the ignited area. For the first stage, ignition experiments of thin PMMA sheets were conducted using a CO2 laser as an external source in the 10 s drop tower. The results of front side surface ignition and of backside surface ignition were observed. The effects of imposed flow velocity, sample thickness, and ambient oxygen concentration on ignition are obtained. Numerical study was conducted to investigate to understand and predict ignition behavior observed in the experiments. For the second stage, numerical study is being conducted to describe the effects of gravity on heat release rate of a PMMA sheet. The gravity level was varied from zero to normal gravity. The preliminary results show that the maximum heat release occurs at around 0.02 g.

Description: We describe a ground-based apparatus that allows the cancellation of gravity on a fluid using magnetic forces. The present system was designed for liquid oxygen studies over the range 0.001 - 5 g s. This fluid is an essential component of any flight mission using substantial amounts of liquid propellant, especially manned missions. The apparatus has been used to reduce the hydrostatic compression near the oxygen critical point and to demonstrate inverted phase separation. It could also be used to study pool boiling and two-phase heat transfer in Martian, Lunar or near-zero gravity, as well as phenomena such as Marangoni flow and convective instabilities. These studies would contribute directly to the reliability and optimization of the Moon and Mars flight programs.

Description: Spray cooling has high potential in thermal management and life support systems by overcoming the deleterious effect of microgravity upon two-phase heat transfer. In particular spray cooling offers several advantages in heat flux removal that include the following: 1) By maintaining a wetted surface, spray droplets impinge upon a thin fluid film rather than a dry solid surface; 2. Most heat transfer surfaces will not be smooth but rough. Roughness can enhance conductive cooling, aid liquid removal by flow channeling; and 3. Spray momentum can be used to a) substitute for gravity delivering fluid to the surface, b) prevent local dryout and potential thermal runaway and c) facilitate liquid and vapor removal. Yet high momentum results in high We and Re numbers characterizing the individual spray droplets. Beyond an impingement threshold, droplets splash rather than spread. Heat flux declines and spray cooling efficiency can markedly decrease. Accordingly we are investigating droplet impingement upon a) dry solid surfaces, b) fluid films, c) rough surfaces and determining splashing thresholds and relationships for both dry surfaces and those covered by fluid films. We are presently developing engineering correlations delineating the boundary between splashing and non-splashing regions. Determining the splash/non-splash boundary is important for many practical applications. Coating and cooling processes would each benefit from near-term empirical relations and subsequent models. Such demonstrations can guide theoretical development by providing definitive testing of its predictive capabilities. Thus, empirical relations describing the boundary between splash and non-splash are given for drops impinging upon a dry solid surface and upon a thin fluid film covering a similar surface. Analytical simplification of the power laws describing the boundary between the splash and non-splash regions yields insight into the engineering parameters governing the splash and non-splash outcomes of the fluid droplets. The power law correlation is shown separating the splashing versus non-splashing regions as developed for droplets impinging upon a dry solid surface. Splashing upon a dry surface is reasonably described by Ca greater than 0.85, reflecting the competing roles of surface tension and viscosity. The power law correlation is shown separating the splashing versus non-splashing regions as developed for droplets impinging upon a thin fluid film covering the solid surface. Splashing upon a thin fluid film, as described by v (pd/s) greater than 63, is governed by fluid density and surface tension, but is rather independent of viscosity. Finally, the data presented here suggests that a more direct dependence upon the surface tension and viscosity, given a better understanding of their interplay, would allow accurate description of the droplet-surface impacts for more complicated situations involving non-Newtonian fluids, specifically those exhibiting viscoelastic behavior.

Description: The dynamics and stability of systems of interfaces is central to a range of technologies related to the Human Exploration and Development of Space (HEDS). Our premise is that dramatic shape changes can be manipulated to advantage with minimal input, if the system is near instability. The primary objective is to develop the science base to allow novel approaches to liquid management in low-gravity based on this premise. HEDS requires efficient, reliable and lightweight technologies. Our poster will highlight our progress toward this goal using the capillary switch as an example. A capillary surface is a liquid/liquid or liquid/gas interface whose shape is determined by surface tension. For typical liquids (e.g., water) against gas on earth, capillary surfaces occur on the millimeterscale and smaller where shape deformation due to gravity is unimportant. In low gravity, they can occur on the centimeter scale. Capillary surfaces can be combined to make a switch a system with multiple stable states. A capillary switch can generate motion or effect force. To be practical, the energy barriers of such a switch must be tunable, its switching time (kinetics) short and its triggering mechanism reliable. We illustrate these features with a capillary switch that consists of two droplets, coupled by common pressure. As long as contact lines remained pinned, motions are inviscid, even at sub-millimeter scales, with consequent promise of low-power consumption at the device level. Predictions of theory are compared to experiment on i) a soap-film prototype at centimeter scale and ii) a liquid droplet switch at millimeter-scale.

Description: The present paper reports ongoing work to develop numerical and modeling tools used to design efficient and effective spray cooling processes and to determine characteristic non-dimensional parametric dependence for practical fluids and conditions. In particular, we present data that will delineate conditions towards control of the impingement dynamics of droplets upon a heated substrate germane to practical situations.

Description: A series of unsteady Reynolds-averaged Navier-Stokes computations are performed for the flow of a synthetic jet issuing into a turbulent boundary layer through a circular orifice. This is one of the validation test cases from a synthetic jet validation workshop held in March 2004. Several numerical parameters are investigated, and the effects of three different turbulence models are explored. Both long-time-averaged and time-dependent phase-averaged results are compared to experiment. On the whole, qualitative comparisons of the mean flow quantities are fairly good. There are many differences evident in the quantitative comparisons. The calculations do not exhibit a strong dependence on the type of turbulence model employed.

Description: One of the greatest uncertainties affecting the design of multiphase flow technologies for space exploration is the spatial distribution of phases that will arise in microgravity or reduced gravity. On Earth, buoyancy-driven motion predominates whereas the shearing of the bubble suspension controls its behavior in microgravity. We are conducting a series of ground-based experiments and a flight experiment spanning the full range of ratios of buoyancy to shear. These include: (1) bubbles rising in a quiescent liquid in a vertical channel; (2) weak shear flow induced by slightly inclining the channel; (3) moderate shear flow in a terrestrial vertical pipe flow; and (4) shearing of a bubble suspension in a cylindrical Couette cell in microgravity. We consider nearly monodisperse suspensions of 1 to 1.8 mm diameter bubbles in aqueous electrolyte solutions. The liquid velocity disturbance produced by bubbles in this size range can often be described using an inviscid analysis. Electrolytic solutions lead to hydrophilic repulsion forces that stabilize the bubble suspension without causing Marangoni stresses. We will discuss the mechanisms that control the flow behavior and phase distribution in the ground-based experiments and speculate on the factors that may influence the suspension flow and bubble volume fraction distribution in the flight experiment.

Description: Man has strived to make objects fly faster, first from subsonic to supersonic and then to hypersonic speeds. Spacecraft and high-speed missiles routinely fly at hypersonic Mach numbers, M greater than 5. In defense applications, aircraft reach hypersonic speeds at high altitude and so may civilian aircraft in the future. Hypersonic flight, while presenting opportunities, has formidable challenges that have spurred vigorous research and development, mainly by NASA and the Air Force in the USA. Although NASP, the premier hypersonic concept of the eighties and early nineties, did not lead to flight demonstration, much basic research and technology development was possible. There is renewed interest in supersonic and hypersonic flight with the HyTech program of the Air Force and the Hyper-X program at NASA being examples of current thrusts in the field. At high-subsonic to supersonic speeds, fluid compressibility becomes increasingly important in the turbulent boundary layers and shear layers associated with the flow around aerospace vehicles. Changes in thermodynamic variables: density, temperature and pressure, interact strongly with the underlying vortical, turbulent flow. The ensuing changes to the flow may be qualitative such as shocks which have no incompressible counterpart, or quantitative such as the reduction of skin friction with Mach number, large heat transfer rates due to viscous heating, and the dramatic reduction of fuel/oxidant mixing at high convective Mach number. The peculiarities of compressible turbulence, so-called compressibility effects, have been reviewed by Fernholz and Finley. Predictions of aerodynamic performance in high-speed applications require accurate computational modeling of these "compressibility effects" on turbulence. During the course of the project we have made fundamental advances in modeling the pressure-strain correlation and developed a code to evaluate alternate turbulence models in the compressible shear layer.

Description: An upwind scheme is presented for solving the three-dimensional Euler equations on unstructured tetrahedral meshes. Spatial discretization is accomplished by a cell-centered finite-volume formulation using flux-difference splitting. Higher-order differences are formed by a novel cell reconstruction process which results in computational times per cell comparable to those of structured codes. The approach yields highly resolved solutions in regions of smooth flow while avoiding oscillations across shocks without explicit limiting. Solutions are advanced in time by a 3-stage Runge-Kutta time-stepping scheme with convergence accelerated to steady state by local time stepping and implicit residual smoothing. Solutions are presented for a range of configurations in the transonic speed regime to demonstrate code accuracy, speed, and robustness. The results include an assessment of grid sensitivity and convergence acceleration by mesh sequencing.

Description: For DES, RANS is responsible for predicting boundary layer growth and separation. LES is responsible for predicting the geometry dependent turbulent flow features. Grid adaptation done using NASA Langley s RefineMesh program. Adaptation on time average of vorticity

Description: Flow cytometry is a powerful technique for obtaining quantitative information from fluorescence in cells. Quantitation is achieved by assuring a high degree of uniformity in the optical excitation and detection, generally by using a highly controlled flow such as is obtained via hydrodynamic focusing. In this work, we demonstrate a two-beam, two- channel detection and two-photon excitation flow cytometry (T(sup 3)FC) system that enables multi-dye analysis to be performed very simply, with greatly relaxed requirements on the fluid flow. Two-photon excitation using a femtosecond near-infrared (NIR) laser has the advantages that it enables simultaneous excitation of multiple dyes and achieves very high signal-to-noise ratio through simplified filtering and fluorescence background reduction. By matching the excitation volume to the size of a cell, single-cell detection is ensured. Labeling of cells by targeted nanoparticles with multiple fluorophores enables normalization of the fluorescence signal and thus ratiometric measurements under nonuniform excitation. Quantitative size measurements can also be done even under conditions of nonuniform flow via a two-beam layout. This innovative detection scheme not only considerably simplifies the fluid flow system and the excitation and collection optics, it opens the way to quantitative cytometry in simple and compact microfluidics systems, or in vivo. Real-time detection of fluorescent microbeads in the vasculature of mouse ear demonstrates the ability to do flow cytometry in vivo. The conditions required to perform quantitative in vivo cytometry on labeled cells will be presented.

Description: Work performed over the last three years has resulted in the addition of several new algorithms to the VULCAN code, NASA's standard for Navier-Stokes calculations in high-speed aeropropulsion devices. This final report describes the new techniques in brief and presents sample results from their use.

Description: The performance of five thermal mass flow meters, MKS Instruments 179A and 258C, Unit Instruments UFM-8100, Sierra Instruments 830L, and Hastings Instruments HFM-200, were tested on the KC-135 Reduced Gravity Aircraft in orthogonal, coparallel, and counterparallel orientations relative to gravity. Data was taken throughout the parabolic trajectory where the g-level varied from 0.01 to 1.8 times normal gravity. Each meter was calibrated in normal gravity in the orthogonal position prior to flight followed by ground testing at seven different flow conditions to establish a baseline operation. During the tests, the actual flow rate was measured independently using choked-flow orifices. Gravitational acceleration and attitude had a unique effect on the performance of each meter. All meters operated within acceptable limits at all gravity levels in the calibrated orthogonal position. However, when operated in other orientations, the deviations from the reference flow became substantial for several of the flow meters. Data analysis indicated that the greatest source of error was the effect of orientation, followed by the gravity level. This work emphasized that when operating thermal flow meters in a variable gravity environment, it is critical to orient the meter in the same direction relative to gravity in which it was calibrated. Unfortunately, there was no test in normal gravity that could predict the performance of a meter in reduced gravity. When operating in reduced gravity, all meters indicated within 5 percent of the full scale reading at all flow conditions and orientations.

Description: This slide presentation reviews the transport and decay of wake vortices in ground effect and cites a need for a physics-based parametric model. The encounter of a vortex with a solid body is always a complex event involving turbulence enhancement, unsteadiness, and very large gradients of velocity and pressure. Wake counter in ground effect is the most dangerous of them all. The interaction of diverging, area-varying, and decaying aircraft wake vortices with the ground is very complex because both the vortices and the flow field generated by them are altered to accommodate the presence of the ground (where there is very little room to maneuver) and the background turbulent flow. Previous research regarding vortex models, wake vortex decay mechanisms, time evolution within in ground effect of a wake vortex pair, laminar flow in ground effect, and the interaction of the existing boundary layer with a convected vortex are reviewed. Additionally, numerical simulations, 3-dimensional large-eddy simulations, a probabilistic 2-phase wake vortex decay and transport model and a vortex element method are discussed. The devising of physics-based, parametric models for the prediction of (operational) real-time response, mindful of the highly three-dimensional and unsteady structure of vortices, boundary layers, atmospheric thermodynamics, and weather convective phenomena is required. In creating a model, LES and field data will be the most powerful tools.

Description: This is the report of a Scientific Working Group (SWG) formed by NASA to determine the feasibility of using a liquid metal cooled nuclear reactor and Rankine energy conversion cycle for dual purpose power and propulsion in space. This is a high level technical report which is intended for use by NASA management in program planning. The SWG was composed of a team of specialists in nuclear energy and multiphase flow and heat transfer technology from academia, national laboratories, NASA and industry. The SWG has identified the key technology issues that need to be addressed and have recommended an integrated short term (approx. 2 years) and a long term (approx. 10 year) research and development (R&D) program to qualify a Rankine cycle power plant for use in space. This research is ultimately intended to give NASA and its contractors the ability to reliably predict both steady and transient multiphase flow and heat transfer phenomena at reduced gravity, so they can analyze and optimize designs and scale-up experimental data on Rankine cycle components and systems. In addition, some of these results should also be useful for the analysis and design of various multiphase life support and thermal management systems being considered by NASA.

Description: An analytic model for predicting the effect of unsteady local surface injection on the flow separating from a streamlined body at angle of attack is proposed. The model uses the premise that separation control results from enhanced mixing along the shear layer that develops between the main stream and the fluid in the underlying recirculation zone. High-Reynolds-number asymptotic methods are used to connect the unsteady surface injection to an instability wave propagating on the separating shear layer and then to the large-scale coherent structures that produce the increased mixing. The results is a tool that can guide the choice of fluid-actuator parameters to maximize flow-control effectiveness and may also facilitate computer-based numerical experiments.

Description: The safe and reliable operation of high pressure test stands for rocket engine and component testing places an increased emphasis on the performance of control valves and flow metering devices. In this paper, we will present a series of high fidelity computational analyses of systems ranging from cryogenic control valves and pressure regulator systems to cavitating venturis that are used to support rocket engine and component testing at NASA Stennis Space Center. A generalized multi-element framework with sub-models for grid adaption, grid movement and multi-phase flow dynamics has been used to carry out the simulations. Such a framework provides the flexibility of resolving the structural and functional complexities that are typically associated with valve-based high pressure feed systems and have been difficult to deal with traditional CFD methods. Our simulations revealed a rich variety of flow phenomena such as secondary flow patterns, hydrodynamic instabilities, fluctuating vapor pockets etc. In the paper, we will discuss performance losses related to cryogenic control valves, and provide insight into the physics of the dominant multi-phase fluid transport phenomena that are responsible for the choking like behavior in cryogenic control elements. Additionally, we will provide detailed analyses of the modal instability that is observed in the operation of the dome pressure regulator valve. Such instabilities are usually not localized and manifest themselves as a system wide phenomena leading to an undesirable chatter at high flow conditions.

Description: A two-dimensional multi-block topology generation technique has been developed. Very general configurations are addressable by the technique. A configuration is defined by a collection of non-intersecting closed curves, which will be referred to as loops. More than a single loop implies that holes exist in the domain, which poses no problem. This technique requires only the medial vertices and the touch points that define each vertex. From the information about the medial vertices, the connectivity between medial vertices is generated. The physical shape of the medial edge is not required. By applying a few simple rules to each medial edge, the multiblock topology is generated with no user intervention required. The resulting topologies contain only the level of complexity dictated by the configurations. Grid lines remain attached to the boundary except at sharp concave turns where a change in index family is introduced as would be desired. Keeping grid lines attached to the boundary is especially important in the area of computational fluid dynamics where highly clustered grids are used near no-slip boundaries. This technique is simple and robust and can easily be incorporated into the overall grid generation process.

Description: In this paper we consider the robust control of a thermal mixer using multivariable Sliding Mode Control (SMC). The mixer consists of a mixing chamber, hot and cold fluid valves, and an exit valve. The commanded positions of the three valves are the available control inputs, while the controlled variables are total mass flow rate, chamber pressure and the density of the mixture inside the chamber. Unsteady thermodynamics and linear valve models are used in deriving a 5th order nonlinear system with three inputs and three outputs, An SMC controller is designed to achieve robust output tracking in the presence of unknown energy losses between the chamber and the environment. The usefulness of the technique is illustrated with a simulation.

Description: This paper provides a discussion of the history of Carbon Cloth Phenolic (CCP) ply lifting in the Redesigned Solid Rocket Motor (RSRM) Program, a brief presentation of theoretical methods used for analytical evaluation, and results of parametric analyses of CCP material subject to test conditions of the Laser Hardened Material Evaluation Laboratory. CCP ply lift can occur in regions of the RSRM nozzle where ply angle to flame surface is generally less than about 20 degrees. There is a heat rate dependence on likelihood and severity of the condition with the higher heating rates generally producing more ply lift. The event occurs in-depth, near the heated surface, where the load necessary to mechanically separate the CCP plies is produced by the initial stages of pyrolysis gas generation due to the thermal decomposition of the phenolic resin matrix. Due to the shallow lay-up angle of the composite, normal components of the indepth mechanical load, due to "pore pressure", are imparted primarily as a cross-ply tensile force on the interlaminar ply boundaries. Tensile capability in the cross-ply (out of plane) direction is solely determined by the matrix material capability. The elevated temperature matrix material capabilities are overcome by pressure induced mechanical normal stress and ply-lift occurs. A theoretical model used for CCP in-depth temperature, pressure, and normal stress prediction, based on first principles, is briefly discussed followed by a parametric evaluation of response variables subject to boundary conditions typical of on-going test programs at the LHMEL facility. Model response demonstrates general trends observed in test and provides insight into the interactivity of material properties and constitutive relationships.

Description: From 1994 to 1996, NASA s Marshall Space Flight Center conducted a Center Director's Discretionary Fund research effort to apply artificial intelligence technologies to the health management of plant equipment and space propulsion systems. Through this effort, NASA established a business relationship with Quality Monitoring and Control (QMC), of Kingwood, Texas, to provide hardware modeling and artificial intelligence tools. Very detailed and accurate Space Shuttle Main Engine (SSME) analysis and algorithms were jointly created, which identified several missing, critical instrumentation needs for adequately evaluating the engine health status. One of the missing instruments was a liquid oxygen (LOX) flow measurement. This instrument was missing since the original SSME included a LOX turbine flow meter that failed during a ground test, resulting in considerable damage for NASA. New balanced flow meter technology addresses this need with robust, safe, and accurate flow metering hardware.

Description: A system and method uses a wireless tether comprising a transmitter and a receiver to alert a caregiver that an object or person has been left unattended. A detector Senses the presence of the object, usually a child, located in a position such as a safety seat. The detector couples to the transmitter, which is located near the object. The transmitter transmits at least one wireless signal when the object is in the position. The receiver, which is remotely located from the transmitter, senses the at least one signal as long as the receiver is within a prescribed range of transmission. By performing a timing function, the receiver monitors the proximity of the caregiver, who maintains possession of the receiver, to the transmitter. The system communicates an alarm to the caregiver when the caregiver ventures outside the range of transmission without having removed the object/child from the position.

Description: A test apparatus and method of its use for evaluating various performance aspects of a piping segment locates a piping segment between two cold boxes. A first cold box conditions test fluid before providing the fluid into the piping segment- The first and second cold boxes both significantly reduce, if not eliminate, any heat transfer from the ends of the piping so that accurate measurements of heat leak rates from the sides of the piping segment may be determined.

Description: The use of detailed CFD modeling for the description of cooling in rocket chambers is discussed. The overall analysis includes a complete three-dimensional analysis of the flow in the regenerative cooling passages, conjugate heat transfer in the combustor walls, and the effects of film cooling on the inside chamber. The results in the present paper omit the effects of film cooling and include only regen cooling and the companion conjugate heat transfer. The hot combustion gases are replaced by a constant temperature wall boundary condition. Load balancing for parallel cluster computations is ensured by using single-block unstructured grids for both fluids and solids, and by using a 'multiple physical zones' to account for differences in the number of equations. Validation of the method is achieved by comparing simple two-dimensional solutions with analytical results. Representative results for cooling passages are presents showing the effects of heat conduction in the copper walls with tube aspect ratios of 1.5:l.

Description: The main objective of this work is to determine the unsteady hydrodynamic characteristics of coaxial swirl atomizers of interest in oxidizer-rich staged combustion (ORSC) liquid rocket engines. To this end, the pseudo-density (homogeneous flow) treatment combined with the Marker-and-Cell (MAC) numerical algorithm has been used to develop an axisymmetric with swirl, two-phase, unsteady model. The numerical model is capable of assessing the time-dependent orifice exit conditions and internal mixing for arbitrary fuel and oxidizer gas injection conditions. Parametric studies have been conducted to determine the effect of geometry, gas properties, and liquid properties on the exit massflow rate and velocity. It has been found that the frequency at which the liquid film oscillates increases as the density ratio and thickness increase, decreases as film thickness and liquid swirl velocity increase, and is unaffected by the mixing length. Additionally, it has been determined that the variation in the massflow rate increases as the liquid swirl velocity and liquid film thickness increase, and decreases as the density ratio, collar thickness, and mixing length increase.

Description: The two most common approaches used to formulate thermodynamic properties of pure substances are fundamental (or characteristic) equations of state (Helmholtz and Gibbs functions) and a piecemeal approach that is described in Adebiyi and Russell (1992). This paper neither presents a different method to formulate thermodynamic properties of pure substances nor validates the aforementioned approaches. Rather its purpose is to present a method to generate property tables from existing property packages and a method to facilitate the accurate interpretation of fluid thermodynamic property data from those tables. There are two parts to this paper. The first part of the paper shows how efficient and usable property tables were generated, with the minimum number of data points, using an aerospace industry standard property package. The second part describes an innovative interpolation technique that has been developed to properly obtain thermodynamic properties near the saturated liquid and saturated vapor lines.

Description: Multiscale modeling of materials requires simulations of multiple levels of structural hierarchy. The computational efficiency of numerical methods becomes a critical factor for simulating large physical systems with highly desperate length scales. Multigrid methods are known for their superior efficiency in representing/resolving different levels of physical details. The efficiency is achieved by employing interactively different discretizations on different scales (grids). To assist optimization of manufacturing conditions for materials processing with numerous particles (e.g., dispersion of particles, controlling flow viscosity and clusters), a new multigrid algorithm has been developed for a case of multiscale modeling of flows with small particles that have various length scales. The optimal efficiency of the algorithm is crucial for accurate predictions of the effect of processing conditions (e.g., pressure and velocity gradients) on the local flow fields that control the formation of various microstructures or clusters.

Description: An unsteady, three dimensional Navier-Stokes solution in rotating frame formulation for turbomachinery applications is presented. Casting the governing equations in a rotating frame enabled the freezing of grid motion and resulted in substantial savings in computer time. The turbine blade was computationally simulated and probabilistically evaluated in view of several uncertainties in the aerodynamic, structural, material and thermal variables that govern the turbine blade. The interconnection between the computational fluid dynamics code and finite element structural analysis code was necessary to couple the thermal profiles with the structural design. The stresses and their variations were evaluated at critical points on the Turbine blade. Cumulative distribution functions and sensitivity factors were computed for stress responses due to aerodynamic, geometric, mechanical and thermal random variables.

Description: This report consists of two abstracts submitted for possible presentation at the AIAA Aerospace Science Meeting to be held in January 2005. Since the submittal of these abstracts we are continuing refinement of the model coefficients derived for the case of a variable Turbulent Prandtl number. The test cases being investigated are a Mach 9.2 flow over a degree ramp and a Mach 8.2 3-D calculation of crossing shocks. We have developed an axisymmetric code for treating axisymmetric flows. In addition the variable Schmidt number formulation was incorporated in the code and we are in the process of determining the model constants.

Description: In 1999, Stolz and Adams unveiled a subgrid-scale model for LES based upon approximately inverting (defiltering) the spatial grid-filter operator and termed .the approximate deconvolution model (ADM). Subsequently, the utility and accuracy of the ADM were demonstrated in a posteriori analyses of flows as diverse as incompressible plane-channel flow and supersonic compression-ramp flow. In a prelude to the current paper, a parameterized temporal ADM (TADM) was developed and demonstrated in both a priori and a posteriori analyses for forced, viscous Burger's flow. The development of a time-filtered variant of the ADM was motivated-primarily by the desire for a unifying theoretical and computational context to encompass direct numerical simulation (DNS), large-eddy simulation (LES), and Reynolds averaged Navier-Stokes simulation (RANS). The resultant methodology was termed temporal LES (TLES). To permit exploration of the parameter space, however, previous analyses of the TADM were restricted to Burger's flow, and it has remained to demonstrate the TADM and TLES methodology for three-dimensional flow. For several reasons, plane-channel flow presents an ideal test case for the TADM. Among these reasons, channel flow is anisotropic, yet it lends itself to highly efficient and accurate spectral numerical methods. Moreover, channel-flow has been investigated extensively by DNS, and a highly accurate data base of Moser et.al. exists. In the present paper, we develop a fully anisotropic TADM model and demonstrate its utility in simulating incompressible plane-channel flow at nominal values of Re(sub tau) = 180 and Re(sub tau) = 590 by the TLES method. The TADM model is shown to perform nearly as well as the ADM at equivalent resolution, thereby establishing TLES as a viable alternative to LES. Moreover, as the current model is suboptimal is some respects, there is considerable room to improve TLES.

Description: Work performed over the last three years has resulted in the addition of several new algorithms to the VULCAN code, NASA's standard for Navier-Stokes calculations in high-speed aeropropulsion devices. This final report describes the new techniques in brief and presents sample results from their use.

Description: Successful use of CFD to provide aerodynamics for stability and control (S&C) applications will require that the traditional time and costs associated with CFD be reduced and that the errors and uncertainties currently associated with CFD be better understood. CFD will be required to work under a wide range of flow conditions and provide fast and reliable aerodynamics if it is to contribute to this next generation of S&C analyses. CFD solutions have errors and uncertainties due to poorly converged solutions, solution anomalies caused by grids, turbulence models, and parameter selection, and other manifold reasons. In addition to the above problems, there will be a requirement for communications between the CFD expert and the S&C expert and possibly experts from other related disciplines. The CFD expert may not understand the technical problems associated with S&C, and it is almost certain the converse is true.

Description: Analysis of high spectral resolution observations of the lambda6614 DIB line profile show systematic variations in the positions of the peaks in the substructure of the profile. These variations can only be understood in the framework of rotational contours of large molecules, where the variations are caused by changes in the rotational excitation temperature. We show that the rotational excitation temperature for the DIB carrier is of the order 10-40 K - much lower than the gas kinetic temperature - indicating that for this particular DIB carrier angular momentum buildup is not very efficient. The rotational constant indicates that the carrier of this DIB is smaller than previously assumed:7-22 C atoms, depending on the geometry.

Description: Computational Fluid Dynamics (CFD) simulations are routinely performed as part of the design process of most fluid handling devices. In order to efficiently and effectively use the results of a CFD simulation, visualization tools are often used. These tools are used in all stages of the CFD simulation including pre-processing, interim-processing, and post-processing, to interpret the results. Each of these stages requires visualization tools that allow one to examine the geometry of the device, as well as the partial or final results of the simulation. An engineer will typically generate a series of contour and vector plots to better understand the physics of how the fluid is interacting with the physical device. Of particular interest are detecting features such as shocks, recirculation zones, and vortices (which will highlight areas of stress and loss). As the demand for CFD analyses continues to increase the need for automated feature extraction capabilities has become vital. In the past, feature extraction and identification were interesting concepts, but not required in understanding the physics of a steady flow field. This is because the results of the more traditional tools like; iso-surface, cuts and streamlines, were more interactive and easily abstracted so they could be represented to the investigator. These tools worked and properly conveyed the collected information at the expense of a great deal of interaction. For unsteady flow-fields, the investigator does not have the luxury of spending time scanning only one "snapshot" of the simulation. Automated assistance is required in pointing out areas of potential interest contained within the flow. This must not require a heavy compute burden (the visualization should not significantly slow down the solution procedure for (co-processing environments). Methods must be developed to abstract the feature of interest and display it in a manner that physically makes sense.

Description: Non-coalescence of two bodies of the same liquid and the suppression of contact between liquid drops and solid surfaces (nonwetting) has been studied through a pair of parallel investigations being conducted at the Georgia Institute of Technology and the Microgravity Research and Support (MARS) Center in Naples, Italy. Both non- coalescence and nonwetting may achieved by either: i) exploiting the mechanism of thermocapillary convection to drive a lubricating film of surrounding gas (air) into the space between the two liquid free surfaces (non-coalescence) or between the drop free surface and the solid (nonwetting); or ii) using the otherwise forced motion of the solid (or one liquid) surface to generate the lubricating film. These, and other types of related phenomena are discussed in our recent review article (Neitzel GP, Dell'Aversana P. 2002. Noncoalescence and nonwetting behavior of liquids. Annual Review of Fluid Mechanics 34: 267-89). This project was a continuation of earlier research on this subject and includes the possible uses of non-coalescing/wetting systems in technological applications.

Description: A control valve that can throttle high pressure cryogenic fluid embodies several design features that distinguish it over conventional valves designed for similar applications. Field and design engineers worked together to create a valve that would simplify installation, trim changes, and maintenance, thus reducing overall cost. The seals and plug stem packing were designed to perform optimally in cryogenic temperature ranges. Unlike conventional high-pressure cryogenic valves, the trim size can be changed independent of the body. The design feature that provides flexibility for changing the trim is a split body. The body is divided into an upper and a lower section with the seat ring sandwiched in between. In order to maintain the plug stem packing at an acceptable sealing temperature during cryogenic service, heat-exchanging fins were added to the upper body section (see figure). The body is made of stainless steel.

Description: We have employed a Hamiltonian model based on a self-consistent Gaussian appoximation to examine the unfolding process of proteins in external - both mechanical and thermal - force elds. The motivation was to investigate the unfolding pathways of proteins by including only the essence of the important interactions of the native-state topology. Furthermore, if such a model can indeed correctly predict the physics of protein unfolding, it can complement more computationally expensive simulations and theoretical work. The self-consistent Gaussian approximation by Micheletti et al. has been incorporated in our model to make the model mathematically tractable by signi cantly reducing the computational cost. All thermodynamic properties and pair contact probabilities are calculated by simply evaluating the values of a series of Incomplete Gamma functions in an iterative manner. We have compared our results to previous molecular dynamics simulation and experimental data for the mechanical unfolding of the giant muscle protein Titin (1TIT). Our model, especially in light of its simplicity and excellent agreement with experiment and simulation, demonstrates the basic physical elements necessary to capture the mechanism of protein unfolding in an external force field.

Description: A workshop was organized on the topic of the title and held on August 17-20, 2003 at the Syracuse University Minnowbrook Conference Center in Blue Mountain Lake, New York. Attendance was by invitation only, 47 guests attended and 30 presentations were made. Support was received from NASA Glenn Research Center, the US Air Force Office of Scientific Research, the European Office of Aeronautical Research and Development, the Asian Office of Aeronautical Research and Development and Syracuse University. This workshop was the fourth in a trienniel series beginning in 1993. A publication under a NASA CP 2004-212913 will be issued and include all abstracts. No full written papers were required. This report includes a list of attendees and the program of presentations. The next workshop is scheduled for August 20-23, 2006.

Description: The performance of a heat pipe system is greatly improved by the use of a dilute aqueous solution of about 0.0005 and about 0.005 moles per liter of a long chain alcohol as the working fluid. The surface tension-temperature gradient of the long-chain alcohol solutions turns positive as the temperature exceeds a certain value, for example about 40.degree. C. for n-heptanol solutions. Consequently, the Marangoni effect does not impede, but rather aids in bubble departure from the heating surface. Thus, the bubble size at departure is substantially reduced at higher frequencies and, therefore, increases the boiling limit of heat pipes. This feature is useful in microgravity conditions. In addition to microgravity applications, the heat pipe system may be used for commercial, residential and vehicular air conditioning systems, micro heat pipes for electronic devices, refrigeration and heat exchangers, and chemistry and cryogenics.

Description: EUPDF-II provides the solution for the species and temperature fields based on an evolution equation for PDF (Probability Density Function) and it is developed mainly for application with sprays, combustion, parallel computing, and unstructured grids. It is designed to be massively parallel and could easily be coupled with any existing gas-phase CFD and spray solvers. The solver accommodates the use of an unstructured mesh with mixed elements of either triangular, quadrilateral, and/or tetrahedral type. The manual provides the user with an understanding of the various models involved in the PDF formulation, its code structure and solution algorithm, and various other issues related to parallelization and its coupling with other solvers. The source code of EUPDF-II will be available with National Combustion Code (NCC) as a complete package.

Description: One of the power systems under consideration for nuclear electric propulsion or as a planetary surface power source is a heatpipe-cooled reactor coupled to a Brayton cycle. In this system, power is transferred from the heatpipes to the Brayton gas via a heat exchanger attached to the heatpipes. This paper discusses the fluid, thermal and structural analyses that were performed in support of the design of the heat exchanger to be tested in the SAFE-100 experimental program at the Marshall Space Flight Center: An important consideration throughout the design development of the heat exchanger w its capability to be utilized for higher power and temperature applications. This paper also discusses this aspect of the design and presents designs for specific applications that are under consideration.

Description: An overview of the data acquisition, reduction, and uncertainty of experimental measurements of the flowfield created by the interaction of an isolated synthetic jet and a turbulent boundary layer is presented. The experimental measurements were undertaken to serve as the second of three computational fluid dynamics validation databases for Active Flow Control. The validation databases were presented at the NASA Langley Research Center Workshop on CFD Validation of Synthetic Jets and Turbulent Separation Control in March, 2004. Detailed measurements were made to document the boundary conditions for the flow and also for the phase-averaged flowfield itself. Three component Laser-Doppler Velocimetry, 2-D Particle Image Velocimetry, and Stereo Particle Image Velocimetry were utilized to document the phase averaged velocity field and the turbulent stresses.

Description: An approach to predicting the effects of free stream turbulence on turbine vane and blade heat transfer is described. Four models for predicting the effects of free stream turbulence were in incorporated into a Navier-Stokes CFD analysis. Predictions were compared with experimental data in order to identify an appropriate model for use across a wide range of flow conditions. The analyses were compared with data from five vane geometries and from four rotor geometries. Each of these nine geometries had data for different Reynolds numbers. Comparisons were made for twenty four cases. Steady state calculations were done because all experimental data were obtained in steady state tests. High turbulence levels often result in suction surface transition upstream of the throat, while at low to moderate Reynolds numbers the pressure surface remains laminar. A two-dimensional analysis was used because the flow is predominately two-dimensional in the regions where free stream turbulence significantly augments surface heat transfer. Because the evaluation of models for predicting turbulence effects can be affected by other factors, the paper discusses modeling for transition, relaminarization, and near wall damping. Quantitative comparisons are given between the predictions and data.

Description: Heated jets in a wide range of temperature ratios (TR), and acoustic Mach numbers (Ma) were investigated experimentally using far field microphones and a molecular Rayleigh scattering technique. The latter provided density fluctuations measurements. Two sets of operating conditions were considered: (1) TR was varied between 0.84 and 2.7 while Ma was fixed at 0.9; (2) Ma was varied between 0.6 and 1.48, while TR was fixed at 2.27. The implementation of the molecular Rayleigh scattering technique required dust removal and usage of a hydrogen combustor to avoid soot particles. Time averaged density measurements in the first set of data showed differences in the peripheral density shear layers between the unheated and heated jets. The nozzle exit shear layer showed increased turbulence level with increased plume temperature. Nevertheless, further downstream the density fluctuations spectra are found to be nearly identical for all Mach number and temperature ratio conditions. To determine noise sources a correlation study between plume density fluctuations and far field sound pressure fluctuations was conducted. For all jets the core region beyond the end of the potential flow was found to be the strongest noise source. Except for an isothermal jet, the correlations did not differ significantly with increasing temperature ratio. The isothermal jet created little density fluctuations. Although the far field noise from this jet did not show any exceptional trend, the flow-sound correlations were very low. This indicated that the density fluctuations only acted as a "tracer parameter" for the noise sources.

Description: In order to enhance the fundamental understanding of thin film evaporation and thereby improve the critical design concept for two-phase heat transfer devices, microscale heat and mass transport is to be investigated for the transition film region using state-of-the-art optical diagnostic techniques. By utilizing a microgravity environment, the length scales of the transition film region can be extended sufficiently, from submicron to micron, to probe and measure the microscale transport fields which are affected by intermolecular forces. Extension of the thin film dimensions under microgravity will be achieved by using a conical evaporator made of a thin silicon substrate under which concentric and individually controlled micro-heaters are vapor-deposited to maintain either a constant surface temperature or a controlled temperature variation. Local heat transfer rates, required to maintain the desired wall temperature boundary condition, will be measured and recorded by the concentric thermoresistance heaters controlled by a Wheatstone bridge circuit, The proposed experiment employs a novel technique to maintain a constant liquid volume and liquid pressure in the capillary region of the evaporating meniscus so as to maintain quasi-stationary conditions during measurements on the transition film region. Alternating use of Fizeau interferometry via white and monochromatic light sources will measure the thin film slope and thickness variation, respectively. Molecular Fluorescence Tracking Velocimetry (MFTV), utilizing caged fluorophores of approximately 10-nm in size as seeding particles, will be used to measure the velocity profiles in the thin film region. An optical sectioning technique using confocal microscopy will allow submicron depthwise resolution for the velocity measurements within the film for thicknesses on the order of a few microns. Digital analysis of the fluorescence image-displacement PDFs, as described in the main proposal, can further enhance the depthwise resolution.

Description: An implementation of premixed equilibrium chemistry has been completed for the OVERFLOW code, a chimera capable, complex geometry flow code widely used to predict transonic flowfields. The implementation builds on the computational efficiency and geometric generality of the solver.

Description: The major accomplishments of the experimental portion of the research were documented in Ling Zheng's doctoral dissertation. Using Pentane, he obtained a considerable amount of data on the stability and heat transfer characteristics of an evaporating meniscus. The important points are that experimental equipment to obtain data on the stability and heat transfer characteristics of an evaporating meniscus were built and successfully operated. The data and subsequent analyses were accepted by the Journal of Heat Transfer for publication in 2004 [PU4]. The work was continued by a new graduate student using HFE-7000 [PU3] and then Pentane at lower heat fluxes. The Pentane results are being analyzed for publication. The experimental techniques are currently being used in our other NASA Grant. The oscillation of the contact line observed in the experiments involves evaporation (retraction part) and spreading. Since both processes occur with finite contact angles, it is important to derive a precise equation of the intermolecular forces (disjoining pressure) valid for non-zero contact angles. This theoretical derivation was accepted for publication by Journal of Fluid Mechanics [PU5]. The evaporation process near the contact line is complicated, and an idealized micro heat pipe has been proposed to help in elucidating the detailed evaporation process [manuscripts in preparation].

Description: This work represents preliminary thermal stability results for liquid hydrocarbon fuels. High Reynolds Number Thermal Stability experiments with Jet A and RP-1 resulted in a quantitative measurement of the thermal stability. Each fuel flowed through a heated capillary tube that held the outlet temperature at 290 C. An optical pyrometer measured the surface temperature of the tube at 12 locations as a function of time. The High Reynolds Number Thermal Stability number was then determined using standards published by the American Society for Testing and Materials. The results for Jet A showed lower thermal stability than similar tests conducted at another facility. The RP-1 results are the first reported using this technique. Because the temperature rise on the capillary tube during testing for the RP-1 fuels was not significant, a new standard for the testing conditions should be developed for these types of fuels.

Description: The Hyper-X (X-43A) program is a flight experiment to demonstrate scramjet performance and operability under controlled powered free-flight conditions at Mach 7 and 10. The Mach 7 flight was successfully completed on March 27, 2004. Thermocouple instrumentation in the hot structures (nose, horizontal tail, and vertical tail) recorded the flight thermal response of these components. Preflight thermal analysis was performed for design and risk assessment purposes. This paper will present a comparison of the preflight thermal analysis and the recorded flight data.

Description: This Minnowbrook IV 2003 workshop on Transition and Unsteady Aspects of Turbomachinery Flows includes the following topics: 1) Current Issues in Unsteady Turbomachinery Flows; 2) Global Instability and Control of Low-Pressure Turbine Flows; 3) Influence of End Wall Leakage on Secondary Flow Development in Axial Turbines; 4) Active and Passive Flow Control on Low Pressure Turbine Airfoils; 5) Experimental and Numerical Investigation of Transitional Flows as Affected by Passing Wakes; 6) Effects of Freestream Turbulence on Turbine Blade Heat Transfer; 7) Bypass Transition Via Continuous Modes and Unsteady Effects on Film Cooling; 8) High Frequency Surface Heat Flux Imaging of Bypass Transition; 9) Skin Friction and Heat Flux Oscillations in Upstream Moving Wave Packets; 10) Transition Mechanisms and Use of Surface Roughness to Enhance the Benefits of Wake Passing in LP Turbines; 11) Transient Growth Approach to Roughness-Induced Transition; 12) Roughness- and Freestream-Turbulence-Induced Transient Growth as a Bypass Transition Mechanism; 13) Receptivity Calculations as a Means to Predicting Transition; 14) On Streamwise Vortices in a Curved Wall Jet and Their Effect on the Mean Flow; 15) Plasma Actuators for Separation Control of Low Pressure Turbine Blades; 16) Boundary-Layer Separation Control Under Low-Pressure-Turbine Conditions Using Glow-Discharge Plasma Actuators; 17) Control of Separation for Low Pressure Turbine Blades: Numerical Simulation; 18) Effects of Elevated Free-Stream Turbulence on Active Control of a Separation Bubble; 19) Wakes, Calming and Transition Under Strong Adverse Pressure Gradients; 20) Transitional Bubble in Periodic Flow Phase Shift; 21) Modelling Spots: The Calmed Region, Pressure Gradient Effects and Background; 22) Modeling of Unsteady Transitional Flow on Axial Compressor Blades; 23) Challenges in Predicting Component Efficiencies in Turbomachines With Low Reynolds Number Blading; 24) Observations on the Causal Relationship Between Blade Count and Developing Rotating Stall in a Four Stage Axial Compressor; 25) Experimental and Numerical Study of Non-Linear Interactions in Transonic Nozzle Flow; 26) Clocking Effects on a Modern Stage and One-Half Transonic Turbine; 27) DNS and LES of Transition on Turbine Blades; 28) The Use of Cellular Automata in Modeling the Transition; 29) Predicting Unsteady Buffet Onset Using RANS Solutions; 30) Transition Modelling With the SST Turbulence Model and an Intermittency Transport; and 31) Equation Workshop Summary Transcript

Description: Marshall Space Flight Center has developed and demonstrated a measurement device for sensing and resolving the hydrodynamic loads on fluid machinery. The device - a derivative of the six-component wind tunnel balance - senses the forces and moments on the rotating device through a weakened shaft section instrumented with a series of strain gauges. This rotating balance was designed to directly measure the steady and unsteady hydrodynamic loads on an inducer, thereby defining the amplitude and frequency content associated with operating in various cavitation modes. The rotating balance was calibrated statically using a dead-weight load system in order to generate the 6 x 12 calibration matrix later used to convert measured voltages to engineering units. Structural modeling suggested that the rotating assembly first bending mode would be significantly reduced with the balance s inclusion. This reduction in structural stiffness was later confirmed experimentally with a hammer-impact test. This effect, coupled with the relatively large damping associated with the rotating balance waterproofing material, limited the device s bandwidth to approximately 50 Hertz Other pre-test validations included sensing the test article rotating assembly built-in imbalance for two configurations and directly measuring the assembly mass and buoyancy while submerged under water. Both tests matched predictions and confirmed the device s sensitivity while stationary and rotating. The rotating balance was then demonstrated in a water test of a full-scale Space Shuttle Main Engine high-pressure liquid oxygen pump inducer. Experimental data was collected a scaled operating conditions at three flow coefficients across a range of cavitation numbers for the single inducer geometry and radial clearance. Two distinct cavitation modes were observed symmetric tip vortex cavitation and alternate-blade cavitation. Although previous experimental tests on the same inducer demonstrated two additional cavitation modes at lower inlet pressures, these conditions proved unreachable with the rotating balance installed due to the intense dynamic environment. The sensed radial load was less influenced by flow coefficient than by cavitation number or cavitation mode although the flow coefficient range was relatively narrow. Transition from symmetric tip vortex to alternate-blade cavitation corresponded to changes in both radial load magnitude and radial load orientation relative to the inducer. Sensed moments indicated that the effective load center moved downstream during this change in cavitation mode. An occurrence of "higher+rdex cavitation" was also detected in both the stationary pressures and the rotating balance data although the frequency of the phenomena was well above the reliable bandwidth of the rotating balance. In summary the experimental tests proved both the concept and device s capability despite the limitations and confirmed that hydrodynamically-induced forces and moments develop in response to the unbalanced pressure field, which is, in turn, a product of the cavitation environment.

Description: A baseline solution for CFD Point 1 (Mach 24) in the STS-107 accident investigation was modified to include effects of holes through the leading edge into a vented cavity. The simulations were generated relatively quickly and early in the investigation by making simplifications to the leading edge cavity geometry. These simplifications in the breach simulations enabled: 1) A very quick grid generation procedure; 2) High fidelity corroboration of jet physics with internal surface impingements ensuing from a breach through the leading edge, fully coupled to the external shock layer flow at flight conditions. These simulations provided early evidence that the flow through a 2 inch diameter (or larger) breach enters the cavity with significant retention of external flow directionality. A normal jet directed into the cavity was not an appropriate model for these conditions at CFD Point 1 (Mach 24). The breach diameters were of the same order or larger than the local, external boundary-layer thickness. High impingement heating and pressures on the downstream lip of the breach were computed. It is likely that hole shape would evolve as a slot cut in the direction of the external streamlines. In the case of the 6 inch diameter breach the boundary layer is fully ingested.

Description: This paper presents computational results obtained with the direct simulation Monte Carlo (DSMC) method for several biconic test cases in which shock interactions and flow separation-reattachment are key features of the flow. Recent ground-based experiments have been performed for several biconic configurations, and surface heating rate and pressure measurements have been proposed for code validation studies. The present focus is to expand on the current validating activities for a relatively new DSMC code called DS2V that Bird (second author) has developed. Comparisons with experiments and other computations help clarify the agreement currently being achieved between computations and experiments and to identify the range of measurement variability of the proposed validation data when benchmarked with respect to the current computations. For the test cases with significant vibrational nonequilibrium, the effect of the vibrational energy surface accommodation on heating and other quantities is demonstrated.

Description: An oscillatory zero net mass flow jet was generated by a cavity-pumping device, namely a synthetic jet actuator. This basic oscillating jet flow field was selected as the first of the three test cases for the Langley workshop on CFD Validation of Synthetic Jets and Turbulent Separation Control. The purpose of this workshop was to assess the current CFD capabilities to predict unsteady flow fields of synthetic jets and separation control. This paper describes the characteristics and flow field database of a synthetic jet in a quiescent fluid. In this experiment, Particle Image Velocimetry (PIV), Laser Doppler Velocimetry (LDV), and hot-wire anemometry were used to measure the jet velocity field. In addition, the actuator operating parameters including diaphragm displacement, internal cavity pressure, and internal cavity temperature were also documented to provide boundary conditions for CFD modeling.

Description: Computational Fluid Dynamics (CFD) is now routinely used to analyze isolated points in a design space by performing steady-state computations at fixed flight conditions (Mach number, angle of attack, sideslip), for a fixed geometric configuration of interest. This "point analysis" provides detailed information about the flowfield, which aides an engineer in understanding, or correcting, a design. A point analysis is typically performed using high fidelity methods at a handful of critical design points, e.g. a cruise or landing configuration, or a sample of points along a flight trajectory.

Description: The report summarizes the results of the redesign efforts directed towards the gas-turbine combustor rapid-injector flow diagnostic probe developed under sponsorship of NASA-GRC and earlier reported in NASA-CR-2003-212540. Lessons learned during the theoretical development, developmental testing and field-testing in the previous phase of this research were applied to redesign of both the probe sensing elements and of the rapid injection device. This redesigned probe (referred to herein as Turboprobe) has been fabricated and is ready, along with the new rapid injector, for field-testing. The probe is now designed to capture both time-resolved and mean total temperatures, total pressures and, indirectly, one component of turbulent fluctuations.

Description: A study of active-flow control in a small-scale boundary layer ingestion inlet was conducted at the NASA Langley Basic Aerodynamic Research Tunnel (BART). Forty MEMS pressure sensors, in a rake style configuration, were used to examine both the mean (DC) and high frequency (AC) components of the total pressure across the inlet/engine interface plane. The mean component was acquired and used to calculate pressure distortion. The AC component was acquired separately, at a high sampling rate, and is used to study the unsteady effects of the active-flow control. An identical total pressure rake, utilizing an Electronically Scanned Pressure (ESP) system, was also used to calculate distortion; a comparison of the results obtained using the two rakes is presented.

Description: A second international AIAA Drag Prediction Workshop (DPW-II) was organized and held in Orlando Florida on June 21-22, 2003. The primary purpose was to inves- tigate the code-to-code uncertainty. address the sensitivity of the drag prediction to grid size and quantify the uncertainty in predicting nacelle/pylon drag increments at a transonic cruise condition. This paper presents an in-depth analysis of the DPW-II computational results from three state-of-the-art unstructured grid Navier-Stokes flow solvers exercised on similar families of tetrahedral grids. The flow solvers are USM3D - a tetrahedral cell-centered upwind solver. FUN3D - a tetrahedral node-centered upwind solver, and NSU3D - a general element node-centered central-differenced solver. For the wingbody, the total drag predicted for a constant-lift transonic cruise condition showed a decrease in code-to-code variation with grid refinement as expected. For the same flight condition, the wing/body/nacelle/pylon total drag and the nacelle/pylon drag increment predicted showed an increase in code-to-code variation with grid refinement. Although the range in total drag for the wingbody fine grids was only 5 counts, a code-to-code comparison of surface pressures and surface restricted streamlines indicated that the three solvers were not all converging to the same flow solutions- different shock locations and separation patterns were evident. Similarly, the wing/body/nacelle/pylon solutions did not appear to be converging to the same flow solutions. Overall, grid refinement did not consistently improve the correlation with experimental data for either the wingbody or the wing/body/nacelle pylon configuration. Although the absolute values of total drag predicted by two of the solvers for the medium and fine grids did not compare well with the experiment, the incremental drag predictions were within plus or minus 3 counts of the experimental data. The correlation with experimental incremental drag was not significantly changed by specifying transition. Although the sources of code-to-code variation in force and moment predictions for the three unstructured grid codes have not yet been identified, the current study reinforces the necessity of applying multiple codes to the same application to assess uncertainty.

Description: This paper describes and analyzes a series of nearly 90 CFD test cases performed as a contribution to the second Drag Prediction Workshop, held in association with the AIAA in June 2003. Two configurations are included: DLR-F6 wing-body and wing-body-nacelle-pylon. The ability of CFD to predict the drag, lift, and pitching moment from experiment-including the "delta" arising from the addition of the nacelle and pylon-is assessed. In general, at a fixed angle of attack CFD overpredicts lift, but predicts the delta C (sub L) reasonably well. At low lift levels (C (sub L) less than 0.3)), delta C (sub D) is 20-30 drag counts (30-45%) high. At the target lift coefficient of C(sub L) = 0.5, delta C (sub D) is overpredicted by between 11-16 counts. However, the primary contribution of this paper is mot so much the assessment of CFD against experiment, but rather a detailed assessment and analysis of CFD variation. The series of test cases are designed to determine the sensitivity/variability of CFD to a variety of factors, including grid, turbulence model, transition code, and viscous model. Using medium-level grids (6-11 million points) at the target lift coefficient, the maximum variation in drag due to different grids is 5-11 drag counts, due to code is 5-10 counts, due to turbulence model is 7-15 counts, due to transition is 10-11 counts, and due to viscous model is 4-5 counts. Other specific variations are described in the paper.

Description: An agglomeration multigrid scheme has been implemented into the sequential version of the NASA code USM3Dns, tetrahedral cell-centered finite volume Euler/Navier-Stokes flow solver. Efficiency and robustness of the multigrid-enhanced flow solver have been assessed for three configurations assuming an inviscid flow and one configuration assuming a viscous fully turbulent flow. The inviscid studies include a transonic flow over the ONERA M6 wing and a generic business jet with flow-through nacelles and a low subsonic flow over a high-lift trapezoidal wing. The viscous case includes a fully turbulent flow over the RAE 2822 rectangular wing. The multigrid solutions converged with 12%-33% of the Central Processing Unit (CPU) time required by the solutions obtained without multigrid. For all of the inviscid cases, multigrid in conjunction with an explicit time-stepping scheme performed the best with regard to the run time memory and CPU time requirements. However, for the viscous case multigrid had to be used with an implicit backward Euler time-stepping scheme that increased the run time memory requirement by 22% as compared to the run made without multigrid.

Description: A heat pipe cooled reactor is one of several candidate reactor cores being considered for advanced space power and propulsion systems to support future space exploration applications. Long life heat pipe modules, with designs verified through a combination of theoretical analysis and experimental lifetime evaluations, would be necessary to establish the viability of any of these candidates, including the heat pipe reactor option. A hardware-based program was initiated to establish the infrastructure necessary to build heat pipe modules. This effort, initiated by Los Alamos National Laboratory and referred to as the Safe Affordable Fission Engine (SAFE) project, set out to fabricate and perform non-nuclear testing on a modular heat pipe reactor prototype that can provide 100 kilowatt from the core to an energy conversion system at 700 C. Prototypic heat pipe hardware was designed, fabricated, filled, closed-out and acceptance tested.

Description: A CFD validation workshop for synthetic jets and turbulent separation control (CFDVAL2004) was held in Williamsburg, Virginia in March 2004. Three cases were investigated: synthetic jet into quiescent air, synthetic jet into a turbulent boundary layer crossflow, and flow over a hump model with no-flow-control, steady suction, and oscillatory control. This paper is a summary of the CFD results from the workshop. Although some detailed results are shown, mostly a broad viewpoint is taken, and the CFD state-of-the-art for predicting these types of flows is evaluated from a general point of view. Overall, for synthetic jets, CFD can only qualitatively predict the flow physics, but there is some uncertainty regarding how to best model the unsteady boundary conditions from the experiment consistently. As a result. there is wide variation among CFD results. For the hump flow, CFD as a whole is capable of predicting many of the particulars of this flow provided that tunnel blockage is accounted for, but the length of the separated region compared to experimental results is consistently overpredicted.

Description: The use of differential equations such as Eikonal, Hamilton-Jacobi and Poisson for the economical calculation of the nearest wall distance d, which is needed by some turbulence models, is explored. Modifications that could palliate some turbulence-modeling anomalies are also discussed. Economy is of especial value for deforming/adaptive grid problems. For these, ideally, d is repeatedly computed. It is shown that the Eikonal and Hamilton-Jacobi equations can be easy to implement when written in implicit (or iterated) advection and advection-diffusion equation analogous forms, respectively. These, like the Poisson Laplacian term, are commonly occurring in CFD solvers, allowing the re-use of efficient algorithms and code components. The use of the NASA CFL3D CFD program to solve the implicit Eikonal and Hamilton-Jacobi equations is explored. The re-formulated d equations are easy to implement, and are found to have robust convergence. For accurate Eikonal solutions, upwind metric differences are required. The Poisson approach is also found effective, and easiest to implement. Modified distances are not found to affect global outputs such as lift and drag significantly, at least in common situations such as airfoil flows.

Description: Turbomachines for rocket propulsion applications operate with many different working fluids and flow conditions. Oxidizer boost turbines often operate in liquid oxygen, resulting in an incompressible flow field. Vortex shedding from airfoils in this flow environment can have adverse effects on both turbine performance and durability. In this study the effects of vortex shedding in a low-pressure oxidizer turbine are investigated. Benchmark results are also presented for vortex shedding behind a circular cylinder. The predicted results are compared with available experimental data.

Description: This paper describes the design and testing of a miniature LHP having a 7 mm O.D. evaporator with an integral CC. The vapor line and liquid line are made of 1.6mm stainless steel tubing. The evaporator and the CC are connected on the outer surface by a copper strap and a thermoelectric (TEC) is installed on the strap. The TEC is used to control the CC temperature by applying an electrical current for heating or cooling. Tests performed in ambient included start-up, power cycle, sink temperature cycle, and CC temperature control using TEC. The LHP demonstrated very robust operation in all tests where the heat load varied between 0.5W and 1OOW, and the sink temperature varied between 243K and 293K. The heat leak from the evaporator to the CC was extremely small. The TEC was able to control the CC temperature within +/-0.3K under all test conditions, and the required control heater power was less than 1W.

Description: In our previous computational studies of a generic high-lift configuration, quasi-laminar (as opposed to fully turbulent) treatment of the slat cove region proved to be an effective approach for capturing the unsteady dynamics of the cove flow field. Combined with acoustic propagation via Ffowes Williams and Hawkings formulation, the quasi-laminar simulations captured some important features of the slat cove noise measured with microphone array techniques. However. a direct assessment of the computed cove flow field was not feasible due to the unavailability of off-surface flow measurements. To remedy this shortcoming, we have undertaken a combined experiment and computational study aimed at characterizing the flow structures and fluid mechanical processes within the slat cove region. Part I of this paper outlines the experimental aspects of this investigation focused on the 30P30N high-lift configuration; the present paper describes the accompanying computational results including a comparison between computation and experiment at various angles of attack. Even through predictions of the time-averaged flow field agree well with the measured data, the study indicates the need for further refinement of the zonal turbulence approach in order to capture the full dynamics of the cove's fluctuating flow field.

Description: As a structural entity of turbulence, hairpin vortices are believed to play a major role in developing and sustaining the turbulence process in the near wall region of turbulent boundary layers and may be regarded as the simplest conceptual model that can account for the essential features of the wall pressure fluctuations. In this work we focus on fully developed typical hairpin vortices and estimate the associated surface pressure distributions and their corresponding spectra. On the basis of the attached eddy model, we develop a representation of the overall surface pressure spectra in terms of the eddy size distribution. Instantaneous wavenumber spectra and spatial correlations are readily derivable from this representation. The model is validated by comparison of predicted wavenumber spectra and cross-correlations with existing emperical models and experimental data.

Description: A comprehensive computational and experimental study has been performed at the NASA Langley Research Center as part of the Quiet Aircraft Technology (QAT) Program to investigate the unsteady flow near a leading-edge slat of a two-dimensional, high-lift system. This paper focuses on the experimental effort conducted in the NASA Langley Basic Aerodynamics Research Tunnel (BART) where Particle Image Velocimetry (PIV) data was acquired in the slat cove and at the slat trailing edge of a three-element, high-lift model at 4, 6, and 8 degrees angle of attack and a freestream Mach Number of 0.17. Instantaneous velocities obtained from PIV images are used to obtain mean and fluctuating components of velocity and vorticity. The data show the recirculation in the cove, reattachment of the shear layer on the slat lower surface, and discrete vortical structures within the shear layer emanating from the slat cusp and slat trailing edge. Detailed measurements are used to examine the shear layer formation at the slat cusp, vortex shedding at the slat trailing edge, and convection of vortical structures through the slat gap. Selected results are discussed and compared with unsteady, Reynolds-Averaged Navier-Stokes (URANS) computations for the same configuration in a companion paper by Khorrami, Choudhari, and Jenkins (2004). The experimental dataset provides essential flow-field information for the validation of near-field inputs to noise prediction tools.

Description: The research described in this viewgraph presentation investigates the ascent of STS-107 and foam-debris impact, and contributes to understanding of the STS-107 accident using CFD tools. The goals of the research are to: 1) Quantify loads on foam bipod ramp during ascent; 2) Provide steady-state flow-fields to debris-transport simulations; 3) Simulate flight of foam debris using unsteady six-degree-of-freedom calculations; 4) Provide estimates of foam mass, velocity, and impact angle which correlate with video and film evidence.

Description: Inlets and exhaust nozzles are often omitted or fared over in aerodynamic simulations of aircraft due to the complexities involving in the modeling of engine details such as complex geometry and flow physics. However, the assumption is often improper as inlet or plume flows have a substantial effect on vehicle aerodynamics. A tool for specifying inlet and exhaust plume conditions through the use of high-energy boundary conditions in an established inviscid flow solver is presented. The effects of the plume on the flow fields near the inlet and plume are discussed.

Description: This paper describes the methods used to produce an interactive visualization of a 2 TB computational fluid dynamics (CFD) data set using particle tracing (streaklines). We use the method introduced by Bruckschen et al. [2001] that pre-computes a large number of particles, stores them on disk using a space-filling curve ordering that minimizes seeks, and then retrieves and displays the particles according to the user's command. We describe how the particle computation can be performed using a PC cluster, how the algorithm can be adapted to work with a multi-block curvilinear mesh, and how the out-of-core visualization can be scaled to 296 billion particles while still achieving interactive performance on PG hardware. Compared to the earlier work, our data set size and total number of particles are an order of magnitude larger. We also describe a new compression technique that allows the lossless compression of the particles by 41% and speeds the particle retrieval by about 30%.

Description: Thermal control is a generic need for all spacecraft. In response to ever more demanding science and exploration requirements, spacecraft are becoming ever more complex, and hence their thermal control systems must evolve. This paper briefly discusses the process of technology development, the state-of-the-art in thermal control, recent experiences with on-orbit two-phase systems, and the emerging thermal control technologies to meet these evolving needs. Some "lessons learned" based on experience with on-orbit systems are also presented.

Description: Analysis tools are needed to guide the development and evaluate the performance of multigrid solvers for the fluid flow equations. Classical analysis tools, such as local mode analysis, often fail to accurately predict performance. Two-grid analysis tools, herein referred to as Idealized Coarse Grid and Idealized Relaxation iterations, have been developed and evaluated within a pilot multigrid solver. These new tools are applicable to general systems of equations and/or discretizations and point to problem areas within an existing multigrid solver. Idealized Relaxation and Idealized Coarse Grid are applied in developing textbook-efficient multigrid solvers for incompressible stagnation flow problems.

Description: Low speed flow separation over a wall-mounted hump, and its control using steady suction, were studied experimentally in order to generate a data set for a workshop aimed at validating CFD turbulence models. The baseline and controlled data sets comprised static and dynamic surface pressure measurements, flow field measurements using Particle Image Velocimetry (PIV) and wall shear stress obtained via oil-film interferometry. In addition to the specific test cases studied, surface pressures for a wide variety of conditions were reported for different Reynolds numbers and suction rates. Stereoscopic PIV and oil-film flow visualization indicated that the baseline separated flow field was mainly two-dimensional. With the application of control, some three-dimensionality was evident in the spanwise variation of pressure recovery, reattachment location and spanwise pressure fluctuations. Part 2 of this paper, under preparation for the AIAA Meeting in Reno 2005, considers separation control by means of zero-efflux oscillatory blowing.

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

Description: The synthesis of physical models for gas chemistry and turbulence from the structured grid codes LAURA and VULCAN into the unstructured grid code FUN3D is described. A directionally Symmetric, Total Variation Diminishing (STVD) algorithm and an entropy fix (eigenvalue limiter) keyed to local cell Reynolds number are introduced to improve solution quality for hypersonic aeroheating applications. A simple grid-adaptation procedure is incorporated within the flow solver. Simulations of flow over an ellipsoid (perfect gas, inviscid), Shuttle Orbiter (viscous, chemical nonequilibrium) and comparisons to the structured grid solvers LAURA (cylinder, Shuttle Orbiter) and VULCAN (flat plate) are presented to show current capabilities. The quality of heating in 3D stagnation regions is very sensitive to algorithm options in general, high aspect ratio tetrahedral elements complicate the simulation of high Reynolds number, viscous flow as compared to locally structured meshes aligned with the flow.

Description: The present study examines the impact of large-scale cross flow on the creation of ice roughness elements on the leading edge of a swept wing under glaze icing conditions. A three-dimensional triple-deck structure is developed to describe the local interaction of a 3 D air boundary layer with ice sheets and liquid films. A linear stability analysis is presented here. It is found that, as the sweep angle increases, the local icing instabilities enhance and the most linearly unstable modes are strictly three dimensional.

Description: Analysis tools are needed to guide the development and evaluate the performance of multigrid solvers for the fluid flow equations. Classical analysis tools, such as local mode analysis, often fail to accurately predict performance. Two-grid analysis tools, herein referred to as Idealized Coarse Grid and Idealized Relaxation iterations, have been developed and evaluated within a pilot multigrid solver. These new tools are applicable to general systems of equations and/or discretizations and point to problem areas within an existing multigrid solver. Idealized Relaxation and Idealized Coarse Grid are applied in developing textbook-efficient multigrid solvers for incompressible stagnation flow problems.

Description: Boundary layer ingestion (BLI) is explored as means to improve overall system performance for Blended Wing Body configuration. The benefits of BLI for vehicle system performance benefit are assessed with a process derived from first principles suitable for highly-integrated propulsion systems. This performance evaluation process provides framework within which to assess the benefits of an integrated BLI inlet and lays the groundwork for higher-fidelity systems studies. The results of the system study show that BLI provides a significant improvement in vehicle performance if the inlet distortion can be controlled, thus encouraging the pursuit of active flow control (AFC) as a BLI enabling technology. The effectiveness of active flow control in reducing engine inlet distortion was assessed using a 6% scale model of a 30% BLI offset, diffusing inlet. The experiment was conducted in the NASA Langley Basic Aerodynamics Research Tunnel with a model inlet designed specifically for this type of testing. High mass flow pulsing actuators provided the active flow control. Measurements were made of the onset boundary layer, the duct surface static pressures, and the mass flow through the duct and the actuators. The distortion was determined by 120 total pressure measurements located at the aerodynamic interface plane. The test matrix was limited to a maximum freestream Mach number of 0.15 with scaled mass flows through the inlet for that condition. The data show that the pulsed actuation can reduce distortion from 29% to 4.6% as measured by the circumferential distortion descriptor DC60 using less than 1% of inlet mass flow. Closed loop control of the actuation was also demonstrated using a sidewall surface static pressure as the response sensor.

Description: Experimental test results of air flow inside and at the cylindrical cavity located on axisymmetric body are presented. These tests were conducted in the wind tunnel A-7 of Institute of Mechanics at Moscow State University. Pressure distribution along the cavities and optical measurements were obtained. Dependence of these characteristics of length of a cavity in the range: L/D = 0.5 - 14 and free stream Mach in the range: M(sub infinity) = 0.6 - 3.0 was determined. Flow structure inside the cavity, cause of flow regime change, separation zones geometry and others were studied. In particular, the flow modes of with open and closed separation zones are determined.

Description: This project was composed of three sub-tasks. The objective of the first task was to use the CFD code INS3D to generate both on- and off-design predictions for the consortium optimized impeller flowfield. The results of the flow simulations are given in the first section. The objective of the second task was to construct a turbomachinery testing database comprised of measurements made on several different impellers, an inducer and a diffuser. The data was in the form of static pressure measurements as well as laser velocimeter measurements of velocities and flow angles within the stated components. Several databases with this information were created for these components. The third subtask objective was two-fold: first, to validate the Enigma CFD code for pump diffuser analysis, and secondly, to perform steady and unsteady analyses on some wide flow range diffuser concepts using Enigma. The code was validated using the consortium optimized impeller database and then applied to two different concepts for wide flow diffusers.

Description: The paper contains some experimental and numerical simulation test results on cylindrical blunt body drag reduction using thin spikes or shell mounted in front of a body against a supersonic flow. Experimental tests were conducted using the Aeromechanics and Gas Dynamics Laboratory facilities at the Institute of Mechanics of Moscow State University (IMMSU). Numerical simulations utilizing NASA and IM/MSU codes were conducted at the Hampton University Fluid Mechanics and Acoustics Laboratory. The main purpose of this research is to examine the efficiency of application of multiple spikes for drag reduction and flow stability at the front of a blunt body in different flight conditions, i.e. Mach number, angle of attack, etc. The principal conclusions of these test results are: multiple spike/needle application leads to decrease of drag reduction benefits by comparison with the case of one central mounted needle at the front of a blunt body, but increase lift benefits.

Description: An implementation of both finite rate and equilibrium chemistry have been completed for the OVERFLOW code, a chimera capable, complex geometry flow code widely used to predict transonic flow fields. The implementation builds on the computational efficiency and geometric generality of the solver.

Description: The Microgravity Fluid Physics Program at NASA has developed a substantial investigator base engaging a broad crosssection of the U.S. scientific community. As a result, it enjoys a rich history of many significant scientific achievements. The research supported by the program has produced many important findings that have been published in prestigious journals such as Science, Nature, Journal of Fluid Mechanics, Physics of Fluids, and many others. The focus of the program so far has primarily been on fundamental scientific studies. However, a recent shift in emphasis at NASA to develop advanced technologies to enable future exploration of space has provided motivation to add a strategic research component to the program. This has set into motion a year of intense planning within NASA including three workshops to solicit inputs from the external scientific community. The planning activities and the workshops have resulted in a prioritized list of strategic research issues along with a corresponding detailed roadmap specific to fluid physics. The results of these activities were provided to NASA s Office of Biological and Physical Research (OBPR) to support the development of the Enterprise Strategy document. This paper summarizes these results while showing how the planned research supports NASA s overall vision through OBPR s organizing questions.

Description: Computational fluid dynamics (CFD) has proven to be an invaluable tool for the design and analysis of high- speed propulsion devices. Massively parallel computing, together with the maturation of robust CFD codes, has made it possible to perform simulations of complete engine flowpaths. Steady-state Reynolds-Averaged Navier-Stokes simulations are now routinely used in the scramjet engine development cycle to determine optimal fuel injector arrangements, investigate trends noted during testing, and extract various measures of engine efficiency. Unfortunately, the turbulence and combustion models used in these codes have not changed significantly over the past decade. Hence, the CFD practitioner must often rely heavily on existing measurements (at similar flow conditions) to calibrate model coefficients on a case- by-case basis. This paper provides an overview of the modeled equations typically employed by commercial- quality CFD codes for high-speed combustion applications. Careful attention is given to the approximations employed for each of the unclosed terms in the averaged equation set. The salient features (and shortcomings) of common models used to close these terms are covered in detail, and several academic efforts aimed at addressing these shortcomings are discussed.

Description: A source term model for an array of vortex generators was implemented into a non-proprietary Navier-Stokes computer code, OVERFLOW. The source term models the side force created by a vortex generator vane. The model is obtained by introducing a side force to the momentum and energy equations that can adjust its strength automatically based on the local flow. The model was tested and calibrated by comparing data from numerical simulations and experiments of a single low profile vortex generator vane on a flat plate. In addition, the model was compared to experimental data of an S-duct with 22 co-rotating, low profile vortex generators. The source term model allowed a grid reduction of about seventy percent when compared with the numerical simulations performed on a fully gridded vortex generator on a flat plate without adversely affecting the development and capture of the vortex created. The source term model was able to predict the shape and size of the stream-wise vorticity and velocity contours very well when compared with both numerical simulations and experimental data. The peak vorticity and its location were also predicted very well when compared to numerical simulations and experimental data. The circulation predicted by the source term model matches the prediction of the numerical simulation. The source term model predicted the engine fan face distortion and total pressure recovery of the S-duct with 22 co-rotating vortex generators very well. The source term model allows a researcher to quickly investigate different locations of individual or a row of vortex generators. The researcher is able to conduct a preliminary investigation with minimal grid generation and computational time.

Description: We have succeeded in measuring a substantial portion of the two-point space-time velocity correlation in hot, high speed turbulent jets. This measurement, crucial in aeroacoustic theory and the prediction of jet noise, has been sought for a long time, but has not been made due to the limitations of anemometry. Particle Image Velocimetry has reached a stage of maturity where sufficient measurement density in both time and space allow the computation of space-time correlations. This paper documents these measurements along with lower-order statistics to document the adherence of the jet rig and instrumentation to conventional measures of the turbulence of jets. These measures have been made for a simple round convergent nozzle at acoustic Mach numbers of 0.5, 0.9, both cold and at a static temperature ratio of 2.7, allowing some estimation of the changes in turbulence that take place with changes in jet temperature. Since the dataset described in this paper is very extensive, attention will be focused on validation of the rig and of the measurement systems, and on some of the interesting observations made from studying the statistics, especially as they relate to jet noise. Of note is the effort to study the acoustically relevant part of the space-time correlation by addressing that part of the turbulence kinetic energy that has sonic phase speed.

Description: Experimental and numerical results are presented here for a separate flow nozzle employing chevrons arranged in an alternating pattern on the core nozzle. Comparisons of these results demonstrate that the combination of the WIND/MGBK suite of codes can predict the noise reduction trends measured between separate flow jets with and without chevrons on the core nozzle. Mean flow predictions were validated against Particle Image Velocimetry (PIV), pressure, and temperature data, and noise predictions were validated against acoustic measurements recorded in the NASA Glenn Aeroacoustic Propulsion Lab. Comparisons are also made to results from the CRAFT code. The work presented here is part of an on-going assessment of the WIND/MGBK suite for use in designing the next generation of quiet nozzles for turbofan engines.

Description: For rime ice - where the ice buildup has only rough and jagged surfaces but no protruding horns - this study shows two dimensional CFD analysis based on the one-equation Spalart-Almaras (S-A) turbulence model to predict accurately the lift, drag, and pressure coefficients up to near the stall angle. For glaze ice - where the ice buildup has two or more protruding horns near the airfoil's leading edge - CFD predictions were much less satisfactory because of the large separated region produced by the horns even at zero angle of attack. This CFD study, based on the WIND and the Fluent codes, assesses the following turbulence models by comparing predictions with available experimental data: S-A, standard k-epsilon, shear-stress transport, v(exp 2)-f, and differential Reynolds stress.

Description: Heat pipes are among the most promising technologies for space radiator systems. Water heat pipes are explored in the intermediate temperature range of 400 to above 500 K. The thermodynamic and thermo-physical properties of water are reviewed in this temperature range. Test Data are reported for a copper-water heat pipe. The heat pipe was tested under different orientations. Water heat pipes show promise in this temperature range.Fabrication and testing issues are being addressed.

Description: The newly developed adaptive numerical dissipation control in spatially high order filter schemes for the compressible Euler and Navier-Stokes equations has been recently extended to the ideal and non-ideal magnetohydrodynamics (MHD) equations. These filter schemes are applicable to complex unsteady MHD high-speed shock/shear/turbulence problems. They also provide a natural and efficient way for the minimization of Div(B) numerical error. The adaptive numerical dissipation mechanism consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free from numerical dissipation contamination. The numerical dissipation considered consists of high order linear dissipation for the suppression of high frequency oscillation and the nonlinear dissipative portion of high-resolution shock-capturing methods for discontinuity capturing. The applicable nonlinear dissipative portion of high-resolution shock-capturing methods is very general. The objective of this paper is to investigate the performance of three commonly used types of nonlinear numerical dissipation for both the ideal and non-ideal MHD.

Description: The objectives of this research project were: 1) To study the fluid dynamics of sheared particle-liquid suspensions and the impact of differential particle-fluid inertia; 2) To develop new techniques for observing suspension particle contact and deposition upon solid surfaces. Dr. Yoda was supported by the NASA Office of Biological and Physical Research on a four-year grant from March 2000 through November 2004 for a ground-based study on the fluid dynamics of sheared particle-liquid suspensions and the impact of differential particle-fluid inertia on such flows. Such inertial effects can only be observed in reduced-gravity environments since they are overwhelmed by buoyancy effects on Earth. Moreover, these inertial effects will have a significant impact upon suspension flows in microgravity. Suspension dynamics are of importance in a wide variety of advanced life systems applications, including water reclamation and dust mitigation in confined habitats.

Description: Turbulent separated flow over a two-dimensional hump is computed by solving the RANS equations with k-omega (SST) turbulence model for the baseline, steady suction and oscillatory blowing/suction flow control cases. The flow equations and the turbulent model equations are solved using a fifth-order accurate weighted essentially nonoscillatory (WENO) scheme for space discretization and a third order, total variation diminishing (TVD) Runge-Kutta scheme for time integration. Qualitatively the computed pressure distributions exhibit the same behavior as they are observed in the experiments. The computed separation regions are much longer than that are observed. However, the percentage reduction in the separation region in the steady suction case is closer to that was measured in the experiment. The computations did not predict the expected reduction in the separation length in the oscillatory case. The predicted turbulent quantities are two to three times smaller than that are measured and it points towards the deficiencies in the existing turbulent models when they are applied to strong steady/unsteady separated flows.

Description: A new computer program to design nozzles with non-uniform inflow has been developed using the rotational method of characteristics (MOC). This program has been used to design a nozzle for the NASA's HYPULSE shock-expansion tunnel for use in scramjet engine tests at a Mach-15 flight-enthalpy condition. The nozzle has an area ratio of 9.5:1 that expands the inflow from Mach 6 along the centerline to Mach 8.7. Although the density and Mach number vary radially at the exit due to the non-uniformities of the inflow, the MOC procedure produces exit flow that is parallel and has uniform static pressure. The design has been verified with CFD which compares favorably with the MOC solution.

Description: We used a permanent-magnet MRI system to estimate the integral and spatially- and/or temporally-resolved void-fraction distributions and flow patterns in gas-liquid two-phase flows. Air was introduced at the bottom of the stagnant liquid column using an accurate and programmable syringe pump. Air flow rates were varied between 1 and 200 ml/min. The cylindrical non-conducting test tube in which two-phase flow was measured was placed in a 2.67 kGauss MRI with MRT spectrometer/imager. Roughly linear relationship has been obtained for the integral void-fraction, obtained by volume-averaging of the spatially-resolved signals, and the air flow rate in upward direction. The time-averaged spatially-resolved void fraction has also been obtained for the quasi-steady flow of air in a stagnant liquid column. No great accuracy is claimed as this was an exploratory proof-of-concept type of experiment. Preliminary results show that MRI a non-invasive and non-intrusive experimental technique can indeed provide a wealth of different qualitative and quantitative data and is especially well suited for averaged transport processes in adiabatic and diabatic multi-phase and/or multi-component flows.

Description: Active flow control of boundary-layer separation using glow-discharge plasma actuators is studied experimentally. Separation is induced on a flat plate installed in a closed-circuit wind tunnel by a shaped insert on the opposite wall. The flow conditions represent flow over the suction surface of a modern low-pressure-turbine airfoil. The Reynolds number, based on wetted plate length and nominal exit velocity, is varied from 50,000 to 300,000, covering cruise to takeoff conditions. Low (0.2 percent) and high (2.5 percent) free-stream turbulence intensities are set using passive grids. A spanwise-oriented phased-plasma-array actuator, fabricated on a printed circuit board, is surface-flush-mounted upstream of the separation point and can provide forcing in a wide frequency range. Static surface pressure measurements and hot-wire anemometry of the base and controlled flows are performed and indicate that the glow-discharge plasma actuator is an effective device for separation control.

Description: The object of this paper is to investigate the feasibility of applying CFD methods to aerodynamic analyses for aircraft stability and control. The integrated aerodynamic parameters used in stability and control, however, are not necessarily those extensively validated in the state of the art CFD technology. Hence, an exploratory study of such applications and the comparison of the solutions to available experimental data will help to assess the validity of the current computation methods. In addition, this study will also examine issues related to wind tunnel measurements such as measurement uncertainty and support interference effects. Several sets of experimental data from the NASA Langley 14x22-Foot Subsonic Tunnel and the National Transonic Facility are presented. Two Navier-Stokes flow solvers, one using structured meshes and the other unstructured meshes, were used to compute longitudinal static stability derivatives for an advanced Blended Wing Body configuration over a wide range of angles of attack. The computations were performed for two different Reynolds numbers and the resulting forces and moments are compared with the above mentioned wind tunnel data.

Description: A rapid, structured volume grid smoothing and adaption technique, based on signal processing methods, was developed and applied to the Shuttle Orbiter at hypervelocity flight conditions in support of the Columbia Accident Investigation. Because of the fast pace of the investigation, computational aerothermodynamicists, applying hypersonic viscous flow solving computational fluid dynamic (CFD) codes, refined and enhanced a grid for an undamaged baseline vehicle to assess a variety of damage scenarios. Of the many methods available to modify a structured grid, most are time-consuming and require significant user interaction. By casting the grid data into different coordinate systems, specifically two computational coordinates with arclength as the third coordinate, signal processing methods are used for filtering the data [Taubin, CG v/29 1995]. Using a reverse transformation, the processed data are used to smooth the Cartesian coordinates of the structured grids. By coupling the signal processing method with existing grid operations within the Volume Grid Manipulator tool, problems related to grid smoothing are solved efficiently and with minimal user interaction. Examples of these smoothing operations are illustrated for reduction in grid stretching and volume grid adaptation. In each of these examples, other techniques existed at the time of the Columbia accident, but the incorporation of signal processing techniques reduced the time to perform the corrections by nearly 60%. This reduction in time to perform the corrections therefore enabled the assessment of approximately twice the number of damage scenarios than previously possible during the allocated investigation time.