ALBERT

Description: Under SBIR (Small Business Innovative Research) contracts with Lewis Research Center, Nektonics, Inc., developed coating process simulation tools, known as Nekton. This powerful simulation software is used specifically for the modeling and analysis of a wide range of coating flows including thin film coating analysis, polymer processing, and glass melt flows. Polaroid, Xerox, 3M, Dow Corning, Mead Paper, BASF, Mitsubishi, Chugai, and Dupont Imaging Systems are only a few of the companies that presently use Nekton.

Description: A major portion of the energy produced in the world today comes from the burning of liquid hydrocarbon fuels in the form of droplets. Understanding the fundamental physical processes involved in droplet combustion is not only important in energy production but also in propulsion, in the mitigation of combustion-generated pollution, and in the control of the fire hazards associated with handling liquid combustibles. Microgravity makes spherically symmetric combustion possible, allowing investigators to easily validate their droplet models without the complicating effects of gravity. The Fiber-Supported Droplet Combustion (FSDC-2) investigation was conducted in the Microgravity Glovebox facility of the shuttles' Spacelab during the reflight of the Microgravity Science Laboratory (MSL- 1R) on STS-94 in July 1997. FSDC-2 studied fundamental phenomena related to liquid fuel droplet combustion in air. Pure fuels and mixtures of fuels were burned as isolated single and duo droplets with and without forced air convection. FSDC-2 is sponsored by the NASA Lewis Research Center, whose researchers are working in cooperation with several investigators from industry and academia. The rate at which a droplet burns is important in many commercial applications. The classical theory of droplet burning assumes that, for an isolated, spherically symmetric, single-fuel droplet, the gas-phase combustion processes are much faster than the droplet surface regression rate and that the liquid phase is at a uniform temperature equal to the boiling point. Recent, more advanced models predict that both the liquid and gas phases are unsteady during a substantial portion of the droplet's burning history, thus affecting the instantaneous and average burning rates, and that flame radiation is a dominant mechanism that can extinguish flames in a microgravity environment. FSDC-2 has provided well-defined, symmetric droplet burning data including radiative emissions to validate these theoretical models for heptane, decane, ethanol, and methanol fuels. Since most commercial combustion systems burn droplets in a convective environment, data were obtained without and with convective flow over the burning droplet (see the following photos).

Description: The paper summarizes Active Flow Control projects currently underway at the NASA Langley Research Center. Technology development is being pursued within a multidisciplinary, cooperative approach, involving the classical disciplines of fluid mechanics, structural mechanics, material science, acoustics, and stability and control theory. Complementing the companion papers in this session, the present paper will focus on projects that have the goal of extending the state-of-the-art in the measurement, prediction, and control of unsteady, nonlinear aerodynamics. Toward this goal, innovative actuators, micro and macro sensors, and control strategies are considered for high payoff flow control applications. The target payoffs are outlined within each section below. Validation of the approaches range from bench-top experiments to wind-tunnel experiments to flight tests. Obtaining correlations for future actuator and sensor designs are implicit in the discussion. The products of the demonstration projects and design tool development from the fundamental NASA R&D level technology will then be transferred to the Applied Research components within NASA, DOD, and US Industry. Keywords: active flow control, separation control, MEMS, review

Description: As part of a contract with the NASA Lewis Research Center, a simple, accurate method of predicting the performance characteristics of a nozzle design has been developed for use in conceptual design studies. The Nozzle Performance Analysis Code (NPAC) can predict the on- and off-design performance of axisymmetric or two-dimensional convergent and convergent-divergent nozzle geometries. NPAC accounts for the effects of overexpansion or underexpansion, flow divergence, wall friction, heat transfer, and small mass addition or loss across surfaces when the nozzle gross thrust and gross thrust coefficient are being computed. NPAC can be used to predict the performance of a given nozzle design or to develop a preliminary nozzle system design for subsequent analysis. The input required by NPAC consists of a simple geometry definition of the nozzle surfaces, the location of key nozzle stations (entrance, throat, exit), and the nozzle entrance flow properties. NPAC performs three analysis "passes" on the nozzle geometry. First, an isentropic control volume analysis is performed to determine the gross thrust and gross thrust coefficient of the nozzle. During the second analysis pass, the skin friction and heat transfer losses are computed. The third analysis pass couples the effects of wall shear and heat transfer with the initial internal nozzle flow solutions to produce a system of equations that is solved at steps along the nozzle geometry. Small mass additions or losses, such as those resulting from leakage or bleed flow, can be included in the model at specified geometric sections. A final correction is made to account for divergence losses that are incurred if the nozzle exit flow is not purely axial.

Description: Several schemes for introducing an artificial dissipation into a central difference approximation to the Euler and Navier Stokes equations are considered. The focus of the paper is on the convective upwind and split pressure (CUSP) scheme, which is designed to support single interior point discrete shock waves. This scheme is analyzed and compared in detail with scalar dissipation and matrix dissipation (MATD) schemes. Resolution capability is determined by solving subsonic, transonic, and hypersonic flow problems. A finite-volume discretization and a multistage time-stepping scheme with multigrid are used to compute solutions to the flow equations. Numerical solutions are also compared with either theoretical solutions or experimental data. For transonic airfoil flows the best accuracy on coarse meshes for aerodynamic coefficients is obtained with a simple MATD scheme. The coarse-grid accuracy for the original CUSP scheme is improved by modifying the limiter function used with the scheme, giving comparable accuracy to that obtained with the MATD scheme. The modifications reduce the background dissipation and provide control over the regions where the scheme can become first order.

Description: Aircraft laminar flow control (LFC) from the 1930's through the 1990's is reviewed and the current status of the technology is assessed. Examples are provided to demonstrate the benefits of LFC for subsonic and supersonic aircraft. Early studies related to the laminar boundary-layer flow physics, manufacturing tolerances for laminar flow, and insect-contamination avoidance are discussed. LFC concept studies in wind-tunnel and flight experiments are the major focus of the paper. LFC design tools are briefly outlined for completeness.

Description: Spacecraft solar dynamic power systems typically use high-temperature phase-change materials to efficiently store thermal energy for heat engine operation in orbital eclipse periods. Lithium fluoride salts are particularly well suited for this application because of their high heat of fusion, long-term stability, and appropriate melting point. Considerable attention has been focused on the development of thermal energy storage (TES) canisters that employ either pure lithium fluoride (LiF), with a melting point of 1121 K, or eutectic composition lithium-fluoride/calcium-difluoride (LiF-20CaF2), with a 1040 K melting point, as the phase-change material. Primary goals of TES canister development include maximizing the phase-change material melt fraction, minimizing the canister mass per unit of energy storage, and maximizing the phase-change material thermal charge/discharge rates within the limits posed by the container structure.

Description: Results from a recently completed enhanced mixing program are summarized in the two technical papers. These studies were parts of a High Speed Research (HSR)-supported joint Government/industry/university program that involved, in addition to the NASA Lewis Research Center, researchers at United Technologies Research Center, Allison Engine Company, CFD Research Corporation, and the University of California, Irvine. The studies investigated the mixing of jets injected normal to a confined subsonic mainsteam in both rectangular and cylindrical ducts. Experimental and computational studies were performed in both nonreacting and reacting flows. The orifice geometries and flow conditions were selected as typical of the complex three-dimensional flows in the combustion chambers in low-emission gas turbine engines. The principal conclusion from both the experiments and modeling was that the momentum-flux ratio J and orifice spacing S/H were the most significant flow and geometry variables, respectively. Conserved scalar distributions were similar-independent of reaction, orifice diameter H/d, and shape-when the orifice spacing and the square root of the momentum-flux ratio were inversely proportional. Jet penetration was critical, and penetration decreased as either momentum-flux ratio or orifice spacing decreased. We found that planar averages must be considered in context with the distributions. The mass-flow ratios and the orifices investigated were often very large. The jet-to-mainstream mass-flow ratio was varied from significantly less than 1 to greater than 1. The orifice-area to mainstream-cross-sectional-area was varied from approx. 0 to 0.5, and the axial planes of interest were often just downstream of the orifice trailing edge. Three-dimensional flow was a key part of efficient mixing and was observed for all configurations. As an example of the results, the accompanying figure shows the effects of different rates of mass addition on the opposite walls of a rectangular duct.

Description: Capillary-pumped loops (CPL's) are devices that are used to transport heat from one location to another--specifically to transfer heat away from something. In low-gravity applications, such as satellites (and possibly the International Space Station), CPL's are used to transfer heat from electrical devices to space radiators. This is accomplished by evaporating one liquid surface on the hot side of the CPL and condensing the vapor produced onto another liquid surface on the cold side. Capillary action, the phenomenon that causes paper towels to absorb spilled liquids, is used to "pump" the liquid back to the evaporating liquid surface (hot side) to complete the "loop." CPL's require no power to operate and can transfer heat over distances as large as 30 ft or more. Their reliance upon evaporation and condensation to transfer heat makes them much more economical in terms of weight than conventional heat transfer systems. Unfortunately, they have proven to be unreliable in space operations, and the explanation for this unreliability has been elusive. The Capillary-Driven Heat Transfer (CHT) experiment is investigating the fundamental fluid physics phenomena thought to be responsible for the failure of CPL's in low-gravity operations. If the failure mechanism can be identified, then appropriate design modifications can be developed to make capillary phase-change heat-transport devices a more viable option in space applications. CHT was conducted onboard the Space Shuttle Columbia during the first Microgravity Science Laboratory (MSL-1) mission, STS-94, which flew from July 1 to 17, 1997. The CHT glovebox investigation, which was conceived by Dr. Kevin Hallinan and Jeffrey Allen of the University of Dayton, focused on studying the dynamics associated with the heating and cooling at the evaporating meniscus within a capillary phase-change device in a low-gravity environment. The CHT experimental hardware was designed by a small team of engineers from Aerospace Design & Fabrication (ADF), the NASA Lewis Research Center, and the University of Dayton. The hardware consisted of two experiment modules that each contained an instrumented test loop (idealized capillary-pumped loop), a base unit for power conversion and backlighting, a display unit with 15 LED's (light-emitting diodes) to display temperatures, pressure, heater power, and time, a control unit to select heaters and heater settings, a cooling fan, and associated cables.

Description: Applying smart materials in aeropropulsion systems may improve the performance of aircraft engines through a variety of vibration, noise, and shape-control applications. To facilitate the experimental characterization of these smart structures, researchers have been focusing on developing analytical models to account for the coupled mechanical, electrical, and thermal response of these materials. One focus of current research efforts has been directed toward incorporating a comprehensive thermal analysis modeling capability. Typically, temperature affects the behavior of smart materials by three distinct mechanisms: Induction of thermal strains because of coefficient of thermal expansion mismatch 1. Pyroelectric effects on the piezoelectric elements; 2. Temperature-dependent changes in material properties; and 3. Previous analytical models only investigated the first two thermal effects mechanisms. However, since the material properties of piezoelectric materials generally vary greatly with temperature (see the graph), incorporating temperature-dependent material properties will significantly affect the structural deflections, sensory voltages, and stresses. Thus, the current analytical model captures thermal effects arising from all three mechanisms through thermopiezoelectric constitutive equations. These constitutive equations were incorporated into a layerwise laminate theory with the inherent capability to model both the active and sensory response of smart structures in thermal environments. Corresponding finite element equations were formulated and implemented for both the beam and plate elements to provide a comprehensive thermal effects modeling capability.

Description: In our world-class, optically accessible combustion facility at the NASA Lewis Research Center, we have developed the unique capability of making three-dimensional fuel distribution measurements of aviation gas turbine fuel injectors at actual operating conditions. These measurements are made in situ at the actual operating temperatures and pressures using the JP-grade fuels of candidate next-generation advanced aircraft engines for the High Speed Research (HSR) and Advanced Subsonics Technology (AST) programs. The inlet temperature and pressure ranges used thus far are 300 to 1100 F and 80 to 250 psia. With these data, we can obtain the injector spray angles, the fuel mass distributions of liquid and vapor, the degree of fuel vaporization, and the degree to which fuel has been consumed. The data have been used to diagnose the performance of injectors designed both in-house and by major U.S. engine manufacturers and to design new fuel injectors with overall engine performance goals of increased efficiency and reduced environmental impact. Mie scattering is used to visualize the liquid fuel, and laser-induced fluorescence is used to visualize both liquid and fuel vapor.

Description: Sprays occur in a wide variety of industrial and power applications and in the processing of materials. A liquid spray is a phase flow with a gas as the continuous phase and a liquid as the dispersed phase (in the form of droplets or ligaments). Interactions between the two phases, which are coupled through exchanges of mass, momentum, and energy, can occur in different ways at different times and locations involving various thermal, mass, and fluid dynamic factors. An understanding of the flow, combustion, and thermal properties of a rapidly vaporizing spray requires careful modeling of the rate-controlling processes associated with the spray's turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates, as well as other phenomena. In an attempt to advance the state-of-the-art in multidimensional numerical methods, we at the NASA Lewis Research Center extended our previous work on sprays to unstructured grids and parallel computing. LSPRAY, which was developed by M.S. Raju of Nyma, Inc., is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo probability density function (PDF) solver. The LSPRAY solver accommodates the use of an unstructured mesh with mixed triangular, quadrilateral, and/or tetrahedral elements in the gas-phase solvers. It is used specifically for fuel sprays within gas turbine combustors, but it has many other uses. The spray model used in LSPRAY provided favorable results when applied to stratified-charge rotary combustion (Wankel) engines and several other confined and unconfined spray flames. The source code will be available with the National Combustion Code (NCC) as a complete package.

Description: The Interface Configuration Experiment (ICE) is part of a multifaceted study that is exploring the often striking behavior of liquid-vapor interfaces in low-gravity environments. Although the experiment was posed largely as a test of current mathematical theory, applications of the results should be manifold. In space almost every fluid system is affected, if not dominated, by capillarity (the effects of surface tension). As a result, knowledge of fluid interface behavior, in particular an equilibrium interface shape from which any analysis must begin, is fundamental--from the control of liquid fuels and oxygen in storage tanks to the design and development of inspace thermal systems, such as heat pipes and capillary pumped loops. ICE has increased, and should continue to increase, such knowledge as it probes the specific peculiarities of current theory upon which our present understanding rests. Several versions of ICE have been conducted in the drop towers at the NASA Lewis Research Center, on the space shuttles during the first and second United States Microgravity Laboratory missions (USML-1 and USML-2), and most recently aboard the Russian Mir space station. These studies focused on interfacial problems concerning the existence, uniqueness, configuration, stability, and flow characteristics of liquid-vapor interfaces. Results to date have clearly demonstrated the value of the present theory and the extent to which it can predict the behavior of capillary systems.

Description: Particle Imaging Velocimetry (PIV) is an optical technique whereby a pulsed laser sheet is used to illuminate particles entrained in a fluid across an extended planar cross section of a flow field. Electronic recording of the particle positions at two closely timed laser pulses permits the computation of the flow velocity. PIV captures the instantaneous flow field, permitting the study of unsteady flow phenomena. Mean flow statistics can be computed by acquiring several hundred images and averaging the results. The first-ever successful application of PIV to acquire measurements in a high-speed rotating turbomachinery blade row was completed in NASA Lewis Research Center's W-8 Single Stage Axial Compressor Facility. Measurements were acquired in a 20-in.-diameter transonic compressor rotor operating at 17,188 rpm. A custom-designed light-sheet generating probe was used to insert the high-energy, pulsed light-sheet illumination required for recording the unblurred images of particles entrained in the fluid. Measurements of the shock wave formed within the rotor blade passage and of unsteady structures within the blade wakes were acquired. These measurements provide insight into unsteady spatial structures in the flow field which cannot be measured with the more conventional laser anemometry technique. The PIV technique provides both instantaneous and average velocity data in a transonic compressor in an order of magnitude less time than required for other conventional optical diagnostic techniques.

Description: Ceramic thermal barrier coatings have received increasing attention for advanced gas turbine and diesel engine applications because of their ability to provide thermal insulation to engine components. However, the durability of these coatings under the severe thermal cycling conditions encountered in a diesel engine (ref. 1) still remains a major issue. In this research at the NASA Lewis Research Center, a high-power laser was used to investigate the thermal fatigue behavior of a yttria-stabilized zirconia coating system under simulated diesel engine conditions. The mechanisms of fatigue crack initiation and propagation, and of coating failure under complex thermal low-cycle fatigue (LCF, representing stop/start cycles) and thermal high-cycle fatigue (HCF, representing operation at 1300 rpm) are described. Continuous wave and pulse laser modes were used to simulate pure LCF and combined LCF/HCF, respectively (ref. 2). The LCF mechanism was found to be closely related to the coating sintering and creep at high temperatures. These creep strains in the ceramic coating led to a tensile stress state during cooling, thus providing the major driving force for crack growth under LCF conditions. The combined LCF/HCF tests induced more severe coating surface cracking, microspallation, and accelerated crack growth than did the pure LCF test. HCF thermal loads also facilitated lateral crack branching and ceramic/bond coat interface delaminations. HCF is associated with the cyclic stresses originating from the high-frequency temperature fluctuation at the ceramic coating surface. The HCF thermal loads act on the crack by a wedging mechanism (ref. 1), resulting in continuous crack growth at temperature. The HCF stress intensity factor amplitude increases with the interaction depth and temperature swing, and decreases with the crack depth. HCF damage also increases with the thermal expansion coefficient and the Young's modulus of the ceramic coating (refs. 1 and 3).

Description: One of the most challenging areas of research in aerodynamics is the reduction of skin friction, especially for turbulent flow. Reduced skin friction means less drag. For aircraft, less drag can lead to less fuel burned or to a greater flight range for a fixed amount of fuel. Many techniques and methods have been tried; however, none of them has significantly reduced skin friction in the flight environment. An innovative skin-friction reduction technique, the Microblowing Technique (MBT), was invented in 1993. This is a unique concept in which an extremely small amount of air is blown vertically at a surface through very small holes. It can be used for aircraft or marine vehicles, such as submarines (where water is blown through the holes instead of air). As shown in the figure, the outer layer, which controls vertical flow, is a plate with high-vertical holes. The inner layer, which produces evenly distributed flow, is a low-permeability porous plate. Microblowing reduces the surface roughness and changes the flow velocity profile on the surface, thereby reducing skin friction.

Description: The "CFD General Notation System" (CGNS) consists of a collection of conventions, and conforming software, for the storage and retrieval of Computational Fluid Dynamics (CFD) data. It facilitates the exchange of data between sites and applications, and helps stabilize the archiving of aerodynamic data. This effort was initiated in order to streamline the procedures in exchanging data and software between NASA and its customers, but the goal is to develop CGNS into a National Standard for the exchange of aerodynamic data. The CGNS development team is comprised of members from Boeing Commercial Airplane Group, NASA-Ames, NASA-Langley, NASA-Lewis, McDonnell-Douglas Corporation (now Boeing-St. Louis), Air Force-Wright Lab., and ICEM-CFD Engineering. The elements of CGNS address all activities associated with the storage of data on external media and its movement to and from application programs. These elements include: - The Advanced Data Format (ADF) Database manager, consisting of both a file format specification and its I/O software, which handles the actual reading and writing of data from and to external storage media; - The Standard Interface Data Structures (SIDS), which specify the intellectual content of CFD data and the conventions governing naming and terminology; - The SIDS-to-ADF File Mapping conventions, which specify the exact location where the CFD data defined by the SIDS is to be stored within the ADF file(s); and - The CGNS Mid-level Library, which provides CFD-knowledgeable routines suitable for direct installation into application codes. The CGNS Mid-level Library was designed to ease the implementation of CGNS by providing developers with a collection of handy I/O functions. Since knowledge of the ADF core is not required to use this library, it will greatly facilitate the task of interfacing with CGNS. There are currently 48 user callable functions that comprise the Mid-level library and are described in the Users Guide. The library is written in C, but each function has a FORTRAN counterpart.

Description: The flow fields of vortices, whether bouyancy-driven or lift-generated, are fascinating fluid-dynamic phenomena which often possess intense swirl velocities and complex time-dependent behavior. As part of the on-going study of vortex behavior, this paper presents a historical overview of the research conducted on the structure and modification of the vortices generated by the lifting surfaces of subsonic transport aircraft. It is pointed out that the characteristics of lift-generated vortices are related to the aerodynamic shapes that produce them and that various arrangements of surfaces can be used to produce different vortex structures. The primary purpose of the research to be described is to find a way to reduce the hazard potential of lift-generated vortices shed by subsonic transport aircraft in the vicinity of airports during landing and takeoff operations. It is stressed that lift-generated vortex wakes are so complex that progress towards a solution requires application of a combined theoretical and experimental research program because either alone often leads to incorrect conclusions. It is concluded that a satisfactory aerodynamic solution to the wake-vortex problem at airports has not yet been found but a reduction in the impact of the wake-vortex hazard on airport capacity may become available in the foreseeable future through wake-vortex avoidance concepts currently under study. The material to be presented in this overview is drawn from aerospace journals that are available publicly.

Description: The "CFD General Notation System" (CGNS) consists of a collection of conventions, and conforming software, for the storage and retrieval of Computational Fluid Dynamics (CFD) data. It facilitates the exchange of data between sites and applications, and helps stabilize the archiving of aerodynamic data. This effort was initiated in order to streamline the procedures in exchanging data and software between NASA and its customers, but the goal is to develop CGNS into a National Standard for the exchange of aerodynamic data. The CGNS development team is comprised of members from Boeing Commercial. Airplane Group, NASA-Ames, NASA-Langley, NASA-Lewis, McDonnell-Douglas Corporation (now Boeing-St. Louis), Air Force-Wright Lab., and ICEM-CFD Engineering. The elements of CGNS address all activities associated with the storage of data on external media and its movement to and from application programs. These elements include: 1) The Advanced Data Format (ADF) Database manager, consisting of both a file format specification and its 1/0 software, which handles the actual reading and writing of data from and to external storage media; 2) The Standard Interface Data Structures (SIDS), which specify the intellectual content of CFD data and the conventions governing naming and terminology; 3) The SIDS-to-ADF File Mapping conventions, which specify the exact location where the CFD data defined by the SIDS is to be stored within the ADF file(s); and 4) The CGNS Mid-level Library, which provides CFD-knowledgeable routines suitable for direct installation into application codes. The ADF is a generic database manager with minimal intrinsic capability. It was written for the purpose of storing large numerical datasets in an efficient, platform independent manner. To be effective, it must be used in conjunction with external agreements on how the data will be organized within the ADF database such defined by the SIDS. There are currently 34 user callable functions that comprise the ADF Core library and are described in the Users Guide. The library is written in C, but each function has a FORTRAN counterpart.

Description: Experiments and computations were carried out on the adiabatic laminar boundary layer developing along the surfaces of a two-dimensional supersonic nozzle, consisting of upper and lower contoured nozzle blocks and flat sidewalls. Two- and three-dimensional Navier-Stokes codes, as well as two-dimensional boundary-layer codes were employed. These codes were 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 were made for the boundary-layer growth on the nozzle contoured surfaces, 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 laminar boundary-layer flow on both the sidewalls and the contoured surfaces. Theory and experiment both indicated that the sidewall flow is highly three-dimensional, with non-uniform shear, corner vortices and a boundary layer strongly distorted by cross flows induced by lateral pressure gradients.

Description: This summary presents the motivation for the Special Section on the credibility of computational fluid dynamics (CFD) simulations, its objective, its background and context, its content, and its major conclusions. Verification and validation (V&V) are the processes for establishing the credibility of CFD simulations. Validation assesses whether correct things are performed and verification assesses whether they are performed correctly. Various aspects of V&V are discussed. Progress is made in verification of simulation models. Considerable effort is still needed for developing a systematic validation method that can assess the credibility of simulated reality.

Description: Workload characterization is used for modeling and evaluating of computing systems at different levels of detail. We present workload characterization for a class of Computational Fluid Dynamics (CFD) applications that solve Partial Differential Equations (PDEs). This workload characterization focuses on three high performance computing platforms: SGI Origin2000, EBM SP-2, a cluster of Intel Pentium Pro bases PCs. We execute extensive measurement-based experiments on these platforms to gather statistics of system resource usage, which results in workload characterization. Our workload characterization approach yields a coarse-grain resource utilization behavior that is being applied for performance modeling and evaluation of distributed high performance metacomputing systems. In addition, this study enhances our understanding of interactions between PDE solver workloads and high performance computing platforms and is useful for tuning these applications.

Description: A review of the algorithmic features and capabilities of the unstructured-grid flow solver USM3Dns is presented. This code, along with the tetrahedral grid generator, VGRIDns, is being extensively used throughout the U.S. for solving the Euler and Navier-Stokes equations on complex aerodynamic problems. Spatial discretization is accomplished by a tetrahedral cell-centered finite-volume formulation using Roe's upwind flux difference splitting. The fluxes are limited by either a Superbee or MinMod limiter. Solution reconstruction within the tetrahedral cells is accomplished with a simple, but novel, multidimensional analytical formula. Time is advanced by an implicit backward-Euler time-stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one-equation model, which is coupled with a wall function to reduce the number of cells in the near-wall region of the boundary layer. The issues of accuracy and robustness of USM3Dns Navier-Stokes capabilities are addressed for a flat-plate boundary layer, and a full F-16 aircraft with external stores at transonic speed.

Description: The use of conservation laws in nonconservative form for deriving shock jump conditions by Schwartz distribution theory leads to ambiguous products of generalized functions. Nonstandard analysis is used to define a class of Heaviside functions where the jump from zero to one occurs on an infinitesimal interval. These Heaviside functions differ by their microstructure near x = 0, i.e., by the nature of the rise within the infinitesimal interval it is shown that the conservation laws in nonconservative form can relate the different Heaviside functions used to define jumps in different flow parameters. There are no mathematical or logical ambiguities in the derivation of the jump conditions. An important result is that the microstructure of the Heaviside function of the jump in entropy has a positive peak greater than one within the infinitesimal interval where the jump occurs. This phenomena is known from more sophisticated studies of the structure of shock waves using viscous fluid assumption. However, the present analysis is simpler and more direct.

Description: An experimental and computational study was conducted on an exhaust nozzle with fluidic injection for yaw thrust-vector control. The nozzle concept was tested experimentally in the NASA Langley Jet Exit Test Facility (JETF) at nozzle pressure ratios up to 4 and secondary fluidic injection flow rates up to 15 percent of the primary flow rate. Although many injection-port geometries and two nozzle planforms (symmetric and asymmetric) were tested experimentally, this paper focuses on the computational results of the more successful asymmetric planform with a slot injection port. This nozzle concept was simulated with the Navier-Stokes flow solver, PAB3D, invoking the Shih, Zhu, and Lumley algebraic Reynolds stress turbulence model (ASM) at nozzle pressure ratios (NPRs) of 2,3, and 4 with secondary to primary injection flow rates (w(sub s)/w(sub p)) of 0, 2, 7 and 10 percent.

Description: The use of a high molecular weight test gas to increase the Reynolds number range of transonic wind tunnels is explored. Modifications to a small transonic wind tunnel are described and the real gas properties of the example heavy gas (sulfur hexafluoride) are discussed. Sulfur hexafluoride is shown to increase the test Reynolds number by a factor of more than 2 over air at the same Mach number. Experimental and computational pressure distributions on an advanced supercritical airfoil configuration at Mach 0.7 in both sulfur hexafluoride and nitrogen are presented. Transonic similarity theory is shown to be partially successful in transforming the heavy gas results to equivalent nitrogen (air) results, provided the correct definition of gamma is used.

Description: This paper presents grid generation techniques available in the Volume Grid Manipulation (VGM) code. The VGM code is designed to manipulate existing line, surface and volume grids to improve the quality of the data. It embodies an easy to read rich language of commands that enables such alterations as topology changes, grid adaption and smoothing. Additionally, the VGM code can be used to construct simplified straight lines, splines, and conic sections which are common curves used in the generation and manipulation of points, lines, surfaces and volumes (i.e., grid data). These simple geometric curves are essential in the construction of domain discretizations for computational fluid dynamic simulations. By comparison to previously established methods of generating these curves interactively, the VGM code provides control of slope continuity and grid point-to-point stretchings as well as quick changes in the controlling parameters. The VGM code offers the capability to couple the generation of these geometries with an extensive manipulation methodology in a scripting language. The scripting language allows parametric studies of a vehicle geometry to be efficiently performed to evaluate favorable trends in the design process. As examples of the powerful capabilities of the VGM code, a wake flow field domain will be appended to an existing X33 Venturestar volume grid; negative volumes resulting from grid expansions to enable flow field capture on a simple geometry, will be corrected; and geometrical changes to a vehicle component of the X33 Venturestar will be shown.

Description: The heavy gas test medium has recently been changed in the Transonic Dynamics Tunnel (TDT) at the NASA Langley Research Center. A NASA Construction of Facilities project has converted the TDT heavy gas from dichlorodifluoromethane (R12) to 1,1,1,2 tetrafluoroethane (R134a). The facility s heavy gas processing system was extensively modified to implement the conversion to R134a. Additional system modifications have improved operator interfaces, hardware reliability, and quality of the research data. The facility modifications included improvements to the heavy gas compressor and piping, the cryogenic heavy gas reclamation system, and the heavy gas control room. A series of wind tunnel characterization and calibration tests are underway. Results of the flow characterization tests show the TDT operating envelope in R134a to be very similar to the previous operating envelope in R12.

Description: Detailed aeroheating information is critical to the successful design of a thermal protection system (TPS) for an aerospace vehicle. This report describes NASA Langley Research Center's (LaRC) two-color relative-intensity phosphor thermography method and the IHEAT software package which is used for the efficient data reduction and analysis of the phosphor image data. Development of theory is provided for a new weighted two-color relative-intensity fluorescence theory for quantitatively determining surface temperatures on hypersonic wind tunnel models; an improved application of the one-dimensional conduction theory for use in determining global heating mappings; and extrapolation of wind tunnel data to flight surface temperatures. The phosphor methodology at LaRC is presented including descriptions of phosphor model fabrication, test facilities and phosphor video acquisition systems. A discussion of the calibration procedures, data reduction and data analysis is given. Estimates of the total uncertainties (with a 95% confidence level) associated with the phosphor technique are shown to be approximately 8 to 10 percent in the Langley's 31-Inch Mach 10 Tunnel and 7 to 10 percent in the 20-Inch Mach 6 Tunnel. A comparison with thin-film measurements using two-inch radius hemispheres shows the phosphor data to be within 7 percent of thin-film measurements and to agree even better with predictions via a LATCH computational fluid dynamics solution (CFD). Good agreement between phosphor data and LAURA CFD computations on the forebody of a vertical takeoff/vertical lander configuration at four angles of attack is also shown. In addition, a comparison is given between Mach 6 phosphor data and laminar and turbulent solutions generated using the LAURA, GASP and LATCH CFD codes. Finally, the extrapolation method developed in this report is applied to the X-34 configuration with good agreement between the phosphor extrapolation and LAURA flight surface temperature predictions. The phosphor process outlined in the paper is believed to provide the aerothermodynamic community with a valuable capability for rapidly obtaining (4 to 5 weeks) detailed heating information needed in TPS design.

Description: Under an SBIR (Small Business Innovative Research) contract with Goddard Space Flight Center, (NAS5-32730), Frederick A. Costello, Inc., has developed a low cost computer program for integrating thermal modeling with other engineering disciplines. TCON(tm) Concurrent Engineering Software enables graphical user-interface programs to create three dimensional models. This technology reduces time, and facilitates accuracy while performing these applications. TCON(tm) has been sold to NASA and major US industries.

Description: This report presents work performed by The Boeing Company to satisfy the deliverable "Evaluation of aeroservoelastic Effects on Symmetric Flutter" for Subtask 7 of Reference 1. The objective of this report is to incorporate the improved methods for studying the effects of a closed-loop control system on the aeroservoelastic behavior of the airplane planned under NASA HSR technical Integration Task 20 work. Also, a preliminary evaluation of the existing pitch control laws on symmetric flutter of the TCA configuration was addressed."The goal is to develop an improved modeling methodology and perform design studies that account for the aero-structures-systems interaction effects.

Description: The objective of this study is to investigate the effect of nonlinear aerodynamics, especially at high angles-of-attack with leading-edge separation, on the TCA flutter properties at transonic speeds. In order to achieve that objective, flutter simulations with Navier-Stokes CFD must be performed. To this end, time-marching Navier-Stokes solutions are computed for the TCA wing/body configuration at high angles-of-attack in transonic flight regimes. The approach is to perform non-linear flutter calculations on the TCA at two angles-of-attack, the first one being a case with attached flow (a=2.8 degrees) and the second one being a high angle-of-attack case with a wing leading edge vortex (a=12.11 degrees). Comparisons of the resulting histories and frequency damping information for both angles-of-attack will evaluate the impact of high-alpha aerodynamics on flutter.

Description: We present results of inviscid simulations, in three dimensions, of Richtmyer-Meshkov instability for high incident shock Mach number. The growth rate of a single harmonic perturbation is quantified and compared with the results of a 2D calculation. Upon re-shock, the perturbation amplitude undergoes a phase reversal while the mean velocity of the interface is zero. Before re-shock the normalized growth rate of a 2D and 3D interface are nearly the same, but the growth rate after re-shock is significantly larger for the 3D than the 2D case. We also examine the evolution of multiple harmonic perturbations. Computational and parallelization issues of the simulation code will also be briefly discussed. The computations were done on the T3E at Pittsburgh Supercomputing Center.

Description: A series of flight tests over the ocean of a four engine turboprop airplane in the cruise configuration have provided a data set for improved understanding of wake vortex physics and atmospheric interaction. An integrated database has been compiled for wake characterization and validation of wake-vortex computational models. This paper describes the wake-vortex flight tests, the data processing, the database development and access, and results obtained from preliminary wake-characterization analysis using the data sets.

Description: A non-invasive flow analysis system and method wherein a sensor, such as an acoustic sensor, is coupled to a conduit for transmitting a signal which varies depending on the characteristics of the flow in the conduit. The signal is amplified and there is a filter, responsive to the sensor signal, and tuned to pass a narrow band of frequencies proximate the resonant frequency of the sensor. A demodulator generates an amplitude envelope of the filtered signal and a number of flow indicator quantities are calculated based on variations in amplitude of the amplitude envelope. A neural network, or its equivalent, is then used to determine the flow rate of the flow in the conduit based on the flow indicator quantities.

Description: A flying wire system was utilized in conjunction with a rake of nine cross-wire probes to obtain simultaneous velocity measurements in an axisymmetric sudden expansion at a Reynolds number of 41,000. From these measurements, the correlation tensor could be calculated. Knowledge of the two-point correlation tensor reveals more in-depth information of the physical attributes of this flow. The two point correlation tensor allowed for calculation of the integrated length scales in both the radial and axial directions. This gives insight into the growth of structures with increasing downstream distance and at different radial locations through out the sudden expansion. The length scales were calculated by integrating the two-point correlation tensor in the radial direction from the centerline to the outer pipe wall and by integrating between several step heights for the axial direction. Calculated correlations at z/h = 8 and 9 at r/R = 0.46 showed a correlation length of 1/3 step height for the radial direction. It was found that length scales in the radial direction became larger with increasing radius with peaks at 0.70 e r/R c 0.81 but then decreased slightly towards the wall. Length scales in the axial direction yielded a recirculating bubble on the order of 3 step heights in the recirculating region. After the recirculating region, the length scales decreased to 1/4 of a step height.

Description: A jet in cross-flow (JIFC) consists of a jet exhausting at a large angle into a freestream flow. It is a flow field which is relevant to a wide variety of technologies and applications. Despite the nearly 65 years of JIFC research there are few results available for laminar hypersonic flows, a combination which will be encountered by re-entry and high altitude vehicles over some portion of their flight path. This research consists of developing a numerical model to investigate the interaction of a normal sonic jet exhausting into a hypersonic cross-flow. The model was validated by comparing experimental measurements with corresponding numerical results generated by the model.

Description: Chapter 1 briefly reviews several related topics associated with the symmetrization of systems of conservation laws and quasi-conservation laws: (1) Basic Entropy Symmetrization Theory; (2) Symmetrization and eigenvector scaling; (3) Symmetrization of the compressible Navier-Stokes equations; and (4) Symmetrization of the quasi-conservative form of the magnetohydrodynamic (MHD) equations. Chapter 2 describes one of the best known tools employed in the study of differential equations, the maximum principle: any function f(x) which satisfies the inequality f(double prime)〉0 on the interval [a,b] attains its maximum value at one of the endpoints on the interval. Chapter three examines the upwind finite volume schemes for scalar and system conservation laws. The basic tasks in the upwind finite volume approach have already been presented: reconstruction, flux evaluation, and evolution. By far, the most difficult task in this process is the reconstruction step.

Description: A capillary pumped loop for transferring heat from one body part to another body part, the capillary pumped loop comprising a capillary evaporator for vaporizing a liquid refrigerant by absorbing heat from a warm body part, a condenser for turning a vaporized refrigerant into a liquid by transferring heat from the vaporized liquid to a cool body part, a first tube section connecting an output port of the capillary evaporator to an input of the condenser, and a second tube section connecting an output of the condenser to an input port of the capillary evaporator. A wick may be provided within the condenser. A pump may be provided between the second tube section and the input port of the capillary evaporator. Additionally, an esternal heat source or heat sink may be utilized.

Description: Description, capabilities, recent upgrades, and utilization of the NASA Langley Research Center (LaRC) Aerothermodynamic Facilities Complex (AFC) are presented. The AFC consists of five hypersonic, blow-down-to-vacuum wind tunnels that collectively provide a range of Mach number from 6 to 20, unit Reynolds number from 0.04 to 22 million per foot and, most importantly for blunt configurations, normal shock density ratio from 4 to 12. These wide ranges of hypersonic simulation parameters are due, in part, to the use of three different test gases (air, helium, and tetrafluoromethane), thereby making several of the facilities unique. The Complex represents nearly three-fourths of the conventional (as opposed to impulse)-type hypersonic wind tunnels operational in this country. AFC facilities are used to assess and optimize the hypersonic aerodynamic performance and aeroheating characteristics of aerospace vehicle concepts and to provide benchmark aerodynamic/aeroheating data fr generating the flight aerodynamic databook and final design of the thermal protection system (TPS) (e.g., establishment of flight limitations not to exceed TPS design limits). Modifications and enhancements of AFC hardware components and instrumentation have been pursued to increase capability, reliability, and productivity in support of programmatic goals. Examples illustrating facility utilization in recent years to generate essentially all of the experimental hypersonic aerodynamic and aeroheating information for high-priority, fast-paced Agency programs are presented. These programs include Phase I of the Reusable Launch Vehicle (RLV) Advanced Technology Demonstrator, X-33 program, PHase II of the X-33 program, X-34 program, the Hyper-X program ( a Mach 5,7, and 10 airbreathing propulsion flight experiment), and the X-38 program (Experimental Crew Return Vehicle, X-CRV). Current upgrades/enchancements and future plans for the AFC are discussed.

Description: The X-33 experimental program is a cooperative program between industry and NASA, managed by Lockheed-Martin Skunk Works to develop an experimental vehicle to demonstrate new technologies for a single-stage-to-orbit, fully reusable launch vehicle (RLV). One of the new technologies to be demonstrated is an advanced Thermal Protection System (TPS) being designed by BF Goodrich (formerly Rohr, Inc.) with support from NASA. The calculation of an aerothermal database is crucial to identifying the critical design environment data for the TPS. The NASA Ames X-33 team has generated such a database using Computational Fluid Dynamics (CFD) analyses, engineering analysis methods and various programs to compare and interpolate the results from the CFD and the engineering analyses. This database, along with a program used to query the database, is used extensively by several X-33 team members to help them in designing the X-33. This paper will describe the methods used to generate this database, the program used to query the database, and will show some of the aerothermal analysis results for the X-33 aircraft.

Description: The shear-sensitive liquid crystal coating (SSLCC) method is an image-based technique for both visualizing dynamic surface-flow phenomena, such as transition and separation, and for measuring the continuous shear-stress vector distribution acting on an aerodynamic surface. Under proper lighting and viewing conditions (discussed below), the coating changes color in response to an applied aerodynamic shear. This color-change response is continuous and reversible, with a response time of milliseconds, and is a function of both the shear magnitude and the shear vector orientation relative to the observer. The liquid crystal phase of matter is a weakly-ordered, viscous, non-Newtonian fluid state that exists between the nonuniform liquid phase and the ordered solid phase of certain organic compounds. Cholesteric liquid crystal compounds possess a helical molecular arrangement that selectively scatters white light, incident along the helical axis, as a three-dimensional spectrum. This property is linked to the helical pitch length, which is within the range of wavelengths in the visible spectrum. The pitch length, and hence the wavelength of the scattered light, is influenced by shear stress normal to the helical axis. This unique optical property produces a measurable color change in response to an applied shearing force. The full-surface shear stress vector measurement method, developed at NASA-Ames, is schematically illustrated. As with the visualization method, the coated test surface is illuminated from the normal direction with white light and the camera is positioned at an above-plane view angle of approximately 30 deg. Experiments have been initiated at NASA Ames to begin the process of quantifying surface-inclination (surface-curvature) effects on shear vector measurement accuracy. In preliminary experiments, surface-inclination angles theta(sub x), theta(sub y) of 0, +/-5, +/-10, and +/-15 deg were employed. In this arrangement, white-light illumination was positioned normal to the untilted test surface, and the camera above-plane view angle was set at 30 deg relative to the untilted test surface. As can be seen, vector-aligned lambda(sub d) values showed no dependence on theta(sub x) or theta(sub y) for absolute values of these tilt angles is less than or equal to 15 deg. Acquisition and analyses of full-surface color images are presently underway to definitively document the insensitivity limits of the shear vector measurement methodology to surface-slope variations.

Description: Experimental results for the spreading and centerline velocity decay rates for round, compressible jets, from a convergent and a convergent-divergent nozzle, are presented. The spreading rate is determined from the variation of streamwise mass flux obtained from Pitot probe surveys. Results for the far asymptotic region show that both spreading and centerline velocity decay rates, when nondimensionalized by parameters at the nozzle exit, decrease with increasing "jet Mach number" M(sub j). Dimensional analysis with the assumption of momentum conservation, together with compressible flow calculations for the conditions at the nozzle exit, predict this Mach number dependence well. The analysis also demonstrates that an increase in the "potential core length" of the jet occurring with increasing M(sub j), a commonly observed trend, is largely accounted for simply by the variations in the density and static pressure at the nozzle exit. The effect of decreasing mixing efficiency with increasing compressibility is inferred to contribute only partially to the latter trend.

Description: The purpose of this paper is to address the problems of Aircraft Engine Technology, Airframe Technology, and Rotorcraft Technology. In addition, several applied problems for domestic industry are also studied using knowledge and experience from Aerospace Sciences. The reduction of aircraft noise is a significant driver in the success of the NASA AST and HSR programs as they attempt to meet stringent international environment regulations on noise for commercial aircraft. In accordance with the project fulfillment under this NASA grant the (Fluid Mechanics and Acoustics Laboratory) FM&AL investigates novel and promising concepts for reduction of noise and improvement of propulsion efficiency in jet exhaust nozzles and fans.

Description: Experimental results for the spreading and centerline velocity decay rates are presented for round, compressible, submerged jets, with total temperature same as ambient temperature, covering a 'jet Mach number' (M(sub J)) range of 0.3 - 2.0. Data from a convergent and a convergent-divergent nozzle are considered. The spreading rate is determined from the total longitudinal mass fluxes measured by Pilot probe survey. Results for the far asymptotic regions show that both spreading and centerline velocity decay rates, when normalized by conditions at the nozzle exit, decrease with increasing M(sub J). Dimensional analysis and isentropic calculations predict this Mach number dependence well. The analysis also demonstrates that an increase in the 'potential core length' of the jet occurring with increasing M(sub J), a commonly observed trend, is largely accounted for simply by the variations in the density and static pressure at the nozzle exit. The effect of decreasing mixing efficiency with increasing compressibility is inferred to contribute only partially to the latter trend.

Description: This paper presents an overview and summary of the many different research work related to tip vortex flows and wake/trailing vortices as applied to practical engineering problems. As a literature survey paper, it outlines relevant analytical, theoretical, experimental and computational study found in literature. It also discusses in brief some of the fundamental aspects of the physics and its complexities. An appendix is also included. The topics included in this paper are: 1) Analytical Vortices; 2) Experimental Studies; 3) Computational Studies; 4) Wake Vortex Control and Management; 5) Wake Modeling; 6) High-Lift Systems; 7) Issues in Numerical Studies; 8) Instabilities; 9) Related Topics; 10) Visualization Tools for Vertical Flows; 11) Further Work Needed; 12) Acknowledgements; 13) References; and 14) Appendix.

Description: Flat plate skin friction calculations over a range of Mach numbers from 0.4 to 3.5 at Reynolds numbers from 16 million to 492 million using a Navier Stokes method with advanced turbulence modeling are compared with incompressible skin friction coefficient correlations. The semi-empirical correlation theories of van Driest; Cope; Winkler and Cha; and Sommer and Short T' are used to transform the predicted skin friction coefficients of solutions using two algebraic Reynolds stress turbulence models in the Navier-Stokes method PAB3D. In general, the predicted skin friction coefficients scaled well with each reference temperature theory though, overall the theory by Sommer and Short appeared to best collapse the predicted coefficients. At the lower Reynolds number 3 to 30 million, both the Girimaji and Shih, Zhu and Lumley turbulence models predicted skin-friction coefficients within 2% of the semi-empirical correlation skin friction coefficients. At the higher Reynolds numbers of 100 to 500 million, the turbulence models by Shih, Zhu and Lumley and Girimaji predicted coefficients that were 6% less and 10% greater, respectively, than the semi-empirical coefficients.

Description: This paper presents an approach for the generation of closed manifold surface triangulations from CAD geometry. CAD parts and assemblies are used in their native format, without translation, and a part's native geometry engine is accessed through a modeler-independent application programming interface (API). In seeking a robust and fully automated procedure, the algorithm is based on a new physical space manifold triangulation technique which was developed to avoid robustness issues associated with poorly conditioned mappings. In addition, this approach avoids the usual ambiguities associated with floating-point predicate evaluation on constructed coordinate geometry in a mapped space, The technique is incremental, so that each new site improves the triangulation by some well defined quality measure. Sites are inserted using a variety of priority queues to ensure that new insertions will address the worst triangles first, As a result of this strategy, the algorithm will return its 'best' mesh for a given (prespecified) number of sites. Alternatively, the algorithm may be allowed to terminate naturally after achieving a prespecified measure of mesh quality. The resulting triangulations are 'CFD-ready' in that: (1) Edges match the underlying part model to within a specified tolerance. (2) Triangles on disjoint surfaces in close proximity have matching length-scales. (3) The algorithm produces a triangulation such that no angle is less than a given angle bound, alpha, or greater than Pi - 2alpha This result also sets bounds on the maximum vertex degree, triangle aspect-ratio and maximum stretching rate for the triangulation. In addition to tile output triangulations for a variety of CAD parts, tile discussion presents related theoretical results which assert the existence of such all angle bound, and demonstrate that maximum bounds of between 25 deg and 30 deg may be achieved in practice.

Description: Experiments were conducted on a 76 deg. swept cylinder to establish the conditions for the attachment-line transition process with, and without, surface roughness in a low-disturbance ("quiet"), Mach number 1.6 flow. Local flow parameters were estimated from pressure measurements. These were in good agreement with predictions from an Euler code (CFL3D) and a boundary layer code (BL3D). Hot-wires and Schlieren photography were used to determine the state of the boundary layer. It was found that, for a near-adiabatic wall condition and a smooth surface, the attachment-line, boundary-layer remained laminar up to the highest attainable Reynolds number (R-bar of 790). Transition under the influence of trip wires was found to depend on wind-tunnel disturbance levels and the onset conditions have been established. Results suggest that current design practice, which is based upon data from conventional ("noisy") tunnels, may be highly conservative.

Description: A parallel CFD code solving 3D time accurate Navier-Stokes equations with multipartitioning parallel Methodology is being developed in collaboration with Ohio State University within the Air Vehicle Directorate, at Wright Patterson Air Force Base. The advantage of the multipartitioning parallel method is that the domain decomposition will not introduce domain boundaries for the implicit operators. A ring structure data communication is employed so that the implicit time accurate method can be implemented for multi-processors with the same accuracy as for the single processor. No sub-iteration is needed at the domain boundaries. The code has been validated for some typical unsteady flows, which include Coutte Flow, flow passing a cylinder. The code now is being employed for a large scale time accurate wall jet transient flow computation. 'ne preliminary results are promising. The mesh has been refined to capture more details of the flow field. The mesh refinement computation is in progress and would be difficult to successfully implement without the parallel computation techniques used. A modified version of the code with more efficient inversion of the diagonalized block matrix is currently being tested.

Description: The state of the art in multidimensional combustor modeling as evidenced by the level of sophistication employed in terms of modeling and numerical accuracy considerations, is also dictated by the available computer memory and turnaround times afforded by present-day computers. With the aim of advancing the current multi-dimensional computational tools used in the design of advanced technology combustors, a solution procedure is developed that combines the novelty of the coupled CFD/spray/scalar Monte Carlo PDF (Probability Density Function) computations on unstructured grids with the ability to run on parallel architectures. In this approach, the mean gas-phase velocity and turbulence fields are determined from a standard turbulence model, the joint composition of species and enthalpy from the solution of a modeled PDF transport equation, and a Lagrangian-based dilute spray model is used for the liquid-phase representation. The gas-turbine combustor flows are often characterized by a complex interaction between various physical processes associated with the interaction between the liquid and gas phases, droplet vaporization, turbulent mixing, heat release associated with chemical kinetics, radiative heat transfer associated with highly absorbing and radiating species, among others. The rate controlling processes often interact with each other at various disparate time 1 and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and liquid phase evaporation in many practical combustion devices.

Description: A spacecraft includes heat-generating payload equipment, and a heat transport system with a cold plate thermally coupled to the equipment and a capillary-wick evaporator, for evaporating coolant liquid to cool the equipment. The coolant vapor is coupled to a condenser and in a loop back to the evaporator. A heated coolant reservoir is coupled to the loop for pressure control. If the wick is not wetted, heat transfer will not begin or continue. A pair of check valves are coupled in the loop, and the heater is cycled for augmentation pumping of coolant to and from the reservoir. This augmentation pumping, in conjunction with the check valves, wets the wick. The wick liquid storage capacity allows the augmentation pump to provide continuous pulsed liquid flow to assure continuous vapor transport and a continuously operating heat transport system. The check valves are of the ball type to assure maximum reliability. However, any type of check valve can be used, including designs which are preloaded in the closed position. The check valve may use any ball or poppet material which resists corrosion. For optimum performance during testing on Earth, the ball or poppet would have neutral buoyancy or be configured in a closed position when the heat transport system is not operating. The ball may be porous to allow passage of coolant vapor.

Description: Microgravity processes must rely on mechanisms other than buoyancy to move bubbles or droplets from one region to another in a continuous liquid phase. One suggested method is thermocapillary migration in which a temperature gradient is applied to the continuous phase. A significant and as yet unresolved impediment to the use of thermocapillary migration to direct bubble or drop motion is that these migrations can be significantly retarded by the adsorption onto the fluid particle surface of surface active impurities unavoidably present in the continuous or (if the particle is a liquid) droplet phases. The focus of our research was to develop a theory for remobilizing fluid particle interfaces retarded by a surfactant impurity in an effort to make more viable the use of thermocapillary migrations for the management of bubbles and drops in microgravity. We postulated that a surfactant at high bulk concentration which kinetically exchanges rapidly with the surface can restore interface mobility. The scaling arguments along with a discussion of the previous literature is reviewed in the context of the scaling framework. The specific objectives of the research were twofold. The first was to prove the remobilization theory by studying a model problem. As the mechanism for remobilization is independent of the force which drives the particle, the fluid particle shape and the presence of fluid inertia, we chose the simplest model consisting of a spherical bubble rising steadily by buoyancy in creeping flow. We solved the hydrodynamic and surfactant transport equations for rapid kinetic exchange to demonstrate that as the concentration increases, the Marangoni retardation at first increases (the low k behavior) and then decreases (the high k behavior). The second objective was to develop a method to determine the kinetic rate constants of a surfactant molecule, since this information is necessary to select surfactants which will exchange rapidly enough relative to the convective rate in the thermocapillary process of interest. To measure the kinetic rate, we measure the dynamic tension change accompanying adsorption onto an initially clean interface, or the re-equilibration in tension when an equilibrium interface is compressed. The dynamic tension measurements are made by a pendant bubble method, in which surfactant adsorbs onto a pendant bubble, and the tension is measured by analyzing the shape change in the bubble. We conclude this report by detailing the publications, presentations and doctoral thesis completed under the auspices of this grant.

Description: This viewgraph presentation gives an overview of the parallel domain decomposition preconditioning for computational fluid dynamics. Details are given on some difficult fluid flow problems, stabilized spatial discretizations, and Newton's method for solving the discretized flow equations. Schur complement domain decomposition is described through basic formulation, simplifying strategies (including iterative subdomain and Schur complement solves, matrix element dropping, localized Schur complement computation, and supersparse computations), and performance evaluation.

Description: This paper reports computational analysis of radiation emission experiments in a high enthalpy arc-jet wind tunnel at NASA Ames Research Center. Recently, as part of ongoing arc-jet characterization work, spectroscopic radiation emission experiments have been conducted at the 20 MW NASA Ames arc-jet facility. The emission measurements were obtained from the arc-jet freestream and from a shock layer formed in front of flatfaced models. Analysis of these data is expected to provide valuable information about the thermodynamic state of the gas in the arc-jet freestream and in the shock layer as well as thermochemical equilibration processes behind the shock in arc-jet flows. Knowledge of the thermodynamic state of the gas in arc-jet test flows and especially within the shock layer is essential to interpret the heat transfer measurements such as in surface catalysis experiments. The present work is a continuation of previous work and focuses on analysis of the emission data obtained at relatively low-pressure conditions for which the arc-jet shock layer is expected to be in thermal and chemical nonequilibrium. Building blocks of the present computational analysis are: (1) simulation of nonequilibrium expanding flow in the converging-diverging conical nozzle and supersonic jet; (2) simulation of nonequilibrium shock layer formed in front of the flat-faced cylinder model; and (3) prediction of line-of-sight radiation from the computed flowfield. For computations of the nonequilibrium flow in the conical nozzle and shock layer, multi-temperature nonequilibrium codes with the axisymmetric formulation are used. For computations of line-of-sight radiation. a nonequilibrium radiation code (NEQAIR) is used to predict emission spectra from the computed flowfield. Computed line-of-sight averaged flow properties such as vibrational and rotational temperatures, species number densities within the shock layer will be compared with those deduced from the experimental spectra. Detailed comparisons of computational and experimental spectra will also be presented.

Description: Much of the ground based testing of advanced thermal protection system (TPS) components for the X33 program is done in arc-heated wind tunnels such as those located in the Arc-Jet Complex at NASA Ames Research Center. These facilities are capable of simulating the high temperature, chemically reacting flow environment experienced by the vehicle during flight. This allows one to test critical design issues such as maximum reuse temperatures, seals, gaps, and increases in heating due to interfaces between different materials. Computational fluid dynamics (CFD) has evolved to the point where it now can be used in the vehicle design process for accurate and timely prediction of trajectory based aerothermal heating environments for re-entry vehicles. It can also be used for simulation of the flow environments in ground based facilities such as arcjets. By utilization of the same CFD code and solution methodology, the important differences between ground test and flight may be quantified. The goal of this paper is to utilize CFD to provide validated simulations of the flow environment in the NASA-Ames semi elliptic nozzle arcjet facilities. The validation of the ground simulations will come From comparison to existing calibration data. Specific tests in support of the X33 TPS test program will ilso be simulated. In this manner, the differences between the ground test simulation and the flight environment can be identified for a measure of ground test to flight traceability.

Description: Researchers from the NASA Lewis Research Center have obtained the first combustion/emissions data under extreme future engine operating conditions. In Lewis' new world-class 60-atm combustor research facility--the Advanced Subsonic Combustion Rig (ASCR)--a flametube was used to conduct combustion experiments in environments as extreme as 900 psia and 3400 F. The greatest challenge for combustion researchers is the uncertainty of the effects of pressure on the formation of nitrogen oxides (NOx). Consequently, U.S. engine manufacturers are using these data to guide their future combustor designs. The flametube's metal housing has an inside diameter of 12 in. and a length of 10.5 in. The flametube can be used with a variety of different flow paths. Each flow path is lined with a high-temperature, castable refractory material (alumina) to minimize heat loss. Upstream of the flametube is the injector section, which has an inside diameter of 13 in. and a length of 0.5-in. It was designed to provide for quick changeovers. This flametube is being used to provide all U.S. engine manufacturers early assessments of advanced combustion concepts at full power conditions prior to engine production. To date, seven concepts from engine manufacturers have been evaluated and improved. This collaborated development can potentially give U.S. engine manufacturers the competitive advantage of being first in the market with advanced low-emission technologies.

Description: Detailed flow field measurements have been carried out for a turbulent circular jet perturbed by tabs and artificial excitation. Two "delta tabs" were placed at the nozzle exit at diametricall opposite y locations. The excitation condition involved subharmonic resonance that manifested in a periodic vortex pairing in the near flow field. While the excitation and the tabs independently increased jet spreading, a combination of the two diminished the effect. The jet spreading was most pronounced with the tabs but was reduced when excitation was applied to the tabbed jet. The tabs generated streamwise vortex pairs that caused a lateral spreading of the jet in a direction perpendicular to the plane containing the tabs. ne excitation, on the other hand, organized the azimuthal vorticity into coherent ring structures whose evolution and pairing also increased entrainment by the jet. In the tabbed case, the excitation produced coherent azimuthal structures that were distorted and asymmetric in shape. The self-induction of these structures produced an effect that opposed the tendency for the lateral spreading of the streamwise vortex pairs. The passage of the distorted vortices, and their pairing, also had a cancellation effect on the time-averaged streamwise vorticity field. These led to the reduction in jet spreading.

Description: The proposed paper will describe the state of the art associated with numerical solution of the full or exact velocity potential equation for solving transonic, external-aerodynamic flows. The presentation will begin with a review of the literature emphasizing research activities of the past decade. Next, the various forms of the full or exact velocity potential equation, the equation's corresponding mathematical characteristics, and the derivation assumptions will be presented and described in detail. Impact of the derivation assumptions on simulation accuracy, especially with respect to shock wave capture, will be presented and discussed relative to the more complete Euler or Navier-Stokes formulations. The technical presentation will continue with a description of recently developed full potential numerical approach characteristics. This description will include governing equation nondimensionalization, physical-to-computational-domain mapping procedures, a limited description of grid generation requirements, the spatial discretization scheme, numerical implementation of boundary conditions, and the iteration scheme. The next portion of the presentation will present and discuss numerical results for several two- and three-dimensional aerodynamic applications. Included in the results section will be a discussion and demonstration of a typical grid refinement analysis for determining spatial convergence of the numerical solution and level of solution accuracy. Computer timings for a variety of full potential applications will be compared and contrasted with similar results for the Euler equation formulation. Finally. the presentation will end with concluding remarks and recommendations for future work.

Description: Small Business Innovation Research contracts from Goddard Space Flight Center to Thermacore Inc. have fostered the company work on devices tagged "heat pipes" for space application. To control the extreme temperature ranges in space, heat pipes are important to spacecraft. The problem was to maintain an 8-watt central processing unit (CPU) at less than 90 C in a notebook computer using no power, with very little space available and without using forced convection. Thermacore's answer was in the design of a powder metal wick that transfers CPU heat from a tightly confined spot to an area near available air flow. The heat pipe technology permits a notebook computer to be operated in any position without loss of performance. Miniature heat pipe technology has successfully been applied, such as in Pentium Processor notebook computers. The company expects its heat pipes to accommodate desktop computers as well. Cellular phones, camcorders, and other hand-held electronics are forsible applications for heat pipes.

Description: This in-house experiment on the gust response of an annular turbine cascade was particularly designed to provide data to compare with the results of a typical, linearized gust-response analysis. Reduced frequency, Mach number, and incidence were varied independently. Except for the lowest reduced frequency, the gust velocity distribution was nearly sinusoidal. For the high inlet-velocity series of tests, the cascade was near choking. The mean flow was documented by measuring blade surface pressures and the cascade exit flow, and high-response pressure transducers were used to measure the unsteady pressure distribution. Inlet-velocity components and turbulence parameters were measured using hot wire anemometry. In addition to the synchronous time-averaged pressure spectra, typical power spectra are included for several representative conditions.

Description: Through a combination of experimental and computational studies, researchers at the NASA Lewis Research Center investigated the effect of upstream blade-row wake passing on the showerhead (leading edge) film cooling of a downstream turbine. The experiments were performed in a steady-flow annular turbine cascade facility equipped with an upstream rotating row of cylindrical rods to produce a periodic wake field similar to that found in an actual turbine. Spanwise, chordwise, and temporal resolution of the blade surface temperature were achieved through the use of an array of nickel thin-film surface gauges covering one unit cell of a showerhead film hole pattern. Film effectiveness and Nusselt numbers were determined for a test matrix of various injectants, injectant blowing ratios, and wake Strouhal numbers (St).

Description: The NASA Lewis Research Center is developing analytical methods and software tools to create a bridge between the controls and computational fluid dynamics (CFD) disciplines. Traditionally, control design engineers have used coarse nonlinear simulations to generate information for the design of new propulsion system controls. However, such traditional methods are not adequate for modeling the propulsion systems of complex, high-speed vehicles like the High Speed Civil Transport. To properly model the relevant flow physics of high-speed propulsion systems, one must use simulations based on CFD methods. Such CFD simulations have become useful tools for engineers that are designing propulsion system components. The analysis techniques and software being developed as part of this effort are an attempt to evolve CFD into a useful tool for control design as well. One major aspect of this research is the generation of linear models from steady-state CFD results. CFD simulations, often used during the design of high-speed inlets, yield high resolution operating point data. Under a NASA grant, the University of Akron has developed analytical techniques and software tools that use these data to generate linear models for control design. The resulting linear models have the same number of states as the original CFD simulation, so they are still very large and computationally cumbersome. Model reduction techniques have been successfully applied to reduce these large linear models by several orders of magnitude without significantly changing the dynamic response. The result is an accurate, easy to use, low-order linear model that takes less time to generate than those generated by traditional means. The development of methods for generating low-order linear models from steady-state CFD is most complete at the one-dimensional level, where software is available to generate models with different kinds of input and output variables. One-dimensional methods have been extended somewhat so that linear models can also be generated from two- and three-dimensional steady-state results. Standard techniques are adequate for reducing the order of one-dimensional CFD-based linear models. However, reduction of linear models based on two- and three-dimensional CFD results is complicated by very sparse, ill-conditioned matrices. Some novel approaches are being investigated to solve this problem.

Description: The entropy associated with absolute equilibrium ensemble theories of ideal, homogeneous, fluid and magneto-fluid turbulence is discussed and the three-dimensional fluid case is examined in detail. A sigma-function is defined, whose minimum value with respect to global parameters is the entropy. A comparison is made between the use of global functions sigma and phase functions H (associated with the development of various H-theorems of ideal turbulence). It is shown that the two approaches are complimentary though conceptually different: H-theorems show that an isolated system tends to equilibrium while sigma-functions allow the demonstration that entropy never decreases when two previously isolated systems are combined. This provides a more complete picture of entropy in the statistical mechanics of ideal fluids.

Description: The year 1971 saw the publication of one of the landmark papers in computational aerodynamics, that of Murman and Cole. As with many seminal works, its significance lies not so much in the specific problem that it addressed| small disturbance, plane transonic flow - but in the identification of a general approach to the solution of a technically important and theoretically difficult problem. The key features of Murman and Cole's work were the use of type- dependent differencing to correctly account for the proper domain of dependence of a mixed elliptic/hyperbolic equation, and the introduction of line relaxation to solve the steady flow equation. All subsequent work in transonic potential flows was based on these concepts. Jameson extended Murman and Cole's ideas to the full potential equation with two important contributions. First, he introduced the rotated difference stencil, which generalized the Murman and Cole type-dependent difference operator to general coordinates. Second, he used the interpretation, introduced by Garabedian, of relaxation as an iteration in artificial time to construct stable relaxation schemes, generalizing the original line relaxation method of Reference. The decade of the 1970s saw an explosion of activity in the solution of transonic potential flows, which has been summarized in the review article of Caughey.

Description: A computational method for accurately predicting the static aeroelastic deformations of typical transonic transport wind tunnel models is described. The method utilizes a finite element method (FEM) for predicting the deformations. Extensive calibration/validation of this method was carried out using a novel wind-off wind tunnel model static loading experiment and wind-on optical wing twist measurements obtained during a recent wind tunnel test in the National Transonic Facility (NTF) at NASA LaRC. Further validations were carried out using a Navier-Stokes computational fluid dynamics (CFD) flow solver to calculate wing pressure distributions about several aeroelastically deformed wings and comparing these predictions with NTF experimental data. Results from this aeroelastic deformation method are in good overall agreement with experimentally measured values. Including the predicted deformations significantly improves the correlation between CFD predicted and experimentally measured wing & pressures.

Description: Boundary-layer bleed is an invaluable tool for controlling the airflow in supersonic aircraft engine inlets. Incoming air is decelerated to subsonic speeds prior to entering the compressor via a series of oblique shocks. The low momentum flow in the boundary layer interacts with these shocks, growing in thickness and, under some conditions, leading to flow separation. To remedy this, bleed holes are strategically located to remove mass from the boundary layer, reducing its thickness and helping to maintain uniform flow to the compressor. The bleed requirements for any inlet design are unique and must be validated by extensive wind tunnel testing to optimize performance and efficiency. To accelerate this process and reduce cost, researchers at the NASA Lewis Research Center initiated an experimental program to study the flow phenomena associated with bleed holes. Knowledge of these flow properties will be incorporated into computational fluid dynamics (CFD) models that will aid engine inlet designers in optimizing bleed configurations before any hardware is fabricated. This ongoing investigation is currently examining two hole geometries, 90 and 20 (both with 5-mm diameters), and various flow features.

Description: A new turbomachinery analysis code called SWIFT has been developed at the NASA Lewis Research Center. SWIFT solves thin-layer Navier-Stokes equations with the Baldwin-Lomax turbulence model and an explicit finite-difference scheme. Preconditioning allows the code to be used for all speed ranges, from incompressible to supersonic flows. Multiblock capability allows three types of grids to be patched together to simulate many types of turbomachinery geometries, including hub and tip clearances. In addition, code users can analyze multistage turbomachinery by using a steady averaging-plane approach. This approach uses linearized characteristic boundary conditions to pass information accurately between the stages.

Description: The Surface Tension Driven Convection Experiment (STDCE) was designed to study basic fluid mechanics and heat transfer on thermocapillary flows generated by temperature variations along the free surfaces of liquids in microgravity. STDCE first flew on the USML-1 mission in July 1992 and was rebuilt for the USML-2 mission that was launched in October 1995. This was a collaborative project with principal investigators from Case Western Reserve University (CWRU), Professors Simon Ostrach and Yasuhiro Kamotani, along with a team from the NASA Lewis Research Center composed of civil servants and contractors from Aerospace Design & Fabrication, Inc. (ADF), Analex, and NYMA, Inc.

Description: The NASA Lewis Research Center is managing a task to numerically simulate overnight, on a parallel computing testbed, the aerodynamic flow in the complete low-pressure subsystem (LPS) of a gas turbine engine. The model solves the three-dimensional Navier- Stokes flow equations through all the components within the LPS, as well as the external flow around the engine nacelle. The LPS modeling task is being performed by Allison Engine Company under the Small Engine Technology contract. The large computer simulation was evaluated on networked computer systems using 8, 16, and 32 processors, with the parallel computing efficiency reaching 75 percent when 16 processors were used.

Description: Renewed interest in hypersonic propulsion systems has led to research programs investigating combined cycle engines that are designed to operate efficiently across the flight regime. The Rocket-Based Combined Cycle Engine is a propulsion system under development at the NASA Lewis Research Center. This engine integrates a high specific impulse, low thrust-to-weight, airbreathing engine with a low-impulse, high thrust-to-weight rocket. From takeoff to Mach 2.5, the engine operates as an air-augmented rocket. At Mach 2.5, the engine becomes a dual-mode ramjet; and beyond Mach 8, the rocket is turned back on. One Rocket-Based Combined Cycle Engine variation known as the "Strut-Jet" concept is being investigated jointly by NASA Lewis, the U.S. Air Force, Gencorp Aerojet, General Applied Science Labs (GASL), and Lockheed Martin Corporation. Work thus far has included wind tunnel experiments and computational fluid dynamics (CFD) investigations with the NPARC code. The CFD method was initiated by modeling the geometry of the Strut-Jet with the GRIDGEN structured grid generator. Grids representing a subscale inlet model and the full-scale demonstrator geometry were constructed. These grids modeled one-half of the symmetric inlet flow path, including the precompression plate, diverter, center duct, side duct, and combustor. After the grid generation, full Navier-Stokes flow simulations were conducted with the NPARC Navier-Stokes code. The Chien low-Reynolds-number k-e turbulence model was employed to simulate the high-speed turbulent flow. Finally, the CFD solutions were postprocessed with a Fortran code. This code provided wall static pressure distributions, pitot pressure distributions, mass flow rates, and internal drag. These results were compared with experimental data from a subscale inlet test for code validation; then they were used to help evaluate the demonstrator engine net thrust.

Description: An experimental and numerical investigation of the thermochemical state of arcjet flows is currently being conducted at NASA Ames Research Center. The experimental approach relies on the use of laser- and emission-spectroscopic diagnostic techniques in three regions of the flow. A fiber optic sensor is used to record spectrally resolved emission signals from the electrode package region, where the flow is most likely to be in thermochemical equilibrium. A second emission diagnostic measurement is made in the shock layer formed over a blunt-body test article placed in the stream, and a CCD camera is used to simultaneously record spectral emission from several measurement locations along the stagnation streamline. Downstream of the nozzle exit, but upstream of the test article, Laser-Induced Fluorescence (LIF) of atomic nitrogen is used to assess the nonequilibrium distribution of flow enthalpy in the free stream. Results from the measurements are compared with predictions from a two-temperature, axisymmetric flow model that solves the nozzle and shock-layer flows.

Description: Implicit methods have been the workhorse for the Euler and Navier-Stokes equations for the last 25 years. The ground breaking work of Dr. Joe Steger in implementing such techniques in practical Euler and Navier-Stokes codes provided the basis for all the success in this area. This presentation will highlight his contribution and technical excellence in the area of implicit methods for CFD.

Description: The use of Direct numerical simulations (DNS) data in developing and testing turbulence models is reviewed. The data is used to test turbulence models at all levels: algebraic, one-equation, two-equation and full Reynolds stress models were tested. Particular examples on the development of models for the dissipation rate equation are presented. Homogeneous flows are used to test new scaling arguments for the various terms in the dissipation rate equation. The channel flow data is used to develop modifications to the equation model that take into account near-wall effects. DNS of compressible flows under mean compression are used in testing new compressible modifications to the two-equation models.

Description: A numerical scheme utilizing a chimera zonal grid approach for solving the three dimensional full potential equation is described. Special emphasis is placed on describing the spatial differencing algorithm around the chimera interface. Results from two spatial discretization variations are presented; one using a hybrid first-order/second-order-accurate scheme and the second using a fully second-order-accurate scheme. The presentation is highlighted with a number of transonic wing flow field computations.

Description: Grid generation issues relating to the simulation of the X33 aerothermal environment using the GASP code are explored. Required grid densities and normal grid stretching are discussed with regards to predicting the fluid dynamic and heating environments with the desired accuracy. The generation of volume grids is explored and includes discussions of structured grid generation packages such as GRIDGEN, GRIDPRO and HYPGEN. Volume grid manipulation techniques for obtaining desired outer boundary and grid clustering using the OUTBOUND code are examined. The generation of the surface grid with the required surface grid with the required surface grid topology is also discussed. Utilizing grids without singular axes is explored as a method of avoiding numerical difficulties at the singular line.

Description: The Internal Fluid Mechanics Division conducts both basic research and technology, and system technology research for aerospace propulsion systems components. The research within the division, which is both computational and experimental, is aimed at improving fundamental understanding of flow physics in inlets, ducts, nozzles, turbomachinery, and combustors. This article and the following three articles highlight some of the work accomplished in 1996. A multidisciplinary combustor design system is critical for optimizing the combustor design process. Such a system should include sophisticated computer-aided design (CAD) tools for geometry creation, advanced mesh generators for creating solid model representations, a common framework for fluid flow and structural analyses, modern postprocessing tools, and parallel processing. The goal of the present effort is to develop some of the enabling technologies and to demonstrate their overall performance in an integrated system called the National Combustion Code.

Description: The receptivity of a laminar boundary layer to an isolated three-dimensional convected disturbance is investigated in a low-speed wind tunnel experiment. The disturbance is created by the short-duration pulsed displacement of a small low-aspect-ratio wing located upstream of a flat plate. The height of the wing is set so that the convected disturbance grazes the edge of the flat-plate boundary layer. A receptivity site is provided by a two-dimensional roughness strip on the surface of the plate. The different propagation speeds of acoustic, convected and instability waves cause the various wave packets from the pulsed displacement to arrive at a downstream measurement station at different times, separating the phenomena and allowing them to be studied independently. Ensemble- averaged measurements are made with and without roughness on the plate. Preliminary analysis of the measurements suggest the presence of a two-dimensional T-S wave packet arising from an interaction between an acoustic wave and the roughness, and a three-dimensional T-S wave packet arising from an interaction between the localized convected disturbance and the roughness strip. The growth rates and spatial characteristics of the disturbances and the instability wave packets are measured as they propagate downstream.

Description: This paper presents two unstructured mesh adaptation schemes for problems in computational fluid dynamics. The procedures allow localized grid refinement and coarsening to efficiently capture aerodynamic flow features of interest. The first procedure is for purely tetrahedral grids; unfortunately, repeated anisotropic adaptation may significantly deteriorate the quality of the mesh. Hexahedral elements, on the other hand, can be subdivided anisotropically without mesh quality problems. Furthermore, hexahedral meshes yield more accurate solutions than their tetrahedral counterparts for the same number of edges. Both the tetrahedral and hexahedral mesh adaptation procedures use edge-based data structures that facilitate efficient subdivision by allowing individual edges to be marked for refinement or coarsening. However, for hexahedral adaptation, pyramids, prisms, and tetrahedra are used as buffer elements between refined and unrefined regions to eliminate hanging vertices. Computational results indicate that the hexahedral adaptation procedure is a viable alternative to adaptive tetrahedral schemes.

Description: This paper reports computational comparisons with experimental studies of a nonequilibrium blunt body shock layer in a high enthalpy arc-jet wind tunnel at NASA Ames Research Center. The primary objective of this work is to investigate the existence of a thermochemical equilibrium region in the shock layer. The existence of such an equilibrium region is of interest for following reasons: (1) to understand the equilibration process behind the shock in an arc-jet flow environment; (2) to interpret measured surface heat transfer data for purpose of determining surface catalytic efficiency, and (3) to determine the total enthalpy from the spectroscopic measurements. The paper will present an analysis of the experimental data obtained in the arc-jet wind tunnel. Experimental data includes measurements of emission spectra of radiation emanating from a shock layer formed in front of a 6-inch flat-faced cylinder. The measurements, obtained using a two dimensional CCD camera mounted on a spectrograph, provide spatially resolved spectra along the stagnation streamline of the model. Computational analysis includes simulation of nonequilibrium flow in the arc-jet facility (flow in the conical nozzle and shock layer in front of a flat-faced cylinder) using 2-D/axisymmetric Navier-Stokes codes and prediction of the radiation spectra from the axisymmetric flowfield using NEQAIR radiation code. Various line-of-sight averaged flow properties such as vibrational and rotational temperatures, species number densities within the shock layer are deduced from the experimental spectra. Comparison of the computed and experimental line-of-sight averaged flow properties provides assessment of thermochemical equilibration processes in an arc-jet shock layer.

Description: The thermal protection system of the windward surface of the X-33 vehicle consists of metallic honeycomb sandwich panels. Thermal gradients experienced during the descent phase of the trajectory result in a different rate of thermal expansion between the inner and outer face sheets of the metallic panels. This causes the panels to bow outward when the temperature of the outer face sheet is larger than that of the inner face sheet and inward when the temperature of the outer face sheet is less than that of he inner face sheet. This results in a quilted-type body surface. Using computational fluid dynamic analysis, this study will determine the effect the metallic TPS panel bowing has on the surface heating.

Description: The usefulness of vortex-fitting in the computational fluid dynamics (CFD) methods to preserve the vortex strength and structure while convecting in a uniform free stream is demonstrated through the numerical simulations of two- and three-dimensional blade-vortex interactions. The fundamental premise of the formulation is the velocity and pressure field of the interacting vortex are unaltered either in the presence of an airfoil or a rotor blade or by the resulting nonlinear interactional flowfield. Although, the governing Euler and Navier-Stokes equations are nonlinear and independent solutions cannot be superposed, the interactional flowfield can be accurately captured by adding and subtracting the flowfield of the convecting vortex at each instant. The aerodynamics and aeroacoustics of two- and three-dimensional blade-vortex interactions have been calculated in Refs. 1-6 using this concept. Some of the results from these publications and similar other published material will be summarized in this paper.

Description: This paper will use high-resolution Navier-Stokes computational fluid dynamics (CFD) simulations to model the near-wake vortex roll-up behind rotor blades. The locations and strengths of the trailing vortices will be determined from newly-developed visualization and analysis software tools applied to the CFD solutions. Computational results for rotor nearwake vortices will be used to study the near-wake vortex roll up for highly-twisted tiltrotor blades. These rotor blades typically have combinations of positive and negative spanwise loading and complex vortex wake interactions. Results of the computational studies will be compared to vortex-lattice wake models that are frequently used in rotorcraft comprehensive codes. Information from these comparisons will be used to improve the rotor wake models in the Tilt-Rotor Acoustic Code (TRAC) portion of NASA's Short Haul Civil Transport program (SHCT). Accurate modeling of the rotor wake is an important part of this program and crucial to the successful design of future civil tiltrotor aircraft. The rotor wake system plays an important role in blade-vortex interaction noise, a major problem for all rotorcraft including tiltrotors.

Description: Knowledge from dynamical systems theory is used to study numerical uncertainties in direct numerical simulation of transition from laminar to turbulent flows. Spurious behavior due to underresolved grids and temporal discretizations is illustrated with 2-D CFD simulations.

Description: The integration of high-fidelity Computational Fluid Dynamics (CFD) analysis tools with the industrial design process benefits greatly from the robust implementations that are transportable across a wide range of computer architectures. In the present work, a hybrid domain-decomposition and parallelization concept was developed and implemented into the widely-used NASA multi-block Computational Fluid Dynamics (CFD) packages implemented in ENSAERO and OVERFLOW. The new parallel solver concept, PENS (Parallel Euler Navier-Stokes Solver), employs both fine and coarse granularity in data partitioning as well as data coalescing to obtain the desired load-balance characteristics on the available computer platforms. This multi-level parallelism implementation itself introduces no changes to the numerical results, hence the original fidelity of the packages are identically preserved. The present implementation uses the Message Passing Interface (MPI) library for interprocessor message passing and memory accessing. By choosing an appropriate combination of the available partitioning and coalescing capabilities only during the execution stage, the PENS solver becomes adaptable to different computer architectures from shared-memory to distributed-memory platforms with varying degrees of parallelism. The PENS implementation on the IBM SP2 distributed memory environment at the NASA Ames Research Center obtains 85 percent scalable parallel performance using fine-grain partitioning of single-block CFD domains using up to 128 wide computational nodes. Multi-block CFD simulations of complete aircraft simulations achieve 75 percent perfect load-balanced executions using data coalescing and the two levels of parallelism. SGI PowerChallenge, SGI Origin 2000, and a cluster of workstations are the other platforms where the robustness of the implementation is tested. The performance behavior on the other computer platforms with a variety of realistic problems will be included as this on-going study progresses.

Description: The widely held view is that transition to turbulence in the Blasius boundary layer occurs via amplification and eventual nonlinear breakdown of initially small amplitude instabilities i.e. Tollmien-Schlichting (TS) waves. However this scenario is only observed for low amplitude free-stream turbulence levels, i.e. u/U 〈 0.1%. Bypass of linear TS instability mechanism occurs for higher EST levels, yet considerable differences exist between the few experiments carefully designed to assess the effect of EST on transition. The consensus is that EST leads to longitudinal streaks that form near the leading edge in the boundary layer . These streaks appeal to be regions of concentrated streamwise vorticity and they are often referred to as Klebanoff modes. The importance of mean flow free-stream nonuniformity (FSN) is not as widely appreciated as EST for characterizing wind tunnel flow quality. Here it is shown that, although the v like generated by a d=50micron wire located upstream of the contraction (Re(sub d)=6.6, x/d=45,000) is immeasurably small by the time it interacts with the leading edge in the test section, it is responsible for generation of a pair of weak streamwise vortices in the boundary layer downstream. The characteristics of these wake-induced vortices and their effect on TS waves are demonstrated. Small remnant FSN variations are also shown to exist downstream of a turbulence grid. The question arises Are the adverse effects introduced by the turbulence grid caused by FST or by small remnant FSN variations?

Description: Challenges in applying CFD to real world problems will be discussed using incompressible flow examples. Despite the remarkable progress made in CFD technology during the past two decades, new challenges are ahead of us. Some of these will be discussed in conjunction with the incompressible Navier-Stokes solver development and their applications. Computed examples include wing tip vortex formation and propagation, flow simulation of an advanced rocket pump, and an extension of the similar technology to biofluid analysis and design. Numerical issues dealing with these applications will be discussed both from algorithm development and from application point of view. The material will be printed in a bound volume after the symposium.

Description: An incompressible flow analysis code, INS3D, has been applied to the development of a mechanical heart assist device. The solution method is based on the artificial compressibility approach and uses an implicit-upwind differencing scheme together with a Gauss-Seidel line relaxation method. The equations are solved in steadily rotating reference frames and the centrifugal and the Coriolis force terms are included as source terms. The resulting computational procedure is validated for liquid rocket engine analysis and applied subsequently to analyze a Ventricular Assist Device (VAD). A new design configuration is developed which includes an inducer upstream of the impeller main blades, and substantial improvement is observed in the performance of the VAD.

Description: Laser Scanning Confocal Microscopy (LSCM) has been used to obtain digital images of the complicated 3-D (three-dimensional) microstructures of rigid, fibrous thermal protection system (TPS) materials. These orthotropic materials are comprised of refractory ceramic fibers with diameters in the range of 1 to 10 microns and have open porosities of 0.8 or more. Algorithms are being constructed to extract quantitative microstructural information from the digital data so that it may be applied to specific heat and mass transport modeling efforts; such information includes, for example, the solid and pore volume fractions, the internal surface area per volume, fiber diameter distributions, and fiber orientation distributions. This type of information is difficult to obtain in general, yet it is directly relevant to many computational efforts which seek to model macroscopic thermophysical phenomena in terms of microscopic mechanisms or interactions. Two such computational efforts for fibrous TPS materials are: i) the calculation of radiative transport properties; ii) the modeling of gas permeabilities.

Description: The talk will present some initial results from the direct numerical simulation (DNS) of compressible turbulent boundary layers. We solve numerically the compressible Navier-Stokes equations using a method based on Spalart's transformation for the incompressible turbulent boundary layer. This allows the spatially developing boundary layer to be transformed to a calculation with periodic boundary conditions in the streamwise and spanwise directions. The equations are solved using Fourier expansions in the horizontal directions and B-splines in the wall-normal direction. The first simulation is at Mach 2.5 with a momentum thickness Reynolds number based on wall viscosity of R(sub theta(sup 1)) = 825. We are examining the physics of the compressible boundary layer using turbulence statistics and budget equations. The turbulence statistics include: rms (root mean square) and mean profiles, energy spectra, and two-point correlations. It is found that there are large density gradients which require significantly more resolution than the incompressible case.

Description: Several stabilized discretization procedures for conservation law equations on triangulated domains will be considered. Specifically, numerical schemes based on upwind finite volume, fluctuation splitting, Galerkin least-squares, and space discontinuous Galerkin discretization will be considered in detail. A standard energy analysis for several of these methods will be given via entropy symmetrization. Next, we will present some relatively new theoretical results concerning congruence relationships for left or right symmetrized equations. These results suggest new variants of existing FV, DG, GLS and FS methods which are computationally more efficient while retaining the pleasant theoretical properties achieved by entropy symmetrization. In addition, the task of Jacobian linearization of these schemes for use in Newton's method is greatly simplified owing to exploitation of exact symmetries which exist in the system. These variants have been implemented in the "ELF" library for which example calculations will be shown. The FV, FS and DG schemes also permit discrete maximum principle analysis and enforcement which greatly adds to the robustness of the methods. Some prevalent limiting strategies will be reviewed. Next, we consider embedding these nonlinear space discretizations into exact and inexact Newton solvers which are preconditioned using a nonoverlapping (Schur complement) domain decomposition technique. Elements of nonoverlapping domain decomposition for elliptic problems will be reviewed followed by the present extension to hyperbolic and elliptic-hyperbolic problems. Other issues of practical relevance such the meshing of geometries, code implementation, turbulence modeling, global convergence, etc. will be addressed as needed.

Description: The Lighthill acoustic analogy, as embodied in the Ffowcs Williams-Hawkings (FW-H) equation, is compared with the Kirchhoff formulation for moving surfaces. A comparison of the two governing equations reveals that the main Kirchhoff advantage (namely nonlinear flow effects are included in the surface integration) is also available to the FW-H method if the integration surface used in the FW-H equation is not assumed impenetrable. The FW-H equation is analytically superior for aeroacoustics because it is based upon the conservation laws of fluid mechanics rather than the wave equation. This means that the FW-H equation is valid even if the integration surface is in the nonlinear region. This is demonstrated numerically in the paper. The Kirchhoff approach can lead to substantial errors if the integration surface is not positioned in the linear region. These errors may be hard to identify. Finally, new metrics based on the Sobolev norm are introduced which may be used to compare input data for both quadrupole noise calculations and Kirchhoff noise predictions.

Description: This paper presents a set of surface and volume grid generation techniques which reuse existing surface and volume grids. These methods use combinations of data manipulations to reduce grid generation time, improve grid characteristics, and increase the capabilities of existing domain discretization software. The manipulation techniques utilize physical and computational domains to produce basis function on which to operate and modify grid character and smooth grids using Trans-Finite Interpolation, a vector interpolation method and parametric re-mapping technique. With these new techniques, inviscid grids can be converted to viscous grids, multiple zone grid adaption can be performed to improve CFD solver efficiency, and topological changes to improve modeling of flow fields can be done simply and quickly. Examples of these capabilities are illustrated as applied to various configurations.

Description: This study attempts t o improve the modeling and computational prediction of high- speed transitional wake flows. The recently developed kappa - zeta (Enstrophy) turbulence model is coupled with a newly developed transition prediction method and implemented in an implicit flow solver well-suited to hypersonic flows. In this model, transition onset is determined as part of the solution. Results obtained using the new model for a 70- deg blunted cone/sting geometry demonstrate better agreement with experimental heat- transfer measurements when compared to laminar calculations as well as solutions using the kappa - omega model. Results are also presented for the situation where transition onset is preselected. It is shown that, in this case, results are quite sensitive to location of the transition point.

Description: A set of compressible flow relations describing flow properties across oblique shock waves, derived for a thermally perfect, calorically imperfect gas, is applied within the existing thermally perfect gas (TPG) computer code. The relations are based upon the specific heat expressed as a polynomial function of temperature. The updated code produces tables of compressible flow properties of oblique shock waves, as well as the original properties of normal shock waves and basic isentropic flow, in a format similar to the tables for normal shock waves found in NACA Rep. 1135. The code results are validated in both the calorically perfect and the calorically imperfect, thermally perfect temperature regimes through comparisons with the theoretical methods of NACA Rep. 1135. The advantages of the TPG code for oblique shock wave calculations, as well as for the properties of isentropic flow and normal shock waves, are its ease of use and its applicability to any type of gas (monatomic, diatomic, triatomic, polyatomic, or any specified mixture thereof).

Description: This paper summarizes the results of a series of experimental studies in the LENS shock tunnel and computations with DSMC and Navier Stokes codes which have been made to examine the aerothermal and flowfield characteristics of the flow over a sting-supported planetary probe configuration in hypervelocity air and nitrogen flows. The experimental program was conducted in the LENS hypervelocity shock tunnel at total enthalpies of 5and 10 MJkg for a range of reservoir pressure conditions from 70 to 500 bars. Heat transfer and pressure measurements were made on the front and rear face of the probe and along the supporting sting. High-speed and single shot schlieren photography were also employed to examine the flow over the model and the time to establish the flow in the base recirculation region. Predictions of the flowfield characteristics and the distributions of heat transfer and pressure were made with DSMC codes for rarefied flow conditions and with the Navier-Stokes solvers for the higher pressure conditions where the flows were assumed to be laminar. Analysis of the time history records from the heat transfer and pressure instrumentation on the face of the probe and in the base region indicated that the base flow was fully established in under 4 milliseconds from flow initiation or between 35 and 50 flow lengths based on base height. The measurements made in three different tunnel entries with two models of identical geometries but with different instrumentation packages, one prepared by NASA Langley and the second prepared by CUBRC, demonstrated good agreement between heat transfer measurements made with two different types of thin film and coaxial gage instrumentation. The measurements of heat transfer and pressure to the front face of the probe were in good agreement with theoretical predictions from both the DSMC and Navier Stokes codes. For the measurements made in low density flows, computations with the DSMC code were found to compare well with the pressure and heat transfer measurements on the sting, although the computed heat transfer rates in the recirculation region did not exhibit the same characteristics as the measurements. For the 10MJkg and 500 bar reservoir match point condition, the measurements and heat transfer along the sting from the first group of studies were in agreement with the Navier Stokes solutions for laminar conditions. A similar set of measurements made in later tests where the model was moved to a slightly different position in the test section indicated that the boundary layer in the reattachment compression region was close to transition or transitional where small changes in the test environment can result in larger than laminar heating rates. The maximum heating coefficients on the sting observed in the present studies was a small fraction of similar measurements obtained at nominally the same conditions in the HEG shock tunnel, where it is possible for transition to occur in the base flow, and in the low enthalpy studies conducted in the NASA Langley high Reynolds number Mach 10 tunnel where the base flow was shown to be turbulent. While the hybrid Navier- StokedDMSC calculations by Gochberg et al. (Reference 1) suggested that employing the Navier- Stokes calculations for the entire flowfield could be seriously in error in the base region for the 10 MJkg, 500 bar test case, similar calculations performed by Cornell, presented here, do not.

Description: The large and small scale vortical motions produced by delta tabs in a two-stream shear layer have been studied experimentally. The tabs were placed at the trailing edge of a splitter plate, mounted in the test section of a wind tunnel, that produced the shear layer. Hot-wire measurements detailed the velocity, vorticity and small scale population distributions in the 3-D space behind the tabs. Small scale structures, which represent the peak in the dissipation spectrum, were identified and counted using the Peak-Valley-Counting technique. The optimal tab inclination angle, theta, and the optimal spanwise spacing, S, for a multiple tab array were identified. The optimal condition was that which maximized the small scale population and generated the strongest streamwise vorticity. Increasing the tab pitch angle from flush to normal to the splitter plate increased the streamwise vorticity as well as the overall small scale population. In fact, the increase in the overall small scale population correlated strongly with the increase in streamwise vorticity intensity. The optimal conditions for peak streamwise vorticity and peak small scale population, with minimum area blockage, were determined to be theta = 45deg and S = 1.5 tab widths.

Description: While there has been much recent progress in simulating nonlinear aeroelastic systems, and in predicting many of the aeroelastic phenomena of concern in transport aircraft design (i.e. transonic flutter buckets), the utility of a simulation in generating an understanding of the flutter behavior is limited. This is due in part to the high cost of generating these simulations; and the implied limitation on the number of conditions that can be analyzed, but there are also some difficulties introduced by the very nature of a simulation. Flutter engineers have traditionally worked in the frequency domain, and are accustomed to describing the flutter behavior of an airplane in terms of its V-G and V-F (or Q-G and Q-F) plots and flutter mode shapes. While the V-G and V-F plots give information about how the dynamic response of an airplane changes as the airspeed is increased, the simulation only gives information about one isolated condition (Mach, airspeed, altitude, etc.). Therefore, where a traditional flutter analysis can let the engineer determine an airspeed at which an airplane becomes unstable, while a simulation only serves as a binary check: either the airplane is fluttering at this condition, or it is not. In this document, a new technique is described in which system identification is used to easily extract modal frequencies and damping ratios from simulation time histories, and shows how the identified parameters can be used to determine the variation in frequency and dampin,o ratio as the airspeed is changed. This technique not only provides the flutter engineer with added insight into the aeroelastic behavior of the airplane, but it allows calculation of flutter mode shapes, and allows estimation of flutter boundaries while minimizing the number of simulations required.

Description: In this lecture we trace the historical developments of alternating direction implicit methods. In particular, we emphasize contributions originating in the Computational Fluid Dynamics Branch at Ames Research Center in the 1970's and early 1980's. Joe Steger played a seminal role in demonstrating the practicality of using an efficient, vectorized, implicit code for solving the compressible Navier-Stokes equations. Numerous discussions with Joe had a significant impact on our own research and it is a pleasure to dedicate this lecture to honor his memory.

Description: We consider preconditioning methods for nonself-adjoint advective-diffusive systems based on a non-overlapping Schur complement procedure for arbitrary triangulated domains. The ultimate goal of this research is to develop scalable preconditioning algorithms for fluid flow discretizations on parallel computing architectures. In our implementation of the Schur complement preconditioning technique, the triangulation is first partitioned into a number of subdomains using the METIS multi-level k-way partitioning code. This partitioning induces a natural 2X2 partitioning of the p.d.e. discretization matrix. By considering various inverse approximations of the 2X2 system, we have developed a family of robust preconditioning techniques. A computer code based on these ideas has been developed and tested on the IBM SP2 and the SGI Power Challenge array using MPI message passing protocol. A number of example CFD calculations will be presented to illustrate and assess various Schur complement approximations.