ALBERT

Description: 〈br /〉〈span class="detail_caption"〉Notes: 〈/span〉Ajayi, Thomas -- Sherman, Kenneth -- Tang, Qisheng -- New York, N.Y. -- Science. 2002 Aug 2;297(5582):772.〈br /〉〈span class="detail_caption"〉Record origin:〈/span〉 〈a href="http://www.ncbi.nlm.nih.gov/pubmed/12162321" target="_blank"〉PubMed〈/a〉

Description: There is debate concerning the most effective conservation of marine biodiversity, especially regarding the appropriate location, size, and connectivity of marine reserves. We describe a means of establishing marine reserve networks by using optimization algorithms and multiple levels of information on biodiversity, ecological processes (spawning, recruitment, and larval connectivity), and socioeconomic factors in the Gulf of California. A network covering 40% of rocky reef habitat can fulfill many conservation goals while reducing social conflict. This quantitative approach provides a powerful tool for decision-makers tasked with siting marine reserves.〈br /〉〈span class="detail_caption"〉Notes: 〈/span〉Sala, Enric -- Aburto-Oropeza, Octavio -- Paredes, Gustavo -- Parra, Ivan -- Barrera, Juan C -- Dayton, Paul K -- New York, N.Y. -- Science. 2002 Dec 6;298(5600):1991-3.〈br /〉〈span class="detail_caption"〉Author address: 〈/span〉Center for Marine Biodiversity and Conservation, Scripps Institution of Oceanography, La Jolla, CA 92093, USA. esala@ucsd.edu〈br /〉〈span class="detail_caption"〉Record origin:〈/span〉 〈a href="http://www.ncbi.nlm.nih.gov/pubmed/12471258" target="_blank"〉PubMed〈/a〉

Description: Various alternative formulations of the LES equations have been explored in which additional evolution equations for variables such as the acceleration, the subgrid-scale stress tensor, or the subgrid-scale force are explicitly carried. Statistics of the velocity field obtained from the equation for the acceleration are shown to depend strongly on the initial conditions. This feature, which is independent of LES modeling issues, seems to prove that the velocity-acceleration formulation of the Navier-Stokes is not useful for numerical simulation. Equations for the subgrid-scale quantities appear to be much more stable. However, models required by this formulation of the LES problem still require additional study.

Description: Two approaches for the identification of internal gravity waves in sheared and unsheared homogeneous stratified turbulence are investigated. First, the phase angle between the vertical velocity and density fluctuations is considered. It is found, however, that a continuous distribution of the phase angle is present in weakly and strongly stratified flow. Second, a projection onto the solution of the linearized inviscid equations of motion of unsheared stratified flow is investigated. It is found that a solution of the fully nonlinear viscous Navier-Stokes equations can be represented by the linearized inviscid solution. The projection yields a decomposition into vertical wave modes and horizontal vortical modes.

Description: The bubble coalescence and interfacial instabilities that are important to modeling critical heat flux (CHF) in reduced-gravity systems can be sensitive to even minute body forces. Understanding these complex phenomena is vital to the design and safe implementation of two-phase thermal management loops proposed for space and planetary-based thermal systems. While reduced gravity conditions cannot be accurately simulated in 1g ground-based experiments, such experiments can help isolate the effects of the various forces (body force, surface tension force and inertia) which influence flow boiling CHF. In this project, the effects of the component of body force perpendicular to a heated wall were examined by conducting 1g flow boiling experiments at different orientations. FC-72 liquid was boiled along one wall of a transparent rectangular flow channel that permitted photographic study of the vapor-liquid interface at conditions approaching CHF. High-speed video imaging was employed to capture dominant CHF mechanisms. Six different CHF regimes were identified: Wavy Vapor Layer, Pool Boiling, Stratification, Vapor Counterflow, Vapor Stagnation, and Separated Concurrent Vapor Flow. CHF showed great sensitivity to orientation for flow velocities below 0.2 m/s, where very small CHF values where measured, especially with downflow and downward-facing heated wall orientations. High flow velocities dampened the effects of orientation considerably. Figure I shows representative images for the different CHF regimes. The Wavy Vapor Layer regime was dominant for all high velocities and most orientations, while all other regimes were encountered at low velocities, in the downflow and/or downward-facing heated wall orientations. The Interfacial Lift-off model was modified to predict the effects of orientation on CHF for the dominant Wavy Vapor Layer regime. The photographic study captured a fairly continuous wavy vapor layer travelling along the heated wall while permitting liquid contact only in wetting fronts, located in the troughs of the interfacial waves. CHF commenced when wetting fronts near the outlet were lifted off the wall. The Interfacial Lift-off model is shown to be an effective tool for predicting the effects of body force on CHF at high velocities.

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

Description: Parallel Plate Plastometer (PPP) is a device commonly used for measuring the viscosity of high polymers at low rates of shear in the range 10(exp 4) to 10(exp 9) poises. This device is being validated for use in measuring the viscosity of liquid glasses at high temperatures having similar ranges for the viscosity values. PPP instrument consists of two similar parallel plates, both in the range of 1 inch in diameter with the upper plate being movable while the lower one is kept stationary. Load is applied to the upper plate by means of a beam connected to shaft attached to the upper plate. The viscosity of the fluid is deduced from measuring the variation of the plate separation, h, as a function of time when a specified fixed load is applied on the beam. Operating plate speeds measured with the PPP is usually in the range of 10.3 cm/s or lower. The flow field within the PPP can be simulated using the equations of motion of fluid flow for this configuration. With flow speeds in the range quoted above the flow field between the two plates is certainly incompressible and laminar. Such flows can be easily simulated using numerical modeling with computational fluid dynamics (CFD) codes. We present below the mathematical model used to simulate this flow field and also the solutions obtained for the flow using a commercially available finite element CFD code.

Description: The Physics of Colloids in Space (PCS) experiment was accommodated within International Space Station (ISS) EXpedite the PRocessing of Experiments to Space Station (EXPRESS) Rack 2 and was remotely operated from early June 2001 until February 2002 from NASA Glenn Research Center's Telescience Support Center in Cleveland, Ohio and from a remote site at Harvard University in Cambridge, Massachusetts. PCS is an experiment conceived by Professor David A. Weitz of Harvard University (the Principal Investigator), focusing on the behavior of three different classes of colloid mixtures. The sophisticated light scattering instrumentation comprising PCS is capable of color imaging, and dynamic and static light scattering from 11 to 169 degrees, Bragg scattering over the range from 10 to 60 degrees, and laser light scattering at low angles from 0.3 to 6.0 degrees. The PCS instrumentation performed remarkably well, demonstrating a flexibility that enabled experiments to be performed that had not been envisioned prior to launch. While on-orbit, PCS accomplished 2400 hours of science operations, and was declared a resounding success. Each of the eight sample cells worked well and produced interesting and important results. Crystal nucleation and growth and the resulting structures of two binary colloidal crystal alloys were studied, with the long duration microgravity environment of the ISS facilitating extended studies on the growth and coarsening characteristics of the crystals. In another experiment run, the de-mixing of the colloid-polymer critical-point sample was studied as it phase-separates into two phases, one that resembles a gas and one that resembles a liquid. This process was studied over four decades of length scale, from 1 micron to 1 centimeter, behavior that cannot be observed in this sample on Earth because sedimentation would cause the colloids to fall to the bottom of the cell faster than the de-mixing process could occur. Similarly, the study of gelation and aging of another colloid-polymer sample, the colloid-polymer gel, also provided valuable information on gelation mechanisms, as did investigations on the extremely the low concentration silica and polystyrene fractal gel samples.

Description: The overall objective of this work is to study nucleate boiling heat transfer under microgravity conditions in such a way that while providing basic knowledge of the phenomena, it also leads to development of simulation models and correlations that can be used as design tools for a wide range of gravity levels. In the study a building block type of approach is used and both pool and low velocity flow boiling are investigated. Starting with experiments using a single bubble, the complexity of the experiments is increased to two or three inline bubbles, to five bubbles placed on a two-dimensional grid. Finally, experiments are conducted where a large number of prescribed cavities nucleate on the heater and when a commercial surface is used. So far experiments have been conducted at earth normal gravity and in the reduced gravity environment of the KC-135 aircraft whereas experiments on the space station are planned. Modeling/complete numerical simulation of the boiling process is an integral part of the total effort. Experiments conducted with single bubbles formed on a nucleation site microfabricated on a polished silicon wafer show that for gravity levels (g) varying from 1.5g(sub e) to 0.01g(sub e), the bubble diameter at departure varies approximately as (g(sub e)/g)(exp 1/2) and the growth period as (g(sub e)/g). When bubbles merge either inline or in a plane, the bubble diameter at departure is found to be smaller than that obtained for a single bubble and shows a weaker dependence on the level of gravity. The possible reason is that as the bubbles merge they create fluid circulation around the bubbles, which in turn induces a lift force that is responsible for the earlier departure of the bubbles. The verification of this proposition is being sought through numerical simulations. There is a merger of two inline, three inline, and several bubbles in a plane in the low gravity environment of the KC-135 aircraft. After merger and before departure, a mushroom type of bubble with several stems attached to the heater surface is clearly evident. Local heat fluxes during growth and departure of a single bubble were also measured. It was found that during most of the growth period of the bubble, generally the wall heat flux decreased with time because of the increased dry area under the bubble. However, the heat flux increased rapidly just prior to departure of the bubble because of the transient conduction into the cold liquid rushing to fill the space vacated by the bubble as the bubble base shrinks. The measured heat fluxes at various radial locations are found to be in qualitative agreement with the numerical predictions. Single bubble studies at earth normal gravity have also been performed on surfaces oriented at different angles to the gravitational acceleration with flow parallel to the surface. It is found that in all cases the bubbles slide along the surface before lift-off from the surface. The lift force generated as a result of the relative motion between the sliding bubbles and the imposed flow is found to play an important role when the normal force due to buoyancy is reduced. An experimental apparatus for the study of the bubble behavior with imposed flow under reduced gravity conditions has been developed and will soon be employed for experiments in the KC-135 aircraft.

Description: Dynamic wetting, the displacement of one fluid by another immiscible fluid on a surface, controls many natural and technological phenomena, such as coating, printing, spray painting and lubricating. Particularly in coating and spraying applications, contact lines advance across pre-existing fluid films. Most previous work has focused on contact lines advancing across films sufficiently thick that they behave as simple Newtonian fluids. Ultrathin films, where the film thickness may impinge on fundamental length scales in the fluid, have received less attention. In this talk, we will discuss the effects of ultrathin polymer films on dynamic wetting. We measure the interface shape within microns of moving contact lines advancing across preexisting films and compare the measurements to existing models of viscous bending for interfaces advancing across dry surfaces and 'thick' (in the sense that they behave as liquids) films. In the experiments, we advance a contact line of 10-poise and 1-poise polydimethylsiloxane (silicone oil) across pre-coated films of the same fluid with thickness from a single chain thickness (approx. 10 A) through a couple of radii of gyration (100-200 A) to films so thick they are likely bulk in behavior (10(exp 3) A). All films are physisorbed, i.e. they readily rinse from the surface. Thus, molecules in the film are not anchored to the surface and can move within the film if the hydrodynamics dictate such motion. For films of the thickness of a single chain (approx. 10 A), our experiments indicate that the advancing fluid behaves just as it would if it advanced over a dry surface. For the thicker films (10(exp 3) A), we find behavior indicating that the molecules in the film are acting as a fluid with the bulk properties. In this regime, results for the two different fluids are identical when the experiments are performed at the same pre-existing film thickness and advancing capillary number, Ca. For film of thickness of a few radii of gyration (approx. 100-200 A), the behavior depends on Ca of the advancing meniscus. At low Ca, the viscous bending of the interface near the contact line does not behave as it would on a dry surface. It has a lower curvature than expected. However, at higher Ca, the viscous bending is described by the model for spreading over a dry surface. These results show that the fluid flow in the film does behave differently than bulk as the film thickness becomes comparable to molecular length scale. But even more intriguing is the unusual velocity dependence of that behavior where the film behaves more solid-like at higher contact line speeds. We will discuss these results in terms of the properties of confined polymer melts.

Description: An important component in biotechnology, particularly in the area of protein engineering and rational drug design is the knowledge of the precise three-dimensional molecular structure of proteins. The quality of structural information obtained from X-ray diffraction methods is directly dependent on the degree of perfection of the protein crystals. As a consequence, the growth of high quality macromolecular crystals for diffraction analyses has been the central focus for biochemists, biologists, and bioengineers. Macromolecular crystals are obtained from solutions that contain the crystallizing species in equilibrium with higher aggregates, ions, precipitants, other possible phases of the protein, foreign particles, the walls of the container, and a likely host of other impurities. By changing transport modes in general, i.e., reduction of convection and sedimentation, as is achieved in "microgravity", researchers have been able to dramatically affect the movement and distribution of macromolecules in the fluid, and thus their transport, formation of crystal nuclei, and adsorption to the crystal surface. While a limited number of high quality crystals from space flights have been obtained, as the recent National Research Council (NRC) review of the NASA microgravity crystallization program pointed out, the scientific approach and research in crystallization of proteins has been mainly empirical yielding inconclusive results. We postulate that we can reduce convection in ground-based experiments and we can understand the different aspects of convection control through the use of strong magnetic fields and field gradients. Whether this limited convection in a magnetic field will provide the environment for the growth of high quality crystals is still a matter of conjecture that our research will address. The approach exploits the variation of fluid magnetic susceptibility with concentration for this purpose and the convective damping is realized by appropriately positioning the crystal growth cell so that the magnetic susceptibility force counteracts terrestrial gravity. The general objective is to test the hypothesis of convective control using a strong magnetic field and magnetic field gradient and to understand the nature of the various forces that come into play. Specifically we aim to delineate causative factors and to quantify them through experiments, analysis and numerical modeling. Once the basic understanding is obtained, the study will focus on testing the hypothesis on proteins of pyruvate dehydrogenase complex (PDC), proteins E1 and E3. Obtaining high crystal quality of these proteins is of great importance to structural biologists since their structures need to be determined. Specific goals for the investigation are: 1. To develop an understanding of convection control in diamagnetic fluids with concentration gradients through experimentation and numerical modeling. Specifically solutal buoyancy driven convection due to crystal growth will be considered. 2. To develop predictive measures for successful crystallization in a magnetic field using analyses and numerical modeling for use in future protein crystal growth experiments. This will establish criteria that can be used to estimate the efficacy of magnetic field flow damping on crystallization of candidate proteins. 3. To demonstrate the understanding of convection damping by high magnetic fields to a class of proteins that is of interest and whose structure is as yet not determined. 4. To compare quantitatively, the quality of the grown crystals with and without a magnetic field. X-ray diffraction techniques will be used for the comparative studies. In a preliminary set of experiments, we studied crystal dissolution effects in a 5 Tesla magnet available at NASA Marshall Space Flight Center (MSFC). Using a Schlieren setup, a 1mm crystal of Alum (Aluminum-Potassium Sulfate) was introduced in a 75% saturated solution and the resulting dissolution plume was observed. The experiment was conducted both in the presence and absence of a magnetic field gradient. The magnet produces a gradient field of approx. 1 Tesla2/cm. Image analysis of the recorded images indicated an enhanced plume velocity that was of the order of the measurement limit. For this experiment, both the gradient and gravity fields are in the same direction resulting in an enhanced effective gravity that tends to accelerate the observed plume velocity. While the results are not conclusive, pending further tests, it clearly points out the inadequacy of the MSFC magnet for conducting protein crystallization experiments and the need for a stronger magnet. In spacebased experiments, however, where the gravitational effects are small, only a weak magnetic field will be required to control or mitigate the effects of convective contamination.

Description: The wall-matching methodology of Wilcox is modified to include a solid-wall, thermal-conduction model. This coupled fluid-thermal-structure model is derived assuming that the wall thermal-structure behavior is locally one-dimensional and that structural deformations, due to thermally induced stresses, are not significant. The one-dimensional coupled fluid-thermal-structure model is derived such that the wall temperature is removed as an independent boundary condition variable. The one-dimensional coupled fluid-thermal-structure model is also derived for the general case of an arbitrary mixture of thermally prefect gases and a wall of arbitrary thickness and conductivity by using a compressible, streamwise-pressure-gradient-corrected, wall-matching function and Fourier's law of heat conduction. The resulting model was implemented in the VULCAN CFD code as a new boundary condition type. VULCAN was then used to simulate a two-dimensional Mach 6 wind tunnel facility nozzle flow to demonstrate/validate the one-dimensional coupled fluid-thermal-structure model. The nozzle internal-wall surface temperature and heat transfer distributions computed using the one-dimensional coupled fluid-thermal-structure model are compared to wall temperature and heat transfer distributions from an iterative multi-dimensional analysis obtained by coupling the VULCAN CFD code and the MSC/NASTRAN-thermal code. The one-dimensional coupled fluid-thermal-structure model analysis is shown to be very robust and in excellent agreement with the multi-dimensional iteratively coupled analysis. It is also shown that the one-dimensional analysis can be used as an initial guess for the multi-dimensional iteratively coupled analysis.

Description: Results from Direct Numerical Simulations of temporal, supercritical mixing layers for two species systems are analyzed to elucidate species-specific turbulence aspects. The two species systems, O2/H2 and C7HG16/N2, have different thermodynamic characteristics; thus, although the simulations are performed at similar reduced pressure (ratio of the pressure to the critical pressure), the former system is dose to mixture ideality and has a relatively high solubility with respect to the latter, which exhibits strong departures from mixture ideality Due to the specified, smaller initial density stratification, the C7H16/N2 layers display higher growth and increased global molecular mixing as well as larger turbulence levels. However, smaller density gradients at the transitional state for the O2/H2 system indicate that on a local basis, the layer exhibits an enhanced mixing, this being attributed to the increased solubility and to mixture ideality. These thermodynamic features are shown to affect the irreversible entropy production (i.e. the dissipation), which is larger for the O2/H2 layer and is primarily concentrated in high density-gradient magnitude regions that are distortions of the initial density stratification boundary. In contrast, the regions of largest dissipation in the C7H16/N2 layer are located in high density-gradient magnitude regions resulting from the mixing of the two fluids.

Description: High temperature composite heat exchangers are an enabling technology for a number of aeropropulsion applications. They offer the potential for mass reductions of greater than fifty percent over traditional metallics designs and enable vehicle and engine designs. Since they offer the ability to operate at significantly higher operating temperatures, they facilitate operation at reduced coolant flows and make possible temporary uncooled operation in temperature regimes, such as experienced during vehicle reentry, where traditional heat exchangers require coolant flow. This reduction in coolant requirements can translate into enhanced range or system payload. A brief review of the approaches and challengers to exploiting this important technology are presented, along with a status of recent government-funded projects.

Description: Inspired by the rapid elongation of air columns injected into vortices by dolphins, we present an exact inviscid solution for the axial speed (assumed steady) of propagation of the tip of a semi-infinite cylindrical bubble along the axis of a cylindrical vortex. The bubble is assumed to be held at constant pressure by being connected to a reservoir, the lungs of the dolphin, say. For a given bubble pressure, there is a modest critical rotation rate above which steadily propagating bubbles exist. For a bubble at ambient pressure, the propagation speed of the bubble (relative to axial velocity within the vortex) varies between 0.5 and 0.6 of the maximum rotational speed of the vortex. Surprisingly, the bubble tip can propagate (almost as rapidly) even when the pressure minimum in the vortex core is greater than the bubble pressure; in this case, solutions exhibit a dimple on the nose of the bubble. A situation important for incipient vortex cavitation, and one which dolphins also demonstrate, is elongation of a free bubble, i.e., one whose internal pressure may vary. Under the assumption that the acceleration term is small (checked a posteriori), the steady solution is applied at each instant during the elongation. Three types of behavior are then possible depending on physical parameters and initial conditions: (A) Unabated elongation with slowly increasing bubble pressure, and nearly constant volume. Volume begins to decrease in the late stages. (B1) Elongation with decreasing bubble pressure. A limit point of the steady solution is encountered at a finite bubble length. (B2) Unabated elongation with decreasing bubble pressure and indefinite creation of volume. This is made possible by the existence of propagating solutions at bubble pressures below the minimum vortex pressure. As the bubble stretches, its radius initially decreases but then becomes constant; this is also observed in experiments on incipient vortex cavitation.

Description: This paper reviews recent progress made in incompressible Navier-Stokes simulation procedures and their application to problems of engineering interest. Discussions are focused on the methods designed for complex geometry applications in three dimensions, and thus are limited to primitive variable formulation. A summary of efforts in flow solver development is given followed by numerical studies of a few example problems of current interest. Both steady and unsteady solution algorithms and their salient features are discussed. Solvers discussed here are based on a structured-grid approach using either a finite -difference or a finite-volume frame work. As a grand-challenge application of these solvers, an unsteady turbopump flow simulation procedure has been developed which utilizes high performance computing platforms. In the paper, the progress toward the complete simulation capability of the turbo-pump for a liquid rocket engine is reported. The Space Shuttle Main Engine (SSME) turbo-pump is used as a test case for evaluation of two parallel computing algorithms that have been implemented in the INS3D code. The relative motion of the grid systems for the rotorstator interaction was obtained using overact grid techniques. Unsteady computations for the SSME turbo-pump, which contains 114 zones with 34.5 million grid points, are carried out on SCSI Origin 3000 systems at NASA Ames Research Center. The same procedure has been extended to the development of NASA-DeBakey Ventricular Assist Device (VAD) that is based on an axial blood pump. Computational, and clinical analysis of this device are presented.

Description: We study the no-slip boundary conditions for water at a hydrophobic (graphite) surface using non-equilibrium molecular-dynamics simulations. For the planar Couette flow, we find a slip length of 64 nm at 1 bar and 300 K, decreasing with increasing system pressure to a value of 31 nm at 1000 bar. Changing the properties of the interface to from hydrophobic to strongly hydrophilic reduces the slip to 14 nm. Finally, we study the flow of water past an array of carbon nanotubes mounted in an inline configuration with a spacing of 16.4 x 16.4 nm. For tube diameters of 1.25 and 2.50 nm we find drag coefficients in good agreement with the macroscopic, Navier-Stokes values. For carbon nanotubes, the no-slip condition is valid to within the definition of the position of the interface.

Description: In this work we analyze the issue of the formation of condensation trails ("contrails") in the near-field of an aircraft wake. The basic configuration consists in an exhaust engine jet interacting with a wing-tip training vortex. The procedure adopted relies on a mixed Eulerian/Lagrangian two-phase flow approach; a simple micro-physics model for ice growth has been used to couple ice and vapor phases. Large eddy simulations have carried out at a realistic flight Reynolds number to evaluate the effects of turbulent mixing and wake vortex dynamics on ice-growth characteristics and vapor thermodynamic properties.

Description: Mixing driven by buoyancy-induced flows is of particular interest to microgravity processes, as the body force that governs the intensity of flow fields can be directly controlled. We consider a model experimental system to explore the dynamics of mixing which employs two miscible liquids inside a cavity separated initially by a divider. The two liquids are oriented vertically inside a rectangular cavity with constant width and height, and varying depths to span the range of a Hele-Shaw cell to a 3-D configuration. The two miscible liquids can be sufficiently diluted and died, for example water and deuterium oxide, such that a distinct interface exists across the divider. The transient mixing characteristic of the two fluids is addressed by following the Lagrangian history of the interface for various aspect ratios in the z-plane (depth variation) as well as a range of pulling velocities of the divider. The mixing characteristics of the two fluids are quantified from measurement of the length stretch of the interface and its flow field using respectively image processing techniques and Particle Imaging Velocimetry. Scaling analysis shows that the length stretch depends on four governing parameters, namely the Grashof number (Gr), Schmidt number (Sc), aspect ratio (Ar), and Reynolds number (Re). Variation of the Schmidt number is taken into account through thermophysical property variation. Thus our problem reduces to a codimension three bifurcation in parametric space for Gr, Ar, and Re.

Description: Richtmyer-Meshkov (RM) instability occurs when two different density fluids are impulsively accelerated in the direction normal to their nearly planar interface. The instability causes small perturbations on the interface to grow and eventually become a turbulent flow. It is closely related to Rayleigh-Taylor instability, which is the instability of a planar interface undergoing constant acceleration, such as caused by the suspension of a heavy fluid over a lighter one in the earth's gravitational field. Like the well-known Kelvin-Helmholtz instability, RM instability is a fundamental hydrodynamic instability which exhibits many of the nonlinear complexities that transform simple initial conditions into a complex turbulent flow. Furthermore, the simplicity of RM instability (in that it requires very few defining parameters), and the fact that it can be generated in a closed container, makes it an excellent test bed to study nonlinear stability theory as well as turbulent transport in a heterogeneous system. However, the fact that RM instability involves fluids of unequal densities which experience negligible gravitational force, except during the impulsive acceleration, requires RM instability experiments to be carried out under conditions of microgravity. This experimental study investigates the instability of an interface between incompressible, miscible liquids with an initial sinusoidal perturbation. The impulsive acceleration is generated by bouncing a rectangular tank containing two different density liquids off a retractable vertical spring. The initial perturbation is produced prior to release by oscillating the tank in the horizontal direction to produce a standing wave. The instability evolves in microgravity as the tank travels up and then down the vertical rails of a drop tower until hitting a shock absorber at the bottom. Planar Laser Induced Fluorescence (PLIF) is employed to visualize the flow. PLIF images are captured by a video camera that travels with the tank. Figure 1 is as sequence of images showing the development of the instability from the initial sinusoidal disturbance far into the nonlinear regime which is characterized by the appearance of mushroom structures resulting from the coalescence of baroclinic vorticity produced by the impulsive acceleration. At later times in this sequence the vortex cores are observed to become unstable showing the beginnings of the transition to turbulence in this flow. The amplitude of the growing disturbance after the impulsive acceleration is measured and found to agree well with theoretical predictions. The effects of Reynolds number (based on circulation) on the development of the vortices and the transition to turbulence are also determined.

Description: Three flow regimes have been identified for gas-liquid flow in a microgravity environment: Bubble, Slug, and Annular. For the slug and annular flow regimes, the behavior observed in vertical upflow in normal gravity is similar to microgravity flow with a thin, symmetrical annular film wetting the tube wall. However, the motion and behavior of this film is significantly different between the normal and low gravity cases. Specifically, the liquid film will slow and come to a stop during low frequency wave motion or slugging. In normal gravity vertical upflow, the film has been observed to slow, stop, and actually reverse direction until it meets the next slug or wave.

Description: One of the longest standing unsolved problems in physics relates to the behavior of fluids that are driven far from equilibrium such as occurs when they become turbulent due to fast flow through a grid or tidal motions. In turbulent flows the distribution of vortex energy as a function of the inverse length scale [or wavenumber 'k'] of motion is proportional to 1/k(sup 5/3) which is the celebrated law of Kolmogorov. Although this law gives a good description of the average motion, fluctuations around the average are huge. This stands in contrast with thermally activated motion where large fluctuations around thermal equilibrium are highly unfavorable. The problem of turbulence is the problem of understanding why large fluctuations are so prevalent which is also called the problem of 'intermittency'. Turbulence is a remarkable problem in that its solution sits simultaneously at the forefront of physics, mathematics, engineering and computer science. A recent conference [March 2002] on 'Statistical Hydrodynamics' organized by the Los Alamos Laboratory Center for Nonlinear Studies brought together researchers in all of these fields. Although turbulence is generally thought to be described by the Navier-Stokes Equations of fluid mechanics the solution as well as its existence has eluded researchers for over 100 years. In fact proof of the existence of such a solution qualifies for a 1 M$ millennium prize. As part of our NASA funded research we have proposed building a bridge between vortex turbulence and wave turbulence. The latter occurs when high amplitude waves of various wavelengths are allowed to mutually interact in a fluid. In particular we have proposed measuring the interaction of ripples [capillary waves] that run around on the surface of a fluid sphere suspended in a microgravity environment. The problem of ripple turbulence poses similar mathematical challenges to the problem of vortex turbulence. The waves can have a high amplitude and a strong nonlinear interaction. Furthermore, the steady state distribution of energy again follows a Kolmogorov scaling law; in this case the ripple energy is distributed according to 1/k (sup 7/4). Again, in parallel with vortex turbulence ripple turbulence exhibits intermittency. The problem of ripple turbulence presents an experimental opportunity to generate data in a controlled, benchmarked system. In particular the surface of a sphere is an ideal environment to study ripple turbulence. Waves run around the sphere and interact with each other, and the effect of walls is eliminated. In microgravity this state can be realized for over 2 decades of frequency. Wave turbulence is a physically relevant problem in its own right. It has been studied on the surface of liquid hydrogen and its application to Alfven waves in space is a source of debate. Of course, application of wave turbulence perspectives to ocean waves has been a major success. The experiment which we plan to run in microgravity is conceptually straightforward. Ripples are excited on the surface of a spherical drop of fluid and then their amplitude is recorded with appropriate photography. A key challenge is posed by the need to stably position a 10cm diameter sphere of water in microgravity. Two methods are being developed. Orbitec is using controlled puffs of air from at least 6 independent directions to provided the positioning force. This approach has actually succeeded to position and stabilize a 4cm sphere during a KC 135 segment. Guigne International is using the radiation pressure of high frequency sound. These transducers have been organized into a device in the shape of a dodecahedron. This apparatus 'SPACE DRUMS' has already been approved for use for combustion synthesis experiments on the International Space Station. A key opportunity presented by the ripple turbulence data is its use in driving the development of codes to simulate its properties.

Description: An in-depth review of boundary-layer flow-separation control by a passive method using low-profile vortex generators is presented. The generators are defined as those with a device height between 10% and 50% of the boundary layer thickness. Key results are presented for several research efforts, all of which were performed within the past decade and a half where the majority of these works emphasize experimentation with some recent efforts on numerical simulations. Topics of discussion consist of both basic fluid dynamics and applied aerodynamics research. The fluid dynamics research includes comparative studies on separation control effectiveness as well as device-induced vortex characterization and correlation. The comparative studies cover the controlling of low-speed separated flows in adverse pressure gradient and supersonic shock-induced separation. The aerodynamics research includes several applications for aircraft performance enhancement and covers a wide range of speeds. Significant performance improvements are achieved through increased lift and/or reduced drag for various airfoils-low-Reynolds number, high-lift, and transonic-as well as highly swept wings. Performance enhancements for non-airfoil applications include aircraft interior noise reduction, inlet flow distortion alleviation inside compact ducts, and a more efficient overwing fairing. The low-profile vortex generators are best for being applied to applications where flow-separation locations are relatively fixed and the generators can be placed reasonably close upstream of the separation. Using the approach of minimal near-wall proturbances through substantially reduced device height, these devices can produce streamwise vortices just strong enough to overcome the separation without unnecessarily persisting within the boundary layer once the flow-control objective is achieved. Practical advantages of low-profile vortex generators, such as their inherent simplicity and low device drag, are demonstrated to be critically important for many applications as well.

Description: This is an old issue in computational fluid dynamics (CFD). What is the so-called contravariant velocity or contravariant velocity component? In the article, we review the basics of tensor analysis and give the contravariant velocity component a rigorous explanation. For a given coordinate system, there exist two uniquely determined sets of base vector systems - one is the covariant and another is the contravariant base vector system. The two base vector systems are reciprocal. The so-called contravariant velocity component is really the contravariant component of a velocity vector for a time-independent coordinate system, or the contravariant component of a relative velocity between fluid and coordinates, for a time-dependent coordinate system. The contravariant velocity components are not physical quantities of the velocity vector. Their magnitudes, dimensions, and associated directions are controlled by their corresponding covariant base vectors. Several 2-D (two-dimensional) linear examples and 2-D mass-conservation equation are used to illustrate the details of expressing a vector with respect to the covariant and contravariant base vector systems, respectively.

Description: The complexity of International Space Station (ISS) systems modeling often necessitates the concurrence of various dissimilar, parallel analysis techniques to validate modeling. This was the case with a feasibility and performance study of the ISS Node 3 Regenerative Heat Exchanger (RHX). A thermo-hydraulic network model was created and analyzed in SINDA/FLUINT. A less complex, closed form solution of the systems dynamics was created using an Excel Spreadsheet. The purpose of this paper is to provide a brief description of the modeling processes utilized, the results and benefits of each to the ISS Node 3 RHX study.

Description: Knowledge of aerothermally induced convective heat transfer is important in the design of thermal protection systems for launch vehicles. Aerothermal models are typically calibrated via the data from circular, in-flight, flush-mounted surface heat flux gauges exposed to the thermal and velocity boundary layers of the external flow. Typically, copper or aluminum Schmidt- Boelter gauges, which take advantage of the one-dimensional Fourier's law of heat conduction, are used to measure the incident heat flux. This instrumentation, when surrounded by low-conductivity insulation, has a wall temperature significantly lower than the insulation. As a result of this substantial disturbance to the thermal boundary layer, the heat flux incident on the gauge tends to be considerably higher than it would have been on the insulation had the calorimeter not been there. In addition, radial conductive heat transfer from the hotter insulation can cause the calorimeter to indicate heat fluxes higher than actual. An overview of an effort to develop and calibrate gauge correction techniques for both of these effects will be presented.

Description: This study provides development and verification of analysis methods used to assess performance of a carbon fiber rope (CFR) thermal barrier system that is currently being qualified for use in Reusable Solid Rocket Motor (RSRM) nozzle joint-2. Modeled geometry for flow calculations considers the joint to be vented with the porous CFR barriers placed in the 'open' assembly gap. Model development is based on a 1-D volume filling approach where flow resistances (assembly gap and CFRs) are defined by serially connected internal flow and the porous media 'Darcy' relationships. Combustion gas flow rates are computed using the volume filling code by assuming a lumped distribution total joint fill volume on a per linear circumferential inch basis. Gas compressibility, friction and heat transfer are included in the modeling. Gas-to-wall heat transfer is simulated by concurrent solution of the compressible flow equations and a large thermal 2-D finite element (FE) conduction grid. The derived numerical technique loosely couples the FE conduction matrix with the compressible gas flow equations. Free constants that appear in the governing equations are calibrated by parametric model comparison to hot fire subscale test results. The calibrated model is then used to make full-scale motor predictions using RSRM aft dome environments. Model results indicate that CFR thermal barrier systems will provide a thermally benign and controlled pressurization environment for the RSRM nozzle joint-2 primary seal activation.

Description: Pure water is distilled from waste water in the International Space Station. The distillation assembly consists of an evaporator, a compressor and a condenser. Vapor is periodically purged from the condenser to avoid vapor accumulation. Purged vapor is condensed in a tube by coolant water prior to entering the purge pump. The paper presents a condensation model of purged vapor in a tube. This model is based on the Finite Volume Method. In the Finite Volume Method, the flow domain is discretized into multiple control volumes and a simultaneous analysis is performed.

Description: This viewgraph presentation evaluates CFD (Computational Fluid Dynamics) tools for solving stage separation problems. The demonstration and validation of the tools is for a second generation RLV (Reusable Launch Vehicle) stage separation. The flow solvers are: Cart3D; Overflow/Overflow-D; Unic.

Description: This viewgraph presentation provides information on computerized simulations of flow fields in a Simplex turbine. The motivations for the simulation were: Determining the effects of partial admission flow on rotor performance as a function of circumferential location and on unsteady rotor loading; Providing an efficient technique for determining turbine performance. The simulation used the flow code CORSAIR.

Description: Contents include the following: Mission Overview. Thermal Environments. Driving Thermal Requirements. Thermal Design Approach. Thermal Control Block Diagram. Thermal Design Description. Thermal Analysis Results Summary. Testing Plans. Issues & Concerns.

Description: This presentation provides an overview of thermal technology development activities carried out at NASA's Goddard Space Flight Center during 2001. Specific topics covered include: two-phase systems (heat pipes, capillary pumped loops, vapor compression systems and phase change materials), variable emittance systems, advanced coatings, high conductivity materials and electrohydrodynamic (EHD) thermal coatings. The application of these activities to specific space missions is also discussed.

Description: An overview of the NASA Glenn Low-Pressure Turbine (LPT) Flow Physics Program will be presented. The flow in the LPT is unique for the gas turbine. It is characterized by low Reynolds number and high freestream turbulence intensity and is dominated by interplay of three basic mechanisms: transition, separation and wake interaction. The flow of most interest is on the suction surface, where large losses are generated due to separation. The LPT is a large, multistage, heavy, jet engine component that suffers efficiency degradation between takeoff and cruise conditions due to decrease in Reynolds number with altitude. The performance penalty is around 2 points for large commercial bypass engines and as much as 7 points for small, high cruise altitude, military engines. The gas-turbine industry is very interested in improving the performance of the LPT and in reducing its weight, part count and cost. Many improvements can be accomplished by improved airfoil design, mainly by increasing the airfoil loading that can yield reduction of airfoils and improved performance. In addition, there is a strong interest in reducing the design cycle time and cost. Key enablers of the needed improvements are computational tools that can accurately predict LPT flows. Current CFD tools in use cannot yet satisfactorily predict the unsteady, transitional and separated flow in the LPT. The main reasons are inadequate transition & turbulence models and incomplete understanding of the LPT flow physics. NASA Glenn has established its LPT program to answer these needs. The main goal of the program is to develop and assess models for unsteady CFD of LPT flows. An approach that consists of complementing and augmenting experimental and computational work elements has been adopted. The work is performed in-house and by several academic institutions, in cooperation and interaction with industry. The program was reviewed at the Minnowbrook II meeting in 1997. This review will summarize the progress that was made since and will introduce newly started projects. The LPT program is focused on three areas: acquisition of experimental and numerical databases and on modeling and computation. Priority was initially given to experiments. There are three classes of experiments: simulated LPT passages, linear cascade, both with and without wakes, and low-speed rotating rig. They are being conducted as follows: At NASA GRC on a flat surface with blade pressure distribution, at the US Naval Academy on a curved surface. The addition of wakes is studied at the University of Minnesota in a curved passage with a retractable wake generator, and at Texas A&M University in a linear cascade with continuously running wake generator. The pressure distribution of the Pratt & Whitney blade 'Pak B' is used in all these experiments. Experiments have been performed also in the GEAE Low-Speed Rotating Turbine (LSRT) rig with GE-designed airfoils. Work on numerically generated database is in progress at the University of Kentucky, using the DNS/LES code LESTool developed there. Turbulence/transition model assessment and development is performed also at the University of Kentucky, where a new intermittency transport model was developed and many experimental test cases have been numerically computed. Assessments of models using simulations of multistage LPT experiments were performed at Virginia Commonwealth University using the Corsair code. Work on suction surface separation delay, using passive and active flow-control, has also been initiated. Following the overview, Principal Investigators attending the workshop will present in detail several of the projects supported by NASA.

Description: Recent measurements on a turbine rotor showed significant relaminarization effects. These effects were evident on the pressure surface heat transfer measurements. The character of the heat transfer varied with Reynolds number. Data were obtained for exit Reynolds numbers between 500,000 and 880,000. Tests were done with a high level of inlet turbulence, 7.5%. At lower Reynolds numbers the heat transfer was similar to that for laminar flow, but at a level higher than for laminar flow. At higher Reynolds numbers the heat transfer was similar to turbulent flow, when the acceleration parameter, K, was sufficiently small. The proposed paper discusses the experimental results, and also discusses approaches to calculating the surface heat transfer for the blade surface. Calculations were done using a three-dimensional Navier-Stokes CFD analysis. The results of these tests, when compared with previous blade tests in the same facility, illustrate modeling difficulties that were encountered in CFD predictions. The two blades were in many ways similar. However, the degree of agreement between the same analysis and the experimental data was significantly different. These differences are highlighted to illustrate where improvements in modeling approaches are needed for transitional flows.

Description: The MCI Engine turbopump supplied the propellants to the main injector. The turbopump consisted of four parts; lox pump, interpropellant seal package (IPS), RP pump and turbine. The thermal analysis was divided into two 2D finite element models; Housing or stationary parts and rotor or rotating parts. Both models were analyzed at the same boundary conditions using SINDA. The housing model consisted of; lox pump housing, ips housing, RP housing, turbine inlet housing, turbine housing, exit guide vane, heat shield and both bearing outer races. The rotor model consisted of the lox impeller; lox end bearing and id race, RP impeller, and RP bearing and id race, shaft and turbine disk. The objectives of the analysis were to (1) verified the original design and recommend modifications to it, (2) submitted a thermal environment to support the structural analysis, (3) support the component and engine test program and (4) to support the X34 vehicle program.

Description: This viewgraph presentation provides information on software development tools to facilitate the use of CFD (Computational Fluid Dynamics) codes. The specific CFD codes FDNS and CORSAIR are profiled, and uses for software development tools with these codes during pre-processing, interim-processing, and post-processing are explained.

Description: The need to identify spikes in the concentration of hazardous gases during countdowns to space shuttle launches has led Kennedy Space Center to acquire considerable expertise in the design, construction, and operation of special-purpose gas analyzers of mass-spectrometer type. If such devices could be miniaturized so as to fit in a small airborne package or backpack them their potential applications would include integrated vehicle health monitoring in later-generation space shuttles and in hazardous material detection in airports, to name two examples. The bulkiest components of such devices are vacuum pumps, particularly those that function in the low vacuum range. Now some pumps that operate in the high vacuum range (e.g. molecular-drag and turbomolecular pumps) are already small and rugged. The present work aims to determine whether, on physical grounds, one may or may not adopt the molecular-drag principle to the low-vacuum range (in which case viscous-drag principle is the appropriate term). The deliverable of the present effort is the derivation and justification of some key formulas and calculation methods for the preliminary design of a single-spool, spiral-channel viscous-drag pump.

Description: A detonation is a supersonic combustion wave. A Pulse Detonation Engine (PDE) repetitively creates a series of detonation waves to take advantage of rapid burning and high peak pressures to efficiently produce thrust. NASA Glenn Research Center's Combustion Branch has developed a PDE test bed that can reproduce the operating conditions that might be encountered in an actual engine. It allows the rapid and cost-efficient evaluation of the technical issues and technologies associated with these engines. The test bed is modular in design. It consists of various length sections of both 2- and 2.6- in. internal-diameter combustor tubes. These tubes can be bolted together to create a variety of combustor configurations. A series of bosses allow instrumentation to be inserted on the tubes. Dynamic pressure sensors and heat flux gauges have been used to characterize the performance of the test bed. The PDE test bed is designed to utilize an existing calorimeter (for heat load measurement) and windowed (for optical access) combustor sections. It uses hydrogen as the fuel, and oxygen and nitrogen are mixed to simulate air. An electronic controller is used to open the hydrogen and air valves (or a continuous flow of air is used) and to fire the spark at the appropriate times. Scheduled tests on the test bed include an evaluation of the pumping ability of the train of detonation waves for use in an ejector and an evaluation of the pollutants formed in a PDE combustor. Glenn's Combustion Branch uses the National Combustor Code (NCC) to perform numerical analyses of PDE's as well as to evaluate alternative detonative combustion devices. Pulse Detonation Engine testbed.

Description: Investigators have used infrared devices to accurately measure heated surface temperatures. Several of these applications have been for turbine heat transfer studies involving film cooling and surface roughness, typically, these measurements use an infrared camera positioned externally to the test section. In cascade studies, where several blades are used to ensure periodic flow, adjacent blades block the externally positioned camera's views of the test blade. To obtain a more complete mapping of the surface temperatures, researchers at the NASA Glenn Research Center fabricated a probe with an infrared detector to sense the blade temperatures. The probe size was kept small to minimize the flow disturbance. By traversing and rotating the probe, using the same approach as for total pressure surveys, one can find the blade surface temperatures. Probe mounted infrared detectors are appropriate for measuring surface temperatures where an externally positioned infrared camera is unable to completely view the test object. This probe consists of a 8-mm gallium arsenide (GaAs) lens mounted in front of a mercury-cadmium-zinc-tellurium (HgCdZnTe) detector. This type of photovoltaic detector was chosen because of its high sensitivity to temperature when the detector is uncooled. The particular application is for relatively low surface temperatures, typically ambient to 100 C. This requires a detector sensitive at long wavelengths. The detector is a commercial product enclosed in a 9-mm-diameter package. The GaAs lens material was chosen because of its glass-like hardness and its good long-wavelength transmission characteristics. When assembled, the 6.4-mm probe stem is held in the traversing actuator. Since the entire probe is above the measurement plane, the flow field disturbance in the measurement plane is minimized. This particular probe body is somewhat wider than necessary, because it was designed to have replaceable detectors and lenses. The signal for the detector is fed through the hollow probe body. The detector's signal goes to an externally mounted preamplifier. The detector assembly, along with a preamplifier, is calibrated as a function of the surface temperature for various detector temperatures. The output voltage is a function of both the detector and object temperatures.

Description: The spreading of an evaporating liquid on a solid surface occurs in many practical processes and is of importance in a number of practical situations such as painting, textile dyeing, coating, gluing, and thermal engineering. Typical processes involving heat transfer where the contact angle plays an important role are film cooling, boiling, and the heat transfer through heat pipes. The biological phenomenon of cell spreading also is analogous to a drop spreading (ref. 1). In the study of spreading, the dynamic contact angle describes the interfacial properties on solid substrates and, therefore, has been studied by physicists and fluid mechanics investigators. The dynamic contact angle of a spreading nonvolatile liquid drop provides a simple tool in the study of the free-boundary problem, but the study of the spreading of a volatile liquid drop is of more practical interest because the evaporation of common liquids is inevitable in practical processes. The most common method to measure the contact angle, the contact radius, and the height of a sessile drop on a solid surface is to view the drop from its edge through an optical microscope. However, this method gives only local information in the view direction. Zhang and Yang (ref. 2) developed a laser shadowgraphy method to investigate the evaporation of sessile drop on a glass plate. As described here, Zhang and Chao (refs. 3 and 4) improved the method and suggested a new optical arrangement to measure the dynamic contact angle and the instant evaporation rate of a sessile drop with much higher accuracy (less than 1 percent). With this method, any fluid motion in the evaporating drop can be visualized through shadowgraphy without using a tracer, which often affects the field under investigation.

Description: Sprays occur in a wide variety of industrial and power applications and in materials processing. A liquid spray is a two-phase flow with a gas as the continuous phase and a liquid as the dispersed phase in the form of droplets or ligaments. The interactions between the two phases--which are coupled through exchanges of mass, momentum, and energy--can occur in different ways at disparate time and length scales involving various thermal, mass, and fluid dynamic factors. An understanding of the flow, combustion, and thermal properties of a rapidly vaporizing spray requires careful modeling of the ratecontrolling processes associated with turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates of the spray, among many other factors. With the aim of developing an efficient solution procedure for use in multidimensional combustor modeling, researchers at the NASA Glenn Research Center have advanced the state-of-the-art in spray computations in several important ways.

Description: Microgravity poses many challenges to the designer of spacecraft tanks. Chief among these are the lack of phase separation and the need to supply vapor-free liquid or liquidfree vapor to the spacecraft processes that require fluid. One of the principal problems of phase separation is the creation of liquid jets. A jet can be created by liquid filling, settling of the fluid to one end of the tank, or even closing a valve to stop the liquid flow. Anyone who has seen a fountain knows that jets occur in normal gravity also. However, in normal gravity, the gravity controls and restricts the jet flow. In microgravity, with gravity largely absent, surface tension forces must contain jets. To model this phenomenon, a numerical method that tracks the fluid motion and the surface tension forces is required. Jacqmin has developed a phase model that converts the discrete surface tension force into a barrier function that peaks at the free surface and decays rapidly away. Previous attempts at this formulation were criticized for smearing the interface. This can be overcome by sharpening the phase function, double gridding the fluid function, and using a higher order solution for the fluid function. The solution of this equation can be rewritten as two coupled Poisson equations that also include the velocity. After the code was implemented in axisymmetric form and verified by several test cases at the NASA Glenn Research Center, the drop tower runs of Aydelott were modeled. Work last year with a laminar model was found to overpredict Aydelott's results, except at the lowest Reynolds number conditions of 400. This year, a simple turbulence model was implemented by adding a turbulent viscosity based on the mixing-length hypothesis and empirical measurements of previous works. Predictions made after this change was implemented have been much closer to experimentally observed flow patterns and geyser heights. Two model runs is shown. The first, without any turbulence correction, breaks through the free surface and strikes the far end of the tank. In the second, the turbulence spreads the jet momentum over more of the free surface, enabling the surface tension forces to turn the jet back into the bulk liquid. The model geyser height with the second model is 1.1 cm. This is quite close to the 1.5-cm geyser height measured by Aydelott.

Description: A computational model was developed at the NASA Glenn Research Center to investigate possible uses of vortex generators (VG's) for improving the performance of turbomachinery. A vortex generator is a small, winglike device that generates vortices at its tip. The vortices mix high-speed core flow with low-speed boundary layer flow and, thus, can be used to delay flow separation. VG's also turn the flow near the walls and, thus, can be used to control flow incidence into a turbomachinery blade row or to control secondary flows.

Description: Aerodynamic stability is an important parameter in the design of compressors for aircraft gas turbine engines. Compression system instabilities can cause compressor surge, which may lead to the loss of an aircraft. As a result, engine designers include a margin of safety between the operating line of the engine and the stability limit line of the compressor. The margin of safety is typically referred to as "surge margin." Achieving the highest possible level of surge margin while meeting design point performance objectives is the goal of the compressor designer. However, performance goals often must be compromised in order to achieve adequate levels of surge margin. Techniques to improve surge margin will permit more aggressive compressor designs. Centrifugal compressor surge margin improvement was demonstrated at the NASA Glenn Research Center by injecting air into the vaned diffuser of a 4:1-pressure-ratio centrifugal compressor. Tests were performed using injector nozzles located on the diffuser hub surface of a vane-island diffuser in the vaneless region between the impeller trailing edge and the diffuser-vane leading edge. The nozzle flow path and discharge shape were designed to produce an air stream that remained tangent to the hub surface as it traveled into the diffuser passage. Injector nozzles were located near the leading edge of 23 of the 24 diffuser vanes. One passage did not contain an injector so that instrumentation located in that passage would be preserved. Several orientations of the injected stream relative to the diffuser vane leading edge were tested over a range of injected flow rates. Only steady flow (nonpulsed) air injection was tested. At 100 percent of the design speed, a 15-percent improvement in the baseline surge margin was achieved with a nozzle orientation that produced a jet that was bisected by the diffuser vane leading edge. Other orientations also improved the baseline surge margin. Tests were conducted at speeds below the design speed, and similar results were obtained. In most cases, the greatest improvement in surge margin occurred at fairly low levels of injected flow rate. Externally supplied injection air was used in these experiments. However, the injected flow rates that provided the greatest benefit could be produced using injection air that is recirculating between the diffuser discharge and nozzles located in the diffuser vaneless region. Future experiments will evaluate the effectiveness of recirculating air injection.

Description: An optical measurement technique known as Three-Dimensional Digital Particle Image Velocimetry (3-D DPIV) was used to characterize the tip clearance flow in NASA Glenn Research Center's low-speed axial compressor. 3-D DPIV is a technique in which a stereoscopic imaging system consisting of two cross-correlation cameras is used to record particles entrained in a flow as a laser light sheet is pulsed at two instances in time. Although 3-D DPIV has been used elsewhere, this is the first time it has been used to measure compressor tip clearance flows. In-house modifications of the DPIV system include the use of effective seeding technology and a novel system to perform a priori calibrations at all five measurement planes, greatly reducing facility run time. Computational fluid dynamics predictions, which are used to guide design changes toward improving the efficiency and operating range of turbomachinery, can be verified and improved by comparison with 3-D DPIV measurements of the actual tip clearance flow. This measurement campaign dealt with the characterization of the tip clearance vortex in the first stage of a four-stage axial compressor. The tip clearance vortex is formed in compressors operating with a clearance gap between the moving rotor blade tips and the stationary casing when a leakage flow, forced from the pressure side of the blade over the blade tip, forms a vortical structure on the suction side of the blade. 3-D DPIV is ideally suited to measure the clearance vortex for two reasons: (1) this technique captures the entire blade passage flow at one instant in time, so that wandering of the vortex during the measurement does not smear out velocity gradients in the flow field, and (2) the spanwise component of velocity changes sign across the vortex core, providing a more accurate measurement of the vortex location than was available with previous two-dimensional measurement approaches. These two attributes of the data will enable computational fluid dynamics researchers to validate their predictions to a level of accuracy not previously attainable. In an effort to understand the effects of this flow phenomenon on the operation of the low speed axial compressor, data were acquired at two mass flow coefficients: 0.395 (design operating point) and 0.35 (operating point just above stall). In order to identify the spatial extent, location, and magnitude of the tip clearance vortex for each mass flow condition, data were acquired at five equally spaced spanwise locations, from 90- to 100-percent span. The data presented shows a cross section of the tip clearance vortex at 94-percent span and a mass flow coefficient of 0.395, where the color contours represent the radial velocity component. The blue (flow towards the hub) and red (flow towards the casing) regions illustrate the extent of the tip clearance vortex. The interface between the blue and red regions delineates the location of the vortex core.

Description: An advanced high-speed fan was recently designed under a cooperative effort between the NASA Glenn Research Center and Honeywell Engines & Systems. The principal design goals were to improve performance and to reduce fan noise at takeoff. Scale models of the Quiet High-Speed Fan were tested for operability, performance, and acoustics. During testing, the fan showed significantly improved noise characteristics, but a self-excited aeroelastic vibration known as flutter was encountered in the operating range. Flutter calculations were carried out for the Quiet High-Speed Fan using a three-dimensional, unsteady aerodynamic, Reynolds-averaged Navier-Stokes turbomachinery code named "TURBO." The TURBO code can accurately model the viscous flow effects that can play an important role in various aeroelastic problems such as flutter with flow separation, flutter at high loading conditions near the stall line (stall flutter), and flutter in the presence of shock and boundary-layer interaction. Initially, calculations were performed with no blade vibrations. These calculations were at a constant rotational speed and a varying mass flow rate. The mass flow rate was varied by changing the backpressure at the exit boundary of the computational domain. These initial steady calculations were followed by aeroelastic calculations in which the blades were prescribed to vibrate harmonically in a natural mode, at a natural frequency, and with a fixed interblade phase angle between adjacent blades. The AE-prep preprocessor was used to interpolate the in-vacuum mode shapes from the structural dynamics mesh onto the computational fluid dynamics mesh and to smoothly propagate the grid deformations from the blade surface to the interior points of the grid. The aeroelastic calculations provided the unsteady aerodynamic forces on the blade surface due to blade vibrations. These forces were vector multiplied with the structural dynamic mode shape to calculate the work done on the blade during one vibration period, then this result was converted to an aerodynamic damping. Flutter occurs when the aerodynamic damping becomes negative, if structural damping is ignored. The results of these aeroelastic calculations are summarized in this plot of aerodynamic damping versus mass flow rate at a constant rotational speed. As the backpressure is increased, the mass flow rate through the fan decreases and the fan operating point moves towards the stall line. The aeroelastic calculations showed that the aerodynamic damping decreases as the stall line is approached, as observed during testing. In addition, the aeroelastic calculations with the TURBO code correctly predicted the aeroelastic parameters: the most unstable vibration mode and interblade phase angle, as observed during testing. The Quiet High-Speed Fan demonstrated significant noise reductions during testing, but flutter imposed limits on its operating range. The accurate calculation of the aeroelastic characteristics using the TURBO code is a significant step toward eliminating flutter from the operating range and toward realizing the benefits of reduced fan noise. The aeroelastic calculations described here were performed under a grant by University of Toledo researchers in collaboration with Glenn's researchers.

Description: An analytical solution to the steady-state fluid temperature for 1-D (one dimensional) transpiration cooling has been derived. Transpiration cooling has potential use in the aerospace industry for protection against high heating environments for re-entry vehicles. Literature for analytical treatments of transpiration cooling has been largely confined to the assumption of thermal equilibrium between the porous matrix and fluid. In the present analysis, the fundamental fluid and matrix equations are coupled through a volumetric heat transfer coefficient and investigated in non-thermal equilibrium. The effects of varying the thermal conductivity of the solid matrix and the heat transfer coefficient are investigated. The results are also compared to existing experimental data.

Description: During a post-test inspection of a Booster Separation Motor (BSM) from a Lot Acceptance Test (LAT), a crack was noticed in the graphite throat. Since this was an out-of-family occurrence, an investigation team was formed to determine the cause of the crack. This paper will describe thermal analysis techniques used in support of this investigation. Models were generated to predict gradients in nominal motor conditions, as well as potentially anomalous conditions. Analysis was also performed on throats that were tested in the Laser Hardened Material Evaluation Laboratory (LHMEL). Some of these throats were pre-cracked, while others represented configurations designed to amplify effects of thermal stresses. Results from these analyses will be presented in this paper.

Description: Thermogravimetric analysis (TGA) is widely employed in the thermal characterization of non-metallic materials, yielding valuable information on decomposition characteristics of a sample over a wide temperature range. However, a potential wealth of chemical information is lost during the process, with the evolving gases generated during thermal decomposition escaping through the exhaust line. Fourier Transform-Infrared spectroscopy (FT-IR) is a powerful analytical technique for determining many chemical constituents while in any material state, in this application, the gas phase. By linking these two techniques, evolving gases generated during the TGA process are directed into an appropriately equipped infrared spectrometer for chemical speciation. Consequently, both thermal decomposition and chemical characterization of a material may be obtained in a single sample run. In practice, a heated transfer line is employed to connect the two instruments while a purge gas stream directs the evolving gases into the FT-IR. The purge gas can be either high purity air or an inert gas such as nitrogen to allow oxidative and pyrolytic processes to be examined, respectively. The FT-IR data is collected realtime, allowing continuous monitoring of chemical compositional changes over the course of thermal decomposition. Using this coupled technique, an array of diverse materials has been examined, including composites, plastics, rubber, fiberglass epoxy resins, polycarbonates, silicones, lubricants and fluorocarbon materials. The benefit of combining these two methodologies is of particular importance in the aerospace community, where newly developing materials have little available data with which to refer. By providing both thermal and chemical data simultaneously, a more definitive and comprehensive characterization of the material is possible. Additionally, this procedure has been found to be a viable screening technique for certain materials, with the generated data useful in the selection of other appropriate analytical procedures for further material characterization.

Description: Contents include the following: Quench Module Insert (QMI) science requirements. QMI interfaces. QMI design layout. QMI thermal analysis and design methodology. QMI bread board testing and instrumentation approach. QMI thermal probe design parameters. Design features for gradient measurement. Design features for heated zone measurements. Thermal gradient analysis results. Heated zone analysis results. Bread board thermal probe layout. QMI bread board correlation and performance. Summary and conclusions.

Description: The Marshall Space Flight Center (MSFC) is hosting the Thermal and Fluids Analysis Workshop (TFAWS) during the week of September 10, 2001. Included in this year's TFAWS is a panel session on Multidisciplinary Analysis techniques. The intent is to provide an opportunity for the users to gain information as to what product may be best suited for their applications environment and to provide feedback to you, the developers, on future desired developments. Potential users of multidisciplinary analysis (MDA) techniques are often overwhelmed by the number of choices available to them via commercial products and by the pace of new developments in this area. The purpose of this panel session is to provide a forum wherein MDA tools available and under development can be discussed, compared, and contrasted. The intent of this panel is to provide the end-user with the information necessary to make educated decisions on how to proceed with selecting their MDA tool. It is anticipated that the discussions this year will focus on MDA techniques that couple discipline codes or algorithms (as opposed to monolithic, unified MDA approaches). The MDA developers will be asked to prepare a product overview presentation addressing specific questions provided by the panel organizers. The purpose of these questions will be to establish the method employed by the particular MDA technique for communication between the discipline codes, to establish the similarities and differences amongst the various approaches, and to establish the range of experience and applications for each particular MDA approach.

Description: The prestart thermal conditioning of the hardware in LOX (liquid oxygen) systems involve heat transfer between LOX and metal where boiling plays a large role. Information is easily found on nucleate boiling, maximum heat flux, minimum heat flux and film boiling for common fluids like water. After looking at these standard correlations it was felt more data was needed for the cool down side transition boiling for the LN2 and LOX. In particular interest is the film boiling values, the temperature at which transition begins and the slope as peak heat flux is approached. The ultimate goal is an array of boiling heat transfer coefficient as a function of surface temperature which can be used in the chilldown model of the feed system, engine and bleed system for X-34. The first experiment consisted of an actual MC-1 LOX Impeller which had been machined backwards, that was instrumented with 17 surface thermocouples and submerged in liquid nitrogen. The thermocouples were installed on metal thicknesses varying from the thin inducer to the thick hub.

Description: As part of efforts to design a regeneratively cooled composite nozzle ramp for use on the reusable vehicle (RLV) rocket engine, an C-SiC composites heat exchanger concept was proposed for thermal performance evaluation. To test the feasibility of the concept, sample heat exchanger panels were made to fit the Glenn Research Center's cell 22 for testing. Operation of the heat exchanger was demonstrated in a combustion environment with high heat fluxes similar to the RLV Aerospike Ramp. Test measurements were reviewed and found to be valuable for the on going fluid and thermal analysis of the actual RLV composite ramp. Since the cooling fluid for the heat exchanger is water while the RLV Ramp cooling fluid is LH2, fluid and thermal models were constructed to correlate to the specific test set-up. The knowledge gained from this work will be helpful for analyzing the thermal response of the actual RLV Composite Ramp. The coolant thermal properties for the models are taken from test data. The heat exchanger's cooling performance was analyzed using the Generalized Fluid System Simulation Program (GFSSP). Temperatures of the heat exchanger's structure were predicted in finite element models using Patran and Sinda. Results from the analytical models and the tests show that RSC's heat exchanger satisfied the combustion environments in a series of 16 tests.

Description: Two ISS (International Space Station) experiment payloads will vent a volume of gas overboard via either the ISS Vacuum Exhaust System or the Vacuum Resource System. A system of ducts, valves and sensors, under design, will connect the experiments to the ISS systems. The following tasks are required: Create an analysis tool that will verify the rack vacuum system design with respect to design requirements, more specifically approximate pressure at given locations within the vacuum systems; Determine the vent duration required to achieve desired pressure within the experiment modules; Update the analysis as systems and operations definitions mature.

Description: Unsteady flow simulations for RLV (Reusable Launch Vehicles) 2nd Generation baseline turbopump for one and half impeller rotations have been completed by using a 34.3 Million grid points model. MLP (Multi-Level Parallelism) shared memory parallelism has been implemented in INS3D, and benchmarked. Code optimization for cash based platforms will be completed by the end of September 2001. Moving boundary capability is obtained by using DCF module. Scripting capability from CAD (computer aided design) geometry to solution has been developed. Data compression is applied to reduce data size in post processing. Fluid/Structure coupling has been initiated.

Description: Modern low-pressure turbine airfoils are subject to increasingly stronger pressure gradients as designers impose higher loading in an effort to improve efficiency and to reduce part count. The adverse pressure gradients on the suction side of these airfoils can lead to boundary-layer separation, particularly under cruise conditions. Separation bubbles, notably those which fail to reattach, can result in a significant degradation of engine efficiency. Accurate prediction of separation and reattachment is hence crucial to improved turbine design. This requires an improved understanding of the transition flow physics. Transition may begin before or after separation, depending on the Reynolds number and other flow conditions, has a strong influence on subsequent reattachment, and may even eliminate separation. Further complicating the problem are the high free-stream turbulence levels in a real engine environment, the strong pressure gradients along the airfoils, the curvature of the airfoils, and the unsteadiness associated with wake passing from upstream stages. Because of the complicated flow situation, transition in these devices can take many paths that can coexist, vary in importance, and possibly also interact, at different locations and instances in time. The present work was carried out in an attempt to systematically sort out some of these issues. Detailed velocity measurements were made along a flat plate subject to the same nominal dimensionless pressure gradient as the suction side of a modern low-pressure turbine airfoil ('Pak-B'). The Reynolds number based on wetted plate length and nominal exit velocity, Re, was varied from 50;000 to 300; 000, covering cruise to takeoff conditions. Low, 0.2%, and high, 7%, inlet free-stream turbulence intensities were set using passive grids. These turbulence levels correspond to about 0.2% and 2.5% turbulence intensity in the test section when normalized with the exit velocity. The Reynolds number and free-stream turbulence level do not have a significant effect on the location of boundary-layer separation unless they are high enough to induce transition upstream of separation. The location and extent of the transition zone, in contrast, depend strongly on Re and TI. The beginning of reattachment closely follows the onset of transition. Under low free-stream turbulence conditions the boundary layer is laminar at separation and then begins to exhibit fluctuations in a finite frequency band in the shear layer over the separation bubble. These fluctuations are due to instability waves. The fluctuations grow in magnitude, higher harmonics are generated, and finally lead to a breakdown to turbulence. Transition begins in the shear layer, but quickly spreads to the near wall region and causes the boundary layer to reattach. The transition is rapid and the resulting turbulence contains a full range of high and low frequencies. Under high free-stream turbulence conditions, slowly growing low-frequency fluctuations are induced in the pretransitional boundary layer by the free-stream. The separation bubbles are considerably thinner than in the low TI cases, resulting in thinner boundary layers at the end of the test wall. At Re=50,000 and 100,000, the pre-transitional boundary layer separates at about the same location as in the low TI cases. Transition occurs through a bypass mode, begins upstream of the corresponding low-TI location, and proceeds in a manner similar to that of an attached boundary layer. Under high TI at Re=200,000 and 300,000, transition begins before separation. The boundary layer may separate, but if it does the separation bubble is very short and does not significantly affect the downstream development of the boundary layer. A comparison is made to previous work in a simulated cascade.

Description: In the typical operation of a packed-bed reactor, gas and liquid flow simultaneously through a fixed bed of solid particles. Depending on the application, the particles can be of various shapes and sizes and provide for intimate contact and high rates of transport between the phases needed to sustain chemical or biological reactions. The packing may also serve as either a catalyst or as a support for growing biological material. NASA has flown two of these packed-bed systems in a microgravity environment with limited or no success. The goal of this research is to develop models (with scale-up capability) needed for the design of the physicochemical equipment to carry out these unit operations in microgravity. New insight will also lead to improvements in normal gravity operations. Our initial experiment was flown using an existing KC-135 two-phase flow rig with a modified test section. The test section is a clear polycarbonate rectangular column with a depth of 2.54 cm, a width of 5.08 cm, and 60 cm long. The column was randomly packed with spherical glass beads by slowly dropping the beads into the bed. Even though care was taken in handling the column after it was filled with packing, the alternating high and low gravity cycles with each parabola created a slightly tighter packed bed than is typically reported for this type. By the usual method of comparing the weight difference of a completely dry column versus a column filled with water, the void fraction was found to be .345 for both sizes of beads used. Five flush mounted differential pressure transducers are spaced at even intervals with the first location 4 cm from the inlet port and the subsequent pressure transducers spaced at 13 cm intervals along the column. Differential pressure data was acquired at 1000 Hz to adequately observe pulse formation and characteristics. Visual images of the flow were recorded using a high-speed SVHS system at 500 frames per second. Over 250 different test conditions were evaluated along with a companion set of tests in normal gravity. The flow rates, fluid properties and packing properties were selected to provide a range of several orders-of-magnitude for the important dimensionless parameters. Additional information is included in the original extended abstract.

Description: Aluminum-Lithium (AL-Li) alloys offer significant performance benefits for aerospace structural applications due to their higher specific properties compared with conventional aluminum alloys. For example, the application of an Al-Li alloy to the space shuttle external cryogenic fuel tank contributed to the weight savings that enabled successful deployment of International Space Station components. The composition and heat treatment of this alloy were optimized specifically for strength-toughness considerations for an expendable cryogenic tank. Time dependent properties related to reliability, such as thermal stability, fatigue, and corrosion, will be of significant interest when materials are evaluated for a reusable cryotank structure. As most aerospace structural hardware is weight sensitive, a reusable cryotank will be designed to the limits of the materials mechanical properties. Therefore, this effort was designed to establish the effects of thermal exposure on the mechanical properties and microstructure of one relatively production mature alloy and two developmental alloys C458 and L277. Tensile and fracture toughness behavior was evaluated after exposure to temperatures as high as 3oooF for up to IO00 hrs. Microstructural changes were also evaluated to correlate with the observed data trends. The ambient temperature parent metal data showed an increase in strength and reduction in elongation after exposure at lower temperatures. Strength reached a peak with intermediate temperature exposure followed by a decrease at highest exposure temperature. Characterizing the effect of thermal exposure on the properties of Al-Li alloys is important to defining a service limiting temperature, exposure time, and end-of-life properties.

Description: In a propulsion system, the propellant flow and mixture ratio could be controlled either by variable area flow control valves or by passive flow control elements such as cavitating venturies. Cavitating venturies maintain constant propellant flowrate for fixed inlet conditions (pressure and temperature) and wide range of outlet pressures, thereby maintain constant, engine thrust and mixture ratio. The flowrate through the venturi reaches a constant value and becomes independent of outlet pressure when the pressure at throat becomes equal to vapor pressure. In order to develop a numerical model of propulsion system, it is necessary to model cavitating venturies in propellant feed systems. This paper presents a finite volume model of flow network of a cavitating venturi. The venturi was discretized into a number of control volumes and mass, momentum and energy conservation equations in each control volume are simultaneously solved to calculate one-dimensional pressure, density, and flowrate and temperature distribution. The numerical model predicts cavitations at the throat when outlet pressure was gradually reduced. Once cavitation starts, with further reduction of downstream pressure, no change in flowrate is found. The numerical predictions have been compared with test data and empirical equation based on Bernoulli's equation.

Description: The properties of molten semiconductors are good indicators for material structure transformation and hysteresis under temperature variations. Viscosity, as one of the most important properties, is difficult to measure because of high temperature, high pressure, and vapor toxicity of melts. Recently, a novel method was developed by applying a rotating magnetic field to the melt sealed in a suspended quartz ampoule, and measuring the transient torque exerted by rotating melt flow on the ampoule wall. The method was designed to measure viscosity in short time period, which is essential for evaluating temperature hysteresis. This paper compares the theoretical prediction of melt flow and ampoule oscillation with the experimental data. A theoretical model was established and the coupled fluid flow and ampoule torsional vibration equations were solved numerically. The simulation results showed a good agreement with experimental data. The results also showed that both electrical conductivity and viscosity could be calculated by fitting the theoretical results to the experimental data. The transient velocity of the melt caused by the rotating magnetic field was found reach equilibrium in about half a minute, and the viscosity of melt could be calculated from the altitude of oscillation. This would allow the measurement of viscosity in a minute or so, in contrast to the existing oscillation cup method, which requires about an hour for one measurement.

Description: This paper reports the recent research results of shockwave structure predictions using a new developed code. The modified Rankine-Hugoniot relations across a standing normal shock wave are discussed and adopted to obtain jump conditions. Coupling a electrostatic body force to the Burnett equations, the weakly ionized flow field across the shock wave was solved. Results indicated that the Modified Rankine-Hugoniot equations for shock wave are valid for a wide range of ionization fraction. However, this model breaks down with small free stream Mach number and with large ionization fraction. The jump conditions also depend on the value of free stream pressure, temperature and density. The computed shock wave structure with ionization provides results, which indicated that shock wave strength may be reduced by existence of weakly ionized gas.

Description: This paper reports on ground based work conducted to support the Spaceflight Definition project SHIVA (Spaceflight Holography Investigation in a Virtual Apparatus). SHIVA will advance our understanding of the movement of a particle in a fluid. Gravity usually dominates the equations of motion, but in microgravity as well as on earth other terms can become important. Before two members of our team found an analytical solution of the equations, numerical methods and/or neglecting terms were required. The general solution predicts that the usually neglected history term becomes important when the characteristic viscous time is in the same order as the vibration period and peaks when the two times are equal. In this case three force terms, the Stokes drag, the added mass, and the history drag must all be included in predicting particle movement. We also developed diagnostic recording methods using holography to save all of the particle field data, allowing the experiment to essentially be transferred from space back to earth in what we call the "virtual apparatus". Using state-of-the-art methods in holography we will quantify the three-dimensional motion of sets of particles, allowing us to test and apply the new analytical solutions. The motion of particles up to 4 mm in diameter in a fluid that oscillates at frequencies up to 100 Hz with amplitudes up to 200 microns is being examined. Ground studies to support the flight development program have employed various schemes to simulate microgravity. One of the most reliable and meaningful methods uses spheres tethered to a fine hair suspended in the fluid. We have also investigated particles with nearly neutral buoyancy. Recordings are made at the peak amplitudes of vibration of the cell providing a measure of the ratio of fluid to particle amplitude. The experiment requires precise location of the particle at the time of recording. The hologram of the particle provides microscopic images of the particle that are used for finding the position with an accuracy of a few microns. To make the experiment more versatile, the spaceflight system will record holograms both on film and electronically. The electronic holograms can be downlinked providing real time data. Results of the ground experiments, the flight experiment design, and data analysis procedures are reported.

Description: While there is a general understanding of the acceleration environment onboard an orbiting spacecraft, past research efforts in the modeling and analysis area have still not produced a general theory that predicts the effects of multi-spectral periodic accelerations on a general class of experiments nor have they produced scaling laws that a prospective experimenter can use to assess how an experiment might be affected by this acceleration environment. Furthermore, there are no actual flight experimental data that correlates heat or mass transport with measurements of the periodic acceleration environment. The present investigation approaches this problem with carefully conducted terrestrial experiments and rigorous numerical modeling for better understanding the effect of residual gravity and gentler on experiments. The approach is to use magnetic fluids that respond to an imposed magnetic field gradient in much the same way as fluid density responds to a gravitational field. By utilizing a programmable power source in conjunction with an electromagnet, both static and dynamic body forces can be simulated in lab experiments. The paper provides an overview of the technique and includes recent results from the experiments.

Description: The measurement of 3-D three-component velocity fields is of great importance in both ground and space experiments for understanding materials processing and fluid physics. Here, we present the investigation results of stereoscopic tracking velocimetry (STV) for measuring 3-D velocity fields. The effort includes diagnostic technology development, experimental velocity measurement, and comparison with analytical and numerical computation. The advantages of STV stems from the system simplicity for building compact hardware and in software efficiency for continual near-real-time process monitoring. It also has illumination flexibility for observing volumetric flow fields from arbitrary directions. STV is based on stereoscopic CCD observations of particles seeded in a flow. Neural networks are used for data analysis. The developed diagnostic tool is tested with a simple directional solidification apparatus using Succinonitrile. The 3-D velocity field in the liquid phase is measured and compared with results from detailed numerical computations. Our theoretical, numerical, and experimental effort has shown STV to be a viable candidate for reliably quantifying the 3-D flow field in materials processing and fluids experiments.

Description: A localized precipitation network has been established over a 4000 sq km region of northern Alabama in support of local weather and climate research at the Global Hydrology and Climate Center (GHCC) in Huntsville. This Cooperative Huntsville-Area Rainfall Measurement (CHARM) network is comprised of over 80 volunteers who manually take daily rainfall measurements from 85 sites. The network also incorporates 20 automated gauges that report data at 1-5 minute intervals on a 24 h a day basis. The average spacing of the gauges in the network is about 6 kin, however coverage in some regions benefit from gauges every 1-2 km. The 24 h rainfall totals from the CHARM network have been used to validate Stage III rainfall estimates of daily and storm totals derived from the WSR-88D radars that cover northern Alabama. The Stage III rainfall product is produced by the Lower Mississippi River Forecast Center (LMRFC) in support of their daily forecast operations. The intercomparisons between the local rain gauge and the radar estimates have been useful to understand the accuracy and utility of the Stage III data. Recently, the Stage III and CHARM rainfall measurements have been combined to produce an hourly rainfall dataset at each CHARM observation site. The procedure matches each CHARM site with a time sequence of Stage III radar estimates of precipitation. Hourly stage III rainfall estimates were used to partition the rain gauge values to the time interval over which they occurred. The new hourly rain gauge dataset is validated at selected points where 1-5 minute rainfall measurements have been made. This procedure greatly enhances the utility of the CHARM data for local weather and hydrologic modeling studies. The conference paper will present highlights of the Stage III intercomparison and some examples of the combined radar / rain gauge product demonstrating its accuracy and utility in deriving an hourly rainfall product from the 24 h CHARM totals.

Description: This paper treats the axisymmetric flow and mass transport in a cylindrical liquid bridge between the melting end of a feed rod and the solidifying end of an alloyed semiconductor crystal. There is a strong, uniform, steady, axial magnetic field. The surface tension depends on the temperature and the concentration of the species, while variations of the concentration occur because one species is rejected into the liquid during solidification. The thermocapillary and solutocapillary convections tend to cancel- over part of the liquid bridge. For certain parameter ranges, there are two different stable steady solutions: one where the concentration gradient along the free surface leads to dominance by the solutocapillary convection and one where the mass transport due to the thermocapillary convection makes the concentration gradient along the free surface small, so that the thermocapillary convection is dominant.

Description: The Convection And Moisture Experiment (CAMEX) is a series of field research investigations sponsored by the Earth Science Enterprise of the National Aeronautics and Space Administration (NASA). The fourth field campaign in the CAMEX series (CAMEX-4) was recently conducted during 16 August - 24 September 2001 using the Jacksonville Naval Air Station in Florida as the main base of operations. CAMEX-4 focused on the study of tropical cyclone (hurricane) development, tracking, intensification, and landfalling impacts using NASA-funded aircraft and surface remote sensing instrumentation. The results of this study will be used to address key issues pertinent to a larger NASA ESE study of the global water cycle as well as to provide synergistic contributions to the research goals of the Hurricane Research Division (HRD) of the National Oceanic and Atmospheric Administration (NOAA) and the Hurricanes At Landfall Initiative of the United States Weather Research Program. All CAMEX-4 aircraft missions were planned and jointly conducted with NOAA aircraft to insure comprehensive sampling. An overview of preliminary observations of Tropical Storms Chantal and Gabrielle as well as Hurricanes Erin and Humberto will be presented.

Description: Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes for the compressible Euler and Navier-Stokes equations has been developed and verified by the authors and collaborators. These schemes are suitable for the problems in question. Basically, the scheme consists of sixth-order or higher non-dissipative spatial difference operators as the base scheme. To control the amount of numerical dissipation, multiresolution wavelets are used as sensors to adaptively limit the amount and to aid the selection and/or blending of the appropriate types of numerical dissipation to be used. Magnetohydrodynamics (MHD) waves play a key role in drag reduction in highly maneuverable high speed combat aircraft, in space weather forecasting, and in the understanding of the dynamics of the evolution of our solar system and the main sequence stars. Although there exist a few well-studied second and third-order high-resolution shock-capturing schemes for the MHD in the literature, these schemes are too diffusive and not practical for turbulence/combustion MHD flows. On the other hand, extension of higher than third-order high-resolution schemes to the MHD system of equations is not straightforward. Unlike the hydrodynamic equations, the inviscid MHD system is non-strictly hyperbolic with non-convex fluxes. The wave structures and shock types are different from their hydrodynamic counterparts. Many of the non-traditional hydrodynamic shocks are not fully understood. Consequently, reliable and highly accurate numerical schemes for multiscale MHD equations pose a great challenge to algorithm development. In addition, controlling the numerical error of the divergence free condition of the magnetic fields for high order methods has been a stumbling block. Lower order methods are not practical for the astrophysical problems in question. We propose to extend our hydrodynamics schemes to the MHD equations with several desired properties over commonly used MHD schemes.

Description: The rapid increase in available computational power over the last decade has enabled higher resolution flow simulations and more widespread use of unstructured grid methods for complex geometries. While much of this effort has been focused on steady-state calculations in the aerodynamics community, the need to accurately predict off-design conditions, which may involve substantial amounts of flow separation, points to the need to efficiently simulate unsteady flow fields. Accurate unsteady flow simulations can easily require several orders of magnitude more computational effort than a corresponding steady-state simulation. For this reason, techniques for improving the efficiency of unsteady flow simulations are required in order to make such calculations feasible in the foreseeable future. The purpose of this work is to investigate possible reductions in computer time due to the choice of an efficient time-integration scheme from a series of schemes differing in the order of time-accuracy, and by the use of more efficient techniques to solve the nonlinear equations which arise while using implicit time-integration schemes. This investigation is carried out in the context of a two-dimensional unstructured mesh laminar Navier-Stokes solver.

Description: The results from the first AIAA CFD Drag Prediction Workshop are summarized. The workshop was designed specifically to assess the state-of-the-art of computational fluid dynamics methods for force and moment prediction. An impartial forum was provided to evaluate the effectiveness of existing computer codes and modeling techniques, and to identify areas needing additional research and development. The subject of the study was the DLR-F4 wing-body configuration, which is representative of transport aircraft designed for transonic flight. Specific test cases were required so that valid comparisons could be made. Optional test cases included constant-C(sub L) drag-rise predictions typically used in airplane design by industry. Results are compared to experimental data from three wind tunnel tests. A total of 18 international participants using 14 different codes submitted data to the workshop. No particular grid type or turbulence model was more accurate, when compared to each other, or to wind tunnel data. Most of the results overpredicted C(sub Lo) and C(sub Do), but induced drag (dC(sub D)/dC(sub L)(exp 2)) agreed fairly well. Drag rise at high Mach number was underpredicted, however, especially at high C(sub L). On average, the drag data were fairly accurate, but the scatter was greater than desired. The results show that well-validated Reynolds-Averaged Navier-Stokes CFD methods are sufficiently accurate to make design decisions based on predicted drag.

Description: Two-dimensional CFD analysis for iced airfoils can be a labor-intensive task. The software toolkit SmaggIce 2D is being developed to help streamline the CFD process and provide the unique features needed for icing. When complete, it will include a combination of partially automated and fully interactive tools for all aspects of the tasks leading up to the flow analysis: geometry preparation, domain decomposition, block boundary discretization. gridding, and linking with a flow solver. It also includes tools to perform ice shape characterization, an important aid in determining the relationship between ice characteristics and their effects on aerodynamic performance. Completed tools, work-in-progress, and planned features of the software toolkit are presented here.

Description: Laminar and turbulent perfect-gas air, Navier-Stokes computations have been performed for a proposed Mars Smart Lander entry vehicle at Mach 6 over a free stream Reynolds number range of 6.9 x 10(exp 6/m to 2.4 x 10(exp 7)m(2.1 x 10(exp 6)/ft to 7.3 x 10(exp 6)ft) for angles-of-attack of 0-deg, 11-deg, 16-deg, and 20-deg, and comparisons were made to wind tunnel heating data obtained at the same conditions. Boundary layer edge properties were extracted from the solutions and used to correlate experimental data on the effects of heat-shield penetrations (bolt-holes where the entry vehicle would be attached to the propulsion module during transit to Mars) on boundary-layer transition. A non-equilibrium Martian-atmosphere computation was performed for the peak heating point on the entry trajectory in order to determine if the penetrations would produce boundary-layer transition by using this correlation.

Description: An experimental study of the application of discrete-time, linear quadratic control design methods to the cavity tone problem is described. State space models of the dynamics from a synthetic jet actuator at the leading edge of the cavity to two pressure sensors in the cavity were computed from experimental data. Variations in model order, control order, control bandwidth, and properties of a Kalman state estimator were studied. Feedback control reduced the levels of multiple cavity tones at Mach 0.275, 0.35, and 0.45. Closed loop performance was often limited by excitation of sidebands of cavity tones, and creation of new tones in the spectrum. State space models were useful for explaining some of these limitations, but were not able to account for non-linear dynamics, such as interactions between tones at different frequencies.

Description: Stochastic versions of stability equations are developed in order to develop integrated models of transition and turbulence and to understand the effects of uncertain initial conditions on disturbance growth. Stochastic forms of the resonant triad equations, a high Reynolds number asymptotic theory, and the parabolized stability equations are developed.

Description: As a result of CFD analysis and remachining of the nozzle, the flow quality of the Mach 10 Hypersonic Wind Tunnel at NAL Chofu, Japan was improved. The subsequent test results validated the CFD analytical predictions by NASA and MHL.

Description: Numerical simulations of a single low-profile vortex generator vane, which is only a small fraction of the boundary-layer thickness, and a vortex generating jet have been performed for flows over a flat plate. The numerical simulations were computed by solving the steady-state solution to the Reynolds-averaged Navier-Stokes equations. The vortex generating vane results were evaluated by comparing the strength and trajectory of the streamwise vortex to experimental particle image velocimetry measurements. From the numerical simulations of the vane case, it was observed that the Shear-Stress Transport (SST) turbulence model resulted in a better prediction of the streamwise peak vorticity and trajectory when compared to the Spalart-Allmaras (SA) turbulence model. It is shown in this investigation that the estimation of the turbulent eddy viscosity near the vortex core, for both the vane and jet simulations, was higher for the SA model when compared to the SST model. Even though the numerical simulations of the vortex generating vane were able to predict the trajectory of the stream-wise vortex, the initial magnitude and decay of the peak streamwise vorticity were significantly under predicted. A comparison of the positive circulation associated with the streamwise vortex showed that while the numerical simulations produced a more diffused vortex, the vortex strength compared very well to the experimental observations. A grid resolution study for the vortex generating vane was also performed showing that the diffusion of the vortex was not a result of insufficient grid resolution. Comparisons were also made between a fully modeled trapezoidal vane with finite thickness to a simply modeled rectangular thin vane. The comparisons showed that the simply modeled rectangular vane produced a streamwise vortex which had a strength and trajectory very similar to the fully modeled trapezoidal vane.

Description: Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes with incremental studies was initiated. Here we further refine the analysis on, and improve the understanding of the adaptive numerical dissipation control strategy. Basically, the development of these schemes focuses on high order nondissipative schemes and takes advantage of the progress that has been made for the last 30 years in numerical methods for conservation laws, such as techniques for imposing boundary conditions, techniques for stability at shock waves, and techniques for stable and accurate long-time integration. We concentrate on high order centered spatial discretizations and a fourth-order Runge-Kutta temporal discretizations as the base scheme. Near the bound-aries, the base scheme has stable boundary difference operators. To further enhance stability, the split form of the inviscid flux derivatives is frequently used for smooth flow problems. To enhance nonlinear stability, linear high order numerical dissipations are employed away from discontinuities, and nonlinear filters are employed after each time step in order to suppress spurious oscillations near discontinuities to minimize the smearing of turbulent fluctuations. Although these schemes are built from many components, each of which is well-known, it is not entirely obvious how the different components be best connected. For example, the nonlinear filter could instead have been built into the spatial discretization, so that it would have been activated at each stage in the Runge-Kutta time stepping. We could think of a mechanism that activates the split form of the equations only at some parts of the domain. Another issue is how to define good sensors for determining in which parts of the computational domain a certain feature should be filtered by the appropriate numerical dissipation. For the present study we employ a wavelet technique introduced in as sensors. Here, the method is briefly described with selected numerical experiments.

Description: For the last several years, Glenn-HT, a three-dimensional (3D) Computational Fluid Dynamics (CFD) computer code for the analysis of gas turbine flow and convective heat transfer has been evolving at the NASA Glenn Research Center. The code is unique in the ability to give a highly detailed representation of the flow field very close to solid surfaces in order to get accurate representation of fluid beat transfer and viscous shear stresses. The code has been validated and used extensively for both internal cooling passage flow and for hot gas path flows, including detailed film cooling calculations and complex tip clearance gap flow and heat transfer. In its current form, this code has a multiblock grid capability and has been validated for a number of turbine configurations. The code has been developed and used primarily as a research tool, but it can be useful for detailed design analysis. In this presentation, the code is described and examples of its validation and use for complex flow calculations are presented, emphasizing the applicability to turbomachinery.

Description: The Low Reynolds number version of the Stress-omega model and the two equation k-omega model of Wilcox were used for the calculation of turbulent heat transfer in a 180 degree turn simulating an internal coolant passage. The Stress-omega model was chosen for its robustness. The turbulent thermal fluxes were calculated by modifying and using the Generalized Gradient Diffusion Hypothesis. The results showed that using this Reynolds Stress model allowed better prediction of heat transfer compared to the k-omega two equation model. This improvement however required a finer grid and commensurately more CPU time.

Description: Comparison of turbulence data taken in three separate flow nozzles, two with mixing enhancement features on their core nozzle, shows how the mixing enhancement features modify turbulence to reduce jet noise. The three nozzles measured were the baseline axisymmetric nozzle 3BB, the alternating chevron nozzle, 3A12B, with 6-fold symmetry, and the flipper tab nozzle 3T24B also with 6-fold symmetry. The data presented show the differences in turbulence characteristics produced by the geometric differences in the nozzles, with emphasis on those characteristics of interest in jet noise. Among the significant findings: the enhanced mixing devices reduce turbulence in the jet mixing region while increasing it in the fan/core shear layer, the ratios of turbulence components are significantly altered by the mixing devices, and the integral lengthscales do not conform to any turbulence model yet proposed. These findings should provide guidance for modeling the statistical properties of turbulence to improve jet noise prediction.

Description: The Sixth Microgravity Fluid Physics and Transport Phenomena Conference provides the scientific community the opportunity to view the current scope of the Microgravity Fluid Physics and Transport Phenomena Program, current research opportunities, and plans for the near future. The conference focuses not only on fundamental research but also on applications of this knowledge towards enabling future space exploration missions. A whole session dedicated to biological fluid physics shows increased emphasis that the program has placed on interdisciplinary research. The conference includes invited plenary talks, technical paper presentations, poster presentations, and exhibits. This TM is a compilation of abstracts of the papers and the posters presented at the conference. Web-based proceedings, including the charts used by the presenters, will be posted on the web shortly after the conference.

Description: Three examples of response surface modeling with CFD are presented for combined cycle propulsion components. The examples include a mixed-compression-inlet during hypersonic flight, a hydrogen-fueled scramjet combustor during hypersonic flight, and a ducted-rocket nozzle during all-rocket flight. Three different experimental strategies were examined, including full factorial, fractionated central-composite, and D-optimal with embedded Plackett-Burman designs. The response variables have been confined to integral data extracted from multidimensional CFD results. Careful attention to uncertainty assessment and modeling bias has been addressed. The importance of automating experimental setup and effectively communicating statistical results are emphasized.

Description: This paper discusses the effect of the pressure drop on the operating temperature in a loop heat pipe (LHP). Because the evaporator and the compensation chamber (CC) both contain two-phase fluid, a thermodynamic constraint exists between the temperature difference and the pressure drop for these two components. As the pressure drop increases, so will the temperature difference. The temperature difference in turn causes an increase of the heat leak from the evaporator to the CC, resulting in a higher CC temperature. Furthermore, the heat leak strongly depends on the vapor void fraction inside the evaporator core. Tests were conducted by installing a valve on the vapor line so as to vary the pressure drop, and by charging the LHP with various amounts of fluid. Test results verify that the LHP operating temperature increases with an increasing differential pressure, and the temperature increase is a strong function of the fluid inventory in the loop.

Description: Capillary dominated fluid dynamics will be examined in a reduced-gravity environment onboard the KC-135; in particular, the behavior of the lower portion of the meniscus in triangular tank geometries. Seven clear acrylic tanks were constructed to view seven angles of the four geometries. Silicon oil with two different viscosities, 2cs and 5cs silicon oil, were used on different days of the flight. Six tanks and one control tank are filled with a certain viscosity fluid for each flight day. During each parabola, three tanks are tested at time. The experimental tanks are exchanged between parabola sets on the KC-135. The 60deg -60deg -60deg control tank is viewed throughout the flight. To gather data, two digital video cameras and one digital still camera are placed perpendicular the viewing surface. To provide a greater contrast in the meniscus, an EL backlighting sheet was used to backlight the tanks. These images and video are then digitized, passed through NASA's mini-tracker software, and compared to a theory published my M. M. Weislogel, "Fluid Interface Phenomena in a Low-Gravity Environment: Recent Results from Drop Tower Experimentation." By focusing on a lower portion of the meniscus and using longer periods of reduced gravity, this experiment may confirm that a stationary point exists on the fluid surface. This information will enable better designing of propellant management devices, especially satellite propellant refilling and gas venting. Also, biological and material processing systems in reduced gravity environments will benefit from this data.

Description: A thermal analysis of the compressible CO2 flow for the Portable Fire Extinguisher (PFE) system has been performed. A SINDA/FLUINT model has been developed for this analysis. The model includes the PFE tank and the TeSS nozzle, and both have initial temperature of 72 F. In order to investigate the thermal effect on the nozzle due to discharging CO2, the PFE TeSS nozzle pipe has been divided into three segments. This model also includes heat transfer predictions for PFE tank inner and outer wall surfaces. The simulation results show that the CO2 discharge rates have fulfilled the minimum flow requirements that the PFE system discharges 3.0 Ibm CO2 in 10 seconds and 5.5 Ibm of CO2 in 45 seconds during its operation. At 45 seconds, the PFE tank wall temperature is 63 F, and the TeSS nozzle cover wall temperatures for the three segments are 47 F, 53 F and 37 F, respectively. Thermal insulation for personal protection is used for the first two segments of the TeSS nozzle. The simulation results also indicate that at 50 seconds, the remaining CO2 in the tank may be near the triple point (gas, liquid and solid) state and, therefore, restricts the flow.

Description: Experimental results from a study of the effects of passing wakes upon laminar-to-turbulent transition in a low-pressure turbine passage are presented. The test section geometry is designed to simulate the effects of unsteady wakes resulting from rotor-stator interaction upon laminar-to-turbulent transition in turbine blade boundary layers and separated flow regions over suction surfaces. Single-wire, thermal anemometry techniques were used to measure time- resolved and phase-averaged, wall-normal profiles of velocity, turbulence intensity, and intermittency at multiple streamwise locations over the turbine airfoil suction surface. These data are compared to steady state, wake-free data collected in the same geometry to identify the effects of wakes upon laminar-to-turbulent transition. Results are presented for flows with a Reynolds number based on suction surface length and exit velocity of 50,000 and an approach flow turbulence intensity of 2.5 percent. From these data, the effects of passing wakes and associated increased turbulence levels and varying pressure gradients on transition and separation in the near-wall flow are presented. The results show that the wakes affect transition both by virtue of their difference in turbulence level from that of the free-stream but also by virtue of their velocity deficit relative to the freestream velocity, and the concomitant change in angle of attack and temporal pressure gradients. The results of this study seem to support the theory that bypass transition is a response of the near-wall viscous layer to pressure fluctuations imposed upon it from the free-stream flow. The data also show a significant lag between when the wake is present over the surface and when transition begins. The accompanying CD-ROM includes tabulated data, animations, higher resolution plots, and an electronic copy of this report.

Description: This paper reports on ground based work conducted to support the Spaceflight Definition project SHIVA (Spaceflight Holography Investigation in a Virtual Apparatus). SHIVA will advance our understanding of the movement of a particle in a fluid. Gravity usually dominates the equations of motion, but in microgravity as well as on earth other terms can become important. Before two members of our team found an analytical solution of the equations, numerical methods and/or neglecting terms were required. The general solution predicts that the usually neglected history term becomes important when the characteristic viscous time is in the same order as the vibration period and peaks when the two times are equal. In this case three force terms, the Stokes drag, the added mass, and the history drag must all be included in predicting particle movement. We also developed diagnostic recording methods using holography to save all of the particle field data, allowing the experiment to essentially be transferred from space back to earth in what we call the "virtual apparatus". Using state-of-the-art methods in holography we will quantify the three-dimensional motion of sets of particles, allowing us to test and apply the new analytical solutions. The motion of particles up to 4 mm in diameter in a fluid that oscillates at frequencies up to 100 Hz with amplitudes up to 200 microns is being examined. Ground studies to support the flight development program have employed various schemes to simulate microgravity. One of the most reliable and meaningful methods uses spheres tethered to a fine hair suspended in the fluid. We have also investigated particles with nearly neutral buoyancy. Recordings are made at the peak amplitudes of vibration of the cell providing a measure of the ratio of fluid to particle amplitude. The experiment requires precise location of the particle at the time of recording. The hologram of the particle provides microscopic images of the particle that are used for finding the position with an accuracy of a few microns. To make the experiment more versatile, the spaceflight system will record holograms both on film and electronically. The electronic holograms can be downlinked providing real time data. Results of the ground experiments, the flight experiment design, and data analysis procedures are reported.

Description: When an alternating magnetic field interacts with induced eddy currents in a conducting body, there will be a repulsive force between the body and the driving coil system generating the field. This repulsive force is the basis of electromagnetic levitation, which allows containerless processing of different materials. The eddy currents in the conducting body also generate Joule heating. Axial rotation of electromagnetically levitated objects is a common observation in levitation systems and often an undesirable side effect of such experiments on 1-g and -g. There have been recent efforts to use magnetic damping and suppress this tendency of body rotation. The first report of rotation in EML drops was attributed to a slight asymmetry of the shape and location of the levitation coils could change the axis and speed of rotation. Other theories of sample rotation include a frequency difference in the traveling electromagnetic waves and a phase difference in two different applied fields of the same frequency. All of these different mechanisms share the following characteristics: the torque is small, constant for constant field strength, and very weakly dependent on the sample's temperature and phase (solid or liquid). During experiments on the MSL-1 (First Microgravity Science Laboratory) mission of the Space Shuttle (STS-83 and STS-94, April and July 1997), a droplet of palladium-silicon alloy was electromagnetically levitated for viscosity measurements. For the non-deforming droplet, the resultant MHD flow inside the drop is inferred from motion of impurities on the surface. These observations indicate formation of a pair of co-rotating toroidal flow structures inside the spheroidal levitated drop that undergo secondary flow instabilities. As rise in the fluid temperature rises, the viscosity falls and the internal flow accelerates and becomes oscillatory; and beyond a point in the experiments, the surface impurities exhibit non-coherent chaotic motion signifying emergence of turbulence inside the drop. In this work, a background of these set of observations will be given followed by a presentation of our results on the digital particle tracking analysis that has been performed on a number of available videos. The analysis indicates that the levitated drop attains a constant rotational speed during the melting phase and formation of the co-rotating axi-symmetric laminar toroidal structures. However, the rate of axial rotation increases dramatically during the deformation of the toroidal structures anti their breakup into chaotic entities. This new data suggests an interaction between the flow inside the levitated molten drop and the driving coils in the experiments. Possible mechanisms for this interaction will be reviewed. The data will also be used to make an assessment of existing theories on droplet rotation.

Description: Oscillating droplets are of interest in a number of disciplines. A practical application is the oscillating drop method, which is a technique for measuring surface tension and viscosity of liquid metals. It is especially suited to undercooled and highly reactive metals, because it is usually performed by electromagnetic levitation. The natural oscillation frequency a f the droplets is related to the surface tension of the material, and the decay of oscillations is related to its viscosity. The fluid flow inside the droplet must be laminar in order for this technique to yield good results. Because no experimental method has yet been developed to visualize flow in electromagnetically-levitated oscillating metal droplets, mathematical modeling is required to determine whether or not turbulence occurs. Three mathematical models of the flow: (1) assuming laminar conditions, (2) using the turbulence model, and (3) using the RNG turbulence model, respectively, are compared and contrasted to determine the physical characteristics of the flow. It is concluded that the RNG model is the best suited for describing this problem. The goal of the presented work was to characterize internal flow in an oscillating droplet of liquid metal, and to verify the accuracy of the characterization by comparing calculated surface tension and viscosity values to experimentally-obtained results.

Description: During experiments on the MSL-1 (First Microgravity Science Laboratory) mission of the Space Shuttle (STS-83 and STS-94, April and July 1997), a droplet of palladium-silicon alloy was electromagnetically levitated for viscosity measurements. For the non-deforming droplet, the resultant MHD flow inside the drop can be inferred from motion of impurity particulates on the surface. In the experiments, subsequent to melting, Joule heating produces a continuous reduction of viscosity of the fluid resulting in an acceleration of the flow with time. These observations indicate formation of a pair of co-rotating toroidal flow structures inside the spheroidal drop that undergo flow instabilities. As the fluid temperature rises, the amplitude of the secondary flow increases; and beyond a point, the tracers exhibit non-coherent chaotic motion signifying emergence of turbulence inside the drop. Assuming that the observed laminar-turbulent transition is shear-layer type, the internal structure of the toroidal loops is used to develop a semi- empirical correlation for the onset of turbulence. Our calculations indicate that the suggested correlation is in modest agreement with the experimental data, with the transition occurring at a Reynolds number of 600.

Description: I will discuss recent experiments from my lab, which use surface templates to induce ordered colloidal structures. Particle assembly driven by entropic depletion, fluid convection, and sedimentation will be described. Confocal microscopy was used to visualize most of these samples.

Description: The widespread use of electro-hydrodynamic devices and processes emphasizes a critical need for developing a comprehensive predictive theory capable of improving our fundamental understanding of the behavior of a suspension subject to an AC electric field and shear, and of facilitating the design and optimization of such devices. The currently favored approach to the qualitative interpretation of the AC field driven manipulation of suspensions is based on a model which considers only the force exerted on a single particle by an external field and neglects the field-induced and hydrodynamic interparticle interactions both being inversely proportional to the interparticle distance raised to the power three. On the other hand, the purpose of the field-induced separation is to concentrate particles in certain regions of a device. This clearly raises the fundamental question regarding the extent to which we can neglect these slow decaying electrical and hydrodynamic collective interactions and rely on the predictions of a single-particle model. Another important issue that still remains open is how to characterize the polarization of a particle exposed to a strong electric field. The presentation will address both these questions. Experiments were conducted in a parallel-plate channel in which a 10(exp -3) (v/v) suspension of heavy, positively polarized Al2O3 spheres was exposed to an AC field under conditions such that the field lines were arranged in the channel cross-section perpendicular to the streamlines of the main flow. To reduce the effects of the gravitational settling of the particles, the channel was slowly rotated (4 rpm) around a horizontal axis. Following the application of a high-gradient strong AC field (approx. kV/mm), the particles were found to move towards both the high-voltage (HV) and grounded (GR) electrodes and to form 'bristles' along their edges.

Description: The development of a unifying framework among direct numerical simulations, large-eddy simulations, and statistically averaged formulations of the Navier-Stokes equations, is of current interest. Toward that goal, the properties of the residual (subgrid-scale) stress of the temporally filtered Navier-Stokes equations are carefully examined. Causal time-domain filters, parameterized by a temporal filter width 0 less than Delta less than infinity, are considered. For several reasons, the differential forms of such filters are preferred to their corresponding integral forms; among these, storage requirements for differential forms are typically much less than for integral forms and, for some filters, are independent of Delta. The behavior of the residual stress in the limits of both vanishing and in infinite filter widths is examined. It is shown analytically that, in the limit Delta to 0, the residual stress vanishes, in which case the Navier-Stokes equations are recovered from the temporally filtered equations. Alternately, in the limit Delta to infinity, the residual stress is equivalent to the long-time averaged stress, and the Reynolds-averaged Navier-Stokes equations are recovered from the temporally filtered equations. The predicted behavior at the asymptotic limits of filter width is further validated by numerical simulations of the temporally filtered forced, viscous Burger's equation. Finally, finite filter widths are also considered, and a priori analyses of temporal similarity and temporal approximate deconvolution models of the residual stress are conducted.

Description: A computer program performs calculations for refinement or coarsening of computational grids of the type called "structured" (signifying that they are geometrically regular and/or are specified by relatively simple algebraic expressions). This program is designed to facilitate analysis of the numerical effects of changing structured grids utilized in computational fluid dynamics (CFD) software. Unlike prior grid-refinement and -coarsening programs, this program is not limited to doubling or halving: the user can specify any refinement or coarsening ratio, which can have a noninteger value. In addition to this ratio, the program accepts, as input, a grid file and the associated restart file, which is basically a file containing the most recent iteration of flow-field variables computed on the grid. The program then refines or coarsens the grid as specified, while maintaining the geometry and the stretching characteristics of the original grid. The program can interpolate from the input restart file to create a restart file for the refined or coarsened grid. The program provides a graphical user interface that facilitates the entry of input data for the grid-generation and restart-interpolation routines.

Description: The present series of three consecutive papers develops a general theory for linear and finite solid viscoelasticity. Because the most important object for nonlinear studies are rubber-like materials, the general approach is specified in a form convenient for solving problems important for many industries that involve rubber-like materials. General linear and nonlinear theories for non-isothermal deformations of viscoelastic solids are developed based on the quasi-linear approach of non-equilibrium thermodynamics. In this, the first paper of the series, we analyze non-isothermal linear viscoelasticity, which is applicable in a range of small strains not only to all synthetic polymers and bio-polymers but also to some non-polymeric materials. Although the linear case seems to be well developed, there still are some reasons to implement a thermodynamic derivation of constitutive equations for solid-like, non-isothermal, linear viscoelasticity. The most important is the thermodynamic modeling of thermo-rheological complexity , i.e. different temperature dependences of relaxation parameters in various parts of relaxation spectrum. A special structure of interaction matrices is established for different physical mechanisms contributed to the normal relaxation modes. This structure seems to be in accord with observations, and creates a simple mathematical framework for both continuum and molecular theories of the thermo-rheological complex relaxation phenomena. Finally, a unified approach is briefly discussed that, in principle, allows combining both the long time (discrete) and short time (continuous) descriptions of relaxation behaviors for polymers in the rubbery and glassy regions.

Description: A discontinuous Galerkin (DG) method is formulated, implemented, and tested for simulation of compressible turbulent flows. The method is applied to turbulent channel flow at low Reynolds number, where it is found to successfully predict low-order statistics with fewer degrees of freedom than traditional numerical methods. This reduction is achieved by utilizing local hp-refinement such that the computational grid is refined simultaneously in all three spatial coordinates with decreasing distance from the wall. Another advantage of DG is that Dirichlet boundary conditions can be enforced weakly through integrals of the numerical fluxes. Both for a model advection-diffusion problem and for turbulent channel flow, weak enforcement of wall boundaries is found to improve results at low resolution. Such weak boundary conditions may play a pivotal role in wall modeling for large-eddy simulation.

Description: Experimental results from a study of the effects of passing wakes upon laminar-to-turbulent transition in a low-pressure turbine passage are presented. The test section geometry is designed to simulate the effects of unsteady wakes resulting from rotor-stator interaction upon laminar-to-turbulent transition in turbine blade boundary layers and separated flow regions over suction surfaces. Single-wire, thermal anemometry techniques were used to measure time-resolved and phase-averaged, wall-normal profiles of velocity, turbulence intensity, and intermittency at multiple streamwise locations over the turbine airfoil suction surface. These data are compared to steady state, wake-free data collected in the same geometry to identify the effects of wakes upon laminar-to-turbulent transition. Results are presented for flows with a Reynolds number based on suction surface length and exit velocity of 50,000 and an approach flow turbulence intensity of 2.5 percent. From these data, the effects of passing wakes and associated increased turbulence levels and varying pressure gradients on transition and separation in the near-wall flow are presented. The results show that the wakes affect transition both by virtue of their difference in turbulence level from that of the free-stream but also by virtue of their velocity deficit relative to the freestream velocity, and the concomitant change in angle of attack and temporal pressure gradients. The results of this study seem to support the theory that bypass transition is a response of the near-wall viscous layer to pressure fluctuations imposed upon it from the free-stream flow. The data also show a significant lag between when the wake is present over the surface and when transition begins. The accompanying CD-ROM includes tabulated data, animations, higher resolution plots, and an electronic copy of this report.

Description: Whipped cream stays in place even when turned upside down. Yet it readily flows through the nozzle of a spray can to reach the dessert plate. This demonstrates the phenomenon of shear thinning that is important to many industrial and physical processes. Paints, film emulsions, and other complex solutions that are highly viscous under normal conditions but become thin and flow easily under shear forces. A simple fluid, such as water, does not exhibit shear thinning under normal conditions. Very close to the liquid-vapor critical point, where the distinction between liquid and vapor disappears, the fluid becomes more complex and is predicted to display shear thinning. At the critical point, xenon atoms interact over long distances in a classical model of cooperative phenomena. Physicists rely on this system to learn how long-range order arises. The Critical Viscosity of Xenon Experiment (CVX-2) will measure the viscous behavior of xenon, a heavy inert gas used in flash lamps and ion rocket engines, at its critical point. Although it does not easily combine with other chemicals, its viscosity at the critical point can be used as a model for a range of fluids. Viscosity originates from the interactions of individual molecules. It is so complicated that, except for the simplest gas, it cannot be calculated accurately from theory. Tests with critical fluids can provide key data, but are limited on Earth because critical fluids are highly compressed by gravity. CVX-2 employs a tiny metal screen vibrating between two electrodes in a bath of critical xenon. The vibrations and how they dampen are used to measure viscosity. CVX flew on STS-85 (1997), where it revealed that, close to the critical point, the xenon is partly elastic: it can 'stretch' as well as flow. For STS-107, the hardware has been enhanced to determine if critical xenon is a shear-thinning fluid.

Description: We have synthesized two deuterated boomerang-shaped liquid crystals based on 2,5-bis(4-hydroxyphenyl)-1,3,4-oxadiazole (ODBP). Deuterium was introduced in the rigid 2,5-diphenyl-1,3,4-oxadiazole core and in the aromatic ring of the terminal 4-dodecyloxyphenyl moiety using standard acid catalyzed deuterium exchange conditions. Both compounds, ([4,4'(1,3,4-oxadiazole-2,5-diyl-d4)] di-4-dodecyloxybenzoate: ODBP-d4-Ph-O-C12) and ([4,4'(1,3,4-oxadiazole-2,5-diyl)] di-4-dodecyloxy-benzoate-d4; ODBP-Ph-d4-O-C12) were investigated by nuclear magnetic resonance, optical microscopy and differential scanning calorimetry. The optical textures and thermal behavior of both compounds were found to be identical to the non-deuterated analog [4,4(1,3,4-oxadiazole-2,5-diyl)] di-4-dodecyloxybenzoate (ODBP-Ph-O-C12) which we reported earlier. These compounds exhibit behavior indicative of a biaxial nematic liquid crystal phase, which we hope to confirm using deuterium NMR spectroscopy in the next phase of this study.

Description: This paper presents an extension of a numerical algorithm for network flow analysis code to perform multi-dimensional flow calculation. The one dimensional momentum equation in network flow analysis code has been extended to include momentum transport due to shear stress and transverse component of velocity. Both laminar and turbulent flows are considered. Turbulence is represented by Prandtl's mixing length hypothesis. Three classical examples (Poiseuille flow, Couette flow and shear driven flow in a rectangular cavity) are presented as benchmark for the verification of the numerical scheme.