ALBERT

Description: Microgravity research at NASA has been an undertaking that has included both science and commercial approaches since the late 80s and early 90s. The Fluid Physics and Transport Phenomena community has been developed, through NASA's science grants, into a valuable base of expertise in microgravity science. This was achieved through both ground and flight scientific research. Commercial microgravity research has been primarily promoted thorough NASA sponsored Centers for Space Commercialization which develop cost sharing partnerships with industry. As an example, the Center for Advanced Microgravity Materials Processing (CAMMP)at Northeastern University has been working with cost sharing industry partners in developing Zeolites and zeo-type materials as an efficient storage medium for hydrogen fuel. Greater commercial interest is emerging. The U.S. Congress has passed the Commercial Space Act of 1998 to encourage the development of a commercial space industry in the United States. The Act has provisions for the commercialization of the International Space Station (ISS). Increased efforts have been made by NASA to enable industrial ventures on-board the ISS. A Web site has been established at http://commercial/nasa/gov which includes two important special announcements. One is an open request for entrepreneurial offers related to the commercial development and use of the ISS. The second is a price structure and schedule for U.S. resources and accommodations. The purpose of the presentation is to make the Fluid Physics and Transport Phenomena community, which understands the importance of microgravity experimentation, aware of important aspects of ISS commercial development. It is a desire that this awareness will be translated into a recognition of Fluid Physics and Transport Phenomena application opportunities coordinated through the broad contacts of this community with industry.

Description: This overview presents in viewgraph form, the NASA Program organization regarding fluid physics, physical sciences research in space and the connection to biology, the dual thrust of the fluid physics program, and the immediate and future plans of the physical science research division.

Description: On a twelve-month voyage to Mars, one astronaut will require at least two tons of potable water and two tons of pure oxygen. Efficient, reliable fluid reclamation is therefore necessary for manned space exploration. Space habitats require a compact, flexible, and robust apparatus capable of solid-fluid mechanical separation over a wide range of fluid and particle densities and particle sizes. In space, centrifugal filtration, where particles suspended in fluid are captured by rotating fixed-fiber mat filters, is a logical candidate for mechanical separation. Non-colloidal particles are deposited on the fibers due to inertial impaction or direct interception. Since rotation rates are easily adjustable, inertial effects are the most practical way to control separation rates for a wide variety of multiphase mixtures in variable gravity environments. Understanding how fluid inertia and differential fluid-particle inertia, characterized by the Reynolds and Stokes numbers, respectively, affect deposition is critical in optimizing filtration in a microgravity environment. This work will develop non-intrusive optical diagnostic techniques for directly visualizing where and when non-colloidal particles deposit upon, or contact, solid surfaces: 'particle proximity sensors'. To model particle deposition upon a single filter fiber, these sensors will be used in ground-based experiments to study particle dynamics as in the vicinity of a large (compared with the particles) cylinder in a simply sheared (i.e., linearly-varying, zero-mean velocity profile) neutrally-buoyant, refractive-index matched solid-liquid suspension.

Description: The present work analyses the dynamics of a suspension of heavy particles in shear flow. The magnitude of the particle inertia is given by the Stokes number St = m(gamma/6(pi)a, which is the ratio of the viscous relaxation time of a particle tau(sub p) = m=6pi(eta)a to the flow time gamma(sup -1). Here, m is the mass of the particle, a is its size, eta is the viscosity of the suspending fluid and gamma is the shear rate. The ratio of the Stokes number to the Reynolds number, Re = (rho)f(gamma)a(exp 2)/eta, is the density ratio rho(sub p)/rho(sub f). Of interest is to understand the separate roles of particle (St) and fluid (Re) inertia in the dynamics of suspensions. In this study we focus on heavy particles, rho(sub p)/rho(sub f) much greater than 1, for which the Stokes number is finite, but the Reynolds number is sufficiently small for inertial forces in the fluid to be neglected; thus, the fluid motion is governed by the Stokes equations. On the other hand, the probability density governing the statistics of the suspended particles satisfies a Fokker-Planck equation that accounts for both configuration and momentum coordinates, the latter being essential for finite St. The solution of the Fokker-Planck equation is obtained to O(St) via a Chapman-Enskog type-procedure, and the conditional velocity distribution so obtained is used to derive a configuration-space Smoluchowski equation with inertial corrections. The inertial effects are responsible for asymmetry in the relative trajectories of two spheres in shear flow, in contrast to the well known symmetric structure in the absence of inertia. Finite St open trajectories in the plane of shear suffer a downward lateral displacement resulting from the inability of a particle of finite mass to follow the curvature of the zero-Stokes-number pathlines. In addition to the induced asymmetry, the O(St) inertial perturbation dramatically alters the nature of the near-field trajectories. The stable closed orbits (for St = 0) in the plane of shear now spiral in, approaching particle-particle contact in the limit. All trajectories starting from an initial offset of O(St(sup 1/2) or less (which remain open for St = 0) also spiral in. The asymmetry of the trajectories leads to a non-Newtonian rheology and diffusive behavior. The latter because a given particle (moving along a finite St open trajectory) suffers a net displacement in the transverse direction after a single interaction. A sequence of such uncorrelated displacements leads to the particle executing a random walk. The inertial diffusivity tensor is anisotropic on account of differing strengths of interaction in the gradient and vorticity directions. Since the entire region (constituting an in finite area) of closed orbits in the plane of shear spirals onto contact for #finite St, the latter represents a singular surface for the pair-distribution function. The exact form of the pair-distribution function at contact is still, however, indeterminate in the absence of non-hydrodynamic effects. It should also be noted that finite St non-rectilinear flows do not support a spatially uniform number density owing to the cross-streamline inertial migration of particles.

Description: Experiments as well as accompanying simulations are described that serve in preparation of a space flight experiment to study the dynamics of miscible interfaces. The investigation specifically addresses the importance of both nonsolenoidal effects as well as nonconventional Korteweg stresses in flows that give rise to steep but finite concentration gradients. The investigation focuses on the flow in which a less viscous fluid displaces one of higher viscosity and different density within a narrow capillary tube. The fluids are miscible in all proportions. An intruding finger forms that occupies a fraction of the total tube diameter. Depending on the flow conditions, as expressed by the Peclet number, a dimensionless viscosity ratio, and a gravity parameter, this fraction can vary between approximately 0.9 and 0.2. For large Pe values, a quasi-steady finger forms, which persists for a time of O(Pe) before it starts to decay, and Poiseuille flow and Taylor dispersion are approached asymptotically. Depending on the specific flow conditions, we observe a variety of topologically different streamline patterns, among them some that leak fluid from the finger tip. For small Pe values, the flow decays from the start and asymptotically reaches Taylor dispersion after a time of O(Pe). Comparisons between experiments and numerical simulations based on the 'conventional' assumption of solenoidal velocity fields and without Korteweg stresses yield poor agreement as far as the Pe value is concerned that distinguishes these two regimes. As one possibility, we attribute this lack of agreement to the disregard of these terms. An attempt is made to use scaling arguments in order to evaluate the importance of the Korteweg stresses and of the assumption of solenoidality. While these effects should be strongest in absolute terms when steep concentration fronts exist, i.e., at large Pe, they may be relatively most important at lower values of Pe. We subsequently compare these conventional simulations to more complete simulations that account for nonvanishing divergence as well as Korteweg stresses. While the exact value of the relevant stress coefficients are not known, ballpark numbers do exist, and their use in the simulations indicates that these stresses may indeed be important. We plan to evaluate these issues in detail by means of comparing a space experiment with corresponding simulations, in order to extract more accurate Korteweg stress coefficients, and to confirm or deny the importance of such stresses.

Description: The vibratory conveyor, routinely employed for normal-gravity transport of granular materials, usually consists of a continuous open trough vibrated sinusoidally to induce axial movement of a granular material. Motivated in part by a hypothetical application in zero gravity, we propose a novel modification of the vibratory conveyor based on a closed 2d trough operating in a "slide-conveying" mode, with the granular mass remaining permanently in contact with the trough walls. We present a detailed analysis of the mechanics of transport, based on a rigid-slab model for the granular mass with frictional (Coulomb) slip at the upper and lower walls. The form of the vibration cycle plays a crucial role, and the optimal conveying cycle is not the commonly assumed rectilinear sinusoidal motion. The conveying efficiency for the novel slide conveyor will be presented for several simple vibration cycles, including one believed to represent the theoretical optimum.

Description: The research carried out in the Heat Transfer Laboratory of the Johns Hopkins University was motivated by previous studies indicating that in terrestrial applications nucleate boiling heat transfer can be increased by a factor of 50 when compared to values obtained for the same system without electric fields. Imposing an external electric field holds the promise to improve pool boiling heat transfer in low gravity, since a phase separation force other than gravity is introduced. The influence of electric fields on bubble formation has been investigated both experimentally and theoretically.

Description: The existing RocketWeb(TradeMark) Internet Analysis System (httr)://www.iohnsonrockets.com/rocketweb) provides an integrated set of advanced analysis tools that can be securely accessed over the Internet. Since these tools consist of both batch and interactive analysis codes, the system includes convenient methods for creating input files and evaluating the resulting data. The RocketWeb(TradeMark) system also contains many features that permit data sharing which, when further developed, will facilitate real-time, geographically diverse, collaborative engineering within a designated work group. Adding work group management functionality while simultaneously extending and integrating the system's set of design and analysis tools will create a system providing rigorous, controlled design development, reducing design cycle time and cost.

Description: A methodology for efficiently including three-dimensional effects on aircraft plume signature is presented. First, exploratory work on the use of passive mixing enhancement devices, namely chevrons and tabs, in IR signature reduction for external turbofan plumes is demonstrated numerically and experimentally. Such small attachments, when properly designed, cause an otherwise axisymmetric plume to have significant 3D structures, affecting signature prediction. Second, an approach for including non-axisymmetric and installation effects in plume signature prediction is discussed using unstructured methodology. Unstructured flow solvers, using advanced turbulence modeling and plume thermochemistry, facilitate the modeling of aircraft effects on plume structure that previously have been neglected due to gridding complexities. The capabilities of the CRUNCH unstructured Navier-Stokes solver for plume modeling is demonstrated for a passively mixed turbofan nozzle, a generic fighter nozzle, and a complete aircraft.

Description: The report documents the recent effort to enhance a transient linear heat transfer code so as to solve nonlinear problems. The linear heat transfer code was originally developed by Dr. Kim Bey of NASA Largely and called the Structure-Compatible Heat Transfer (SCHT) code. The report includes four parts. The first part outlines the formulation of the heat transfer problem of concern. The second and the third parts give detailed procedures to construct the nonlinear finite element equations and the required Jacobian matrices for the nonlinear iterative method, Newton-Raphson method. The final part summarizes the results of the numerical experiments on the newly enhanced SCHT code.

Description: Recent technology developments and a systems engineering design approach have led to the development of a practical battery-free solar refrigerator as a spin-off of NASA's aerospace refrigeration research. Off-grid refrigeration is a good application of solar photovoltaic (PV) power if thermal storage is incorporated and a direct connection is made between the cooling system and the PV panel. This was accomplished by integrating water as a phase-change material into a well insulated refrigerator cabinet and by developing a microprocessor based control system that allows direct connection of a PV panel to a variable speed compressor. This second innovation also allowed peak power-point tracking from the PV panel and elimination of batteries from the system. First a laboratory unit was developed to prove the concept and then a commercial unit was produced and deployed in a field test. The laboratory unit was used to test many different configurations including thermoelectric, Stirling and vapor compression cooling systems. The final configuration used a vapor compression cooling cycle, vacuum insulation, a passive condenser, an integral evaporator/ thermal storage tank, two 77 watt PV panels and the novel controller mentioned above. The system's only moving part was the variable speed BD35 compressor made by Danfoss. The 365 liter cabinet stayed cold with as little as 274 watt-hours per day average PV power. Battery-free testing was conducted for several months with very good results. The amount of thermal storage, size of compressor and power of PV panels connected can all be adjusted to optimize the design for a given application and climate. In the commercial unit, the high cost of the vacuum insulated refrigerator cabinet and the stainless steel thermal storage tank were addressed in an effort to make the technology commercially viable. This unit started with a 142 liter, mass-produced chest freezer cabinet that had the evaporator integrated into its inner walls. Its compressor was replaced with a Danfoss DC compressor slightly larger than the one used in the laboratory unit. The control system was integrated onto a single electronics card and packaged with its starting capacitors. The water for thermal storage was placed behind a liner that was made to fit inside the original factory liner. The original condenser was also augmented with additional surface area to improve performance. PV panels with a total rated power of 180 watts were used. The unit was tested with very successful results in an outside ambient environment, demonstrating its potential for widespread use in many off-grid applications for solar refrigeration.

Description: The parabolized stability equations (PSE) are used to investigate issues of nonlinear flow development and mixing in compressible reacting shear layers. Particular interest is placed on investigating the change in flow structure that occurs when compressibility and heat release are added to the flow. These conditions allow the 'outer' instability modes- one associated with each of the fast and slow streams-to dominate over the 'central', Kelvin-Helmholtz mode that unaccompanied in incompressible nonreacting mixing layers. Analysis of scalar probability density functions in flows with dominant outer modes demonstrates the ineffective, one-sided nature of mixing that accompany these flow structures. Colayer conditions, where two modes have equal growth rate and the mixing layer is formed by two sets of vortices, offer some opportunity for mixing enhancement. Their extent, however, is found to be limited in the mixing layer's parameter space. Extensive validation of the PSE technique also provides a unique perspective on central- mode vortex pairing, further supporting the view that pairing is primarily governed perspective sheds insight on how linear stability theory is able to provide such an accurate prediction of experimentally-observed, fully nonlinear flow phenomenon.

Description: The flow structure of a pulsating helium jet was investigated using quantitative rainbow schlieren deflectometry. The operating parameters included the tube inside diameter, the jet Reynolds number and the jet Richardson number. The jet structure was characterized by the frequency spectrum, temporal evolution of the oscillations and mean and root-mean-square profiles of the species mole fraction, which in this case, was helium. Experiments were conducted using a variable nozzle facility. Angular deflection data were obtained using rainbow schlieren deflectometry across full field of color images taken at a temporal resolution of 60Hz. The flicker cycle was analyzed by instantaneous rainbow schlieren images, contour plots of deflection angle and helium mole fraction placed sequentially in time. Two observe the flow downstream, the tube was lowered with respect to the optical setup and images were taken at different axial planes. In this way, the laminar, transition and turbulent regions of the jet flow were observed. Abel inversion algorithm was utilized to reconstruct the refractive index field from the measurements of the beam deflection angle. The concentration field was then generated from the refractive index field. The phenomenon of vortex initiation and propagation in the flicker cycle was described by correlating the ray deflection angle and concentration contour plots. Experiments in the flickering jet reveal global oscillation in the flow field. The effect of jet exit Reynolds number and jet Richardson number on the flickering frequency was analyzed. The effect of jet Richardson number was more pronounced than that of the jet exit Reynolds number. The flow field was studies quantitatively in terms of temporal evolution and statistical description of helium mole fraction.

Description: Richtmyer-Meshkov (R-M) 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 possibly become turbulent given the proper initial conditions. R-M instability is similar to the Rayleigh-Taylor (R-T) instability, which is generated when the two fluids undergo a constant acceleration. R-M instability is a fundamental fluid instability that is important to fields ranging from astrophysics to high-speed combustion. For example, R-M instability is currently the limiting factor in achieving a net positive yield with inertial confinement fusion. The experiments described here utilize a novel technique that circumvents many of the experimental difficulties previously limiting the study of the R-M instability. A Plexiglas tank contains two unequal density liquids and is gently oscillated horizontally to produce a controlled initial fluid interface shape. The tank is mounted to a sled on a high speed, low friction linear rail system, constraining the main motion to the vertical direction. The sled is released from an initial height and falls vertically until it bounces off of a movable spring, imparting an impulsive acceleration in the upward direction. As the sled travels up and down the rails, the spring retracts out of the way, allowing the instability to evolve in free-fall until impacting a shock absorber at the end of the rails. The impulsive acceleration provided to the system is measured by a piezoelectric accelerometer mounted on the tank, and a capacitive accelerometer measures the low-level drag of the bearings. Planar Laser-Induced Fluorescence is used for flow visualization, which uses an Argon ion laser to illuminate the flow and a CCD camera, mounted to the sled, to capture images of the interface. This experimental study investigates the instability of an interface between incompressible, miscible liquids with an initial sinusoidal perturbation. The amplitude of the disturbance during the experiment is measured and compared to theory. The results show good agreement (within 10%) with linear stability theory up to nondimensional amplitude ka = 0.7 (wavenumber x amplitude). These results hold true for an initial ka (before acceleration) of -0.7 less than ka less than -0.06, while the linear theory was developed for absolute value of ka much less than 1. In addition, a third order weakly nonlinear perturbation theory is shown to be accurate for amplitudes as large as ka = 1.3, even though the interface becomes double-valued at ka = 1.1. As time progresses, the vorticity on the interface concentrates, and the interface spirals around the alternating sign vortex centers to form a mushroom pattern. At higher Reynolds Number (based on circulation), an instability of the vortex cores has been observed. While time limitations of the apparatus prevent determination of a critical Reynolds Number, the lowest Reynolds Number this vortex instability has been observed at is 5000.

Description: Recent experiments have shown that low-density gas jets injected into a high-density ambient gas undergo an instability mode, leading to highly-periodic oscillations in the flow-field for certain conditions. The transition from laminar to turbulent flow in these jets is abrupt, without the gradual change in scales. Even the fine scale turbulent structure repeats itself with extreme regularity from cycle to cycle. Similar observations were obtained in buoyancy-dominated and momentum-dominated jets characterized by the Richardson numbers, Ri = [gD(rho(sub a)-rho(sub j))/rho(sub j)U(sub j)(exp 2) ] where g is the gravitational acceleration, D is the jet diameter, rho(sub a) and rho(sub a) are, respectively, the free-stream and jet densities, and U(sub j) is the mean jet exit velocity. At high Richardson numbers, the instability is presumably caused by buoyancy since the flow-oscillation frequency (f) or the Strouhal number, St = [fD/U(sub j)] scales with Ri. In momentum-dominated jets, however, the Strouhal number of the oscillating flow is relatively independent of the Ri. In this case, a local absolute instability is predicted in the potential core of low-density jets with S [= rho(sub j)/rho(sub a)] 〈 0.7, which agrees qualitatively with experiments. Although the instability in gas jets of high Richardson numbers is attributed to buoyancy, direct physical evidence has not been acquired in experiments. If the instability is indeed caused by buoyancy, the near-field flow structure of the jet will change significantly when the buoyancy is removed, for example, in the microgravity environment. Thus, quantitative data on the spatial and temporal evolutions of the instability, length and time scale of the oscillating mode and its effects on the mean flow and breakdown of the potential core are needed in normal and microgravity to delineate gravitational effects in buoyant jets. In momentum dominated low-density jets, the instability is speculated to originate in the potential core. However, experiments have not succeeded in identifying the direct physical cause of the instability. For example, the theory predicts an oscillating mode for S〈0.62 in the limit of zero momentum thickness, which contradicts with the experimental findings of Kyle and Sreenivasan. The analyses of momentum-dominated jets neglect buoyancy effects because of the small Richardson number. Although this assumption is appropriate in the potential core, the gravitational effects are important in the annular region surrounding the jet, where the density and velocity gradients are large. This reasoning provides basis for the hypothesis that the instability in low Richardosn number jets studied by Kyle and Sreenivasan and Monkewitz et al. is caused by buoyancy. The striking similarity in characteristics of the instability and virtually the identical conclusions reached by Subbarao and Cantwell in buoyant (Ri〉0.5) helium jets on one hand and by Kyle and Sreenivasan in momentum-dominated (Ri〈1x10(exp -3)) helium jets on the other support this hypothesis. However, quantitative experiments in normal and microgravity are necessary to obtain direct physical evidence of buoyancy effects on the flow instability and structure of momentum-dominated low-density jets. The primary objective of this new research project is to quantify how buoyancy affects the flow instability and structure in the near field of low-density jets. The flow will be described by the spatial and temporal evolutions of the instability, length and time scales of the oscillating mode, and the mean and fluctuating concentration fields. To meet this objective, concentration measurements will be obtained across the whole field using quantitative Rainbow Schlieren Deflectometry, providing spatial resolution of 0.1mm and temporal resolution of 0.017s to 1ms. The experimental effort will be supplemented with linear stability analysis of low-density jets by considering buoyancy. The first objective of this research is to investigate the effects of gravity on the flow instability and structure of low-density jets. The flow instability in these jets has been attributed to buoyancy. By removing buoyancy in our experiments, we seek to obtain the direct physical evidence of the instability mechanism. In the absence of the instability, the flow structure will undergo a significant change. We seek to quantify these changes by mapping the flow field (in terms of the concentration profiles) of these jets at non-buoyant conditions. Such information is presently lacking in the existing literature. The second objective of this research is to determine if the instability in momentum-driven, low-density jets is caused by buoyancy. At these conditions, the buoyancy effects are commonly ignored because of the small Richardson based on global parameters. By eliminating buoyancy in our experiments, globally as well as locally, we seek to examine the possibility that the instability mechanism in self-excited, buoyant or momentum-driven jets is the same. To meet this objective, we would quantify the jet flow in normal and microgravity, while systematically decreasing the Richardson number from buoyancy-driven to momentum driven flow regime. The third objective of this research is to perform a linear stability analysis of low-density gas jets by including the gravitational effects. The flow oscillations in these jets are attributed to an absolute instability, whereby the disturbance grows exponentially at the site to ultimately contaminate the entire flow field. We seek to study the characteristics of both convective and absolute instabilities and demarcate the boundary between them.

Description: Most of the second year of our research program focused on exploring the potentially favorable the effects of expansion waves on homogeneous and isotropic turbulence, which is formed downstream of a grid. Expansion waves are associated with compressible flows and may reduce the drag over airfoils by suppressing turbulence. In the very few previous investigations of interactions of turbulence with expansion waves the effects due to stabilizing streamline curvature substantially masked the effects of turbulence suppression due to flow expansion though the waves. In the present flow configuration planar expansion waves interact with grid generated turbulence in our high-resolution shock tube research facility. This approach will assess directly the effects of the interaction on turbulence. The first objective of our study was to identify the nature of expansion waves present in our shock tube facility. Our time-dependent numerical simulations of the flow in our facility indicated the existence of two regions of traveling expansion waves. The system of expansion waves utilized in this investigation is generated by the exiting shock wave and the induced flow behind it at the end of the driver. Several new measuring techniques are being developed which are capable of providing velocity-gradient-related quantities in compressible flows for the first time.

Description: Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is applied to the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow which propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillarity to direct the movement of bubbles in space is the fact that surfactant impurities which are unavoidably present in the continuous phase can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature tension gradient, and decreases to near zero the thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have been retarded by the adsorption onto the bubble surface of a surfactant impurity, Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition the rapid kinetic exchange and high bulk concentration maintain a saturated surface with a uniform surface concentrations. This prevents retarding surface tension gradients and keeps the velocity high. In our first report last year, we detailed experimental results which verified the theory of remobilization in ground based experiments in which the steady velocity of rising bubbles was measured in a continuous phase consisting of a glycerol/water mixture containing a polyethylene glycol surfactant C12E6 (CH3(CH2)11(OCH2CH2)6OH). In our report this year, we detail our efforts to describe theoretically the remobilization observed. We construct a model in which a bubble rises steadily by buoyancy in a continuous (Newtonian) viscous fluid containing surfactant with a uniform far field bulk concentration. We account for the effects of inertia as well as viscosity in the flow in the continuous phase caused by the bubble motion (order one Reynolds number), and we assume that the bubble shape remains spherical (viscous and inertial forces are smaller than capillary forces, i e. small Weber and capillary numbers). The surfactant distribution is calculated by solving the mass transfer equations including convection and diffusion in the bulk, and finite kinetic exchange the bulk and the surface. Convective effects dominate diffusive mass transfer in the bulk of the liquid (high Peclet numbers) except in a thin boundary layer near the surface. A finite volume method is used to numerically solve the hydrodynamic and mass transfer equations on a staggered grid which accounts specifically for the thin boundary layer. We present the results of the nondimensional drag as a function of the bulk concentration of surfactant for different rates of kinetic exchange, from which we develop criteria for the concentration necessary to develop a prescribed degree of remobilization. The criteria compare favorably with the experimental results.

Description: When a fluid interface with surfactants is at rest, the interfacial stress is isotropic (as given by the equilibrium interfacial tension), and is described by the equation of state which relates the surface tension to the surfactant surface concentration. When surfactants are subjected to shear and dilatational flows, flow induced interaction of the surfactants; can create interfacial stresses apart from the equilibrium surface tension. The simplest relationship between surface strain rate and surface stress is the Boussinesq-Scriven constitutive equation completely characterized by three coefficients: equilibrium interfacial tension, surface shear viscosity, and surface dilatational viscosity Equilibrium interfacial tension and surface shear viscosity measurements are very well established. On the other hand, surface dilatational viscosity measurements are difficult because a flow which change the surface area also changes the surfactant surface concentration creating changes in the equilibrium interfacial tension that must be also taken into account. Surface dilatational viscosity measurements of existing techniques differ by five orders of magnitude and use spatially damped surface waves and rapidly expanding bubbles. In this presentation we introduce a new technique for measuring the surface dilatational viscosity by contracting an aqueous pendant drop attached to a needle tip and having and insoluble surfactant monolayer at the air-water interface. The isotropic total tension on the surface consists of the equilibrium surface tension and the tension due to the dilation. Compression rates are undertaken slow enough so that bulk hydrodynamic stresses are small compared to the surface tension force. Under these conditions we show that the total tension is uniform along the surface and that the Young-Laplace equation governs the drop shape with the equilibrium surface tension replaced by the constant surface isotropic stress. We illustrate this technique using DPPC as the insoluble surfacant monolayer and measured for it a surface dilatational viscosity in the LE phase that is 20 surface poise.

Description: Over the last contract year, a numerical procedure for combined conduction-radiation heat transfer using unstructured grids has been developed. As a result of this research, one paper has been published in the Numerical Heat Transfer Journal. One paper has been accepted for presentation at the International Center for Heat and Mass Transfer's International Symposium on Computational Heat Transfer to be held in Australia next year. A journal paper is under review by my NASA's contact. A conference paper for the ASME National Heat Transfer conference is under preparation. In summary, a total of four (4) papers (two journal and two conference) have been published, accepted or are under preparation. There are two (2) to three (3) more papers to be written for the project. In addition to the above publications, one book chapter, one journal paper and six conference papers have been published as a result of this project. Over the last contract year, the research project resulted in one Ph.D. thesis and partially supported another Ph.D. student. My NASA contact and myself have formulated radiation heat transfer procedures for materials with different indices of refraction and for combined conduction-radiation heat transfer. We are trying to find other applications for the procedures developed under this grant.

Description: Several beta-adrenergic receptor (bAR) agonists are known to cause hypertrophy of skeletal muscle tissue. Because it seems logical that these agonists exert their action on muscle through stimulation of cAMP synthesis, five bAR agonists encompassing a range in activity from strong to weak were evaluated for their ability to stimulate cAMP accumulation in embryonic chicken skeletal muscle cells in culture. Two strong agonists (epinephrine and isoproterenol), one moderate agonist (albuterol), and two weak agonists known to cause hypertrophy in animals (clenbuterol and cimaterol) were studied. Dose response curves were determined over six orders of magnitude in concentration for each agonist, and values were determined for their maximum stimulation of cAMP synthesis rate (Bmax) and the agonist concentration at which 50% stimulation of cAMP synthesis (EC50) occurred. Bmax values decreased in the following order: isoproterenol, epinephrine, albuterol, cimaterol, clenbuterol. Cimaterol and clenbuterol at their Bmax levels were approximately 15-fold weaker than isoproterenol in stimulating the rate of cAMP synthesis. In addition, the EC50 values for isoproterenol, cimaterol, clenbuterol, epinephrine, and albuterol were 360 nM, 630 nM, 900 nM, 2,470 nM, and 3,650 nM, respectively. Finally, dose response curves show that the concentrations of cimaterol and clenbuterol in culture media at concentrations known to cause significant muscle hypertrophy in animals had no detectable effect on stimulation of CAMP accumulation in chicken skeletal muscle cells.

Description: Electromagnetic levitation (EML) is an important tool in materials research. Because a sample can be processed without contact with a container, experiments may be performed on high temperature, highly reactive, and undercooled liquid metals. Many of these experiments are affected by fluid flow in the sample, driven by the electromagnetic positioning force. Despite the importance of convection in these experiments, the transition to turbulence is not well understood in this system. However, we have observed a transition from laminar to turbulent flow in EML droplets in the course of microgravity experiments in TEMPUS on the Space Shuttle (STS-94). The transition occurs repeatably and over a narrow range of conditions. These experimental observations are compared with two competing theories about the transition to turbulence. Also, the results of a particle tracking study of the instabilities leading up to the transition to turbulence are presented.

Description: The influence of gravity on the two-phase flow dynamics is obvious.As the gravity level is reduced,there is a new balance between inertial and interfacial forces, altering the behavior of the flow. In bubbly flow,the absence of drift velocity leads to spherical-shaped bubbles with a rectilinear trajectory.Slug flow is a succession of long bubbles and liquid slug carrying a few bubbles. There is no flow reversal in the thin liquid film as the long bubble and liquid slug pass over the film. Although the flow structure seems to be simpler than in normal gravity conditions,the models developed for the prediction of flow behavior in normal gravity and extended to reduced gravity flow are unable to predict the flow behavior correctly.An additional benefit of conducting studies in microgravity flows is that these studies aide the development of understanding for normal gravity flow behavior by removing the effects of buoyancy on the shape of the interface and density driven shear flows between the gas and the liquid phases. The proposal calls to study specifically the following: 1) The dynamics of isolated bubbles in microgravity liquid flows will be analyzed: Both the dynamics of spherical isolated bubbles and their dispersion by turbulence, their interaction with the pipe wall,the behavior of the bubbles in accelerated or decelerated flows,and the dynamics of isolated cylindrical bubbles, their deformation in accelerated/decelerated flows (in converging or diverging channels), and bubble/bubble interaction. Experiments will consist of the use of Particle Image Velocimetry (PIV) and Laser Doppler Velocimeters (LDV) to study single spherical bubble and single and two cylindrical bubble behavior with respect to their influence on the turbulence of the surrounding liquid and on the wall 2) The dynamics of bubbly and slug flow in microgravity will be analyzed especially for the role of the coalescence in the transition from bubbly to slug flow (effect of fluid properties and surfactant), to identify clusters that promote coalescence and transition the void fraction distribution in bubbly and slug flow,to measure the wall friction in bubbly flow. These experiments will consist of multiple bubbles type flows and will utilize hot wire and film anemometers to measure liquid velocity and wall shear stress respectively and double fiber optic probes to measure bubble size and velocity as a function of tube radius and axial location.

Description: Mixing of two fluids generated by steady and particularly g-jitter acceleration is fundamental towards the understanding of transport phenomena in a microgravity environment. We propose to carry out flight and ground-based experiments to quantify flow fields due to g-jitter type of accelerations using Stereo Imaging Velocimetry (SIV), and measure the concentration field using laser fluorescence. The understanding of the effects of g-jitter on transport phenomena is of great practical interest to the microgravity community and impacts the design of experiments for the Space Shuttle as well as the International Space Station. The aim of our proposed research is to provide quantitative data to the community on the effects of g-jitter on flow fields due to mixing induced by buoyancy forces. The fundamental phenomenon of mixing occurs in a broad range of materials processing encompassing the growth of opto-electronic materials and semiconductors, (by directional freezing and physical vapor transport), to solution and protein crystal growth. In materials processing of these systems, crystal homogeneity, which is affected by the solutal field distribution, is one of the major issues. The understanding of fluid mixing driven by buoyancy forces, besides its importance as a topic in fundamental science, can contribute towards the understanding of how solutal fields behave under various body forces. The body forces of interest are steady acceleration and g-jitter acceleration as in a Space Shuttle environment or the International Space Station. Since control of the body force is important, the flight experiment will be carried out on a tunable microgravity vibration isolation mount, which will permit us to precisely input the desired forcing function to simulate a range of body forces. To that end, we propose to design a flight experiment that can only be carried out under microgravity conditions to fully exploit the effects of various body forces on fluid mixing. Recent flight experiments, by the P.I. through collaboration with the Canadian Space Agency (STS-85, August 1997), aimed at determining the stability of the interface between two miscible liquids inside an enclosure show that a long liquid column (5 cm) under microgravity isolation conditions can be stable, i.e. the interface remains sharp and vertical over a short time scale; thus transport occurs by molecular mass diffusion. On the other hand, when the two liquids were excited from a controlled vibration source (Microgravity Vibration Isolation Mount) two to four mode large amplitude quasi-stationary waves were observed. The data was limited to CCD recording of the dynamics of the interface between the two fluids. We propose to carry out flight experiments to quantify the dynamics of the flow field using Stereo Imaging Velocimetry and measure the concentration field using laser fluorescence. The results will serve as a basis to understand effects of g-jitter on transport phenomena, in this case mass diffusion. As the measurement of the kinematics of the flow field will shed light on the instability mechanism. The research will allow measurement of the flow field in microgravity environment to prove two hypotheses: (1) Maxwell's hypothesis: finite convection always exists in diffusing systems, and (2) Quasi-stationary waves inside a bounded enclosure in a microgravity environment is generated by Kelvin-Helmholtz instability; resonance of the interface which produces incipient mixing is due to Rayleigh-Taylor instability. The first hypothesis can be used as a benchmark experiment to illustrate diffusive mixing. The second hypothesis will lead to the understanding of g-jitter effects on buoyancy driven flow fields which occur in many situations involving materials processing, and other basic fluid physics phenomena. In addition, the second hypothesis will also provide insight in how Rayleigh-Taylor and Kelvin-Helmholtz instabilities propagate concentration fronts during mixing. Measurement of the flow field using SIV is important because it is the flow field which causes instability at the interface between the two fluids. Mixing driven by buoyancy induced flow fields will be addressed both experimentally and computationally. The experimental effort will address the kinematics of mixing: stretching, transport and chaos. Quantification of the mechanisms of mixing will consists of measuring the flow field using the SIV system at Glenn and capturing the dynamics of the interface, to measure mass transport, using a CCD camera. These experiments will be carried out within the framework of Earth's gravity and g-jitter microgravity acceleration as in a Space Shuttle environment or the International Space Station. The g-jitter will be induced and controlled using a tunable vibration isolation platform to isolate against vibration as well as input periodic and random vibration to the system. The parametric range of the microgravity experiment will be extended from the experiments on STS-85 to investigate higher mode quasi-stationary waves (8 to 12), as well as resonance regions which leads to chaos and turbulence. Ground-based experiments will focus on effects of vibration on stably stratified fluid layers in order to scale for possible scenarios in a microgravity environment. These vibrations will be subjected perpendicular to the concentration field on the ground since the parallel case can only be carried out in a microgravity environment. The concept of dynamical similarity will be applied to tune the experiments as closely as possible to a Space Shuttle environment or the International Space Station. The computational effort will take advantage of the Computational Laboratory at Glenn to corroborate the experimental findings with predictions of the dynamics of the flow field using the codes FLUENT (finite difference based) and FIDAP (finite element based). We will investigate two important cases, single-fluid model to address dilute systems with negligible jump in viscosity and the more general two-fluid model which accounts for finite jump in viscosity. Apart from its microgravity relevance, this experiment is well suited to study dynamics in nonlinear systems.

Description: Our recent discovery of the spontaneous formation of chiral domains in fluid smectic phases of achiral bow-shaped molecules opens up a wide variety of possibilities for new liquid crystal phases and phenomena. The basic, spontaneously chiral layer structure of the highest temperature fluid smectic phases, the B2 and B7, are shown. One of the most intriguing aspects of this structure is the plethora of possible phases coming from different stacking sequences of the polar ordering and tilt directions. The four possibilities of next-nearest neighbor alternation are shown. In the original material studied, NOBOW, the ground states found are antiferroelectric, either the racemic SmC(sub S)P(sub A) or the chiral SmC(sub A)P(sub A). We are currently studying MHOBOW, synthesized by D. Walba which, by virtue of its methyl hexyloxy tail has a tendency to form anticlinic layer interfaces, in the hope of finding a phase with a ferroelectric ground state, either SmC(sub A)P(sub S) or SmC(sub S)P(sub A), which can be obtained in NOBOW only by applying a field. Preliminary observations of MHO-BOW have made its study, from the point of view of understanding novel LC structures, extremely high priority. The following truly remarkable characteristics have been revealed: (i) The smectic phase grows out of the isotropic in the form of helical ribbons. The resulting planar aligned textures of focal conics with layers normal to glass plates exhibit bizarre modulations, including stripes and checker-boards. These have also been seen in other materials suggesting that this is a new phase (tentatively called B7), which is a fluid smectic with some kind of in-layer structure. (ii) It is virtually impossible to make freely suspended films of MHOBOW. Rather it makes the freely suspended filaments which preliminary x-ray scattering experiments reveal to have the nested cylinder layer structure indicated; (iii) The powder x-ray diffraction exhibits four resolution-limited smectic layering peaks, very close in layer spacing, which vary continuously with T. This is further evidence for a more complex three dimensional structure than NOBOW, which has a typical single layering reflection. (iv) The x-ray structure factor of the layering peak of the filaments is extraordinarily complex and rich. Varying in qL (the scattering vector component along the filament axis) from a double slit-like pattern to modulated layer-like patterns, as qH (the scattering vector component normal to the filament axis) is varied over the range where the four powder peaks are located. These results suggest some kind of mosaic structure, perhaps with different layer spacings corresponding to the different stacking sequences. Recent x-ray diffraction experiments show that the peaks are modulated in intensity upon translation along a filament, in domains of several hundred microns dimension. These preliminary experiments suggest that the B7 is a fluid smectic with extremely unusual and fascinating structures. Of all of the many hundreds of fluid smectic materials we have attempted to study in the freely suspended film geometry over the years, only a few have failed to form films, and none showed any great tendency to form filaments, although this clearly should be a possible freely suspended smectic LC morphology. On several occasions in the past we have intentionally tried to make filaments from a variety of smectics without success. Thus the smectic filament formation property makes the B7 phase unique. It seems quite likely that the stability of filaments is related to the in-plane structure. The filaments exhibit other interesting structural and optical features. They are birefringent with a local optic axis which is oblique and which can vary continuously along filament and which can be manipulated with an electric field applied normal to the fiber, as if the field were causing a rotation of the optic axis about the fiber axis. Rapid displacement of the ends of the fiber toward one another causes a macroscopic helixing at low T and causes thick regions to transiently appear at high T, a 1D analog of island formation on a rapidly compressed film.

Description: Microgravity technologies often require aqueous phases to spread over nonwetting hydrophobic solid/surfaces. At a hydrophobic surface, the air/hydrophobic solid tension is low, and the solid/aqueous tension is high. A large contact angle forms as the aqueous/air tension acts together with the solid/air tension to balance the large solid/aqueous tension. The aqueous phase, instead of spreading, is held in a meniscus by the large angle. Surfactants facilitate the wetting of water on hydrophobic surfaces by adsorbing on the water/air and hydrophobic solid/water interfaces and lowering the surface tensions of these interfaces. The tension reductions decrease the contact angle, which increases the equilibrium wetted area. Hydrocarbon surfactants (i.e. amphiphiles with a hydrophobic chain of methylene groups attached to a large polar group to give aqueous solubility) do not reduce significantly the contact angles of the very hydrophobic surfaces such as parafilm or polyethylene. Trisiloxane surfactants (amphiphiles with a hydrophobe consisting of methyl groups linked to a trisiloxane backbone in the form of a disk ((CH3)3-Si-O-Si-O-Si(CH3)3)) and an extended ethoxylate (-(OCH2CH2)n-) polar group in the form of a chain with seven or eight units) can significantly reduce the contact angle of water on a very hydrophobic surface and cause rapid and complete (or nearly complete) spreading (lermed superspreading). The overall goal of the research described in this proposal is to establish and verify a theory for how trisiloxanes cause superspreading, and then use this knowledge as a guide to developing more general hydrocarbon based surfactant systems which superspread and can be used in microgravity. We propose that the trisiloxane surfactants superspread when the siloxane adsorbs, the hydrophobic disk parts of the molecule adsorb onto the surface removing the surface water. Since the cross sectional area of the disk is larger than that of the extended ethoxylate chain, the disks can form a space filling mat on the surface which removes a significant amount of the surface water. The water adjacent to the hydrophobic solid surface is of high energy due to incomplete hydrogen bonding; its removal significantly lowers the tension and reduces the contact angle. Hydrocarbon surfactants cannot remove as much surface water because their large polar groups prevent the chains from cohering lengthwise. In our report last year we presented a poster describing the preparation of model very hydrophobic surfaces which are homogeneous and atomically smooth using self assembled monolayers of octadecyl trichlorosilane (OTS). In this poster we will use these surfaces as test substrates in developing hydrocarbon based surfactant systems which superspread. We studied a binary hydrocarbon surfactant systems consisting of a very soluble large polar group polyethylene oxide surfactant (C12E6 (CH3(CH2)11(OCH2CH2)6OH) and a long chain alcohol dodecanol. By mixing the alcohol with this soluble surfactant we have found that the contact angle of the mixed system on our test hydrophobic surfaces is very low. We hypothesize that the alcohol fills in the gaps between adjacent adsorbed chains of the large polar group surfactant. This filling in removes the surface water and effects the decrease in contact angle. We confirm this hypothesis by demonstrating that at the air/water interface the mixed layer forms condensed phases while the soluble large polar group surfactant by itself does not. We present drop impact experiments which demonstrate that the dodecanol/C12E6 mixture is effective in causing impacting drops to spread on the very hydrophobic model OTS surfaces.

Description: Development of HPF versions of NPB and ARC3D has shown that HPF provides an efficient, concise way to express parallelism and to organize data traffic. The use of HPF, as noted in the papers, requires an intimate knowledge of the applications and a detailed analysis of data affinity, data movement, and data granularity. To simplify and accelerate the task of developing HPF versions of existing CFD applications we have designed and implemented ADAPT (Automatic Data Alignment and Placement Tool). ADAPT analyzes a CFD application working on a single structured grid and generates HPF TEMPLATE, (RE)DISTRIBUTION, ALIGNMENT, and INDEPENDENT directives. The directives can be generated on the nest level, subroutine level, application level, or on the application interface level. ADAPT annotates an existing CFD FORTRAN application, performing computations on single or multiple grids. On each grid the application is considered as a sequence of operators, each applied to a set of variables defined in a particular grid domain. ADAPT automatically detects implicit operators (i.e., having data dependences) and explicit operators (without data dependences). For parallelization of an explicit operator ADAPT creates a template for the operator domain, aligns arrays used in the operator with the template, distributes the template, and declares the loops over the distributed dimensions as INDEPENDENT. For parallelization of an implicit operator, the distribution of the operator's domain should be consistent with the operator's dependences. Any dependence between sections distributed on different processors would preclude parallelization if the compiler does not have an ability to pipeline computations. If a data distribution is "orthogonal" to the dependences of an implicit operator, then the loop which implements the operator can be declared as INDEPENDENT. ADAPT starts with an analysis of array index expressions of the loop nests. For each pair of arrays referenced in an assignment statement, it generates an arc in the alignment graph and annotates it with an affinity relation. The template, alignment, and distribution directives for a particular loop nest are then derived from a transitive closure of the affinity relation. A compromise of data distributions in different nests and subroutines is achieved by merging annotated alignment graphs for adjacent nests/stibroutine calls in the nest/call graph of the application in the process called distribution lifting. ADAPT has been implemented as a C++ program running in conjunction with a parallelization tool called CAPTools. ADAPT uses the parse tree, interprocedural analysis and application database generated by CAPTools. It also uses the Directed Graph class, initially implemented in p2d2 (parallel debugger oi distributed programs), and some other classes supporting symbolic computations. ADAPT uses data distribution techniques described. ADAPT was tested with ARC3D and the FT benchmark and has demonstrated a code performance within a factor of 1.5 of handwritten versions.

Description: During hypersonic flight, high temperatures and high heat fluxes are generated on the surfaces of vehicles. The Flight Loads Laboratory (FLL) at Dryden Flight Research Center (DFRC) is equipped with a calibration furnace, capable of calibrating heat flux gages up to 1100kW per square meters, and temperature sensors up to 2600 C. One heating configuration of the calibration furnace is a cylindrical blackbody cavity. Throughout the blackbody there are temperature gradients due to various boundary conditions. These boundary conditions include resistance heating, radiant heat transfer, and conduction to water-cooled electrodes. Also, an inert gas is purged through the graphite blackbody to prevent it from oxidizing. Consequently, the various modes of heat transfer present during operation of the blackbody cavity must be well understood in order to produce accurate heat flux gage and temperature sensor calibrations for use in ground testing or flight testing of hypersonic vehicles. The first step towards understanding the heat transfer in the blackbody cavity was to perform experiments at 1100 C, with and without outer surface insulation, while taking detailed temperature measurements inside the blackbody cavity. Steady state thermal models of the blackbody cavity were then developed. These models included detailed thermal analysis using commercial thermal analysis software. Conduction, radiation, and convection were considered in the thermal models for two cases: one with the outside of the blackbody cavity insulated and the second without insulation. This paper describes the experimental and numerical efforts used to characterize the steady state operation of the blackbody cavity. It describes the analysis of the test measurements, the boundary conditions used in the numerical models, and how the models were calibrated to fit the experimental data. Effects of various uncertainties, such as material properties, and convection are discussed.Initial thermal models predicted temperatures in the deepest part of the blackbody cavity within 7 C of the measured value and produced trends comparable to the experimental data, throughout the models. Adjustment of the boundary conditions, which were included in the thermal models, produced good agreement with measured temperatures. Free and forced convection of the purge gas inside the blackbody was found to be insignificant.

Description: There are innumerable two-component systems in which two very different liquid phases co-exist in equilibrium over a range of temperature and composition, e.g., oil and water, salt fluxes and solders, aluminum and lead. Often it is of practical concern to fabricate a solid component consisting of a uniform dispersion of one phase in the other. Unfortunately, uniform microstructural development during solidification of two immiscible liquids is hampered by inherent, often large, density differences between the phases that lead to severe segregation. Uniformity is also compromised by preferential wetting and coalescence phenomena. It is, however, well known that ultrasonic energy can initiate and maintain a fine liquid-liquid dispersion. The work presented here extends that observation by application of ultrasonic energy to promote uniform phase incorporation during controlled directional solidification. To this end experiments with the transparent organic, immiscible, succinonitrile-glycerol system were conducted and the numerous processing parameters associated with this technique were evaluated in view of optimizing dispersion uniformity. In view of the initial experimental results a model that predicts the dispersed liquid droplet size as a function of material properties, sample geometry, and applied energy has been developed. In the mathematical model we consider the ultrasonic field in an experimental ampoule of length L and diameter D induced by a probe having a vibration frequency of f=2OKhz (circular frequency omega = 2 pi f). The amplitude is adjustable from A=65 to 13Omicrons. The probe tip diameter is d, the liquid has a density of p, in which the speed of sound and surface tension are, respectively, c and sigma. The mathematical model and numerical investigation for the experiments [1] is done using the following assumptions: (i) The droplet size is small in comparison to the sound wave length; (ii) The forces between droplets are neglected (relative concentration is small); (iii) The droplet is stable if the kinetic energy, E(sub K), of the liquid motion due to ultrasonic field influence is less then the binding energy, E(sub S), due to the surface tension (it is easy to show that the surface energy of two droplets resulting from one is larger by about a factor of two.); (iv) The stability limit is characterized by E(sub S) to approx. E(sub K).

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 performed by electromagnetic levitation. The natural oscillation frequency of 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 k-epsilon 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.

Description: Nucleate boiling, especially near the critical heat flux (CHF), can provide excellent economy along with high efficiency of heat transfer. However, the performance of nucleate boiling may deteriorate in a reduced gravity environment and the nucleate boiling usually has a potentially dangerous characteristic in CHF regime. That is, any slight overload can result in burnout of the boiling surface because the heat transfer will suddenly move into the film-boiling regime. Therefore, enhancement of nucleate boiling heat transfer becomes more important in reduced gravity environments. Enhancing nucleate boiling and critical heat flux can be reached using micro-configured metal-graphite composites as the boiling surface. Thermocapillary force induced by temperature difference between the graphite-fiber tips and the metal matrix, which is independent of gravity, will play an important role in bubble detachment. Thus boiling heat transfer performance does not deteriorate in a reduced-gravity environment. Based on the existing experimental data, and a two-tier theoretical model, correlation formulas are derived for nucleate boiling on the copper-graphite and aluminum-graphite composite surfaces, in both the isolated and coalesced bubble regimes. Experimental studies were performed on nucleate pool boiling of pentane on cooper-graphite (Cu-Gr) and aluminum-graphite (Al-Gr) composite surfaces with various fiber volume concentrations for heat fluxes up to 35 W per square centimeter. It is revealed that a significant enhancement in boiling heat transfer performance on the composite surfaces is achieved, due to the presence of micro-graphite fibers embedded in the matrix. The onset of nucleate boiling (the isolated bubble regime) occurs at wall superheat of about 10 C for the Cu-Gr surface and 15 C for the Al-Gr surface, much lower than their respective pure metal surfaces. Transition from an isolated bubble regime to a coalesced bubble regime in boiling occurs at a superheat of about 14 C on Cu-Gr surface and 19 C on Al-Gr surface.

Description: The molecular structure of biological macromolecules is important in understanding how these molecules work and has direct application to rational drug design for new medicines and for the improvement and development of industrial enzymes. In order to obtain the molecular structure, large, well formed, single macromolecule crystals are required. The growth of macromolecule crystals is a difficult task and is often hampered on the ground by fluid flows that result from the interaction of gravity with the crystal growth process. One such effect is the bulk movement of the crystal through the fluid due to sedimentation. A second is buoyancy driven convection close to the crystal surface. On the ground the crystallization process itself induces both of these flows. Buoyancy driven convection results from density differences between the bulk solution and fluid close to the crystal surface which has been depleted of macromolecules due to crystal growth. Schlieren photograph of a growing lysozyme crystal illustrating a 'growth plume' resulting from buoyancy driven convection. Both sedimentation and buoyancy driven convection have a negative effect on crystal growth and microgravity is seen as a way to both greatly reduce sedimentation and provide greater stability for 'depletion zones' around growing crystals. Some current crystal growth hardware however such as those based on a vapor diffusion techniques, may also be introducing unwanted Marangoni convection which becomes more pronounced in microgravity. Negative effects of g-jitter on crystal growth have also been observed. To study the magnitude of fluid flows around growing crystals we have attached a number of different fluorescent probes to lysozyme molecules. At low concentrations, less than 40% of the total protein, the probes do not appear to effect the crystal growth process. By using these probes we expect to determine not only the effect of induced flows due to crystal growth hardware design but also hope to optimize crystallization hardware so that destructive flows are minimized both on the ground and in microgravity.

Description: Microgravity fluid physics experiments frequently measure concentration and temperature. Interferometers such as the Twyman Green illustrated have performed full-field measurement of these quantities. As with most such devices, this interferometer uses a reference path that is not common with the path through the test section. Recombination of the test and reference wavefronts produces interference fringes. Unfortunately, in order to obtain stable fringes, the alignment of both the test and reference paths must be maintained to within a fraction of the wavelength of the light being used for the measurement. Otherwise, the fringes will shift and may disappear. Because these interferometers are extremely sensitive to bumping, jarring and transmitted vibration, they are typically mounted on optical isolation tables. Schlieren deflectometers or the more recent Shack-Hartmann wavefront sensors also measure concentration and temperature in laboratory fluid flows. Ray optics describe the operation of both devices. In a schlieren system, an expanded, collimated beam passes through a test section where refractive index gradients deflect rays. A lens focuses the beam to a filter placed in the rear focal plane of the decollimating lens. In a quantitative color schlieren system, gradients in the index of refraction appear as colors in the field of view due to the action of the color filter. Since sensitivity is a function of the focal length of the decollimating lens, these systems are rather long and filter fabrication and calibration is rather difficult. A Shack-Hartmann wavefront sensor is an array of small lenslets. Typical diameters are on the order of a few hundred microns. Since these lenslets divide the test section into resolution elements, the spatial resolution can be no smaller than an individual lenslet. Such a device was recently used to perform high-speed tomography of heated air exiting a 1.27 cm diameter nozzle. While these wavefront sensors are very compact, the limited spatial resolution and the methods required for data reduction suggest that a more useful instrument needs to be developed. The category of interferometers known as common path interferometers can eliminate much of the vibration sensitivity associated with traditional interferometry as described above. In these devices, division of the amplitude of the wavefront following the test section produces the reference beam. Examples of these instruments include shearing and point diffraction interferometers. In the latter case, shown schematically, a lens focuses light passing through the test section onto a small diffracting object. Such objects are typically either a circle of material on a high quality glass plate or a small sphere in a glass cell. The size of the focused spot is several times larger than the object so that the light not intercepted by the diffracting object forms the test beam while the diffracted light generates a spherical reference beam. While this configuration is mechanically stable, phase shifting one beam with respect to the other is difficult due to the common path. Phase shifting enables extremely accurate measurements of the phase of the interferogram using only gray scale intensity measurements and is the de facto standard of industry. Mercer and Creath 2 demonstrated phase shifting in a point diffraction interferometer using a spherical spacer in a liquid crystal cell as the diffracting object. By changing the voltage across the cell, they were able to shift the phase of the undiffracted beam relative to the reference beam generated by diffraction from the sphere. While they applied this technology to fluid measurements, the device shifted phase so slowly that it was not useful for studying transient phenomena. We have identified several technical problems that precluded operation of the device at video frame rates and intend to solve them to produce a phase-shifting liquid crystal point-diffraction interferometer operating at video frame rates. The first task is to produce high contrast fringes. Since the diffracted beam is much weaker than the transmitted beam, interferograms have poor contrast unless a dye is added to the liquid crystal to reduce the intensity of the undiffracted light. Dyes previously used were not rigorously characterized and suffered from hysteresis in both the initial alignment state of the device and the electro-optic switching characteristics. Hence, our initial effort will identify and characterize dyes that do not suffer from these difficulties and are readily soluble in the liquid crystal host. Since the ultimate goal of this research is to produce interferometers capable of phase shifting at video frame rates, we will quantify the difference in switching times between ferroelectric and nematic liquid crystals. While we have more experience with nematic crystals, they typically switch more slowly than ferroelectric cells. As part of that effort, we will investigate the difference in the modulation of the interferograms as a function of the type of liquid crystal in the cell. Because the temporal switching response of a liquid crystal cell is directly related its thickness, we intend to explore techniques required to produce cells that are as thin as possible. However, the cells must still produce a total phase shift of two pi radians.

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

Description: Our previous float zone experiments 1,2 with NaNO3 revealed that steady thermocapillary flow (TC flow) can be balanced/offset by the controlled surface streaming flow (CSS flow), induced by end-wall vibration. In the current experiments, we are examining the effects of surface streaming flow on steadying/stabilizing oscillatory thermocapillary flow. To this effect, we have set up a controlled NaNO3 half-zone experiment, where the processing parameters like zone dimensions and temperature gradients scan be easily varied to achieve oscillatory TC flow. In the present paper, we discuss the thermal signature of the TC flow, and how it is affected by imposition of CSS flow. The results will also include a comparison of the microstructure of a NaNO3- BaNO3 eutectic, processed under oscillatory TC conditions, with and without imposed CSS flow.

Description: The distribution of insoluble particles in a metal casting depends primarily on the interaction of the particles with the solid/liquid interface (SLI) during the solidification process. The balance of the forces acting on the particle essentially determines whether a particle will be engulfed or pushed by the SLI. An important component of this force balance is the drag force generated by the particle motion in front of the SLI. Previously developed mathematical models for particle/SLI interaction made use of steady-state solutions of this force provided by the lubrication theory. However, our numerical model based on the SLI tracking approach shows that the steady-state approach is inappropriate to model the interaction process and that at steady-state the theoretical solution underestimates the drag force. It was found that regression analysis of steady-state numerical solutions for cylindrical particles moving normal to a flat SLI gives a relationship of the form: Abstract The distribution of insoluble particles in a metal casting depends primarily on the interaction of the particles with the solid/liquid interface (SLI) during the solidification process. The balance of the forces acting on the particle essentially determines whether a particle will be engulfed or pushed by the SLI. An important component of this force balance is the drag force generated by the particle motion in front of the SLI. Previously developed mathematical models for particle/SLI interaction made use of steady-state solutions of this force provided by the lubrication theory. However, our numerical model based on the SLI tracking approach shows that the steady-state approach is inappropriate to model the interaction process and that at steady-state the theoretical solution underestimates the drag force. It was found that regression analysis of steady-state numerical solutions for cylindrical particles moving normal to a flat SLI gives a relationship of the form: F(sub D, sup num) =sqoare root of 3(pi)(eta)V(sub p)(R(sub p)/d)(sup 10(gamma)/3). This is to compared to the solution provided by the lubrication theory: F(sub D), sup theor) = 3 square root of 2(pi)(eta)V(sub p)(R(sub p)/d)(sub 3/2), where F(sub D) is the drag force, eta is the dynamic viscosity of the fluid, V(sub p) is the particle velocity, R(sub p) is the particle radius, gamma is Euler's constant (gamma approx. = 0.577), and d is the width of the gap between the particle and the SLI. The domain on which the equations proposed by the lubrication theory are relevant is clearly identified in the paper. The numerical model was then validated against the classical lubrication theory within the domain of its validity.

Description: Rotating magnetic fields (RMF) have found applications in modern metallurgy. Examples are casting and semiconductor crystal growth technology. Rotating convective flows can also be induced in ionic fluids. Due to optical transparency, these fluids offer a great opportunity to visualize convection. In this work, we investigate the dynamics of convection induced by a high frequency (100 kHz range) RMF. A dye injected into a cylindrical column of salty water serves as an indicator of the flow. The developed technique has been used to study mixing phenomena induced by a RMF both when the direc'tion of the field rotation is constant and when it is alternated. Optical recording has been used to analyze this process. A numerical model describing RMF mixing in the laminar regime will also be presented.

Description: Most models that describe the interaction of an insoluble particle with an advancing solid-liquid interface are based on the assumption of steady state. However, as demonstrated by experimental work, the process does not reach steady state until the particle is pushed for a while by the interface. In this work, a dynamic mathematical model was developed. The dynamic model demonstrates that this interaction is essentially non-steady state and that steady state eventually occurs only when solidification is conducted at sub-critical velocities. The model was tested for three systems: aluminum-zirconia particles, succinonitrilepolystyrene particles, and biphenyl-glass particles. The calculated values for critical velocity of the pushing/engulfment transition were in same range with the experimental ones.

Description: A procedure that models coupled thermo-mechanical deformations of viscoelastic rubber cylinders by employing the ABAQUS finite element code is described. Computational simulations of hysteretic heating are presented for several tall and short rubber cylinders both with and without a steel disk at their centers. The cylinders are compressed axially and are then cyclically loaded about the compressed state. The non-uniform hysteretic heating of the rubber cylinders containing a steel disk is presented. The analyses performed suggest that the coupling procedure should be considered for further development as a design tool for rubber degradation studies.

Description: The study of the shock structure in a viscous heat conducting fluid is an old problem. We study this problem from a novel mathematical point of view. A new class of generalized functions is defined where multiplication of any two functions is allowed with the usual properties. A Heaviside function in this class has the unit jump at occurring on an infinitesimal interval of the nonstandard analysis (NSA) in the halo of . This jump has a smooth microstructure over the infinitesimal interval . From this point of view, we have a new class of Heaviside functions, and their derivatives the Dirac delta functions, which are equivalent when viewed as continuous linear functionals over the test function space of Schwartz. However, they differ in their microstructures which in applications are determined from physics of the problem as shown in our presentation.

Description: Unseeded molecular tagging methods based on single-photon processes that produce long tag lines (〉50 mm) have been recently developed and demonstrated by the Combustion Laser Diagnostics Group (Mechanical Engineering Department) at Vanderbilt University [1,2]. In Ozone Tagging Velocimetry (OTV) a line of ozone (O3) is produced by a single photon from a pulsed narrowband argon fluoride (ArF) excimer laser operating at - 193 nm. After a known time delay, t, the position of the displaced (convected in the flow field) O3 tag line is revealed by photodissociation of O3 and subsequent fluorescence of O2, caused by a pulsed laser sheet from a krypton fluoride (KrF) excimer laser operating at - 248 nm. Intensified CCD camera images of the fluorescence are taken from the initial and final tag line locations thus providing unobtrusive means of establishing a velocity profile in the interrogated flow field. The O3 lines are "written" and subsequently "read" by the following reactions:

Description: A stable high order numerical scheme for direct numerical simulation (DNS) of shock-free compressible turbulence is presented. The method is applicable to general geometries. It contains no upwinding, artificial dissipation, or filtering. Instead the method relies on the stabilizing mechanisms of an appropriate conditioning of the governing equations and the use of compatible spatial difference operators for the interior points (interior scheme) as well as the boundary points (boundary scheme). An entropy splitting approach splits the inviscid flux derivatives into conservative and non-conservative portions. The spatial difference operators satisfy a summation by parts condition leading to a stable scheme (combined interior and boundary schemes) for the initial boundary value problem using a generalized energy estimate. A Laplacian formulation of the viscous and heat conduction terms on the right hand side of the Navier-Stokes equations is used to ensure that any tendency to odd-even decoupling associated with central schemes can be countered by the fluid viscosity. A special formulation of the continuity equation is used, based on similar arguments. The resulting methods are able to minimize spurious high frequency oscillation producing nonlinear instability associated with pure central schemes, especially for long time integration simulation such as DNS. For validation purposes, the methods are tested in a DNS of compressible turbulent plane channel flow at a friction Mach number of 0.1 where a very accurate turbulence data base exists. It is demonstrated that the methods are robust in terms of grid resolution, and in good agreement with incompressible channel data, as expected at this Mach number. Accurate turbulence statistics can be obtained with moderate grid sizes. Stability limits on the range of the splitting parameter are determined from numerical tests.

Description: In a series of papers, Olsson (1994, 1995), Olsson & Oliger (1994), Strand (1994), Gerritsen Olsson (1996), Yee et al. (1999a,b, 2000) and Sandham & Yee (2000), the issue of nonlinear stability of the compressible Euler and Navier-Stokes Equations, including physical boundaries, and the corresponding development of the discrete analogue of nonlinear stable high order schemes, including boundary schemes, were developed, extended and evaluated for various fluid flows. High order here refers to spatial schemes that are essentially fourth-order or higher away from shock and shear regions. The objective of this paper is to give an overview of the progress of the low dissipative high order shock-capturing schemes proposed by Yee et al. (1999a,b, 2000). This class of schemes consists of simple non-dissipative high order compact or non-compact central spatial differencings and adaptive nonlinear numerical dissipation operators to minimize the use of numerical dissipation. The amount of numerical dissipation is further minimized by applying the scheme to the entropy splitting form of the inviscid flux derivatives, and by rewriting the viscous terms to minimize odd-even decoupling before the application of the central scheme (Sandham & Yee). The efficiency and accuracy of these scheme are compared with spectral, TVD and fifth- order WENO schemes. A new approach of Sjogreen & Yee (2000) utilizing non-orthogonal multi-resolution wavelet basis functions as sensors to dynamically determine the appropriate amount of numerical dissipation to be added to the non-dissipative high order spatial scheme at each grid point will be discussed. Numerical experiments of long time integration of smooth flows, shock-turbulence interactions, direct numerical simulations of a 3-D compressible turbulent plane channel flow, and various mixing layer problems indicate that these schemes are especially suitable for practical complex problems in nonlinear aeroacoustics, rotorcraft dynamics, direct numerical simulation or large eddy simulation of compressible turbulent flows at various speeds including high-speed shock-turbulence interactions, and general long time wave propagation problems. These schemes, including entropy splitting, have also been extended to freestream preserving schemes on curvilinear moving grids for a thermally perfect gas (Vinokur & Yee 2000).

Description: In the past, feature extraction and identification were interesting concepts, but not required in understanding the physics of a steady flow field. This is because the results of the more traditional tools like iso-surfaces, cuts and streamlines, were more interactive and easily abstracted so they could be represented to the investigator. These tools worked and properly conveyed the collected information at the expense of a great deal of interaction. For unsteady flow-fields, the investigator does not have the luxury of spending time scanning only one 'snap-shot' of the simulation. Automated assistance is required in pointing out areas of potential interest contained within the flow. This must not require a heavy compute burden (the visualization should not significantly slow down the solution procedure for co-processing environments like pV3). And methods must be developed to abstract the feature and display it in a manner that physically makes sense.

Description: The objectives of this viewgraph presentation are to develop Chimera-based potential methodology which is compatible with overflow and overflow infrastructure, creating options for an advanced problem solving environment and to significantly reduce turnaround time for aerodynamic analysis and design (primarily cruise conditions).

Description: This paper discusses the current technologies for on-orbit refueling of spacecraft. The findings of 55 references are reviewed and summarized. Highlights include: (1) the Russian Progress system used by the International Space Station; (2) a flight demonstration of superfluid helium transfer; and (3) ground tests of large cryogenic systems. Key technologies discussed include vapor free liquid outflow, control of fluid inflow to prevent liquid venting, and quick disconnects for on-orbit mating of transfer lines.

Description: The coolant flow characteristics at the hole exits of a film-cooled blade are derived from an earlier analysis where the hole pipes and coolant plenum were also discretized. The blade chosen is the VKI rotor with three staggered rows of shower-head holes. The present analysis applies these flow characteristics at the shower-head hole exits. A multi-block three-dimensional Navier-Stokes code with Wilcox's k-omega model is used to compute the heat transfer coefficient on the film-cooled turbine blade. A reasonably good comparison with the experimental data as well as with the more complete earlier analysis where the hole pipes and coolant plenum were also gridded is obtained. If the 1/7th power law is assumed for the coolant flow characteristics at the hole exits, considerable differences in the heat transfer coefficient on the blade surface, specially in the leading-edge region, are observed even though the span-averaged values of h (heat transfer coefficient based on T(sub o)-T(sub w)) match well with the experimental data. This calls for span-resolved experimental data near film-cooling holes on a blade for better validation of the code.

Description: This paper describes the results of a numerical study of interacting hypersonic flows at conditions that can be produced in ground-based test facilities. The computations are made with the direct simulation Monte Carlo (DSMC) method of Bird. The focus is on Mach 10 flows about flared axisymmetric configurations, both hollow cylinder flares and double cones. The flow conditions are those for which experiments have been or will be performed in the ONERA R5Ch low-density wind tunnel and the Calspan-University of Buffalo Research Center (CUBRC) Large Energy National Shock (LENS) tunnel. The range of flow conditions, model configurations, and model sizes provides a significant range of shock/shock and shock/boundary layer interactions at low Reynolds number conditions. Results presented will highlight the sensitivity of the calculations to grid resolution, contrast the differences in flow structure for hypersonic cold flows and those of more energetic but still low enthalpy flows, and compare the present results with experimental measurements for surface heating, pressure, and extent of separation.

Description: Loop heat pipes (LHPs) are versatile two-phase heat transfer devices that have gained increasing acceptance for space and terrestrial applications. The operating temperature of an LHP is a function of its operating conditions. The LHP usually reaches a steady operating temperature for a given heat load and sink temperature. The operating temperature will change when the heat load and/or the sink temperature changes, but eventually reaches another steady state in most cases. Under certain conditions, however, the loop operating temperature never really reaches a true steady state, but instead becomes oscillatory. This paper discusses the temperature oscillation phenomenon using test data from a miniature LHP.

Description: An adaptive unstructured grid refinement technique has been developed and successfully applied to several three dimensional inviscid flow test cases. The method is based on a combination of surface mesh subdivision and local remeshing of the volume grid Simple functions of flow quantities are employed to detect dominant features of the flowfield The method is designed for modular coupling with various error/feature analyzers and flow solvers. Several steady-state, inviscid flow test cases are presented to demonstrate the applicability of the method for solving practical three-dimensional problems. In all cases, accurate solutions featuring complex, nonlinear flow phenomena such as shock waves and vortices have been generated automatically and efficiently.

Description: The effects of rotating cavitation and cavitation surges on the Fastrac Engine Turbopump are described in a viewgraph presentation format. The bent inducer blade dilemma and observations of unsteady data and oscillation components are discussed. The pump-feed system stability modeling assessment is outlined. Recommendations are made urging further investigation.

Description: This paper describes a series of wind tunnel experiments carried out with the aim of providing data suitable for evaluating the performance of Computational Fluid Dynamics (CFD) codes. The configurations and flow conditions studied are most relevant to slender supersonic missiles. However, the data obtained, which includes forces and moments, surface pressures, flowfield surveys and a selection of flow visualization images, should he of interest to other CFD practitioners. Results for three test cases are presented and discussed in this paper. These cases have been the subject of a collaborative study concerned with the evaluation of Navier-Stokes solvers for missiles, carried out under the auspices of The Technical Cooperation Programme (TTCP).

Description: The objective of this effort is to better understand the physics of evaporation, condensation, and fluid flow as they affect the heat transfer processes in a constrained vapor bubble heat exchanger (CVBHX). This CVBHX consists of a small enclosed container with a square cross section (inside dimensions. 3 x 3 x 40 mm) partially filled with a liquid. The major portion of the liquid is in the corners, which act as arteries. When a temperature difference is applied to the ends of the CVBHX, evaporation occurs at the hot end and condensation at the cold end resulting in a very effective heat transfer device with great potential in space applications. Liquid is returned by capillary flow in the corners. A complete description of the system and the results obtained to date are given in the papers listed.

Description: Computational Fluid Dynamics (CFD) analysis was used to compute effective nozzle discharge coefficients for subscale sharp-edged converging/diverging nozzles, with a variety of convergence half-angles, motor operating conditions, and two propellants with different ballistics. Convergence half-angles ranged from 0 to 80 deg. Analysis was conducted at total temperatures from 2946K (5303R) to 3346K (6023R) and over total pressures ranged from 2.72 MPa (395 psia) to 20.68 MPa (3000 psia). Area ratios (A(sub e)/A*) ranged from 7.43 to 9.39. Ratio of specific heats (gamma) ranged from 1.13 to 1.18. Throat and exit Reynolds numbers were calculated to be 8.26 x 10(exp 5) and 5.51 x 10(exp 5), respectively. Present results of nozzle discharge coefficients are reported and correlated as a function of nozzle convergence half-angle (theta(sub c)) and area ratios (A(sub e)/A*) for a constant divergence half-angle (theta(sub d)) of 15 deg. Computed discharge coefficients ranged from 0.88 to 0.97. They are compared with theory and experimental data available in literature. Available turbulence models with respect to grid refinements and heat transfer are discussed.

Description: The effect of accelerating forces on the performance of loop heat pipes (LHP) is of interest and importance to terrestrial and space applications. They are being considered for cooling of military combat vehicles and for spinning spacecraft. In order to investigate the effect of an accelerating force on LHP operation, a miniature LHP was installed on a spin table. Variable accelerating forces were imposed on the LHP by spinning the table at different angular speeds. Several patterns of accelerating forces were applied, i.e. continuous spin at different speeds and periodic spin at different speeds and frequencies. The resulting accelerations ranged from 1.17 g's to 4.7 g's. This paper presents the first part of the experimental study, i.e. the effects of a centrifugal force on the LHP start-up. Tests were conducted by varying the heat load to the evaporator, sink temperature, magnitude and frequency of centrifugal force, and LHP orientation relative to the direction of the accelerating force. The accelerating force seems to have little effect on the loop start-up in terms of temperature overshoot and superheat at boiling incipience. Changes in these parameters seem to be stochastic with or without centrifugal accelerating forces. The LHP started successfully in all tests.

Description: The effect of accelerating forces on the performance of loop heat pipes (LHP) is of interest and importance to terrestrial and space applications. LHP's are being considered for cooling of military combat vehicles and for spinning spacecraft. In order to investigate the effect of an accelerating force on LHP operation, a miniature LHP was installed on a spin table. Variable accelerating forces were imposed on the LHP by spinning the table at different angular speeds. Several patterns of accelerating forces were applied, i.e. continuous spin at different speeds and periodic spin at different speeds and frequencies. The resulting accelerations ranged from 1.17 g's to 4.7 g's. This paper presents the second part of the experimental study, i.e. the effect of an accelerating force on the LHP operating temperature. It has been known that in stationary tests the LHP operating temperature is a function of the evaporator power and the condenser sink temperature when the compensation temperature is not actively controlled. Results of this test program indicate that any change in the accelerating force will result in a chance in the LHP operating temperature through its influence on the fluid distribution in the evaporator, condenser and compensation chamber. However, the effect is not universal, rather it is a function of other test conditions. A steady, constant acceleration may result in an increase or decrease of the operating temperature, while a periodic spin will lead to a quasi-steady operating temperature over a sufficient time interval. In addition, an accelerating force may lead to temperature hysteresis and changes in the temperature oscillation. In spite of all these effects, the LHP continued to operate without any problems in all tests.

Description: A wind tunnel experiment was conducted in the NASA Langley Research Center (LaRC) Unitary Plan Wind Tunnel (UPW7) to determine the effects of passive surface porosity on vortex flow interactions about a general research fighter configuration at supersonic speeds. Optical flow measurement and flow visualization techniques were used and included pressure-sensitive paint (PSP), schlieren, and laser vapor screen (LVS). These techniques were combined with force and moment and conventional electronically-scanned pressure (ESP) measurements to quantify and to visualize the effects flow-through porosity applied to a wing leading-edge extension (LEX) mounted to a 65 deg cropped delta wing model.

Description: In this paper, parallel processing is used to analyze the mixing, and combustion behavior of hypersonic flow. Preliminary work for a sonic transverse hydrogen jet injected from a slot into a Mach 4 airstream in a two-dimensional duct combustor has been completed [Moon and Chung, 1996]. Our aim is to extend this work to three-dimensional domain using multithreaded domain decomposition parallel processing based on the flowfield-dependent variation theory. Numerical simulations of chemically reacting flows are difficult because of the strong interactions between the turbulent hydrodynamic and chemical processes. The algorithm must provide an accurate representation of the flowfield, since unphysical flowfield calculations will lead to the faulty loss or creation of species mass fraction, or even premature ignition, which in turn alters the flowfield information. Another difficulty arises from the disparity in time scales between the flowfield and chemical reactions, which may require the use of finite rate chemistry. The situations are more complex when there is a disparity in length scales involved in turbulence. In order to cope with these complicated physical phenomena, it is our plan to utilize the flowfield-dependent variation theory mentioned above, facilitated by large eddy simulation. Undoubtedly, the proposed computation requires the most sophisticated computational strategies. The multithreaded domain decomposition parallel processing will be necessary in order to reduce both computational time and storage. Without special treatments involved in computer engineering, our attempt to analyze the airbreathing combustion appears to be difficult, if not impossible.

Description: Computational Fluid Dynamics (CFD) analysis results are compared with experimental data from the Pennsylvania State University's (PSU) Propulsion Engineering Research Center (PERC) rocket based combined cycle (RBCC) rocket-ejector experiments. The PERC RBCC experimental hardware was in a direct-connect configuration in diffusion and afterburning (DAB) operation. The objective of the present work was to validate the Finite Difference Navier Stokes (FDNS) CFD code for the rocket-ejector mode internal fluid mechanics and combustion phenomena. A second objective was determine the best application procedures to use FDNS as a predictive/engineering tool. Three-dimensional CFD analysis was performed. Solution methodology and grid requirements are discussed. CFD results are compared to experimental data for static pressure, Raman Spectroscopy species distribution data and RBCC net thrust and specified impulse.

Description: Recent experiences in utilizing the global optimization methodology, based on polynomial and neural network techniques for fluid machinery design are summarized. Global optimization methods can utilize the information collected from various sources and by different tools. These methods offer multi-criterion optimization, handle the existence of multiple design points and trade-offs via insight into the entire design space can easily perform tasks in parallel, and are often effective in filtering the noise intrinsic to numerical and experimental data. Another advantage is that these methods do not need to calculate the sensitivity of each design variable locally. However, a successful application of the global optimization method needs to address issues related to data requirements with an increase in the number of design variables and methods for predicting the model performance. Examples of applications selected from rocket propulsion components including a supersonic turbine and an injector element and a turbulent flow diffuser are used to illustrate the usefulness of the global optimization method.

Description: The air density fluctuations in the plumes of fully-expanded, unheated free jets were investigated experimentally using a Rayleigh scattering based technique. The point measuring technique used a continuous wave laser, fiber-optic transmission and photon counting electronics. The radial and centerline profiles of time-averaged density and root-mean-square density fluctuation provided a comparative description of jet growth. To measure density fluctuation spectra a two-Photomultiplier tube technique was used. Crosscorrelation between the two PMT signals significantly reduced electronic shot noise contribution. Turbulent density fluctuations occurring up to a Strouhal number (Sr) of 2.5 were resolved. A remarkable feature of density spectra, obtained from the same locations of jets in 0.5〈 M〈1.5 range, is a constant Strouhal frequency for peak fluctuations. A detailed survey at Mach numbers M = 0.95, 1.4 and 1.8 showed that, in general, distribution of various Strouhal frequency fluctuations remained similar for the three jets. In spite of the similarity in the flow fluctuation the noise characteristics were found to be significantly different. Spark schlieren photographs and near field microphone measurements confirmed that the eddy Mach wave radiation was present in Mach 1.8 jet, and was absent in Mach 0.95 jet. To measure correlation between the flow and the far field sound pressure fluctuations, a microphone was kept at a distance of 50 diameters, 30 deg. to the flow direction, and the laser probe volume was moved from point to point in the flow. The density fluctuations in the peripheral shear layer of Mach 1.8 jet showed significant correlation up to the measurement limit of Sr = 2.5, while for Mach 0.95 jet no correlation was measured. Along the centerline measurable correlation was found from the end of the potential core and at the low frequency range (Sr less than 0.5). Usually the normalized correlation values increased with an increase of the jet Mach number. The experimental data point out eddy Mach waves as a strong source of sound generation in supersonic jets and fail to locate the primary noise mechanism in subsonic jets.

Description: This paper presents two means for enhancing nucleate boiling and critical heat flux under microgravity conditions: using micro-configured metal-graphite composites as the boiling surface and dilute aqueous solutions of long-chain alcohols as the working fluid. In the former, thermocapillary force induced by temperature difference between the graphite-fiber tips and the metal matrix plays an important role in bubble detachment. Thus boiling-heat transfer performance does not deteriorate in a reduced-gravity environment. In the latter cases, the surface tension-temperature gradient of the long-chain alcohol solutions turns positive as the temperature exceeds a certain value. Consequently, the Marangoni effect does not impede, but rather aids in bubble departure from the heating surface. This feature is most favorable in microgravity. As a result, the bubble size of departure is substantially reduced at higher frequencies. Based on the existing experimental data, and a two-tier theoretical model, correlation formulas are derived for nucleate boiling on the copper-graphite and aluminum-graphite composite surfaces, in both the isolated and coalesced bubble regimes. In addition, performance equations for nucleate boiling and critical heat flux in dilute aqueous solutions of long-chain alcohols are obtained.

Description: The second-order factorizable discretization of the compressible Euler equations developed by Sidilkover is extended to conservation form on general curvilinear body-fitted grids. The discrete equations are solved by symmetric collective Gauss-Seidel relaxation and FAS multigrid. Solutions for flow in a channel with Mach numbers ranging from 0.0001 to a supercritical Mach number are shown, demonstrating uniform convergence rates and no loss of accuracy in the incompressible limit. A solution for the flow around the leading edge of a semi-infinite parabolic body demonstrates that the scheme maintains rapid convergence for a flow containing a stagnation point.

Description: Detailed velocity measurements were made along a flat plate subject to the same dimensionless pressure gradient as the suction side of a modern low-pressure turbine airfoil. Reynolds numbers based on wetted plate length and nominal exit velocity were varied from 50,000 to 300,000, covering cruise to takeoff conditions. Low and high inlet free-stream turbulence intensities (0.2% and 7%) were set using passive grids. The location of boundary-layer separation does not depend strongly on the free-stream turbulence level or Reynolds number, as long as the boundary layer remains non-turbulent prior to separation. Strong acceleration prevents transition on the upstream part of the plate in all cases. Both free-stream turbulence and Reynolds number have strong effects on transition in the adverse pressure gradient region. Under low free-stream turbulence conditions transition is induced by instability waves in the shear layer of the separation bubble. Reattachment generally occurs at the transition start. At Re = 50,000 the separation bubble does not close before the trailing edge of the modeled airfoil. At higher Re, transition moves upstream, and the boundary layer reattaches. With high free-stream turbulence levels, transition appears to occur in a bypass mode, similar to that in attached boundary layers. Transition moves upstream, resulting in shorter separation regions. At Re above 200,000, transition begins before separation. Mean velocity, turbulence and intermittency profiles are presented.

Description: Jet spreading enhancement with a certain coannular nozzle configuration has been explored. When the outer nozzle is flaired (i.e., made convergent-divergent) the ensuing jet spreads faster than the case where the outer nozzle is convergent. The spreading enhancement is most pronounced when the outer flow is run near 'transonic' condition, in an overexpanded state. Under this condition, the increased spreading takes place regardless of the operating conditions of the inner jet. This observation, first made in a small scale facility, has been confirmed and studied in some detail in a larger-scale facility. Results of the latter experiment are presented in this paper. The spreading increase is shown to be substantial and comparable to or better than that achieved by a lobed nozzle. Estimates based on idealized flow indicate that there is an accompanying thrust penalty - the actual penalty is expected to be less than the estimate but remains undetermined at this time. In both the earlier and the present experiments, the spreading increase has often been found to accompany a flow resonance. The nature of this resonance is addressed in this paper. It is shown that the spreading increase takes place even if the resonance is absent. Thus, flow excitation due to the resonance is ruled out as the underlying mechanism. While the complete mechanism remains unclear, it is conjectured that pressure gradients near the nozzle, characteristic of overexpanded flow, are at the root of the phenomenon.

Description: Pressure gradients act differently on liquid particles and suspended bubbles and are, therefore, capable of inducing a relative motion between the phases even when no relative velocity initially exists. As a consequence of the enhanced heat transfer in the presence of convection, this fact may have a major impact on the evolution of a vapor bubble. The effect is particularly strong in the case of a collapsing bubble for which, due to the conservation of the system's impulse, the induced relative velocity tends to be magnified when the bubble volume shrinks. A practical application could be, for instance, the enhancement of the condensation rate of bubbles downstream of a heated region, thereby reducing the quality of a flowing liquid-vapor mixture. A simple model of the process, in which the bubble is assumed to be spherical and the flow potential, is developed in the paper.

Description: A combined experimental and computational study has been performed to investigate the detailed heat transfer coefficient distributions within a complex blade trailing edge passage. The experimental measurements are made using a steady liquid crystal thermography technique applied to one major side of the passage. The geometry of the trailing edge passage is that of a two-pass serpentine circuit with a sharp 180-degree turning region at the tip. The upflow channel is split by interrupted ribs into two major subchannels, one of which is turbulated. This channel has an average aspect ratio of roughly 14:1. The spanwise extent of the channel geometry includes both area convergence from root to tip, as well as taper towards the trailing edge apex. The average section Reynolds numbers tested in this upflow channel range from 55,000 to 98,000. The tip section contains a turning vane near the extreme comer. The downflow channel has an aspect ratio of about 5:1, and also includes convergence and taper. Turbulators of varying sizes are included in this channel also. Both detailed heat transfer and pressure distribution measurements are presented. The pressure measurements are incorporated into a flow network model illustrating the major loss contributors.

Description: This paper discusses a new technology for variable emissivity (vari-e) radiator surfaces, which has significant advantages over traditional radiators and promises an alternative design technique for future spacecraft thermal control systems. All spacecraft rely on radiative surfaces to dissipate waste heat. These radiators have special coatings, typically with a low solar absorptivity and a high infrared-red emissivity, that are intended to optimize performance under the expected heat load and thermal sink environment. The dynamics of the heat loads and thermal environment make it a challenge to properly size the radiator and often require some means of regulating the heat rejection rate of the radiators in order to achieve proper thermal balance. Specialized thermal control coatings, which can passively or actively adjust their emissivity offer an attractive solution to these design challenges. Such systems would allow intelligent control of the rate of heat loss from a radiator in response to heat load and thermal environmental variations. Intelligent thermal control through variable emissivity systems is well suited for nano and pico spacecraft applications where large thermal fluctuations are expected due to the small thermal mass and limited electric resources. Presently there are three different types of vari-e technologies under development: Micro ElectroMechanical Systems (MEMS) louvers, Electrochromic devices, and Electrophoretic devices. This paper will describe several prototypes of micromachined (MEMS) louvers and experimental results for the emissivity variations measured on theses prototypes. It will further discuss possible actuation mechanisms and space reliability aspects for different designs. Finally, for comparison parametric evaluations of the thermal performances of the new vari-e technology and standard thermal control systems are presented in this paper.

Description: Observations regarding rotating cavitation and cavitation surge experienced during the development of the Fastrac engine turbopump are discussed. Detailed observations acquired from the analysis of both water flow and liquid oxygen test data are offered in this paper. Scaling and general comparison of rotating cavitation between water flow and liquid oxygen testing are discussed. Complex data features linking the localized rotating cavitation mechanism of the inducer to system surge components are described in detail. Finally a description of a lumped-parameter hydraulic system model developed to better understand observed data is given.

Description: A method to obtain coupled Computational Fluid Dynamics-Direct Simulation Monte Carlo (CFD-DSMC), 3-D flow field solutions for highly blunt bodies at low incidence is presented and applied to one concept of the Mars Sample Return Orbiter vehicle as a demonstration of the technique. CFD is used to solve the high-density blunt forebody flow defining an inflow boundary condition for a DSMC solution of the afterbody wake flow. By combining the two techniques in flow regions where most applicable, the entire mixed flow field is modeled in an appropriate manner.

Description: This paper advances the state-of-the-art in spray computations with some of our recent contributions involving scalar Monte Carlo PDF (Probability Density Function), unstructured grids and parallel computing. It provides a complete overview of the scalar Monte Carlo PDF and Lagrangian spray computer codes developed for application with unstructured grids and parallel computing. Detailed comparisons for the case of a reacting non-swirling spray clearly highlight the important role that chemistry/turbulence interactions play in the modeling of reacting sprays. The results from the PDF and non-PDF methods were found to be markedly different and the PDF solution is closer to the reported experimental data. The PDF computations predict that some of the combustion occurs in a predominantly premixed-flame environment and the rest in a predominantly diffusion-flame environment. However, the non-PDF solution predicts wrongly for the combustion to occur in a vaporization-controlled regime. Near the premixed flame, the Monte Carlo particle temperature distribution shows two distinct peaks: one centered around the flame temperature and the other around the surrounding-gas temperature. Near the diffusion flame, the Monte Carlo particle temperature distribution shows a single peak. In both cases, the computed PDF's shape and strength are found to vary substantially depending upon the proximity to the flame surface. The results bring to the fore some of the deficiencies associated with the use of assumed-shape PDF methods in spray computations. Finally, we end the paper by demonstrating the computational viability of the present solution procedure for its use in 3D combustor calculations by summarizing the results of a 3D test case with periodic boundary conditions. For the 3D case, the parallel performance of all the three solvers (CFD, PDF, and spray) has been found to be good when the computations were performed on a 24-processor SGI Origin work-station.

Description: A series of heated tube experiments was performed to investigate fluid instabilities that occur during heating of supercritical fluids. In these tests, JP-7 flowed vertically through small diameter tubes at supercritical pressures. Test section heated length, diameter, mass flow rate, inlet temperature, and heat flux were varied in an effort to determine the range of conditions that trigger the instabilities. Heat flux was varied up to 4 BTU/sq in./s, and test section wall temperatures reached as high as 1950 F. A statistical model was generated to explain the trends and effects of the control variables. The model included no direct linear effect of heat flux on the occurrence of the instabilities. All terms involving inlet temperature were negative, and all terms involving mass flow rate were positive. Multiple tests at conditions that produced instabilities provided inconsistent results. These inconsistencies limit the use of the model as a predictive tool. Physical variables that had been previously postulated to control the onset of the instabilities, such as film temperature, velocity, buoyancy, and wall-to-bulk temperature ratio, were evaluated here. Film temperatures at or near critical occurred during both stable and unstable tests. All tests at the highest velocity were stable, but there was no functional relationship found between the instabilities and velocity, or a combination of velocity and temperature ratio. Finally, all of the unstable tests had significant buoyancy at the inlet of the test section, but many stable tests also had significant buoyancy forces.

Description: A study is conducted on high heat-flux pool boiling of pentane on micro-configured composite surfaces. The boiling surfaces are copper-graphite (Cu-Gr) and aluminum-graphite (Al-Gr) composites with a fiber volume concentration of 50%. The micro-graphite fibers embedded in the matrix contribute to a substantial enhancement in boiling heat-transfer performance. Correlation equations are obtained for both the isolated and coalesced bubble regimes, utilizing a mathematical model based on a metal-graphite, two-tier configuration with the aid of experimental data. A new model to predict the critical heat flux (CHF) on the composites is proposed to explain the fundamental aspects of the boiling phenomena. Three different factors affecting the CHF are considered in the model. Two of them are expected to become the main agents driving vapor volume detachment under microgravity conditions, using the metal-graphite composite surfaces as the heating surface and using liquids with an unusual Marangoni effect as the working fluid.

Description: An experimental investigation was conducted to examine acoustic receptivity and subsequent boundary-layer instability evolution for a Blasius boundary layer formed on a flat plate in the presence of two-dimensional (2-D) and oblique (3-D) surface waviness. The effect of the non-localized surface roughness geometry and acoustic wave amplitude on the receptivity process was explored. The surface roughness had a well defined wavenumber spectrum with fundamental wavenumber k (sub w). A planar downstream traveling acoustic wave was created to temporally excite the flow near the resonance frequency of an unstable eigenmode corresponding to k (sub ts) = k (sub w). The range of acoustic forcing levels, epsilon, and roughness heights, DELTA h, examined resulted in a linear dependence of receptivity coefficients; however, the larger values of the forcing combination epsilon dot DELTA h resulted in subsequent nonlinear development of the Tollmien-Schlichting (T-S) wave. This study provided the first experimental evidence of a marked increase in the receptivity coefficient with increasing obliqueness of the surface waviness in excellent agreement with theory. Detuning of the 2-D and oblique disturbances was investigated by varying the streamwise wall-roughness wavenumber a,, and measuring the T-S response. For the configuration where laminar-to-turbulent breakdown occurred, the breakdown process was found to be dominated by energy at the fundamental and harmonic frequencies, indicative of K-type breakdown.

Description: Experimental data from jet-engine tests have indicated that turbine efficiencies at takeoff can be as much as two points higher than those at cruise conditions. Recent studies have shown that Reynolds number effects contribute to the lower efficiencies at cruise conditions. In the current study numerical simulations have been performed to study the boundary layer development in a two-stage low-pressure turbine, and to evaluate the models available for low Reynolds number flows in turbomachinery. In a previous study using the same geometry the predicted time-averaged boundary layer quantities showed excellent agreement with the experimental data, but the predicted unsteady results showed only fair agreement with the experimental data. It was surmised that the blade count approximation used in the numerical simulations generated more unsteadiness than was observed in the experiments. In this study a more accurate blade approximation has been used to model the turbine, and the method of post-processing the boundary layer information has been modified to more closely resemble the process used in the experiments. The predicted results show improved agreement with the unsteady experimental data.

Description: This paper presents viewgraphs of the functionality of a loop heat pipe that was subjected to variable accelerating forces. The topics include: 1) Summary of LHP (Loop Heat Pipe) Design Parameters; 2) Picture of the LHP; 3) Schematic of Test Setup; 4) Test Configurations; 5) Test Profiles; 6) Overview of Test Results; 7) Start-up; 8) Typical Start-up without Temperature Overshoot; 9) Start-up with a Large Temperature Overshoot; 10) LHP Operation Under Stationary Condition; 11) LHP Operation Under Continuous Acceleration; 12) LHP Operation Under Periodic Acceleration; 13) Effects of Acceleration on Temperature Oscillation and Hysteresis; 14) Temperature Oscillation/Hysteresis vs Spin Rate; and 15) Summary.

Description: A review of some of the basic gas turbine technology being developed at the NASA John H. Glenn Research Center at Lewis Field, which may have the potential to be applied to ground-based systems, is presented in this paper. Only a sampling of the large number of research activities underway at the Glenn Research Center can be represented here. The items selected for presentation are those that may lead to increased power and efficiency, reduced cycle design time and cost, improved thermal design, reduced fatigue and fracture, reduced mechanical friction and increased operating margin. The topic of improved material will be presented in this conference and shall not be discussed here. The topics selected for presentation are key research activities at the Glenn Center of Excellence on Turbo-machinery. These activities should be of interest and utility to this ISABE (International Symposium on Air Breathing Engines) Special Forum on Aero-Derivative Land-Based Gas Turbines and to the power industry.

Description: Orbital experiments on the behavior of liquid in spacecraft are planned. The Sloshsat free-flyer is described. Preparation of the experiments, and later evaluation, are supported by models of varying complexity. The characteristics of the models are discussed. Particular attention is given to the momentum transfer between the liquid and the spacecraft, in connection with the liquid impact that may occur at the end of a reorientation maneuver of the spacecraft.

Description: Infrared thermography was used to obtain data on the state of the boundary layer of a natural laminar flow airfoil in supersonic flight. In addition to the laminar-to-turbulent transition boundary, the infrared camera was able to detect shock waves and present a time dependent view of the flow field. A time dependent heat transfer code was developed to predict temperature distributions on the test subject and any necessary surface treatment. A commercially available infrared camera was adapted for airborne use in this application. Readily available infrared technology has the capability to provide detailed visualization of various flow phenomena in subsonic to hypersonic flight regimes.

Description: A general purpose, one dimensional fluid flow code has been interfaced with the thermal analysis program SINDA/G. The flow code, GFSSP, is capable of analyzing steady state and transient flow in a complex network. The flow code is capable of modeling several physical phenomena including compressibility effects, phase changes, body forces (such as gravity and centrifugal) and mixture thermodynamics for multiple species. The addition of GFSSP to SINDA/G provides a significant improvement in convective heat transfer modeling for SINDA/G. The interface development was conducted in two phases. This paper describes the first (which allows for steady and quasi-steady - unsteady solid, steady fluid - conjugate heat transfer modeling). The second (full transient conjugate heat transfer modeling) phase of the interface development will be addressed in a later paper. Phase 1 development has been benchmarked to an analytical solution with excellent agreement. Additional test cases for each development phase demonstrate desired features of the interface. The results of the benchmark case, three additional test cases and a practical application are presented herein.

Description: For long-duration space missions, the life support and In-Situ Resource Utilization (ISRU) systems necessary to lower the mass and volume of consumables carried from Earth will require more sophisticated chemical processing technologies involving gas-liquid two-phase flows. This paper discusses some preliminary two-phase flow work in packed columns and generation of bubbly suspensions, two types of flow systems that can exist in a number of chemical processing devices. The experimental hardware for a co-current flow, packed column operated in two ground-based low gravity facilities (two-second drop tower and KC- 135 low-gravity aircraft) is described. The preliminary results of this experimental work are discussed. The flow regimes observed and the conditions under which these flow regimes occur are compared with the available co-current packed column experimental work performed in normal gravity. For bubbly suspensions, the experimental hardware for generation of uniformly sized bubbles in Couette flow in microgravity conditions is described. Experimental work was performed on a number of bubbler designs, and the capillary bubble tube was found to produce the most consistent size bubbles. Low air flow rates and low Couette flow produce consistent 2-3 mm bubbles, the size of interest for the "Behavior of Rapidly Sheared Bubbly Suspension" flight experiment. Finally the mass transfer implications of these two-phase flows is qualitatively discussed.

Description: Generation of a uniform monodisperse bubbly suspension in low gravity is a rather difficult task because bubbles do not detach as easily as on Earth. Under microgravity, the buoyancy force is not present to detach the bubbles as they are formed from the nozzles. One way to detach the bubbles is to establish a detaching force that helps their detachment from the orifice. The drag force, established by flowing a liquid in a cross or co-flow configuration with respect to the nozzle direction, provides this additional force and helps detach the bubbles as they are being formed. This paper is concerned with studying the generation of a bubbly suspension in low gravity in support of a flight definition experiment titled "Behavior of Rapidly Sheared Bubbly Suspension." Generation of a bubbly suspension, composed of 2 and 3 mm diameter bubbles with a standard deviation 〈10% of the bubble diameter, was identified as one of the most important engineering/science issues associated with the flight definition experiment. This paper summarizes the low gravity experiments that were conducted to explore various ways of making the suspension. Two approaches were investigated. The first was to generate the suspension via a chemical reaction between the continuous and dispersed phases using effervescent material, whereas the second considered the direct injection of air into the continuous phase. The results showed that the reaction method did not produce the desired bubble size distribution compared to the direct injection of bubbles. However, direct injection of air into the continuous phase (aqueous salt solution) resulted in uniform bubble-diameter distribution with acceptable bubble-diameter standard deviation.

Description: Detailed heat transfer measurements and predictions are given for a power generation turbine rotor with 129 deg of nominal turning and an axial chord of 137 mm. Data were obtained for a set of four exit Reynolds numbers comprised of the design point of 628,000, -20%, +20%, and +40%. Three ideal exit pressure ratios were examined including the design point of 1.378, -10%, and +10%. Inlet incidence angles of 0 deg and +/-2 deg were also examined. Measurements were made in a linear cascade with highly three-dimensional blade passage flows that resulted from the high flow turning and thick inlet boundary layers. Inlet turbulence was generated with a blown square bar grid. The purpose of the work is the extension of three-dimensional predictive modeling capability for airfoil external heat transfer to engine specific conditions including blade shape, Reynolds numbers, and Mach numbers. Data were obtained by a steady-state technique using a thin-foil heater wrapped around a low thermal conductivity blade. Surface temperatures were measured using calibrated liquid crystals. The results show the effects of strong secondary vortical flows, laminar-to-turbulent transition, and also show good detail in the stagnation region.

Description: An experimental study was made to obtain heat transfer and air temperature data for a simple 3-leg serpentine test section that simulates a turbine blade internal cooling passage with trip strips and bleed holes. The objectives were to investigate the interaction of ribs and various bleed conditions on internal cooling and to gain a better understanding of bulk air temperature in an internal passage. Steady state heat transfer measurements were obtained using a transient technique with thermochromic liquid crystals. Trip strips were attached to one wall of the test section and were located either between or near the bleed holes. The bleed holes, used for film cooling, were metered to simulate the effect of external pressure on the turbine blade. Heat transfer enhancement was found to be greater for ribs near bleed holes compared to ribs between holes, and both configurations were affected slightly by bleed rates upstream. Air temperature measurements were taken at discreet locations along one leg of the model. Average bulk air temperatures were found to remain fairly constant along one leg of the model.

Description: This paper, the first of a two-part series, presents results of an unsteady rotordynamic analysis of a long-bearing squeeze film damper executing circular centered orbits using a fluid circuit approach. A series of nodes and branches represent the geometry of the flow circuit. The mass and momentum conservation equations are solved to predict the pressure distribution in the squeeze film. The motion of the bearing is simulated by the variation of geometry within the flow path. The modeling methodology is benchmarked against published experimental long-bearing squeeze film damper test results. The model provides good agreement with the experimental damping coefficient.

Description: Experimental and complementary CFD results from the study of the rocket-ejector mode of a Rocket Based Combined Cycle (RBCC) engine are presented and discussed. The experiments involved systematic flowfield measurements in a two-dimensional, variable geometry rocket-ejector system. The rocket-ejector system utilizes a single two-dimensional, gaseous oxygen/gaseous hydrogen rocket as the ejector. To gain a thorough understanding of the rocket-ejector's internal fluid mechanic/combustion phenomena, experiments were conducted with both direct-connect and sea-level static configurations for a range of rocket operating conditions. Overall system performance was obtained through global measurements of wall static pressure profiles, heat flux profiles and engine thrust, whereas detailed mixing and combustion information was obtained through Raman spectroscopy measurements of major species (oxygen, hydrogen, nitrogen and water vapor). The experimental results for both the direct-connect and sea-level static configurations are compared with CFD predictions of the flowfield.

Description: This pper provides a validation summary of the spray computations performed as a part of the NCC (National Combustion Code) development activity. NCC is being developed with the aim of advancing the current prediction tools used in the design of advanced technology combustors based on the multidimensional computational methods. The solution procedure combines the novelty of the application of the scalar Monte Carlo PDF (Probability Density Function) method to the modeling of turbulent spray flames with the ability to perform the computations on unstructured grids with parallel computing. The calculation procedure was applied to predict the flow properties of three different spray cases. One is a nonswirling unconfined reacting spray, the second is a nonswirling unconfined nonreacting spray, and the third is a confined swirl-stabilized spray flame. The comparisons involving both gas-phase and droplet velocities, droplet size distributions, and gas-phase temperatures show reasonable agreement with the available experimental data. The comparisons involve both the results obtained from the use of the Monte Carlo PDF method as well as those obtained from the conventional computational fluid dynamics (CFD) solution. Detailed comparisons in the case of a reacting nonswirling spray clearly highlight the importance of chemistry/turbulence interactions in the modeling of reacting sprays. The results from the PDF and non-PDF methods were found to be markedly different and the PDF solution is closer to the reported experimental data. The PDF computations predict that most of the combustion occurs in a predominantly diffusion-flame environment. However, the non-PDF solution predicts incorrectly that the combustion occurs in a predominantly vaporization-controlled regime. The Monte Carlo temperature distribution shows that the functional form of the PDF for the temperature fluctuations varies substantially from point to point. The results also bring to the fore some of the deficiencies associated with the use of assumed-shape PDF methods in spray computations.

Description: The proposed paper compares predicted turbine vane heat transfer for a rough surface over a wide range of test conditions with experimental data. Predictions were made for the entire vane surface. However, measurements were made only over the suction surface of the vane, and the leading edge region of the pressure surface. Comparisons are shown for a wide range of test conditions. Inlet pressures varied between 3 and 15 psia, and exit Mach numbers ranged between 0.3 and 0.9. Thus, while a single roughened vane was used for the tests, the effective rougness,(k(sup +)), varied by more than a factor of ten. Results were obtained for freestream turbulence levels of 1 and 10%. Heat transfer predictions were obtained using the Navier-Stokes computer code RVCQ3D. Two turbulence models, suitable for rough surface analysis, are incorporated in this code. The Cebeci-Chang roughness model is part of the algebraic turbulence model. The k-omega turbulence model accounts for the effect of roughness in the application of the boundary condition. Roughness causes turbulent flow over the vane surface. Even after accounting for transition, surface roughness significantly increased heat transfer compared to a smooth surface. The k-omega results agreed better with the data than the Cebeci-Chang model. However, the low Reynolds number k-omega model did not accurately account for roughness when the freestream turbulence level was low. The high Reynolds number version of this model was more suitable when the freestream turbulence was low.

Description: This paper details the design and logic of an experimental investigation to study axial reorientation in low gravity. The Sloshsat free-flyer is described. The planned axial reorientation experiments and test matrixes are presented. Existing analytical tools are discussed. Estimates for settling range from 64 to 1127 seconds. The planned experiments are modelled using computational fluid dynamics. These models show promise in reducing settling estimates and demonstrate the ability of pulsed high thrust settling to emulate lower thrust continuous firing.

Description: In the setting of the conservation element-solution element (CE/SE) method, a new and unified wall boundary treatment for the Navier-Stokes and Euler Equations is proposed. In essence, the shear stress exerted on the fluid by a wall is modeled as a source term as a part of local spacetime flux conservation in the vicinity of a wall boundary. When the fluid is inviscid, the source term vanishes and the boundary condition reduces to the usual 'slip' condition. On the other hand, when the fluid is viscous, the source-term effect is consistent with the traditional no-slip condition. Numerical results show that the new treatment is robust, efficient, and accurate for viscous and inviscid flows.

Description: Turbine vane heat transfer distributions obtained using an infrared camera technique are described. Infrared thermography was used because noncontact surface temperature measurements were desired. Surface temperatures were 80 C or less. Tests were conducted in a three vane linear cascade, with inlet pressures between 0.14 and 1.02 atm., and exit Mach numbers of 0.3, 0.7, and 0.9, for turbulence intensities of approximately 1 and 10%. Measurements were taken on the vane suction side, and on the pressure side leading edge region. The designs for both the vane and test facility are discussed. The approach used to account for conduction within the vane is described. Midspan heat transfer distributions are given for the range of test conditions.

Description: A blackbody calibration furnace at the NASA Dryden Flight Research Center is used to calibrate heat flux gages. These gages are for measuring the aerodynamic heat flux on hypersonic flight vehicle surfaces. The blackbody is a graphite tube with a midplane partition which divides the tube into two compartments (dual cavities). Electrical resistance heating is used to heat the graphite tube. This heating and the boundary conditions imposed on the graphite tube result in temperature gradients along the walls of the blackbody cavity. This paper describes measurements made during steady-state operation and development of finite-difference thermal models of the blockbody furnace at 1100 C. Two configurations were studied, one with the blackbody outer surface insulated and the other without insulation. The dominant modes of heat transfer were identified for each configuration and the effect of variations in material properties and electric current that was passed through the blackbody were quantified.

Description: This paper presents recent work on developing methods for analyzing radiation heat transfer between diffuse-gray surfaces using p-version finite elements. The work was motivated by a thermal analysis of a High Speed Civil Transport (HSCT) wing structure which showed the importance of radiation heat transfer throughout the structure. The analysis also showed that refining the finite element mesh to accurately capture the temperature distribution on the internal structure led to very large meshes with unacceptably long execution times. Traditional methods for calculating surface-to-surface radiation are based on assumptions that are not appropriate for p-version finite elements. Two methods for determining internal radiation heat transfer are developed for one and two-dimensional p-version finite elements. In the first method, higher-order elements are divided into a number of sub-elements. Traditional methods are used to determine radiation heat flux along each sub-element and then mapped back to the parent element. In the second method, the radiation heat transfer equations are numerically integrated over the higher-order element. Comparisons with analytical solutions show that the integration scheme is generally more accurate than the sub-element method. Comparison to results from traditional finite elements shows that significant reduction in the number of elements in the mesh is possible using higher-order (p-version) finite elements.

Description: The ability of three types of turbulence models to accurately predict the effects of curvature on the flow in a U-duct is studied. An explicit algebraic stress model performs slightly better than one- or two-equation linear eddy viscosity models, although it is necessary to fully account for the variation of the production-to-dissipation-rate ratio in the algebraic stress model formulation. In their original formulations, none of these turbulence models fully captures the suppressed turbulence near the convex wall, whereas a full Reynolds stress model does. Some of the underlying assumptions used in the development of algebraic stress models are investigated and compared with the computed flowfield from the full Reynolds stress model. Through this analysis, the assumption of Reynolds stress anisotropy equilibrium used in the algebraic stress model formulation is found to be incorrect in regions of strong curvature. By the accounting for the local variation of the principal axes of the strain rate tensor, the explicit algebraic stress model correctly predicts the suppressed turbulence in the outer part of the boundary layer near the convex wall.

Description: Scientific visualizations of two-dimensional compressible flow of a gas with discontinuities are presented. The numerical analogue to experimental techniques such as schlieren imaging, shadowgraphs, and interferograms are discussed. Edge detection techniques are utilized to identify the discontinuities. In particular, the zero crossing of the Laplacian of a field (usually density) is recommended for extracting the discontinuities. An algorithm to extract and quantify the discontinuities is presented. To illustrate the methods developed in the report, the example chosen is that of an unsteady interaction of a shock wave with a contact discontinuity.

Description: The Fluid Mechanics and Acoustics Laboratory (FM&AL) was established At Hampton University in June of 1996. In addition, the FM&AL jointly conducted research with the Central AeroHydrodynamics Institute (TsAGI, Moscow) in Russia under a 2.5 year Civilian Research and Development Foundation (CRDF). The goals of the FM&AL programs are two fold: 1) to improve the working efficiency of the FM&AL team in generating new innovative ideas and in conducting research in the field of fluid dynamics and acoustics, basically for improvement of supersonic and subsonic aircraft engines, and 2) to attract promising minority students to this research and training and, in cooperation with other HU departments, to teach them basic knowledge in Aerodynamics, Gas Dynamics, and Theoretical and Experimental Methods in Aeroacoustics and Computational Fluid Dynamics (CFD). The research at the FM&AL supports reduction schemes associated with the emission of engine pollutants for commercial aircraft and concepts for reduction of observables for military aircraft. These research endeavors relate to the goals of the NASA Strategic Enterprise in Aeronautics concerning the development of environmentally acceptable aircraft. It is in this precise area, where the US aircraft industry, academia, and Government are in great need of trained professionals and which is a high priority goal of the Minority University Research and Education (MUREP) Program, that the HU FM&AL can make its most important contribution. This project already benefits NASA and HU because: First, the innovation, testing, and further development of new techniques for advanced propulsion systems are necessary for the successful attainment of the NASA Long Term Goals in Aeronautics and Space Transportation Technology (ASTT) including Global Civil Aviation, Revolutionary Technology Leaps, Access to Space, R&D Services, and the economic competitiveness of the US Aircraft Industry in the 2 1 st century. Secondly, the joint theoretical and experimental research and training by the GRC-HU Teams aids: using advanced methods and experience in Aerospace Engineering for domestic industries and training of HU students for interesting innovative work in the numerical simulation field as well as engineering and experimental research. HU students use and modify existing numerical codes for the solution of actual applied problems of the NASA Langley Research Center (LaRC)

Description: A review of the research conducted at the National Aeronautics and Space Administration (NASA), Langley Research Center (LaRC) into high-speed vortex flows during the 1970s, 1980s, and 1990s is presented. The data reviewed is for flat plates, cavities, bodies, missiles, wings, and aircraft. These data are presented and discussed relative to the design of future vehicles. Also presented is a brief historical review of the extensive body of high-speed vortex flow research from the 1940s to the present in order to provide perspective of the NASA LaRC's high-speed research results. Data are presented which show the types of vortex structures which occur at supersonic speeds and the impact of these flow structures to vehicle performance and control is discussed. The data presented shows the presence of both small- and large scale vortex structures for a variety of vehicles, from missiles to transports. For cavities, the data show very complex multiple vortex structures exist at all combinations of cavity depth to length ratios and Mach number. The data for missiles show the existence of very strong interference effects between body and/or fin vortices and the downstream fins. It was shown that these vortex flow interference effects could be both positive and negative. Data are shown which highlights the effect that leading-edge sweep, leading-edge bluntness, wing thickness, location of maximum thickness, and camber has on the aerodynamics of and flow over delta wings. The observed flow fields for delta wings (i.e. separation bubble, classical vortex, vortex with shock, etc.) are discussed in the context of' aircraft design. And data have been shown that indicate that aerodynamic performance improvements are available by considering vortex flows as a primary design feature. Finally a discussing of a design approach for wings which utilize vortex flows for improved aerodynamic performance at supersonic speed is presented.

Description: A new bias flow impedance model is developed for perforated plates from basic principles using as little empiricisms as possible. A quality experimental database was used to determine the predictive validity of the model. Results show that the model performs better for higher (15%) rather than lower (5%) percent open area (POA) samples. Based on the least squares ratio of numerical vs. experimental results, model predictions were on average within 20% and 30% for the higher and lower (POA), respectively. It is hypothesized on the work of other investigators that at lower POAs the higher fluid velocities in the perforate's orifices start forming unsteady vortices, which is not accounted for in our model. The numerical model, in general also underpredicts the experiments. It is theorized that the actual acoustic C(sub D) is lower than the measured raylometer C(sub D) used in the model. Using a larger C(sub D) makes the numerical model predict lower impedances. The frequency domain model derived in this paper shows very good agreement with another model derived using a time domain approach.

Description: This paper, the second of a two-part series, presents results of an unsteady rotordynamic analysis of a long-bearing squeeze film damper executing orbits about an off center position using a fluid circuit approach. A series of nodes and branches represent the geometry of the flow circuit. The mass and momentum conservation equations are solved to predict the pressure distribution in the squeeze film. The motion of the bearing is simulated by the variation of geometry within the flow path. This effort represents the first modeling approach which allows for an arbitrary orbit size about an arbitrary position.

Description: An experimental investigation of the flow over the rear end of a 0.16 scale notchback automobile configuration has been conducted in the NASA Langley Basic Aerodynamics Research Tunnel (BART). The objective of this work was to investigate the flow separation that occurs behind the backlight and obtain experimental data that can be used to understand the physics and time-averaged structure of the flow field. A three-component laser velocimeter was used to make non-intrusive, velocity measurements in the center plane and in a single cross-flow plane over the decklid. In addition to off-body measurements, flow conditions on the car surface were documented via surface flow visualization, boundary layer measurements, and surface pressures. The experimental data show several features previously identified by other researchers, but also reveal differences between the flow field associated with this particular configuration and the generally accepted models for the flow over a notchback rear end.

Description: This viewgraph presentation gives an overview of the future of computational fluid dynamics (CFD), which in the past has pioneered the field of flow simulation. Over time CFD has progressed as computing power. Numerical methods have been advanced as CPU and memory capacity increases. Complex configurations are routinely computed now and direct numerical simulations (DNS) and large eddy simulations (LES) are used to study turbulence. As the computing resources changed to parallel and distributed platforms, computer science aspects such as scalability (algorithmic and implementation) and portability and transparent codings have advanced. Examples of potential future (or current) challenges include risk assessment, limitations of the heuristic model, and the development of CFD and information technology (IT) tools.