ALBERT

Description: This volume and its accompanying CD-ROM contain materials presented at the Minnowbrook III-2000 Workshop on Boundary Layer Transition and Unsteady Aspects of Turbomachinery Flows held at the Syracuse University Minnowbrook Conference Center, Blue Mountain Lake, New York, August 20-23, 2000. Workshop organizers were John E. LaGraff (Syracuse University), Terry V Jones (Oxford University), and J. Paul Gostelow (University of Leicester). The workshop followed the theme, venue, and informal format of two earlier workshops: Minnowbrook I (1993) and Minnowbrook II (1997). The workshop was focused on physical understanding the late stage (final breakdown) boundary layer transition, separation, and effects of unsteady wakes with the specific goal of contributing to engineering application of improving design codes for turbomachinery. The workshop participants included academic researchers from the USA and abroad, and representatives from the gas-turbine industry and government laboratories. The physical mechanisms discussed included turbulence disturbance environment in turbomachinery, flow instabilities, bypass and natural transition, turbulent spots and calmed regions, wake interactions with attached and separated boundary layers, turbulence and transition modeling and CFD, and DNS. This volume contains abstracts and copies of the viewgraphs presented, organized according to the workshop sessions. The viewgraphs are included on the CD-ROM only. The workshop summary and the plenary-discussion transcripts clearly highlight the need for continued vigorous research in the technologically important area of transition, separated and unsteady flows in turbomachines.

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

Description: This paper presents viewgraphs of physical sciences research priorities and plans at the Office of Biological and Physical Sciences Research (OBPR). The topics include: 1) Sixth Microgravity Fluid Physics and Transport Phenomena Conference; 2) Beneficial Characteristics of the Space Environment; 3) Windows of Opportunity for Research Derived from Microgravity; 4) Physical Sciences Research Program; 5) Fundamental Research: Space-based Results and Ground-based Applications; 6) Nonlinear Oscillations; and 7) Fundamental Research: Applications to Mission-Oriented Research.

Description: This paper presents viewgraphs of NASA's strategic and fundamental research program at the Office of Biological and Physical Research (OBPR). The topics include: 1) Colloid-Polymer Samples; 2) Pool Boiling Experiment; 3) The Dynamics of Miscible Interfaces: A Space Flight Experiment (MIDAS); and 4) ISS and Ground-based Facilities.

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: The stability of cylindrical liquid bridges in reduced gravity is affected by ambient vibrations of the spacecraft. Such vibrations are expected to excite capillary modes of the bridge. The lowest-order unstable mode is particularly susceptible to vibration as the length of the bridge approaches the stability limit. This low-order mode is known as the (2,0) mode and is an axisymmetric varicose mode of one wavelength in the axial direction. In this work, an optical system is used to detect the (2,0)-mode amplitude. The derivative of the error signal produced by this detector is used to produce the appropriate voltages on a pair of ring electrodes which are concentric with the bridge. A mode-coupled Maxwell stress profile is thus generated in proportional to the modal velocity. Depending on the sign of the gain, the damping of the capillary oscillation can be either increased or decreased. This effect has been demonstrated in Plateau-tank experiments. Increasing the damping of the capillary modes on free liquid surfaces in space could be beneficial for containerless processing and other novel technologies.

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

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

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

Description: 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 stability of cylindrical liquid bridges in reduced gravity is affected by ambient vibrations of the spacecraft. Such vibrations are expected to excite capillary modes of the bridge. The lowest-order unstable mode is particularly susceptible to vibration as the length of the bridge approaches the stability limit. This low-order mode is known as the (2,0) mode and is an axisymmetric varicose mode of one wavelength in the axial direction. In this work, an optical system is used to detect the (2,0)-mode amplitude. The derivative of the error signal produced by this detector is used to produce the appropriate voltages on a pair of ring electrodes which are concentric with the bridge. A mode-coupled Maxwell stress profile is thus generated in proportional to the modal velocity. Depending on the sign of the gain, the damping of the capillary oscillation can be either increased or decreased. This effect has been demonstrated in Plateau-tank experiments. Increasing the damping of the capillary modes on free liquid surfaces in space could be beneficial for containerless processing and other novel technologies. [work supported by NASA]

Description: Spray cooling has high potential in thermal management and life support systems by overcoming the deleterious effect of microgravity upon two-phase heat transfer. In particular spray cooling offers several advantages in heat flux removal that include the following: 1. By maintaining a wetted surface, spray droplets impinge upon a thin fluid film rather than a dry solid surface 2. Most heat transfer surfaces will not be smooth but rough. Roughness can enhance conductive cooling, aid liquid removal by flow channeling. 3. Spray momentum can be used to a) substitute for gravity delivering fluid to the surface, b) prevent local dryout and potential thermal runaway and c) facilitate liquid and vapor removal. Yet high momentum results in high We and Re numbers characterizing the individual spray droplets. Beyond an impingement threshold, droplets splash rather than spread. Heat flux declines and spray cooling efficiency can markedly decrease. Accordingly we are investigating droplet impingement upon a) dry solid surfaces, b) fluid films, c) rough surfaces and determining splashing thresholds and relationships for both dry surfaces and those covered by fluid films. We are presently developing engineering correlations delineating the boundary between splashing and non-splashing regions.

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: The focus of this project was to study the physical processes that govern tachocline dynamics and structure. Specific features explored included stratification, shear, waves, and toroidal and poloidal background fields. In order to address recent theoretical work on anisotropic mixing and dynamics in the tachocline, we were particularly interested in such anisotropic mixing for the specific tachocline processes studied. Transition to turbulence often shapes the largest-scale features that appear spontaneously in a flow during the development of turbulence. The resulting large-scale straining field can control the subsequent dynamics; therefore, anticipation of the large-scale straining field that results for individual realizations of the transition to turbulence can be important for subsequent dynamics, flow morphology, and transport characteristics. As a result, we paid particular attention to the development of turbulence in the stratified and sheared environment of the tachocline. This is complicated by the fact that the linearly stability of sheared MHD flows is non-self-adjoint, implying that normal asymptotic linear stability theory may not be relevant.

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

Description: 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 gravity-driven flow of non-neutrally buoyant suspensions is shown to be unstable to spanwise perturbations when the shearing motion generates a density profile that increases with height. The instability is simply due to having heavier material over light. The wavelength of the perturbation is found to be on the order of the thickness of the suspension layer. The parameters important to the problem are the angle of inclination of the layer relative to gravity, the relative density difference between the particles and fluid, the ratio of the particle size to the suspension layer, and the bulk volume fraction of particles. An example showing the growth rate as a function of wave number is shown.

Description: Results of an ongoing effort to quantify the role turbulence in scattering sound in jets are reported. Using a direct numerical simulation database to provide the flow data, ray paths traced through the mean flow are compared with those traced through the actual time evolving turbulent flow. Significant scattering by the turbulence is observed. The most notable effect is that upstream traveling waves that are trapped in the potential core by the mean flow, which acts as a wave guide, easily escape in the turbulent flow. A crude statistical estimate based on ray number density suggests that directivity is modified by the turbulence, but no rigorous treatment of non-uniformities in the high-frequency approximation is attempted.

Description: A numerical model of the tropical Atlantic ocean is used to investigate the upper layer pathways of the Meridional Overturning Circulation (MOC) in the tropical Atlantic. The main focus of this thesis is on those parts of the tropical circulation that are thought to be important for the MOC return flow, but whose dynamics have not been understood yet. It is shown how the particular structure of the tropical gyre and the MOO act to inhibit the flow of North Atlantic water into the equatorial thermocline. As a result, the upper layers of the tropical Atlantic are mainly fed by water from the South Atlantic. The processes that carry the South Atlantic water across the tropical Atlantic into the North Atlantic as part of the MOO are described here, and three processes that were hitherto not understood are explained as follows: The North Brazil Current rings are created as the result of the reflection of Rossby waves at the South American coast. These Rossby waves are generated by the barotropically unstable North Equatorial Countercurrent. The deep structure of the rings can be explained by merger of the wave's anticyclones with the deeper intermediate eddies that are generated as the intermediate western boundary current crosses the equator. The bands of strong zonal velocity in intermediate depths along the equator have hitherto been explained as intermediate currents. Here, an alternative interpretation of the observations is offered: The Eulerian mean flow along the equator is negligible and the observations are the signature of strong seasonal Rossby waves. The previous interpretation of the observations can then be explained as aliasing of the tropical wave field. The Tsuchyia Jets are driven by the Eliassen-Palm flux of the tropical instability waves. The equatorial current system with its strong shears is unstable and generates tropical instability waves.

Description: The magnetic Kelvin force has been proposed as an artificial gravity to control the orientation of paramagnetic liquid propellants such as liquid oxygen in a microgravity environment. This paper reports experiments performed in the NASA "Weightless Wonder" KC-135 aircraft, through the Reduced Gravity Student Flight Opportunities Program. The aircraft flies through a series of parabolic arcs providing about 25 s of microgravity in each arc. The experiment was conceived, designed, constructed, and performed by the undergraduate student team and their two faculty advisors. Two types of tanks were tested: square-base prismatic tanks 5 cm x 5 cm x 8.6 cm and circular cylinders 5 cm in diameter and 8.6 cm tall. The paramagnetic liquid was a 3.3 molar solution of MnCl2 in water. Tests were performed with each type of tank filled to depths of 1 cm and 4 cm. Each test compared a pair of tanks that were identical except that the base of one was a pole face of a 0.6 Tesla permanent magnet. The Kelvin force attracts paramagnetic materials toward regions of higher magnetic field. It was hypothesized that the Kelvin force would hold the liquid in the bottom of the tanks during the periods of microgravity. The tanks were installed in a housing that could slide on rails transverse to the flight direction. By manually shoving the housing, an identical impulse could be provided to each tank at the beginning of each period of microgravity. The resulting fluid motions were videotaped for later analysis.

Description: A pulse thermal loop heat transfer system includes a means to use pressure rises in a pair of evaporators to circulate a heat transfer fluid. The system includes one or more valves that iteratively, alternately couple the outlets the evaporators to the condenser. While flow proceeds from one of the evaporators to the condenser, heating creates a pressure rise in the other evaporator, which has its outlet blocked to prevent fluid from exiting the other evaporator. When the flow path is reconfigured to allow flow from the other evaporator to the condenser, the pressure in the other evaporator is used to circulate a pulse of fluid through the system. The reconfiguring of the flow path, by actuating or otherwise changing the configuration of the one or more valves, may be triggered when a predetermined pressure difference between the evaporators is reached.

Description: A NASA grant has been awarded to Cleveland State University (CSU) to develop a multi-dimensional (multi-D) Stirling computer code with the goals of improving loss predictions and identifying component areas for improvements. The University of Minnesota (UMN) and Gedeon Associates are teamed with CSU. Development of test rigs at UMN and CSU and validation of the code against test data are part of the effort. The one-dimensional (1-D) Stirling codes used for design and performance prediction do not rigorously model regions of the working space where abrupt changes in flow area occur (such as manifolds and other transitions between components). Certain hardware experiences have demonstrated large performance gains by varying manifolds and heat exchanger designs to improve flow distributions in the heat exchangers. 1-D codes were not able to predict these performance gains. An accurate multi-D code should improve understanding of the effects of area changes along the main flow axis, sensitivity of performance to slight changes in internal geometry, and, in general, the understanding of various internal thermodynamic losses. The commercial CFD-ACE code has been chosen for development of the multi-D code. This 2-D/3-D code has highly developed pre- and post-processors, and moving boundary capability. Preliminary attempts at validation of CFD-ACE models of MIT gas spring and "two space" test rigs were encouraging. Also, CSU's simulations of the UMN oscillating-flow fig compare well with flow visualization results from UMN. A complementary Department of Energy (DOE) Regenerator Research effort is aiding in development of regenerator matrix models that will be used in the multi-D Stirling code. This paper reports on the progress and challenges of this

Description: A liquid drop present on a solid surface can move because of a gradient in wettability along the surface, as manifested by a gradient in the contact angle. The contact angle at a given point on the contact line between a solid and a liquid in a gaseous medium is the angle between the tangent planes to the liquid and the solid surfaces at that point and is measured within the liquid side, by convention. The motion of the drop occurs in the direction of increasing wettability. The cause of the motion is the net force exerted on the drop by the solid surface because of the variation of the contact angle around the periphery. This force causes acceleration of an initially stationary drop, and leads to its motion in the direction of decreasing contact angle. The nature of the motion is determined by the balance between the motivating force and the resisting hydrodynamic force from the solid surface and the surrounding gaseous medium. A wettability gradient can be chemically induced as shown by Chaudhury and Whitesides who provided unambiguous experimental evidence that drops can move in such gradients. The phenomenon can be important in heat transfer applications in low gravity, such as when condensation occurs on a surface. Daniel et al have demonstrated that the velocity of a drop on a surface due to a wettability gradient in the presence of condensation can be more than two orders of magnitude larger than that observed in the absence of condensation. In the present research program, we have begun to study the motion of a drop in a wettability gradient systematically using a model system. Our initial efforts will be restricted to a system in which no condensation occurs. The experiments are performed as follows. First, a rectangular strip of approximate dimensions 10 x 20 mm is cut out of a silicon wafer. The strip is cleaned thoroughly and its surface is exposed to the vapor from an alkylchlorosilane for a period lasting between one and two minutes inside a desiccator. This is done using an approximate line source of the vapor in the form of a string soaked in the alkylchlorosilane. Ordinarily, many fluids, including water, wet the surface of silicon quite well. This means that the contact angle is small. But the silanized surface resists wetting, with contact angles that are as large as 100 degs. Therefore, a gradient of wettability is formed on the silicon surface. The region near the string is highly hydrophobic, and the contact angle decreases gradually toward a small value at the hydrophilic end away from this region. The change in wettability occurs over a distance of several mm. The strip is placed on a platform within a Plexiglas cell. Drops of a suitable liquid are introduced on top of the strip near the hydrophobic end. An optical system attached to a video camera is trained on the drop so that images of the moving drop can be captured on videotape for subsequent analysis. We have performed preliminary experiments with water as well as ethylene glycol drops. Results from these experiments will be presented in the poster. Future plans include the refinement of the experimental system so as to permit images to be recorded from the side as well as the top, and the conduct of a systematic study in which the drop size is varied over a good range. Experiments will be conducted with different fluids so as to obtain the largest possible range of suitably defined Reynolds and Capillary numbers. Also, an effort will be initiated on theoretical modeling of this motion. The challenges in the development of the theoretical description lie in the proper analysis of the region in the vicinity of the contact line, as well as in the free boundary nature of the problem. It is known that continuum models assuming the no slip condition all the way to the contact line fail by predicting that the stress on the solid surface becomes singular as the contact line is approached. One approach for dealing with this issue has been to relax the no-slip boundary condition using the Navier model. Molecular dynamics simulations of the contact line region show that for a non-polar liquid on a solid surface, the no-slip boundary condition is in fact incorrect near the contact line. Furthermore, the same simulations also show that the usual relationship between stress and the rate of deformation breaks down in the vicinity of the contact line. In developing continuum theoretical models of the system, we shall accommodate this knowledge to the extent possible.

Description: This research addresses turbulent gas flows laden with fine solid particles at sufficiently large mass loading that strong two-way coupling occurs. By two-way coupling we mean that the particle motion is governed largely by the flow, while the particles affect the gas-phase mean flow and the turbulence properties. Our main interest is in understanding how the particles affect the turbulence. Computational techniques have been developed which can accurately predict flows carrying particles that are much smaller than the smallest scales of turbulence. Also, advanced computational techniques and burgeoning computer resources make it feasible to fully resolve very large particles moving through turbulent flows. However, flows with particle diameters of the same order as the Kolmogorov scale of the turbulence are notoriously difficult to predict. Some simple flows show strong turbulence attenuation with reductions in the turbulent kinetic energy by up to a factor of five. On the other hand, some seemingly similar flows show almost no modification. No model has been proposed that allows prediction of when the strong attenuation will occur. Unfortunately, many technological and natural two-phase flows fall into this regime, so there is a strong need for new physical understanding and modeling capability. Our objective is to study the simplest possible turbulent particle-laden flow, namely homogeneous, isotropic turbulence with a uniform dispersion of monodisperse particles. We chose such a simple flow for two reasons. First, the simplicity allows us to probe the interaction in more detail and offers analytical simplicity in interpreting the results. Secondly, this flow can be addressed by numerical simulation, and many research groups are already working on calculating the flow. Our detailed data can help guide some of these efforts. By using microgravity, we can further simplify the flow to the case of no mean velocity for either the turbulence or the particles. In fact the addition of gravity as a variable parameter may help us to better understand the physics of turbulence attenuation. The experiments are conducted in a turbulence chamber capable of producing stationary or decaying isotropic turbulence with nearly zero mean flow and Taylor microscale Reynolds numbers up to nearly 500. The chamber is a 410 mm cubic box with the corners cut off to make it approximately spherical. Synthetic jet turbulence generators are mounted in each of the eight corners of the box. Each generator consists of a loudspeaker forcing a plenum and producing a pulsed jet through a 20 mm diameter orifice. These synthetic jets are directed into ejector tubes pointing towards the chamber center. The ejector tubes increase the jet mass flow and decrease the velocity. The jets then pass through a turbulence grid. Each of the eight loudspeakers is forced with a random phase and frequency. The resulting turbulence is highly Isotropic and matches typical behavior of grid turbulence. Measurements of both phases are acquired using particle image velocimetry (PIV). The gas is seeded with approximately 1 micron diameter seeding particles while the solid phase is typically 150 micron diameter spherical glass particles. A double-pulsed YAG laser and a Kodak ES-1.0 10-bit PIV camera provide the PIV images. Custom software is used to separate the images into individual images containing either gas-phase tracers or large particles. Modern high-resolution PIV algorithms are then used to calculate the velocity field. A large set of image pairs are acquired for each case, then the results are averaged both spatially and over the ensemble of acquired images. The entire apparatus is mounted in two racks which are carried aboard NASA's KC-135 Flying Microgravity Laboratory. The rack containing the turbulence chamber, the laser head, and the camera floats freely in the airplane cabin (constrained by competent NASA personnel) to minimize g-jitter.

Description: The present paper studies the numerical simulation of flows with shock/boundary-layer upstream interaction, under conditions of symmetry in geometry, boundary conditions, and grid. For this purpose, a series of two- and three-dimensional numerical test-cases were carried out. The tests showed that standard numerical schemes, which appear to be symmetry preserving under most flow configurations, produce nonsymmetric perturbations when large separated regions are present. These perturbations are amplified when the core flow is under compression. If the flow-blockage due to separation is sufficiently large, the symmetry of the flow may collapse altogether. Experimental evidence of this numerical behavior is also considered.

Description: Late-time dynamics and morphology of a stratified turbulent shear layer are examined using 1) Reynolds-stress and heat-flux budgets, 2) the single-point structure tensors introduced by Kassinos et al. (2001), and 3) flow visualization via 3D volume rendering. Flux reversal is observed during restratification in the edges of the turbulent layer. We present a first attempt to quantify the turbulence-mean-flow interaction and to characterize the predominant flow structures. Future work will extend this analysis to earlier times and different values of the Reynolds and Richardson numbers.

Description: Boiling is a complex phenomenon where hydrodynamics, heat transfer, mass transfer, and interfacial phenomena are tightly interwoven. An understanding of boiling and critical heat flux in microgravity environments is of importance to space based hardware and processes such as heat exchange, cryogenic fuel storage and transportation, electronic cooling, and material processing due to the large amounts of heat that can be removed with relatively little increase in temperature. Although research in this area has been performed in the past four decades, the mechanisms by which heat is removed from surfaces in microgravity are still unclear. In earth gravity, buoyancy is an important parameter that affects boiling heat transfer through the rate at which bubbles are removed from the surface. A simple model describing the bubble departure size based on a quasistatic force balance between buoyancy and surface tension is given by the Fritz [I] relation: Bo(exp 1/2) = 0.0208 theta where Bo is the ratio between buoyancy and surface tension forces. For small, rapidly growing bubbles, inertia associated with the induced liquid motion can also cause bubble departure. In microgravity, the magnitude of effects related to natural convection and buoyancy are small and physical mechanisms normally masked by natural convection in earth gravity such as Marangoni convection can substantially influence the boiling and vapor bubble dynamics. CHF (critical heat transfer) is also substantially affected by microgravity. In 1 g environments, Bo has been used as a correlating parameter for CHF. Zuber's CHF model for an infinite horizontal surface assumes that vapor columns formed by the merger of bubbles become unstable due to a Helmholtz instability blocking the supply of liquid to the surface. The jets are spaced lambda(sub D) apart, where lambda(sub D) = 2pi square root of 3[(sigma)/(g(rho(sub l) - rho(sub v)](exp 1/2) = 2pi square root of 3 L Bo(exp -1/2) = square root of 3 lambda(sub c) and is the wavelength that amplifies most rapidly. The critical wavelength, lambda(sub c), is the wavelength below which a vapor layer underneath a liquid layer is stable. For heaters with Bo smaller than about 3 (heaters smaller than lambda(sub D)), the above model is not applicable, and surface tension effects dominate. Bubble coalescence is thought to be the mechanism for CHF under these conditions. Small Bo can result by decreasing the size of a heater in earth gravity, or by operating a large heater in a lower gravity environment. In the microgravity of space, even large heaters can have low Bo, and models based on Helmholtz instability should not be applicable. The macrolayer model of Haramura and Katto is dimensionally equivalent to Zuber's model and has the same dependence on gravity, so it should not be applicable as well. The goal of this work is to determine how boiling heat transfer mechanisms in a low-g environment are altered from those at higher gravity levels. Boiling data using a microheater array was obtained under gravity environments ranging from 1.8 g to 0.02 g with heater sizes ranging from 2.7 mm to 1 mm. The boiling behavior for 2.7 mm at 0.02 g looked quite similar to boiling on the 1 mm heater at 1 g-the formation of a large primary bubble surrounded by smaller satellite bubbles was observed under both conditions. The similarity suggests that for heaters smaller than some fraction of I(sub c), coalescence and surface tension dominate boiling heat transfer. It also suggests that microgravity boiling can be studied by studying boiling on very small heaters.

Description: This work focuses on the properties of sheared granular materials near the jamming transition. The project currently involves two aspects. The first of these is an experiment that is a prototype for a planned ISS (International Space Station) flight. The second is discrete element simulations (DES) that can give insight into the behavior one might expect in a reduced-g environment. The experimental arrangement consists of an annular channel that contains the granular material. One surface, say the upper surface, rotates so as to shear the material contained in the annulus. The lower surface controls the mean density/mean stress on the sample through an actuator or other control system. A novel feature under development is the ability to 'thermalize' the layer, i.e. create a larger amount of random motion in the material, by using the actuating system to provide vibrations as well control the mean volume of the annulus. The stress states of the system are determined by transducers on the non-rotating wall. These measure both shear and normal components of the stress on different size scales. Here, the idea is to characterize the system as the density varies through values spanning dense almost solid to relatively mobile granular states. This transition regime encompasses the regime usually thought of as the glass transition, and/or the jamming transition. Motivation for this experiment springs from ideas of a granular glass transition, a related jamming transition, and from recent experiments. In particular, we note recent experiments carried out by our group to characterize this type of transition and also to demonstrate/ characterize fluctuations in slowly sheared systems. These experiments give key insights into what one might expect in near-zero g. In particular, they show that the compressibility of granular systems diverges at a transition or critical point. It is this divergence, coupled to gravity, that makes it extremely difficult if not impossible to characterize the transition region in an earth-bound experiment. In the DE modeling, we analyze dynamics of a sheared granular system in Couette geometry in two (2D) and three (3D) space dimensions. Here, the idea is to both better understand what we might encounter in a reduced-g environment, and at a deeper level to deduce the physics of sheared systems in a density regime that has not been addressed by past experiments or simulations. One aspect of the simulations addresses sheared 2D system in zero-g environment. For low volume fractions, the expected dynamics of this type of system is relatively well understood. However, as the volume fraction is increased, the system undergoes a phase transition, as explained above. The DES concentrate on the evolution of the system as the solid volume fraction is slowly increased, and in particular on the behavior of very dense systems. For these configurations, the simulations show that polydispersity of the sheared particles is a crucial factor that determines the system response. Figures 1 and 2 below, that present the total force on each grain, show that even relatively small (10 %) nonuniformity of the size of the grains (expected in typical experiments) may lead to significant modifications of the system properties, such as velocity profiles, temperature, force propagation, and formation shear bands. The simulations are extended in a few other directions, in order to provide additional insight to the experimental system analyzed above. In one direction, both gravity, and driving due to vibrations are included. These simulations allow for predictions on the driving regime that is required in the experiments in order to analyze the jamming transition. Furthermore, direct comparison of experiments and DES will allow for verification of the modeling assumptions. We have also extended our modeling efforts to 3D. The (preliminary) results of these simulations of an annular system in zero-g environment will conclude the presentation.

Description: Diffusion flame stabilization is of essential importance in both Earth-bound combustion systems and spacecraft fire safety. Local extinction, re-ignition, and propagation processes may occur as a result of interactions between the flame zone and vortices or fire-extinguishing agents. By using a computational fluid dynamics code with a detailed chemistry model for methane combustion, the authors have revealed the chemical kinetic structure of the stabilizing region of both jet and flat-plate diffusion flames, predicted the flame stability limit, and proposed diffusion flame attachment and detachment mechanisms in normal and microgravity. Because of the unique geometry of the edge of diffusion flames, radical back-diffusion against the oxygen-rich entrainment dramatically enhanced chain reactions, thus forming a peak reactivity spot, i.e., reaction kernel, responsible for flame holding. The new results have been obtained for the edge diffusion flame propagation and attached flame structure using various C1-C3 hydrocarbons.

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: This paper addresses some theoretical modelling and control issues for a mixing chamber used in rocket engine testing at NASA Stennis Space Center. The mixer is responsible for combining high pressure LH2 and GH2 to produce a hydrogen flow that meets certain thermodynamic properties before it is fed into a test article. The desired properties are maintained by precise control of the LH2 and GH2 flows. The mixer is modelled as a general multi-flow lumped volume for single constituent fluids using density and internal energy as states. The set of nonlinear differential equations is modelled in the SIMULINK environment including a table look-up feature of the fluid thermodynamic properties. a small-signal (linear) model is developed based on the nonlinear model and simulated as well. Pulse disturbances are introduced to the valve positions and the quality of the linear model is ascertained by comparing its behavior against the nonlinear model simulations. Valve control strategies that simulate an operator-in-the-loop scenario are then explored demonstrating the need for automatic feedback control. Finally, classical optimal single-output and multi-output Proportional/Integral controllers are designed based on the linear model and applied to the nonlinear model with excellent results to track simultaneous, constant setpoint changes in desired exit flow, exit temperature, and mixer pressure, as well as to reject unmeasurable but bounded additive step perturbations in the valve positions.

Description: This paper addresses some modeling and control issues for a mixing chamber used in rocket engine testing at NASA Stennis Space Center. The mixer must combine high pressure liquid hydrogen (LH2) and gaseous hydrogen (GH2) to produce and output flow that meets certain thermodynamic properties before it is fed into a test article. More precisely, this paper considers that the quantities to be tracked and/or regulated are mixer internal pressure, exit mass flow, and exit temperature. The available control inputs are given by three value positions, namely those of the GH2, LH2 and exit valves. The mixer is modelled by a system of two nonlinear ordinary differential equations having density and internal energy as states. The model must be simple enough to lend itself to subsequent feedback controller design, yet its accuracy must be tested against real data. For this reason, the model includes function calls to thermodynamic property data. Some structural properties of the resulting model that pertain to controller design, such as controllability and uniqueness of the equilibrium point are shown to hold. Validation of the model against real data is also provided. As a first control approach, a small-signal (linear) model is developed based on the nonlinear model and simulated as well. Pulse disturbances are introduced to the valve positions and the quality of the linear model is ascertained by comparing its behavior against the nonlinear model simulations. Valve control strategies that simulate an operator-in-the-loop scenario are then explored demonstrating the need for automatic feedback control. Classical optimal single-output and multi-output Proportional/Integral controllers are designed based on the linear model and applied to the nonlinear model with excellent results to track simultaneous, constant setpoint changes in desired exit flow, exit temperature, and mixer pressure, as well as to reject unmeasurable but bounded additive step perturbations in the valve positions. A feedback linearization controller is designed and used to achieve tracking and regulation of the outputs over an extended range of the variables of interest.

Description: A mixing chamber used in rocket engine combustion testing at NASA Stennis Space Center is modeled by a second order nonlinear MIMO system. The mixer is used to condition the thermodynamic properties of cryogenic liquid propellant by controlled injection of the same substance in the gaseous phase. The three inputs of the mixer are the positions of the valves regulating the liquid and gas flows at the inlets, and the position of the exit valve regulating the flow of conditioned propellant. The outputs to be tracked and/or regulated are mixer internal pressure, exit mass flow, and exit temperature. The outputs must conform to test specifications dictated by the type of rocket engine or component being tested downstream of the mixer. Feedback linearization is used to achieve tracking and regulation of the outputs. It is shown that the system is minimum-phase provided certain conditions on the parameters are satisfied. The conditions are shown to have physical interpretation.

Description: Various configurations of characterization systems such as ion mobility spectrometers and mass spectrometers are disclosed that are coupled to an ionization device. The ionization device is formed of a membrane that houses electrodes therein that are located closer to one another than the mean free path of the gas being ionized. Small voltages across the electrodes generate large electric fields which act to ionize substantially all molecules passing therethrough without fracture. Methods to manufacture the mass spectrometer and ion mobility spectrometer systems are also described.

Description: Linear response (specifically, Fourier's Law) in He-4 has been observed to fail in heat flow experiments near the superfluid transition. A detailed analysis of the data suggests that the hydrostatic pressure gradient across the helium column limits the divergence of the correlation length in our earth-based experiments. This is consistent with other observations, such as the surprising lack of mutual friction and hysteresis near the superfluid transition, and a 'rounding' of the transition that appears to be independent of heat flux in the low heat flux limit. I will discuss these unusual results from earth-based measurements, and will show predictions for the very different results that may result when we make our measurements on orbit as part of the M1 Mission of the Low- Temperature, Microgravity Physics Facility. This work has been funded by the Fundamental Physics Discipline within the Physical Sciences Research Office of NASA, and is conducted by the DYNAMX (UNM) and CQ (Caltech) Groups, with assistance from the Low Temperature Science and Quantum Sensors Group at JPL.

Description: Heat pipes are among the most promising technologies for space radiator systems. The paper reports further evaluation of potential heat pipe fluids in the intermediate temperature range of 400 to 700 K in continuation of two recent reports. More thermo-physical property data are examined. Organic, inorganic and elemental substances are considered. The evaluation of surface tension and other fluid properties are examined. Halides are evaluated as potential heat pipe fluids. Reliable data are not available for all fluids and further database development in necessary. Many of the fluids considered are promising candidates as heat pipe fluids. Water is promising as a heat pipe fluid up to 500-550 K. Life test data for thermo-chemical compatibility are almost non-existent.

Description: Are eddies an important contributor to subduction in the eastern subtropical gyre? Here, an adjoint model is used to combine a regional, eddy-resolving numerical model with observations to produce a state estimate of the ocean circulation. The estimate is a synthesis of a variety of in- situ observations from the Subduction Experiment, TOPEX/POSEIDON altimetry, and the MTI General Circulation Model. The adjoint method is successful because the Northeast Atlantic Ocean is only weakly nonlinear. The state estimate provides a physically-interpretable, eddy-resolving information source to diagnose subduction. Estimates of eddy subduction for the eastern subtropical gyre of the North Atlantic are larger than previously calculated from parameterizations in coarse-resolution models. Furthermore, eddy subduction rates have typical magnitudes of 15% of the total subduction rate. Eddies contribute as much as 1 Sverdrup to water-mass transformation, and hence subduction, in the North Equatorial Current and the Azores Current. The findings of this thesis imply that the inability to resolve or accurately parameterize eddy subduction in climate models would lead to an accumulation of error in the structure of the main thermocline, even in the relatively-quiescent eastern subtropical gyre.

Description: The main goal of this work is to better understand foam behavior both on the Earth and in microgravity conditions and to determine the relation between a foam's structure and wetness and its rheological properties. Our experiments focused on the effects of the bubble size distribution (BSD) on the foam behavior under gradual or stepwise in the liquid flow rate and on the onset of the convective instability. We were able to show experimentally, that the BSD affects foam rheology very strongly so any theory must take foam texture into account.

Description: Building on over four decades of research and technology development related to the behavior of fluids in low gravity environments, the current NASA Microgravity Fluid Physics Program continues the quest for knowledge to further understand and design better fluids systems for use on earth and in space. NASA's Biological and Physical Research Enterprise seeks to exploit the space environment to conduct research supporting human exploration of space (strategic research), research of intrinsic scientific importance and impact (fundamental research), and commercial research. The strategic research thrust will build the vital knowledge base needed to enable NASA's mission to explore the Universe and search for life. There are currently five major research areas in the Microgravity Fluid Physics Program: complex fluids, niultiphase flows and phase change, interfacial phenomena, biofluid mechanics, and dynamics and instabilities. Numerous investigations into these areas are being conducted in both ground-based laboratories and facilities and in the flight experiments program. Most of the future NASA- sponsored flight experiments in microgravity fluid physics and transport phenomena will be carried out on the International Space Station (ISS) in the Fluids Integrated Rack (FIR), in the Microgravity Science Glovebox (MSG), in EXPRESS racks, and in other facilities provided by international partners. This paper presents an overview of the near- and long-term visions for NASA's Microgravity Fluid Physics Research Program and brief descriptions of hardware systems planned to enable this research.

Description: The Two-phase Flow, Fluid Stability and Dynamics Workshop was held on May 15, 2003 in Cleveland, Ohio to define a coherent scientific research plan and roadmap that addresses the multiphase fluid problems associated with NASA s technology development program. The workshop participants, from academia, industry and government, prioritized various multiphase issues and generated a research plan and roadmap to resolve them. This report presents a prioritization of the various multiphase flow and fluid stability phenomena related primarily to power, propulsion, fluid and thermal management and advanced life support; and a plan to address these issues in a logical and timely fashion using analysis, ground-based and space-flight experiments.

Description: System and methods are disclosed for fluid measurements which may be utilized to determine mass flow rates such as instantaneous mass flow of a fluid stream. In a preferred embodiment, the present invention may be utilized to compare an input mass flow to an output mass flow of a drilling fluid circulation stream. In one embodiment, a fluid flow rate is determined by utilizing a microwave detector in combination with an acoustic sensor. The acoustic signal is utilized to eliminate 2pi phase ambiguities in a reflected microwave signal. In another embodiment, a fluid flow rate may be determined by detecting a phase shift of an acoustic signal across two different predetermined transmission paths. A fluid density may be determined by detecting a calibrated phase shift of an acoustic signal through the fluid. In another embodiment, a second acoustic signal may be transmitted through the fluid to define a particular 2pi phase range which defines the phase shift. The present invention may comprise multiple transmitters/receivers operating at different frequencies to measure instantaneous fuel levels of cryogenic fuels within containers positioned in zero or near zero gravity environments. In one embodiment, a moveable flexible collar of transmitter/receivers may be utilized to determine inhomogenuities within solid rocket fuel tubes.

Description: Computational Fluid Dynamics (CFD) analyses of axisymmetric circular-arc boattail nozzles have been completed in support of NASA's Next Generation Launch Technology Program to investigate the effects of high-speed nozzle geometries on the nozzle internal flow and the surrounding boattail regions. These computations span the very difficult transonic flight regime, with shock-induced separations and strong adverse pressure gradients. External afterbody and internal nozzle pressure distributions computed with the Wind code are compared with experimental data. A range of turbulence models were examined in Wind, including an Explicit Algebraic Stress model (EASM). Computations on two nozzle geometries have been completed at freestream Mach numbers ranging from 0.6 to 0.9, driven by nozzle pressure ratios (NPR) ranging from 2.9 to 5. Results obtained on converging-only geometry indicate reasonable agreement to experimental data, with the EASM and Shear Stress Transport (SST) turbulence models providing the best agreement. Calculations completed on a converging-diverging geometry involving large-scale internal flow separation did not converge to a true steady-state solution when run with variable timestepping (steady-state). Calculations obtained using constant timestepping (time-accurate) indicate less variations in flow properties compared with steady-state solutions. This failure to converge to a steady-state solution was found to be the result of difficulties in using variable time-stepping with large-scale separations present in the flow. Nevertheless, time-averaged boattail surface pressure coefficient and internal nozzle pressures show fairly good agreement with experimental data. The SST turbulence model demonstrates the best over-all agreement with experimental data.

Description: In 1995, Carlos Jorquera left NASA s Jet Propulsion Laboratory (JPL) to focus on erasing the growing void between high-performance cameras and the requisite software to capture and process the resulting digital images. Since his departure from NASA, Jorquera s efforts have not only satisfied the private industry's cravings for faster, more flexible, and more favorable software applications, but have blossomed into a successful entrepreneurship that is making its mark with improvements in fields such as medicine, weather forecasting, and X-ray inspection. Formerly a JPL engineer who constructed imaging systems for spacecraft and ground-based astronomy projects, Jorquera is the founder and president of the three-person firm, Boulder Imaging Inc., based in Louisville, Colorado. Joining Jorquera to round out the Boulder Imaging staff are Chief Operations Engineer Susan Downey, who also gained experience at JPL working on space-bound projects including Galileo and the Hubble Space Telescope, and Vice President of Engineering and Machine Vision Specialist Jie Zhu Kulbida, who has extensive industrial and research and development experience within the private sector.

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

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

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

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

Description: This viewgraph presentation provides information on activities pertaining to the experimental characterization of gas/gas injector flowfields. An experimental testbed for uni-element gas/gas injector studies at realistic conditions has been fabricated and verified. Experiments for characterizing mixing/combustion of gas/gas injectors with raman spectroscopy have been initiated.

Description: This viewgraph presentation gives an overview on recent improvements in the Finite Difference Navier Stokes (FDNS) computational fluid dynamics (CFD) code and its associated process. The development of a utility, PreViewer, has essentially eliminated the creeping of simple human error into the FDNS Solution process. Extension of PreViewer to encapsulate the Domain Decompression process has made practical the routine use of parallel processing. The combination of CVS source control and ATS consistency validation significantly increases the efficiency of the CFD process.

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

Description: The specific heat at constant volume C(sob V) of a simple fluid diverges near its liquid-vapor critical point. However, gravity-induced density stratification due to the divergence of isothermal susceptibility hinders the direct comparison of the experimental data with the predictions of renormalization group theory. In the past, a microgravity environment has been considered essential to eliminate the density stratification. We propose to perform specific heat measurements of He-3 on the ground using a method to cancel the density stratification. A He-3 fluid layer will be heated from below, using the thermal expansion of the fluid to cancel the hydrostatic compression. A 6% density stratification at a reduced temperature of 10(exp -5) can be cancelled to better than 0.1% with a steady 1.7 micro K temperature difference across a 0.05 cm thick fluid layer. A conventional AC calorimetry technique will be used to determine the heat capacity. The minimized bulk density stratification with a relaxation time 6500 sec at a reduced temperature of 10(exp -5) will stay unchanged during 1 Hz AC heating. The smear of the specific heat divergence due to the temperature difference across the cell is about 0.1% at a reduced temperature of 10(exp -6). The combination of using High Resolution Thermometry with a 0.5 n K temperature resolution in the AC technique and the cancellation of the density stratification will enable C(sub V) to be measured down to a reduced temperature of 10(exp -6) with less than a 1% systematic error.

Description: When a large density stratification is no longer a problem in a microgravity environment, one would like to increase the sample size in order to increase the signal-to-noise ratio for a specific heat measurement. To reduce the equilibration time associated with the large sample size, we designed a cylindrical cell containing a stack of plates that separate the bulk fluid into 60 equally thin layers. To understand the thermal behavior of the whole cell, we analyzed the thermal behavior of a 2-D composite system of a cylindrical near-critical fluid layer in contact with a cylindrical copper plate. In this 2-D analysis, the circumference boundary of the two cylindrical layers is subjected to a step temperature change. The solution of this 2-D composite system includes the piston effect that speeds up the equilibration in the near-critical fluid layer and the pure diffusion in the copper plate. The results of this analysis indicate that the characteristic length for the equilibration of the stacked cell is determined by an effective thickness of a single fluid layer instead of the total height of the cylindrical cell.

Description: We have observed a new temperature-entropy wave that propagates opposite to the direction of a steady heat flux Q when the helium column is heated from above. This propagating mode is due to non-linear thermo-conductance of the helium sample in the self-organized state. Such a mode had been predicted to exist on the self-organized heat transport state for Q less than about 100 nW/sq cm. We confirm that this mode exists in this regime. However, we also observe that it propagates even when the helium is pushed away from the self-organized heat transport state into the normal state.

Description: The shape of the liquid-gas coexistence curve of He-3 very near the critical point (-2x10(exp -6) 〈 t 〈 -5x10(exp -3) was measured using the quasi-static thermogram method. The study was performed in Earth s gravitational field using two different height calorimetry cells, both originally designed for simultaneous measurements of the isochoric heat capacity, isothermal compressibility, and PVT. The heights of two cells were 0.5 mm and 4.8 cm. The uncertainty in measuring the phase transition temperature was typically +/-2 micro-K. The measured coexistence curve near the critical point was strongly affected by the gravitational field. Away from the critical point, the coexistence curve obtained using this technique was also consistent with the earlier work using the local density measurements of Pittman et al. The recent crossover parametric model of the equation-of-state are used to analyze the height-dependent measured coexistence curves. Data analyses have indicated that microgravity will permit measurements within two additional decades in reduced temperatures beyond the best gravity-free data obtained in Earth-bound experiments.

Description: In addition to its primary role of studying non-linear heat transport effects near the lambda transition of He-4, the DYNAMX apparatus is suitable for measurements of the specific heat and the velocity of second sound. We plan to take advantage of available time on orbit to make measurements in these areas near to the lambda transition. The specific heat work would be similar to LPE, aimed at improving our knowledge of the singularity in the bulk heat capacity at the transition, but would provide more accurate results close to the transition. It would focus roughly equally on each side of the transition and would be synergistic with the CQ experiment, providing wider-range data at Q = 0. The second sound measurements are made possible by the fast time constant and high resolution of the DYNAMX thermometers, which allow accurate time-of-flight measurements of second sound pulses. It appears possible to measure the second sound velocity to about 1% at a reduced temperature of t = 5x10(exp -8) by averaging over a moderate number of pulses. The data would complement and extend earlier ground-based measurements, leading to improved tests of the theory of static critical phenomena at the lambda transition.

Description: We present progress on evaporative cooling of Rb-87 atoms in our Holographic Atom Trap (HAT). The HAT is formed by the interference of five intersecting YAG laser beams: atoms are loaded from a vapor-cell MOT into the bright fringes of the interference pattern through the dipole force. The interference pattern is composed of Talbot fringes along the direction of propagation of the YAG beams, prior to evaporative cooling each Talbot fringe contains 300,000 atoms at 50 micro-K and peak densities of 2 x 10(exp 14)/cu cm. Evaporative cooling is achieved through adiabatically decreasing the intensity of the YAG laser. We present data and calculations covering a range of HAT geometries and cooling procedures.

Description: Break-throughs in the study of superfluid He-3 weak links and recent demonstration of Josephson effect in He-4 are a result of significant advances in ultra-sensitive transducer and nanofabrication technology. However, further progress in the performance of superfluid weak links and quantum rotation interferometry devices depends, in part, on reducing the mechanical noise and increasing the effective duty cycle of such devices. In existing devices, the DC Josephson effect is driven by chemical potential difference produced by a pressure applied across the weak link. We propose a novel drive technique, where the chemical potential is due to a controlled temperature difference. This technique promises to eliminate mechanical shock associated with the switch of the direction of applied pressure and to achieve 100% duty cycle. The thermally driven Josephson effect may also answer outstanding questions about dissipation in superfluid weak links.

Description: Over the years, many ground-based studies have been performed near liquid-gas critical points to elucidate the expected divergences in thermodynamic quantities. The unambiguous interpretation of these studies very near the critical point is hindered by a gravity-induced density stratification. However, these ground-based measurements can give insight into the crossover behavior between the asymptotic critical region near the transition and the mean field region farther away. We have completed a detailed analysis of heat capacity, susceptibility and coexistence curve measurements near the He-3 liquid-gas critical point using the minimal-subtraction renormalization (MSR) scheme within the phi(exp 4) model. This MSR scheme, using only two adjustable parameters, provides a reasonable global fit to all of these experimental measurements in the gravity-free region out to a reduced temperature of |t| approx. 2x10(exp -2). Recently this approach has also been applied to the earlier microgravity measurements of Haupt and Straub in SF(sub 6) with surprising results. The conclusions drawn from the MSR analyses will be presented. Measurements in the gravity-affected region closer to the He-3 critical point have also been analyzed using the recent crossover parametric model (CPM) of the equation-of-state. The results of fitting heat capacity measurements to the CPM model along the He-3 critical isochore in the gravity-affected region will also be presented.

Description: We present measurements of the thermal conductivity near the superfluid transition of He-4 in confined geometries. The confinements we have studied include: cylindrical geometries with radii L=.5 and 1.0 microns, and parallel plates with 5 micron spacing. For L=1.0 microns, measurements at six pressures were conducted, whereas only SVP measurements have been done for other geometries. For the 1-D confinement in cylinders, the data are consistent with a universal scaling for all pressures at and above T(sub lambda). There are indications of breakdown of scaling and universality below T(sub lambda). For the 2-D confinement between parallel plates, the preliminary results indicate that the thermal conductivity is finite at the bulk superfluid transition temperature. Further analyses are needed to compare the 2-D results with those in bulk and 1-D confinement.

Description: We present a formal thermodynamic treatment of superfluid flow in a Josephson junction. We show that the current and the phase difference are thermodynamic conjugate variables. We derive quantitative expressions for the rms fluctuations of these variables. Also, we discuss the thermodynamic stability and the thermal activation to the phase slip region. We apply the developed formalism to show why an array of apertures in He-4 can exhibit the Josephson effect near the Lambda transition despite strong thermal fluctuations.

Description: A vortex-loop theory of the superfluid lambda transition has been developed over the last decade, with many results in agreement with experiments. It is a very simple theory, consisting of just three basic equations. When it was first proposed the main uncertainty in the theory was the use Flory scaling to find the fractal dimension of the random-walking vortex loops. Recent developments in high-resolution Monte Carlo simulations have now made it possible to verify the accuracy of this Flory-scaling assumption. Although the loop theory is not yet rigorously proven to be exact, the Monte Carlo results show at the least that it is an extremely good approximation. Recent loop calculations of the critical Casimir effect in helium films in the superfluid phase T 〈 Tc will be compared with similar perturbative RG calculations in the normal phase T 〉 Tc; the two calculations are found to match very nicely right at Tc.

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

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

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

Description: The present research concerns the development of high-frequency pressure and temperature probes and related instrumentation capable of performing spectral characterization of unsteady pressure and temperature fluctuations over the 0.05 20 kHz range, at the exit of a gas turbine combustor operating at conditions close to nominal ones for large power generation turbomachinery. The probes used a transient technique pioneered at Oxford University; in order to withstand exposure to the harsh environment the probes were fitted on a rapid injection and cooling system jointly developed by Centrospazio CPR and Syracuse University. The experimental runs were performed on a large industrial test rig being operated by ENEL Produzione. The achieved results clearly show the satisfactory performance provided by this diagnostic tool, even though the poor location of the injection port prevented the tests from yielding more insight of the core flow turbulence characteristics. The pressure and temperature probes survived several dozen injections in the combustor hot jet, while consistently providing the intended high frequency performance. The apparatus was kept connected to the combustor during long duration firings, operating as an unobtrusive, self contained, piggy-back experiment: high frequency flow samplings were remotely recorded at selected moments corresponding to different combustor operating conditions.

Description: A system and a method is provided for measuring three dimensional velocities at a plurality of points in a fluid employing at least two cameras positioned approximately perpendicular to one another. Image frames captured by the cameras may be filtered using background subtraction with outlier rejection with spike-removal filtering. The cameras may calibrated to accurately represent image coordinates in a world coordinate system using calibration grids modified using warp transformations. The two-dimensional views of the cameras may be recorded fur image processing and particle track determination. The tracer particles may be tracked on a two-dimensional basis and then stereo matched to obtain three-dimensional locations of the particles as a function of time so that velocities can be measured there from.

Description: Heat conduction plays an important role in the efficiency and life span of electronic components. To keep electronic components running efficiently and at a proper temperature, thermal management systems transfer heat generated from the components to thermal surfaces such as heat sinks, heat pipes, radiators, or heat spreaders. Thermal surfaces absorb the heat from the electrical components and dissipate it into the environment, preventing overheating. To ensure the best contact between electrical components and thermal surfaces, thermal interface materials are applied. In addition to having high conductivity, ideal thermal interface materials should be compliant to conform to the components, increasing the surface contact. While many different types of interface materials exist for varying purposes, Energy Science Laboratories, Inc. (ESLI), of San Diego, California, proposed using carbon velvets as thermal interface materials for general aerospace and electronics applications. NASA s Johnson Space Center granted ESLI a Small Business Innovation Research (SBIR) contract to develop thermal interface materials that are lightweight and compliant, and demonstrate high thermal conductance even for nonflat surfaces. Through Phase II SBIR work, ESLI created Vel-Therm for the commercial market. Vel-Therm is a soft, carbon fiber velvet consisting of numerous high thermal conductivity carbon fibers anchored in a thin layer of adhesive. The velvets are fabricated by precision cutting continuous carbon fiber tows and electrostatically flocking the fibers into uncured adhesive, using proprietary techniques.

Description: A NASA-developed aerodynamic simulation tool is ensuring the safety of future space operations while providing designers and engineers with an automated, highly accurate computer simulation suite. Cart3D, co-winner of NASA's 2002 Software of the Year award, is the result of over 10 years of research and software development conducted by Michael Aftosmis and Dr. John Melton of Ames Research Center and Professor Marsha Berger of the Courant Institute at New York University. Cart3D offers a revolutionary approach to computational fluid dynamics (CFD), the computer simulation of how fluids and gases flow around an object of a particular design. By fusing technological advancements in diverse fields such as mineralogy, computer graphics, computational geometry, and fluid dynamics, the software provides a new industrial geometry processing and fluid analysis capability with unsurpassed automation and efficiency.

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

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

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

Description: SLI engine designs will require pumps to throttle over a wide flow range while maintaining high performance. Unsteadiness generated by impeller/diffuser interaction is one of the major factors affecting off-design performance. Initial unsteady simulations are completed for impeller/diffuser stage. The Corsair simulations will continue across a wide flow range and for inducer/impeller/diffuser combinations. Results of unsteady simulations are being used to guide and explore new designs.

Description: A baseline solution for CFD Point 1 (Mach 24) in the STS-107 accident investigation was modified to include effects of: (1) holes through the leading edge into a vented cavity; and (2) a scarfed, conical nozzle directed toward the centerline of the vehicle from the forward, inboard corner of the landing gear door. The simulations were generated relatively quickly and early in the investigation because simplifications were made to the leading edge cavity geometry and an existing utility to merge scarfed nozzle grid domains with structured baseline external domains was implemented. These simplifications in the breach simulations enabled: (1) a very quick grid generation procedure; and (2) high fidelity corroboration of jet physics with internal surface impingements ensuing from a breach through the leading edge, fully coupled to the external shock layer flow at flight conditions. These simulations provided early evidence that the flow through a two-inch diameter (or larger) breach enters the cavity with significant retention of external flow directionality. A normal jet directed into the cavity was not an appropriate model for these conditions at CFD Point 1 (Mach 24). The breach diameters were of the same order or larger than the local, external boundary-layer thickness. High impingement heating and pressures on the downstream lip of the breach were computed. It is likely that hole shape would evolve as a slot cut in the direction of the external streamlines. In the case of the six-inch diameter breach the boundary layer is fully ingested. The intent of externally directed jet simulations in the second scenario was to approximately model aerodynamic effects of a relatively large internal wing pressure, fueled by combusting aluminum, which deforms the corner of the landing gear door and directs a jet across the windside surface. These jet interactions, in and of themselves, were not sufficiently large to explain observed aerodynamic behavior.

Description: The propagation of one-dimensional detonations in ethylene-air mixtures is investigated numerically by solving the one-dimensional Euler equations with detailed finite-rate chemistry. The numerical method is based on a second-order spatially accurate total-variation-diminishing scheme and a point implicit, first-order-accurate, time marching algorithm. The ethylene-air combustion is modeled with a 20-species, 36-step reaction mechanism. A multi-level, dynamically adaptive grid is utilized, in order to resolve the structure of the detonation. Parametric studies over an equivalence ratio range of 0.5 less than phi less than 3 for different initial pressures and degrees of detonation overdrive demonstrate that the detonation is unstable for low degrees of overdrive, but the dynamics of wave propagation varies with fuel-air equivalence ratio. For equivalence ratios less than approximately 1.2 the detonation exhibits a short-period oscillatory mode, characterized by high-frequency, low-amplitude waves. Richer mixtures (phi greater than 1.2) exhibit a low-frequency mode that includes large fluctuations in the detonation wave speed; that is, a galloping propagation mode is established. At high degrees of overdrive, stable detonation wave propagation is obtained. A modified McVey-Toong short-period wave-interaction theory is in excellent agreement with the numerical simulations.

Description: Widely used for the modeling of gas flows through the computation of the motion and collisions of representative molecules, the Direct Simulation Monte Carlo method has become the gold standard for producing research and engineering predictions in the field of rarefied gas dynamics. Direct Simulation Monte Carlo was first introduced in the early 1960s by Dr. Graeme Bird, a professor at the University of Sydney, Australia. It has since proved to be a valuable tool to the aerospace and defense industries in providing design and operational support data, as well as flight data analysis. In 2002, NASA brought to the forefront a software product that maintains the same basic physics formulation of Dr. Bird's method, but provides effective modeling of complex, three-dimensional, real vehicle simulations and parallel processing capabilities to handle additional computational requirements, especially in areas where computational fluid dynamics (CFD) is not applicable. NASA's Direct Simulation Monte Carlo Analysis Code (DAC) software package is now considered the Agency s premier high-fidelity simulation tool for predicting vehicle aerodynamics and aerothermodynamic environments in rarified, or low-density, gas flows.

Description: This viewgraph presentation discusses unsteady flow simulations for a turbopump intended for a reusable launch vehicle (RLV). The simulation process makes use of computational grids and parallel processing. The architecture of the parallel computers used is discussed, as is the scripting of turbopump simulations.

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

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

Description: This paper highlights results from the first of a series of hierarchical simulations aimed at assessing the modeling requirements for application of the large eddy simulation technique to cryogenic injection and combustion processes in liquid rocket engines. The focus is on liquid-oxygen-hydrogen coaxial injectors at a condition where the liquid-oxygen is injected at a subcritical temperature into a supercritical environment. For this situation a diffusion dominated mode of combustion occurs in the presence of exceedingly large thermophysical property gradients. Though continuous, these gradients approach the behavior of a contact discontinuity. Significant real gas effects and transport anomalies coexist locally in colder regions of the flow, with ideal gas and transport characteristics occurring within the flame zone. The current focal point is on the interfacial region between the liquid-oxygen core and the coaxial hydrogen jet where the flame anchors itself.

Description: Marshall Space Flight Center (MSFC) has demonstrated two measurement devices for sensing and resolving the hydrodynamic loads on fluid machinery. The first - a derivative of the six-component wind tunnel balance - senses the forces and moments on the rotating device through a weakened shaft section instrumented with a series of strain gauges. This rotating balance was designed to directly measure the steady and unsteady hydrodynamic loads on an inducer, thereby defining both the amplitude and frequency content associated with operating in various cavitation modes. The second device - a high frequency response pressure transducer surface mounted on a rotating component - was merely an extension of existing technology for application in water. MSFC has recently completed experimental evaluations of both the rotating balance and surface-mount transducers in a water test loop. The measurement bandwidth of the rotating balance was severely limited by the relative flexibility of the device itself, resulting in an unexpectedly low structural bending mode and invalidating the higher-frequency response data. Despite these limitations, measurements confirmed that the integrated loads on the four-bladed inducer respond to both cavitation intensity and cavitation phenomena. Likewise, the surface-mount pressure transducers were subjected to a range of temperatures and flow conditions in a non-rotating environment to record bias shifts and transfer functions between the transducers and a reference device. The pressure transducer static performance was within manufacturer's specifications and dynamic response accurately followed that of the reference.

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

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

Description: Computational Fluid Dynamics (CFD) flow solutions for the M2129 diffusing S-duct with and without vane effectors were computed by the Wind-US flow solver. Both structured and unstructured 3-D grids were used. Without vane effectors, the duct exhibited massive flow separation in both experiment and CFD. With vane effectors installed, the flow remained attached and aerodynamic losses were reduced. Total pressure recovery and distortion near the duct outlet were computed from the solutions and compared favorably to experimental values. These calculations are part of a validation effort for the Wind-US code. They also provide an example case to aid engineers in learning to use the Wind-US software.

Description: This work deals with the documentation and control of flow separation that occurs over turbine blades in the low-pressure turbine stage at low Reynolds numbers that exist at high altitude cruise. We utilize a specially constructed linear cascade that is designed to study the flow field over a generic LPT cascade consisting of Pratt & Whitney 'Pak B' shaped blades. This facility was constructed under a previous one-year NASA Glenn RC initiative. The center blade in the cascade is instrumented to measure the surface pressure coefficient distribution. Optical access allows two-component LDV measurement for boundary layer profiles. Experimental conditions have been chosen to give a range of chord Reynolds numbers from 10 to 100K, and a range of free-stream turbulence levels from u'/U(sub infinity)= 0.08 to 3 percent. The surface pressure measurements were used to define a region of separation and reattachment that depend on the free-stream conditions. The location of separation was found to be relatively insensitive to the experimental conditions. However, reattachment location was very sensitive to the turbulence level and Reynolds number. Excellent agreement was found between the measured pressure distributions and predictions from Euler and RANS simulations. Two-component LDV measurements are presently underway to document the mean and fluctuating velocity components in the boundary layer over the center blade for the range of experimental conditions. The fabrication of the plasma actuator is underway. These are designed to produce either streamwise vortices, or a downstream-directed wall jet. A precursor experiment for the former approach was performed with an array of vortex generators placed just upstream of the separation line. These led to reattachment except for the lowest Reynolds number. Progress has also been made on the proposed concept for a laterally moving wake. This involved constructing a smaller wind tunnel and molding an array of symmetric airfoils to form an array. Following its development, it will be scaled up and used to introduce lateral moving wakes upstream up the Pak-B cascade.

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: A computational and experimental study was conducted to investigate the effects of multiple injection ports in a two-dimensional, convergent-divergent nozzle, for fluidic thrust vectoring. The concept of multiple injection ports was conceived to enhance the thrust vectoring capability of a convergent-divergent nozzle over that of a single injection port without increasing the secondary mass flow rate requirements. The experimental study was conducted at static conditions in the Jet Exit Test Facility of the 16-Foot Transonic Tunnel Complex at NASA Langley Research Center. Internal nozzle performance was obtained at nozzle pressure ratios up to 10 with secondary nozzle pressure ratios up to 1 for five configurations. The computational study was conducted using the Reynolds Averaged Navier-Stokes computational fluid dynamics code PAB3D with two-equation turbulence closure and linear Reynolds stress modeling. Internal nozzle performance was predicted for nozzle pressure ratios up to 10 with a secondary nozzle pressure ratio of 0.7 for two configurations. Results from the experimental study indicate a benefit to multiple injection ports in a convergent-divergent nozzle. In general, increasing the number of injection ports from one to two increased the pitch thrust vectoring capability without any thrust performance penalties at nozzle pressure ratios less than 4 with high secondary pressure ratios. Results from the computational study are in excellent agreement with experimental results and validates PAB3D as a tool for predicting internal nozzle performance of a two dimensional, convergent-divergent nozzle with multiple injection ports.

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

Description: Progress can be reported in two areas related to characterizing the properties of cusp diamagnetic cavities. Laboratory terrella experiments have been conducted for the purpose of using neutral gas excitation as a tracer of trapped electron populations in the presence of two dipoles that are used to develop a magnetic cusp topology. Figure 1 and 2 show top and side views of two configurations. Dipole trapped electron populations appear as the two luminous annular rings. Other populations are the most intense regions are shown. Interspersed between these regions are narrow regions that represent the topological cusps in these configurations. That they contain luminous gas is evidence for cusp trapping similar to what we believe exists in the terrestrial magnetosphere. The asymmetry of these cusp regions as seen in Figure 1 is the result of a relative tilt between the two dipoles suggestive of what would be expected in space. It is in these regions that particle observations were sought, so as to validate the realization of proposed and laboratory achieved trapping in a diamagnetic cusp. Figure 3 shows particle trajectories in a modeled cusp magnetic topology for three particle energies. Blue, green, and red traces correspond to increasing energies. Due to factors discussed outside of this final report, a thorough exploration of relevant satellite observations have not been achieved.

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: Mixing driven by buoyancy-induced flows inside a cavity consists of stretching and folding of an interface. Measurement of the flow field using particle imaging velocimetry shows that during stretching the flow field has a single elliptic point, thus dominated by a single vortex. However, global bifurcation that results in folding introduces a hyperbolic point whereby the flow field degenerates to multiple vortex interactions. The short-lived coherent structure observed during mixing which results in the Rayleigh- Taylor morphology is attributed to vortex interactions. The mixing characteristics of non-homogeneous fluids driven by buoyancy are important towards understanding transport phenomenon in a microgravity environment. Mixing consists of stretching and folding of an interface due to a flow field whose intensity depends on the body force. For miscible liquids, the characteristic of the flow field determines whether mass transport is governed by diffusion or bulk stirring which induces mixing. For technologically important processes, transport of mass is governed by the coupling of the body force to scalar gradients such as concentration and or temperature' 2 3 . In order to lend insight into these classes of problems we consider a model experimental system to study mixing driven by buoyancy-induced flows. The characteristics of mixing is addressed from detail measurements of the flow field using particle imaging velocimetry (PIV), and its corresponding interface dynamics using image processing techniques.

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

Description: Turbulence is a universal, nonlinear phenomenon found in all energetic fluid and plasma motion. In particular. understanding magneto hydrodynamic (MHD) turbulence and incorporating its effects in the computation and prediction of the flow of ionized gases in space, for example, are great challenges that must be met if such computations and predictions are to be meaningful. Although a general solution to the "problem of turbulence" does not exist in closed form, numerical integrations allow us to explore the phase space of solutions for both ideal and dissipative flows. For homogeneous, incompressible turbulence, Fourier methods are appropriate, and phase space is defined by the Fourier coefficients of the physical fields. In the case of ideal MHD flows, a fairly robust statistical mechanics has been developed, in which the symmetry and ergodic properties of phase space is understood. A discussion of these properties will illuminate our principal discovery: Coherent structure and randomness co-exist in ideal MHD turbulence. For dissipative flows, as opposed to ideal flows, progress beyond the dimensional analysis of Kolmogorov has been difficult. Here, some possible future directions that draw on the ideal results will also be discussed. Our conclusion will be that while ideal turbulence is now well understood, real turbulence still presents great challenges.

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: It is well known that numerical warm season quantitative precipitation forecasts lack significant skill for numerous reasons. Some are related to the model--it may lack physical processes required to realistically simulate convection or the numerical algorithms and dynamics employed may not be adequate. Others are related to initialization-mesoscale features play an important role in convective initialization and atmospheric observation systems are incapable of properly depicting the three-dimensional stability structure at the mesoscale. The purpose of this study is to determine if a mesoscale model initialized with a diabatic initialization scheme can improve short-term (0 to 12h) warm season quantitative precipitation forecasts in the Southeastern United States. The Local Analysis and Prediction System (LAPS) developed at the Forecast System Laboratory is used to diabatically initialize the Pennsylvania State University/National center for Atmospheric Research (PSUNCAR) Mesoscale Model version 5 (MM5). The SPORT Center runs LAPS operationally on an hourly cycle to produce analyses on a 15 km covering the eastern 2/3 of the United States. The 20 km National Centers for Environmental Prediction (NCEP) Rapid Update Cycle analyses are used for the background fields. Standard observational data are acquired from MADIS with GOES/CRAFT Nexrad data acquired from in-house feeds. The MM5 is configured on a 140 x 140 12 km grid centered on Huntsville Alabama. Preliminary results indicate that MM5 runs initialized with LAPS produce improved 6 and 12h QPF threat scores compared with those initialized with the NCEP RUC.

Description: Four Sunpower M87N Stirling-cycle cryocoolers will be used to extend the lifetime of the Alpha Magnetic Spectrometer-02 (AMS-02) experiment. The cryocoolers will be mounted to the AMS-02 vacuum case using a structure that will thermally and mechanically decouple the cryocooler from the vacuum case while providing compliance to allow force attenuation using a passive balancer system. The cryocooler drive is implemented using a 60Hz pulse duration modulated square wave. Details of the testing program, mounting assembly and drive scheme will be presented. AMS-02 is a state-of-the-art particle physics detector containing a large superfluid helium-cooled superconducting magnet. Highly sensitive detector plates inside the magnet measure a particle s speed, momentum, charge, and path. The AMS-02 experiment, which will be flown as an attached payload on the International Space Station, will study the properties and origin of cosmic particles and nuclei including antimatter and dark matter. Two engineering model cryocoolers have been under test at NASA Goddard since November 2001. Qualification testing of the engineering model cryocooler bracket assembly is near completion. Delivery of the flight cryocoolers to Goddard is scheduled for September 2003.

Description: Experiments have shown that moderate turbulence levels can nearly double turbine blade stagnation region heat transfer. Data have also shown that heat transfer is strongly affected by the scale of turbulence as well as its level. In addition to the stagnation region, turbulence is often seen to increase pressure surface heat transfer. This is especially evident at low to moderate Reynolds numbers. Vane and rotor stagnation region, and vane pressure surface heat transfer augmentation is often seen in a pre-transition environment. Accurate predictions of transition and relaminarization are critical to accurately predicting blade surface heat transfer. An approach is described which incorporates the effects of both turbulence level and scale into a CFD analysis. The model is derived from experimental data for cylindrical and elliptical leadng edges. Results using this model are compared to experimental data for both vane and rotor geometries. The comparisons are made to illustrate that using a model which includes the effects of turbulence length scale improves agreement with data, and to illustrate where improvements in the modeling are needed.

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