ALBERT

Description: This chapter is concerned with three-dimensional imaging of fluid flows. Although relatively young, this field of research has already yielded an enormous range of techniques. These vary widely in cost and complexity, with the cheapest light sheet systems being within the budgets of most laboratories, and the most expensive Magnetic Resonance Imaging systems available to a select few. Taking the view that the most likely systems to be developed are those using light sheets, the authors will relate their knowledge and experience of such systems. Other systems will be described briefly and references provided. Flows are inherently three-dimensional in structure; even those generated around nominally 2-D surface geometry. It is becoming increasingly apparent to scientists and engineers that the three-dimensionalities, both large and small scale, are important in terms of overall flow structure and species, momentum, and energy transport. Furthermore, we are accustomed to seeing the world in three dimensions, so it is natural that we should wish to view, measure and interpret flows in three-dimensions. Unfortunately, 3-D images do not lend themselves to convenient presentation on the printed page, and this task is one of the challenges facing us.

Description: The objectives of this effort were to determine (if possible) the best method: 1) for forcing the boundary layer to transition, 2) for assessing trip effectiveness, 3) for quantifying trip drag, 4) for testing at Reynolds numbers per foot from 5 million to maximum available rather than I to 5 million, and 5) for boundary layer state determination.

Description: NASA Langley Research Center (LaRC) has field tested a 2.0 gm, 100 Hertz, pulsed coherent lidar to detect and characterize wake vortices and to measure atmospheric winds and turbulence. The quantification of aircraft wake-vortex hazards is being addressed by the Wake Vortex Lidar (WVL) Project as part of Aircraft Vortex Spacing System (AVOSS), which is under the Reduced Spacing Operations Element of the Terminal Area Productivity (TAP) Program. These hazards currently set the minimum, fixed separation distance between two aircraft and affect the number of takeoff and landing operations on a single runway under Instrument Meteorological Conditions (IMC). The AVOSS concept seeks to safely reduce aircraft separation distances, when weather conditions permit, to increase the operational capacity of major airports. The current NASA wake-vortex research efforts focus on developing and validating wake vortex encounter models, wake decay and advection models, and wake sensing technologies. These technologies will be incorporated into an automated AVOSS that can properly select safe separation distances for different weather conditions, based on the aircraft pair and predicted/measured vortex behavior. The sensor subsystem efforts focus on developing and validating wake sensing technologies. The lidar system has been field-tested to provide real-time wake vortex trajectory and strength data to AVOSS for wake prediction verification. Wake vortices, atmospheric winds, and turbulence products have been generated from processing the lidar data collected during deployments to Norfolk (ORF), John F. Kennedy (JFK), and Dallas/Fort Worth (DFW) International Airports.

Description: Subgrid analysis of a transitional temporal mixing layer with evaporating droplets has been performed using three sets of results from a Direct Numerical Simulation (DNS) database, with Reynolds numbers (based on initial vorticity thickness) as large as 600 and with droplet mass loadings as large as 0.5. In the DNS, the gas phase is computed using a Eulerian formulation, with Lagrangian droplet tracking. The Large Eddy Simulation (LES) equations corresponding to the DNS are first derived, and key assumptions in deriving them are first confirmed by computing the terms using the DNS database. Since LES of this flow requires the computation of unfiltered gas-phase variables at droplet locations from filtered gas-phase variables at the grid points, it is proposed to model these by assuming the gas-phase variables to be the sum of the filtered variables and a correction based on the filtered standard deviation; this correction is then computed from the Subgrid Scale (SGS) standard deviation. This model predicts the unfiltered variables at droplet locations considerably better than simply interpolating the filtered variables. Three methods are investigated for modeling the SGS standard deviation: the Smagorinsky approach, the Gradient model and the Scale-Similarity formulation. When the proportionality constant inherent in the SGS models is properly calculated, the Gradient and Scale-Similarity methods give results in excellent agreement with the DNS.

Description: A three-dimensional Navier Stokes code has been used to compute the heat transfer coefficient on a film-cooled, rotating turbine blade. The blade chosen is the ACE rotor with five rows containing 93 film cooling holes covering the entire span. This is the only film-cooled rotating blade over which experimental data is available for comparison. Over 2.278 million grid points are used to compute the flow over the blade including the tip clearance region. using Wilcox's k-omega model, Coakley's q-omega model, and the zero-equation Baldwin-Lomax (B-L) model. A reasonably good comparison with the experimental data is obtained on the suction surface for all the turbulence models. At the leading edge, the B-L model yields a better comparison than tile two-equation models. On the pressure surface however the comparison between the experimental data and the prediction from the k-omega model is much better than from the other two models. Overall, the k-omega model provides the best comparison with the experimental data. However, the two-equation models require at least 40% more computational resources than the B-L model.

Description: We report measurements of the depression of the superfluid transition temperature by a heat current (1 less than or = Q less than or = 100 microW/sq cm) along the lambda-line (SVP less than or = P less than or = 21.6 bar). At P = 21.6 bar, measurements were also performed in a reduced gravity (0.2g). Experimental results show that the pressure dependence of the depression and the gravity effect on the measurements are small, in qualitative agreement with theoretical predictions. Keywords: superfluid helium; Lambda transition; heat current

Description: The differences between subcritical liquid drop and supercritical fluid drop behavior are discussed. Under subcritical, evaporative high emission rate conditions, a film layer is present in the inner part of the drop surface which contributes to the unique determination of the boundary conditions; it is this film layer which contributes to the solution's convective-diffusive character. In contrast, under supercritical condition as the boundary conditions contain a degree of arbitrariness due to the absence of a surface, and the solution has then a purely diffusive character. Results from simulations of a free fluid drop under no-gravity conditions are compared to microgravity experimental data from suspended, large drop experiments at high, low and intermediary temperatures and in a range of pressures encompassing the sub-and supercritical regime. Despite the difference between the conditions of the simulations and experiments (suspension vs. free floating), the time rate of variation of the drop diameter square is remarkably well predicted in the linear curve regime. The drop diameter is determined in the simulations from the location of the maximum density gradient, and agrees well with the data. It is also shown that the classical calculation of the Lewis number gives qualitatively erroneous results at supercritical conditions, but that an effective Lewis number previously defined gives qualitatively correct estimates of the length scales for heat and mass transfer at all pressures.

Description: We report here, preliminary data from an experiment studying flow of superfluid helium through a slit orifice (of sub-micron width) very close to T(sub lambda). Critical supercurrent (I(sub c)) data is obtained from a step function drive to the diaphragm in a Helmholtz resonator cell. The superfluid density (rho(sub s)) data can be obtained from the resonant frequency of the Helmholtz oscillator, as determined by transfer function of the resonator or from the free ringing after the step function excitation. Preliminary data shows that I(sub c) is proportional to (rho(sub s))(exp 1.27) and rho(sub s)) is proportional to tau(exp 0.73), where tau is the reduced temperature. However, the magnitude of I(sub c) is much larger than expected, indicating a possible parallel flow path. Further investigations are in progress. Keywords: superfluid; hydrodynamics; critical exponent

Description: This paper describes and discusses the textbook, Fundamentals of Computational Fluid Dynamics by Lomax, Pulliam, and Zingg, which is intended for a graduate level first course in computational fluid dynamics. This textbook emphasizes fundamental concepts in developing, analyzing, and understanding numerical methods for the partial differential equations governing the physics of fluid flow. Its underlying philosophy is that the theory of linear algebra and the attendant eigenanalysis of linear systems provides a mathematical framework to describe and unify most numerical methods in common use in the field of fluid dynamics. Two linear model equations, the linear convection and diffusion equations, are used to illustrate concepts throughout. Emphasis is on the semi-discrete approach, in which the governing partial differential equations (PDE's) are reduced to systems of ordinary differential equations (ODE's) through a discretization of the spatial derivatives. The ordinary differential equations are then reduced to ordinary difference equations (O(Delta)E's) using a time-marching method. This methodology, using the progression from PDE through ODE's to O(Delta)E's, together with the use of the eigensystems of tridiagonal matrices and the theory of O(Delta)E's, gives the book its distinctiveness and provides a sound basis for a deep understanding of fundamental concepts in computational fluid dynamics.

Description: Buoyant flows play an important role in various technological and environmental issues. For example, dispersal of pollutants, smoke, or volcano exhaust in the atmosphere, vertical motion of air, formation of clouds and other weather systems, and flows in cooling towers and fires are all determined primarily by buoyancy effects. The buoyancy force in such flows can originate from either a heat source or due to different densities between a fluid and its surroundings. Whatever the cause, the flow can be understood by studying the effects of the tight coupling between the thermal and the velocity fields since density differences can be characterized as temperature differences.

Description: Our objective is to predict droplet size distributions created by fuel injector nozzles in Jet turbines. These results will be used to determine the initial conditions for numerical simulations of the combustion process in gas turbine combustors. To predict the droplet size distribution, we are currently constructing a numerical model to understand the instability and breakup of thin conical liquid sheets. This geometry serves as a simplified model of the liquid jet emerging from a real nozzle. The physics of this process is difficult to study experimentally as the time and length scales are very short. From existing photographic data, it does seem clear that three-dimensional effects such as the formation of streamwise ligaments and the pulling back of the sheet at its edges under the action of surface tension are important.

Description: A detailed investigation of the flow physics occurring on the suction side of a simulated Low Pressure Turbine (LPT) blade was performed. A contoured upper wall was designed to simulate the y pressure distribution of an actual LPT blade onto a flat plate. The experiments were carried out at Reynolds numbers of 100,000 and 250,000 with three levels of freestream turbulence. The main emphasis in this paper is placed on flow field surveys performed at a y Reynolds number of 100,000 with levels of freestream turbulence ranging from 0.8% to 3%. Smoke-wire flow visualization data was used to confirm that the boundary layer was separated and formed a bubble. The transition process over the separated flow region is observed to be similar to a laminar free shear layer flow with the formation of a large coherent eddy structure. For each condition, the locations defining the separation bubble were determined by careful examination of pressure and mean velocity profile data. Transition onset location and length determined from intermittency profiles decrease as freestream turbulence levels increase. Additionally, the length and height of the laminar separation bubbles were observed to be inversely proportional to the levels of freestream turbulence.

Description: A distinct boundary between turbulent and non-turbulent regions in a fluid of otherwise constant properties is found in many laboratory and engineering turbulent flows, including jets, mixing layers, boundary layers and wakes. Generally, the flow has mean shear in at least one direction within t he turbulent zone, but the non-turbulent zones have no shear (adjacent laminar shear is a different case, e.g. transition in a boundary layer). There may be purely passive differences between the turbulent and non-turbulent zones, e.g. small variations in temperature or scalar concentration, for which turbulent mixing is an important issue. The boundary has several major characteristics of interest for the present study. Firstly, the boundary advances into the non-turbulent fluid, or in other words, nonturbulent fluid is entrained. Secondly, the change in turbulence properties across the boundary is remarkably abrupt; strong turbulent motions come close to the nonturbulent fluid, promoting entrainment. Thirdly, the boundary is irregular with a continually changing convoluted shape, which produces statistical intermittency. Its shape is contorted at all scales of the turbulent motion.

Description: Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This report demonstrates that a good compromise is possible with the Advected Grid Explicit (AGE) method. AGE has proven to be both efficient and accurate for simulating turbulent free-shear flows, including planar mixing layers and planar jets. Its efficiency results from its localized fully explicit finite difference formulation (Bisset 1998a,b) that is very straightforward to compute, outweighing the need for a fairly small timestep. Also, most of the successful simulations were slightly under-resolved, and therefore they were, in effect, large-eddy simulations (LES) without a sub-grid-scale (SGS) model, rather than direct numerical simulations (DNS). The principle is that the role of the smallest scales of turbulent motion (when the Reynolds number is not too low) is to dissipate turbulent energy, and therefore they do not have to be simulated when the numerical method is inherently dissipative at its resolution limits. Such simulations are termed 'auto-LES' (LES with automatic SGS modeling) in this report.

Description: A common problem occurs when refrigerant fluids wick inside the insulation of thermocouple wires through a compressor's casing feedthrough and then leak into the adjacent disconnect box outside the casing. Leaking fluids create an unfavorable situation inside the disconnect box and may contaminate the fluids. To address this problem, NASA Lewis Research Center s Manufacturing Engineering Division developed a customized hermetic feedthrough for a bank of Worthington compressors. In these compressors, bearing temperatures are measured by internal thermocouples embedded in bearings located inside the compressor casings. The thermocouple wires need to be routed outside the casing and read at another location. These wires are short and are terminated to a disconnect strip inside the casing. The bearings operate at about 170 F, but because the casing is filled with R12 refrigerant oil, the casing has a maximum temperature of about 100 F. The operating conditions of these compressors permit the use of an epoxy that is compatible with the R12 fluid. The desired finished product is a stainless steel tube that has been filled solid with epoxy after thermocouple wires bonded and sealed by epoxy have been inserted through its length. Shrink tubing extends from both ends of the tube. The process that was developed to isolate the thermocouple wires from the R12 fluid follows. For this application, use an 8-in.-long piece of 0.500-in. 304 stainless steel tube with six pairs of 24-gauge stranded, PTFE-insulated (polytetrafluoroethylene) type "T" thermocouple wires for each feedthrough. Use shrink tubing to strain relief the insulated wires at their exit from the stainless steel tube. Cut the wire to length and identify the location of the stainless steel tube sleeve with masking tape. Then, remove the outer insulation from a 2-in. section of wire that will be inside the tube, and carefully strip to bare wire a 1-in. section in the middle of the section with the outer insulation removed. For an effective seal, the epoxy must penetrate between the strands when stranded conductors are used. Make the seal with epoxy bond on the bare wire. The bare wire must be encapsulated with a thin layer of the epoxy that leaves only a very low profile. These encapsulated wires must cure before the assembly can be continued. Then, inspect the cured wires for complete encapsulation before going to the next step. Insert the wires in the stainless steel tube and orient them so that the epoxied stripped sections are staggered within the tube; then, apply shrink tubing to one end of the cleaned wires, positioning it inside the edge of the tube. The small gaps between the wires on the other end will be used to inject the epoxy into the tube. Let the epoxy cure inside the tube, free of any voids. Then, continue to fill the tube until the entire 8-in. length is nearly filled, allowing room for the other strain-relieving shrink tubing. Since this first design, the process has been adjusted to fit many needs and situations. Customized feedthroughs have been assembled from various wire types, wire gauges, and/or stainless steel tube passages. The fittings selected to mount these feedthroughs allow their use in other areas, such as pressure or vacuum systems.

Description: High-Speed Civil Transport (HSCT) engine combustor liners are subjected to complex thermal environments and have to endure these for thousands of hours with assured reliability. In the past, several deterministic analyses have been performed, including detailed heat transfer analyses to obtain thermal profiles and deterministic stress analyses to identify critical locations of high stresses. Actual rig tests also have been performed for segments by simulating these loading situations as closely as possible. However, it is well known that many uncertainties exist in loading (primarily thermal loads due to heat transfer), boundary conditions (end fixity unknowns), and material properties (moduli, thermal-expansion coefficients, and conductivities). The present in-house effort at the NASA Lewis Research Center is directed toward accounting for these in a formal way to assess the performance of liner components under complex and uncertain loading conditions as well as subject to other geometry- and material-related uncertainties.

Description: This paper presents the results of a numerical study of shock/shock interactions that include both the Edney type IV and type III interactions, with emphasis on the type IV interactions. Computations are made using the direct simulation Monte Carlo (DSMC) method of Bird for Mach 10 air flow, as produced in the ONERA R5Ch low-density wind tunnel. The simulations include the flow about a shock generator which creates a relatively weak oblique shock that impinges on a much stronger cylinder bow shock. The sensitivity and characteristics of the interactions are examined by varying the horizontal distance separating the shock generator leading edge and cylinder. Results of the simulation for one separation distance are compared with wind tunnel measurements. Comparisons are made for surface heating and pressure and for flow-field values of density and rotational temperatures, as obtained with the Dual-line Coherent Anti-Stokes Scattering (DL-CARS) technique. The comparisons between experiment and calculation yield a consistent description of the shock interaction features and a consistent description of the surface heating and pressure distributions, with the exception of the peak values-the computed values being greater than the measured values.

Description: Problems concerning laminar-turbulent transition are addressed by solving a series of initial value problems. Solutions to the temporal, initial-value problem .with an inhomogeneous forcing term imposed upon the flow are sought. It is shown that: (1) A transient disturbance lying located outside of the boundary layer can lead to the growth of an unstable Tollmein-Schlicting wave; (2) A resonance with the continuous spectrum may provide a mechanism for bypass transition; and (3) The continuum modes of a disturbance feed directly into the Tollmein-Schlicting wave downstream through non-parallel effects.

Description: At the NASA Lewis Research Center, a comprehensive assessment was made of the predictive capability of the average passage flow model as applied to multistage axial-flow compressors. This model, which describes the time-averaged flow field within a typical passage of a blade row embedded in a multistage configuration, is being widely used throughout U.S. aircraft industry as an integral part of their design systems. Rotor flow-angle deviation. In this work, detailed data taken within a four and one-half stage large low-speed compressor were used to assess the weaknesses and strengths of the predictive capabilities of the average passage flow model. The low-speed compressor blading is of modern design and employs stator end-bends. Measurements were made with slow- and high response instrumentation. The high-response measurements revealed the velocity components of both the rotor and stator wakes. From the measured wake profiles, we found that the flow exiting the rotors deviated from the rotor exit metal angle to a lesser degree than was predicted by the average passage flow model. This was found to be due to blade boundary layer transition, which recently has been shown to exist on multistage axial compressor rotor and stator blades, but was not accounted for in the average passage model. Consequently, a model that mimics the effects of blade boundary layer transition, Shih k-epsilon model, was incorporated into the average passage model. Simulations that incorporated this transition model showed a dramatic improvement in agreement with data. The altered model thus improved predictive capability for multistage axial-flow compressors, and this was verified by detailed experimental measurement.

Description: A thorough understanding of dynamic interactions between inlets and compressors is extremely important to the design and development of propulsion control systems, particularly for supersonic aircraft such as the High-Speed Civil Transport (HSCT). Computational fluid dynamics (CFD) codes are routinely used to analyze individual propulsion components. By coupling the appropriate CFD component codes, it is possible to investigate inlet-compressor interactions. The objectives of this work were to gain a better understanding of inlet-compressor interaction physics, formulate a more realistic compressor-face boundary condition for time-accurate CFD simulations of inlets, and to take a first step toward the CFD simulation of an entire engine by coupling multidimensional component codes. This work was conducted at the NASA Lewis Research Center by a team of civil servants and support service contractors as part of the High Performance Computing and Communications Program (HPCCP).

Description: NASA Lewis Research Center's Advanced Controls and Dynamics Technology Branch is investigating active control strategies to mitigate or eliminate the combustion instabilities prevalent in lean-burning, low-emission combustors. These instabilities result from coupling between the heat-release mechanisms of the burning process and the acoustic flow field of the combustor. Control design and implementation require a simulation capability that is both fast and accurate. It must capture the essential physics of the system, yet be as simple as possible. A quasi-one-dimensional, computational fluid dynamics (CFD) based simulation has been developed which may meet these requirements. The Euler equations of mass, momentum, and energy have been used, along with a single reactive species transport equation to simulate coupled thermoacoustic oscillations. A very simple numerical integration scheme was chosen to reduce computing time. Robust boundary condition procedures were incorporated to simulate various flow conditions (e.g., valves, open ends, and choked inflow) as well as to accommodate flow reversals that may arise during large flow-field oscillations. The accompanying figure shows a sample simulation result. A combustor with an open inlet, a choked outlet, and a large constriction approximately two thirds of the way down the length is shown. The middle plot shows normalized, time-averaged distributions of the relevant flow quantities, and the bottom plot illustrates the acoustic mode shape of the resulting thermoacoustic oscillation. For this simulation, the limit cycle peak-to-peak pressure fluctuations were 13 percent of the mean. The simulation used 100 numerical cells. The total normalized simulation time was 50 units (approximately 15 oscillations), which took 26 sec on a Sun Ultra2.

Description: Pulse Phase Thermography (PPT) has been reported as a novel powerful technique of the thermal NDE. It employs application of the Discrete Fourier Transform (DFT) to thermal images obtained following flash heating of the front surface of a specimen. The computed phasegrams are excellent for defect visualization in a wide range of materials. This is in part due to their low sensitivity to uneven heating. This work is an attempt to analyze advantages and limitations of PPT. Results of application of the DFT to simulated temperature decays are presented. The temperature evolution on a surface has been simulated based on an analytical solution of the one-dimensional heat diffusion problem. A more sophisticated study has been done for different sizes of defects using numerical solution of the three-dimensional mathematical model. Capabilities of PPT for in-depth scanning and for monitoring of the material loss are discussed. The recommendations for the practical application of the PPT are presented. Experimental results obtained following these recommendations are reported.

Description: Inertial and gravitational forces determine a particle's motion in a turbulent flow field. Gravity plays the dominant role in this motion by pulling the particles through adjacent regions of fluid turbulence. To better understand and model how a particle's inertia effects its displacement, one must examine the dispersion in a turbulent flow in the absence of gravity. In this paper, we present the particle experiments planned for NASA's KC-135 Reduced-Gravity Aircraft, which generates microgravity conditions for about 20 seconds. We also predict the particle behavior using simulation and ground-based experiments. We will release particles with Stokes numbers of 0.1, 1, and 10 into an enclosed tank of near-isotropic, stationary, and homogenous turbulence. These particle Stoke numbers cover a broad range of flow regimes of interest. Two opposed grids oscillating back and forth generate the turbulent field in the tank with a range of turbulence scales that covers about three orders of magnitude and with turbulence intensities of about ten times the mean velocity. The motion of the particles will be tracked using a stereo image velocimetry technique.

Description: Experimental results are presented for drops and bubbles levitated in a liquid host, with particular attention given to the effect of shape oscillations and capillary waves on the local flow fields. Some preliminary results are also presented on the use of streaming flows for the control of evaporation rate and rotation of electrostatically levitated droplets in 1 g. The results demonstrate the potential for the technological application of acoustic methods to active control of forced convection in microgravity.

Description: We have developed a new instrument that can measure fast transient birefringence and polymer chain orientation angle in complex fluids. The instrument uses a dual-crystal transverse electro-optic modulator with the second crystal's modulation voltage applied 180 deg out of phase from that of the first crystal. In this manner, the second crystal compensates for the intrinsic static birefringence of the first crystal, and it doubles the modulation depth. By incorporating a transverse electro-optic modulator with two lithium-niobate (LiNbO3) crystals oriented orthogonal to each other with a custom-designed optical system, we have produced a very small robust instrument capable of fast transient retardation measurements. By measuring the sample thickness or optical path length through the sample, we can calculate the transient birefringence. This system can also measure dichroism. We have compared the calibration results and retardation and orientation angle measurements of this instrument with those of a photoelastic modulator (PEM) based system using a quarter wave plate and a high-precision 1/16-wave plate to simulate a birefringent sample. Transient birefringence measurements on the order of 10(exp -9) can be measured using either modulator.

Description: The bulk behavior of dispersed, fluidized, or undispersed stationary granular systems cannot be fully understood in terms of adhesive/cohesive properties without understanding the role of electrostatic forces acting at the level of the grains themselves. When grains adhere to a surface, or come in contact with one another in a stationary bulk mass, it is difficult to measure the forces acting on the grains, and the forces themselves that induced the cohesion and adhesion are changed. Even if a single grain were to be scrutinized in the laboratory, it might be difficult, perhaps impossible, to define the distribution and character of surface charging and the three-dimensional relationship that charges (electrons, holes) have to one another. The hypothesis that we propose to test in microgravity (for dielectric materials) is that adhesion and cohesion of granular matter are mediated primarily by dipole forces that do not require the presence of a net charge; in fact, nominally electrically neutral materials should express adhesive and cohesive behavior when the neutrality results from a balance of positive and negative charge carriers. Moreover, the use of net charge alone as a measure of the electrical nature of grain-to-grain relationships within a granular mass may be misleading. We believe that the dipole forces arise from the presence of randomly-distributed positive and negative fixed charge carriers on grains that give rise to a resultant dipole moment. These dipole forces have long-range attraction. Random charges are created whenever there is triboelectrical activity of a granular mass, that is, whenever the grains experience contact/separation sequences or friction.

Description: The present work is aimed at the experimental studies and numerical modeling of the bubble growth mechanisms of a single bubble attached to a heating surface and of a bubble sliding along an inclined heated plate. Single artificial cavity of 10 microns in diameter was made on the polished Silicon wafer which was electrically heated at the back side in order to control the surface nucleation superheat. Experiments with a sliding bubble were conducted at different inclination angles of the downward facing heated surface for the purpose of studying the effect of magnitude of components of gravity acting parallel to and normal to the heat transfer surface. Information on the bubble shape and size, the bubble induced liquid velocities as well as the surface temperature were obtained using the high speed imaging and hydrogen bubble techniques. Analytical/numerical models were developed to describe the heat transfer through the micro-macro layer underneath and around a bubble formed at a nucleation site. In the micro layer model the capillary and disjoining pressures were included. Evolution of the bubble-liquid interface along with induced liquid motion was modeled. As a follow-up to the studies at normal gravity, experiments are being conducted in the KC-135 aircraft to understand the bubble growth/detachment under low gravity conditions. Experiments have been defined to be performed under long duration of microgravity conditions in the space shuttle. The experiment in the space shuttle will provide bubble growth and detachment data at microgravity and will lead to validation of the nucleate boiling heat transfer model developed from the preceding studies conducted at normal and low gravity (KC-135) conditions.

Description: In this work, we study steady and oscillatory thermocapillary and natural convective flows generated by a bubble on a heated solid surface. The interaction between gas and vapor bubbles with the surrounding fluid is of interest for both space and ground-based processing. A combined numerical-experimental approach is adopted here. The temperature field is visualized using Mach-Zehnder and/or Wollaston Prism Interferometry and the flow field is observed by a laser sheet flow visualization technique. A finite element numerical model is developed which solves the transient two-dimensional continuity, momentum, and energy equations and includes the effects of temperature-dependent surface tension and bubble surface deformation. Below the critical Marangoni number, the steady state low-g and 1-g temperature and velocity fields predicted by the finite element model are in excellent agreement with both the visualization experiments in our laboratory and recently published experimental results in the literature. Above the critical Marangoni number, the model predicts an oscillatory flow which is also closely confirmed by experiments. It is shown that the dynamics of the oscillatory flow are directly controlled by the thermal and hydrodynamic interactions brought about by combined natural and thermocapillary convection. Therefore, as numerical simulations show, there are considerable differences between the 1-g and low-g temperature and flow fields at both low and high Marangoni numbers. This has serious implications for both materials processing and fluid management in space.

Description: The WIND code is a general-purpose, structured, multizone, compressible flow solver that can be used to analyze steady or unsteady flow for a wide range of geometric configurations and over a wide range of flow conditions. WIND is the latest product of the NPARC Alliance, a formal partnership between the NASA Lewis Research Center and the Air Force Arnold Engineering Development Center (AEDC). WIND Version 1.0 was released in February 1998, and Version 2.0 will be released in February 1999. The WIND code represents a merger of the capabilities of three existing computational fluid dynamics codes--NPARC (the original NPARC Alliance flow solver), NXAIR (an Air Force code used primarily for unsteady store separation problems), and NASTD (the primary flow solver at McDonnell Douglas, now part of Boeing).

Description: About half of all scientific and commercial spacecraft spin during some part of their mission. Although spinning has many benefits (increasing stability, controlling the location of liquid propellants, and distributing solar heat loads), it also creates problems because a precession (or wobble) motion is unavoidable. For modern spacecraft, by far the largest source of wobble is energy dissipation in the liquid of partially filled tanks. The liquid s energy dissipation cannot, however, be quantified adequately by any ground test. Current analytical models are also inadequate because fundamental data on fluid motion in low gravity are needed to validate them. Consequently, spacecraft attitude-control systems are designed and operated very conservatively. Nonetheless, spacecraft often still perform poorly in orbit, and some have been lost because of a rapid increase of the wobble rate. The Liquid Motion Experiment (LME) was designed to provide spacecraft designers accurate data on the wobble dynamics of spacecraft that contain large quantities of mobile liquids. LME, which was flown on the space shuttle mission STS-84, was built under contract to the NASA Lewis Research Center by the Southwest Research Institute of San Antonio, Texas. Major accomplishments for 1998 include reduction of the flight data and publication of the experimental results. LME was essentially a spin table that created a realistic nutation motion of scale-model tanks containing liquid. Two spherical and two cylindrical transparent tanks were tested simultaneously, and three sets of such tanks were employed to vary liquid viscosity, fill level, and propellant management device (PMD) design. All the tanks were approximately 4.5 in. in diameter. The primary test measurements were the radial and tangential torques exerted on the tanks by the liquid. These torques could not be measured on the ground because of the masking effects of gravity.

Description: For the last several years, LeRC-HT, a three-dimensional computational fluid dynamics (CFD) computer code for analyzing gas turbine flow and convective heat transfer, has been evolving at the NASA Lewis Research Center. The code is unique in its ability to give a highly detailed representation of the flow field very close to solid surfaces. This is necessary for an accurate representation of fluid heat transfer and viscous shear stresses. The code has been used extensively for both internal cooling passage flows and hot gas path flows--including detailed film cooling calculations, complex tip-clearance gap flows, and heat transfer. In its current form, this code has a multiblock grid capability and has been validated for a number of turbine configurations. The code has been developed and used primarily as a research tool (at least 35 technical papers have been published relative to the code and its application), but it should be useful for detailed design analysis. We now plan to make this code available to selected users for further evaluation.

Description: The overall performance of a centrifugal compressor depends on the performance of the impeller and diffuser as well as on the interactions occurring between these components. Accurate measurements of the flow fields in each component are needed to develop computational models that can be used in compressor design codes. These measurements must be made simultaneously over an area that covers both components so that researchers can understand the interactions occurring between the two components. Optical measurement techniques are being used at the NASA Lewis Research Center to measure the velocity fields present in both the impeller and diffuser of a 4:1 pressure ratio centrifugal compressor operating at several conditions ranging from design flow to surge. Laser Doppler Velocimetry (LDV) was used to measure the intrablade flows present in the impeller, and the results were compared with analyses obtained from two three-dimensional viscous codes. The development of a region of low throughflow velocity fluid within this high-speed impeller was examined and compared with a similar region first observed in a large low-speed centrifugal impeller at Lewis. Particle Image Velocimetry (PIV) is a relatively new technique that has been applied to measuring the diffuser flow fields. PIV can collect data rapidly in the diffuser while avoiding the light-reflection problems that are often encountered when LDV is used. The Particle Image Velocimeter employs a sheet of pulsed laser light that is introduced into the diffuser in a quasi-radial direction through an optical probe inserted near the diffuser discharge. The light sheet is positioned such that its centerline is parallel to the hub and shroud surfaces and such that it is parallel to the diffuser vane, thereby avoiding reflections from the solid surfaces. Seed particles small enough to follow the diffuser flow are introduced into the compressor at an upstream location. A high-speed charge-coupled discharge (CCD) camera is synchronized to the laser pulse rate; this allows it to capture images of seed particle position that are separated by a small increment in time. A crosscorrelation of a particle's position in two consecutive images provides an estimate of flow velocity and direction. Multiple image pairs obtained in rapid succession at a particular flow condition provide enough measurements for statistical significance. PIV provides simultaneous velocity measurements over the entire plane that is illuminated by the light sheet instead of at a single point, as is the case when LDV is used. PIV has a further advantage in that the laser light pulse can be triggered by an external source such as a high-response pressure transducer. This feature will allow PIV to synchronize flow imaging to physical phenomena such as rotating stall or stall precursor waves. We hope that this technique can be used to obtain images of the flow field during and just prior to stall.

Description: Recently, refractive secondary solar concentrator systems were developed for solar thermal power and propulsion (ref. 1). Single-crystal oxides-such as yttria-stabilized zirconia (Y2O3-ZrO2), yttrium aluminum garnet (Y3Al5O12, or YAG), magnesium oxide (MgO), and sapphire (Al2O3)-are candidate refractive secondary concentrator materials. However, the refractive concentrator system will experience high-temperature thermal cycling in the solar thermal engine during the sun/shade transition of a space mission. The thermal mechanical reliability of these components in severe thermal environments is of great concern. Simulated mission tests are important for evaluating these candidate oxide materials under a variety of transient and steady-state heat flux conditions. In this research at the NASA Lewis Research Center, a controlled heat flux test approach was developed for investigating the thermal mechanical stability of the candidate oxide. This approach used a 3.0-kW continuous-wave (wavelength, 10.6 mm) carbon dioxide (CO2) laser (ref. 2). The CO2 laser is especially well-suited for single-crystal thermal shock tests because it can directly deliver well-characterized heat energy to the oxide surfaces. Since the oxides are opaque at the 10.6-mm wavelength of the laser beam, the light energy is absorbed at the surfaces rather than transmitting into the crystals, and thus generates the required temperature gradients within the specimens. The following figure is a schematic diagram of the test rig.

Description: The number of optical techniques that may potentially be used during a given wind tunnel test is continually growing. These include parameter sensitive paints that are sensitive to temperature or pressure, several different types of off-body and on-body flow visualization techniques, optical angle-of-attack (AoA), optical measurement of model deformation, optical techniques for determining density or velocity, and spectroscopic techniques for determining various flow field parameters. Often in the past the various optical techniques were developed independently of each other, with little or no consideration for other techniques that might also be used during a given test. Recently two optical techniques have been increasingly requested for production measurements in NASA wind tunnels. These are the video photogrammetric (or videogrammetric) technique for measuring model deformation known as the video model deformation (VMD) technique, and the parameter sensitive paints for making global pressure and temperature measurements. Considerations for, and initial attempts at, simultaneous measurements with the pressure sensitive paint (PSP) and the videogrammetric techniques have been implemented. Temperature sensitive paint (TSP) has been found to be useful for boundary-layer transition detection since turbulent boundary layers convect heat at higher rates than laminar boundary layers of comparable thickness. Transition is marked by a characteristic surface temperature change wherever there is a difference between model and flow temperatures. Recently, additional capabilities have been implemented in the target-tracking videogrammetric measurement system. These capabilities have permitted practical simultaneous measurements using parameter sensitive paint and video model deformation measurements that led to the first successful unified test with TSP for transition detection in a large production wind tunnel.

Description: Microscopic spheres suspended in liquid become highly ordered under the proper conditions. Such collections of particles, called colloidal suspensions or colloids, are the subject of a series of ongoing microgravity experiments at the NASA Lewis Research Center. By studying the way these colloidal suspensions order themselves, scientists can better understand how atoms of a liquid become ordered to form a solid. In addition, highly ordered colloids have special properties that may make them useful in future hightech applications. Work is underway at Lewis to develop an optical microscope to view these colloidal suspensions sphere by sphere in microgravity.

Description: Under the Advanced Subsonic Technology (AST) Program, an aeroelastic analysis code (TURBO-AE) based on Navier-Stokes equations is currently under development at NASA Lewis Research Center s Machine Dynamics Branch. For a blade row, aeroelastic instability can occur in any of the possible interblade phase angles (IBPA s). Analyzing small IBPA s is very computationally expensive because a large number of blade passages must be simulated. To reduce the computational cost of these analyses, we used time shifted, or phase-lagged, boundary conditions in the TURBO-AE code. These conditions can be used to reduce the computational domain to a single blade passage by requiring the boundary conditions across the passage to be lagged depending on the IBPA being analyzed. The time-shifted boundary conditions currently implemented are based on the direct-store method. This method requires large amounts of data to be stored over a period of the oscillation cycle. On CRAY computers this is not a major problem because solid-state devices can be used for fast input and output to read and write the data onto a disk instead of storing it in core memory.

Description: A complete "geometry to drag-polar" analysis capability for three-dimensional high-lift configurations is described. The approach is based on the use of unstructured meshes in order to enable rapid turnaround for complicated geometries which arise in high-lift con gurations. Special attention is devoted to creating a capability for enabling analyses on highly resolved grids. Unstructured meshes of several million vertices are initially generated on a work-station, and subsequently refined on a supercomputer. The flow is solved on these refined meshes on large parallel computers using an unstructured agglomeration multigrid algorithm. Good prediction of lift and drag throughout the range of incidences is demonstrated on a transport take-off configuration using up to 24.7 million grid points. The feasibility of using this approach in a production environment on existing parallel machines is demonstrated, as well as the scalability of the solver on machines using up to 1450 processors.

Description: Development of HPF versions of NPB and ARC3D showed that HPF has potential to be a high level language for parallelization of CFD applications. The use of HPF requires an intimate knowledge of the applications and a detailed analysis of data affinity, data movement and data granularity. Since HPF hides data movement from the user even with this knowledge it is easy to overlook pieces of the code causing low performance of the application. In order to simplify and accelerate the task of developing HPF versions of existing CFD applications we have designed and partially implemented ADAPT (Automatic Data Distribution and Placement Tool). The ADAPT analyzes a CFD application working on a single structured grid and generates HPF TEMPLATE, (RE)DISTRIBUTION, ALIGNMENT and INDEPENDENT directives. The directives can be generated on the nest level, subroutine level, application level or inter application level. ADAPT is designed to annotate existing CFD FORTRAN application performing computations on single or multiple grids. On each grid the application can considered as a sequence of operators each applied to a set of variables defined in a particular grid domain. The operators can be classified as implicit, having data dependences, and explicit, without data dependences. In order to parallelize an explicit operator it is sufficient to create a template for the domain of the operator, align arrays used in the operator with the template, distribute the template, and declare the loops over the distributed dimensions as INDEPENDENT. In order to parallelize an implicit operator, the distribution of the operator's domain should be consistent with the operator's dependences. Any dependence between sections distributed on different processors would preclude parallelization if compiler does not have an ability to pipeline computations. If a data distribution is "orthogonal" to the dependences of an implicit operator then the loop which implements the operator can be declared as INDEPENDENT.

Description: In this work the gas volume fraction and the root-mean-squared fluid velocity are measured in buoyancy driven shear flows of bubble suspensions in a tall, inclined, rectangular channel. The experiments are performed under conditions where We 〈〈 1a nd Re 〉〉 1, for which comparisons are made with kinetic theory and numerical simulations. Here Re = gamma(a(exp 2)/nu is the Reynolds number and We = rho(gamma(exp 2))a(exp 3)/sigma is the Weber number; gamma is the shear rate, a is the bubble radius, nu is the kinematic viscosity of the liquid, rho is the density of the liquid, and sigma is the surface tension of the gas/liquid interface. Kang et al. calculated the bubble phase pressure and velocity variance of sheared bubble suspensions under conditions where the bubbles are spherical and the liquid phase velocity field can be approximated using potential flow theory, i.e. We= 0 and Re 〉〉 1. Such conditions can be achieved in an experiment using gas bubbles, with a radius of O(0.5mm), in water. The theory requires that there be no average relative motion of the gas and liquid phases, hence the motivation for an experimental program in microgravity. The necessity of performing preliminary, Earth based experiments, however, requires performing experiments where the gas phase rises in the liquid, which significantly complicates the comparison of experiments with theory. Rather than comparing experimental results with theory for a uniform, homogeneous shear flow, experiments can be compared directly with solutions of the averaged equations of motion for bubble suspensions. This requires accounting for the significant lift force acting on the gas phase when the bubbles rise parallel to the average velocity of the sheared suspension. Shear flows can be produced in which the bubble phase pressure gradient, arising from shear induced collisions amongst the bubbles, balances a body force (centrifugal or gravitational) on the gas phase. A steady, non-uniform gas volume fraction can be measured, from which the bubble phase pressure gradient can be obtained and compared to theory and numerical simulations. The presence of bounding walls further complicates the experiments, since the detailed interactions of the bubbles with bounding walls is not well understood, especially in the presence of gravity, where the momentum and energy exchange depends on the inclination of the wall.

Description: A two-phase test loop has been designed and constructed to generate the necessary data for two-phase pressure drop and Critical Heat Flux (CHF) under reduced gravity conditions. A series of airplane trajectory tests aboard NASA KC-135 were performed and the data was used to evaluate the applicability of the earth gravity models for prediction of the reduced gravity data. Several commonly used correlations for the two-phase friction multiplier and critical heat flux were used to predict the data. It was generally concluded that the two-phase pressure drop can be predicted by the earth gravity correlations. The critical heat flux under reduced gravity conditions did not show a strong dependence on mass flow rate and the measured CHF were generally lower than the equivalent 1g conditions. The earth gravity models need to be modified for application to reduced gravities.

Description: In a wind tunnel facility, the direct measurement of forces and moments induced on the model are performed by a force measurement balance. The measurement balance is a precision-machined device that has strain gages at strategic locations to measure the strain (i.e., deformations) due to applied forces and moments. The strain gages convert the strain (and hence the applied force) to an electrical voltage that is measured by external meters. Thermal gradients can complicate the process, however. Thermal gradients on the balance cause differential expansion (or contraction) of various parts of the balance that induce a strain that is detected by the strain gages and is indistinguishable from an external applied force. The thermal gradients can result when testing is done at elevated temperatures or at cryogenic temperatures such as at the National Transonic Facility (NTF) at NASA Langley Research Center (LaRC).

Description: Boundary layer and aeroheating characteristics of several X-33 configurations have been experimentally examined in the Langley 20-Inch Mach 6 Air Tunnel. Global surface heat transfer distributions, surface streamline patterns, and shock shapes were measured on 0.013-scale models at Mach 6 in air. Parametric variations include angles-of-attack of 20-deg, 30-deg, and 40-deg; Reynolds numbers based on model length of 0.9 to 6.6 million; and body-flap deflections of 0, 10 and 20-deg. The effects of discrete and distributed roughness elements on boundary layer transition, which included trip height, size, location, and distribution, both on and off the windward centerline, were investigated. The discrete roughness results on centerline were used to provide a transition correlation for the X-33 flight vehicle that was applicable across the range of reentry angles of attack. The attachment line discrete roughness results were shown to be consistent with the centerline results, as no increased sensitivity to roughness along the attachment line was identified. The effect of bowed panels was qualitatively shown to be less effective than the discrete trips; however, the distributed nature of the bowed panels affected a larger percent of the aft-body windward surface than a single discrete trip.

Description: A parametric weight assessment of advanced metallic panel, ceramic blanket, and ceramic tile thermal protection systems (TPS) was conducted using an implicit, one-dimensional (1 -D) thermal finite element sizing code. This sizing code contained models to ac- count for coatings, fasteners, adhesives, and strain isolation pads. Atmospheric entry heating profiles for two vehicles, the Access to Space (ATS) rocket-powered single-stage-to-orbit (SSTO) vehicle and a proposed Reusable Launch Vehicle (RLV), were used to ensure that the trends were not unique to a particular trajectory. Eight TPS concepts were compared for a range of applied heat loads and substructural heat capacities to identify general trends. This study found the blanket TPS concepts have the lightest weights over the majority of their applicable ranges, and current technology ceramic tiles and metallic TPS concepts have similar weights. A proposed, state-of-the-art metallic system which uses a higher temperature alloy and efficient multilayer insulation was predicted to be significantly lighter than the ceramic tile systems and approaches blanket TPS weights for higher integrated heat loads.

Description: A new method has been developed to accelerate the convergence of explicit time-marching, laminar, Navier-Stokes codes through the combination of local preconditioning and multi-stage time marching optimization. Local preconditioning is a technique to modify the time-dependent equations so that all information moves or decays at nearly the same rate, thus relieving the stiffness for a system of equations. Multi-stage time marching can be optimized by modifying its coefficients to account for the presence of viscous terms, allowing larger time steps. We show it is possible to optimize the time marching scheme for a wide range of cell Reynolds numbers for the scalar advection-diffusion equation, and local preconditioning allows this optimization to be applied to the Navier-Stokes equations. Convergence acceleration of the new method is demonstrated through numerical experiments with circular advection and laminar boundary-layer flow over a flat plate.

Description: This paper discusses a method for the identification and application of reduced-order models based on linear and nonlinear aerodynamic impulse responses. The Volterra theory of nonlinear systems and an appropriate kernel identification technique are described. Insight into the nature of kernels is provided by applying the method to the nonlinear Riccati equation in a non-aerodynamic application. The method is then applied to a nonlinear aerodynamic model of RAE 2822 supercritical airfoil undergoing plunge motions using the CFL3D Navier-Stokes flow solver with the Spalart-Allmaras turbulence model. Results demonstrate the computational efficiency of the technique.

Description: USM3D is a widely-used unstructured flow solver for simulating inviscid and viscous flows over complex geometries. The current version (version 5.0) of USM3D, however, does not have advanced turbulence models to accurately simulate complicated flows. We have implemented two modified versions of the original Jones and Launder k-epsilon two-equation turbulence model and the Girimaji algebraic Reynolds stress model in USM3D. Tests have been conducted for two flat plate boundary layer cases, a RAE2822 airfoil and an ONERA M6 wing. The results are compared with those of empirical formulae, theoretical results and the existing Spalart-Allmaras one-equation model.

Description: The spreading characteristics of jets from several asymmetric nozzles, and a set of rectangular orifices are compared, covering a jet Mach number range of 0.3-2.0. The effect of 'tabs' for a rectangular and a round nozzle is also included in the comparison. Compared to a round jet, the jets from the asymmetric nozzles spread only slightly more at subsonic conditions whereas at supersonic conditions, when 'screech' occurs, they spread much more. The dynamics of the azimuthal vortical structures of the jet, organized and intensified under the screeching condition, are thought to be responsible for the observed effect at supersonic conditions. Curiously, the jet from a 'lobed' nozzle spreads much less at supersonic condition compared to all other cases; this is due to the absence of screech with this nozzle. Screech stages inducing flapping, rather than varicose or helical, flow oscillation cause a more pronounced jet spreading. At subsonic conditions, only a slight increase in jet spreading with the asymmetric nozzles contrasts previous observations by others. The present results show that the spreading of most asymmetric jets is not much different from that of a round jet. This inference is further supported by data from the rectangular orifices. In fact, jets from the orifices with small aspect ratio (AR) exhibit virtually no increase in the spreading. A noticeable increase commences only when AR is larger than about 10. Thus, 'shear layer perimeter stretching', achieved with a larger AR for a given cross-sectional area of the orifice, by itself, proves to be a relatively inefficient mechanism for increasing jet spreading. In contrast, the presence of streamwise vortices or 'natural excitation' can cause a significant increase - effects that might explain the observations in the previous investigations. Thus far, the biggest increase in jet spreading is observed with the tabs. This is true in the subsonic regime, as well as in the supersonic regime, in spite of the fact that screech is eliminated by the tabs. The characteristic spreading of the tabbed jets is explained by the induced motion of the tab-generated streamwise vortex pairs. The tabs, however, incur thrust loss; the flow blockage and loss in thrust coefficient, vis-a-vis the spreading increase, are evaluated for various configurations.

Description: The Ninth Thermal and Fluids Analysis Workshop (TFAWS 98) was held at the Ohio Aerospace Institute in Cleveland, Ohio from August 31 to September 4, 1998. The theme for the hands-on training workshop and conference was "Integrating Computational Fluid Dynamics and Heat Transfer into the Design Process." Highlights of the workshop (in addition to the papers published herein) included an address by the NASA Chief Engineer, Dr. Daniel Mulville; a CFD short course by Dr. John D. Anderson of the University of Maryland; and a short course by Dr. Robert Cochran of Sandia National Laboratories. In addition, lectures and hands-on training were offered in the use of several cutting-edge engineering design and analysis-oriented CFD and Heat Transfer tools. The workshop resulted in international participation of over 125 persons representing aerospace and automotive industries, academia, software providers, government agencies, and private corporations. The papers published herein address issues and solutions related to the integration of computational fluid dynamics and heat transfer into the engineering design process. Although the primary focus is aerospace, the topics and ideas presented are applicable to many other areas where these and other disciplines are interdependent.

Description: This paper reviews the effects of 'screech', 'asymmetric nozzle shaping', 'tabs' and 'overexpansion' on the spreading of free jets. Corresponding thrust penalty for the tabs and overexpanded condition are also evaluated. The asymmetric shapes include rectangular ones with varying aspect ratio. Tabs investigated are triangular shaped 'delta-tabs' placed at the exit of a convergent circular nozzle. The effect of overexpansion is examined with circular convergent-divergent (C-D) nozzles. Tabs and overexpansion are found to yield the largest increase in jet spreading. Each, however, involves a performance penalty, i.e., a loss in thrust coefficient. Variation of the size of four delta-tabs show that there exists an optimum size for which the gain in jet spreading is the maximum per unit loss in thrust coefficient. With the C-D nozzles, the minimum in thrust coefficient is expected near the beginning of the overexpanded regime based on idealized flow calculations. The maximum increase in jet spreading, however, is found to occur at higher pressure ratios well into the overexpanded regime. The optimum benefit with the overexpanded flow, in terms of gain in spreading for unit penalty, is found to be comparable to the optimum tab case.

Description: Continuous improvement of aerospace product development processes is a driving requirement across much of the aerospace community. As up to 90% of the cost of an aerospace product is committed during the first 10% of the development cycle, there is a strong emphasis on capturing, creating, and communicating better information (both requirements and performance) early in the product development process. The community has responded by pursuing the development of computer-based systems designed to enhance the decision-making capabilities of product development individuals and teams. Recently, the historical foci on sharing the geometrical representation and on configuration management are being augmented: Physics-based analysis tools for filling the design space database; Distributed computational resources to reduce response time and cost; Web-based technologies to relieve machine-dependence; and Artificial intelligence technologies to accelerate processes and reduce process variability. Activities such as the Advanced Design Technologies Testbed (ADTT) project at NASA Ames Research Center study the strengths and weaknesses of the technologies supporting each of these trends, as well as the overall impact of the combination of these trends on a product development event. Lessons learned and recommendations for future activities will be reported.

Description: The overall objective of this work is to characterize the flow to rough showerheads by deriving pressure drop versus velocity correlations to at can be then used in reactor scale simulations where the showerhead is approximated as a porous medium. At relatively low Reynolds numbers (less than 1-10 based on the hole length scale) and in the absence of slip flow, Darcy's Law, grad P = mu U/k, can be used to express the relation between the pressure drop and velocity where @mu@ is the fluid viscosity and it is the permeability that can be theoretically predicted as k= e R^2 /8, where e is the porosity. However, at sufficiently small hole diameters and decreased pressures (less than 5 Torr), the Knudsen number based on showerhead tube radius increases, and the flow may be in a transition regime. Different expressions have been proposed to account for this effect in the permeability by expressing k as a function of either pressure or Knudsen number. But at even higher Knudsen numbers, the pressure drop - velocity dependence is non-linear, and Darcy's Law no longer holds such that a permeability cannot be defined. The direct simulation Monte Carlo method is used along side conventional CFD techniques to determine the extent to which the CFD technique is appropriate and helps to derive correlations for the more rarefied cases of interest in these showerhead flows.

Description: This conference presents information to the scientific community on research results, future directions, and research opportunities in microgravity fluid physics and transport phenomena within NASA's microgravity research program. The conference theme is "The International Space Station." Plenary sessions provide an overview of the Microgravity Fluid Physics Program, the International Space Station and the opportunities ISS presents to fluid physics and transport phenomena researchers, and the process by which researchers may become involved in NASA's program, including information about the NASA Research Announcement in this area. Two plenary lectures present promising areas of research in electrohydrodynamics/electrokinetics in the movement of particles and in micro- and meso-scale effects on macroscopic fluid dynamics. Featured speakers in plenary sessions present results of recent flight experiments not heretofore presented. The conference publication consists of this book of abstracts and the full Proceedings of the 4th Microgravity Fluid Physics and Transport Phenomena Conference on CD-ROM, containing full papers presented at the conference (NASA/CP-1999-208526/SUPPL1).

Description: Following is a brief summary of a theoretical investigation of material (or constitutive) instability associated with shear induced particle migration in dense particulate suspensions or granular media. It is shown that one can obtain a fairly general linear-stability analysis, including the effects of shear-induced anisotropy in the base flow as well as Reynolds dilatancy. A criterion is presented here for simple shearing instability in the absence of inertia and dilatancy.

Description: Data assimilation methods are routinely used in oceanography. The statistics of the model and measurement errors need to be specified a priori. This study addresses the problem of estimating model and measurement error statistics from observations. We start by testing innovation based methods of adaptive error estimation with low-dimensional models in the North Pacific (5-60 deg N, 132-252 deg E) to TOPEX/POSEIDON (TIP) sea level anomaly data, acoustic tomography data from the ATOC project, and the MIT General Circulation Model (GCM). A reduced state linear model that describes large scale internal (baroclinic) error dynamics is used. The methods are shown to be sensitive to the initial guess for the error statistics and the type of observations. A new off-line approach is developed, the covariance matching approach (CMA), where covariance matrices of model-data residuals are "matched" to their theoretical expectations using familiar least squares methods. This method uses observations directly instead of the innovations sequence and is shown to be related to the MT method and the method of Fu et al. (1993). Twin experiments using the same linearized MIT GCM suggest that altimetric data are ill-suited to the estimation of internal GCM errors, but that such estimates can in theory be obtained using acoustic data. The CMA is then applied to T/P sea level anomaly data and a linearization of a global GFDL GCM which uses two vertical modes. We show that the CMA method can be used with a global model and a global data set, and that the estimates of the error statistics are robust. We show that the fraction of the GCM-T/P residual variance explained by the model error is larger than that derived in Fukumori et al.(1999) with the method of Fu et al.(1993). Most of the model error is explained by the barotropic mode. However, we find that impact of the change in the error statistics on the data assimilation estimates is very small. This is explained by the large representation error, i.e. the dominance of the mesoscale eddies in the T/P signal, which are not part of the 21 by 1" GCM. Therefore, the impact of the observations on the assimilation is very small even after the adjustment of the error statistics. This work demonstrates that simult&neous estimation of the model and measurement error statistics for data assimilation with global ocean data sets and linearized GCMs is possible. However, the error covariance estimation problem is in general highly underdetermined, much more so than the state estimation problem. In other words there exist a very large number of statistical models that can be made consistent with the available data. Therefore, methods for obtaining quantitative error estimates, powerful though they may be, cannot replace physical insight. Used in the right context, as a tool for guiding the choice of a small number of model error parameters, covariance matching can be a useful addition to the repertory of tools available to oceanographers.

Description: Turbulence attenuation by greater than a factor of two has been observed in many practical gas flows carrying volume fractions as small as 0.01% of dispersed particles. Particles which cause such attenuation usually are smaller than the smallest scales of the turbulence and have time constants 5 to 10 times greater than the time scale of a typical turbulent eddy. That is, strongly attenuating particles usually have Stokes numbers in the range of 5 to 10, indicating that they do not respond to the turbulent fluctuations, but instead just fall through the flow responding only to the mean flow. There are two mechanisms by which free falling particles may attenuate turbulence. First, the unresponsive particles act as a drag on the turbulent eddies, passing energy from the turbulent eddies to the small scale wakes of the particles where it is quickly dissipated by viscosity. The second mechanism is more complicated. Particles falling under gravity convert gravitational potential energy to turbulent velocity fluctuations. If the particles are large, this mechanism increases the overall turbulence level. However, with moderate size particles, the small scale turbulence generated apparently distorts the turbulent eddies leading to more rapid dissipation. Unfortunately, this conclusion is supported only by circumstantial evidence to date. The objectives of the experiment are to use microgravity to separate the two mechanisms. A region of nearly-isotropic decaying turbulence with zero mean flow will be formed in a box in the microgravity environment. Different sets of particles with Stokes numbers in the range of 2 to 20 will be dispersed in the flow. With zero gravity and no mean fluid velocity the particles will have zero mean velocity. With the large Stokes numbers, the fluctuating velocities will also be small. Therefore, the only attenuation mechanism will be the direct action of the particles on the turbulence. Control experiments will also be done in which the particles fall through the measurement volume. Measurements will be acquired using a high resolution image velocimetry (PIV) system being developed specifically for work in particle-laden flows. The measurements will include the decay of the turbulence kinetic energy under various particle loadings. The spatial spectra of the turbulence will also be measured. In a second set of experiments, the interaction of a single eddy with a collection of nearly stationary particles will be examined. The eddy will be a vortex ring emitted by a jet pulse through an orifice. The distortion of the vortex under the influence of the particles will be examined to gain a better understanding of how fine particles can cause such large reductions in turbulence levels. This experiment could not be conducted in terrestrial gravity because the high particle velocities would overwhelm the relatively low speed motion of the vortex ring. This experimental program is just getting underway. The initial challenge is to build a closed facility containing reasonably homogeneous and isotropic turbulence with zero mean velocity. Our approach is to use a set of synthetic jets mounted on the periphery of a transparent plexiglass box to create the turbulence. A synthetic jet is a plenum chamber with an orifice open to the volume of interest. The volume of the chamber fluctuates periodically so alternately a jet is ejected from the volume or flow is drawn back in as a sink. The asymmetry of this situation results in a net transport of momentum and kinetic energy into the volume of interest. The present apparatus includes eight synthetic jets each powered independently by a six inch loudspeaker. The synthetic jets discharge through ejector tubes to increase the scale of the turbulence. Construction of the apparatus is now complete and preliminary flow visualization studies have been conducted. The PIV system is also under development. A compact dual-pulse YAG laser has been acquired as the light source and special software is under development to allow simultaneous measurements of both the particle phase and the fluid phase (marked by fine tracers).

Description: The Aeroelasticity Task is intended to provide demonstrated technology readiness to predict and improve flutter characteristics of an HSCT configuration. This requires aerodynamic codes that are applicable to the wide range of flight regimes in which the HSCT will operate, and are suitable to provide the higher fidelity required for evaluation of aeroservoelastic coupling effects. Prediction of these characteristics will result in reduced airplane weight and risk associated with a highly flexible, low-aspect ratio supersonic airplane with narrow fuselage, relatively thin wings, and heavy engines. This Task is subdivided into three subtasks. The first subtask includes the design, fabrication, and testing of wind-tunnel models suitable to provide an experimental database relevant to HSCT configurations. The second subtask includes validation of candidate unsteady aerodynamic codes, applicable in the Mach and frequency ranges of interest for the HSCT, through analysis test correlation with the test data. The third subtask includes efforts to develop and enhance these codes for application to HSCT configurations. The wind tunnel models designed and constructed during this program furnished data which were useful for the analysis test correlation work but there were shortcomings. There was initial uncertainty in the proper tunnel configuration for testing, there was a need for higher quality measured model geometry, and there was a need for better measured model displacements in the test data. One of the models exhibited changes in its dynamic characteristics during testing. Model design efforts were hampered by a need for more and earlier analysis support and better knowledge of material properties. Success of the analysis test correlation work was somewhat muted by the uncertainties in the wind tunnel model data. The planned extent of the test data was not achieved, partly due to the delays in the model design and fabrication which could not be extended due to termination of the HSR program.

Description: The local viscous-inviscid interaction field generated by a wall temperature jump on a flat plate in supersonic flow and on the windside of a Reusable Launch Vehicle in hypersonic flow is studied in detail by both a Navier-Stokes numerical code and an analytical triple-deck model. Treatment of the rapid heat transfer changes both upstream and downstream of the jump is included. Closed form relationships derived from the triple-deck theory are presented. The analytically predicted pressure and heating variations including upstream influence are found to be in generally good agreement with the Computational Fluid Dynamic (CFD) predictions. These analyses not only clarify the interactive physics involved but also are useful in preliminary design of thermal protection systems and as an insertable module to improve CFD code efficiency when applied to such small-scale interaction problems. The analyses only require conditions at the wall and boundary-layer edge which are easily extracted from a baseline, constant wall temperature, CFD solution.

Description: The Line Integral Convolution (LIC) algorithm has received a lot of attention and interest. Yet, only a few of the current LIC related algorithms deal specifically with color textures for automatic detection of flow features. This paper provides an overview of research in this area.

Description: The present invention provides an ozone generation and delivery system that lends itself to small scale applications and requires very low maintenance. The system includes an anode reservoir and a cathode phase separator each having a hydrophobic membrane to allow phase separation of produced gases from water. The system may be configured to operate passively with no moving parts or in a self-pressurizing manner with the inclusion of a pressure controlling device or valve in the gas outlet of the anode reservoir. The hydrogen gas, ozone gas and water containing ozone may be delivered under pressure.

Description: Gas bubbles driven in radial oscillations are subject to an instability of the spherical shape that is opposed by surface tension and viscosity. An exact linear formulation for the study of the phenomenon has been available for many years, but its complexity has discouraged a detailed investigation. With the recent theory of sonoluminescence of Lohse and co-workers, there has been a renewed interest in the problem and new data have become available. This paper presents a numerical method for the solution of the pertinent equations and compares the theory with these new data. The coupling of the strong nonlinearity of the bubble radial oscillations with the parametric mechanism of the surface instability results in a very complex structure for the stability boundary. Nevertheless, a good agreement between theory and data is found. A comparison with earlier approximate models is also made.

Description: A steady, two dimensional cellular convection modifies the morphological instability of a binary alloy that undergoes directional solidification. When the convection wavelength is far longer than that of the morphological cells, the behavior of the moving front is described by a slow, spatial-temporal dynamics obtained through a multiple-scale analysis. The resulting system has a "parametric-excitation" structure in space, with complex parameters characterizing the interactions between flow, solute diffusion, and rejection. The convection stabilizes two dimensional disturbances oriented with the flow, but destabilizes three dimensional disturbances in general. When the flow is weak, the morphological instability behaves incommensurably to the flow wavelength, but becomes quantized and forced to fit into the flow-box as the flow gets stronger. At large flow magnitudes the instability is localized, confined in narrow envelopes with cells traveling with the flow. In this case the solutions are discrete eigenstates in an unbounded space. Their stability boundary and asymptotics are obtained by the WKB analysis.

Description: Foams are extremely important in a variety of industrial applications. Foams are widely used in fire-fighting applications, and are especially effective in fighting flammable liquid fires. In fact the Fire Suppression System aboard the Space Shuttle utilizes cylinders of Halon foam, which, when fired, force a rapidly expanding foam into the convoluted spaces behind instrument panels. Foams are critical in the process of enhanced oil recovery, due to their surface-active and highly viscous nature. They are also used as drilling fluids in underpressurized geologic formations. They are used as transport agents, and as trapping agents. They are also used as separation agents, where ore refinement is accomplished by froth flotation of the typically lighter and hydrophobic contaminants. The goal of the proposed investigation is the determination of the mechanical and rheological properties of foams, utilizing the microgravity environment to explore foam rheology for foams which cannot exist, or only exist for a short time, in 1g.

Description: Integrated modeling of spacecraft systems is a rapidly evolving area in which multidisciplinary models are developed to design and analyze spacecraft configurations. These models are especially important in the early design stages where rapid trades between subsystems can substantially impact design decisions. Integrated modeling is one of the cornerstones of two of NASA's planned missions in the Origins Program -- the Next Generation Space Telescope (NGST) and the Space Interferometry Mission (SIM). Common modeling tools for control design and opto-mechanical analysis have recently emerged and are becoming increasingly widely used. A discipline that has been somewhat less integrated, but is nevertheless of critical concern for high precision optical instruments, is thermal analysis and design. A major factor contributing to this mild estrangement is that the modeling philosophies and objectives for structural and thermal systems typically do not coincide. Consequently the tools that are used in these discplines suffer a degree of incompatibility, each having developed along their own evolutionary path. Although standard thermal tools have worked relatively well in the past. integration with other disciplines requires revisiting modeling assumptions and solution methods. Over the past several years we have been developing a MATLAB based integrated modeling tool called IMOS (Integrated Modeling of Optical Systems) which integrates many aspects of structural, optical, control and dynamical analysis disciplines. Recent efforts have included developing a thermal modeling and analysis capability, which is the subject of this article. Currently, the IMOS thermal suite contains steady state and transient heat equation solvers, and the ability to set up the linear conduction network from an IMOS finite element model. The IMOS code generates linear conduction elements associated with plates and beams/rods of the thermal network directly from the finite element structural model. Conductances for temperature varying materials are accommodated. This capability both streamlines the process of developing the thermal model from the finite element model, and also makes the structural and thermal models compatible in the sense that each structural node is associated with a thermal node. This is particularly useful when the purpose of the analysis is to predict structural deformations due to thermal loads. The steady state solver uses a restricted step size Newton method, and the transient solver is an adaptive step size implicit method applicable to general differential algebraic systems. Temperature dependent conductances and capacitances are accommodated by the solvers. In addition to discussing the modeling and solution methods. applications where the thermal modeling is "in the loop" with sensitivity analysis, optimization and optical performance drawn from our experiences with the Space Interferometry Mission (SIM), and the Next Generation Space Telescope (NGST) are presented.

Description: The application of an unstructured grid methodology on a three-dimensional high-lift configuration is presented. The focus of this paper is on the grid generation aspect of an integrated effort for the development of an unstructured-grid computational fluid dynamics (CFD) capability at the NASA Langley Research Center. The meshing approach is based on tetrahedral grids generated by the advancing-front and the advancing-layers procedures. The capability of the method for solving high-lift problems is demonstrated on an aircraft model referred to as the energy efficient transport configuration. The grid generation issues, including the pros and cons of the present approach, are discussed in relation to the high-lift problems. Limited viscous flow results are presented to demonstrate the viability of the generated grids. A corresponding Navier-Stokes solution capability, along with further computations on the present grid, is presented in a companion SAE paper.

Description: The Thermal Radiation Group at Virginia Polytechnic Institute and State University has been working closely with scientists and engineers at NASA's Langley Research Center to develop accurate analytical and numerical models suitable for designing next generation thin-film thermal radiation detectors for earth radiation budget measurement applications. The current study provides an analytical model of the notional thermal radiation detector that takes into account thermal transport phenomena, such as the contact resistance between the layers of the detector, and is suitable for use in parameter estimation. It was found that the responsivity of the detector can increase significantly due to the presence of contact resistance between the layers of the detector. Also presented is the effect of doping the thermal impedance layer of the detector with conducting particles in order to electrically link the two junctions of the detector. It was found that the responsivity and the time response of the doped detector decrease significantly in this case. The corresponding decrease of the electrical resistance of the doped thermal impedance layer is not sufficient to significantly improve the electrical performance of the detector. Finally, the "roughness effect" is shown to be unable to explain the decrease in the thermal conductivity often reported for thin-film layers.

Description: A method and apparatus for cold gas reinjection in through-flow and reverse-flow wave rotors having a plurality of channels formed around a periphery thereof. A first port injects a supply of cool air into the channels. A second port allows the supply of cool air to exit the channels and flow to a combustor. A third port injects a supply of hot gas from the combustor into the channels. A fourth port allows the supply of hot gas to exit the channels and flow to a turbine. A diverting port and a reinjection port are connected to the second and third ports, respectively. The diverting port diverts a portion of the cool air exiting through the second port as reinjection air. The diverting port is fluidly connected to the reinjection port which reinjects the reinjection air back into the channels. The reinjection air evacuates the channels of the hot gas resident therein and cools the channel walls, a pair of end walls of the rotor, ducts communicating with the rotor and subsequent downstream components. In a second embodiment, the second port receives all of the cool air exiting the channels and the diverting port diverts a portion of the cool air just prior to the cool air flowing to the combustor.

Description: A fractional step method for the solution of steady and unsteady incompressible Navier-Stokes equations is outlined. The method is based on a finite volume formulation and uses the pressure in the cell center and the mass fluxes across the faces of each cell as dependent variables. Implicit treatment of convective and viscous terms in the momentum equations enables the numerical stability restrictions to be relaxed. The linearization error in the implicit solution of momentum equations is reduced by using three subiterations in order to achieve second order temporal accuracy for time-accurate calculations. In spatial discretizations of the momentum equations, a high-order (3rd and 5th) flux-difference splitting for the convective terms and a second-order central difference for the viscous terms are used. The resulting algebraic equations are solved with a line-relaxation scheme which allows the use of large time step. A four color ZEBRA scheme is employed after the line-relaxation procedure in the solution of the Poisson equation for pressure. This procedure is applied to a Couette flow problem using a distorted computational grid to show that the method minimizes grid effects. Additional benchmark cases include the unsteady laminar flow over a circular cylinder for Reynolds Numbers of 200, and a 3-D, steady, turbulent wingtip vortex wake propagation study. The solution algorithm does a very good job in resolving the vortex core when 5th-order upwind differencing and a modified production term in the Baldwin-Barth one-equation turbulence model are used with adequate grid resolution.

Description: This paper analytically examines the unsteady fluid dynamics of a vortex filament subsequent to a normal collision of the vortex with a solid body. In particular, the breakdown or reconnection phenomena, post-collision, for a vortex filament is studied. The paper does not investigate the collision dynamics process itself. The derived exact solution is based upon the laminar viscous form of the Helmholtz equations.

Description: In spite of a superficial similarity with gas bubbles, the intimate coupling between dynamical and thermal processes confers to oscillating vapor bubbles some unique characteristics. This paper examines numerically the validity of some asymptotic-theory predictions such as the existence of two resonant radii and a limit size for a given sound amplitude and frequency. It is found that a small vapor bubble in a sound field of sufficient amplitude grows quickly through resonance and continues to grow thereafter at a very slow rate, seemingly indefinitely. Resonance phenomena therefore play a role for a few cycles at most, and reaching a limit size-if one exists at all-is found to require far more than several tens of thousands of cycles. It is also found that some small bubbles may grow or collapse depending on the phase of the sound field. The model accounts in detail for the thermo-fluid-mechanic processes in the vapor. In the second part of the paper, an approximate formulation valid for bubbles small with respect to the thermal penetration length in the vapor is derived and its accuracy examined, The present findings have implications for acoustically enhanced boiling heat transfer and other special applications such as boiling in microgravity.

Description: Condition queries on distributed data ask where particular conditions are satisfied. It is possible to represent condition queries as geometric objects by plotting field data in various spaces derived from the data, and by selecting loci within these derived spaces which signify the desired conditions. Rather simple geometric partitions of derived spaces can represent complex condition queries because much complexity can be encapsulated in the derived space mapping itself A geometric view of condition queries provides a useful conceptual unification, allowing one to intuitively understand many existing vector field feature detection algorithms -- and to design new ones -- as variations on a common theme. A geometric representation of condition queries also provides a simple and coherent basis for computer implementation, reducing a wide variety of existing and potential vector field feature detection techniques to a few simple geometric operations.

Description: This paper describes an experiment in which a large-scale scientific application development for tightly-coupled parallel machines is adapted to the distributed execution environment of the Information Power Grid (IPG). A brief overview of the IPG and a description of the computational fluid dynamics (CFD) algorithm are given. The Globus metacomputing toolkit is used as the enabling device for the geographically-distributed computation. Modifications related to latency hiding and Load balancing were required for an efficient implementation of the CFD application in the IPG environment. Performance results on a pair of SGI Origin 2000 machines indicate that real scientific applications can be effectively implemented on the IPG; however, a significant amount of continued effort is required to make such an environment useful and accessible to scientists and engineers.

Description: A compact modulated birefringence-measurement system has been developed for use in microgravity fluid physics applications with non-Newtonian fluids such as polymer solutions. This instrument uses a dual-crystal transverse electro-optical modulator capable of modulation frequencies in excess of 100 MHz. The two crystals are modulated 180 deg. out of phase from each other, The theoretical framework governing the development of this instrument using the Mueller-Stokes polarization matrices is discussed. Several ground-based experiments are performed to compare this system with the theoretical results. Results from this transverse electro-optical modulator-based birefringence-measurement system agree well with the theory. The instrument is also very stable and robust, making it suitable for the extreme acceleration environment to be encountered in a NASA Black Brandt sounding rocket.

Description: The creeping flow of a dilute (0.025 wt%) monodisperse polystyrene/polystyrene Boger fluid through a 4:1:4 axisymmetric contraction/expansion is experimentally observed for a wide range of Deborah numbers. Pressure drop measurements across the orifice plate show a large extra pressure drop that increases monotonically with Deborah number above the value observed for a similar Newtonian fluid at the same flow rate. This enhancement in the dimensionless pressure drop is not associated with the onset of a flow instability, yet it is not predicted by existing steady-state or transient numerical computations with simple dumbbell models. It is conjectured that this extra pressure drop is the result of an additional dissipative contribution to the polymeric stress arising from a stress-conformation hysteresis in the strong non-homogeneous extensional flow near the contraction plane. Such a hysteresis has been independently measured and computed in recent studies of homogeneous transient uniaxial stretching of PS/PS Boger fluids. Flow visualization and velocity field measurements using digital particle image velocimetry (DPIV) show large upstream growth of the corner vortex with increasing Deborah number. At large Deborah numbers, the onset of an elastic instability is observed, first locally as small amplitude fluctuations in the pressure measurements, and then globally as an azimuthal precessing of the upstream corner vortex accompanied by periodic oscillations in the pressure drop across the orifice.

Description: Scale-similar models employ multiple filtering operations to identify the smallest resolved scales, which have been shown to be the most active in the interaction with the unresolved subgrid scales. They do not assume that the principal axes of the strain-rate tensor are aligned with those of the subgrid-scale stress (SGS) tensor, and allow the explicit calculation of the SGS energy. They can provide backscatter in a numerically stable and physically realistic manner, and predict SGS stresses in regions that are well correlated with the locations where large Reynolds stress occurs. In this paper, eddy viscosity and mixed models, which include an eddy-viscosity part as well as a scale-similar contribution, are applied to the simulation of two flows, a high Reynolds number plane channel flow, and a three-dimensional, nonequilibrium flow. The results show that simulations without models or with the Smagorinsky model are unable to predict nonequilibrium effects. Dynamic models provide an improvement of the results: the adjustment of the coefficient results in more accurate prediction of the perturbation from equilibrium. The Lagrangian-ensemble approach [Meneveau et al., J. Fluid Mech. 319, 353 (1996)] is found to be very beneficial. Models that included a scale-similar term and a dissipative one, as well as the Lagrangian ensemble averaging, gave results in the best agreement with the direct simulation and experimental data.

Description: Experiments were designed, implemented, and evaluated in the thermal and fluid sciences at the NASA Langley Research Center. This research was conducted cooperatively with NASA employees using, where necessary, equipment and facilities provided by the U.S. Government. The research fell within the scope of the University Agreement between the NASA Langley Research Center and The George Washington University for Joint Research and Education Projects dated June 7, 8, 1994, which continues the Joint Institute for the Advancement of Flight Sciences (JIAFS).

Description: In the largest solar energetic particle (SEP) events, acceleration occurs at shock waves driven out from the Sun by coronal mass ejections (CMEs). Peak particle intensities are a strong function of CME speed, although the intensities, spectra, and angular distributions of particles escaping the shock are highly modified by scattering on Alfven waves produced by the streaming particles themselves. Element abundances vary in complex ways because ions with different values of Q/A resonate with different parts of the wave spectrum, which varies with space and time. Just recently, we have begun to model these systematic variations theoretically and to explore other consequences of proton-generated waves.

Description: Computational Fluid Dynamics (CFD) analysis results are compared with benchmark quality test data from the Propulsion Engineering Research Center's (PERC) Rocket Based Combined Cycle (RBCC) experiments to verify fluid dynamic code and application procedures. RBCC engine flowpath development will rely on CFD applications to capture the multi-dimensional fluid dynamic interactions and to quantify their effect on the RBCC system performance. Therefore, the accuracy of these CFD codes must be determined through detailed comparisons with test data. The PERC experiments build upon the well-known 1968 rocket-ejector experiments of Odegaard and Stroup by employing advanced optical and laser based diagnostics to evaluate mixing and secondary combustion. The Finite Difference Navier Stokes (FDNS) code was used to model the fluid dynamics of the PERC RBCC ejector mode configuration. Analyses were performed for both Diffusion and Afterburning (DAB) and Simultaneous Mixing and Combustion (SMC) test conditions. Results from both the 2D and the 3D models are presented.

Description: In this paper we describe the numerical analysis underlying our efforts to develop an accurate and reliable code for simulating flame propagation using complex physical and chemical models. We discuss our spatial and temporal discretization schemes, which in our current implementations range in order from two to six. In space we use staggered meshes to define discrete divergence and gradient operators, allowing us to approximate complex diffusion operators while maintaining ellipticity. Our temporal discretization is based on the use of preconditioning to produce a highly efficient linearly implicit method with good stability properties. High order for time accurate simulations is obtained through the use of extrapolation or deferred correction procedures. We also discuss our techniques for computing stationary flames. The primary issue here is the automatic generation of initial approximations for the application of Newton's method. We use a novel time-stepping procedure, which allows the dynamic updating of the flame speed and forces the flame front towards a specified location. Numerical experiments are presented, primarily for the stationary flame problem. These illustrate the reliability of our techniques, and the dependence of the results on various code parameters.

Description: We present an Abel-inversion algorithm to reconstruct mean and rms refractive-index profiles from spatially resolved statistical measurements of the beam-deflection angle in time-dependent, axisymmetric flows. An oscillating gas-jet diffusion flame was investigated as a test case for applying the algorithm. Experimental data were obtained across the whole field by a rainbow schlieren apparatus. Results show that simultaneous multipoint measurements are necessary to reconstruct the rms refractive index accurately.

Description: A scale-invariant statistical theory of fields is presented that leads to invariant definition of density, velocity, temperature, and pressure, The definition of Boltzmann constant is introduced as k(sub k) = m(sub k)v(sub k)c = 1.381 x 10(exp -23) J x K(exp -1), suggesting that the Kelvin absolute temperature scale is equivalent to a length scale. Two new state variables called the reversible heat Q(sub rev) = TS and the reversible work W(sub rev) = PV are introduced. The modified forms of the first and second law of thermodynamics are presented. The microscopic definition of heat (work) is presented as the kinetic energy due to the random (peculiar) translational, rotational, and pulsational motions. The Gibbs free energy of an element at scale Beta is identified as the total system energy at scale (Beta-1), thus leading to an invariant form of the first law of thermodynamics U(sub Beta) = Q(sub Beta) - W(sub Beta) +N(e3)U(sub Beta-1).

Description: Mainframe supercomputers such as the Cray C90 was invaluable in obtaining large scale computations using several millions of grid points to resolve salient features of a tip vortex flow over a lifting wing. However, real flight configurations require tracking not only of the flow over several lifting wings but its growth and decay in the near- and intermediate- wake regions, not to mention the interaction of these vortices with each other. Resolving and tracking the evolution and interaction of these vortices shed from complex bodies is computationally intensive. Parallel computing technology is an attractive option in solving these flows. In planetary science vortical flows are also important in studying how planets and protoplanets form when cosmic dust and gases become gravitationally unstable and eventually form planets or protoplanets. The current paradigm for the formation of planetary systems maintains that the planets accreted from the nebula of gas and dust left over from the formation of the Sun. Traditional theory also indicate that such a preplanetary nebula took the form of flattened disk. The coagulation of dust led to the settling of aggregates toward the midplane of the disk, where they grew further into asteroid-like planetesimals. Some of the issues still remaining in this process are the onset of gravitational instability, the role of turbulence in the damping of particles and radial effects. In this study the focus will be with the role of turbulence and the radial effects.

Description: This paper provides an overview of the microgravity fluid physics and transport phenomena experiments planned for the International Spare Station. NASA's Office of Life and Microgravity Science and Applications has established a world-class research program in fluid physics and transport phenomena. This program combines the vast expertise of the world research community with NASA's unique microgravity facilities with the objectives of gaining new insight into fluid phenomena by removing the confounding effect of gravity. Due to its criticality to many terrestrial and space-based processes and phenomena, fluid physics and transport phenomena play a central role in the NASA's Microgravity Program. Through widely publicized research announcement and well established peer-reviews, the program has been able to attract a number of world-class researchers and acquired a critical mass of investigations that is now adding rapidly to this field. Currently there arc a total of 106 ground-based and 20 candidate flight principal investigators conducting research in four major thrust areas in the program: complex flows, multiphase flow and phase change, interfacial phenomena, and dynamics and instabilities. The International Space Station (ISS) to be launched in 1998, provides the microgravity research community with a unprecedented opportunity to conduct long-duration microgravity experiments which can be controlled and operated from the Principal Investigators' own laboratory. Frequent planned shuttle flights to the Station will provide opportunities to conduct many more experiments than were previously possible. NASA Lewis Research Center is in the process of designing a Fluids and Combustion Facility (FCF) to be located in the Laboratory Module of the ISS that will not only accommodate multiple users but, allow a broad range of fluid physics and transport phenomena experiments to be conducted in a cost effective manner.

Description: Two-phase flows under microgravity condition find a large number of important applications in fluid handling and storage, and spacecraft thermal management. Specifically, under microgravity condition heat transfer between heat exchanger surfaces and fluids depend critically on the distribution and interaction between different fluid phases which are often qualitatively different from the gravity-based systems. Heat transfer and flow analysis in two-phase flows under these conditions require a clear understanding of the flow pattern transition and development of appropriate dimensionless scales for its modeling and prediction. The physics of this flow is however very complex and remains poorly understood. This has led to various inadequacies in flow and heat transfer modeling and has made prediction of flow transition difficult in engineering design of efficient thermal and flow systems. In the present study the available published data for flow transition under microgravity condition are considered for mapping. The transition from slug to annular flow and from bubbly to slug flow are mapped using dimensionless variable combination developed in a previous study by the authors. The result indicate that the new maps describe the flow transitions reasonably well over the range of the data available. The transition maps are examined and the results are discussed in relation to the presumed balance of forces and flow dynamics. It is suggested that further evaluation of the proposed flow and transition mapping will require a wider range of microgravity data expected to be made available in future studies.

Description: In order to expand the predictive capability of single-point turbulence closure models to account for the early-stage transition regime, a methodology for the formulation and calibration of model equations for the ensemble-averaged disturbance kinetic energy and energy dissipation rate is presented. The calibration is based on homogeneous shear flow where disturbances can be described by rapid distortion theory (RDT). The relationship between RDT and linear stability theory is exploited in order to obtain a closed set of modeled equations. The linear disturbance equations are solved directly so that the numerical simulation yields a database from which the closure coefficients in the ensemble-averaged disturbance equations can be determined. Key words: transition/turbulence modeling, homogeneous shear, rapid distortion theory.

Description: This paper examines the modeling of two-dimensional homogeneous stratified turbulent shear flows using the Reynolds-stress and Reynolds-heat-flux equations. Several closure models have been investigated-, the emphasis is placed on assessing the effect of modeling the dissipation rate tensor in the Reynolds-stress equation. Three different approaches are considered: one is an isotropic approach while the other two are anisotropic approaches. The isotropic approach is based on Kolmogorov's hypothesis and a dissipation rate equation modified to account for vortex stretching. One of the anisotropic approaches is based on an algebraic representation of the dissipation rate tensor, while another relies on solving a modeled transport equation for this tensor. In addition, within the former anisotropic approach, two different algebraic representations are examined one is a function of the Reynolds-stress anisotropy tensor, and the other is a function of' the mean velocity gradients. The performance of these closure models is evaluated against experimental and direct numerical simulation data of pure shear flows. pure buoyant flows and buoyant shear flows. Calculations have been carried out over a range of Richardson numbers (Ri) and two different Prandtl numbers (Pr); thus the effect of Pr on the development of counter-gradient heat flux in a stratified shear flow can be assessed. At low Ri, the isotropic model performs well in the predictions of stratified shear flows; however, its performance deteriorates as Ri increases. At high Ri, the transport equation model for the dissipation rate tensor gives the best result. Furthermore, the results also lend credence to the algebraic dissipation rate model based on the Reynolds stress anisotropy tensor. Finally, it is found that Pr has an effect on the development of counter-gradient heat flux. The calculations show that, under the action of shear, counter-gradient heat flux does not occur even at Ri = 1 in an air flow.

Description: This conference presents information to the scientific community on research results, future directions, and research opportunities in microgravity fluid physics and transport phenomena within NASA's microgravity research program. The conference theme is "The International Space Station." The conference publication consists of the full Proceedings of the 4th Microgravity Fluid Physics and Transport Phenomena Conference on CD-ROM, containing full papers presented at the conference. Ninety papers are presented in 21 technical sessions, and a special exposition session presents 32 posters describing the work of principal investigators new to NASA's program in this discipline. Eighty-eight papers and 25 posters are presented in their entirety on the CD-ROM.

Description: In early nineties Mars was designated as the planet to be explored. The present results of Mars program include the following: Background of Mars program. Thermal control design for Martian spacecraft. Loop heat pipe for Lander and Rover. Test results. and LHP applications for future Martian Lander and Rover.

Description: We examine, with asymptotic analysis and numerical simulation, the iterated stretching dynamics of FENE and Oldroyd-B jets of initial radius r(sub 0), shear viscosity nu, Weissenberg number We, retardation number S, and capillary number Ca. The usual Rayleigh instability stretches the local uniaxial extensional flow region near a minimum in jet radius into a primary filament of radius [Ca(1 - S)/ We](sup 1/2)r(sub 0) between two beads. The strain-rate within the filament remains constant while its radius (elastic stress) decreases (increases) exponentially in time with a long elastic relaxation time 3We(r(sup 2, sub 0)/nu). Instabilities convected from the bead relieve the tension at the necks during this slow elastic drainage and trigger a filament recoil. Secondary filaments then form at the necks from the resulting stretching. This iterated stretching is predicted to occur successively to generate high-generation filaments of radius r(sub n), (r(sub n)/r(sub 0)) = square root of 2[r(sub n-1)/r(sub 0)](sup 3/2) until finite-extensibility effects set in.

Description: A reduced gravity environment offers the ability to remove the effect of buoyancy on two phase flows whereby density differences that normally would promote relative velocities between the phases and also alter the shape of the interface are removed. However, besides being a potent research tool, there are also many space-based technologies that will either utilize or encounter two-phase flow behavior, and as a consequence, several questions must be addressed. This paper presents some of these technologies missions. Finally, this paper gives a description of web-sites for some funding.

Description: The contact angle and the spreading process of sessile droplet are very crucial in many technological processes, such as painting and coating, material processing, film-cooling applications, lubrication, and boiling. Additionally, as it is well known that the surface free energy of polymers cannot be directly, measured for their elastic and viscous restraints. The measurements of liquid contact angle on the polymer surfaces become extremely important to evaluate the surface free energy of polymers through indirect methods linked with the contact angle data. Due to the occurrence of liquid evaporation is inevitable, the effects of evaporation on the contact angle and the spreading become very important for more complete understanding of these processes. It is of interest to note that evaporation can induce Marangoni-Benard convection in sessile drops. However, the impacts of the inside convection on the wetting and spreading processes are not clear. The experimental methods used by previous investigators cannot simultaneously measure the spreading process and visualize the convection inside. Based on the laser shadowgraphic system used by the present author, a very simple optical procedure has been developed to measure the contact angle, the spreading speed, the evaporation rate, and to visualize inside convection of a sessile drop simultaneously. Two CCD cameras were used to synchronously record the real-time diameter of the sessile drop, which is essential for determination of both spreading speed and evaporation rate, and the shadowgraphic image magnified by the sessile drop acting as a thin plano-convex lens. From the shadowgraph, the inside convection of the drop can be observed if any and the image outer diameter, which linked to the drop profile, can be measured. Simple equations have been derived to calculate the drop profile, including the instantaneous contact angle, height, and volume of the sessile drop, as well as the evaporation rate. The influence of the inside convection on the wetting and spreading processes can be figured out through comparison of the drop profiles with and without inside convection when the sessile drop is placed at different evaporation conditions.

Description: The current study is characterized by two distinct phases in the development of the vortex tube (VT) technology as a primary means for in-flight air separation. The purpose of the first phase was to systematically identify parameters that influence oxygen concentration and recovery and to quantify the extent of that influence. To that end, the project team used a series of planned factorial experiments to identify statistically significant variables (factors) and their interactions. These experiments identified a best range of the operating envelope that includes nozzle diameter, orifice diameter, inlet air pressure, and liquid phase content in the inlet air. The best results observed in this envelope were an oxygen content of approximately 68% and a recovery factor of approximately 38%. The primary objectives of the second phase of the current study were to investigate the application effects of the two different air separation efficiency enhancement methods. One of these methods resulted in a concentration increase of 12% and second resulted in a concentration increase of 5%. Several aspects of these methods application are subject to optimize.

Description: The objectives of this research are: (1) to develop and implement a new methodology for large eddy simulation of (LES) of high-speed reacting turbulent flows. (2) To develop algebraic turbulence closures for statistical description of chemically reacting turbulent flows. We have just completed the third year of Phase III of this research. This is the Final Report of our activities on this research sponsored by the NASA LaRC.

Description: An experiment on a cooled swept cylinder in a low-disturbance Mach 1.6 wind tunnel is described. The flow attachment line is disturbed by trip wires of varying size and the laminar/turbulent state of the downstream boundary layer is determined with a hot wire. The results demonstrate that although cooling the wall increases the stability of the boundary layer, it promotes roughness induced transition. Analysis of the data suggests that the attachment- line Reynolds number can account for the effect of wall cooling if the viscosity is evaluated at a particular reference temperature.

Description: This paper describes the unsteady blade surface pressures on the first-stage rotor blades of a two-stage transonic axial flow compressor experiencing inlet flow distortion. This study was conducted to demonstrate the ability of a full annulus unsteady Reynolds-averaged Navier-Stokes numerical technique to predict unsteady pressures on the rotor blades operating in a distorted inflow. A total pressure distortion produced by a variable mesh screen mounted near the inlet was used to excite the unsteady blade loading on the rotor. On-blade pressure transducers were used to measure the unsteady blade surface pressure. These pressures and the resulting transient load on the rotor blades were compared to the numerical prediction. It is important to develop numerical techniques to predict these transient loads to better understand the response of compressor blades to forcing functions. With this enhanced understanding and ability to predict these transient forces, more robust compressors can be developed. In the study, a high definition of the inlet flow distortion was achieved by rotating the distortion screens. In this manner the inlet flow distortion and the distortion at the first stage stator leading edge were measured at approximately every 0.7 degrees. This full annulus high definition of the inlet flow distortion was used as the inlet boundary condition for the numerical technique. The experimental measurements and numerical analyses are highly complementary in this study. Detailed comparisons between the measurements and the numerical analyses indicate that the current numerical procedure calculates the unsteady aerodynamic pressure on the blade surfaces reasonably well. Further, the agreement of the measured and predicted rotor exit flow distortion at the first stage stator leading edge provides verification of the numerical technique.

Description: An energy equation is derived for the general case of a viscous drop suspended in a viscous medium with surfactants contaminating the interface. It contains terms that clearly identify dissipation contributions from the viscous effects in the bulk fluids, surface shear and dilatational viscosity effects at the interface, and surfactant transport. An efficient boundary integral method is developed which incorporates the effects of a constant surface dilatational viscosity in simulations of an oscillating two-dimensional inviscid drop. Surface dilatational viscosity is shown to have a significant damping effect on the otherwise undamped inviscid oscillations.

Description: We discuss results from two parts of our study on the behavior of liquids under low-gravity conditions. The first concerns the Interface Configuration Experiment (ICE) aboard the Space Station Mir on the Mir-21/NASA-2 mission; for a certain 'exotic' container, distinct asymmetric liquid configurations are found as locally stable ones, even though the container itself is rotationally symmetric, in confirmation of mathematical results and numerical computations. The second investigation concerns the behavior of slowly rotating liquids; it is found that a rotating film instability observed previously in a physical experiment in 1-g, scaled to render gravity effects small, does not correspond to mathematical and computational results obtained for low gravity. These latter results are based on the classical equilibrium theory enhanced with a van der Waals potential of adhesion.

Description: The existing theoretical research concerned with thermal creep at fluid/solid interfaces is briefly reviewed, and the importance of microgravity-based experimental data is then discussed. It is noted that the ultimate goal of this research is a rational molecular level theory that predicts the dependence of a dimensionless thermal creep coefficient, Ctc, on relevant dimensionless parameters describing the way fluid molecules interact with the solid surface and how they interact among themselves. The discussion covers thermophoresis of isolated solid spheres and aggregates in gases; solid sphere thermophoresis in liquids and dense vapors; thermophoresis of small immiscible liquid droplets; and applications of the direct simulation Monte Carlo method.

Description: Results of three problems are summarized in this contribution. Each involves the fundamental capillary instability of an interfacial bridge and is an extension of previous work. The first two problems concern equilibrium shapes of liquid bridges near the stability boundary corresponding to maximum length (Plateau-Rayleigh limit). For the first problem, a previously formulated nonlinear theory to account for imposed gravity and interfacial shear disturbances in an isothermal environment is quantitatively tested in experiment. For the second problem, the liquid bridge is subjected to a shear that models the effect of a thermocapillary flow generated by a ring heater in a liquid encapsulated float-zone configuration. In the absence of gravity, this symmetric perturbation can stabilize the bridge to lengths on the order of 30 percent beyond the Plateau-Rayleigh limit, which is on the order of heretofore unexplained Shuttle observations. The third problem considers the dynamics of collapse and pinchoff of a film bridge (no gravity), which happens in the absence of stabilization. Here, we summarize experimental efforts to measure the self-similar cone-and-crater structure predicted by a previous theory.

Description: The use of both acoustic radiation stress and the Maxwell stress to stabilize liquid bridges is reported. Acoustic radiation stress arises from the time-averaged acoustic pressure at the surface of an object immersed in a sound field. Both passive and active acoustic stabilization schemes as well as an active electrostatic method are examined.

Description: Liquid bridges have been the focus of numerous theoretical and experimental investigations since the early work by Plateau more than a century ago. More recently, motivated by interest in their physical behavior and their occurrence in a variety of technological situations, there has been a resurgence of interest in the static and dynamic behavior of liquid bridges. Furthermore, opportunities to carry out experiments in the near weightless environment of a low-Earth-orbit spacecraft have also led to a number of low-gravity experiments involving large liquid bridges. In this paper, we present selected results from our work concerning the stability of nonaxisymmetric liquid bridges, the bifurcation of weightless bridges in the neighborhood of the maximum volume stability limit, isorotating axisymmetric bridges contained between equidimensional disks, and bridges contained between unequal disks. For the latter, we discuss both theoretical and experimental results. Finally, we present results concerning the stability of axisymmetric equilibrium configurations for a capillary liquid partly contained in a closed circular cylinder.