ALBERT

Description: A computational investigation was performed to assess the effectiveness of Detached Eddy Simulation (DES) as a tool for predicting icing effects. The AVUS code, a public domain flow solver, was employed to compute solutions for an iced wing configuration using DES and steady Reynolds Averaged Navier-Stokes (RANS) equation methodologies. The wing section considered here was a business jet airfoil (GLC305) with a 22.5-minute glaze ice accretion (944-ice shape). The section was extruded to form a rectangular planform. The model was mounted between two walls so no tip effects were considered. The numerical results were validated by comparison with experimental data for the same configuration. The time-averaged DES computations showed some improvement in lift and drag results near stall when compared to steady RANS results. However, comparisons of the flow field details did not show the level of agreement suggested by the integrated quantities. More validation is needed to determine what role DES can play as part of an overall icing effects prediction strategy.

Description: We present a Brook streaming language implementation of the 3-D discontinuous Galerkin method for compressible fluid flow on tetrahedral meshes. Efficient implementation of the discontinuous Galerkin method using the streaming model of computation introduces several algorithmic design challenges. Using a cycle-accurate simulator, performance characteristics have been obtained for the Stanford Merrimac stream processor. The current Merrimac design achieves 128 Gflops per chip and the desktop board is populated with 16 chips yielding a peak performance of 2 Teraflops. Total parts cost for the desktop board is less than $20K. Current cycle-accurate simulations for discretizations of the 3-D compressible flow equations yield approximately 40-50% of the peak performance of the Merrimac streaming processor chip. Ongoing work includes the assessment of the performance of the same algorithm on the 2 Teraflop desktop board with a target goal of achieving 1 Teraflop performance.

Description: A series of experiments have been performed to investigate the oscillatory flow pressure and pressure drop through complex geometries. These experiments were conducted at the CSU-SLRE facility which is a horizontally opposed, two-piston, single-acting engine with a split crankshaft driving mechanism. Flow through a rectangular duct, with no insert (obstruction), was studied first. Then four different inserts were examined: Abrupt, Manifold, Diverging Short and Diverging Long. The inserts were mounted in the center of the rectangular duct to represent different type of geometries that could be encountered in Stirling machines. The pressure and pressure drop of the oscillating flow was studied for: 1) different inserts, 2) different phase angle between the two pistons of the engine (zero, 90 lead, 180, and 90 lag), and 3) for different piston frequencies (5, 10, 15, and 20 Hz). It was found that the pressure drop of the oscillatory flow increases with increasing Reynolds number. The pressure drop was shown to be mainly due to the gas inertia for the case of oscillatory flow through a rectangular duct with no insert. On the other hand, for the cases with different inserts into the rectangular duct, the pressure drop has three sources: inertia, friction, and local losses. The friction pressure drop is only a small fraction of the total pressure drop. It was also shown that the dimensionless pressure drop decreases with increasing kinetic Reynolds number.

Description: This is the final report on our NASA grant, Problems in Microgravity Fluid Mechanics NAG3-2513: 12/14/2000 - 11/30/2003, extended through 11/30/2004. This grant was made to Stanford University and then transferred to the University of California at Santa Barbara when the PI relocated there in January 2001. Our main activity has been to conduct both experimental and theoretical studies of instabilities in fluids that are relevant to the microgravity environment, i.e. those that do not involve the action of buoyancy due to a steady gravitational field. Full details of the work accomplished under this grant are given below. Our work has focused on: (i) Theoretical and computational studies of the effect of g-jitter on instabilities of convective states where the convection is driven by forces other than buoyancy (ii) Experimental studies of instabilities during displacements of miscible fluid pairs in tubes, with a focus on the degree to which these mimic those found in immiscible fluids. (iii) Theoretical and experimental studies of the effect of time dependent electrohydrodynamic forces on chaotic advection in drops immersed in a second dielectric liquid. Our objectives are to acquire insight and understanding into microgravity fluid mechanics problems that bear on either fundamental issues or applications in fluid physics. We are interested in the response of fluids to either a fluctuating acceleration environment or to forces other than gravity that cause fluid mixing and convection. We have been active in several general areas.

Description: This paper examines the behavior of flux and slope limiters on non-uniform grids in multiple dimensions. Many slope limiters in standard use do not preserve linear solutions on irregular grids impacting both accuracy and convergence. We rewrite some well-known limiters to highlight their underlying symmetry, and use this form to examine the proper - ties of both traditional and novel limiter formulations on non-uniform meshes. A consistent method of handling stretched meshes is developed which is both linearity preserving for arbitrary mesh stretchings and reduces to common limiters on uniform meshes. In multiple dimensions we analyze the monotonicity region of the gradient vector and show that the multidimensional limiting problem may be cast as the solution of a linear programming problem. For some special cases we present a new directional limiting formulation that preserves linear solutions in multiple dimensions on irregular grids. Computational results using model problems and complex three-dimensional examples are presented, demonstrating accuracy, monotonicity and robustness.

Description: The Lag model has shown great promise in prediction of low speed and transonic separations. The predictions of the model, along with other models (Spalart-Allmaras and Menter SST) are assessed for various high speed flowfields. In addition to skin friction and separation predictions, the prediction of heat transfer are compared among these models, and some fundamental building block flowfields, are investigated.

Description: A joint computational and experimental study has been performed at NASA Langley Research Center to investigate the unsteady flow generated by the components of an aircraft landing gear system. Because the flow field surrounding a full landing gear is so complex, the study was conducted on a simplified geometry consisting of two cylinders in tandem arrangement to isolate and characterize the pertinent flow phenomena. This paper focuses on the experimental effort where surface pressures, 2-D Particle Image Velocimetry, and hot-wire anemometry were used to document the flow interaction around the two cylinders at a Reynolds Number of 1.66 x 10(exp 5), based on cylinder diameter, and cylinder spacing-todiameter ratios, L/D, of 1.435 and 3.70. Transition strips were applied to the forward cylinder to produce a turbulent boundary layer upstream of the flow separation. For these flow conditions and L/D ratios, surface pressures on both the forward and rear cylinders show the effects of L/D on flow symmetry, base pressure, and the location of flow separation and attachment. Mean velocities and instantaneous vorticity obtained from the PIV data are used to examine the flow structure between and aft of the cylinders. Shedding frequencies and spectra obtained using hot-wire anemometry are presented. These results are compared with unsteady, Reynolds-Averaged Navier-Stokes (URANS) computations for the same configuration in a companion paper by Khorrami, Choudhari, Jenkins, and McGinley (2005). The experimental dataset produced in this study provides information to better understand the mechanisms associated with component interaction noise, develop and validate time-accurate computer methods used to calculate the unsteady flow field, and assist in modeling of the radiated noise from landing gears.

Description: Turbulence has two distinctive effects on sonic booms: there is distortion in the form of random perturbations that appear behind the shock waves, and shock rise times are increased randomly. A first scattering theory by S.C. Crow in the late 1960s quantified the random distortions, and Crow's theory was shown to agree with available flight test data. A variety of theories for the shock thickness have been presented, all supporting the role of turbulence in increasing rise time above that of a basic molecular-relaxation structure. The net effect of these phenomena on the loudness of shaped minimized booms is of significant interest. Initial analysis suggests that there would be no change to average loudness, but this had not been experimentally investigated. The January 2004 flight test of the Shaped Sonic Boom Demonstrator (SSBD), together with a reference unmodified F-5E, included a 12500- foot linear ground sensor array with 28 digitally recorded sensor sites. This data set provides an opportunity to re-test Crow's theory for the post-shock perturbations, and to examine the net effect of turbulence on the loudness of shaped sonic booms.

Description: A 2-D impedance eduction methodology is extended to quasi-3-D sound fields in uniform or shearing mean flow. We introduce a nonlocal, nonreflecting boundary condition to terminate the duct and then educe the impedance by minimizing an objective function. The introduction of a parallel, sparse, equation solver significantly reduces the wall clock time for educing the impedance when compared to that of the sequential band solver used in the 2-D methodology. The accuracy, efficiency, and robustness of the methodology is demonstrated using two examples. In the first example, we show that the method reproduces the known impedance of a ceramic tubular test liner. In the second example, we illustrate that the approach educes the impedance of a four-segment liner where the first, second, and fourth segments consist of a perforated face sheet bonded to honeycomb, and the third segment is a cut from the ceramic tubular test liner. The ability of the method to educe the impedances of multisegmented liners has the potential to significantly reduce the amount of time and cost required to determine the impedance of several uniform liners by allowing them to be placed in series in the test section and to educe the impedance of each segment using a single numerical experiment. Finally, we probe the objective function in great detail and show that it contains a single minimum. Thus, our objective function is ideal for use with local, inexpensive, gradient-based optimizers.

Description: This viewgraph presentation discusses spray cooling as it applies to thermal management of the new ESMD (Exploration Systems Mission Directorate) programs. It is concluded that given the targeted heat flux and temperature operating regimes, spray cooling is a viable option, although further technology development will be required.

Description: We are working with Project Scientists R. Balasubramanian and Sang Young Son, and a NASA Projects Team headed by Sue Motil at the Glenn Research Center on the design and development of an experimental system for use on the International Space Station during the year 2006. John Eustace is the coordinator for the flight experiment at Zin-Tech (previously Northrop-Grumman) for the design and development of the Constrained Vapor Bubble Heat Exchanger, CVBHX, cell which will fit into the Light Microscope Module, LMM. Good progress is being made. The CDR for the LMM being developed was held on December 10-1 1,2003. Experimental results obtained under microgravity conditions will be compared with those obtained at Rensselaer. Basic and applied research at Rensselaer continues on the experimental and theoretical details associated with passive phase change heat transfer processes controlled by interfacial forces in the CVBHX. The extensive results of our current research are presented in the 23 external publications listed below. Twenty-two external presentations have been given. Briefly, evaporation/condensation data from both vertical and horizontal CVBHX systems were obtained and analyzed for both polar (wetting) and apolar (partially wetting) fluids. The vertical system is axi-symmetric, but strongly effected by gravity. Whereas, the horizontal system is asymmetric, but weakly effected by gravity. Therefore, there will be significant differences in the operation of the cell in the earth s environment versus the operation under microgravity conditions. Due to its relative large size, the system s performance should be optimum under micro-gravity conditions, where the CVBHX should be a very effective passive heat exchanger. The CVBHX was found to be an ideal experimental setup in which to study the effects of interfacial phenomena on both the evaporation and drop-wise condensation processes. The optical technique (Image Analyzing Interferometry, IAI), which is based on the measurement and analysis of the reflectivity pattern of a thin film, was significantly improved. The accuracy of the IAI system is of critical importance to the success of the mission because it is used to measure the details of the pressure field in the liquid by measuring the film thickness profile. The accuracy was found to be excellent and various publications/presentations documenting these new results were written. Significant new results were also obtained for the effect of the oscillating contact line region on evaporation. Three doctoral students graduated under this grant. All three work in US industry, two for Intel. Another doctoral student is in his third year of study and will finish under an extension of the NASA grant: # NNC05GA27G.

Description: The aeroelastic and aerothermoelastic behavior of three-dimensional configurations in hypersonic flow regime are studied. The aeroelastic behavior of a low aspect ratio wing, representative of a fin or control surface on a generic hypersonic vehicle, is examined using third order piston theory, Euler and Navier-Stokes aerodynamics. The sensitivity of the aeroelastic behavior generated using Euler and Navier-Stokes aerodynamics to parameters governing temporal accuracy is also examined. Also, a refined aerothermoelastic model, which incorporates the heat transfer between the fluid and structure using CFD generated aerodynamic heating, is used to examine the aerothermoelastic behavior of the low aspect ratio wing in the hypersonic regime. Finally, the hypersonic aeroelastic behavior of a generic hypersonic vehicle with a lifting-body type fuselage and canted fins is studied using piston theory and Euler aerodynamics for the range of 2.5 less than or equal to M less than or equal to 28, at altitudes ranging from 10,000 feet to 80,000 feet. This analysis includes a study on optimal mesh selection for use with Euler aerodynamics. In addition to the aeroelastic and aerothermoelastic results presented, three time domain flutter identification techniques are compared, namely the moving block approach, the least squares curve fitting method, and a system identification technique using an Auto-Regressive model of the aeroelastic system. In general, the three methods agree well. The system identification technique, however, provided quick damping and frequency estimations with minimal response record length, and therefore o ers significant reductions in computational cost. In the present case, the computational cost was reduced by 75%. The aeroelastic and aerothermoelastic results presented illustrate the applicability of the CFL3D code for the hypersonic flight regime.

Description: Recent advancement in the molecular Rayleigh scattering based technique allowed for simultaneous measurement of velocity and density fluctuations with high sampling rates. The technique was used to investigate unheated high subsonic and supersonic fully expanded free jets in the Mach number range of 0.8 to 1.8. The difference between the Favre averaged and Reynolds averaged axial velocity and axial component of the turbulent kinetic energy is found to be small. Estimates based on the Morkovin's "Strong Reynolds Analogy" were found to provide lower values of turbulent density fluctuations than the measured data.

Description: The aerodynamic performance of an isolated fan or rotor alone model was measured in the NASA Glenn Research Center 9- by 15- Foot Low Speed Wind Tunnel as part of the Fan Broadband Source Diagnostic Test conducted at NASA Glenn. The Source Diagnostic Test was conducted to identify the noise sources within a wind tunnel scale model of a turbofan engine and quantify their contribution to the overall system noise level. The fan was part of a 1/5th scale model representation of the bypass stage of a current technology turbofan engine. For the rotor alone testing, the fan and nacelle, including the inlet, external cowl, and fixed area fan exit nozzle, were modeled in the test hardware; the internal outlet guide vanes located behind the fan were removed. Without the outlet guide vanes, the velocity at the nozzle exit changes significantly, thereby affecting the fan performance. As part of the investigation, variations in the fan nozzle area were tested in order to match as closely as possible the rotor alone performance with the fan performance obtained with the outlet guide vanes installed. The fan operating performance was determined using fixed pressure/temperature combination rakes and the corrected weight flow. The performance results indicate that a suitable nozzle exit was achieved to be able to closely match the rotor alone and fan/outlet guide vane configuration performance on the sea level operating line. A small shift in the slope of the sea level operating line was measured, which resulted in a slightly higher rotor alone fan pressure ratio at take-off conditions, matched fan performance at cutback conditions, and a slightly lower rotor alone fan pressure ratio at approach conditions. However, the small differences in fan performance at all fan conditions were considered too small to affect the fan acoustic performance.

Description: This paper describes the flow within a prototype actuator, energized by pulsed combustion or detonations, that provides a pulsed jet suitable for flow control in high-speed applications. A high-speed valve, capable of delivering a pulsed stream of reactants a mixture of H2 and air at rates of up to 1500 pulses per second, has been constructed. The reactants burn in a resonant chamber, and the products exit the device as a pulsed jet. High frequency pressure transducers have been used to monitor the pressure fluctuations in the device at various reactant injection frequencies, including both resonant and off-resonant conditions. The combustion chamber has been constructed with windows, and the flow inside it has been visualized using Planar Laser-Induced Fluorescence (PLIF). The pulsed jet at the exit of the device has been observed using schlieren.

Description: The goal of the DARPA Shaped Sonic Boom Demonstration (SSBD) Program was to demonstrate for the first time in flight that sonic booms can be substantially reduced by incorporating specialized aircraft shaping techniques. Although mitigation of the sonic boom via specialized shaping techniques was theorized decades ago, until now, this theory had never been tested with a flight vehicle subjected to actual flight conditions in a real atmosphere. The demonstrative success, which occurred on 27 August 2003 with repeat flights in the supersonic corridor at Edwards Air Force Base, is a critical milestone in the development of next generation supersonic aircraft that could one day fly unrestricted over land and help usher in a new era of time-critical air transport. Pressure measurements obtained on the ground and in the air confirmed that the specific modifications made to a Northrop Grumman F-5E aircraft not only changed the shape of the shock wave signature emanating from the aircraft, but also produced a flat-top signature whose shape persisted, as predicted, as the pressure waves propagated through the atmosphere to the ground. This accomplishment represents a major advance towards reducing the startling and potentially damaging noise of a sonic boom. This paper describes the evolution of the SSBD program, including the rationale for test article selection, and provides an overview of the history making accomplishments achieved during the SSBD effort, as well as, the follow-on NASA Shaped Sonic Boom Experiment (SSBE) Program, whose goal was to further evaluate the characteristics and robustness of shaped boom signatures.

Description: The transition from laminar to turbulent flow in helium jets discharged into air was studied using Rainbow Schlieren Deflectometry technique. In particular, the effects of buoyancy on jet oscillations and flow transition length were considered. Experiments to simulate microgravity were conducted in the 2.2s drop tower at NASA Glenn Research Center. The jet Reynolds numbers varied from 800 to1200 and the jet Richardson numbers ranged between 0.01 and 0.004. Schlieren images revealed substantial variations in the flow structure during the drop. Fast Fourier Transform (FFT) analysis of the data obtained in Earth gravity experiments revealed the existence of a discrete oscillating frequency in the transition region, which matched the frequency in the upstream laminar regime. In microgravity, the transition occurred farther downstream indicating laminarization of the jet in the absence of buoyancy. The amplitude of jet oscillations was reduced by up to an order of magnitude in microgravity. Results suggest that jet oscillations were buoyancy induced and that the brief microgravity period may not be sufficient for the oscillations to completely subside.

Description: The static aerodynamics for the Mars Exploration Rover (MER) aeroshell are presented. This aerodynamic database was an integral part of the end-to-end simulation used in pre- entry analysis for determining the MER entry design requirements for development of the MER entry system, as well as targeting the MER landing sites. The database was constructed using the same approach used for Mars Pathfinder (MPF). However, the MER aerodynamic database is of much higher fidelity and tailored to the MER entry trajectories. This set of data includes direct simulation Monte Carlo calculations covering the transitional regime of the entry trajectory and computational fluid dynamics calculations describing the aerodynamics in the hypersonic and supersonic continuum regimes. An overview of the methodology used to generate the data is given along with comparisons to important features in the MPF aerodynamics and related heritage data. The MER and MPF comparison indicates that trajectory specific data is required to properly model the flight characteristics of a.blunt entry capsule at Mars.

Description: A new computer program to design nozzles with non-uniform idow has been developed using the rotational method of characteristics (MOC). This program has been used to design a nozzle for the NASA's HYPULSE shock-expansion tunnel for use in scramjet engine tests at a Mach-15 flight-enthalpy condition. The nozzle has an area ratio of 9.5:l that expands the inflow from Mach 6 along the centerline to Mach 8.7. Although the density and Mach number vary radially at the exit due to the non-uniformities of the inflow, the MOC procedure produces exit flow that is parallel and has uniform static pressure. The design has been verified with CFD which compares favorably with the MOC solution.

Description: A source term model for steady micro jets was implemented into a non-proprietary Navier-Stokes computer code, OVERFLOW. The source term models the mass flow and momentum created by a steady blowing micro jet. The model is obtained by adding the momentum and mass flow created by the jet to the Navier-Stokes equations. The model was tested by comparing with data from numerical simulations of a single, steady micro jet on a flat plate in two and three dimensions. The source term model predicted the velocity distribution well compared to the two-dimensional plate using a steady mass flow boundary condition, which was used to simulate a steady micro jet. The model was also compared to two three-dimensional flat plate cases using a steady mass flow boundary condition to simulate a steady micro jet. The three-dimensional comparison included a case with a grid generated to capture the circular shape of the jet and a case without a grid generated for the micro jet. The case without the jet grid mimics the application of the source term. The source term model compared well with both of the three-dimensional cases. Comparisons of velocity distribution were made before and after the jet and Mach and vorticity contours were examined. The source term model allows a researcher to quickly investigate different locations of individual or several steady micro jets. The researcher is able to conduct a preliminary investigation with minimal grid generation and computational time.

Description: Turbulent separated flow over a two-dimensional hump is computed by solving the RANS equations with k - omega (SST) turbulence model for the baseline, steady suction and oscillatory blowing/suction flow control cases. The flow equations and the turbulent model equations are solved using a fifth-order accurate weighted essentially. nonoscillatory (WENO) scheme for space discretization and a third order, total variation diminishing (TVD) Runge-Kutta scheme for time integration. Qualitatively the computed pressure distributions exhibit the same behavior as those observed in the experiments. The computed separation regions are much longer than those observed experimentally. However, the percentage reduction in the separation region in the steady suction case is closer to what was measured in the experiment. The computations did not predict the expected reduction in the separation length in the oscillatory case. The predicted turbulent quantities are two to three times smaller than the measured values pointing towards the deficiencies in the existing turbulent models when they are applied to strong steady/unsteady separated flows.

Description: For applied aerodynamic computations using a general purpose Navier-Stokes code, the common practice of treating laminar to turbulent flow transition over a non-slip surface is somewhat arbitrary by either treating the entire flow as turbulent or forcing the flow to undergo transition at given trip locations in the computational domain. In this study, the possibility of using the PAB3D code, standard k-epsilon turbulence model, and the Girimaji explicit algebraic stresses model to mimic natural laminar to turbulent flow transition was explored. The sensitivity of flow transition with respect to two limiters in the standard k-epsilon turbulence model was examined using a flat plate and a 6:1 aspect ratio prolate spheroid for our computations. For the flat plate, a systematic dependence of transition Reynolds number on background turbulence intensity was found. For the prolate spheroid, the transition patterns in the three-dimensional boundary layer at different flow conditions were sensitive to the free stream turbulence viscosity limit, the reference Reynolds number and the angle of attack, but not to background turbulence intensity below a certain threshold value. The computed results showed encouraging agreements with the experimental measurements at the corresponding geometry and flow conditions.

Description: The effect of discrete contour bumps on reducing the transonic drag at off-design conditions on an airfoil have been examined. The research focused on fully-turbulent flow conditions, at a realistic flight chord Reynolds number of 30 million. State-of-the-art computational fluid dynamics methods were used to design a new baseline airfoil, and a family of fixed contour bumps. The new configurations were experimentally evaluated in the 0.3-m Transonic Cryogenic Tunnel at the NASA Langley Research center, which utilizes an adaptive wall test section to minimize wall interference. The computational study showed that transonic drag reduction, on the order of 12% - 15%, was possible using a surface contour bump to spread a normal shock wave. The computational study also indicated that the divergence drag Mach number was increased for the contour bump applications. Preliminary analysis of the experimental data showed a similar contour bump effect, but this data needed to be further analyzed for residual wall interference corrections.

Description: Production of single-walled carbon nanotubes (SWNT) has taken place for a number of years and by a variety of methods such-as laser ablation, chemical vapor deposition, and arc-jet ablation. Yet, little is actually understood about the exact chemical kinetics and processes that occur in SWNT formation. In recent time, NASA Johnson Space Center has devoted a considerable effort to the experimental evaluation of the laser ablation production process for SWNT originally developed at Rice University. To fully understand the nature of the laser ablation process it is necessary to understand the development of the carbon plume dynamics within the laser ablation oven. The present work is a continuation of previous studies into the efforts to model plume dynamics using computational fluid dynamics (CFD). The ultimate goal of the work is to improve understanding of the laser ablation process, and through that improved understanding, refine the laser ablation production of SWNT. Fig. 1 shows a basic schematic of the laser-ablation oven at NASA-JSC. Construction of the facility is simple in concept. Two concentric quartz tubes of 1.5 mm thickness form the inner and outer tubes with inside diameters of 2.2 and 5.08 cm respectively. At one end of the inner tube are located two 60 Hz pulsed lasers operating at 1064 nm and 532 nm wavelength with beam diameters of 5 mm aligned coaxially with the longitudinal axis of the inner quartz tube. For standard nanotube production runs, a 10 ns 532 nm pulse is followed 50 ns later by a 10 ns 1064 nm pulse. Each pulse is of 300 mJ energy. A target of carbon graphite with approximately 1% nickel and cobalt catalysts is located at the other end of the inner quartz tube. In the ordinary processing of SWNT, a base flow of 100 sccm of argon is maintained from the laser location and exits past the carbon target at a pressure of 66.7 kPa. These conditions yield a baseline mass flow through the chamber of 2.723x10(exp -6)kg/s of argon. The whole oven facility is heated to a temperature of 1473 K prior to nanotube production runs. Upon laser irradiation, part of the carbon target ablates immediately and forms a carbon vapor plume that penetrates into the argon base flow towards the laser initially at supersonic velocities. In the previous studies either a single carbon species, C3, was used to model the plume development, or a simplified 11 species carbon model reduced down from the Krestinin and Moravsky model for full fullerene chemistry was used. While both of these options yielded significant results, it was felt that the actual chemistry occ&g in the carbon plume might have a greater affect on the plume than assumed. Indeed, in the earlier attempts at modeling the carbon plume, several thermophysical characteristics could never be matched to experimental observations of plume development - mainly the propagation distance of the plume itself. In the present study, two additional chemistry models will be used to duplicate the previous studies simulations of the carbon plume. The first chemistry model used in this study is again a reduced form of the Krestinin and Moravsky rates. However the highest order carbon species allowed has been increased from C6 to C30 - therefore allowing the simulation of up to a half of the standard C60 fullerene. The second chemistry model investigated is a reduced form of a full carbon nanotube model developed at NASA-JSC. The C30 studies have already been accomplished at the present time, and the reduced SWNT model studies are currently underway. To pursue the current study, one sacrifice had to be made in that the simulation grid spacing had to be increased from 0.5 mm spacing to 1 mm spacing for the sake of computational efficiency since computational effort is proportional to the square of the number of grid points multiplied by the number of species considered. propagation that is far more in line with the experimental results observed by Puretzky et al as shown in Fig. 2. e C6 studies had yielded a far greater propagation in previous studies. In addition, chemical species development with the C30 model indicates that many higher order carbon species are produced outside of the plume proper (indicated by plotting contours of the background argon concentrations in Fig. 3) - this result was not observed in previous studies. In fact, some species primarily occurred outside of the plume itself - as shown for C27 in Fig. 4 when it is compared to Fig. 3. It could be asserted that this has occurred because all of the C27 in the plume had already been consumed in the formation of C30, but this does not seem to be indicated over time. Several other factors that arose in the previous studies have also been made more clear by the use of the higher order chemical models - one being that the use of c6 as an indicator species was mistaken. C6 is the only carbon species in the previous studies that was not injected into the flowfield as a boundary condition; it was therefore hoped that this species would provide insight into the formation of higher order carbon species for comparison to full SWNT production. But, when the plot of total mass in the plume is examined on a species by species basis in Fig. 5, it is seen that Cg was a fairly insignificant contributor to the total carbon mass in the plume and would not provide information on higher order carbon formation. the thermophysical characteristics of the carbon plume as well as simulate the carbon plume using the reduced SWNT model to provide an even better simulation of full chemistry effects upon plume propagation.

Description: A transport equation for the intermittency factor is employed to predict transitional flows under the effects of pressure gradients, freestream turbulence intensities, Reynolds number variations, flow separation and reattachment. and unsteady wake-blade interactions representing diverse operating conditions encountered in low-pressure turbines. The intermittent behaviour of the transitional flows is taken into account and incorporated into computations by modifying the eddy viscosity, Mu(sub t), with the intermittency factor, gamma. Turbulent quantities are predicted by using Menter's two-equation turbulence model (SST). The onset location of transition is obtained from correlations based on boundary-layer momentum thickness, acceleration parameter, and turbulence intensity. The intermittency factor is obtained from a transport model which can produce both the experimentally observed streamwise variation of intermittency and a realistic profile in the cross stream direction. The intermittency transport model is tested and validated against several well documented low pressure turbine experiments ranging from flat plate cases to unsteady wake-blade interaction experiments. Overall, good agreement between the experimental data and computational results is obtained illustrating the predicting capabilities of the model and the current intermittency transport modelling approach for transitional flow simulations.

Description: This viewgraph presentation reviews NASA's basics of flight, controlling flight, and flight systems.

Description: We have performed 3-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations of jet formation from an accretion disk with/without initial perturbation around a rotating black hole. We input a sinusoidal perturbation (m = 5 mode) in the rotation velocity of the accretion disk. The simulation results show the formation of a relativistic jet from the accretion disk. Although the initial perturbation becomes weakened by the coupling among different modes, it survives and triggers lower modes. As a result, complex non-axisymmetric density structure develops in the disk and the jet. Newtonian MHD simulations of jet formation with a non-axisymmetric mode show the growth of the m = 2 mode but GRMHD simulations cannot see the clear growth of the m = 2 mode.

Description: A mathematical model was developed to simulate flow control applications using plasma actuators. The effects of the plasma actuators on the external flow are incorporated into Navier Stokes computations as a body force vector. In order to compute this body force vector, the model solves two additional equations: one for the electric field due to the applied AC voltage at the electrodes and the other for the charge density representing the ionized air. The model is calibrated against an experiment having plasma-driven flow in a quiescent environment and is then applied to simulate a low pressure turbine flow with large flow separation. The effects of the plasma actuator on control of flow separation are demonstrated numerically.

Description: The hydrodynamic force experienced by a spherical-cap drop moving on a solid surface is obtained from two approximate analytical solutions and used to predict the quasi-steady speed of the drop in a wettability gradient. One solution is based on approximation of the shape of the drop as a collection of wedges, and the other is based on lubrication theory. Also, asymptotic results from both approximations for small contact angles, as well as an asymptotic result from lubrication theory that is good when the length scale of the drop is large compared with the slip length, are given. The results for the hydrodynamic force also can be used to predict the quasi-steady speed of a drop sliding down an incline.

Description: Turbulent thin-layer, Reynolds-Averaged Navier-Stokes solutions, based on a multi-block structured grid, are presented for a 65 deg delta wing having either a sharp leading edge (SLE) or blunt leading edge (BLE) geometry. The primary objective of the study is to assess the prediction capability of the method for simulating the leading-edge flow separation and the ensuing vortex flow characteristics. Computational results are obtained for two angles of attack of approximately 13 and 20 deg, at free-stream Mach number of 0.40 and Reynolds number of 6 million based on the wing mean aerodynamic chord. The effects of two turbulence models of Baldwin-Lomax with Degani-Schiff (BL/DS) and the Spalart-Allmaras (SA) on the numerical results are also discussed. The computations also explore the effects of two numerical flux-splitting schemes, i.e., flux difference splitting (fds) and flux vector splitting (fvs), on the solution development and convergence characteristics. The resulting trends in solution sensitivity to grid resolution for the selected leading-edge geometries, angles of attack, turbulence models and flux splitting schemes are also presented. The validity of the numerical results is evaluated against a unique set of experimental wind-tunnel data that was obtained in the National Transonic Facility at the NASA Langley Research Center.

Description: In this grant, experimental, numerical and analytical studies of heat transfer in a thin liquid film flowing over a rotating disk have been conducted. Heat transfer coefficients were measured experimentally in a rotating disk heat transfer apparatus where the disk was heated from below with electrical resistance heaters. The heat transfer measurements were supplemented by experimental characterization of the liquid film thickness using a novel laser based technique. The heat transfer measurements show that the disk rotation plays an important role on enhancement of heat transfer primarily through the thinning of the liquid film. Experiments covered both momentum and rotation dominated regimes of the flow and heat transfer in this apparatus. Heat transfer measurements have been extended to include evaporation and nucleate boiling and these experiments are continuing in our laboratory. Empirical correlations have also been developed to provide useful information for design of compact high efficiency heat transfer devices. The experimental work has been supplemented by numerical and analytical analyses of the same problem. Both numerical and analytical results have been found to agree reasonably well with the experimental results on liquid film thickness and heat transfer Coefficients/Nusselt numbers. The numerical simulations include the free surface liquid film flow and heat transfer under disk rotation including the conjugate effects. The analytical analysis utilizes an integral boundary layer approach from which

Description: An experimental study has been performed to develop a large force and moment aerodynamic data set on a slender axisymmetric missile configuration having cruciform strakes and in-line control tail fins. The data include six-component balance measurements of the configuration aerodynamics and three-component measurements on all four tail fins. The test variables include angle of attack, roll angle, Mach number, model buildup, strake length, nose size, and tail fin deflection angles to provide pitch, yaw, and roll control. Test Mach numbers ranged from 0.60 to 4.63. The entire data set is presented on a CD-ROM that is attached to this paper. The CD-ROM also includes extensive plots of both the six-component configuration data and the three-component tail fin data. Selected samples of these plots are presented in this paper to illustrate the features of the data and to investigate the effects of the test variables.

Description: In a recent paper by Zinoviev and Bies in this Journal, the authors have claimed that the well-known theoretical results of Curle and Ffowcs Williams and Hawkings (FW-H) are incorrect. This claim is categorically refuted below and serious errors are pointed out.

Description: Experiments were conducted in Earth gravity and microgravity to acquire quantitative data on near field flow structure of helium jets injected into air. Microgravity conditions were simulated in the 2.2- second drop tower at NASA Glenn Research Center. The jet flow was observed by quantitative rainbow schlieren deflectometry, a non-intrusive line of site measurement technique for the whole field. The flow structure was characterized by distributions of angular deflection and helium mole percentage obtained from color schlieren images taken at 60 Hz. Results show that the jet in microgravity was up to 70 percent wider than that in Earth gravity. The global jet flow oscillations observed in Earth gravity were absent in microgravity, providing direct experimental evidence that the flow instability in the low density jet was buoyancy induced. The paper provides quantitative details of temporal flow evolution as the experiment undergoes change in gravity in the drop tower.

Description: Flow structure of momentum-dominated helium jets discharged vertically into ambient air was investigated using a high-speed rainbow schlieren deflectometry (RSD) apparatus operated at up to 2000 Hz. The operating parameters, i.e., Reynolds number and Richardson number were varied independently to examine the self-excited, flow oscillatory behavior over a range of experimental conditions. Measurements revealed highly periodic oscillations in the laminar region at a unique frequency as well as high regularity in the flow transition and initial turbulent regions. The buoyancy was shown to affect the oscillation frequency and the distance from the jet exit to the flow transition plane. Instantaneous helium concentration contours across the field of view revealed changes in the jet flow structure and the evolution of the vortical structures during an oscillation cycle. A cross-correlation technique was applied to track the vortices and to find their convection velocity. Time-traces of helium concentration at different axial locations provided detailed information about the oscillating flow.

Description: The detached eddy simulation of film cooling has been utilized for a proprietary GE plate-pipe configuration. The blowing ratio was 2.02, the velocity ratio was 1.26, and the temperature ratio was 1.61. Results indicate that the mixing processes downstream of the hole are highly anisotropic. DES solution shows its ability to depict the dynamic nature of the flow and capture the asymmetry present in temperature and velocity distributions. Further, comparison between experimental and DES time-averaged effectiveness is satisfactory. Numerical values of span-averaged effectiveness show better prediction of the experimental values at downstream locations than a steady state Glenn HT solution. While the DES method shows obvious promise, there are several issues that need further investigation. Despite an accurate prediction in the hole vicinity, the simulation still falls short in the region x = 10d to 100d. This should be investigated. Also the model used flat plate. Actual turbine blade should be modeled in the future if additional finding is available.

Description: The aim of this project was is to advance a microstructural understanding of foam and emulsion flows. The dynamics of individual surfactant-covered drops and well as the collective behavior of dilute and concentrated was explored using numerical simulations. The long-range goal of this work is the formulation of reliable microphysically-based statistical models of emulsion flows.

Description: A brief review of the evolutionary progress in computational aerothermodynamics is presented. The current status of computational aerothermodynamics is then discussed, with emphasis on its capabilities and limitations for contributions to the design process of hypersonic vehicles. Some topics to be highlighted include: (1) aerodynamic coefficient predictions with emphasis on high temperature gas effects; (2) surface heating and temperature predictions for thermal protection system (TPS) design in a high temperature, thermochemical nonequilibrium environment; (3) methods for extracting and extending computational fluid dynamic (CFD) solutions for efficient utilization by all members of a multidisciplinary design team; (4) physical models; (5) validation process and error estimation; and (6) gridding and solution generation strategies. Recent experiences in the design of X-33 will be featured. Computational aerothermodynamic contributions to Mars Path finder, METEOR, and Stardust (Comet Sample return) will also provide context for this discussion. Some of the barriers that currently limit computational aerothermodynamics to a predominantly reactive mode in the design process will also be discussed, with the goal of providing focus for future research.

Description: Effects of buoyancy on transition from laminar to turbulent flow are presented for momentum-dominated helium jet injected into ambient air. The buoyancy was varied in a 2.2-sec drop tower facility without affecting the remaining operating parameters. The jet flow in Earth gravity and microgravity was visualized using the rainbow schlieren deflectometry apparatus. Results show significant changes in the flow structure and transition behavior in the absence of buoyancy.

Description: The objective of this study was to determine the effects of buoyancy on the absolute instability of low-density gas jets injected into high-density gas mediums. Most of the existing analyses of low-density gas jets injected into a high-density ambient have been carried out neglecting effects of gravity. In order to investigate the influence of gravity on the near-injector development of the flow, a spatio-temporal stability analysis of a low-density round jet injected into a high-density ambient gas was performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The Briggs-Bers criterion was combined with the spatio-temporal stability analysis to determine the nature of the absolute instability of the jet whether absolutely or convectively unstable. The roles of the density ratio, Froude number, Schmidt number, and the lateral shift between the density and velocity profiles on the absolute instability of the jet were determined. Comparisons of the results with previous experimental studies show good agreement when the effects of these variables are combined together. Thus, the combination of these variables determines how absolutely unstable the jet will be.

Description: Using laser-induced thermal acoustics, we demonstrate non-intrusive and remote sound speed and temperature measurements over the range 10 - 45 C in liquid water. Averaged accuracy of sound speed and temperature measurements (10 s) are 0.64 m/s and 0.45 C respectively. Single-shot precisions based on one standard deviation of 100 or greater samples range from 1 m/s to 16.5 m/s and 0.3 C to 9.5 C for sound speed and temperature measurements respectively. The time resolution of each single-shot measurement was 300 nsec.

Description: A methodology is developed and implemented to mitigate the lengthy software development cycle typically associated with constructing a discrete adjoint solver for aerodynamic simulations. The approach is based on a complex-variable formulation that enables straightforward differentiation of complicated real-valued functions. An automated scripting process is used to create the complex-variable form of the set of discrete equations. An efficient method for assembling the residual and cost function linearizations is developed. The accuracy of the implementation is verified through comparisons with a discrete direct method as well as a previously developed handcoded discrete adjoint approach. Comparisons are also shown for a large-scale configuration to establish the computational efficiency of the present scheme. To ultimately demonstrate the power of the approach, the implementation is extended to high temperature gas flows in chemical nonequilibrium. Finally, several fruitful research and development avenues enabled by the current work are suggested.

Description: As part of the Space Shuttle Return to Flight (RTF) program, Marshall Space Flight Center (MSFC) performed computational fluid dynamics (CFD) analysis to define the velocity flowfields around the Shuttle stack at liftoff. These CFD predicted velocity flowfields were used in debris transport analysis (DTA). High speed flows such as plumes induce or 'entrain' mass from the surrounding environment. Previous work had shown that CFD analysis over-predicts plume induced flows. Therefore, the DTA would tend to 1) predict more debris impacts, and 2) the debris velocity (and kinetic energy) of those impacts would be too high. At a November, 2004 peer-review it was recommended that the Liftoff DTA team quantify the uncertainty in the DTA caused by the CFD's over prediction of plume induced flow. To do so, the Liftoff DTA team needed benchmark quality data for plume induced flow to quantify the CFD accuracy and its effect on the DTA. MSFC's Nozzle Test Facility (NTF) conducted the "Nozzle Induced Flows test, P#2456" to obtain experimental data for plume induced flows for nozzle flow exhausting into q quiescent freestream. Planning for the test began in December, 2004 and the experimental data was obtained in February and March of 2005. The funding for this test was provided by MSFC's Space Shuttle Propulsion Systems Integration and Engineering office.

Description: A gas-liquid separator uses a helical passageway to impart a spiral motion to a fluid passing therethrough. The centrifugal fore generated by the spiraling motion urges the liquid component of the fluid radially outward which forces the gas component radially inward. The gas component is then filtered through a gas-permeable, liquid-impervious membrane and discharged through a central passageway.

Description: A dynamic wake avoidance system utilizes aircraft and atmospheric parameters readily available in flight to model and predict airborne wake vortices in real time. A novel combination of algorithms allows for a relatively simple yet robust wake model to be constructed based on information extracted from a broadcast. The system predicts the location and movement of the wake based on the nominal wake model and correspondingly performs an uncertainty analysis on the wake model to determine a wake hazard zone (no fly zone), which comprises a plurality of wake planes, each moving independently from another. The system selectively adjusts dimensions of each wake plane to minimize spatial and temporal uncertainty, thereby ensuring that the actual wake is within the wake hazard zone. The predicted wake hazard zone is communicated in real time directly to a user via a realistic visual representation. In an example, the wake hazard zone is visualized on a 3-D flight deck display to enable a pilot to visualize or see a neighboring aircraft as well as its wake. The system substantially enhances the pilot's situational awareness and allows for a further safe decrease in spacing, which could alleviate airport and airspace congestion.

Description: The computational fluid dynamics code FLUENT was used to study Rayleigh instability at large temperature differences in a sealed gas-filled enclosure with a cold top surface and a heated bottom wall (Benard problem). Both steady state and transient calculations were performed. The results define the boundaries of instability in a system depending on the geometry, temperature and pressure. It is shown that regardless of how fast the bottom-wall temperature can be ramped up to minimize the time spent in the unstable region of fluid motion, the eventual stability of the system depends on the prevailing final pressure after steady state has been reached. Calculations also show that the final state of the system can be different depending on whether the result is obtained via a steady-state solution or is reached by transient calculations. Changes in the slope of the pressure-versus-time curve are found to be a very good indicator of changes in the flow patterns in the system.

Description: Circulation Control technology is a very effective way of achieving high lift forces required by aircraft during take-off and landing. This technology can also directly control the flow field over the wing. Compared to a conventional high-lift system, a Circulation Control Wing (CCW) can generate comparable or higher lift forces during take-off/landing with fewer or no moving parts and much less complexity. In this work, an unsteady three-dimensional Navier-Stokes analysis procedure has been developed and applied to Circulation Control Wing configurations. The effects of 2-D steady jets and 2-D pulsed jets on the aerodynamic performance of CCW airfoils have been investigated. It is found that a steady jet can generate very high lift at zero angle of attack without stall, and that a small amount of blowing can eliminate vortex shedding at the trailing edge, a potential noise source. It is also found that a pulsed jet can achieve the same high lift as a steady jet at lower mass flow rates, especially at a high frequency, and that the Strouhal number has a more dominant effect on the pulsed jet performance than just the frequency or the free-stream velocity.

Description: Closed-form analytic solutions useful for the design of capillary flows in a variety of containers possessing interior corners were recently collected and reviewed. Low-g drop tower and aircraft experiments performed at NASA to date show excellent agreement between theory and experiment for perfectly wetting fluids. The analytical expressions are general in terms of contact angle, but do not account for variations in contact angle between the various surfaces within the system. Such conditions may be desirable for capillary containment or to compute the behavior of capillary corner flows in containers consisting of different materials with widely varying wetting characteristics. A simple coordinate rotation is employed to recast the governing system of equations for flows in containers with interior corners with differing contact angles on the faces of the corner. The result is that a large number of capillary driven corner flows may be predicted with only slightly modified geometric functions dependent on corner angle and the two (or more) contact angles of the system. A numerical solution is employed to verify the new problem formulation. The benchmarked computations support the use of the existing theoretical approach to geometries with variable wettability. Simple experiments to confirm the theoretical findings are recommended. Favorable agreement between such experiments and the present theory may argue well for the extension of the analytic results to predict fluid performance in future large length scale capillary fluid systems for spacecraft as well as for small scale capillary systems on Earth.

Description: An integral analysis of hydrodynamics and heat transfer in a thin liquid film flowing over a rotating disk surface is presented for both constant temperature and constant heat flux boundary conditions. The model is found to capture the correct trends of the liquid film thickness variation over the disk surface and compare reasonably well with experimental results over the range of Reynolds and Rossby numbers covering both inertia and rotation dominated regimes. Nusselt number variation over the disk surface shows two types of behavior. At low rotation rates, the Nusselt number exhibits a radial decay with Nusselt number magnitudes increasing with higher inlet Reynolds number for both constant wall temperature and heat flux cases. At high rotation rates, the Nusselt number profiles exhibit a peak whose location advances radially outward with increasing film Reynolds number or inertia. The results also compare favorably with the full numerical simulation results from an earlier study as well as with the reported experimental results.

Description: Heat transfer associated with turbulent flow on a step-heated or cooled section of a flat plate at zero angle of attack with an insulated starting section was computationally modeled using the GASP Navier-Stokes code. The algebraic eddy viscosity model of Baldwin-Lomax and the turbulent two-equation models, the K- model and the Shear Stress Turbulent model (SST), were employed. The variations from uniformity of the imposed experimental temperature profile were incorporated in the computations. The computations yielded satisfactory agreement with the experimental results for all three models. The Baldwin- Lomax model showed the closest agreement in heat transfer, whereas the SST model was higher and the K-omega model was yet higher than the experiments. In addition to the step temperature distribution case, computations were also carried out for a uniformly heated or cooled plate. The SST model showed the closest agreement with the Von Karman analogy, whereas the K-omega model was higher and the Baldwin-Lomax was lower.

Description: The vaporization behavior of the A1-O system has been studied on numerous occasions but significant uncertainties remain. The origin of this uncertainty must be understood before A1-O vaporization behavior can be accurately determined. The condensation of A12O3 and clogging of the effusion orifice is a difficult problem for the Knudsen effusion technique that influences the measured vaporization behavior but has only received limited attention. This study reconsiders this behavior in detail. A new theory for A12O3 condensation is proposed together with procedures that will improve the measured thermodynamic properties of A1-O vaporization.

Description: Computational Fluid Dynamics (CFD) simulations are routinely performed as part of the design process of most fluid handling devices. In order to efficiently and effectively use the results of a CFD simulation, visualization tools are often used. These tools are used in all stages of the CFD simulation including pre-processing, interim-processing, and post-processing, to interpret the results. Each of these stages requires visualization tools that allow one to examine the geometry of the device, as well as the partial or final results of the simulation. An engineer will typically generate a series of contour and vector plots to better understand the physics of how the fluid is interacting with the physical device. Of particular interest are detecting features such as shocks, re-circulation zones, and vortices (which will highlight areas of stress and loss). As the demand for CFD analyses continues to increase the need for automated feature extraction capabilities has become vital. In the past, feature extraction and identification were interesting concepts, but not required in understanding the physics of a steady flow field. This is because the results of the more traditional tools like; isc-surface, cuts and streamlines, were more interactive and easily abstracted so they could be represented to the investigator. These tools worked and properly conveyed the collected information at the expense of a great deal of interaction. For unsteady flow-fields, the investigator does not have the luxury of spending time scanning only one "snapshot" of the simulation. Automated assistance is required in pointing out areas of potential interest contained within the flow. This must not require a heavy compute burden (the visualization should not significantly slow down the solution procedure for co-processing environments). Methods must be developed to abstract the feature of interest and display it in a manner that physically makes sense.

Description: A computational study of the separated flow through a 2-D asymmetric subsonic diffuser has been performed. The Wind Computational Fluid Dynamics code is used to predict the separation and reattachment behavior for an incompressible diffuser flow. The diffuser inlet flow is a two-dimensional, turbulent, and fully-developed channel flow with a Reynolds number of 20,000 based on the centerline velocity and the channel height. Wind solutions computed with the Menter SST, Chien k-epsilon, Spalart-Allmaras and Explicit Algebraic Reynolds Stress turbulence models are compared with experimentally measured velocity profiles and skin friction along the upper and lower walls. In addition to the turbulence model study, the effects of grid resolution and use of wall functions were investigated. The grid studies varied the number of grid points across the diffuser and varied the initial wall spacing from y(sup +) = 0.2 to 60. The wall function study assessed the applicability of wall functions for analysis of separated flow. The SST and Explicit Algebraic Stress models provide the best agreement with experimental data, and it is recommended wall functions should only be used with a high level of caution.

Description: In Spring of 2005, the NASA Engineering Safety Center (NESC) was engaged by the Space Shuttle Program (SSP) to peer review the suite of analytical tools being developed to support the determination of impact and damage tolerance of the Orbiter Thermal Protection Systems (TPS). The NESC formed an independent review team with the core disciplines of materials, flight sciences, structures, mechanical analysis and thermal analysis. The Math Model Tools reviewed included damage prediction and stress analysis, aeroheating analysis, and thermal analysis tools. Some tools are physics-based and other tools are empirically-derived. Each tool was created for a specific use and timeframe, including certification, real-time pre-launch assessments, and real-time on-orbit assessments. The tools are used together in an integrated strategy for assessing the ramifications of impact damage to tile and RCC. The NESC teams conducted a peer review of the engineering data package for each Math Model Tool. This report contains the summary of the team observations and recommendations from these reviews.

Description: A method for estimating the sonic-boom overpressures from a conceptual aircraft where the lift is carried by both a canard and a wing during supersonic cruise is presented and discussed. Computer codes used for the prediction of the aerodynamic performance of the wing, the canard-wing interference, the nacelle-wing interference, and the sonic-boom overpressures are identified and discussed as the procedures in the method are discussed. A canard-wing supersonic-cruise concept was used as an example to demonstrate the application of the method.

Description: An investigation using a survey rake with 11 five-hole pyramid-head probes has been conducted in the Langley Transonic Dynamics Tunnel (TDT) to measure the test section flow angularity. Flow measurements were made in a 10-ft square grid centered about the test section centerline at a single streamwise location for nine Mach numbers ranging from 0.50 to 1.19 at dynamic pressures of 100 and 225 pounds per square foot. Test section flow angularity was found to be minimal with a generally random flow pattern. Corrections for survey rake induced in-plane flow were determined to be necessary; however, corrections for rake induced lift effects were not required.

Description: A new computer program has been developed called ASP3D (Advanced Small Perturbation 3D), which solves the small perturbation potential flow equation in an advanced form including mass-consistent surface and trailing wake boundary conditions, and entropy, vorticity, and viscous effects. The purpose of the program is for unsteady aerodynamic and aeroelastic analyses, especially in the nonlinear transonic flight regime. The program exploits the simplicity of stationary Cartesian meshes with the movement or deformation of the configuration under consideration incorporated into the solution algorithm through a planar surface boundary condition. The new ASP3D code is the result of a decade of developmental work on improvements to the small perturbation formulation, performed while the author was employed as a Senior Research Scientist in the Configuration Aerodynamics Branch at the NASA Langley Research Center. The ASP3D code is a significant improvement to the state-of-the-art for transonic aeroelastic analyses over the CAP-TSD code (Computational Aeroelasticity Program Transonic Small Disturbance), which was developed principally by the author in the mid-1980s. The author is in a unique position as the developer of both computer programs to compare, contrast, and ultimately make conclusions regarding the underlying formulations and utility of each code. The paper describes the salient features of the ASP3D code including the rationale for improvements in comparison with CAP-TSD. Numerous results are presented to demonstrate the ASP3D capability. The general conclusion is that the new ASP3D capability is superior to the older CAP-TSD code because of the myriad improvements developed and incorporated.

Description: A mesh deformation scheme is developed for a structured multi-block Navier-Stokes code consisting of two steps. The first step is a finite element solution of either user defined or automatically generated macro-elements. Macro-elements are hexagonal finite elements created from a subset of points from the full mesh. When assembled, the finite element system spans the complete flow domain. Macro-element moduli vary according to the distance to the nearest surface, resulting in extremely stiff elements near a moving surface and very pliable elements away from boundaries. Solution of the finite element system for the imposed boundary deflections generally produces smoothly varying nodal deflections. The manner in which distance to the nearest surface has been found to critically influence the quality of the element deformation. The second step is a transfinite interpolation which distributes the macro-element nodal deflections to the remaining fluid mesh points. The scheme is demonstrated for several two-dimensional applications.

Description: A planar laser-induced fluorescence (PLIF) technique is discussed and applied to measurement of time-averaged values of velocity and temperature in an I(sub 2)-seeded N(sub 2) hypersonic free jet facility. Using this technique, a low temperature, non-reacting, hypersonic flow over a simplified model of a reaction control system (RCS) was investigated. Data are presented of rarefied Mach 12 flow over a sharp leading edge flat plate at zero incidence, both with and without an interacting jet issuing from a nozzle built into the plate. The velocity profile in the boundary layer on the plate was resolved. The slip velocity along the plate, extrapolated from the velocity profile data, varied from nearly 100% down to 10% of the freestream value. These measurements are compared with results of a DSMC solution. The velocity variation along the centerline of a jet issuing from the plate was measured and found to match closely with the correlation of Ashkenas and Sherman. The velocity variation in the oblique shock terminating the jet was resolved sufficiently to measure the shock wave thickness.

Description: Experiments were performed on the collision of a solid sphere with a nearly horizontal flat surface covered with a thin layer of viscous liquid. High-speed collisions were obtained by dropping the ball onto the surface from various heights, using gravitational acceleration. Low-speed collisions were obtained using pendulums with long strings or by launching the balls at low velocities in the reduced-gravity environment of parabolic flight. The sphere bounces only when the impact velocity exceeds a critical value. The coefficient of restitution (ratio of rebound velocity to impact velocity) increases with increasing impact velocity above the critical value, indicating the increasing relative importance of elastic deformation to viscous dissipation. The critical impact velocity increases, and the coefficient of restitution decreases, with increasing viscosity or thickness of the liquid layer and with decreasing density or size of the sphere. The ratio of the wet and dry coefficients is expressed as a function of the Stokes number (ratio of particle inertia and viscous forces), showing good agreement between theory and experiment. Similar experiments were performed with the flat surface inclined at various angles to the approaching sphere. A modified Stokes number, which is a measure of the ratio of inertia of the sphere in the normal direction to the viscous forces exerted by the fluid layer, was used for the analysis of oblique collisions. Even for these oblique collisions, it was found that no rebound of the ball was observed below a certain critical Stokes number. The coefficient of normal restitution, defined as a ratio of normal rebound velocity to normal approach velocity, was found to increase beyond the critical Stokes number and even out as it approaches the value for dry restitution at high Stokes numbers. It was also found that, for smooth spheres like steel, the normal restitution at the same modified Stokes number is independent of the angle of impact. The tangential coefficient of restitution, defined as the ratio of tangential rebound velocity to tangential approach velocity, is found to be nearly unity, except for very low approach velocities. Thus, as a first approximation, the theories that predict the coefficient of restitution for head-on wet collisions can be extended to predict the coefficient of normal restitution for oblique wet collisions. Additional experiments were performed with soft surfaces in which a porous cloth or sponge layer was placed over the hard, flat surface. In these experiments, the coefficient of restitution was found to decrease with increasing impact velocity, due to inelastic losses in the soft material. A model combining inelastic deformation and flow through porous media was developed to describe these findings.

Description: Two optically based flame-detection systems have been developed for use in NASA Langley's 8-Foot High-Temperature Tunnel (8-ft HTT). These systems are used to detect the presence and stability of the main-burner and pilot-level flames during facility operation. System design considerations will be discussed, and a detailed description of the system components and circuit diagrams will be provided in the Appendices of this report. A more detailed description of the manufacturing process used in the fabrication of the fiber-optic probes is covered in NASA TM-2001-211233.

Description: The paper experimentally studies the effects of periodic unsteady wake flow and different Reynolds numbers on boundary layer development, separation and re-attachment along the suction surface of a low pressure turbine blade. The experimental investigations were performed on a large scale, subsonic unsteady turbine cascade research facility at Turbomachinery Performance and Flow Research Laboratory (TPFL) of Texas A&M University. The experiments were carried out at Reynolds numbers of 110,000 and 150,000 (based on suction surface length and exit velocity). One steady and two different unsteady inlet flow conditions with the corresponding passing frequencies, wake velocities, and turbulence intensities were investigated. The reduced frequencies chosen cover the operating range of LP turbines. In addition to the unsteady boundary layer measurements, surface pressure measurements were performed. The inception, onset, and the extent of the separation bubble information collected from the pressure measurements were compared with the hot wire measurements. The results presented in ensemble-averaged, and the contour plot forms help to understand the physics of the separation phenomenon under periodic unsteady wake flow and different Reynolds number. It was found that the suction surface displayed a strong separation bubble for these three different reduced frequencies. For each condition, the locations defining the separation bubble were determined carefully analyzing and examining the pressure and mean velocity profile data. The location of the boundary layer separation was dependent of the Reynolds number. It is observed that starting point of the separation bubble and the re-attachment point move further downstream by increasing Reynolds number from 110,000 to 150,000. Also, the size of the separation bubble is smaller when compared to that for Re=110,000.

Description: The present study, which is the first of a series of investigations dealing with specific issues of low pressure turbine (LPT) boundary layer aerodynamics, is aimed at providing detailed unsteady boundary flow information to understand the underlying physics of the inception, onset, and extent of the separation zone. A detailed experimental study on the behavior of the separation zone on the suction surface of a highly loaded LPT-blade under periodic unsteady wake flow is presented. Experimental investigations were performed at Texas A&M Turbomachinery Performance and Flow Research Laboratory using a large-scale unsteady turbine cascade research facility with an integrated wake generator and test section unit. To account for a high flow deflection of LPT-cascades at design and off-design operating points, the entire wake generator and test section unit including the traversing system is designed to allow a precise angle adjustment of the cascade relative to the incoming flow. This is done by a hydraulic platform, which simultaneously lifts and rotates the wake generator and test section unit. The unit is then attached to the tunnel exit nozzle with an angular accuracy of better than 0.05 , which is measured electronically. Utilizing a Reynolds number of 110,000 based on the blade suction surface length and the exit velocity, one steady and two different unsteady inlet flow conditions with the corresponding passing frequencies, wake velocities and turbulence intensities are investigated using hot-wire anemometry. In addition to the unsteady boundary layer measurements, blade surface pressure measurements were performed at Re=50,000, 75,000, 100,000, and 125,000 at one steady and two periodic unsteady inlet flow conditions. Detailed unsteady boundary layer measurement identifies the onset and extent of the separation zone as well as its behavior under unsteady wake flow. The results presented in ensemble-averaged and contour plot forms contribute to understanding the physics of the separation phenomenon under periodic unsteady wake flow. Several physical mechanisms are discussed.

Description: A method for reducing drag upon a blunt-based vehicle by adaptively increasing forebody roughness to increase drag at the roughened area of the forebody, which results in a decrease in drag at the base of this vehicle, and in total vehicle drag.

Description: The IATCS coolant has experienced a number of anomalies in the time since the US Lab was first activated on Flight 5A in February 2001. These have included: 1) a decrease in coolant pH, 2) increases in inorganic carbon, 3) a reduction in phosphate buffer concentration, 4) an increase in dissolved nickel and precipitation of nickel salts, and 5) increases in microbial concentration. These anomalies represent some risk to the system, have been implicated in some hardware failures and are suspect in others. The ISS program has conducted extensive investigations of the causes and effects of these anomalies and has developed a comprehensive program to remediate the coolant chemistry of the on-orbit system as well as provide a robust and compatible coolant solution for the hardware yet to be delivered. The remediation steps include changes in the coolant chemistry specification, development of a suite of new antimicrobial additives, and development of devices for the removal of nickel and phosphate ions from the coolant. This paper presents an overview of the anomalies, their known and suspected system effects, their causes, and the actions being taken to remediate the coolant.

Description: A high Reynolds number wind tunnel test was conducted to assess Reynolds number effects on the aerodynamic performance characteristics of a realistic, second-generation supersonic transport concept. The tests included longitudinal studies at transonic and low-speed, high-lift conditions across a range of chord Reynolds numbers (8 million to 120 million). Results presented focus on Reynolds number and static aeroelastic sensitivities at Mach 0.30 and 0.90 for a configuration without a tail. Static aeroelastic effects, which mask Reynolds number effects, were observed. Reynolds number effects were generally small and the drag data followed established trends of skin friction as a function of Reynolds number. A more nose-down pitching moment was produced as Reynolds number increased because of an outward movement of the inboard leading-edge separation at constant angles of attack. This study extends the existing Reynolds number database for supersonic transports operating at off-design conditions.

Description: An interactive computer code, written with a readily available software program, Microsoft Excel (Microsoft Corporation, Redmond, WA) is presented which displays 3 D oblique plots of a conserved scalar distribution downstream of jets mixing with a confined crossflow, for a single row, double rows, or opposed rows of jets with or without flow area convergence and/or a non-uniform crossflow scalar distribution. This project used a previously developed empirical model of jets mixing in a confined crossflow to create an Microsoft Excel spreadsheet that can output the profiles of a conserved scalar for jets injected into a confined crossflow given several input variables. The program uses multiple spreadsheets in a single Microsoft Excel notebook to carry out the modeling. The first sheet contains the main program, controls for the type of problem to be solved, and convergence criteria. The first sheet also provides for input of the specific geometry and flow conditions. The second sheet presents the results calculated with this routine to show the effects on the mixing of varying flow and geometric parameters. Comparisons are also made between results from the version of the empirical correlations implemented in the spreadsheet and the versions originally written in Applesoft BASIC (Apple Computer, Cupertino, CA) in the 1980's.

Description: A computational method for the analysis of longitudinal-mode liquid rocket combustion instability has been developed based on the unsteady, quasi-one-dimensional Euler equations where the combustion process source terms were introduced through the incorporation of a two-zone, linearized representation: (1) A two-parameter collapsed combustion zone at the injector face, and (2) a two-parameter distributed combustion zone based on a Lagrangian treatment of the propellant spray. The unsteady Euler equations in inhomogeneous form retain full hyperbolicity and are integrated implicitly in time using second-order, high-resolution, characteristic-based, flux-differencing spatial discretization with Roe-averaging of the Jacobian matrix. This method was initially validated against an analytical solution for nonreacting, isentropic duct acoustics with specified admittances at the inflow and outflow boundaries. For small amplitude perturbations, numerical predictions for the amplification coefficient and oscillation period were found to compare favorably with predictions from linearized small-disturbance theory as long as the grid exceeded a critical density (100 nodes/wavelength). The numerical methodology was then exercised on a generic combustor configuration using both collapsed and distributed combustion zone models with a short nozzle admittance approximation for the outflow boundary. In these cases, the response parameters were varied to determine stability limits defining resonant coupling onset.

Description: The influence of vortex ring state (VRS) on rotorcraft flight dynamics is investigated, specifically the vertical velocity drop of helicopters and the roll-off of tiltrotors encountering VRS. The available wind tunnel and flight test data for rotors in vortex ring state are reviewed. Test data for axial flow, non-axial flow, two rotors, unsteadiness, and vortex ring state boundaries are described and discussed. Based on the available measured data, a VRS model is developed. The VRS model is a parametric extension of momentum theory for calculation of the mean inflow of a rotor, hence suitable for simple calculations and real-time simulations. This inflow model is primarily defined in terms of the stability boundary of the aircraft motion. Calculations of helicopter response during VRS encounter were performed, and good correlation is shown with the vertical velocity drop measured in flight tests. Calculations of tiltrotor response during VRS encounter were performed, showing the roll-off behavior characteristic of tiltrotors. Hence it is possible, using a model of the mean inflow of an isolated rotor, to explain the basic behavior of both helicopters and tiltrotors in vortex ring state.

Description: In phase II, additional inducer rotations are simulated in order to understand the root cause of the flowliner crack problem. CFD results confirmed that there is a strong unsteady interaction between the backflow regions caused by the LPFTP inducer and secondary flow regions in the bellows cavity through the flowliner slots. It is observed that the swirl on the duct side of the downstream flowliner is stronger than on the duct side of the upstream flowliner. Due to this swirl, there are more significant unsteady flow interactions through the downstream slots than those observed in the upstream slots. Averaged values of the local velocity at the slots were provided to the NESC-ITA flow physics acoustics team to guide them in designing the acoustics experiment. A parametric study was performed to compare the flow field in the flowliner area when one upstream slot and one corresponding downstream slot were enlarged. No significant differences were observed between the flow field obtained from the enlarged slot configuration when compared with the original configuration. More cases must be analyzed with various enlarged slot configurations to generalize this observation. The flow through the A1 test stand and the flow through the orbiter fuel feedline manifold were simulated without the LPFTP. It was observed that incoming flow to the flowliner and inducer was more uniform in the A1 test stand then in the orbiter manifold. Additionally, each engine LPFTP in the orbiter receives significantly different velocity distributions. Because of the differences observed in the computed results, it is not possible for the A1 test stand to represent the three different engine feedlines simultaneously.

Description: An experimental investigation of a flush-mounted, S-duct inlet with large amounts of boundary layer ingestion has been conducted at Reynolds numbers up to full scale. The study was conducted in the NASA Langley Research Center 0.3-Meter Transonic Cryogenic Tunnel. In addition, a supplemental computational study on one of the inlet configurations was conducted using the Navier-Stokes flow solver, OVERFLOW. Tests were conducted at Mach numbers from 0.25 to 0.83, Reynolds numbers (based on aerodynamic interface plane diameter) from 5.1 million to 13.9 million (full-scale value), and inlet mass-flow ratios from 0.29 to 1.22, depending on Mach number. Results of the study indicated that increasing Mach number, increasing boundary layer thickness (relative to inlet height) or ingesting a boundary layer with a distorted profile decreased inlet performance. At Mach numbers above 0.4, increasing inlet airflow increased inlet pressure recovery but also increased distortion. Finally, inlet distortion was found to be relatively insensitive to Reynolds number, but pressure recovery increased slightly with increasing Reynolds number.This CD-ROM supplement contains inlet data including: Boundary layer data, Duct static pressure data, performance-AIP (fan face) data, Photos, Tunnel wall P-PTO data and definitions.

Description: Wind tunnel experiments will continue to be a primary source of validation data for many types of mathematical and computational models in the aerospace industry. The increased emphasis on accuracy of data acquired from these facilities requires understanding of the uncertainty of not only the measurement data but also any correction applied to the data. One of the largest and most critical corrections made to these data is due to wall interference. In an effort to understand the accuracy and suitability of these corrections, a statistical validation process for wall interference correction methods has been developed. This process is based on the use of independent cases which, after correction, are expected to produce the same result. Comparison of these independent cases with respect to the uncertainty in the correction process establishes a domain of applicability based on the capability of the method to provide reasonable corrections with respect to customer accuracy requirements. The statistical validation method was applied to the version of the Transonic Wall Interference Correction System (TWICS) recently implemented in the National Transonic Facility at NASA Langley Research Center. The TWICS code generates corrections for solid and slotted wall interference in the model pitch plane based on boundary pressure measurements. Before validation could be performed on this method, it was necessary to calibrate the ventilated wall boundary condition parameters. Discrimination comparisons are used to determine the most representative of three linear boundary condition models which have historically been used to represent longitudinally slotted test section walls. Of the three linear boundary condition models implemented for ventilated walls, the general slotted wall model was the most representative of the data. The TWICS code using the calibrated general slotted wall model was found to be valid to within the process uncertainty for test section Mach numbers less than or equal to 0.60. The scatter among the mean corrected results of the bodies of revolution validation cases was within one count of drag on a typical transport aircraft configuration for Mach numbers at or below 0.80 and two counts of drag for Mach numbers at or below 0.90.

Description: An experimental investigation of a flush-mounted, S-duct inlet with large amounts of boundary layer ingestion has been conducted at Reynolds numbers up to full scale. The study was conducted in the NASA Langley Research Center 0.3-Meter Transonic Cryogenic Tunnel. In addition, a supplemental computational study on one of the inlet configurations was conducted using the Navier-Stokes flow solver, OVERFLOW. Tests were conducted at Mach numbers from 0.25 to 0.83, Reynolds numbers (based on aerodynamic interface plane diameter) from 5.1 million to 13.9 million (full-scale value), and inlet mass-flow ratios from 0.29 to 1.22, depending on Mach number. Results of the study indicated that increasing Mach number, increasing boundary layer thickness (relative to inlet height) or ingesting a boundary layer with a distorted profile decreased inlet performance. At Mach numbers above 0.4, increasing inlet airflow increased inlet pressure recovery but also increased distortion. Finally, inlet distortion was found to be relatively insensitive to Reynolds number, but pressure recovery increased slightly with increasing Reynolds number.

Description: AirShow is visualization post-processing software for Computational Fluid Dynamics (CFD). Upon reading binary PLOT3D grid and solution files into AirShow, the engineer can quickly see how hundreds of complex 3-D structured blocks are arranged and numbered. Additionally, chosen grid planes can be displayed and colored according to various aerodynamic flow quantities such as Mach number and pressure. The user may interactively rotate and translate the graphical objects using the mouse. The software source code was written in cross-platform Java, C++, and OpenGL, and runs on Unix, Linux, and Windows. The graphical user interface (GUI) was written using Java Swing. Java also provides multiple synchronized threads. The Java Native Interface (JNI) provides a bridge between the Java code and the C++ code where the PLOT3D files are read, the OpenGL graphics are rendered, and numerical calculations are performed. AirShow is easy to learn and simple to use. The source code is available for free from the NASA Technology Transfer and Partnership Office.

Description: This report presents the use of a time-marching three-dimensional compressible Navier-Stokes equation numerical solver with a one-equation turbulence model to simulate the flow fields developed along concave wall surfaces without and with a downstream extension flat wall surface. The 3-D Navier- Stokes numerical solver came from the NASA Glenn-HT code. The one-equation turbulence model was derived from the Spalart and Allmaras model. The computational approach was first calibrated with the computations of the velocity and Reynolds shear stress profiles of a steady flat plate boundary layer flow. The computational approach was then used to simulate developing boundary layer flows along concave wall surfaces without and with a downstream extension wall. The author investigated the computational results of surface friction factors, near surface velocity components, near wall temperatures, and a turbulent shear stress component in terms of turbulence modeling, computational mesh configurations, inlet turbulence level, and time iteration step. The computational results were compared with existing measurements of skin friction factors, velocity components, and shear stresses of the developing boundary layer flows. With a fine computational mesh and a one-equation model, the computational approach could predict accurately the skin friction factors, near surface velocity and temperature, and shear stress within the flows. The computed velocity components and shear stresses also showed the vortices effect on the velocity variations over a concave wall. The computed eddy viscosities at the near wall locations were also compared with the results from a two equation turbulence modeling technique. The inlet turbulence length scale was found to have little effect on the eddy viscosities at locations near the concave wall surface. The eddy viscosities, from the one-equation and two-equation modeling, were comparable at most stream-wise stations. The present one-equation turbulence model is an effective approach for turbulence modeling in the near solid wall surface region of flow over a concave wall.

Description: The Experiment to Characterize Aircraft Volatile Aerosols and Trace Species Emissions (EXCAVATE) took place at NASA Langley Research Center during January 2002. This ground based study was conducted to examine the role of fuel sulfur content on particulate emissions. Size distributions as a function of engine operating conditions were measured in the exhaust plume of a B-757 at four downstream axial locations (1 m, 10 m, 25 m and 35 m). The engine was run on JP-5 with three different sulfur concentrations, 810 ppm, 1050 ppm, 1820 ppm; and was operated over a range of power settings from idle to near-full power. Zalabsky differential-mobility analyzers DMAS), Met One condensation-nuclei counters (CNCs), and a TSI 3022 condensation-particle counter (CPC) were used to measure the size distributions. The total number-count (particle concentration), number-based Emissions Index (EInumber) and mass-based Emissions Index (E1-J increased with fuel sulfur-content and engine pressure ratio (EPR). Count Mean Diameter (Ch4D) also increased with EPR yet remained fairly constant with fuel sulfur-content for a fixed location in the exhaust plume. Also the mode and CMD both increased with distance in the plume.

Description: A free-to-roll study of the low-speed lateral characteristics of the pre-production F/A-18E was conducted in the NASA Langley 12-Foot Low-Speed Tunnel. In developmental flight tests the F/A-18E unexpectedly experienced uncommanded lateral motions in the power approach configuration. The objective of this study was to determine the feasibility of using the free-to-roll technique for the detection of uncommanded lateral motions for the pre-production F/A-18E in the power approach configuration. The data revealed that this technique in conjunction with static data revealed insight into the cause of the lateral motions. The free-to-roll technique identified uncommanded lateral motions at the same angle-of-attack range as experienced in flight tests. The cause of the uncommanded lateral motions was unsteady asymmetric wing stall. The paper also shows that free-to-roll data or static force and moment data alone are not enough to accurately capture the potential for an aircraft to experience uncommanded lateral motion.

Description: This report presents experiments of a less viscous liquid displaced by a more viscous, miscible liquid in a vertical tube of small diameter. The liquids used are silicone oils of various viscosities. The more viscous fluid has a higher density than the less viscous fluid. The experiments showed that the interface between the liquids can be unstable and exhibits an asymmetric, sinuous shape for downward displacement. The Reynolds number is negligibly small in these experiments and the Peclet number for mass transfer is O(50)-0(l0(exp 5)). The influence of gravity is characterized by the dimensionless parameter F=g(Delta)(rho)d(sup 2)/(mu)U where g is the acceleration due to gravity, Delta(rho) is the density difference between the fluids, d is the tube diameter, mu is the viscosity of the more viscous fluid, and U is the displacement speed. Visualization and PIV measurement of the unstable linger are reported. Experiments were also conducted to study immiscible displacements in tubes in order to compare with previously obtained results for miscible fluids.

Description: The third flight of the Hyper-X Research Vehicle at nearly Mach 10 in November 2004 demonstrated the potential for airbreathing propulsion at hypersonic speeds. The flowpath of the scramjet engine used in the Hyper-X Research Vehicle was developed using computational fluid dynamic based design method in combination with a systematic series of ground tests in the NASA HyPulse Shock Tunnel at conditions duplicating Mach 10 flight. The test model was a true scale height and length of the internal flight engine flowpath that was partial width and truncated fore and aft. Tests were conducted using a mixture of hydrogen augmented with two-percent silane. Silane gas served as an ignition aid during the short duration of the pulse facility tests. The engine inflow values of pressure, temperature, and Mach number were parametrically varied during the ground test entries to broaden the database over the expected uncertainty in the flight conditions. A review of the ground test technique and comparisons of the ground test pressures along with selected data from the third flight are presented.

Description: This study investigates the sidewall effect on flow within the mixing duct downstream of a lobed mixer-ejector nozzle. Simulations which model only one half-chute width of the ejector array are compared with those which model one complete quadrant of the nozzle geometry and with available experimental data. These solutions demonstrate the applicability of the half-chute technique to model the flowfield far away from the sidewall and the necessity of a full-quadrant simulation to predict the formation of a low-energy flow region near the sidewall. The quadrant solutions are further examined to determine the cause of this low-energy region, which reduces the amount of mixing and lowers the thrust of the nozzle. Grid resolution and different grid topologies are also examined. Finally, an assessment of the half-chute and quadrant approaches is made to determine the ability of these simulations to provide qualitative and/or quantitative predictions for this type of complex flowfield.

Description: This is the source listing of the computer code SPIRALI which predicts the performance characteristics of incompressible cylindrical and face seals with or without the inclusion of spiral grooves. Performance characteristics include load capacity (for face seals), leakage flow, power requirements and dynamic characteristics in the form of stiffness, damping and apparent mass coefficients in 4 degrees of freedom for cylindrical seals and 3 degrees of freedom for face seals. These performance characteristics are computed as functions of seal and groove geometry, load or film thickness, running and disturbance speeds, fluid viscosity, and boundary pressures.

Description: The proper understanding and modeling of a microcalorimeter s response requires the accurate knowledge of a handful of parameters, such as C, G, alpha, . . . . While a few of these, such 8s the normal state resistance and the total thermal conductance to the heat bath (G) are directly determined from the DC IV characteristics, some others, notoriously the heat capacity (C) and alpha, appear in degenerate combinations in most measurable quantities. The case of a complex microcalorimeter, i.e. one in which the absorber s heat capacity is connected by a finite thermal impedance to the sensor, and subsequently by another thermal impedance to the heat bath, results in an added ambiguity in the determination of the individual C's and G's. In general, the dependence of the microcalorimeter s complex impedance on these parameters varies with frequency. This variation allows us to determine the individual parameters by fitting the prediction of the microcalorimeter model to the impedance data. We describe in this paper our efforts at characterizing the Goddard X-ray microcalorimeters. Using the parameters determined with this method we them compare the pulse shape and noise spectra predicted by the microcalorimeter model to data taken with the same devices.

Description: The objective of this paper is to discuss the testing of a spray cooling experiment designed to be flown on NASA's KC-135 Reduced Gravity Testing Platform. Spray cooling is an example of a thermal management technique that may be utilized in high flux heat acquisition and high thermal energy transport concepts. Many researchers have investigated the utility of spray cooling for the thermal management of devices generating high heat fluxes. However, there has been little research addressing the physics and ultimate performance of spray cooling in a variable gravity environment. An experimental package, consisting of a spray chamber coupled to a fluid delivery loop system, was fabricated for variable gravity flight tests. The spray chamber contains two opposing nozzles spraying on target Indium Tin Oxide (ITO) heaters. These heaters are mounted on glass pedestals, which are part of a sump system to remove unconstrained liquid from the test chamber. Liquid is collected in the sumps and returned to the fluid delivery loop. Thermocouples mounted in and around the pedestals are used to determine both the heat loss through the underside of the IT0 heater and the heat extracted by the spray. A series of flight tests were carried out aboard the KC-135, utilizing the ability of the aircraft to produce various gravity conditions. During the flight tests, for a fixed flow rate, heat input was varied at 20, 30, 50, and 80W with variable gravities of 0.01, 0.16, 0.36, and 1.8g. Flight test data was compared to terrestrial baseline data in addition to analytical and numerical solutions to evaluate the heat transfer in the heater and support structure . There were significant differences observed in the spray cooling performance as a result of variable gravity conditions and heat inputs. In general, the Nussult number at the heater surface was found to increase with decreasing gravity conditions for heat loads greater than 30W.

Description: Studies were performed in order to characterize the thrust augmentation potential of an ejector in a Pulse Detonation Engine application. A 49-mm diameter tube of 0.914-m length was constructed with one open end and one closed end. Ethylene, oxygen, and nitrogen were introduced into the tube at the closed end through the implementation of a fast mixing injector. The tube was completely filled with a stoichiometric mixture containing a one to one molar ratio of nitrogen to oxygen. Ethylene was selected as the fuel due to its detonation sensitivity and the molar ratio of the oxidizer was chosen for heat transfer purposes. Detonations were initiated in the tube through the use of a spark ignition system. The PDE was operated in a multi-cycle mode at frequencies ranging from 20-Hz to 50-Hz. Baseline thrust measurements with no ejector present were performed while operating the engine at various frequencies and compared to theoretical estimates. The baseline values were observed to agree with the theoretical model at low operating frequencies and proved to be increasingly lower than the predicted values as the operating frequency was increased. The baseline thrust measurements were observed to agree within 15 percent of the model for all operating frequencies. A straight 152-mm diameter ejector was installed and thrust augmentation percentages were measured. The length of the ejector was varied while the overlap percentage (percent of the ejector length which overlapped the tube) was maintained at 25 percent for all tests. In addition, the effect of ejector inlet geometry was investigated by comparing results with a straight inlet to those of a 38-mm inlet diameter. The thrust augmentation of the straight inlet ejector proved to be independent of engine operating frequency, augmenting thrust by 40 percent for the 0.914-m length ejector. In contrast, the rounded lip ejector of the same length seemed to be highly dependent on the engine operating frequency. An optimum operating frequency observed with the rounded inlet occurred at an operating frequency of 30-Hz, resulting in thrust augmentation percentages greater than 100 percent. The effect that the engine operating frequency had on thrust augmentation levels attained with an ejector was characterized and optimum performance parameters were established. Insight into the frequency dependent nature of the ejector performance was pursued. Suggestions for future experiments which are needed to fully understand the means in which thrust augmentation is achieved in a PDE-ejector configuration were noted.

Description: Thrust augmentation results of an ongoing study of pulse detonation engine driven ejectors are presented and discussed. The experiments were conducted using a pulse detonation engine (PDE) setup with various ejector configurations. The PDE used in these experiments utilizes ethylene (C2H4) as the fuel, and an equi-molar mixture of oxygen and nitrogen as the oxidizer at an equivalence ratio of one. High fidelity thrust measurements were made using an integrated spring damper system. The baseline thrust of the PDE engine was first measured and agrees with experimental and modeling results found in the literature. Thrust augmentation measurements were then made for constant diameter ejectors. The parameter space for the study included ejector length, PDE tube exit to ejector tube inlet overlap distance, and straight versus rounded ejector inlets. The relationship between the thrust augmentation results and various physical phenomena is described. To further understand the flow dynamics, shadow graph images of the exiting shock wave front from the PDE were also made. For the studied parameter space, the results showed a maximum augmentation of 40%. Further increase in augmentation is possible if the geometry of the ejector is tailored, a topic currently studied by numerous groups in the field.

Description: The next generation of aircraft engines will need to pass stricter efficiency and emission tests. NASA's Ultra-Efficient Engine Technology (UEET) program has set an ambitious goal of 70% reduction of NO(x) emissions and a 15% increase in fuel efficiency of aircraft engines. We will demonstrate the state-of-the-art combustion tools developed a t Stanford's Center for Turbulence Research (CTR) as part of this program. In the last decade, CTR has spear-headed a multi-physics-based combustion modeling program. Key technologies have been transferred to the aerospace industry and are currently being used for engine simulations. In this demo, we will showcase the next-generation combustion modeling tools that integrate a very high level of detailed physics into advanced flow simulation codes. Combustor flows involve multi-phase physics with liquid fuel jet breakup, evaporation, and eventual combustion. Individual components of the simulation are verified against complex test cases and show excellent agreement with experimental data.

Description: Determinations of the phase formation sequence, the crystal structures and the thermodynamic properties of materials at high temperatures are difficult because of contamination from the sample container and environment. Containerless processing techniques, such as electrostatic (ESL), electromagnetic (EML), aerodynamic, and acoustic levitation, are most suitable these studies. An adaptation of ESL for in-situ structural studies of a wide range of materials, including metals, semiconductors, insulators using high energy (125 keV) synchrotron x-rays is described here. This beamline ESL (BESL) allows the in-situ determination of the atomic structures of equilibrium solid and liquid phases, including undercooled liquids, as well as real-time studies of solid-solid and liquid-solid phase transformations. The use of image plate (MAR345) or GE-Angio detectors enables fast (30 ms - 1s) acquisition of complete diffraction patterns over a wide q-range (4 - 140/mm). The wide temperature range (300 - 2500 K), containerless processing under high vacuum (10(exp -7) - 10(exp -8) torr), and fast data acquisition, make BESL particularly suitable for phase diagram studies of high temperature materials. An additional, critically important, feature of BESL is the ability to also make simultaneous measurement of a host of thermo-physical properties, including the specific heat, enthalpy of transformation, solidus and liquidus temperatures, density, viscosity, and surface tension; all on the same sample and simultaneous with the structural measurements.

Description: Static and dynamic magnetic fields have been used to control convection in many materials processing applications. In most of the applications, convection control (damping or enhancement) is achieved through the Lorentz force that can be tailored to counteract/assist dominant system flows. This technique has been successfully applied to liquids that are electrically conducting, such as high temperature melts of semiconductors, metals and alloys, etc. In liquids with low electrical conductivity such as ionic solutions of salts in water, the Lorentz force is weak and hence not very effective and alternate ways of flow control are necessary. If the salt in solution is paramagnetic then the variation of magnetic susceptibility with temperature and/or concentration can be used for flow control. For thermal buoyancy driven flows this can be accomplished in a temperature range below the Curie point of the salt. The magnetic force is proportional to the magnetic susceptibility and the product of the magnetic field and its gradient. By suitably positioning the experiment cell in the magnet, system flows can be assisted or countered, as desired. A similar approach can be extended to diamagnetic substances and fluids but the required magnetic force is considerably larger than that required for paramagnetic substances. The presentation will provide an overview of work to date on a NASA fluid physics sponsored project that aims to test the hypothesis of convective flow control using strong magnetic fields in protein crystal growth. The objective is to understand the nature of the various forces that come into play, delineate causative factors for fluid flow and to quantify them through experiments, analysis, and numerical modeling. The seminar will report specifically on the experimental results using paramagnetic salts and solutions in magnetic fields and compare them to analytical predictions. Applications of the concept to protein crystallization studies will be discussed. The use of strong magnetic fields for terrestrially simulating variable gravity environments and applications supporting the NASA Exploration Initiative will also be briefly discussed.

Description: The James Webb Space Telescope will include an optical bench known as the integrated science instrument module (ISIM). Candidate structural materials for the ISIM must have low density, high stiffness, high thermal conductivity, and low thermal expansion coefficient at the operating temperature of 30 Kelvin. The specific heat is also important in modeling the on-orbit cooldown. We built two different systems for measuring the thermal conductivity and specific heat of samples between 4 Kelvin and 290 Kelvin. Both experiments were carefully designed to minimize potential errors due to radiative heat transfer. We chose the cooling system and instrumentation to allow long-term unattended operation. Software was developed to automate each experiment. It used an algorithm designed to ensure that each system was in thermal equilibrium before a measurement was taken. We describe the two experiments and present the data.

Description: The paper experimentally and theoretically studies the effects of periodic unsteady wake flow and aerodynamic characteristics on boundary layer development, separation and re-attachment along the suction surface of a low pressure turbine blade. The experiments were carried out at Reynolds number of 110,000 (based on suction surface length and exit velocity). For one steady and two different unsteady inlet flow conditions with the corresponding passing frequencies, intermittency behaviors were experimentally and theoretically investigated. The current investigation attempts to extend the intermittency unsteady boundary layer transition model developed in previously to the LPT cases, where separation occurs on the suction surface at a low Reynolds number. The results of the unsteady boundary layer measurements and the intermittency analysis were presented in the ensemble-averaged and contour plot forms. The analysis of the boundary layer experimental data with the flow separation, confirms the universal character of the relative intermittency function which is described by a Gausssian function.

Description: An analysis of the optimal control by blowing and suction in order to generate stream- wise velocity streaks is presented. The problem is examined using an iterative process that employs the Parabolized Stability Equations for an incompressible uid along with its adjoint equations. In particular, distributions of blowing and suction are computed for both the normal and tangential velocity perturbations for various choices of parameters.

Description: This report's contents focus on making use of NASA Glenn on-site computational facilities,to develop, validate, and apply models for use in advanced 3D Navier-Stokes Computational Fluid Dynamics (CFD) codes to enhance the capability to compute heat transfer and losses in turbomachiney.

Description: An advantage of using a Radioisotope Power System (RPS) for deep space or planetary surface missions is the readily available waste heat, which can be used for a number of beneficial purposes including: maintaining electronic components within a controlled temperature range, warming propulsion tanks and mobility actuators, and maintaining liquid propellants above their freezing temperature. Previous missions using Radioisotope Thermoelectric Generators (RTGs) dissipated large quantities of waste heat due to the low efficiency of the thermoelectric conversion technology. The next generation RPSs, such as the 110-Watt Stirling Radioisotope Generator (SRG110) will have higher conversion efficiencies, thereby rejecting less waste heat at a lower temperature and may require alternate approaches to transferring waste heat to the spacecraft. RTGs, with efficiencies of ~6 to 7 percent, reject their waste heat at the relatively high heat rejection temperature of 200 C. This is an advantage when rejecting heat to space; however, transferring heat to the internal spacecraft components requires a large and heavy radiator heat exchanger. At the same time, sensitive spacecraft instruments must be shielded from the thermal radiation of the RTG. The SRG110, with an efficiency around 22 percent and 50 C nominal housing surface temperature, can readily transfer the available waste heat directly via heat pipes, thermal straps, or fluid loops. The lower temperatures associated with the SRG110 avoid the chances of overheating other scientific components, eliminating the need for thermal shields. This provides the spacecraft designers more flexibility when locating the generator for a specific mission. A common misconception with high-efficiency systems is that there is not enough waste heat for spacecraft thermal management. This paper will dispel this misconception and investigate the use of a high-efficiency SRG110 for spacecraft thermal management and outline potential methods of waste heat utilization in several conceptual missions (Lunar Rover, Mars Rover, and Titan Lander). The advantages associated with the SRG110 as they relate to ease of assembly, less complex interfaces, and overall mass savings for a spacecraft will be highlighted.

Description: Background on thermal effects on cavitation and numerical framework. Validation of numerical model for cryogens in CRUNCH CFD(R). Comparison with subscale test data hord (1973). Simulations of liquid hydrogen inducer at various flow coefficients. 120% of design, Design, and 80% of Design flow rate. Detailed comparison of flow profiles. Sensitivity of backflow to turbulent viscosity noted. Conclusion.

Description: Experiments were conducted to study the effects of enhanced surface structures on heat flux using spray cooling. The surface enhancements consisted of cubic pin fins machined on the top surface of copper heater blocks. The structure height, pitch, and width were parametrically vaned. Each copper block had a projected cross-sectional area of 2.0 sq cm. Measurements were also obtained on a heater block with a flat surface for baseline comparison purposes. A 2 x 2 nozzle array was used with PF-5060 as the working fluid. Thermal performance data were obtained under nominally degassed (chamber pressure of 41.4 kPa) and gassy conditions (chamber with N2 gas at 100.7 kPa) with a bulk fluid temperature of 20.5 C. Results for both the degassed and gassy cases show that structure width and separation distance have a dominant effect upon the heat transfer for the size ranges used. Cubic pin fin height had little impact upon heat flux. The maximum critical heat flux (CHF) attained for any of the surfaces was 121 W/sq cm, giving an enhancement of 51% relative to the flat surface case under nominally degassed conditions. The gassy case had a maximum CHF of 149 W/sq cm, giving an enhancement of 38% relative to the flat surface case.

Description: A new data analysis technique for the identification of static and dynamic aerodynamic stability coefficients from wind tunnel test video data is presented. This new technique was applied to video data obtained during a parachute wind tunnel test program conducted in support of the Mars Exploration Rover Mission. Total angle-of-attack data obtained from video images were used to determine the static pitching moment curve of the parachute. During the original wind tunnel test program the static pitching moment curve had been determined by forcing the parachute to a specific total angle-of -attack and measuring the forces generated. It is shown with the new technique that this parachute, when free to rotate, trims at an angle-of-attack two degrees lower than was measured during the forced-angle tests. An attempt was also made to extract pitch damping information from the video data. Results suggest that the parachute is dynamically unstable at the static trim point and tends to become dynamically stable away from the trim point. These trends are in agreement with limit-cycle-like behavior observed in the video. However, the chaotic motion of the parachute produced results with large uncertainty bands.

Description: The ability to accurately model details of inlet back flow for inducers operating a t low-flow, off-design conditions is evaluated. A sub-scale version of a three-bladed liquid hydrogen inducer tested in water with detailed velocity and pressure measurements is used as a numerical test bed. Under low-flow, off-design conditions the length of the separation zone as well as the swirl velocity magnitude was under predicted with a standard k-E model. When the turbulent viscosity coefficient was reduced good comparison was obtained a t all the flow conditions examined with both the magnitude and shape of the profile matching well with the experimental data taken half a diameter upstream of the leading edge. The velocity profiles and incidence angles a t the leading edge itself were less sensitive to the back flow length predictions indicating that single-phase performance predictions may be well predicted even if the details of flow separation modeled are incorrect. However, for cavitating flow situations the prediction of the correct swirl in the back flow and the pressure depression in the core becomes critical since it leads to vapor formation. The simulations have been performed using the CRUNCH CFD(Registered Trademark) code that has a generalized multi-element unstructured framework and a n advanced multi-phase formulation for cryogenic fluids. The framework has been validated rigorously for predictions of temperature and pressure depression in cryogenic fluid cavities and has also been shown to predict the cavitation breakdown point for inducers a t design conditions.

Description: This paper describes numerical implementation of a newly developed hybrid model, T-blob/T-TAB, into an existing computational fluid dynamics (CFD) program for primary and secondary breakup simulation of liquid jet atomization. This model extend two widely used models, the Kelvin-Helmholtz (KH) instability of Reitz (blob model) and the Taylor-Analogy-Breakup (TAB) secondary droplet breakup by O'Rourke and Amsden to include turbulence effects. In the primary breakup model, the level of the turbulence effect on the liquid breakup depends on the characteristic scales and the initial flow conditions. For the secondary breakup, an additional turbulence force acted on parent drops is modeled and integrated into the TAB governing equation. Several assessment studies are presented and the results indicate that the existing KH and TAB models tend to under-predict the product drop size and spray angle, while the current model provides superior results when compared with the measured data.

Description: High-energy rotating machines, powering liquid propellant rocket engines, are subject to various sources of high and low cycle fatigue generated by unsteady flow phenomena. Given the tremendous need for reliability in a sustainable space exploration program, a fundamental change in the design methodology for engine components is required for both launch and space based systems. A design optimization system based on neural-networks has been applied and demonstrated in the redesign of the Space Shuttle Main Engine (SSME) Low Pressure Oxidizer Turbo Pump (LPOTP) turbine nozzle. One objective of the redesign effort was to increase airfoil thickness and thus increase its strength while at the same time detuning the vane natural frequency modes from the vortex shedding frequency. The second objective was to reduce the vortex shedding amplitude. The third objective was to maintain this low shedding amplitude even in the presence of large manufacturing tolerances. All of these objectives were achieved without generating any detrimental effects on the downstream flow through the turbine, and without introducing any penalty in performance. The airfoil redesign and preliminary assessment was performed in the Exploration Technology Directorate at NASA ARC. Boeing/Rocketdyne and NASA MSFC independently performed final CFD assessments of the design. Four different CFD codes were used in this process. They include WIL DCA T/CORSAIR (NASA), FLUENT (commercial), TIDAL (Boeing Rocketdyne) and, a new family (AardvarWPhantom) of CFD analysis codes developed at NASA MSFC employing LOX fluid properties and a Generalized Equation Set formulation. Extensive aerodynamic performance analysis and stress analysis carried out at Boeing Rocketdyne and NASA MSFC indicate that the redesign objectives have been fully met. The paper presents the results of the assessment analysis and discusses the future potential of robust optimal design for rocket engine components.

Description: Pulse detonation rocket engines (PDREs) have generated research interest in recent years as a chemical propulsion system potentially offering improved performance and reduced complexity compared to conventional rocket engines. The detonative mode of combustion employed by these devices offers a thermodynamic advantage over the constant-pressure deflagrative combustion mode used in conventional rocket engines and gas turbines. However, while this theoretical advantage has spurred considerable interest in building PDRE devices, the unsteady blowdown process intrinsic to the PDRE has made realistic estimates of the actual propulsive performance problematic. The recent review article by Kailasanath highlights some of the progress that has been made in comparing the available experimental measurements with analytical and numerical models. In recent work by the author, a quasi-one-dimensional, finite rate chemistry CFD model was utilized to study the gasdynamics and performance characteristics of PDREs over a range of blowdown pressure ratios from 1-1000. Models of this type are computationally inexpensive, and enable first-order parametric studies of the effect of several nozzle and extension geometries on PDRE performance over a wide range of conditions. However, the quasi-one-dimensional approach is limited in that it cannot properly capture the multidimensional blast wave and flow expansion downstream of the PDRE, nor can it resolve nozzle flow separation if present. Moreover, the previous work was limited to single-pulse calculations. In this paper, an axisymmetric finite rate chemistry model is described and utilized to study these issues in greater detail. Example Mach number contour plots showing the multidimensional blast wave and nozzle exhaust plume are shown. The performance results are compared with the quasi-one-dimensional results from the previous paper. Both Euler and Navier-Stokes solutions are calculated in order to determine the effect of viscous effects in the nozzle flowfield. Additionally, comparisons of the model results to performance data from CalTech, as well as experimental flowfield measurements from Stanford University, are also reported.