ALBERT

Description: A heat source such as a magnetic induction/eddy current generator remotely heats a region of a surface of a test structure to a desired depth. For example, the frequency of the heating source can be varied to heat to the desired depth. A thermal sensor senses temperature changes in the heated region as a function of time. A computer compares these sensed temperature changes with calibration standards of a similar sample having known disbond and/or inclusion geography(ies) to analyze the test structure. A plurality of sensors can be arranged linearly to sense vector heat flow.

Description: Flat or curved micro heat pipe panels are fabricated by arranging essentially parallel filaments in the shape of the desired panel. The configuration of the filaments corresponds to the desired configuration of the tubes that will constitute the heat pipes. A thermally conductive material is then deposited on and around the filaments to fill in the desired shape of the panel. The filaments are then removed, leaving tubular passageways of the desired configuration and surface texture in the material. The tubes are then filled with a working fluid and sealed. Composite micro heat pipe laminates are formed by layering individual micro heat pipe panels and bonding them to each other to form a single structure. The layering sequence of the micro heat pipe panels can be tailored to transport heat preferentially in specific directions as desired for a particular application.

Description: Numerical solutions for hypersonic flows of carbon-dioxide and air around a 70-deg sphere-cone have been computed using an axisymmetric non-equilibrium Navier-Stokes solver. Freestream flow conditions for these computations were equivalent to those obtained in an experimental blunt-body heat-transfer study conducted in a high-enthalpy, hypervelocity expansion tube. Comparisons have been made between the computed and measured surface heat-transfer rates on the forebody and afterbody of the sphere-cone and on the sting which supported the test model. Computed forebody heating rates were within the estimated experimental uncertainties of 10% on the forebody and 15% in the wake except for within the recirculating flow region of the wake.

Description: Heat Pipes were originally developed by NASA and the Los Alamos Scientific Laboratory during the 1960s to dissipate excessive heat build- up in critical areas of spacecraft and maintain even temperatures of satellites. Heat pipes are tubular devices where a working fluid alternately evaporates and condenses, transferring heat from one region of the tube to another. KONA Corporation refined and applied the same technology to solve complex heating requirements of hot runner systems in injection molds. KONA Hot Runner Systems are used throughout the plastics industry for products ranging in size from tiny medical devices to large single cavity automobile bumpers and instrument panels.

Description: We have developed a calculational model that treats all the components of an orifice pulse tube cooler. We base our analysis on 1-dimensional thermodynamic equations for the regenerator and we assume that all mass flows, pressure oscillations and temperature oscillations are small and sinusoidal. Non-linear pressure drop effects are included in the regenerator to account for finite pressure amplitude effects. The resulting mass flows and pressures are matched at the boundaries with the other components of the cooler: compressor, aftercooler, cold heat exchanger, pulse tube, hot heat exchanger, orifice and reservoir. The results of the calculation are oscillating pressures, mass flows and enthalpy flows in the main components of the cooler. By comparing with the calculations of other available models, we show that our model is very similar to REGEN 3 from NIST and DeltaE from Los Alamos National Lab. Our model is much easier to use than other available models because of its simple graphical interface and the fact that no guesses are required for the operating pressures or mass flows. In addition, the model only requires a few minutes of running time allowing many parameters to be optimized in a reasonable time. A version of the model is available for use over the World Wide Web at http://irtek.arc.nasa.gov. Future enhancements include adding a bypass orifice and including second order terms in steady mass streaming and steady heat transfer. A two-dimensional anelastic approximation of the fluid equations will be used as the basis for the latter analysis. Preliminary results are given in dimensionless numbers appropriate for oscillating compressible flows. The model shows how transverse heat transfer reduces enthalpy flow, particularly for small pulse tubes. The model also clearly shows mass recirculation in the open tube on the order of the tube length. They result from the higher order Reynolds stresses. An interesting result of the linearized approach is that the steady mass streaming does not affect the enthalpy flow at second order. The major effect of recirculating mass streaming is to increase transverse temperature gradients, which leads to higher entropy production and reduced efficiency.

Description: The Harmonic Point Source (HPS) is the simplest possible form of 3D disturbance yet even this case is not fully understood. One of the characteristics of HPS generated Tolmien Schlichting (TS) waves is that the maximum rms amplitude downstream of the source occurs away from the centerline. As the waves fan out in the spanwise direction they can encounter Klebanoff modes (weak streamwise vortices) which may be present as background disturbances within the layer. The vortices originate near the leading edge and they appear to be caused by amplification of almost immeasurably small nonuniformities in the free stream. The vortices can be steady but they often appear as very low frequency background unsteadiness (as observed by Klebanoff). The vortices locally distort the mean flow and therefore (by definition) they are nonlinear phenomena e.g. the growth rate of the TS waves is altered. Kendall (private communication) has demonstrated total suppression of TS waves by deliberately introducing strong vortices generated by a delta wing. For the very weak vortices that are often encountered in experiments, their effect is characterized by a local reduction in the rms TS magnitude. The laminar wake behind a fine wire is used to generate Klebanoff modes and the interactions with HPS generated TS waves are explored in a controlled manner. It is shown that a proportion of the reduction in the observed TS magnitude can be attributed to washout of the phase-averaged signals owing to phase jitter.

Description: Calculations have been carried out on the adiabatic laminar boundary layer developing on the surface of a two-dimensional supersonic nozzle, consisting of contoured nozzle blocks and flat sidewalls. Two- and three-dimensional Navier-Stokes codes, as well as two-dimensional boundary-layer codes have been employed. Thee codes have been adapted to the characteristics of a specific wind tunnel nozzle, so that their numerical results could be directly compared with experimental data obtained in the same nozzle. Such comparisons have been made for the boundary-layer growth on the contoured nozzle, and for the boundary-layer growth, surface streamlines and surface shear on the sidewalls. The three-dimensional Navier-Stokes code was found to be the only one to correctly predict the mean boundary-layer flow on both the sidewalls and the contoured nozzle. Theory and experiment both indicate that the sidewall flow is highly three-dimensional, with non-uniform shear, comer vortices and a boundary layer strongly distorted by cross flows induced by lateral pressure gradients.

Description: RANS-MP, a new implementation of a single-grid Navier-Stokes solver using the diagonalized Beam-Warming approximate-factorization scheme, is presented. This first release of the completely rewritten solver employs the following optimizations: (1) Bi-directional multi-partition method for the ADI solver part; this improves granularity and load balance; (2) Improved cache usage through elimination of non-unit-stride array access (possible in part due to multi-partitioning); (3) Preprocessing of communicating boundary conditions to streamline logic during time stepping; (4) Truly parallel, high-performance I/O using the newly-developed MPI-IO library; (5) Elimination of large amounts of redundant operations through efficient use of workspace. Results of some realistic wing computations on the IBM SP2 computer will be presented. We will demonstrate that excellent absolute performance and scalability are obtained with RANS-MP, even for relatively small grid sizes. Besides high performance, an outstanding feature of RANS-MP is its true portability, due to the use of the portable message passing and I/O libraries MPI and MPI-IO.

Description: Considerable progress in understanding nonlinear phenomena in both unbounded and wallbounded shear flow transition has been made through the use of a combination of high- Reynolds-number asymptotic and numerical methods. The objective of this continuing work is to fully understand the nonlinear dynamics so that ultimately (1) an effective means of mixing and transition control can be developed and (2) the source terms in the aeroacoustic noise problem can be modeled more accurately. Two important aspects of the work are that (1) the disturbances evolve from strictly linear instability waves on weakly nonparallel mean flows so that the proper upstream conditions are applied in the nonlinear or wave-interaction streamwise region and (2) the asymptotic formulations lead to parabolic problems so that the question of proper out-flow boundary conditions--still a research issue for direct numerical simulations of convectively unstable shear flows--does not arise. Composite expansion techniques are used to obtain solutions that account for both mean-flow-evolution and nonlinear effects. A previously derived theory for the amplitude evolution of a two-dimensional instability wave in an incompressible mixing layer (which is in quantitative agreement with available experimental data for the first nonlinear saturation stage for a plane-jet shear layer, a circular-jet shear layer, and a mixing layer behind a splitter plate) have been extended to include a wave-interaction stage with a three-dimensional subharmonic. The ultimate wave interaction effects can either give rise to explosive growth or an equilibrium solution, both of which are intimately associated with the nonlinear self-interaction of the three dimensional component. The extended theory is being evaluated numerically. In contrast to the mixing-layer situation, earlier comparisons of theoretical predictions based on asymptotic methods and experiments in wall-bounded shear-flow transition have been somewhat lacking in one aspect or another. The current work strongly suggests that the main weakness is the underlying asymptotic representation of the linear "part" of the problem and not the explicit modeling of the nonlinear/wave-interaction effects. Consequently, the long-wave-length/high-Reynolds-number asymptotic limit for the Blasius boundary-layer stability problem was reexamined, and a new dispersion relationship for the instability waves that is uniformly valid for both the upper- and lower branch regions to the required order of approximation was obtained. A comparison with numerical results, obtained by solving the Orr-Sommerfeld stability problem, shows that the asymptotic formula provides surprisingly good results, even for values of the frequency parameter usually encountered in experimental investigations. This is particularly evident in the dynamically important upper-branch region, where much of the nonlinear interactions in transition experiments are believed to take place. The result is important in that it can be used to greatly improve the accuracy of weakly nonlinear critical-layer-based theories, and a consistent nonlinear theory is currently under evaluation.

Description: An experimental investigation was carried out into the replacement of air in the driven tube of a reflected shock tunnel by an N2O/N2 mixture in order to increase the test time. The incident shock velocities were between 2 and 3 km/sec. Test times were estimated from light emission histories in the driven tube (at distance of L/D = 46.5 from the main diaphragm) and in the nozzle at an area ratio of 27.9 and from pressure histories just upstream of the nozzle entrance (at L/D = 54). The test times estimated from the light emission histories in the driven tube showed that consistent increases of 60-100% were obtained upon substituting N2O/N2 for air in the driven tube. These increases were in very good agreement with theoretical estimates. The test times estimated from the light emission histories in the nozzle or pressure histories at the nozzle inlet showed significant improvements with N2O/N2 only for cases where the facility was operated at substantially overtailored conditions. It is believed that this is due to the greater stability of the driver-driven interface at overtailored operating conditions. At overtailored operating conditions, test times increases of 60-100% with N2O/N2 were observed with all three diagnostic techniques. These increases were in reasonable agreement with theoretical estimates.

Description: In this work we consider solving matrices which arise from the discretization of advection-diffusion field equations on arbitrary triangulated domains using stabilized numerical methods. The talk will discuss several candidate matrix preconditioning algorithms based on the 2 x 2 block factorization induced by an apriori partitioning of the triangulated domain. Application of the 2 x 2 block preconditioner requires the formation and inversion of the Schur complement submatrix. We consider several strategies for simplifying this task: incomplete Schur complement factorizations, drop tolerance element filling, Schur complement probing, and localized Schur complement inversion. Numerical results will be shown comparing performance and efficiency of these approximations. The matrix preconditioner has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The remainder of the talk will show numerous examples in CFD to demonstrate the efficiency and robustness of the techniques.

Description: To provide the proper test conditions for various applications, the characteristics of arcjet test facilities should be known. To determine the operational characteristics of an arcjet facility, emission spectroscopy measurements have been performed and analyzed in a previous study. As an extension of the study, radiation from the free stream and the shock layer of 15 cm diameter blunt-body test articles is measured under different test conditions in the NASA Ames 20 MW Arcjet Facility. The test gas is a mixture of argon and air or N2 and argon. To capture the spatially resolved emission spectra along the stagnation streamline of the blunt-body, a CCD camera (1024 x 256 array) with .275 m Czerny-Turner spectrograph is used. The optical system is calibrated in situ using tungsten and deuterium radiation sources. The emission measurements supplement the existing data base and will help provide better understanding of the equilibration processes in the shock layer region.

Description: A Navier-Stokes flow solver, OVERFLOW, has been developed by researchers at NASA Ames Research Center to use overset (Chimera) grids to simulate the flow about complex aerodynamic shapes. Primary customers of the OVERFLOW flow solver and related software include McDonnell Douglas and Boeing, as well as the NASA Focused Programs for Advanced Subsonic Technology (AST) and High Speed Research (HSR). Code development has focused on customer issues, including improving code performance, ability to run on workstation clusters and the NAS SP2, and direct interaction with industry on accuracy assessment and validation. Significant interaction with NAS has produced a capability tailored to the Ames computing environment, and code contributions have come from a wide range of sources, both within and outside Ames.

Description: The stratospheric observatory for infrared astronomy (SOFIA) is a 2.5 meter aperture Cassegrain telescope with a Nasmyth focus that will be housed in an open cavity in the Boeing 747-SP aircraft and operated at altitudes around 41,000 feet for infrared (IR) viewing of celestial events of astronomical nature. At these altitudes the IR viewing capability of SOFIA far exceeds that of any ground based system. To minimize IR transmission losses, SOFIA will operate with an open cavity. Such an open cavity during flight creates several challenging aerodynamic and aeroacoustic design problems. Foremost of these are: the shear layer over the cavity may cause unwanted resonance if the cavity is untreated; this might give rise to excessive sound pressure levels (SPL) in the cavity and thus affect the unsteady loads on the telescope; the unsteady flow within the cavity produces large dynamic loads and moments that will impact the pointing accuracy of the telescope; the open cavity and the shear layer control devices produce additional drag that will affect directly the time of flight of the mission; the aft location of the cavity down stream of port wing will affect the the flow on the aircraft control surfaces and thus the stability of the aircraft. Also, the highly turbulent shear layer over the cavity and the temperature gradients and 'hot spots' within the cavity can produce a wave front error of the image when it reaches the focal plane of the recorder.

Description: This investigation is a continuation of a previous study on nonequilibrium convective heat transfer to a blunt body. In the previous study, for relatively high Reynolds number flows, it was found that: nonequilibrium convective heat transfer to a blunt body is not strongly dependent on freestream parameters, provided that the thermochemical equilibrium is reached at the edge of boundary layer; and successful testing of convective heat transfer in an arc-jet environment is possible by duplicating the surface pressure and total enthalpy. The nonequilibrium convective heat transfer computations are validated against the results of Fay and Riddell/Goulard theory. Present work investigates low Reynolds number conditions which are typical in an actual arc-jet flow environment. One expects that there will be departures from the Fay and Riddell/Goulard result since certain assumptions of the classical theory are not satisfied. These departures are of interest because the Fay and Riddell/Goulard formulas are extensively used in arc-jet testing (e.g., to determine the enthalpy of the flow and the catalytic efficiency of heat shield materials). For practical sizes of test materials, density of the test flow (and Reynolds number) in an arc-jet is such that thermochemical equilibrium may not be reached at the edge of boundary layer. For blunt body flows of nitrogen and air, computations will be presented to show the effects of thermochemical nonequilibrium at the boundary layer edge on nonequilibrium heat transfer.

Description: We are interested in the study, via numerical simulations, of active vortex generators. Vortex generators may be used to modify the inner part of the boundary layer or to control separation thus enhancing the performance and maneuverability of aerodynamic configurations. We consider generators that consist of a surface cavity elongated in the streamwise direction and partially covered with a moving lid that at rest lies flush with the boundary. Streamwise voracity is generated and ejected due to the oscillatory motion of the lid. The present simulations c Implement relevant experimental investigations of active vortex generators that have been conducted at NASA Ames Research Center and Stanford University. Jacobson and Reynolds used a piezoelectric device in water, allowing for small amplitude high frequency oscillations. They placed the lid asymmetrically on the cavity and observed a strong outward velocity at the small gap of the cavity. Saddoughi used a larger mechanically driven device in air to investigate this flow and observed a jet emerging from the wide gap of the configuration, contrary to the findings of Jacobson and Reynolds. More recently, Lachowiez and Wlezien are investigating the flow generated by an electro-mechanically driven lid to be used for assertion control in aerodynamic applications. We are simulating the flows generated by these devices and we are conducting a parametric study that would help us elucidate the physical mechanisms present in the flow. Conventional computational schemes encounter difficulties when simulating flows around complex configurations undergoing arbitrary motions. Here we present a formulation that achieves this task on a purely Lagrangian frame by extending the formulation presented by Koumoutsakos, Leonard and Pepin. The viscous effects are taken into account by modifying the strength of the particles, whereas fast multipole schemes employing hundreds of thousands ol'particle's allow for high resolution simulations. We shall present simulation results of an oscillating plate at various Reynolds numbers and Strouhal frequencies. Estimates of the forces needed to drive the devices will also be presented.

Description: Recent progress in incompressible Navier-Stokes solution methods will be presented. Discussions are focused on the methods designed for complex geometry applications in three dimensions, and thus are limited to primitive variable formulation. Both steady- and unsteady-solution algorithms and their salient features are discussed. A summary of our recent progress in flow solver development is given followed by numerical studies of a few example problems of our current interest. Solvers discussed here are based on structured-grid approach using finite-difference or finite-volume frame work. This short course will be delivered in three one-our lectures. The material in the course are collected from the work performed by the Incompressible Navier-Stokes group at NASA Ames Research Center over the past several years, and can be found in our publications widely disseminated in the US and abroad. This short course is sponsored by AGARD Consultant and Exchange Program under Support Project P-110.

Description: We are interested in the study, via numerical simulations, of active vortex generators. Vortex generators may be used to modify the inner part of the boundary layer or to control separation thus enhancing the performance and maneuverability of aerodynamic configurations. we consider generators that consist of a surface cavity elongated in the streamwise direction and partially covered with a moving lid that at rest lies flush with the boundary. Streamwise vorticity is generated and ejected due to the oscillatory motion of the lid. The present simulations complement relevant experimental investigations of active vortex generators that have been conducted at NASA Ames Research Center and Stanford University used a piezoelectric device in water, allowing for small amplitude high frequency oscillations. They placed the lid asymmetrically on the cavity and observed a strong outward velocity at the small gap of the cavity. Saddoughi used a larger mechanically driven device in air to investigate this flow and observed a jet emerging from the wide gap of the configuration, contrary to the findings of Jacobson and Reynolds We are simulating the flows generated by these devices and we are conducting a parametric study that would help us elucidate the physical mechanisms present in the flow. Conventional computational schemes encounter difficulties when simulating flows around complex configurations undergoing arbitrary motions. Here we present a formulation that achieves this task on a purely Lagrangian frame by extending the formulation presented by Koumoutsakos, Leonard and Pepin (1994). The viscous effects are taken into account by modifying the strength Of the particles, whereas fast multipole schemes employing hundreds of thousands of particles allow for high resolution simulations. We shall present simulation results of an oscillating plate at various Reynolds numbers and Strouhal frequencies.

Description: The application of nonlinear dynamics to improve the understanding of numerical uncertainties in computational fluid dynamics (CFD) is reviewed. Elementary examples in the use of dynamics to explain the nonlinear phenomena and spurious behavior that occur in numerics are given. The role of dynamics in the understanding of long time behavior of numerical integrations and the nonlinear stability, convergence, and reliability of using time-marching approaches for obtaining steady-state numerical solutions in CFD is explained. The study is complemented with examples of spurious behavior observed in CFD computations.

Description: Over the past 20 years, the use of oil-film interferometry to measure the skin friction coefficient (C(sub f) = tau/q where tau is the surface shear stress and q is the dynamic pressure) has increased. Different forms of this oil-film technique with various levels of accuracy and ease of use have been successfully applied in a wide range of flows. The method's popularity is growing due to its relative ease of implementation and minimal intrusiveness as well as an increased demand for C(sub f) measurements. Nonetheless, the accuracy of these methods has not been rigorously addressed to date. Most researchers have simply shown that the skin-friction measurements made using these techniques compare favorably with other measurements and theory, most of which are only accurate to within 5-20%. The use of skin-friction data in the design of commercial aircraft, whose drag at cruise is 50% skin-friction drag, and in the validation of computational fluid dynamics programs warrants better uncertainty estimates. Additional information is contained in the original extended abstract.

Description: An empty test section flow calibration of the refurbished NASA Ames 12-Foot Pressure Tunnel was recently completed. Distributions of total pressure, dynamic pressure, Mach number, flow angularity temperature, and turbulence are presented along with results obtained prior to facility demolition. Axial static pressure distributions along tunnel centerline are also compared. Test section model support geometric configurations will be presented along with a discussion of the issues involved with different model mounting schemes.

Description: A set of compressible flow relations for a thermally perfect, calorically imperfect gas are derived for a value of c(sub p) (specific heat at constant pressure) expressed as a polynomial function of temperature and developed into a computer program, referred to as the Thermally Perfect Gas (TPG) code. The code is available free from the NASA Langley Software Server at URL http://www.larc.nasa.gov/LSS. The code produces tables of compressible flow properties similar to those found in NACA Report 1135. Unlike the NACA Report 1135 tables which are valid only in the calorically perfect temperature regime the TPG code results are also valid in the thermally perfect, calorically imperfect temperature regime, giving the TPG code a considerably larger range of temperature application. Accuracy of the TPG code in the calorically perfect and in the thermally perfect, calorically imperfect temperature regimes are verified by comparisons with the methods of NACA Report 1135. The advantages of the TPG code compared to the thermally perfect, calorically imperfect method of NACA Report 1135 are its applicability to any type of gas (monatomic, diatomic, triatomic, or polyatomic) or any specified mixture of gases, ease-of-use, and tabulated results.

Description: A fully automated iterative design method has been developed by which an airfoil with a substantial amount of natural laminar flow can be designed, while maintaining other aerodynamic and geometric constraints. Drag reductions have been realized using the design method over a range of Mach numbers, Reynolds numbers and airfoil thicknesses. The thrusts of the method are its ability to calculate a target N-Factor distribution that forces the flow to undergo transition at the desired location; the target-pressure-N-Factor relationship that is used to reduce the N-Factors in order to prolong transition; and its ability to design airfoils to meet lift, pitching moment, thickness and leading-edge radius constraints while also being able to meet the natural laminar flow constraint. The method uses several existing CFD codes and can design a new airfoil in only a few days using a Silicon Graphics IRIS workstation.

Description: A numerical large-eddy simulation model is under modification and testing for application to aircraft wake vortices. The model, having a meteorological framework, permits the interaction of wake vortices with environments characterized by crosswind shear, stratification, and humidity. As part of the validation process, model results are compared with measured field data from the 1990 Idaho Falls and the 1994-1995 Memphis field experiments. Cases are selected that represent different aircraft and a cross section of meteorological environments. Also included is one case with wake vortex generation in ground effect. The model simulations are initialized with the appropriate meteorological conditions and a post roll-up vortex system. No ambient turbulence is assumed in our initial set of experiments, although turbulence can be self generated by the interaction of the model wakes with the ground and environment.

Description: Current time-dependent CFD simulations produce very large multi-dimensional data sets at each time step. The visual analysis of computational results are traditionally performed by post processing the static data on graphics workstations. We present results from an alternate approach in which we analyze the simulation data in situ on each processing node at the time of simulation. The locally analyzed results, usually more economical and in a reduced form, are then combined and sent back for visualization on a graphics workstation.

Description: Topological concepts are used to derive separation conditions for two- and three-dimensional laminar flows. The result for two-dimensional flow reproduces the form of the well-known Stratford criterion. An extension makes the form applicable to the symmetry plane of a three-dimensional flow.

Description: An improved nozzle for reducing overspray in high temperature supersonic plasma spray devices comprises a body defining an internal passageway having an upstream end and a downstream end through which a selected plasma gas is directed. The nozzle passageway has a generally converging/diverging Laval shape with its upstream end converging to a throat section and its downstream end diverging from the throat section. The upstream end of the passageway is configured to accommodate a high current cathode for producing an electrical arc in the passageway to heat and ionize the gas flow to plasma form as it moves along the passageway. The downstream end of the nozzle is uniquely configured through the methodology of this invention to have a contoured bell-shape that diverges from the throat to the exit of the nozzle. Coating material in powder form is injected into the plasma flow in the region of the bell-shaped downstream end of the nozzle and the powder particles become entrained in the flow. The unique bell shape of the nozzle downstream end produces a plasma spray that is ideally expanded at the nozzle exit and thus virtually free of shock phenomena, and that is highly collimated so as to exhibit significantly reduced fanning and diffusion between the nozzle and the target. The overall result is a significant reduction in the amount of material escaping from the plasma stream in the form of overspray and a corresponding improvement in the cost of the coating operation and in the quality and integrity of the coating itself.

Description: This paper presents a robust method for the generation of zonal volume grids of design parametrics for aerodynamic configurations. The process utilizes simple algebraic techniques with parametric splines coupled with elliptic volume grid generation to generate isolated zonal grids for changes in body configuration needed to perform parametric design studies. Speed of the algorithm is maximized through the algebraic methods and reduced number of grid points to be regenerated for each design parametric without sacrificing grid quality and continuity within the volume domain. The method is directly applicable to grid reusability, because it modifies existing ow adapted volume grids and enables the user to restart the CFD solution process with an established flow field. Use of this zonal approach reduces computer usage time to create new volume grids for design parametric studies by an order of magnitude, as compared to current methods which require the regeneration of an entire volume grid. A sample configuration of a proposed Single Stage-to-Orbit Vehicle is used to illustrate an application of this method.

Description: Hypersonic boundary layer measurements were conducted over a flared cone in a quiet wind tunnel. The flared cone was tested at a freestream unit Reynolds number of 2.82x106/ft in a Mach 6 flow. This Reynolds number provided laminar-to-transitional flow over the model in a low-disturbance environment. Point measurements with a single hot wire using a novel constant voltage anemometry system were used to measure the boundary layer disturbances. Surface temperature and schlieren measurements were also conducted to characterize the laminar-to-transitional state of the boundary layer and to identify instability modes. Results suggest that the second mode disturbances were the most unstable and scaled with the boundary layer thickness. The integrated growth rates of the second mode compared well with linear stability theory in the linear stability regime. The second mode is responsible for transition onset despite the existence of a second mode sub-harmonic. The sub-harmonic wavelength also scales with the boundary layer thickness. Furthermore, the existence of higher harmonics of the fundamental suggests that non-linear disturbances are not associated with high free stream disturbance levels.

Description: Time-dependent CFD has been used to predict aerospace vehicle aerodynamics during a subsonic rotation maneuver. The inviscid 3D3U code is employed to solve the 3-D unsteady flow field using an unstructured grid of tetrahedra. As this application represents a challenge to time-dependent CFD, observations concerning spatial and temporal resolution are included. It is shown that even for a benign rotation rate, unsteady aerodynamic effects are significant during the maneuver. Possibly more significant, however, the rotation maneuver creates ow asymmetries leading to yawing moment, rolling moment, and side force which are not present in the quasi-steady case. A series of steady solutions at discrete points in the maneuver are also computed for comparison with wind tunnel measurements and as a means of quantifying unsteady effects.

Description: North American Marine Jet (NAMJ), Inc. received assistance from Marshall Space Flight Center engineers in the Computational Fluid Dynamics (CFD) branch of the Structure and Dynamics Laboratory in improving the proposed design of a new impeller for their jet-propulsion systems. Marshall used advanced CFD techniques, which included creating a three-dimensional computer model of the impeller for analysis. With Marshall input, the company modified the design, then Marshall used a computer model to make a solid polycarbonate model. The rapid prototyping allowed the company to avoid many time- consuming and costly steps in creating the impeller model. NAMJ is now able to compete with Pacific-area and European manufacturers who have traditionally dominated the market.

Description: On a swept wing, contamination along the leading edge, Tollmien-Schlichting waves, stationary or traveling crossflow vortices, and/or Taylor-Gortler vortices can cause the catastrophic breakdown of laminar to turbulent flow, which leads to increased skin-friction drag for the aircraft. The discussion in this Note will be limited to disturbances which evolve along the attachment line (leading edge of swept wing). If the Reynolds number of the attachment-line boundary layer is greater than some critical value, then the complete wing is inevitably engulfed in turbulent flow. Essentially, there are two critical Reynolds number points that must be considered. The first is for small-amplitude disturbances, and the second is for bypass transition. The present study will use direct numerical simulations to validate a linear 2D-eigenvalue prediction method based on parabolized stability equations by Lin and Malik. This method is considered because it suggests that a number of symmetric and asymmetric modes exist and are stable or unstable on the attachment line depending on the Reynolds number. If validated, the approach would predict a number of modes which are linearly damped in the Reynolds number regime 100 to 245; however, these modes may grow nonlinearly and provide an explanation to this region.

Description: The linear and the nonlinear stability of disturbances that propagate along the attachment line of a three-dimensional boundary layer is considered. The spatially evolving disturbances in the boundary layer are computed by direct numerical simulation (DNS) of the unsteady, incompressible Navier-Stokes equations. Disturbances are introduced either by forcing at the in ow or by applying suction and blowing at the wall. Quasi-parallel linear stability theory and a nonparallel theory yield notably different stability characteristics for disturbances near the critical Reynolds number; the DNS results con rm the latter theory. Previously, a weakly nonlinear theory and computations revealed a high wave-number region of subcritical disturbance growth. More recent computations have failed to achieve this subcritical growth. The present computational results indicate the presence of subcritically growing disturbances; the results support the weakly nonlinear theory. Furthermore, an explanation is provided for the previous theoretical and computational discrepancy. In addition, the present results demonstrate that steady suction can be used to stabilize disturbances that otherwise grow subcritically along the attachment line.

Description: The National PARC (NPARC) Alliance was established by the NASA Lewis Research Center and the Air Force Arnold Engineering Development Center to provide the U.S. aeropropulsion community with a reliable Navier-Stokes code for simulating the nonrotating components of propulsion systems. Recent improvements to the turbulence model capabilities of the NPARC code have significantly improved its capability to simulate turbulent flows. Specifically, the Chien k-epsilon and Wilcox k-omega turbulence models were implemented at Lewis. Lewis researchers installed the Chien k-epsilon model into NPARC to improve the code's ability to calculate turbulent flows with attached wall boundary layers and free shear layers. Calculations with NPARC have demonstrated that the Chien k-epsilon model provides more accurate calculations than those obtained with algebraic models previously available in the code. Grid sensitivity investigations have shown that computational grids must be packed against the solid walls such that the first point off of the wall is placed in the laminar sublayer. In addition, matching the boundary layer and momentum thicknesses entering mixing regions is necessary for an accurate prediction of the free shear-layer growth.

Description: In a continuing in-house program at the NASA Lewis Research Center, analytical and numerical methods are being developed to apply radiative analysis to predict transient temperature distributions and heat flows in partially transmitting materials. Results have been obtained for a single plane layer, and a transient analysis is being developed for a two-layer composite where each layer has a different refractive index. Because the ceramic refractive indices are larger than one, internal reflections are produced at the surfaces and at the internal interface. Reflections tend to distribute energy within a layer, and this affects the transient temperature distributions.

Description: A gas-jet diffusion flame is similar to the flame on a Bunsen burner, where a gaseous fuel (e.g., propane) flows from a nozzle into an oxygen-containing atmosphere (e.g., air). The difference is that a Bunsen burner allows for (partial) premixing of the fuel and the air, whereas a diffusion flame is not premixed and gets its oxygen (principally) by diffusion from the atmosphere around the flame. Simple gas-jet diffusion flames are often used for combustion studies because they embody the mechanisms operating in accidental fires and in practical combustion systems. However, most practical combustion is turbulent (i.e., with random flow vortices), which enhances the fuel/air mixing. These turbulent flames are not well understood because their random and transient nature complicates analysis. Normal gravity studies of turbulence in gas-jet diffusion flames can be impeded by buoyancy-induced instabilities. These gravitycaused instabilities, which are evident in the flickering of a candle flame in normal gravity, interfere with the study of turbulent gas-jet diffusion flames. By conducting experiments in microgravity, where buoyant instabilities are avoided, we at the NASA Lewis Research Center hope to improve our understanding of turbulent combustion. Ultimately, this could lead to improvements in combustor design, yielding higher efficiency and lower pollutant emissions. Gas-jet diffusion flames are often researched as model flames, because they embody mechanisms operating in both accidental fires and practical combustion systems (see the first figure). In normal gravity laboratory research, buoyant air flows, which are often negligible in practical situations, dominate the heat and mass transfer processes. Microgravity research studies, however, are not constrained by buoyant air flows, and new, unique information on the behavior of gas-jet diffusion flames has been obtained.

Description: The Interface Configuration Experiment (ICE) is actually a series of experiments that explore the striking behavior of liquid-vapor interfaces (i.e., fluid surfaces) in a low gravity environment under which major shifts in liquid position can arise from small changes in container shape or contact angle. Although these experiments are designed to test current mathematical theory, there are numerous practical applications that could result from these studies. When designing fluid management systems for space-based operations, it is important to be able to predict the locations and configurations that fluids will assume in containers under low-gravity conditions. The increased ability to predict, and hence control, fluid interfaces is vital to systems and/or processes where capillary forces play a significant role both in space and on the Earth. Some of these applications are in general coating processes (paints, pesticides, printing, etc.), fluid transport in porous media (ground water flows, oil recovery, etc.), liquid propellant systems in space (liquid fuel and oxygen), capillary-pumped loops and heat pipes, and space-based life-support systems. In space, almost every fluid system is affected, if not dominated, by capillarity. Knowledge of the liquid-vapor interface behavior, and in particular the interface shape from which any analysis must begin, is required as a foundation to predict how these fluids will react in microgravity and on Earth. With such knowledge, system designs can be optimized, thereby decreasing costs and complexity, while increasing performance and reliability. ICE has increased, and will continue to increase this knowledge, as it probes the specific peculiarities of current theory upon which our current understanding of these effects is based. Several versions of ICE were conducted in NASA Lewis Research Center's drop towers and on the space shuttle during the first and second United States Microgravity Laboratory missions (USML-1 and USML-2). Additional tests are planned for the space shuttle and for the Russian Mir space station. These studies will focus on interfacial problems concerning surface existence, uniqueness, configuration, stability, and flow characteristics.

Description: The development of thermodynamics and statistical mechanics is very important in the history of physics, and it underlines the difficulty in dealing with systems involving many bodies, even if those bodies are identical. Macroscopic systems of atoms typically contain so many particles that it would be virtually impossible to follow the behavior of all of the particles involved. Therefore, the behavior of a complete system can only be described or predicted in statistical ways. Under a grant to the NASA Lewis Research Center, scientists at the Case Western Reserve University have been examining the use of modern computing techniques that may be able to investigate and find the behavior of complete systems that have a large number of particles by tracking each particle individually. This is the study of molecular dynamics. In contrast to Monte Carlo techniques, which incorporate uncertainty from the outset, molecular dynamics calculations are fully deterministic. Although it is still impossible to track, even on high-speed computers, each particle in a system of a trillion trillion particles, it has been found that such systems can be well simulated by calculating the trajectories of a few thousand particles. Modern computers and efficient computing strategies have been used to calculate the behavior of a few physical systems and are now being employed to study important problems such as supersonic flows in the laboratory and in space. In particular, an animated video (available in mpeg format--4.4 MB) was produced by Dr. M.J. Woo, now a National Research Council fellow at Lewis, and the G-VIS laboratory at Lewis. This video shows the behavior of supersonic shocks produced by pistons in enclosed cylinders by following exactly the behavior of thousands of particles. The major assumptions made were that the particles involved were hard spheres and that all collisions with the walls and with other particles were fully elastic. The animated video was voted one of two winning videos in a competition held at the meeting of the American Physical Society's Division of Fluid Dynamics, held in Atlanta, Georgia, in November 1994. Of great interest was the result that in every shock there were a few high-speed precursor particles racing ahead of the shock, carrying information about its impending arrival. Most recently, Dr. Woo has been applying molecular dynamics techniques to the problem of determining the drag produced by the space station truss structure as it flies through the thin residual atmosphere of low-Earth orbit. This problem is made difficult by the complex structure of the truss and by the extreme supersonic nature of the flow. A fully filled section of the truss has already been examined, and drag predictions have been made. Molecular dynamics techniques promise to make realistic drag calculations possible even for very complex partially filled truss segments flying at arbitrary angles.

Description: The design and implementation of Doppler Global Velocimetry (DGV) for testing in the Langley Unitary Plan Wind Tunnel is presented. The discussion begins by outlining the characteristics of the tunnel and the test environment, with potential problem areas highlighted. Modifications to the optical system design to implement solutions for these problems are described. Since this tunnel entry was the first ever use of DGV in a supersonic wind tunnel, the test series was divided into three phases, each with its own goal. Phase I determined if condensation provided sufficient scattered light for DGV applications. Phase II studied particle lag by measuring the flow about an oblique shock above an inclined flat plate. Phase III investigated the supersonic vortical flow field above a 75-degree delta wing at 24-degrees angle of attack. Example results from these tests are presented.

Description: As opposed to previous explanations based on the effects of anharmonicity of simple diatomic molecules, traces of water vapor are suggested to be the most likely cause of the anomalously fast vibrational relaxation of such gases observed in supersonic and hypersonic nozzles. The mechanism is the strong V-VR coupling with H2O molecules that dramatically facilitates the collisional transfer of vibrational energy. Slight moisture content is thus a real world aspect of gas dynamics that must be considered in characterizations of shock tubes, reflected shock tunnels, and expansion tubes.