ALBERT

Description: A 3D CFD code, SCISEAL, has been developed and validated. Its capabilities include cylindrical seals, and it is employed on labyrinth seals, rim seals, and disc cavities. State-of-the-art numerical methods include colocated grids, high-order differencing, and turbulence models which account for wall roughness. SCISEAL computes efficient solutions for complicated flow geometries and seal-specific capabilities (rotor loads, torques, etc.).

Description: The uses of virtual environment technology in the space program are examined with emphasis on training for the Hubble Space Telescope Repair and Maintenance Mission in 1993. Project ScienceSpace at the Virtual Environment Technology Lab is discussed.

Description: A Workshop on "Nitrogen Dynamics in Controlled Systems" was held September 26-28, 1995 at the Lawrence Berkeley National Laboratory. The meetings were sponsored by the NASA Advanced Life Support program and the Lawrence Berkeley National Laboratory, and hosted by Prof. Lester Packer of the University of California at Berkeley, and of the Lawrence Berkeley National Laboratory. The Workshop participants were asked to: 1. summarize current knowledge on the cycling of nitrogen in closed systems; 2. identify the needs that closed systems may have for specific forms of nitrogen; 3. identify possible ways of generating and maintaining (or avoiding) specific forms and concentrations of nitrogen; 4. compare biological and physical/chemical methods of transforming nitrogen.

Description: As part of bioregenerative life support feasibility testing by NASA, crop residues are being used to resupply elemental nutrients to recirculating hydroponic crop production systems. Methods for recovering nutrients from crop residues have evolved from water soaking (leaching) to rapid aerobic bioreactor processing. Leaching residues recovered the majority of elements but it also recovered significant amounts of soluble organics. The high organic content of leachates was detrimental to plant growth. Aerobic bioreactor processing reduced the organic content ten-fold, which reduced or eliminated phytotoxic effects. Wheat and potato production studies were successful using effluents from reactors having with 8- to 1-day retention times. Aerobic bioreactor effluents supplied at least half of the crops elemental mass needs in these studies. Descriptions of leachate and effluent mineral content, biomass productivity, microbial activity, and nutrient budgets for potato and wheat are presented.

Description: Staphylococcus aureus was isolated over 2 years from Space Shuttle mission crewmembers to determine dissemination and retention of bacteria. Samples before and after each mission were from nasal, throat, urine, and feces and from air and surface sampling of the Space Shuttle. DNA fingerprinting of samples by digestion of DNA with SmaI restriction endonuclease followed by pulsed-field gel electrophoresis showed S. aureus from each crewmember had a unique fingerprint and usually only one strain was carried by an individual. There was only one instance of transfer between crewmembers. Strains from interior surfaces after flight matched those of crewmembers, suggesting microbial fingerprinting may have forensic application.

Description: The flow in a hydrostatic pocket is described by a mathematical model that uses the three dimensional Navier-Stokes equations written in terms of the primary variables, u, v, w, and p. Using a conservative formulation, a finite volume multi-block method is applied through a collocated, body fitted grid. The flow is simulated in a shallow pocket with a depth/length ratio of 0.02. The flow structures obtained and described by the authors in their previous two dimensional models are made visible in their three dimensional aspect for the Couette flow. It has been found that the flow regimes formed central and secondary vortical cells with three dimensional corkscrew-like structures that lead the fluid on an outward bound path in the axial direction of the pocket. The position of the central vortical cell center is at the exit region of the capillary restrictor feedline. It has also been determined that a fluid turn around zone occupies all the upstream space between the floor of the pocket and the runner, thus preventing any flow exit through the upstream port. The corresponding pressure distribution under the shaft presented as well. It was clearly established that for the Couette dominated case the pressure varies significantly in the pocket in the circumferential direction, while its variation is less pronounced axially.

Description: An experimental program was conducted at NASA Langley Research Center that included development and evaluation of an operational facility for wall drag measurement of potential scramjet fuel injection or wall cooling configurations. The facility consisted of a supersonic tunnel, with one wall composed of a series of interchangeable aluminum plates attached to an air bearing suspension system. The system was equipped with load cells that measured drag forces of 115 psia (793 kPa). This flow field contained a train of weak, unsteady, reflecting shock waves that were produced in the Mach 2 nozzle flows, the effect of reflecting shocks (which are to be expected in scramjet combustors) in internal flows has not previously been documented.

Description: The wall drag test tunnel at NASA Langley Research Center was used to evaluate simulated scramjet fuel injection into a wall cavity. In this tunnel, one wall consists of interchangeable aluminum plates attached to an air bearing suspension system. The plates were equipped with load cells to measure drag forces and static taps to determine pressure distributions. The plates were exposed to a Mach 2 air stream at a total pressure of 115 psia (793 kPa). This flow field contained a train of weak unsteady, reflecting shock waves that were produced in the nozzle assembly located upstream of the test section.

Description: Life support systems represent one of the most critical aspects of human space exploration. Future long-term missions such as the establishment of human-tended Lunar and Martian bases require closed life support systems. A conceptual approach to an Engineered Closed/Controlled EcoSystem incorporating bioregenerative capabilities by integrating humans, plants, and waste management processes is presented. The integration of physical/chemical and biological waste treatment processes is suitable for supporting plant growth through hydroponics and materially closing the human and plant metabolic loops. This conceptual design separates wastes into individual loops for treatment according to the specific metabolic needs of humans and plants. The means through which an integrated Engineered Closed/Controlled EcoSystem meets the life support objectives of long-term space habitation are summarized.

Description: No abstract available

Description: Here we report on the in situ performance of inexpensive, miniature sensors that have increased our ability to measure mass and energy fluxes from plant canopies in controlled environments: 1. Surface temperature. Canopy temperature measurements indicate changes in stomatal aperture and thus latent and sensible heat fluxes. Infrared transducers from two manufacturers (Exergen Corporation, Newton, MA; and Everest Interscience, Tucson, AZ, USA) have recently become available. Transducer accuracy matched that of a more expensive hand-held infrared thermometer. 2. Air velocity varies above and within plant canopies and is an important component in mass and energy transfer models. We tested commercially-available needle, heat-transfer anemometers (1 x 50 mm cylinder) that consist of a fine-wire thermocouple and a heater inside a hypodermic needle. The needle is heated and wind speed determined from the temperature rise above ambient. These sensors are particularly useful in measuring the low wind speeds found within plant canopies. 3. Accurate measurements of air temperature adjacent to plant leaves facilitates transport phenomena modeling. We quantified the effect of radiation and air velocity on temperature rise in thermocouples from 10 to 500 micrometers. At high radiation loads and low wind speeds, temperature errors were as large as 7 degrees C above air temperature.

Description: Atmospheres of enclosed environments in which 20 m2 stands of wheat, potato, and lettuce were grown were characterized and quantified by gas chromatography-mass spectrometry. A large number (in excess of 90) of volatile organic compounds (VOCs) were identified in the chambers. Twenty eight VOC's were assumed to be of biogenic origin for these were not found in the chamber atmosphere when air samples were analyzed in the absence of plants. Some of the compounds found were unique to a single crop. For example, only 35% of the biogenic compounds detected in the wheat atmosphere were unique to wheat, while 36% were unique to potato and 26% were unique to lettuce. The number of compounds detected in the wheat (20 compounds) atmosphere was greater than that of potato (11) and lettuce (15) and concentration levels of biogenic and non-biogenic VOC's were similar.

Description: The Biomass Production Chamber (BPC) located at Kennedy Space Center, FL, USA provides a large (20 m2 area, 113 m3 vol.), closed environment for crop growth tests for NASA's Controlled Ecological Life Support System (CELSS) program. Since the summer of 1988, the chamber has operated on a near-continuous basis (over 1200 days) without any major failures (excluding temporary power losses). During this time, five crops of wheat (64-86 days each), three crops of soybean (90 to 97 days), five crops of lettuce (28-30 days), and four crops of potato (90 to 105 days were grown, producing 481 kg of dry plant biomass, 196 kg edible biomass, 540 kg of oxygen, 94,700 kg of condensed water, and fixing 739 kg of carbon dioxide. Results indicate that total biomass yields were close to expected values for the given light input, but edible biomass yields and harvest indices were slightly lower than expected. Stand photosynthesis, respiration, transpiration, and nutrient uptake rates were monitored throughout growth and development of the different crops, along with the build-up of ethylene and other volatile organic compounds in the atmosphere. Data were also gathered on system hardware maintenance and repair, as well as person-hours required for chamber operation. Future tests will include long-term crop production studies, tests in which nutrients from waste treatment systems will be used to grow new crops, and multi-species tests.

Description: The Porous Tube Plant Nutrient Delivery System (PTPNDS), a hydrophilic, microporous ceramic tube hydroponic system designed for microgravity, will be tested in a middeck locker of the Space Shuttle. The flight experiment will focus on hardware operation and assess its ability to support seed germination and early seedling growth in microgravity. The water controlling system of the PTPNDS hardware has been successfully tested during the parabolic flight of the KC-135. One challenge to the development of the space flight experiment was to devise a method of holding seeds to the cylindrical porous tube. The seed-holder must provide water and air to the seed, absorb water from the porous tube, withstand sterilization, provide a clear path for shoots and roots to emerge, and be composed of flight qualified materials. In preparation for the flight experiment, a wheat seed-holder has been designed that utilizes a cellulose acetate plug to facilitate imbibition and to hold the wheat seeds in contact with the porous tube in the correct orientation during the vibration of launch and the microgravity environment of orbit. Germination and growth studies with wheat at a range of temperatures showed that optimal moisture was 78% (by weight) in the cellulose acetate seed holders. These and other design considerations are discussed.

Description: Potato (Solanum tuberosum L.) cv. 'Norland', vegetative growth and tuber productivity grown in the porous water and nutrient delivery system (PTNDS) developed by the Wisconsin Center for Space Automation and Robotics were compared with the vegetative growth and tuber productivity of plants grown in a peat:vermiculite potting mixture (PT/VR). The plants were grown at 12, 16, and 24-h light periods, 18 degrees C constant temperature, 70% relative humidity, and 300 micromol m-2 s-1 photosynthetic photon flux. Canopy height of plants grown in the PT/VR system was taller than that of plants grown in the PTNDS system. Canopy height differences were greatest when the plants were grown under a 24-h photoperiod. Leaf and stem dry masses were similar for plants grown in the two systems under the 12-h photoperiod. Under the 24-h photoperiod, leaf and stem dry masses of plants grown in the PT/VR system were more than 3 times those of plants grown in the PTNDS system. Tuber dry masses were similar for plants grown in the two systems under the 12-h photoperiod. Under the 24 h-photoperiod, tuber dry weights of plants grown in the PT/VR system were more than twice those of plants grown in the PTNDS system. A slightly higher harvest index (ratio of tuber weight to leaf plus stem weight) was noted for the plants grown in the PTNDS than for the plants grown in the PT/VR system. Plants grown in the PTNDS system at the 24-h photoperiod matured earlier than plants grown at this photoperiod in the PT/VR system. Vegetative growth and tuber productivity of plants grown under the 16-h photoperiod generally were intermediate to those noted for plants grown under the 12 and 24-h photoperiods. These results indicate that potato plants grown in a PTNDS system may require less plant growing volume, mature in a shorter time, and likely produce more tubers per unit area compared with plants grown in the PT/VR system. These plant characteristics are a distinct advantage for a plant growing unit of a CELSS.

Description: Assuming that crops grown in controlled ecological life-support systems (CELSS) should provide a basis for meals that are both nutritious and attractive (to taste and vision), and that CELSS diets on the moon or Mars or in space-craft during long voyages will have to be mostly vegetarian, a workshop was convened at the Johnson Space Center, Houston, Texas, U.S.A. on 19 to 21 January, 1994. Participants consisted of trained nutritionists and others; many of the approximately 18 presenters who discussed possible diets were practicing vegetarians, some for more than two decades. Considering all the presentations, seven conclusions (or points for discussion) could be formulated: nutritious vegetarian diets are relatively easily to formulate, vegetarian diets are healthy, variety is essential in vegetarian diets, some experiences (e.g., Bios-3 and Biosphere 2) are relevant to planning of CELSS diets, physical constraints will limit the choice of crops, a preliminary list of recommended crops can be formulated, and this line of research has some potential practical spinoffs. The list of crops and the reasons for including specific crops might be of interest to professionals in the field of health and nutrition as well as to those who are designing closed ecological systems.

Description: Purdue University, as well as the Johnson and Kennedy Space Centers and NASA Ames Research Center, are investigating approximately 5-10 plants that will be grown hydroponically to provide not only the energy and nutrients, but also the oxygen for humans habitating in Mars and lunar bases. The growth and nutritional status of rats fed either a control diet (adequate in all macro- and micronutrients) or a strict vegetarian diet consisting of 5 (vegan-5) or 10 (vegan-10) candidate crop species were investigated. In addition, vegan-10 diets were supplemented with mineral and/or vitamin mix at a level similar to the control diets to assess the effect of supplementation on nutrient status. The assessment of inedible plant material as an alternative food source was also investigated. Results of this study demonstrated that consumption of the vegan-10 diet significantly improved weight gain of rats compared to that for rats fed the vegan-5 diet. Mineral supplementation, at a level present in the control diet, to the vegan-10 diet improved growth and nutrient status, but growth was significantly lower compared to the control-fed rats. Inclusion of inedible plant material, high in ash content, improved some indices of nutrient status, without improving growth.

Description: Edible biomass from four crops of wheat (Triticum aestivum L.), four crops of lettuce (Lactuca sativa L.), four crops of potato (Solanum tuberosum L.), and three crops of soybean (Glycine max (L.) Merr.) grown in NASA's CELSS Biomass Production Chamber were analyzed for proximate composition. All plants were grown using recirculating nutrient (hydroponic) film culture with pH and electrical conductivity automatically controlled. Temperature and humidity were controlled to near optimal levels for each species and atmospheric carbon dioxide partial pressures were maintained near 100 Pa during the light cycles. Soybean seed contained the highest percentage of protein and fat, potato tubers and wheat seed contained the highest levels of carbohydrate, and lettuce leaves contained the highest level of ash. Analyses showed values close to data published for field-grown plants with several exceptions: In comparison with field-grown plants, wheat seed had higher protein levels; soybean seed had higher ash and crude fiber levels; and potato tubers and lettuce leaves had higher protein and ash levels. The higher ash and protein levels may have been a result of the continuous supply of nutrients (e.g., potassium and nitrogen) to the plants by the recirculating hydroponic culture.

Description: A perspective of the process of obtaining performance robustness of space exploration life support systems is presented. Some useful definitions are made, some relevant issues are addressed, and a procedure for establishing performance robustness, so far as it now known, is explained. An example is given to illustrate the procedure.

Description: This paper examines mission simulation as an approach to develop requirements for automation and robotics for Advanced Life Support Systems (ALSS). The focus is on requirements and applications for command and control, control and monitoring, situation assessment and response, diagnosis and recovery, adaptive planning and scheduling, and other automation applications in addition to mechanized equipment and robotics applications to reduce the excessive human labor requirements to operate and maintain an ALSS. Based on principles of systems engineering, an approach is proposed to assess requirements for automation and robotics using mission simulation tools. First, the story of a simulated mission is defined in terms of processes with attendant types of resources needed, including options for use of automation and robotic systems. Next, systems dynamics models are used in simulation to reveal the implications for selected resource allocation schemes in terms of resources required to complete operational tasks. The simulations not only help establish ALSS design criteria, but also may offer guidance to ALSS research efforts by identifying gaps in knowledge about procedures and/or biophysical processes. Simulations of a planned one-year mission with 4 crewmembers in a Human Rated Test Facility are presented as an approach to evaluation of mission feasibility and definition of automation and robotics requirements.

Description: The Regenerative Life Support Systems (RLSS) Test Bed at NASA's Johnson Space Center is an atmospherically closed, controlled environment facility for human testing of regenerative life support systems using higher plants in conjunction with physicochemical life support systems. The facility supports NASA's Advanced Life Support (ALS) Program. The facility is comprised of two large scale plant growth chambers, each with approximately 11 m2 growing area. The root zone in each chamber is configurable for hydroponic or solid media plant culture systems. One of the two chambers, the Variable Pressure Growth Chamber (VPGC), is capable of operating at lower atmospheric pressures to evaluate a range of environments that may be used in a planetary surface habitat; the other chamber, the Ambient Pressure Growth Chamber (APGC) operates at ambient atmospheric pressure. The air lock of the VPGC is currently being outfitted for short duration (1 to 15 day) human habitation at ambient pressures. Testing with and without human subjects will focus on 1) integration of biological and physicochemical air and water revitalization systems; 2) effect of atmospheric pressure on system performance; 3) planetary resource utilization for ALS systems, in which solid substrates (simulated planetary soils or manufactured soils) are used in selected crop growth studies; 4) environmental microbiology and toxicology; 5) monitoring and control strategies; and 6) plant growth systems design. Included are descriptions of the overall design of the test facility, including discussions of the atmospheric conditioning, thermal control, lighting, and nutrient delivery systems.

Description: Future NASA missions to explore the solar system will be long-duration missions, requiring human life support systems which must operate with very high reliability over long periods of time. Such systems must be highly regenerative, requiring minimum resupply, to enable the crews to be largely self-sufficient. These regenerative life support systems will use a combination of higher plants, microorganisms, and physicochemical processes to recycle air and water, produce food, and process wastes. A key step in the development of these systems is establishment of a human-rated test facility specifically tailored to evaluation of closed, regenerative life supports systems--one in which long-duration, large-scale testing involving human test crews can be performed. Construction of such a facility, the Advanced Life Support Program's (ALS) Human-Rated Test Facility (HRTF), has begun at NASA's Johnson Space Center, and definition of systems and development of initial outfitting concepts for the facility are underway. This paper will provide an overview of the HRTF project plan, an explanation of baseline configurations, and descriptive illustrations of facility outfitting concepts.

Description: This study compared the growth of potato plants on nutrients recycled from inedible potato biomass. Plants were grown for 105 days in recirculating, thin-film hydroponic systems containing four separate nutrient solution treatments: (1) modified half-strength Hoagland's (control), 2) liquid effluent from a bioreactor containing inedible potato biomass, 3) filtered (0.2 micrometer) effluent, and 4) the water soluble fraction of inedible potato biomass (leachate). Approximately 50% of the total nutrient requirement in treatments 2-4 were provided (recycled) from the potato biomass. Leachate had an inhibitory effect on leaf conductance, photosynthetic rate, and growth (50% reduction in plant height and 60% reduction in tuber yield). Plants grown on bioreactor effluent (filtered or unfiltered) were similar to the control plants. These results indicated that rapidly degraded, water soluble organic material contained in the inedible biomass, i.e., material in leachate, brought about phytotoxicity in the hydroponic culture of potato. Recalcitrant, water soluble organic material accumulated in all nutrient recycling treatments (650% increase after 105 days), but no increase in rhizosphere microbial numbers was observed.

Description: As space missions become longer in duration, the need to recycle waste into useful compounds rises dramatically. This problem can be addressed by the development of Controlled Ecological Life Support Systems (CELSS) (i.e., Engineered Closed/Controlled Eco-Systems (ECCES)), consisting of human and plant modules. One of the waste streams leaving the human module is urine. In addition to the reclamation of water from urine, recovery of the nitrogen is important because it is an essential nutrient for the plant module. A 3-step biological process for the recycling of nitrogenous waste (urea) is proposed. A packed-bed bioreactor system for this purpose was modeled, and the issues of reaction step segregation, reactor type and volume, support particle size, and pressure drop were addressed. Based on minimization of volume, a bioreactor system consisting of a plug flow immobilized urease reactor, a completely mixed flow immobilized cell reactor to convert ammonia to nitrite, and a plug flow immobilized cell reactor to produce nitrate from nitrite is recommended. It is apparent that this 3-step bioprocess meets the requirements for space applications.

Description: An anaerobic reactor seeded with organisms from an anaerobic lagoon was used to study the degradation of inedible crop residues from potato and wheat crops grown in a closed environment. Conversion of this biomass into other products was also evaluated. Degradation of wheat volatile solids was about 25% where that of potato was about 50%. The main product of the anaerobic fermentation of both crops was acetic acid with smaller quantities of propionate and butyrate produced. Nitrate, known to be high in concentration in inedible potato and wheat biomass grown hydroponically, was converted to ammonia in the anaerobic reactor. Both volatile fatty acid and ammonia production may have implications in a crop production system.

Description: No abstract available

Description: Growth of plants in a Controlled Ecological Life Support System (CELSS) may involve the use of hypobaric pressures enabling lower mass requirements for atmospheres and possible enhancement of crop productivity. A controlled environment plant growth chamber with hypobaric capability designed and built at Ames Research Center was used to determine if reduced pressures influence the rates of photosynthesis (Ps) and dark respiration (DR) of hydroponically grown lettuce plants. The chamber, referred to as a plant volatiles chamber (PVC), has a growing area of about 0.2 m2, a total gas volume of about 0.7 m3, and a leak rate at 50 kPa of 〈0.1%/day. When the pressure in the chamber was reduced from ambient to 51 kPa, the rate of net Ps increased by 25% and the rate of DR decreased by 40%. The rate of Ps increased linearly with decreasing pressure. There was a greater effect of reduced pressure at 41 Pa CO2 than at 81 Pa CO2. This is consistent with reports showing greater inhibition of photorespiration (Pr) in reduced O2 at low CO2 concentrations. When the partial pressure of O2 was held constant but the total pressure was varied between 51 and 101 kPa, the rate of CO2 uptake was nearly constant, suggesting that low pressure enhancement of Ps may be mainly attributable to lowered partial pressure of O2 and the accompanying reduction in Pr. The effects of lowered partial pressure of O2 on Ps and DR could result in substantial increases in the rates of biomass production, enabling rapid throughput of crops or allowing flexibility in the use of mass and energy resources for a CELSS.

Description: The effect of photoperiod (PP) on net carbon assimilation rate (Anet) and starch accumulation in newly mature canopy leaves of 'Norland' potato (Solanum tuberosum L.) was determined under high (412 varies as mol m-2s-1) and low (263 varies as mol m-2s-1) photosynthetic photon flux (PPF) conditions. The Anet decreased from 13.9 to 11.6 and 9.3 micromoles m-2s-1, and leaf starch increased from 70 to 129 and 118 mg g-1 drymass (DM) as photoperiod (PP) was increased from 12/12 to 18/6, and 24/0, respectively. Longer PP had a greater effect with high PPF conditions than with low PPF treatments, with high PPF showing greater decline in Anet. Photoperiod did not affect either the CO2 compensation point (50 micromoles mol-1) or CO2 saturation point (1100-1200 micromoles mol-1) for Anet. These results show an apparent limit to the amount of starch that can be stored (approximately 15% DM) in potato leaves. An apparent feedback mechanism exists for regulating Anet under high PPF, high CO2, and long PP, but there was no correlation between Anet and starch concentration in individual leaves. This suggests that maximum Anet cannot be sustained with elevated CO2 conditions under long PP (〉 or = 12 hours) and high PPF conditions. If a physiological limit exists for the fixation and transport of carbon,then increasing photoperiod and light intensity under high CO2 conditions is not the most appropriate means to maximize the yield of potatoes.

Description: A system and methodology were developed for the nondestructive qualitative and quantitative analysis of volatile emissions from hydroponically grown 'Waldmann's Green' leaf lettuce (Lactuca sativa L.). Photosynthetic photon flux (PPF), photoperiod, and temperature were automatically controlled and monitored in a growth chamber modified for the collection of plant volatiles. The lipoxygenase pathway products (Z)-3-hexenal, (Z)-3-hexenol, and (Z)-3-hexenyl acetate were emitted by lettuce plants after the transition from the light period to the dark period. The volatile collection system developed in this study enabled measurements of volatiles emitted by intact plants, from planting to harvest, under controlled environmental conditions.

Description: The flow through a plane asymmetric diffuser is a good test case for assessing the capability of LES since it contains features such as large scale unsteady separation and strong intermittency which are difficult to capture using conventional modeling approaches. Previous attempts to simulate this flow (Kaltenbach, 1994) have significantly underpredicted the extent of separation. The objective of the present research is to understand why the previous simulations did not predict the flow separation correctly. This study focuses on mesh refinement and matching of the inlet velocity profile. In order to perform this study, the flow solver of Kaltenbach (1994) was modified to increase its accuracy and efficiency. The improved algorithm allows for better resolution at affordable CPU cost.

Description: The objective of the present research is to assess the usefulness of large-eddy simulation (LES) methodology for flows in complex geometries. Flow past a circular cylinder has been calculated using a central-difference based solver, and the results have been compared to those obtained by a solver that employs higher-order upwind biased schemes (Beaudan & Moin, 1994). This comparison allows us to assess the suitability of these schemes for LES in complex geometry flows.

Description: Accurate SGS models must be capable of correctly representing the energy transfer between GS and SGS. Recent direct assessment of the energy transfer carried out using direct numerical simulation (DNS) data for wall-bounded flows revealed that the energy exchange is not unidirectional. Although GS kinetic energy is transferred to the SGS (forward scatter (F-scatter) on average, SGS energy is also transferred to the GS. The latter energy exchange (backward scatter (B-scatter) is very significant, i.e., the local energy exchange can be backward nearly as often as forward and the local rate of B-scatter is considerably higher than the net rate of energy dissipation.

Description: An efficient scheme for the direct numerical simulation of 3D transitional and developed turbulent flow is presented. Explicit and implicit time integration schemes for the compressible Navier-Stokes equations are compared. The nonlinear system resulting from the implicit time discretization is solved with an iterative method and accelerated by the application of a multigrid technique. Since we use central spatial discretizations and no artificial dissipation is added to the equations, the smoothing method is less effective than in the more traditional use of multigrid in steady-state calculations. Therefore, a special prolongation method is needed in order to obtain an effective multigrid method. This simulation scheme was studied in detail for compressible flow over a flat plate. In the laminar regime and in the first stages of turbulent flow the implicit method provides a speed-up of a factor 2 relative to the explicit method on a relatively coarse grid. At increased resolution this speed-up is enhanced correspondingly.

Description: The role of unsteady flow processes in establishing the performance of axial flow turbomachinery was investigated. The development of the flow model, as defined by the time average flow equations associated with the passage of the blade row embedded in a multistage configuration, is presented. The mechanisms for unsteady momentum and energy transport is outlined. The modeling of the unsteady momentum and energy transport are discussed. The procedure for simulating unsteady multistage turbomachinery flows is described.

Description: A mathematical model for closing or mathematically completing the system of equations is proposed. The model describes the time average flow field through the blade passages of multistage turbomachinery. These average-passage equation systems govern a conceptual model useful in turbomachinery aerodynamic design and analysis. The closure model was developed to insure a consistency between these equations and the axisymmetric through-flow equations. The closure model was incorporated into a calculation code for use in the simulation of the flow field about a high-speed counter rotating propeller and a high-speed fan stage.

Description: The objective of this study is to provide careful qualification and detailed measurements in a re-creation of the Obi experiment. The work will include extensive documentation of the flow two-dimensionality and detailed measurements required for testing of flow computations. Also important to this study is the close interaction of the experimental and computational groups to improve the utility of the data obtained and the accuracy of computation.

Description: The control of vortical flows is gaining significance in the design of aeronautical and marine structures. While passive devices have been used effectively in the past, active control strategies have the potential of allowing a leap in the performance of future configurations. The efficiency of control schemes is strongly dependent on the development of accurate flow models that can be devised using information that is available not only from numerical solutions of the governing Navier-Stokes equations but also can be measured experimentally. In that context it is desirable to construct adaptive control schemes using information that can be measured at the wall.

Description: Stability of the elliptic vortex attracted interest in the past decade. Cambon (1982), and Cambon, Teissedre and Jeandel (1985) have studied the stability of such flows with spatially uniform velocity gradient, and have provided RDT solutions for a wide range of the parameter S/Omega (where the strain rate S and the vorticity 2 Omega define the velocity gradient matrix). The range studied included those of hyperbolic streamlines (strain dominated, S/Omega is greater than 1), linear streamlines (simple shear, S/Omega = 1), and elliptical streamlines (vorticity dominated, S/Omega is less than 1). The latter class has more recently attracted interest and several studies appeared (Pierrehumbert 1986, Bayly 1986, Craik and Criminale 1986). These studies will be collectively referred to as PBCC.

Description: The objective is the development of a new theory which enables the algorithmic computation of all self-similar mean velocity profiles. The theory is based on Liegroup analysis and unifies a large set of self-similar solutions for the mean velocity of stationary parallel turbulent shear flows. The results include the logarithmic law of the wall, an algebraic law, the viscous sublayer, the linear region in the middle of a Couette flow and in the middle of a rotating channel flow, and a new exponential mean velocity profile not previously reported. Experimental results taken in the outer parts of a high Reynolds number flat-plate boundary layer, strongly support the exponential profile. From experimental as well as from DNS data of a turbulent channel flow the algebraic scaling law could be confirmed in both the center region and in the near wall region. In the case of the logarithmic law of the wall, the scaling with the wall distance arises as a result of the analysis and has not been assumed in the derivation. The crucial part of the derivation of all the different mean velocity profiles is to consider the invariance of the equation for the velocity fluctuations at the same time as the invariance of the equation for the velocity product equations. The latter is the dyad product of the velocity fluctuations with the equation for the velocity fluctuations. It has been proven that all the invariant solutions are also consistent with similarity of all velocity moment equations up to any arbitrary order.

Description: Nonpremixed combustion is the process whereby fuel and oxidizer species, which are each nonflammable in isolation, concurrently mix to burn a flammable mixture, and chemically react in the flammable mixture. In cases of practical industrial interest, the bulk of nonpremixed combustion occurs in a turbulent mixing regime where enhanced mass transfer rates flow the maximum power density to be achieved in any given thermochemical device.

Description: Direct numerical simulation of turbulent reacting flows places extreme demands on computational resources. At the present time, simulations can be performed only for greatly simplified reaction systems and for very low Reynolds numbers. Direct simulation of more realistic cases occurring at higher Reynolds number and including multiple species and numerous chemical reactions will exceed available computational resources far into the future. Because of this, there is a clear need to develop the technique of large eddy simulation for reacting flows. Unfortunately this task is complicated by the fact that combustion arises from chemical reactions that occur at the smallest scales of the flow. Capturing the large-scale behavior without resolving the small-scale details is extremely difficult in combustion problems. Thus LES modeling for turbulent combustion encounters difficulties not present in modeling momentum transport, in which the main effect of the small scales is to provide dissipation. The difficulty is more pronounced in premixed combustion, where detailed chemistry plays an essential role in determining the flame speed (or overall burning rate); in nonpremixed combustion infinite rate chemistry can be assumed, eliminating the small scale features to a first approximation.

Description: Many industrial applications in such fields as aeronautical, mechanical, thermal, and environmental engineering involve complex turbulent flows containing global separations and subsequent reattachment zones. Accurate prediction of this phenomena is very important because separations influence the whole fluid flow and may have an even bigger impact on surface heat transfer. In particular, reattaching flows are known to be responsible for large local variations of the local wall heat transfer coefficient as well as modifying the overall heat transfer. For incompressible, non-buoyant situations, the fluid mechanics have to be accurately predicted in order to have a good resolution of the temperature field.

Description: In simple flows, where the mean deformation rates are mild and the turbulence has time to come to equilibrium with the mean flow, the Reynolds stresses are determined by the applied strain rate. Hence in these flows, it is often adequate to use an eddy-viscosity representation. The modern family of kappa-epsilon models has been very useful in predicting near equilibrium turbulent flows, where the rms deformation rate S is small compared to the reciprocal time scale of the turbulence (epsilon/kappa). In modern engineering applications, turbulence models are quite often required to predict flows with very rapid deformations (large S kappa/epsilon). In these flows, the structure takes some time to respond and eddy viscosity models are inadequate. The response of turbulence to rapid deformations is given by rapid distortion theory (RDT). Under RDT the nonlinear effects due to turbulence-turbulence interactions are neglected in the governing equations, but even when linearized in this fashion, the governing equations are unclosed at the one-point level due to the non-locality of the pressure fluctuations.

Description: Direct numerical simulation (DNS) was used to study homogeneous turbulence with a mean flow having elliptic streamlines. This flow combines the effects of rotation and strain on the turbulence. There are many important reasons for studying the elliptic streamline flow. This flow contains the effects of both rotation and strain and is therefore similar to the mean flow in a vortex strained in the plane perpendicular to its axis. Such flows provide insight into the fundamental vortical interactions within turbulence. A strained vortex also occurs in airplane wakes, in which each wingtip vortex induces a strain field on the other. The strain field can affect the stability of these vortices and thereby their turbulent structure downstream. The ability to understand and predict the turbulent structure of the vortices is important to the wake hazard problem which is of major concern for the safety of commercial aircraft.

Description: A combined experimental and computational effort to test the performance of turbulence models in prediction of three-dimensional separation is under way. The current phase pf the project is the acquisition of flow field data using laser Doppler velocimetry. Accomplishments for summer 1996 include development of the instrument and signal processing system and acquisition of mean velocity profiles.

Description: The objective of this study is to investigate the effects of cold only, commercially available gloves only, and the combination of gloves and cold on the blood flow and surface (skin) temperature of the medial and proximal phalanxes of digit 3, the metacarpal region of the hand, and the forearm.

Description: The mean and phase averaged pressure and wall shear stress distributions were measured on the stator wall of a 50% eccentric annular seal which was whirling in a circular orbit at the same speed as the shaft rotation. The shear stresses were measured using flush mounted hot-film probes. Four different operating conditions were considered consisting of Reynolds numbers of 12,000 and 24,000 and Taylor numbers of 3,300 and 6,600. At each of the operating conditions the axial distribution (from Z/L = -0.2 to 1.2) of the mean pressure, shear stress magnitude, and shear stress direction on the stator wall were measured. Also measured were the phase averaged pressure and shear stress. These data were combined to calculate the force distributions along the seal length. Integration of the force distributions result in the net forces and moments generated by the pressure and shear stresses. The flow field inside the seal operating at a Reynolds number of 24,000 and a Taylor number of 6,600 has been measured using a 3-D laser Doppler anemometer system. Phase averaged wall pressure and wall shear stress are presented along with phase averaged mean velocity and turbulence kinetic energy distributions located 0.16c from the stator wall where c is the seal clearance. The relationships between the velocity, turbulence, wall pressure and wall shear stress are very complex and do not follow simple bulk flow predictions.

Description: Computational programs developed for the thermal analysis of tilting and flexure-pad hybrid bearings, and the unsteady flow and transient response of a point mass rotor supported on fluid film bearings are described. The motion of a cryogenic liquid on the thin film annular region of a fluid film bearing is described by a set of mass and momentum conservation, and energy transport equations for the turbulent bulk-flow velocities and pressure, and accompanied by thermophysical state equations for evaluation of the fluid material properties. Zeroth-order equations describe the fluid flow field for a journal static equilibrium position, while first-order (linear) equations govern the fluid flow for small amplitude-journal center translational motions. Solution to the zeroth-order flow field equations provides the bearing flow rate, load capacity, drag torque and temperature rise. Solution to the first-order equations determines the rotordynamic force coefficients due to journal radial motions.

Description: Turbulent flows that exhibit inhomogeneities in the streamwise direction pose a particular challenge to numerical simulation approaches due to the need to prescribe time-dependent turbulent inflow conditions. In most cases the flow downstream is more or less 'driven' by the conditions at the inlet, making it necessary to specify realistic turbulent fluctuations that are in equilibrium with the assumed mean flow. This requirement often dictates that the inflow data should satisfy the Navier-Stokes equations, which in turn implies that an independent simulation be used to generate them. Detailed simulations for the purpose of creating inflow conditions can be costly and thus certain levels of approximation are desirable. In this paper we shall focus on an approximate yet accurate method for generating inflow conditions for spatially-developing boundary layer simulations. The proposed method is essentially a simplification of the method of Spalart and Leonard (1985), who devised an ingenious transformation that allows for the calculation of spatially evolving boundary layers in conjunction with periodic boundary conditions applied in the streamwise direction. While this method is elegant and highly accurate, it is more complicated than is necessary for the purpose of generating inflow data. A few key approximations are used in this work to arrive at a 'modified Spalart method' that is very easy to implement and efficient to use. The new method is shown to yield results that compare well with the computations of Spalart (1988). When used as a means of generating inflow data, the modified Spalart method is shown to be superior to existing approaches.

Description: Accurate large eddy simulation (LES) of a wall-bounded flow generally requires a near-wall resolution comparable to that in direct numerical simulation (DNS). As much as 50% of the total grid points and computational costs are expended in the near-wall regions in a typical simulation. This limits LES to fairly low Reynolds numbers on current computers. To perform practical flow applications at realistically high Reynolds numbers, such as flow over an airfoil, it is desirable to replace very thin, near-wall regions in the LES with easily and inexpensively computed wall models to specify the near-wall boundary conditions.

Description: Reducing the infinite degrees of freedom inherent in fluid motion into a manageable number of modes to analyze fluid motion is presented. The concepts behind the center manifold technique are used. Study of the Blasius boundary layer and a precise description of stability within the flow field are discussed.

Description: Regulations require oxygen usage by commercial airliners during check out and during certain aircraft configurations. This oxygen is drawn from a high pressure on-board cylinder storage system. In a typical aircraft, oxygen cylinder removal for oxygen ground servicing is conducted every 4 to 6 weeks. In the early 1990's, it was recognized that an on-board oxygen generating system (OBOGS) could provide an economic advantage for the airlines. An in-flight service evaluation (ISE) of the SPE-OBOGS by United Technologies Corporate is in the planning stage.

Description: A key ingredient in the simulation of flow in porous media is the accurate determination of the velocities that drive the flow. The large scale irregularities of the geology, such as faults, fractures, and layers suggest the use of irregular grids in the simulation. Work has been done in applying the finite volume element (FVE) methodology as developed by McCormick in conjunction with mixed methods which were developed by Raviart and Thomas. The resulting mixed finite volume element discretization scheme has the potential to generate more accurate solutions than standard approaches. The focus of this paper is on a multilevel algorithm for solving the discrete mixed FVE equations. The algorithm uses a standard cell centered finite difference scheme as the 'coarse' level and the more accurate mixed FVE scheme as the 'fine' level. The algorithm appears to have potential as a fast solver for large size simulations of flow in porous media.

Description: The effect of the unsteady flow field in a multistage compressor on the time-averaged performance was assessed. The energy transport by the unsteady deterministic flow field was taken into account. The magnitude of the body force resulting from the aerodynamic loading on a blade row was compared to the gradient of the stress tensor associated with the unsteady time-resolved flow field generated by the blade row. The magnitude of the work performed by these forces was compared to the divergence of the energy correlations produced by the unsteady time-resolved flow field. The stress tensor and the energy correlations are non-negligible in the end wall regions. The results suggest that the turbulence is the primary source of flow mixing in the midspan region.

Description: Turbulent wakes are known to develop self-similarly sufficiently far downstream from obstacles that generate them. It has long been assumed that the spreading rate of the wake in the self-similar regime is independent of the details of the body generating the wake, being dependent only on the total drag (or momentum deficit). This assumption seems to be in contradiction with some recent experiments. In this study we attempt to complement these experimental investigations through a numerical study of a time-developing wake. A numerical study has the advantage of eliminating many of the uncontrolled factors present in experiments and allowing precise control of initial conditions. Large-eddy simulations employing the recently developed dynamic localization model are used to extend previous results from direct numerical simulations. The large-eddy simulation results are compared to the direct numerical simulation database, wherever such comparisons are feasible, as a check of the method. Like the experiments, the large-eddy simulations suggest that non-unique self-similar states, characterized by different spreading rates and turbulent statistics, are possible and that they can be maintained for significant time periods. The study also demonstrates the predictive capability of the dynamic localization subgrid model.

Description: Turbulence modeling plays an important role in the study of high-speed flows in engineering and aerodynamic problems; they include flows in supersonic combustion engines and over hypersonic transport aircraft. The enhancement of the kinetic energy dissipation by the dilatational terms is one of the typical compressibility effects. Zeman (1990) and Sarkar et al. (1991) proposed that the dilatation dissipation is proportional to the solenoidal dissipation and is a function of the turbulent Mach number. Sarkar (1992) also modeled the pressure-dilatation correlation using the turbulent Mach number. Zeman (1991) related the correlation to the rate of change of the pressure variance.

Description: Combustion is generally categorized as either premixed, where flames propagate into homogeneous mixtures of reactants, or as nonpremixed, where initially separated reactants diffuse into the reaction zones. Although these approaches are applicable to many combustion devices, there are cases not in either of these two limiting regimes. Under such circumstances, one must consider partially premixed combustion. In partially premixed combustion, mechanisms from both premixed and nonpremixed regimes coexist and, as a result, some interesting phenomena arise. One such phenomenon is flame stabilization in laminar mixing layers by triple flames.

Description: A new concept, Multiple Scale Simulation (MSS), is presented in this paper. The basic idea is that the flow is decomposed into several component groups according to spatial and temporal length scales. Each group has its own subdomain, governing system, mesh size, and discretization method. The simulation is then performed groupwise. This approach has been successfully applied in combination with the intergrid dissipation technique for simulation of transitional and turbulent flow in 3-D boundary layers, and it is feasible for 3-D airfoils and other more complex configurations. MSS should prove to ameliorate the scale problems associated with conventional direct numerical simulation.

Description: The flow in a hydrostatic journal bearing (HJB) is described by a mathematical model that uses the three dimensional non-orthogonal form of the Navier-Stokes equations. Using the u, v, w, and p, as primary variables, a conservative formulation, finite volume multi-block method is applied through a collocated, body fitted grid. The HJB has four shallow pockets with a depth/length ratio of 0.067. This paper represents a natural extension to the two and three dimensional studies undertaken prior to this project.

Description: A steady, three dimensional average-passage equation system was derived. The purpose was to simulate multistage turbomachinery flows. These equations describe a steady, viscous flow that is periodic from blade passage to blade passage. Moreover, these equations have a closure problem that is similar to that of the Reynolds-average Navier-Stokes equations. A scaled form of the average-passage equation system could provide an improved mathematical model for simulating the flow in the design and in the off-design conditions of a multistage machine.

Description: CFD has played a major role in the resurgence of hypersonic flight, on the premise that numerical methods will allow us to perform simulations at conditions for which no ground test capability exists. Validation of CFD methods is being established using the experimental data base available, which is below Mach 8. It is important, however, to realize the limitations involved in the extrapolation process as well as the deficiencies that exist in numerical methods at the present time. Current features of CFD codes are examined for application to propulsion system components. The shortcomings in simulation and modeling are identified and discussed.

Description: Implementation issues associated with the application of Krylov subspace iterative methods, such as Newton-GMRES, are presented within the framework of practical computational fluid dynamic (CFD) applications. This paper categorizes, evaluates, and contrasts the major ingredients (function evaluations, matrix-vector products, and preconditioners) of Newton-GMRES Krylov subspace methods in terms of their effect on the local linear and global nonlinear convergence, memory requirements, and accuracy. The discussion focuses on Newton-GMRES in both a structured multi-zone incompressible Navier-Stokes solver and an unstructured mesh finite-volume Navier-Stokes solver. Approximate versus exact matrix-vector products, effective preconditioners, and other pertinent issues are addressed.

Description: This paper reviews recent research conducted in the Flight Management and Human Factors Division of NASA Ames Research Center on superimposed symbology (as found on HUDs and HMDs). We first identify various performance problems which suggest that superimposed symbology impairs pilots' ability to maintain simultaneous awareness of instrument information and information in the forward visual scene. Results of experiments supporting an attentional account of the impairment are reported. A design solution involving the concept of 'scene-linked' symbology is developed, and experiments testing the design solution are reported. An application of the scene-linking concept, in the form of a candidate HUD to support ground taxi operations for civil transport, is described.

Description: A human's ability to perform physical tasks is limited by physical strength, not by intelligence. We define "extenders" as a class of robot manipulators worn by humans to augment human mechanical strength, while the wearer's intellect remains the central control system for manipulating the extender. Our research objective is to determine the ground rules for the design and control of robotic systems worn by humans through the design, construction, and control of several prototype experimental direct-drive/non-direct-drive multi-degree-of-freedom hydraulic/electric extenders. The design of extenders is different from the design of conventional robots because the extender interfaces with the human on a physical level. Two sets of force sensors measure the forces imposed on the extender by the human and by the environment (i.e., the load). The extender's compliances in response to such contact forces were designed by selecting appropriate force compensators. This paper gives a summary of some of the selected research efforts related to Extender Technology, carried out during 1980s. The references, at the end of this article, give detailed description of the research efforts.

Description: Ultrasonic and electrostatic levitation techniques have allowed the experimental investigation of the nonlinear oscillatory dynamics of free droplets with diameter between 0.1 and 0.4 cm. The measurement of the resonance frequencies of the first three normal modes of large amplitude shape oscillations in an electric field of varying magnitude has been carried out with and without surface charges for weakly conducting liquids in air. These oscillations of nonspherical levitated drops have been driven by either modulating the ultrasonic field or by using a time-varying electric field, and the free decay from the oscillatory state has been recorded. A decrease in the resonance frequency of the driven fundamental quadrupole mode has been measured for increasing oblate deformation in the absence of an electric field. Similarly, a decrease in this frequency has also been found for increasing DC electric field magnitude. A soft nonlinearity exists in the amplitude dependence of the resonant mode frequencies for freely decaying as well as ultrasonically and electrically driven uncharged drops. This decrease in resonance frequency is accentuated by the presence of free surface charge on the drop. Subharmonic resonance excitation has been observed for drops in a time-varying electric field, and hysteresis exists for resonant modes driven to large amplitude. Mode coupling from lower-order resonances to higher-order modes has been found to be very weak, even for fairly large amplitude shape oscillations. Most of these results are in general agreement with predictions from recent analytical and numerical investigations.

Description: The Liquid Propellant Program (LPP) computer code is a super-set of the industry standard Two Dimensional Kinetics (TDK) computer code. The TDK code uses a two dimensional method of characteristics solution with fully coupled finite rate kinetics for axially symmetric nozzles. The chemical reactions are modeled with a generalized reaction package that include three dimensional body efficiencies and four reaction rate forms. The code performs optional solutions for frozen or equilibrium flow. TDK evaluates discrete shocks, both attached or induced. The Transonic module models variable mixture ratio profiles from the combustion chamber injector. The Mass Addition Boundary Layer module (MABL) calculates the boundary parameters with the same chemistry options, and includes transpiration or tangential slot injection of gas at the wall. The LPP upgrades include: planar nozzle, scarfed nozzles, plug nozzles, and scramjet nozzle configurations. The code evaluates both upper and lower wall flow simulation, and includes the interaction with the external flow. The MABL module evaluates equilibrium radiation heat transfer for both upper and lower walls. In addition, LPP code models combustion effects due to injector inefficiencies with the Spray Combustion Analysis Program (SCAP) module. The LPP package provides extensive post plotting capabilities for flow visualization. The LPP is sufficiently fast and robust to provide performance predictions for extensive parametric studies and sufficiently accurate to provide flow field and performance solutions for detailed studies.

Description: To better understand the Space Shuttle Main Engine (SSME) startup/shutdown tansients, an airflow test of a two dimensional nozzle was conducted at Marshall Space Flight Center's trisonic wind tunnel. Photographic and other instrumentation show during an SSME start large nozzle shell distortions occur as the Mach disk is passing through the nozzle. During earlier develop of the SSME, this startup transient resulted in low cycle fatigue failure of one of the LH2 feedlines. The two dimensional SSME nozzle test was designed to measure the static and fluctuating pressure environment and color Schlieren video during the startup and shutdown phases of the run profile.

Description: The objective of the present research is to develop improved turbulence models for the computation of complex flows through turbomachinery passages, including the effects of streamline curvature, heat transfer and secondary flows. Advanced turbulence models are crucial for accurate prediction of rocket engine flows, due to existance of very large extra strain rates, such as strong streamline curvature. Numerical simulation of the turbulent flows in strongly curved ducts, including two 180-deg ducts, one 90-deg duct and a strongly concave curved turbulent boundary layer have been carried out with Reynolds stress models (RSM) and algebraic Reynolds stress models (ARSM). An improved near-wall pressure-strain correlation has been developed for capturing the anisotropy of turbulence in the concave region. A comparative study of two modes of transition in gas turbine, the by-pass transition and the separation-induced transition, has been carried out with several representative low-Reynolds number (LRN) k-epsilon models. Effects of blade surface pressure gradient, freestream turbulence and Reynolds number on the blade boundary layer development, and particularly the inception of transition are examined in detail. The present study indicates that the turbine blade transition, in the presence of high freestream turbulence, is predicted well with LRN k-epsilon models employed. The three-dimensional Navier-Stokes procedure developed by the present authors has been used to compute the three-dimensional viscous flow through the turbine nozzle passage of a single stage turbine. A low Reynolds number k-epsilon model and a zonal k-epsilon/ARSM (algebraic Reynolds stress model) are utilized for turbulence closure. An assessment of the performance of the turbulence models has been carried out. The two models are found to provide similar predictions for the mean flow parameters, although slight improvement in the prediction of some secondary flow quantities has been obtained by the ARSM model. It's found that the wake profiles inside the endwall boundary layers are predicted better than those near the mid-span.

Description: The effective performance of modular thrusters in an aerospike configuration is difficult to determine. Standard analytical tools are applicable to conventional nozzle shapes, but are limited when applied to an aerospike nozzle (An aerospike nozzle is an altitude compensating external nozzle). Three baseline nozzle shapes are derived using standard analytical procedures. The baseline nozzle sizes are restricted to fill a volume envelope. The three shapes are an axi-symmetric round nozzle, a two dimensional planar square exit nozzle, and a super elliptic round to nearly square nozzle. The integrated (thruster/aerospike) performance of the three nozzles is determined through the use of three dimensional viscous computational fluid dynamic (CFD) calculation where complex features of the flow field can be accurately captured. The resulting installed performance is then used to evaluate the efficiency of these nozzle shapes for aerospike applications. The determination of effective performance of a thruster nozzle integrated into an aerospike nozzle require the solution of the three dimensional turbulent Navier-Stokes equations. The model used in this study consisted of two zones; one of the upstream thruster cowl surface so freestream conditions can be accurately predicted, and two, the aerospike surface beginning with the thruster outflow and extending to the end of the aerospike surface. The numerical grid consisted of over 120,000 nodes and used symmetry on the thruster centerline and edge. A two species non-reacting chemistry model was used to capture the variation of fluid properties between the hot plume base and freestream air. From the results of the three baseline nozzle aerospike calculations, the effictive performance of the nozzle was determined. The flow fields of these calculations do show some variation between the cases. Recirculation zones on the cowl surface is predicted for the two dimensional planar nozzle and a smaller one for the super elliptic nozzle. The recirculation is caused by the strong pressure gradient between the plume and freestream flows. The axi-symmetric nozzle results indicate recirculation zones on the thruster face. These recirculation zones smooth the pressure gradient between the plume and freestream flow limiting the formation of recirculation on the cowl surface. Thruster to thruster interaction is evident for the axi-symmetric and supper elliptic calculation while the two dimensional planar nozzle did not have any lateral expansion in the nozzle, so thruster to thruster interaction is limited. The integrated performance results, at the altitude choosen, show very little variation between the three thruster shapes. This result allows for nozzle shape determination based on additional considerations (thermal, structural, weight) besides performance.

Description: The Hydroblast Research Cell at MSFC is both a research and a processing facility. The cell is used to investigate fundamental phenomena associated with waterjets as well as to clean hardware for various NASA and contractor projects. In the area of research, investigations are made regarding the use of high pressure waterjets to strip paint, grease, adhesive and thermal spray coatings from various substrates. Current industrial methods of cleaning often use ozone depleting chemicals (ODC) such as chlorinated solvents, and high pressure waterjet cleaning has proven to be a viable alternative. Standard methods of waterjet cleaning use hand held or robotically controlled nozzles. The nozzles used can be single-stream or multijet nozzles, and the multijet nozzles may be mounted in a rotating head or arranged in a fan-type shape. We consider in this paper the use of a rotating, multijet, high pressure water nozzle which is robotically controlled. This method enables rapid cleaning of a large area, but problems such as incomplete coverage (e.g. the formation of 'islands' of material not cleaned) and damage to the substrate from the waterjet have been observed. In addition, current stripping operations require the nozzle to be placed at a standoff distance of approximately 2 inches in order to achieve adequate performance. This close proximity of the nozzle to the target to be cleaned poses risks to the nozzle and the target in the event of robot error or the striking of unanticipated extrusions on the target surface as the nozzle sweeps past. Two key motivations of this research are to eliminate the formation of 'coating islands' and to increase the allowable standoff distance of the nozzle.

Description: Accurate prediction of all physical phenomena in a combustion chamber is essential for better understanding of the system performance. Atomization, evaporation, combustion, chemical kinetics, and turbulence are those processes of great importance that need to be well understood. Processes involving the liquid phase in a combustion chamber will be further complicated under supercritical conditions. More advanced and accurate numerical techniques are required to extend our understanding of the above phenomena. A computer program for multi-species/multi-phase flow was developed for NASA/MSFC in 1992. This code, called Liquid Thrust Chamber Performance (LTCP) program takes an Eulerian- Eulerian approach and is based on the Total Variation Diminishing (TVD) technique with Lax-Friedrichs upwind method. Under the NASA/ASEE SFFP the LTCP code was used to predict the performance characteristics of several engines that were of particular interest to NASA. This code was also successful in a combustion detonation study. Converting the program to the PC platform was accomplished which extends usability and makes it available to a wider range of users. The Eulerian formulation of the liquid phase provides a suitable model that can be extended to include combustion modeling under supercritical conditions. The results have been compared against the ones of other codes and available measured data. The algorithm proved to be robust and efficient for problems with stiff source terms.

Description: One of the more mission-critical tasks performed in space is extravehicular activity (EVA) which requires the astronaut to be external to the station or spacecraft, and subsequently at risk from the many threats posed by space. These threats include, but are not limited to: no significant atmosphere, harmful electromagnetic radiation, micrometeoroids, and space debris. To protect the astronaut from this environment, a special EVA suit is worn which is designed to maintain a sustainable atmosphere (at 1/3 atmosphere) and provide protection against the hazards of space. While the EVA suit serves these functions well, it does impose limitations on the astronaut as a consequence of the safety it provides. Since the astronaut is in a virtual vacuum, any atmospheric pressure inside the suit serves to pressurize the suit and restricts mobility of flexible joints (such as fabric). Although some of the EVA suit joints are fixed, rotary-style joints, most of the mobility is achieved by the simple flexibility of the fabric. There are multiple layers of fabric, each of which serves a special purpose in the safety of the astronaut. These multiple layers add to the restriction of motion the astronaut experiences in the space environment. Ground-based testing is implemented to evaluate the capability of EVA-suited astronauts to perform the various tasks in space. In addition to the restriction of motion imposed by the EVA suit, most EVA activity is performed in a micro-gravity (weight less) environment. To simulate weightlessness EVA-suited testing is performed in a neutral buoyancy simulator (NBS). The NBS is composed of a large container of water (pool) in which a weightless environment can be simulated. A subject is normally buoyant in the pressurized suit; however he/she can be made neutrally buoyant with the addition of weights. In addition, most objects the astronaut must interface with in the NBS sink in water and flotation must be added to render them "weightless". The implementation of NBS testing has proven to invaluable in the assessment of EVA activities performed with the Orbiter and is considered to be a key step in the construction of the International Space Station (ISS). While the NBS testing is extremely valuable, it does require considerable overhead to maintain and operate. It has been estimated that the cost of utilizing the facility is approximately $10,000 per day. Therefore it is important to maximize the utility of NBS testing for optimal results. One important aspect to consider in any human/worksite interface is the considerable wealth of anthropometric and ergonomic data available. A subset of this information specific to EVA activity is available in NASA standard 3000. The difficulty in implementing this data is that most of the anthropometric information is represented in a two-dimensional format. This poses some limitations in complete evaluation of the astronaut's capabilities in a three-dimensional environment. Advances in computer hardware and software have provided for three-dimensional design and implementation of hardware with the advance of computer aided design (CAD) software. There are a number of CAD products available and most companies and agencies have adopted CAD as a fundamental aspect of the design process. Another factor which supports the use of CAD is the implementation of computer aided manufacturing (CAM) software and hardware which provides for rapid prototyping and decreases the time to product in the design process. It is probable that most hardware to be accessed by astronauts in EVA or IVA (intravehicular activity) has been designed by a CAD system, and is therefore represented in three-dimensional space for evaluation. Because of the implementation of CAD systems and the movement towards early prototyping, a need has arisen in industry and government for tools which facilitate the evaluation of ergonomic consideration in a three-dimensional environment where the hardware has been designed by the CAD tools. One such product is Jack which was developed by the University of Pennsylvania with funding from several government agencies, including NASA. While the primary purpose of Jack is to model human figures in a ground-based (gravity) environment, it can be utilized to evaluate EVA-suited activities as well. The effects of simulated gravity must be turned off by turning off "behaviors". Although Jack provides human figures for manipulation, the primary instrument to be evaluated for EVA mobility is the work envelope provided by the EVA suit. An EVA Jack suit model has been developed by NASA-JSC and was utilized in this study. This suit model provided a more restrictive motion environment as expected for an EVA suited subject. As part of this study, the anthropometric dimensions for a 50th percentile male were compared with basic anthropometric data and were found to be representative for the population group expected in the NASA flight program. The joints for the suit were created in a manner which provided consistent performance with EVA reach envelopes published in NASA standard #3000.

Description: This is a preliminary study devoted to verifying whether or not direct simulations of turbulent Magneto-Hydro-Dynamic (MHD) flows in liquid metals reproduce experimental observations of drag reduction. Two different cases have been simulated by a finite difference scheme which is second order accurate in space and time. In the first case, an external azimuthal magnetic field is imposed. In this case, the magnetic field acts on the mean axial velocity and complete laminarization of the flow at N(sub a) = 30 has been achieved. In the second case, an axial magnetic field is imposed which affects only fluctuating velocities, and thus the action is less efficient. This second case is more practical, but comparison between numerical and experimental results is only qualitative.

Description: Direct numerical simulations are performed for the elliptic streamline flow, which is a homogeneous turbulent flow that combines the effects of solid body rotation and strain. Simulations are run over a range of parameters in order to determine the effect of changing rotation and strain separately. For early times the nonlinear cascade is suppressed, but then is re-established at later times. The growth rate of turbulent kinetic energy agrees at early times with the trends from linear theory, but at later times the flow seems to approach an asymptotic state that is independent of the ratio of mean flow rotation rate to strain rate. A comparison with standard Reynolds stress turbulence models is made. It is found that for strong rotation rates, the models predict decay of the turbulence, while the simulations show exponential growth. Close examination of the simulation results shows that they are affected by excessively low Reynolds numbers. Suggestions for reducing low Reynolds number effects in future simulations is given.

Description: We report preliminary results obtained from the large eddy simulation of a backward facing step at a Reynolds number of 5100. The numerical platform is based on a high order Legendre spectral element spatial discretization and a least squares time integration scheme. A non-reflective outflow boundary condition is in place to minimize the effect of downstream influence. Smagorinsky model with Van Driest near wall damping is used for sub-grid scale modeling. Comparisons of mean velocity profiles and wall pressure show good agreement with benchmark data. More studies are needed to evaluate the sensitivity of this method on numerical parameters before it is applied to complex engineering problems.

Description: A wall pressure frequency spectrum model (Blake 1971 ) has been evaluated using databases from Direct Numerical Simulations (DNS) of a turbulent boundary layer (Na & Moin 1996). Good agreement is found for moderate to strong adverse pressure gradient flows in the absence of separation. In the separated flow region, the model underpredicts the directly calculated spectra by an order of magnitude. The discrepancy is attributed to the violation of the model assumptions in that part of the flow domain. DNS computed coherence length scales and the normalized wall pressure cross-spectra are compared with experimental data. The DNS results are consistent with experimental observations.

Description: A variable explicit/implicit characteristics-based advection scheme that is second-order accurate in space and time has been developed recently for unstructured deforming meshes (O'Rourke & Sahota 1996a). To explore the suitability of this methodology for Large-Eddy Simulation (LES), three subgrid-scale turbulence models have been implemented in the CHAD CFD code (O'Rourke & Sahota 1996b): a constant-coefficient Smagorinsky model, a dynamic Smagorinsky model for flows having one or more directions of statistical homogeneity, and a Lagrangian dynamic Smagorinsky model for flows having no spatial or temporal homogeneity (Meneveau et al. 1996). Computations have been made for three canonical flows, progressing towards the intended application of in-cylinder flow in a reciprocating engine. Grid sizes were selected to be comparable to the coarsest meshes used in earlier spectral LES studies. Quantitative results are reported for decaying homogeneous isotropic turbulence, and for a planar channel flow. Computations are compared to experimental measurements, to Direct-Numerical Simulation (DNS) data, and to Rapid-Distortion Theory (RDT) where appropriate. Generally satisfactory evolution of first and second moments is found on these coarse meshes; deviations are attributed to insufficient mesh resolution. Issues include mesh resolution and computational requirements for a specified level of accuracy, analytic characterization of the filtering implied by the numerical method, wall treatment, and inflow boundary conditions. To resolve these issues, finer-mesh simulations and computations of a simplified axisymmetric reciprocating piston-cylinder assembly are in progress.

Description: Subgrid-scale models for Large Eddy Simulation (LES) in both the velocity-pressure and the vorticity-velocity formulations were evaluated and compared in a priori tests using spectral Direct Numerical Simulation (DNS) databases of isotropic turbulence: 128(exp 3) DNS of forced turbulence (Re(sub(lambda))=95.8) filtered, using the sharp cutoff filter, to both 32(exp 3) and 16(exp 3) synthetic LES fields; 512(exp 3) DNS of decaying turbulence (Re(sub(Lambda))=63.5) filtered to both 64(exp 3) and 32(exp 3) LES fields. Gaussian and top-hat filters were also used with the 128(exp 3) database. Different LES models were evaluated for each formulation: eddy-viscosity models, hyper eddy-viscosity models, mixed models, and scale-similarity models. Correlations between exact versus modeled subgrid-scale quantities were measured at three levels: tensor (traceless), vector (solenoidal 'force'), and scalar (dissipation) levels, and for both cases of uniform and variable coefficient(s). Different choices for the 1/T scaling appearing in the eddy-viscosity were also evaluated. It was found that the models for the vorticity-velocity formulation produce higher correlations with the filtered DNS data than their counterpart in the velocity-pressure formulation. It was also found that the hyper eddy-viscosity model performs better than the eddy viscosity model, in both formulations.

Description: A general discussion on the structure of the eddy viscosity tensor in anisotropic flows is presented. The systematic use of tensor symmetries and flow symmetries is shown to reduce drastically the number of independent parameters needed to describe the rank 4 eddy viscosity tensor. The possibility of using Onsager symmetries for simplifying further the eddy viscosity is discussed explicitly for the axisymmetric geometry.

Description: In numerical simulation of complex flows, it is important to identify different length scales of the flow and treat them differently. In this report, we introduce a new multilevel scheme for simulating turbulent channel flows. Two different versions of the scheme, namely the spectral and finite difference versions, are presented. The spectral version of the scheme is based on a spectral-Galerkin formulation which provides a natural decomposition of the flow into small and large wavelength parts, and which leads to linear systems that can be solved with quasi-optimal computational complexity. In the finite difference version, the Incremental Unknown (IU) is used to separate the length scales. Preliminary numerical results indicate that the scheme is well suited for turbulence computations and provides results which are comparable to that by Direct Numerical Simulation (DNS) but with significantly less CPU time.

Description: The possibility of using the information from simultaneous equivalent Large Eddy Simulations (LAS) for improving the subgrid scale modeling is investigated. An ensemble average dynamic model is proposed as an alternative to the usual spatial average versions. It is shown to be suitable independently of the existence of any homogeneity directions, and its formulation is thus universal. The ensemble average dynamic model is shown to give very encouraging results for as few as 16 simultaneous LES's.

Description: The equations describing the first two terms of an asymptotic expansion of the solution of the planar turbulent mixing layer for values of the velocity ratio close to one are obtained. The first term of this expansion is the solution of the well-known time-evolving problem and the second, which includes the effects of the increase of the turbulence scales in the stream-wise direction, obeys a linear system of equations. Numerical solutions of these equations for a two-dimensional reacting mixing layer show that the correction to the time-evolving solution may explain the asymmetry of the entrainment and the differences in product generation observed in flip experiments.

Description: A four-equation model of turbulence is applied to the numerical simulation of flows with massive separation induced by a sudden expansion. The model constants are a function of the flow parameters, and two different formulations for these functions are tested. The results are compared with experimental data for a high Reynolds-number case and with experimental and DNS data for a low Reynolds-number case. The computations prove that the recovery region downstream of the massive separation is properly modeled only for the high Re case. The problems in this case stem from the gradient diffusion hypothesis, which underestimates the turbulent diffusion.

Description: The beta-assumed-pdf approximation of (Cook & Riley 1994) is tested as a subgrid model for the LES computation of nonpremixed turbulent reacting flows, in the limit of cold infinitely fast chemistry, for two plane turbulent mixing layers with different degrees of intermittency. Excellent results are obtained for the computation of integrals properties such as product mass fraction, and the model is applied to other quantities such as powers of the temperature and the pdf of the scalar itself. Even in these cases the errors are small enough to be useful in practical applications. The analysis is extended to slightly out of equilibrium problems such as the generation of radicals, and formulated in terms of the pdf of the scalar gradients. It is shown that the conditional gradient distribution is universal in a wide range of cases whose limits are established. Within those limits, engineering approximations to the radical concentration are also possible. It is argued that the experiments in this paper are essentially in the limit of infinite Reynolds number.

Description: A Second Moment Closure computation (SMC) is compared in detail with the Direct Numerical Simulation (DNS) data of Le and Moin for the backstep flow at Re = 5,000 in an attempt to understand why the intensity of the backflow and, consequently, the friction coefficient in the recirculation bubble are severely underestimated. The data show that this recirculation bubble is far from being laminar except in the very near wall layer. A novel 'differential a priori' procedure was used, in which the full transport equation for one isolated component of the Reynolds stress tensor was solved using DNS data as input. Conclusions are then different from what would have been deduced by comparing a full simulation to a DNS. One cause of discrepancy was traced back to insufficient transfer of energy to the normal stress by pressure strain, but was not cured. A significant finding, confirmed by the DNS data in the core region of a channel flow, is that the coefficient that controls destruction of dissipation, C epsilon(sub 2), should be decreased by a factor of 2 when production is vanishing. This is also the case in the recirculation bubble, and a new formulation has cured 25% of the backflow discrepancy.

Description: New models are proposed for the 'slow and 'rapid' parts of the pressure diffusive transport based on the examination of DNS databases for plane mixing layers and wakes. The model for the 'slow' part is non-local, but requires the distribution of the triple-velocity correlation as a local source. The latter can be computed accurately for the normal component from standard gradient diffusion models, but such models are inadequate for the cross component. More work is required to remedy this situation.

Description: A new approach to Reynolds averaged turbulence modeling is proposed which has a computational cost comparable to two equation models but a predictive capability approaching that of Reynolds stress transport models. This approach isolates the crucial information contained within the Reynolds stress tensor, and solves transport equations only for a set of 'reduced' variables. In this work, Direct Numerical Simulation (DNS) data is used to analyze the nature of these newly proposed turbulence quantities and the source terms which appear in their respective transport equations. The physical relevance of these quantities is discussed and some initial modeling results for turbulent channel flow are presented.

Description: A second order effect associated with high intensity sound field, acoustic streaming has been historically investigated to gain a fundamental understanding of its controlling mechanisms and to apply it to practical aspects of heat and mass transfer enhancement. The objectives of this new research project are to utilize a unique experimental technique implementing ultrasonic standing waves in closed cavities to study the details of the generation of the steady-state convective streaming flows and of their interaction with the boundary of ultrasonically levitated near-spherical solid objects. The goals are to further extend the existing theoretical studies of streaming flows and sample interactions to higher streaming Reynolds number values, for larger sample size relative to the wavelength, and for a Prandtl and Nusselt numbers parameter range characteristic of both gaseous and liquid host media. Experimental studies will be conducted in support to the theoretical developments, and the crucial impact of microgravity will be to allow the neglect of natural thermal buoyancy. The direct application to heat and mass transfer in the absence of gravity will be emphasized in order to investigate a space-based experiment, but both existing and novel ground-based scientific and technological relevance will also be pursued.

Description: It is well known that radiative heat transport influences many types of buoyant flows due to its effect on the temperature and thus density field in the fluid medium. It is of interest to study gaseous flows driven solely by radiation in the absence of buoyancy, particularly because of its application to astrophysical flows that are well known from astronomical observations and numerical simulation. However, no laboratory-scale experiments of this phenomenon have ever been conducted. To study the possibility of obtaining such flows in the laboratory, an apparatus was built to produce large temperature differences (Delta (T)) up to 300 K in a gas confined between flat parallel plates. SF6 was used as the radiatively-active gas because its Planck absorption length is much shorter than that of any other common non-reactive gas. The NASA-Lewis 2.2 second drop tower was used to obtain reduced gravity in order to suppress buoyancy effects. To image the resulting flows, a laser shearing interferometer was employed. Initial results indicate the presence of flow that does not appear to be attributable to the residual flow resulting from buoyancy influences before the drop. For Delta(T) greater than 70 K, slight deformations in the interferometer fringes seen at lower Delta(T) became large unsteady swirls. Such behavior did not occur for radiatively-inactive gases, suggesting that a flow driven solely by radiation was obtained in SF6 and to a lesser extent in CO2. This was more pronounced at higher pressures and plate spacings, consistent with our scaling predictions.

Description: The first objective of this experiment is to build an experimental system in which, in analogy to a geophysical system, a compressible fluid in a spherical annulus becomes radially stratified in density through an A.C. electric field. When this density gradient is demonstrated, the system will be augmented so that the fluid can be driven by heating and rotation and tested in preparation for a microgravity experiment. This apparatus consists of a spherical capacitor filled with critical fluid in a temperature controlled environment. To make the fluid critical, the apparatus will be operated near the critical pressure, critical density, and critical temperature of the fluid. This will result in a highly compressible fluid because of the properties of the fluid near its critical point. A high voltage A.C. source applied across the capacitor will create a spherically symmetric central force because of the dielectric properties of the fluid in an electric field gradient. This central force will induce a spherically symmetric density gradient that is analogous to a geophysical fluid system. To generate such a density gradient the system must be small (approx. 1 inch diameter). This small cell will also be capable of driving the critical fluid by heating and rotation. Since a spherically symmetric density gradient can only be made in microgravity, another small cell, of the same geometry, will be built that uses incompressible fluid. The driving of the fluid by rotation and heating in these small cells will be developed. The resulting instabilities from the driving in these two systems will then be studied. The second objective is to study the pattern forming instabilities (bifurcations) resulting from the well controlled experimental conditions in the critical fluid cell. This experiment will come close to producing conditions that are geophysically similar and will be studied as the driving parameters are changed.

Description: We are beginning laboratory experiments using magnetically active ferrofluids to study surface waves in novel geometries. Terrestrial gravity is eliminated from the dynamics, and the magnetic body force felt by ferrofluid in the presence of a magnetic field gradient is used to create a geopotential field which is a section of or an entire sphere or cylinder. New optical, electromagnetic and ultrasonic diagnostic techniques are under development to initially study capillary-gravity wave propagation and interaction in such geometries.

Description: In a recent review it is said that free-surface flows 'represent some of the difficult remaining challenges in computational fluid dynamics'. There has been progress with the development of new approaches to treating interfaces, such as the level-set method and the improvement of older methods such as the VOF method. A common theme of many of the new developments has been the regularization of discontinuities at the interface. One example of this approach is the continuum surface force (CSF) formulation for surface tension, which replaces the surface stress given by Laplace's equation by an equivalent volume force. Here, we describe how CSF formulation might be made more useful. Specifically, we consider a derivation of the CSF equations from a minimization of surface energy as outlined by Jacqmin (1996). This reformulation suggests that if one eliminates the computation of curvature in terms of a unit normal vector, parasitic currents may be eliminated. For this reformulation to work, it is necessary that transition region thickness be controlled. Various means for this, in addition to the one discussed by Jacqmin (1996), are discussed.

Description: The technique of Fluorescence Recovery After Photobleaching (FRAP) has been applied to the measurement of interfacial transport in two-phase systems. FRAP exploits the loss of fluorescence exhibited by certain fluorophores when over-stimulated (photobleached), so that a two-phase system, originally at equilibrium, can be perturbed without disturbing the interface by strong light from an argon-ion laser and its recovery monitored by a microscope-mounted CCD camera as it relaxes to a new equilibrium. During this relaxation, the concentration profiles of the probe solute are measured on both sides of the interface as a function of time, yielding information about the transport characteristics of the system. To minimize the size of the meniscus between the two phases, a photolithography technique is used to selectively treat the glass walls of the cell in which the phases are contained. This allows concentration measurements to be made very close to the interface and increases the sensitivity of the FRAP technique.

Description: Analytical and numerical studies are to be carried out to examine time-averaged thermal effects which are induced by the interaction of strong acoustic fields with a rigid boundary (thermoacoustic streaming). Also of interest is the significance of a second-order thermal expansion coefficient that emerges from this analysis. The model problem to be considered is that of a sphere that is acoustically levitated such that it is effectively isolated in a high-intensity standing acoustic field. The solution technique involves matched asymptotic analysis along with numerical solution of the boundary layer equations. The objective of this study is to predict the thermoacoustic streaming behavior and fully understand the role of the associated second-order thermodynamic modulus.

Description: The motion of a fluid-fluid-solid contact line on a rough surface is well known to display hysteresis in the contact angle vs. velocity relationship. In order to understand the phenomenon at a fundamental microscopic level, we have conducted molecular dynamics computer simulations of a Wilhelmy plate experiment in which a solid surface is dipped into a liquid bath, and the force-velocity characteristics are measured. We directly observe a systematic variation of force and contact angle with velocity, which is single-valued for the case of an atomically smooth solid surface. In the microscopically rough case, however, we find (as intuitively expected) an open hysteresis loop. Further characterization of the interface dynamics is in progress.

Description: We will pursue an improved physical understanding and mathematical model for the boundary condition at an oscillating contact line at high Reynolds number. We expect that the body force is locally unimportant for earth-based systems, and that the local behavior may dominate the mechanics of partially-filled reservoirs in the microgravity environment. One important space-based application for this contact-line study is for Faraday-waves. Oscillations in the direction of gravity (or acceleration) can dominate the fluid motion during take-off and reentry with large steady-state accelerations and in orbit, where fluctuations on the order of 10(exp -4)g occur about a zero mean. Our experience with Faraday waves has shown them to be 'cleaner' than those produced by vertical or horizontal oscillation of walls. They are easier to model analytically or computationally, and they do not have strong vortex formation at the bottom of the plate. Hence many, if not most, of the experiments will be performed in this manner. The importance of contact lines in the microgravity environment is well established. We will compare high resolution measurements of the velocity field (lO micro-m resolution) using particle-tracking and particle-image velocimetry as the fluid/fluid interface is approached from the lower fluid. The spatial gradients in the deviation provide additional means to determine an improved boundary condition and a measure of the slip region. Dissipation, the size of the eddy near the contact line, and hysteresis will be measured and compare to linear and nonlinear models of viscous and irrotational but dissipative models.

Description: The oscillation frequency and damping rate for small-amplitude axisymmetric shape modes of a gas bubble in an ideal liquid are obtained, in the limit when the bubble interface possesses Newtonian interfacial rheology with constant surface shear and dilatational viscosities. Such results permit the latter surface properties to be measured by analyzing experimental data on frequency shift and damping rate of specific shape modes of suspended bubbles in the presence of surfactants.

Description: Single bubble dynamics are investigated using acoustic techniques for isolation and manipulation. The goal of the investigations is to understand the dynamic origin of the various phenomena that bubbles exhibit: light emission, enhanced mass transport, chaotic and quasiperiodic oscillations, and translations. Once understood, acoustically manipulated bubbles can serve as platforms for materials effects on free surfaces, using surfactants to alter surface rheology and observing how that affects both dynamics and also mass transport. The effects of gravity on the problem will be shown to be significant. The first set of observations from 1g experimentation are presented. These observations are of the onset conditions for instability of the spherical shape of the bubble. For the size range 55-90 microns in diameter we observe instability governed by resonant mode coupling, which is significantly affected by the buoyant force and its effects.

Description: The proposed research is divided into three components concerned with molecular structure, molecular orientation, and continuum averages of discrete systems. In the experimental program, we propose exploring how changes in interfacial molecular structure generate contact line motion. Rather than rely on the electrostatic and electrokinetic fields arising from the molecules themselves, we augment their interactions by an imposed field at the solid/liquid interface. By controling the field, we can manipulate the molecular structure at the solid/liquid interface. In response to controlled changes in molecular structure, we observe the resultant contact line motion. In the analytical portion of the proposed research we seek to formulate a system of equations governing fluid motion which accounts for the orientation of fluid molecules. In preliminary work, we have focused on describing how molecular orientation affects the forces generated at the moving contact line. Ideally, as assumed above, the discrete behavior of molecules can be averaged into a continuum theory. In the numerical portion of the proposed research, we inquire whether the contact line region is, in fact, large enough to possess a well-defined average. Additionally, we ask what types of behavior distinguish discrete systems from continuum systems. Might the smallness of the contact line region, in itself, lead to behavior different from that in the bulk? Taken together, our proposed research seeks to identify and accurately account for some of the molecular dynamics of the moving contact line, and attempts to formulate a description from which one can compute the forces at the moving contact line.

Description: Self-diffusion coefficients that describe cross-flow migration of non-Brownian drops in a dilute sheared emulsion were obtained by trajectory calculations. A boundary integral formulation was used to describe pairwise interactions between deformable drops; interactions between undeformed drops were described with mobility functions for spherical drops. The results indicate that drops have large anisotropic self-diffusivities which depend strongly on the drop viscosity and modestly on the shear-rate. Pairwise interactions between drops in shear-flow do not appreciably promote drop breakup.