ALBERT

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Description: It is well established that residual flows exist in contained liquid metal processes. In 1-g processing, buoyancy forces often drive these flows and their magnitudes can be substantial. It is also known that residual flows can exist during microgravity processing, and although greatly reduced in magnitude, they can influence the properties of the processed materials. Unfortunately, there are very few techniques to visualize flows in opaque, high temperature liquid metals, and those available are not easily adapted to flight investigation. In this study, a novel technique is developed that uses liquid tin as the model fluid and solid-state electrochemical cells constructed from Yttria-Stabilized Zirconia (YSZ) to establish and measure dissolved oxygen boundary conditions. The melt serves as a common electrode for each of the electrochemical cells in this design, while independent reference electrodes are maintained at the outside surfaces of the electrolyte. By constructing isolated electrochemical cells at various locations along the container walls, oxygen is introduced or extracted by imposing a known electrical potential or passing a given current between the melt and the reference electrode. This programmed titration then establishes a known oxygen concentration boundary condition at the selected electrolyte-melt interface. Using the other cells, the concentration of oxygen at the electrolyte-melt interface is also monitored by measuring the open-circuit potentials developed between the melt and reference electrodes. Thus the electrochemical cells serve to both establish boundary conditions for the passive tracer and sense its path. Rayleigh-Benard convection was used to validate the electrochemical approach to flow visualization. Thus, a numerical characterization of the second critical Rayleigh numbers in liquid tin was conducted for a variety of Cartesian aspect ratios. The extremely low Prandtl number of tin represents the lowest value studied numerically. Additionally, flow field oscillations are visualized and the effect of tilt on convecting systems is quantified. Experimental studies of the effect of convection in liquid tin are presented. Three geometries are studied: (1) double electrochemical cell with vertical concentration gradients; (2) double cell with horizontal concentration gradients; and (3) multiple cells with vertical temperature gradients. The first critical Rayleigh number transition is detected with geometry (1) and it is concluded that current measurements are not as affected by convection as EMF measurements. The system is compared with numerical simulations in geometry (2), and oscillating convection is detected with geometry (3).

Description: A magnetohydrodynamic model that examines the effect of rotating an electrically conducting cylinder with a uniform external magnetic field applied orthogonal to its axis is presented. Noting a simple geometry, it can be classified as a fundamental dynamo problem. For the case of an infinitely long cylinder, an analytical solution is obtained and analyzed in detail. A semi-analytical model was developed that considers a finite cylinder. Experimental data from a spinning brass wheel in the presence of Earth's magnetic field were compared to the proposed theory and found to fit well.

Description: A theoretical and experimental study is presented on the stability of solutal convection of a magnetized fluid in the presence of a magnetic field. The total force on the fluid is derived and equilibrium positions are computed establishing the field necessary to counter fluid buoyancy. The requirements for stability are developed and compared with experiments with a paramagnetic fluid. The experiments are in good agreement not only with the theoretical predictions for equilibrium but also verify the stability theory which predicts both horizontal and vertical stability. Analogous to results for levitation, the theory indicates that solutal convection in paramagnetic fluids cannot be completely stabilized while that in diamagnetic liquid are possible.

Description: The purpose of our project is to develop, analyze, and test novel numerical technologies central to the long term goal of direct simulations of subsonic jet noise. Our current focus is on two issues: accurate, near-field domain truncations and high-order, single-step discretizations of the governing equations. The Direct Numerical Simulation (DNS) of jet noise poses a number of extreme challenges to computational technique. In particular, the problem involves multiple temporal and spatial scales as well as flow instabilities and is posed on an unbounded spatial domain. Moreover, the basic phenomenon of interest, the radiation of acoustic waves to the far field, involves only a minuscule fraction of the total energy. The best current simulations of jet noise are at low Reynolds number. It is likely that an increase of one to two orders of magnitude will be necessary to reach a regime where the separation between the energy-containing and dissipation scales is sufficient to make the radiated noise essentially independent of the Reynolds number. Such an increase in resolution cannot be obtained in the near future solely through increases in computing power. Therefore, new numerical methodologies of maximal efficiency and accuracy are required.

Description: Development of hydrogen in sealed silica glass ampoules during annealing at elevated temperatures was investigated. The dependence of hydrogen pressure in the ampoules as a function of time, for different temperatures and ampoule parameters was measured. The process was modeled assuming chemical solution of hydrogen according to the reaction: silica + H2 = H- Si= + H-O-Si=. The equilibrium constant of the reaction was determined by fitting the theoretical curves to the experimental data. The Gibbs function for this reaction was estimated at deltaG = -25.8 + 54T.

Description: Volumetric forces resulting from local density variations and gravitational acceleration cause buoyancy induced convective motion in melts and solutions. Solutal buoyancy is a result of concentration differences in an otherwise isothermal fluid. If the fluid also exhibits variations in magnetic susceptibility with concentration then convection control by external magnetic fields can be hypothesized. Magnetic control of thermal buoyancy induced convection in ferrofluids (dispersions of ferromagnetic particles in a carrier fluid) and paramagnetic fluids have been demonstrated. Here we show the nature of magnetic control of solutal buoyancy driven convection of a paramagnetic fluid, an aqueous solution of Manganese Chloride hydrate. We predict the critical magnetic field required for balancing gravitational solutal buoyancy driven convection and validate it through a simple experiment. We demonstrate that gravity driven flow can be completely reversed by a magnetic field but the exact cancellation of the flow is not possible. This is because the phenomenon is unstable. The technique can be applied to crystal growth processes in order to reduce convection and to heat exchanger devices for enhancing convection. The method can also be applied to impose a desired g-level in reduced gravity applications.

Description: Immiscible fluid layers are found in a host of applications ranging from materials processing, for example the use of encapsulants in float zone crystal growth technique and a buffer layer in industrial Czochralski growth of crystals to prevent Marangoni convection, to heat transfer phenomena in day-to-day processes like the presence of air pockets in heat exchangers. In the microgravity and space processing realm, the exploration of other planets requires the development of enabling technologies in several fronts. The reduction in the gravity level poses unique challenges for fluid handling and heat transfer applications. The present work investigates the efficacy of controlling thermal convective flow using magnetic fluids and magnetic fields. The setup is a two-layer immiscible liquids system with one of the fluids being a diluted ferrofluid (super paramagnetic nano particles dispersed in carrier fluid). Using an external magnetic field one can essentially dial in a volumetric force - gravity level, on the magnetic fluid and thereby affect the system thermo-fluid behavior. The paper will describe the experimental and numerical modeling approach to the problem and discuss results obtained to date.

Description: Re-circulating flow in molten metal columns can be conveniently induced by the axisymmetric traveling magnetic field. A number of applications can benefit from this technique, such as mixing under microgravity environment, or crysta1 growth from metallic melts. For small magnetic field excitations, the flow is laminar and stationary. As the imposed field increases, a more complex flow will set up in the cylindrical column. Conditions for stable laminar flow are of importance for practical applications. In this work, a linear stability analysis is performed in order to determine the onset of the bifurcation in the system. Here the analysis is restricted to the axisymmetric modes and the low-frequency regime.

Description: Experiments to suppress thermocapillary oscillations using high-frequency vibrations were performed on float-zones. Such a float-zone is formed by melting one end of a vertically held sodium nitrate-barium nitrate crystal rod in contact with a hot surface at the top. In the experiments, when thermocapillary oscillation occurred, the bottom end of the rod was vibrated at a high frequency to generate fine ripples on the melt surface, driving a streaming flow in the opposite direction to that of the thermocapillary convection. It was observed that by generating a sufficiently strong streaming flow the thermocapillary flow can be offset enough such that the associated thermocapillarity oscillations can be quenched.

Description: Imagine a sequence of equidistant steps parallel to one another on a crystal face growing from solution. Steps may be one or part of a lattice spacing high. With such ideal morphology, the most perfect material is expected to be produced. Sometimes this happens. Much more often, however, steps form bunches, acquire wavy shapes and interlace. With time, these spontaneous perturbations may be increased in amplitude, i.e. the stepped surface becomes morphologically unstable. Interferometric studies of Y.G. Kuznetsov, L.N. Rashkovich, P.G. Vekilov and N.A. Booth and the author demonstrated that the step bunching depends at least on direction and rate of solution flow over the interface, supersaturation and presence of impurities. In particular, if solution and step flow are parallel to one another, the steps convene into bunches. If these flows are antiparallel, the bunches disappear. As the step train propagates, the overall bunch height often increases infinitely. However, we have found recently that if the flow rate is large enough, (approximately 1 meter per second), the bunch height may be limited. All these phenomena come from interaction between steps. The interaction mechanisms are still not fully understood. I plan to overview major approaches to the still unresolved problem on how these dissipative structures on growing crystal face appear and evolve.

Description: Turbulent fluxes of momentum, moisture, and heat at the air-sea interface are essential for climate studies. Version 2 Goddard Satellite-based Surface Turbulent Fluxes (GSSTF2) has been derived from the Special Sensor Microwave/Imager (SSM/I) radiance measurements. This dataset, covering the period July 1987-December 2000 over global oceans, has a spatial resolution of 1 deg x 1 deg lat-long and a temporal resolution of 1 day. Turbulent fluxes are derived from the SSM/I surface winds and surface air humidity, as well as the 2-m air and sea surface temperatures (SST) of the NCEP/NCAR reanalysis, using a bulk aerodynamic algorithm based on the surface layer similarity theory. The GSSTF2 bulk flux model, and retrieved daily wind stress, latent heat flux, wind speed, and surface air humidity validate well with ship observations of ten field experiments over the tropical and midlatitude oceans during 1991-99. The global distributions of 1988-2000 annual- and seasonal-mean turbulent fluxes show reasonable patterns related to the atmospheric general circulation and seasonal variations. Latent heat fluxes and related input parameters over global oceans during 1992-93 have been compared among GSSTF1 (version 1), GSSTF2, HOAPS (Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite Data), NCEP/NCAR reanalysis, and one based on COADS (Comprehensive Ocean-Atmosphere Data Set). Our analyses suggest that the GSSTF2 latent heat flux, surface air humidity, surface wind, and SST are quite realistic compared to the other four flux datasets examined. However, significant differences are found among these five flux datasets. The GSSTF2, available at http://daac.gsfc.nasa.gov/CAMPAIGN_DOCS/hydrology/hd_gsstf2.O.html, is useful for climate studies.

Description: Containerless processing is an important tool for materials research. The freedom from a crucible allows processing of liquid materials in a metastable undercooled state, as well as allowing processing of high temperature and highly reactive melts. Electrostatic levitation (ESL) is a containerless method which provides a number of unique advantages, including the ability to process non-conducting materials, the ability to operate in ultra-high vacuum or at moderate gas pressure (approx. = 5 atm), and the decoupling of positioning force from sample heating. ESL also has the potential to reduce internal flow velocities below those possible with electromagnetic, acoustic, or aero-acoustic techniques. In electrostatic levitation, the acceleration of gravity (or residual acceleration in reduced gravity) is opposed by the action of an applied electric field on a charged sample. Microgravity allows electrostatic levitation to work even more effectively. The ESL facility at NASA s Marshall Space Flight Center is in use for materials research and thermophysical property measurement by a number of different internal and external investigators. Results from the recent CDDF studies on the high energy X-ray beamline at the Advanced Photon Source of Argonne National Laboratory will be presented. The Microgravity Research Program supports the facility.

Description: We investigate the convective dynamics of evaporating pools of volatile liquids using an ultra-sensitive thermal imaging camera. During evaporation, there are significant convective flows inside the liquid due to Marangoni forces. We find that Marangoni convection during evaporation can dramatically affect the evaporation rates of volatile liquids. A simple heat balance model connects the convective velocities and temperature gradients to the evaporation rates.

Description: Low-dissipative high order filter finite difference methods for shock/turbulence/combustion compressible viscous MHD flows has been constructed. Several variants of the filter approach that cater to different flow types are proposed. These filters provide a natural and efficient way for the minimization of the divergence of the magnetic field (del (raided dot) B) numerical error in the sense that no standard divergence cleaning is required. For certain 2-D MHD test problems, divergence free preservation of the magnetic fields of these filter schemes has been achieved.

Description: Low-dissipative high order filter finite difference methods for long time wave propagation of shock/turbulence/combustion compressible viscous MHD flows has been constructed. Several variants of the filter approach that cater to different flow types are proposed. These filters provide a natural and efficient way for the minimization of the divergence of the magnetic field (Delta . B) numerical error in the sense that no standard divergence cleaning is required. For certain 2-D MHD test problems, divergence free preservation of the magnetic fields of these filter schemes has been achieved.

Description: The capabilities of the Dense Media Radiative Transfer model using the Quasi Crystalline Approximation with Coherent Potential (QCA-CP) to reproduce measured radiometric data were tested using the University of Tokyo Ground Based Microwave Radiometer (GBMR-7) during the third Intensive Observation Period (IOP3) of the NASA Cold-land Processes Field Experiment (CLPX). The data were collected at the Local-Scale Observation Site (LSOS), a 0.8-ha study site consisting of two open meadows separated by trees. Intensive measurements were also made of snow depth and temperature, density, and grain size profiles. A DMRT model is needed to describe radiative transfer in a medium such as snow because the assumption of independent scattering used in classical radiative transfer theory (CRT) is not valid. Validation of the DMRT approach requires a relationship between measured snow grain size and the DMRT approximation of snow grain radius as spherical particles with a mean radius of the log-normal particle-size distribution. This relationship is very important for a better understanding of snow modelling and for practical applications. DMRT simulations were compared with observations of microwave brightnesses at 18.7, 36.5 and 89-GHz (V and H polarizations) collected on February-1 9-25, 2003. Observation angles ranged from 30\deg to 70\deg. Model inputs included measured snow parameters except mean grain size. The average snow temperature, fractional volume and depth were held constant, together with the ice and soil permittivities. The minimum and maximum measured mean grain sizes were used to test the capabilities of the DMRT to reproduce the brightnesses as upper and lower limits. The sensitivity to the largest and smallest measured grain size in the three classes of minimum, medium and maximum observed grain sizes was also investigated. DMRT particle sizes yielding a best-fit to the experimental data for each date were computed. Results show that the measured brightnesses fall within the range of simulated brightnesses using the smallest and largest measured grain size values. The DMRT best-fit radii are comparable to the average radii for the medium observed grain sizes.

Description: A new, interferometer-based system for measuring thermal expansion to an absolute accuracy of 20 ppb or better at cryogenic temperatures has been developed. Data from NIST Copper SRM 736 measured from room temperature to 15 K will be presented along with data from many other materials including beryllium, ULE, Zerodur, and composite materials. Particular attention will be given to a study by the Space Optics Manufacturing Technology Center (SOMTC) investigating the variability of ULE and beryllium materials used in the AMSD program Approximately 20 samples of each material, tested from room temperature to below 30 K are compared as a function of billet location.

Description: Experiments on solutal convection in a paramagnetic fluid were conducted in a strong magnetic field gradient using a dilute solution of Manganese Chloride. The observed flows indicate that the magnetic field can completely counter the settling effects of gravity locally and are consistent with the theoretical predictions presented.

Description: One of the power systems under consideration for nuclear electric propulsion or as a planetary surface power source is a heatpipe-cooled reactor coupled to a Brayton cycle. In this system, power is transferred from the heatpipes to the Brayton gas via a heat exchanger attached to the heatpipes. This paper discusses the fluid, thermal and structural analyses that were performed in support of the design of the heat exchanger to be tested in the SAFE-100 experimental program at Marshall Space Flight Center. A companion paper, "Mechanical Design and Fabrication of a SAFE-100 Heat Exchanger for use in NASA s Advanced Propulsion Thermal-hydraulic Simulator", presents the fabrication issues and prototyping studies that, together with these analyses, led to the development of this heat exchanger. An important consideration throughout the design development of the heat exchanger was its capability to be utilized for higher power and temperature applications. This paper also discusses this aspect of the design and presents designs for specific applications that are under consideration.

Description: Re-circulating flow in molten metal columns can be conveniently induced by the axisymmetric traveling magnetic field. A number of applications can benefit from this technique, such as mixing under microgravity environment, or.crysta1 growth from metallic melts. For small magnetic field excitations, the flow is laminar and stationary. As the imposed field increases, a more complex flow will set up in the cylindrical column. Conditions for stable laminar flow are of importance for practical applications. In this work, a linear stability analysis is performed in order to determine the onset of the bifurcation in the system. Here the analysis is restricted to the axisymmetric modes and the low-frequency regime.

Description: Velocity fluctuations and particle concentrations are studied in a liquid fluidized bed to investigate the origin of steady state sedimentation. Both the velocity fluctuations and the particle concentrations are found to strongly depend on height. A flux balance model shows why the bed is stable: velocity fluctuations drive a downward particle flux that is compensated by an upward particle flux stemming from an excess flow velocity due to the stratification in concentration. Our results show that in steady state the magnitudes of the fluctuations are related to the degree of stratification.

Description: The study examined the effect of natural convection in applying the modified laser flash method to measure thermal properties of semiconductor melts. Common laser flash method uses a laser pulse to heat one side of a thin circular sample and measures the temperature response of the other side. Thermal diffusivity can be calculations based on a heat conduction analysis. For semiconductor melt, the sample is contained in a specially designed quartz cell with optical windows on both sides. When laser heats the vertical melt surface, the resulting natural convection can introduce errors in calculation based on heat conduction model alone. The effect of natural convection was studied by CFD simulations with experimental verification by temperature measurement. The CFD results indicated that natural convection would decrease the time needed for the rear side to reach its peak temperature, and also decrease the peak temperature slightly in our experimental configuration. Using the experimental data, the calculation using only heat conduction model resulted in a thermal diffusivity value is about 7.7% lower than that from the model with natural convection. Specific heat capacity was about the same, and the difference is within 1.6%, regardless of heat transfer models.

Description: The Lorentz force induced by a traveling magnetic field (TMF) in a cylindrical container has been calculated. The force can be used to control flow in dectrically conducting melts and the direction of the magnetic field and resulting flow can be reversed. A TMF can be used to partially cancel flow driven by buoyancy. The penetration of the field into the cylinder decreases as the frequency increases, and there exists an optimal value of frequency for which the resulting force is a maximum. Expressions for the Lorentz force in the limiting cases of low frequency and infinite cylinder are also given and compared to the numerical calculations.

Description: The ancillary components for Unitized Regenerative Fuel Cell (URFC) Energy Storage System are being developed at the NASA Glenn Research Center. This URFC system is unique in that it uses the surface area of the hydrogen and oxygen storage tanks as radiating heat surfaces for overall thermal control of the system. The waste heat generated by the URFC stack during charging and discharging is transferred from the cell stack to the surface of each tank by loop heat pipes. The heat pipes are coiled around each tank and covered with a thin layer of thermally conductive layer of carbon composite. The thin layer of carbon composite acts as a fin structure that spreads the heat away from the heat pipe and across the entire tank surface. Two different sized commercial grade composite tanks were constructed with integral heat pipes and tested in a thermal vacuum chamber to examine the feasibility of using the storage tanks as system radiators. The storage radiators were subjected to different steady-state heat loads and varying heat load profiles. The surface emissivity and specific heat capacity of each tank were calculated. The results were incorporated into a model that simulates the performance of similar radiators using lightweight, space rated carbon composite tanks.

Description: The movement of a particle in a fluid is one of the most fundamental processes in physics and plays an important role in materials science. For example, the physics of crystal growth from a nucleation site in solution can be dominated by movement of the solution over the crystal surface. Gravity usually dominates the equations of motion, but in microgravity other terms can dominate, making the equation much more complex. Until recently, the equations were solved only by numerical methods and/or by neglecting terms. During this study, we discovered an exact solution to the equations, which shows that the usually neglected terms become extremely important in microgravity. We also developed diagnostic recording methods using holography to save all of the particle field data, allowing the experiment to essentially be transferred from space back to earth in what we call the virtual apparatus . We will record holograms of particle distributions in motion in microgravity and bring them back to earth for analysis, allowing the study of the full three-dimensional motion of sets of particles, allowing us to test the new analytical solutions. The experiment will also provide accurate measurements of the quasi-steady acceleration of the space platform and other interesting microgravity and g-jitter effects. This program will produce the flight definition for an experiment in the microgravity environment of space to validate the theoretical model. We will design an experiment with the help of the theoretical model that is optimized for testing the model, measuring g, g-jitter, and other microgravity phenomena.

Description: Modeling flows with moving contact lines poses the formidable challenge that the usual assumptions of Newtonian fluid and no-slip condition give rise to a well-known singularity. This singularity prevents one from satisfying the contact angle condition to compute the shape of the fluid-fluid interface, a crucial calculation without which design parameters such as the pressure drop needed to move an immiscible 2-fluid system through a solid matrix cannot be evaluated. Some progress has been made for low Capillary number spreading flows. Combining experimental measurements of fluid-fluid interfaces very near the moving contact line with an analytical expression for the interface shape, we can determine a parameter that forms a boundary condition for the macroscopic interface shape when Ca much les than l. This parameter, which plays the role of an "apparent" or macroscopic dynamic contact angle, is shown by the theory to depend on the system geometry through the macroscopic length scale. This theoretically established dependence on geometry allows this parameter to be "transferable" from the geometry of the measurement to any other geometry involving the same material system. Unfortunately this prediction of the theory cannot be tested on Earth.

Description: Direct and large-eddy simulations of the interaction between the wake of a circular cylinder and a flat-plate boundary layer are conducted. Two Reynolds numbers are examined. The simulations indicate that at the lower Reynolds number the boundary layer is buffeted by the unsteady Karman vortex street shed by the cylinder. The fluctuations, however, cannot be self-sustained due to the low Reynolds-number, and the flow does not reach a turbulent state within the computational domain. In contrast, in the higher Reynolds-number case, boundary-layer fluctuations persist after the wake has decayed (due, in part, to the higher values of the local Reynolds number Re(sub theta) achieved in this case); some evidence could be observed that a self-sustaining turbulence generation cycle was beginning to be established.

Description: We analyze the velocity and temperature fields at steady state due to thermocapillary convection around a gas bubble that is stationary in a liquid. A linear temperature field is imposed in the undisturbed liquid. Our interest is in investigating the effect of convective transport of momentum and energy on the velocity and temperature fields. We assume the pertinent physical properties to be constant, and that buoyant convection is negligible. Suitably defined Reynolds and Marangoni numbers are assumed to be small compared with unity. When both the Reynolds and Marangoni numbers are set equal to zero, a solution can be found. In this solution, far from the bubble, the velocity field decays as the inverse of the distance from the bubble, and the disturbance temperature field decays as the inverse of the square of this distance. We now attempt to obtain a solution when the Reynolds number is zero, but the Marangoni number is small, but non-zero, by a perturbation expansion in the Marangoni number. When the temperature field is expanded in a regular perturbation series in the Marangoni number, we show that the problem for the first correction field is ill-posed. The governing equation for this perturbation field contains an inhomogeneity, and the corresponding particular solution neither decays far from the bubble, nor can be canceled by a homogeneous solution. Additional information is included in the original extended abstract.

Description: A new parabolized Navier-Stokes (PNS) algorithm has been developed to efficiently compute magnetohydrodynamic (MHD) flows in the low magnetic Reynolds number regime. In this regime, the electrical conductivity is low and the induced magnetic field is negligible compared to the applied magnetic field. The MHD effects are modeled by introducing source terms into the PNS equation which can then be solved in a very efficient manner. To account for upstream (elliptic) effects, the flowfields are computed using multiple streamwise sweeps with an iterated PNS algorithm. Turbulence has been included by modifying the Baldwin-Lomax turbulence model to account for MHD effects. The new algorithm has been used to compute both laminar and turbulent, supersonic, MHD flows over flat plates and supersonic viscous flows in a rectangular MHD accelerator. The present results are in excellent agreement with previous complete Navier-Stokes calculations.

Description: This report discusses the: 1. Development coupled electro/fluid/structural lumped-element model (LEM) of a prototypical flow-control actuator. 2. Validation the coupled electro/fluid/structural dynamics lumped-element models. 3. Development simple, yet effective, design tools for actuators. 4. Development structural dynamic models that accurately characterize the dynamic response of piezoelectric flap actuators using the Finite Element Method (FEW as well as analytical methods. 5. Perform a parametric study of a piezo-composite flap actuator. 6.Develop an optimization scheme for maximizing the actuator performance.

Description: The objective of this research is to improve and implement the filtered mass density function (FDF) methodology for large eddy simulation (LES) of high-speed reacting turbulent flows. We have just completed Year 1 of this research. This is the Final Report on our activities during the period: January 1, 2003 to December 31, 2003. 2002. In the efforts during the past year, LES is conducted of the Sandia Flame D, which is a turbulent piloted nonpremixed methane jet flame. The subgrid scale (SGS) closure is based on the scalar filtered mass density function (SFMDF) methodology. The SFMDF is basically the mass weighted probability density function (PDF) of the SGS scalar quantities. For this flame (which exhibits little local extinction), a simple flamelet model is used to relate the instantaneous composition to the mixture fraction. The modelled SFMDF transport equation is solved by a hybrid finite-difference/Monte Carlo scheme.

Description: We evaluate the applicability of a production computational fluid dynamics code for conducting detached eddy simulation for unsteady flows. A second-order accurate Navier-Stokes code developed at NASA Langley Research Center, known as TLNS3D, is used for these simulations. We focus our attention on high Reynolds number flow (Re = 5 x 10(sup 4) - 1.4 x 10(sup 5)) past a circular cylinder to simulate flows with large-scale separations. We consider two types of flow situations: one in which the flow at the separation point is laminar, and the other in which the flow is already turbulent when it detaches from the surface of the cylinder. Solutions are presented for two- and three-dimensional calculations using both the unsteady Reynolds-averaged Navier-Stokes paradigm and the detached eddy simulation treatment. All calculations use the standard Spalart-Allmaras turbulence model as the base model.

Description: This paper presents the results obtained from a recent ethane heat pipe program. Three identical ethane heat pipes were tested individually, and then two selected heat pipes were tested collectively in their system configuration. Heat transport, thermal conductance, and non-condensable gas tests were performed on each heat pipe. To gain insight into the reflux operation as seen at spacecraft level ground testing, the test fixture was oriented in a vertical configuration. The system level test included a computer-controlled heater designed to emulate the heat load generated at the thermoelectric cooler interface. The system performance was successfully characterized for a wide range of environmental conditions while staying within the operating limits.

Description: A silicon-based microhotplate tin oxide (SnO2) gas sensor integrated into a polymer-based microfluidic system for monitoring of contaminants in water systems is presented. This device is designed to sample a water source, control the sample vapor pressure within a microchannel using integrated resistive heaters, and direct the vapor past the integrated gas sensor for analysis. The sensor platform takes advantage of novel technology allowing direct integration of discrete silicon chips into a larger polymer microfluidic substrate, including seamless fluidic and electrical interconnects between the substrate and silicon chip.

Description: Flutter-free operation of advanced transonic fan designs continues to be a challenging task for the designers of aircraft engines. In order to meet the demands of increased performance and lighter weight, these modern fan designs usually feature low-aspect ratio shroudless rotor blade designs that make the task of achieving adequate flutter margin even more challenging for the aeroelastician. This is especially true for advanced forward swept designs that encompass an entirely new design space compared to previous experience. Fortunately, advances in unsteady computational fluid dynamic (CFD) techniques over the past decade now provide an analysis capability that can be used to quantitatively assess the aeroelastic characteristics of these next generation fans during the design cycle. For aeroelastic applications, Mississippi State University and NASA Glenn Research Center have developed the CFD code TURBO-AE. This code is a time-accurate three-dimensional Euler/Navier-Stokes unsteady flow solver developed for axial-flow turbomachinery that can model multiple blade rows undergoing harmonic oscillations with arbitrary interblade phase angles, i.e., nodal diameter patterns. Details of the code can be found in Chen et al. (1993, 1994), Bakhle et al. (1997, 1998), and Srivastava et al. (1999). To assess aeroelastic stability, the work-per-cycle from TURBO-AE is converted to the critical damping ratio since this value is more physically meaningful, with both the unsteady normal pressure and viscous shear forces included in the work-per-cycle calculation. If the total damping (aerodynamic plus mechanical) is negative, then the blade is unstable since it extracts energy from the flow field over the vibration cycle. TURBO-AE is an integral part of an aeroelastic design system being developed at Honeywell Engines, Systems & Services for flutter and forced response predictions, with test cases from development rig and engine tests being used to validate its predictive capability. A recent experimental program (Sanders et al., 2002) was aimed at providing the necessary unsteady aerodynamic and vibratory response data needed to validate TURBO-AE for fan flutter predictions. A comparison of numerical TURBO-AE simulations with the benchmark flutter data is given in Sanders et al. (2003), with the data used to guide the validation of the code and define best practices for performing accurate unsteady simulations. The agreement between the analyses and the predictions was quite remarkable, demonstrating the ability of the analysis to accurately model the unsteady flow processes driving stall-side flutter.

Description: The space-time conservation-element and solution-element method is employed to numerically study the near-field screech-tone noise of a typical underexpanded circular jet issuing from a sonic nozzle. Both axisymmetric and fully three-dimensional computations are carried out. The self-sustained feedback loop is properly simulated. The computed shock-cell structure, acoustic wave length, screech-tone frequency, and sound-pressure levels are in good agreement with existing experimental results.

Description: The flow fields of synthetic jets in a cross-flow from orifices of different geometry are investigated. The geometries include a straight, a tapered, a pitched and a cluster of nine orifices, all having the same cross-sectional area through which the perturbation is discharged into the cross-flow. The strength of the jet from the tapered orifice in comparison to that from the straight one is found to be only slightly enhanced. The flow field from the cluster of orifices, when viewed a few equivalent diameters downstream, is similar to that from the single orifice. However, the penetration is somewhat lower in the former case due to the increased mixing of the distributed jets with the cross-flow. The penetration for the pitched configuration is the lowest, as expected. The jet trajectories for the straight and pitched orifices are well represented by correlation equations available for steady jets-in-cross-flow. Distributions of streamwise velocity, vorticity as well as turbulence intensity are documented for various cases. In addition, distributions of phase-averaged velocity and vorticity for the cylindrical and the clustered orifices are presented providing an insight into the flow dynamics.

Description: Contents include the following: 1. High frequency, low amplitude temperature oscillations: LHP operation - governing equations; interactions among LHP components; factors affecting low amplitude temperature oscillations. 2. Test results. 3. Conclusions.

Description: This paper presents viewgraphs on the low frequency high amplitude temperature oscillations observed in loop heat pipe operations. The topics include: 1) Proposed Theory; 2) Test Loop and Test Results; and 3) Effects of Various Parameters. The author also presents a short summary on the conditiions that must be met in order to sustain a low frequency high amplitude temperature oscillation.

Description: This work examined the wake aerodynamics of a single helicopter rotor blade with several tip shapes operating on a hover test stand. Velocity field measurements were conducted using three-component laser Doppler velocimetry (LDV). The objective of these measurements was to document the vortex velocity profiles and then extract the core properties, such as the core radius, peak swirl velocity, and axial velocity. The measured test cases covered a wide range of wake-ages and several tip shapes, including rectangular, tapered, swept, and a subwing tip. One of the primary differences shown by the change in tip shape was the wake geometry. The effect of blade taper reduced the initial peak swirl velocity by a significant fraction. It appears that this is accomplished by decreasing the vortex strength for a given blade loading. The subwing measurements showed that the interaction and merging of the subwing and primary vortices created a less coherent vortical structure. A source of vortex core instability is shown to be the ratio of the peak swirl velocity to the axial velocity deficit. The results show that if there is a turbulence producing region of the vortex structure, it will be outside of the core boundary. The LDV measurements were supported by laser light-sheet flow visualization. The results provide several benchmark test cases for future validation of theoretical vortex models, numerical free-wake models, and computational fluid dynamics results.

Description: A finite volume based network analysis procedure has been applied to model unsteady flow without and with heat transfer. Liquid has been modeled as compressible fluid where the compressibility factor is computed from the equation of state for a real fluid. The modeling approach recognizes that the pressure oscillation is linked with the variation of the compressibility factor; therefore, the speed of sound does not explicitly appear in the governing equations. The numerical results of chilldown process also suggest that the flow and heat transfer are strongly coupled. This is evident by observing that the mass flow rate during 90-second chilldown process increases by factor of ten.

Description: Numerical simulations have been completed for a variety of designs for a 90 deg elbow duct. The objective is to identify a design that minimizes the dynamic load entering a LOX turbopump located at the elbow exit. Designs simulated to date indicate that simpler duct geometries result in lower losses. Benchmark simulations have verified that the compressible flow codes used in this study are applicable to these incompressible flow simulations.

Description: Reynolds-averaged Navier-Stokes calculations are used to investigate porous side-edge treatment as a passive means for flap noise reduction. Steady-state simulations are used to infer effects of the treatment on acoustically relevant features of the mean flow near the flap side edge. Application of the porous treatment over a miniscule fraction of the wetted flap area (scaling with the flap thickness) results in significantly weaker side-edge vortex structures via modification of the vortex initiation and roll-up processes. At high flap deflections, the region of axial flow reversal associated with the breakdown of the side-edge vortex is also eliminated, indicating an absence of vortex bursting in the presence of the treatment. Potential ramifications of the mean-flow modifications for flap-noise reduction are examined in the light of lessons learned from recent studies on flap noise. Computations confirm that any noise reduction benefit via the porous treatment would be achieved without compromising the aerodynamic effectiveness of the flap. Results of the parameter study contribute additional insight into the measured data from the 7x10 wind tunnel at NASA Ames and provide preliminary guidance for specifying optimal treatment characteristics in terms of treatment location, spatial extent, and flow resistance of the porous skin.

Description: As part of the aero-thermodynamics team supporting the Columbia Accident Investigation Board (CAB), the Marshall Space Flight Center was asked to perform engineering analyses of internal flows in the port wing. The aero-thermodynamics team was split into internal flow and external flow teams with the support being divided between shorter timeframe engineering methods and more complex computational fluid dynamics. In order to gain a rough order of magnitude type of knowledge of the internal flow in the port wing for various breach locations and sizes (as theorized by the CAB to have caused the Columbia re-entry failure), a bulk venting model was required to input boundary flow rates and pressures to the computational fluid dynamics (CFD) analyses. This paper summarizes the modeling that was done by MSFC in Thermal Desktop. A venting model of the entire Orbiter was constructed in FloCAD based on Rockwell International s flight substantiation analyses and the STS-107 reentry trajectory. Chemical equilibrium air thermodynamic properties were generated for SINDA/FLUINT s fluid property routines from a code provided by Langley Research Center. In parallel, a simplified thermal mathematical model of the port wing, including the Thermal Protection System (TPS), was based on more detailed Shuttle re-entry modeling previously done by the Dryden Flight Research Center. Once the venting model was coupled with the thermal model of the wing structure with chemical equilibrium air properties, various breach scenarios were assessed in support of the aero-thermodynamics team. The construction of the coupled model and results are presented herein.

Description: This paper presents the Turbine Air-Flow Test (TAFT) rig computational fluid dynamics (CFD) results for test matrix. The topics include: 1) TAFT Background; 2) Design Point CFD; 3) TAFT Test Plan and Test Matrix; and 4) CFD of Test Points. This paper is in viewgraph form.

Description: This report consists of two published papers, 'Computation of Magnetohydrodynamic Flows Using an Iterative PNS Algorithm' and 'Numerical Simulation of Turbulent MHD Flows Using an Iterative PNS Algorithm'.

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

Description: During the past thirty years, numerical methods and simulation tools for incompressible flows have been advanced as a subset of CFD discipline. Even though incompressible flows are encountered in many areas of engineering, simulation of compressible flow has been the major driver for developing computational algorithms and tools. This is probably due to rather stringent requirements for predicting aerodynamic performance characteristics of flight vehicles, while flow devices involving low speed or incompressible flow could be reasonably well designed without resorting to accurate numerical simulations. As flow devices are required to be more sophisticated and highly efficient, CFD tools become indispensable in fluid engineering for incompressible and low speed flow. This paper is intended to review some of the successes made possible by advances in computational technologies during the same period, and discuss some of the current challenges.

Description: A generalized multi-phase formulation for cavitation in fluids operating at temperatures elevated relative to their critical temperatures is presented. The thermal effects and the accompanying property variations due to phase change are modeled rigorously. Thermal equilibrium is assumed and fluid thermodynamic properties are specified along the saturation line using the NIST-12 databank. Fundamental changes in the physical characteristics of the cavity when thermal effects become pronounced are identified; the cavity becomes more porous, the interface less distinct, and has increased entrainment when temperature variations are present. Quantitative estimates of temperature and pressure depressions in both liquid nitrogen and liquid hydrogen were computed and compared with experimental data of Hord for hydrofoils. Excellent estimates of the leading edge temperature and pressure depression were obtained while the comparisons in the cavity closure region were reasonable. Liquid nitrogen cavities were consistently found to be in thermal equilibrium while liquid hydrogen cavities exhibited small, but distinct, non-equilibrium effects.

Description: A new software tool, AeroDB, is used to compute thousands of Euler and Navier-Stokes solutions for a 2nd generation glide-back booster in one week. The solution process exploits a common job-submission grid environment using 13 computers located at 4 different geographical sites. Process automation and web-based access to the database greatly reduces the user workload, removing much of the tedium and tendency for user input errors. The database consists of forces, moments, and solution files obtained by varying the Mach number, angle of attack, and sideslip angle. The forces and moments compare well with experimental data. Stability derivatives are also computed using a monotone cubic spline procedure. Flow visualization and three-dimensional surface plots are used to interpret and characterize the nature of computed flow fields.

Description: In this study 3-D Computational Fluid Dynamic (CFD) runs have been made for the Space Shuttle RSRM using 2 different grids and 4 different turbulent models, which were the Standard KE, the RNG KE, the Realizable KE, and the Reynolds stress model. The RSRM forward segment consists of 11 fins. By taking advantage of the forward fin symmetry only half of one fin along the axis had to be used in making the grid. This meant that the 3-D model consisted of a pie slice that encompassed 1/22nd of the motor circumference and went along the axis of the entire motor. The 3-D flow patterns in the forward fin region are of particular interest. Close inspection of these flow patterns indicate that 2 counter-rotating axial vortices emerge from each submerged solid propellant fin. Thus, the 3-D CFD analysis allows insight into complicated internal motor flow patterns that are not available from the simpler 2-D axi-symmetric studies. In addition, a comparison is made between the 3-D bore pressure drop and the 2-D axi-symmetric pressure drop.

Description: A four-parameter, three-level, central composite experiment design has been used to optimize the configuration of an annular scramjet injector geometry using computational fluid dynamics. The computational fluid dynamic solutions played the role of computer experiments, and response surface methodology was used to capture the simulation results for mixing efficiency and total pressure recovery within the scramjet flowpath. An optimization procedure, based upon the response surface results of mixing efficiency, was used to compare the optimal design configuration against the target efficiency value of 92.5%. The results of three different optimization procedures are presented and all point to the need to look outside the current design space for different injector geometries that can meet or exceed the stated mixing efficiency target.

Description: A concerted effort is underway at NASA Langley Research Center to create a benchmark for Computational Fluid Dynamic (CFD) codes. both unstructured and structured, against a data set for the hump model with actuation. The hump model was tested in the NASA Langley 0.3-m Transonic Cryogenic Tunnel. The CFD codes used for the analyses are the FUN2D (Full Unstructured Navier-Stokes 2-Dimensional) code, the structured TLNS3D (Thin-Layer Navier-Stokes 3-Dimensional) code, and the structured CFL3D code, all developed at NASA Langley. The current investigation uses the time-accurate Reynolds-Averaged Navier-Stokes (RANS) approach to predict aerodynamic performance of the active flow control experimental database for the hump model. Two-dimensional computational results verified that steady blowing and suction and oscillatory suction/blowing can be used to significantly reduce the separated flow region on the model. Discrepancies do exist between the CFD results and experimental data in the region downstream of the slot with the largest differences in the oscillatory cases. Overall, the structured CFD codes exhibited similar behavior with each other for a wide range of control conditions, with the unstructured FUN2D code showing moderately different results in the separated flow region for the suction and oscillatory cases.

Description: A new high Reynolds number test capability for boundary layer ingesting inlets has been developed for the NASA Langley Research Center 0.3-Meter Transonic Cryogenic Tunnel. Using this new capability, an experimental investigation of four S-duct inlet configurations with large amounts of boundary layer ingestion (nominal boundary layer thickness of about 40% of inlet height) was conducted at realistic operating conditions (high subsonic Mach numbers and full-scale Reynolds numbers). The objectives of this investigation were to 1) provide a database for CFD tool validation on boundary layer ingesting inlets operating at realistic conditions and 2) provide a baseline inlet for future inlet flow-control studies. Tests were conducted at Mach numbers from 0.25 to 0.83, Reynolds numbers (based on duct exit diameter) from 5.1 million to a full-scale value of 13.9 million, and inlet mass-flow ratios from 0.39 to 1.58 depending on Mach number. Results of this investigation indicate that inlet pressure recovery generally decreased and inlet distortion generally increased with increasing Mach number. Except at low Mach numbers, increasing inlet mass-flow increased pressure recovery and increased distortion. Increasing the amount of boundary layer ingestion (by decreasing inlet throat height) or ingesting a boundary layer with a distorted (adverse) profile decreased pressure recovery and increased distortion. Finally, increasing Reynolds number had almost no effect on inlet distortion but increased inlet recovery by about one-half percent at a Mach number near cruise.

Description: The initial efforts to develop the capability to model valves used in rocket engine component testing at Stennis Space Center are documented. An axisymmetric model of a control valve with LN2 as the working fluid was developed. The goal was to predict the effect of change in the plug/sear region of the valve prior to testing. The valve flow coefficient was predicted for a range of plug positions. Verification of the calculations was carried out to quantify the uncertainty in the numerical answer. The modeled results compared well qualitatively to experimental trends. Additionally, insights into the flow processes in the valve were obtained. Benefits from the verification process included the ability to use coarser grids and insight into ways to reduce computational time by using double precision accuracy and non-integer grid ratios. Future valve modeling activities will include shape optimization of the valve/seat region and dynamic grid modeling.

Description: The YO1.5-ZrO2 system consists of five solid solutions, one liquid solution, and one intermediate compound. A thermodynamic description of this system is developed, which allows calculation of the phase diagram and thermodynamic properties. Two different solution models are used-a neutral species model with YO1.5 and ZrO2 as the components and a charged species model with Y(+3), Zr(+4), O(-2), and vacancies as components. For each model, regular and sub-regular solution parameters are derived fiom selected equilibrium phase and thermodynamic data.

Description: In this work, a new type of boundary condition for time-accurate Computational Aeroacoustics solvers is described. This boundary condition is designed to complement the existing nonreflective boundary conditions while ensuring that the correct mean flow conditions are maintained throughout the flow calculation. Results are shown for a loaded 2D cascade, started with various initial conditions.

Description: A three dimensional linear plug nozzle of area ratio 12.79 was designed by EADS Space Transportation (former Astrium Space Infrastructure). The nozzle was tested within the German National Technology Program 'LION' in a cold air wind tunnel by TU Dresden. The experimental hardware and test conditions are described. Experimental data was obtained for the nozzle without plug side wall fences at a nozzle pressure ratio of 116 and then with plug side wall fences at NPR 110. Schlieren images were recorded and axial profiles of plug wall static pressures were measured at several spanwise locations and on the plug base. Detailed CFD analysis was performed for these nozzle configurations at NPR 116 by NASA MSFC. The CFD exhibits good agreement with the experimental data. A detailed comparison of the CFD results and the experimental plug wall pressure data are given. Comparisons are made for both the without and with plug side wall fence configurations. Numerical results for density gradient are compared to experimental Schlieren images. Experimental nozzle thrust efficiencies are calculated based on the CFD results. The CFD results are used to illustrate the plug nozzle fluid dynamics. The effect of the plug side wall is emphasized.

Description: Pulse detonation rocket engines (PDREs) offer potential performance improvements over conventional designs, but represent a challenging modeling task. A quasi 1-D, finite-rate chemistry CFD model for a PDRE & described and implemented. A parametric study of the effect of blowdown pressure ratio on the performance of an optimized, fixed PDRE nozzle configuration is reported. The results are compared to a steady-state rocket system using similar modeling assumptions.

Description: A 2D parallel Euler code based on the space-time conservation element and solution element (CE/SE) method is validated by solving the benchmark problem I in Category 3 of the Third CAA Workshop. This problem concerns the acoustic field generated by the interaction of a convected harmonic vortical gust with a single airfoil. Three gust frequencies, two gust configurations, and three airfoil geometries are considered. Numerical results at both near and far fields are presented and compared with the analytical solutions, a frequency-domain solver GUST3D solutions, and a time-domain high-order Discontinuous Spectral Element Method (DSEM) solutions. It is shown that the CE/SE solutions agree well with the GUST3D solution for the lowest frequency, while there are discrepancies between CE/SE and GUST3D solutions for higher frequencies. However, the CE/SE solution is in good agreement with the DSEM solution for these higher frequencies. It demonstrates that the CE/SE method can produce accurate results of CAA problems involving complex geometries by using unstructured meshes.

Description: This viewgraph presentation provides an overview of future research issues and objectives of the Marshall Space Flight Center's Turbomachinery Fluid Dynamics team. Research focuses on developing and enhancing aerodynamic and hydrodynamic design tools and computational fluid dynamics (CFD) analysis tools, applying these tools to development of hardware technology and concepts, and defining cold flow experiments to provide tool validation, technology and concept verification, and to further expand fluid dynamics knowledge and understanding.

Description: Twin thruster tests have been conducted with the Penn State RBCC test article operating at sea- level static conditions. Significant differences were observed in the performance characteristics for two different thruster centerline spacings. Changing the thruster spacing from 2.50 to 1.75 in. reduced the entrained air velocity (-17%) and the thrust (-7%) for tests at a thruster chamber pressure of 200 psia and MR = 8. In addition, significant differences were seen in the static pressure profiles, the Raman spectroscopy profiles, and the acoustic power spectrum for these two configurations.

Description: To develop future Reusable Launch Vehicle concepts, we have conducted design optimization for a single element rocket injector, with overall goals of improving reliability and performance while reducing cost. Computational solutions based on the Navier-Stokes equations, finite rate chemistry, and the k-E turbulence closure are generated with design of experiment techniques, and the response surface method is employed as the optimization tool. The design considerations are guided by four design objectives motivated by the consideration in both performance and life, namely, the maximum temperature on the oxidizer post tip, the maximum temperature on the injector face, the adiabatic wall temperature, and the length of the combustion zone. Four design variables are selected, namely, H2 flow angle, H2 and O2 flow areas with fixed flow rates, and O2 post tip thickness. In addition to establishing optimum designs by varying emphasis on the individual objectives, better insight into the interplay between design variables and their impact on the design objectives is gained. The investigation indicates that improvement in performance or life comes at the cost of the other. Best compromise is obtained when improvements in both performance and life are given equal importance.