ALBERT

Description: In the design of a combustor, information is necessary for the mixing of the fuel and air in order to determine the optimum combustor length. In scramjet combustors the mixing often takes place in a shear layer that is formed between the fuel and air. This research was an experimental study of shear layers in supersonic flows aimed at determining what mechanisms affect the shear layer so that the mixing could be better predicted. A second goal was to provide sufficient instream information for use in checking existing Computational Fluid Dynamic (CFD) codes. The shear layer between a supersonic two-dimensional air stream (M = 2 or M = 3) was mixed with a near sonic two-dimensional air stream (M = 1.2). Instream measurements of pitot pressure and cone static pressure were used to determine mean velocity profiles at various axial locations. These velocity profiles were used to determine the shear layer spreading rate and are compared with various predictions. Wall measurements of static pressure, temperature and skin friction were also taken and are presented. The instream measurements were also used for comparison with an existing CFD code. The upstream velocity, pressure and temperature profiles were used as a starting profile and the code was used to calculate downstream profiles for comparison with the experimental results. Reasonable agreement between the measured and calculated results was obtained.

Description: This paper presents a set of second-order closure models for low-Reynolds-number turbulence near the wall. Existing closure models for the Reynolds-stress equations were modified to show proper near-wall behavior. A dissipation-rate equation for the turbulent kinetic energy is also reformulated. The proposed models satisfy realizability and will not produce unphysical behavior. Fully developed channel flows are used for model testing. The equations are solved for the mean velocity, the Reynolds stresses, and the dissipation rate of the turbulent kinetic energy. The calculations are compared with both direct numerical simulations and with measurements. It is shown that the present models perform well in predicting the behavior of the turbulence near a wall. Significant improvements over previous models in predicting the components of the Reynolds stress tensor are obtained in the present models.

Description: The effect of rapid mean compression on compressible turbulence at a range of turbulent Mach numbers is investigated. Rapid dist'ortion theory (RDT) and direct numerical simulation results for the case of axial (one-dimensional) compression are used to illustrate the existence of two distinct rapid compression regimes. These regimes - the nearly solenoidal and the 'pressure-released' - are defined by a single parameter involving the timescales of the mean distortion, the turbulence, and the speed of sound. A general RDT formulation is developed and is proposed as a means of improving turbulence models for compressible flows. In contrast to the well-documented observation that 'compressibility' (measured, for example, by the turbulent Mach number) is often associated with a decrease in the growth rate of turbulent kinetic energy, we find that under rapid distortion compressibility can produce an amplification of the kinetic energy growth rate. We also find that as the compressibility increases, the magnitude of the pressure-dilation correlation increases, in absolute terms, but its relative importance decreases compared to the magnitude of the kinetic energy production.

Description: Hi-Tech Inc., a company which manufactures water jetting equipment, needed a high pressure rotating swivel, but found that available hardware for the system was unsatisfactory. They were assisted by Marshall, which had developed water jetting technology to clean the Space Shuttles. The result was a completely automatic water jetting system which cuts rock and granite and removes concrete. Labor costs have been reduced; dust is suppressed and production has been increased.

Description: Stars and Stripes racing yacht brought the American's Cup back to the United States. Originating from NASA's drag reduction technology, the boats "secret weapon" was that the hull's underside was coated with riblets. Riblets are small, barely visible grooves on the surface of an airplane intended to reduce skin friction by smoothing the turbulent airflow next to the skin. Grooves are V-shaped with the angle pointing in the direction of the airflow. No deeper than a scratch, they have a pronounced beneficial influence on air turbulence. *No longer commercially available.

Description: Developed under a Small Business Innovation Research (SBIR) contract, RAMPANT is a CFD software package for computing flow around complex shapes. The package is flexible, fast and easy to use. It has found a great number of applications, including computation of air flow around a Nordic ski jumper, prediction of flow over an airfoil and computation of the external aerodynamics of motor vehicles.

Description: Macrodyne, Inc.'s laser velocimeter (LV) is a system used in wind tunnel testing of aircraft, missiles and spacecraft employing electro optical techniques to probe the flow field as the tunnel blows air over a model of flight vehicle and to determine velocity of air and its direction at many points around the model. However, current state-of-the-art minicomputers cannot handle the massive flow of real time data from several sources simultaneously. Langley developed instrument Laser Velocimeter Autocovariance Buffer Interface (LVABI). LVABI is interconnecting instrument between LV and computer. It acquires data from as many as six LV channels at high real time data rates, stores it in memory and sends it to computer on command. LVABI has application in variety of research, industrial and defense functions requiring precise flow measurement.

Description: A dual-frequency acoustic levitator containing water was developed for studying bubble and drop dynamics in low gravity. It was flown on USML-1 where it was used in the Glovebox facility. High frequency (21 or 63 kHz) ultrasonic waves were modulated by low frequencies to excite shape oscillations on bubbles and oil drops ultrasonically trapped in the water. Bubble diameters were typically close to 1 cm or larger. When such large bubbles are acoustically trapped on the Earth, the acoustic radiation pressure needed to overcome buoyancy tends to shift the natural frequency for quadrupole (n = 2) oscillations above the prediction of Lamb's equation. In low gravity, a much weaker trapping force was used and measurements of n = 2 and 3 mode frequencies were closer to the ideal case. Other video observations in low gravity include: (i) the transient reappearance of a bulge where a small bubble has coalesced with a large one, (ii) observations of the dynamics of bubbles coated by oil indicating that shape oscillations can shift a coated bubble away from the oil-water interface of the coating giving a centering of the core, and (iii) the agglomeration of bubbles induced by the sound field.

Description: A study of instabilities in incompressible boundary-layer flow on a flat plate is conducted by spatial direct numerical simulation (DNS) of the Navier-Stokes equations. Here, the DNS results are used to critically evaluate the results obtained using parabolized stability equations (PSE) theory and to study mechanisms associated with breakdown from laminar to turbulent flow. Three test cases are considered: two-dimensional Tollmien-Schlichting wave propagation, subharmonic instability breakdown, and oblique-wave break-down. The instability modes predicted by PSE theory are in good quantitative agreement with the DNS results, except a small discrepancy is evident in the mean-flow distortion component of the 2-D test problem. This discrepancy is attributed to far-field boundary- condition differences. Both DNS and PSE theory results show several modal discrepancies when compared with the experiments of subharmonic breakdown. Computations that allow for a small adverse pressure gradient in the basic flow and a variation of the disturbance frequency result in better agreement with the experiments.

Description: A new measurement technique is being developed by NASA to measure off-surface flow fields. This method, Doppler global velocimetry, will allow quantification of complex three-dimensional flow fields at video camera rates. The entire flow field structure within a selected plane is measured simultaneously rather than by scanned, point-by-point measurements using conventional laser velocimetry. Data obtained using this technique will be used to correlate with other data sets for verification, and following verification, provide a quantified, highly detailed definition of the flow field. This will help to improve the understanding of fluid physics, supplement and broaden the database required to validate and refine computational fluid dynamics (CFD) models, and improve aircraft design methodology. To assess the capability of the technique, velocity measurements of the vortical flow field above a thin 75-degree delta wing were made in the NASA - Langley Basic Aerodynamics Research Tunnel. Preliminary comparisons of the results were made with similar measurements obtained using a three component laser velocimeter indicate that this technique is capable of describing the entire three - component velocity flow field simultaneously within a measurement plane in real time.

Description: The spatial evolution of cross flow-vortex packets in a laminar boundary layer on a swept wing are computed by the direct numerical simulation of the incompressible Navier- Stokes equations. A wall-normal velocity distribution of steady suction and blowing at the wing surface is used to generate a strip of equally spaced and periodic disturbances along the span. Three simulations are conducted to study the effect of initial amplitude on the disturbance evolution, to determine the role of traveling cross ow modes in transition, and to devise a correlation function to guide theories of transition prediction. In each simulation, the vortex packets first enter a chordwise region of linear independent growth, then, the individual packets coalesce downstream and interact with adjacent packets, and, finally, the vortex packets nonlinearly interact to generate inflectional velocity profiles. As the initial amplitude of the disturbance is increased, the length of the evolution to breakdown decreases. For this pressure gradient, stationary modes dominate the disturbance evolution. A two-coeffcient function was devised to correlate the simulation results. The coefficients, combined with a single simulation result, provide sufficient information to generate the evolution pattern for disturbances of any initial amplitude.

Description: A new nonintrusive flow diagnostics instrumentation system, Doppler global velocimetry, is presented. The system is capable of making simultaneous, three-component velocity measurements within a selected measurement plane at video camera rates. These velocity images can provide the researcher with spatial and temporal information about the flow field in a global sense. The investigation of a vortical flow above a 75-degree delta wing comparing standard three-component, fringe-type laser velocimetry measurements with Doppler global velocimetry measurements is presented.

Description: This paper documents the development of the National Aeronautics and Space Administration s (NASA) Langley Research Center ( LaRC) Coherent Antistokes Raman Spectroscopy (CARS) systems for measurements of temperature in a turbulent subsonic or supersonic reacting hydrogen-air environment. Spectra data provides temperature data when compared to a precalculated library of nitrogen CARS spectra. Library validity was confirmed by comparing CARS temperatures derived through the library with three different techniques for determination of the temperature in hydrogen-air combustion and an electrically heated furnace. The CARS system has been used to survey temperature profiles in the simulated flow of a supersonic combustion ramjet (scramjet) model. Measurement results will be discussed.

Description: An implicit, Navier-Stokes solution algorithm is presented for the computation of turbulent flow on unstructured grids. The inviscid fluxes are computed using an upwind algorithm and the solution is advanced in time using a backward-Euler time-stepping scheme. At each time step, the linear system of equations is approximately solved with a point-implicit relaxation scheme. This methodology provides a viable and robust algorithm for computing turbulent flows on unstructured meshes. Results are shown for subsonic flow over a NACA 0012 airfoil and for transonic flow over a RAE 2822 airfoil exhibiting a strong upper-surface shock. In addition, results are shown for 3 element and 4 element airfoil configurations. For the calculations, two one equation turbulence models are utilized. For the NACA 0012 airfoil, a pressure distribution and force data are compared with other computational results as well as with experiment. Comparisons of computed pressure distributions and velocity profiles with experimental data are shown for the RAE airfoil and for the 3 element configuration. For the 4 element case, comparisons of surface pressure distributions with experiment are made. In general, the agreement between the computations and the experiment is good.

Description: A technique is presented for triangulation of NURBS surfaces. This technique is built upon an advancing front technique combined with grid point projection. This combined approach has been successfully implemented for structured and unstructured grids.

Description: An approximate method for development of flow and thermal boundary layers in laminar regime on cylinders with arbitrary cross section and transpiration-cooled walls is obtained by use of Karman's integrated momentum equation and an analogous heat-flow equation. Incompressible flow with constant property values throughout boundary layer is assumed. Shape parameters for approximated velocity and temperature profiles and functions necessary for solution of boundary-layer equations are presented as charts, reducing calculations to a minimum. The method is applied to determine local heat-transfer coefficients and surface temperature-cooled turbine blades for a given flow rate. Coolant flow distributions necessary for maintaining uniform blade temperatures are also determined.

Description: The presence of radomes and instruments that are sensitive to water films or ice formations in the nose section of all-weather aircraft and missiles necessitates a knowledge of the droplet impingement characteristics of bodies of revolution. Because it is possible to approximate many of these bodies with an ellipsoid of revolution, droplet trajectories about an ellipsoid of revolution with a fineness ratio of 10 were computed for incompressible axisymmetric air flow. From the computed droplet trajectories, the following impingement characteristics of the ellipsoid surface were obtained and are presented in terms of dimensionless parameters: (1) total rate of water impingement, (2) extent of droplet impingement zone, and (3) local rate of water impingement. These impingement characteristics are compared briefly with those previously reported for an ellipsoid of revolution with a fineness ratio of 5.

Description: An investigation at a free-stream Mach number of 2.02 was made to determine the effects of a propulsive jet on a wing surface located in the vicinity of a choked convergent nozzle. Static-pressure surveys were made on a flat surface that was located in the vicinity of the propulsive jet. The nozzle was operated over a range of exit pressure ratios at different fixed vertical distances from the flat surface. Within the scope of this investigation, it was found that shock waves, formed in the external flow because of the presence of the propulsive jet, impinged on the flat surface and greatly altered the pressure distribution. An integration of this pressure distribution, with the location of the propulsive jet exit varied from 1.450 propulsive-jet exit diameters to 3.392 propulsive-jet exit diameters below the wing, resulted in an incremental lift for all jet locations that was equal to the gross thrust at an exit pressure ratio of 2.86. This incremental lift increased with increase in exit pressure ratio, but not so rapidly as the thrust increased, and was approximately constant at any given exit pressure ratio.

Description: Numerical solutions of the differential equation obtained from the momentum theorem for the development of a turbulent boundary layer along a thermally insulated surface in two-dimensional and in radial shock-free flow are presented in tabular form for a range of Mach numbers from 0.100 to 10. The solution can be used in a step-wise procedure with any given distribution of favorable pressure gradients and for zero pressure gradients. Solutions are also given for use with moderate adverse pressure gradients. The mean velocity in the boundary layer is approximated by a power-law profile. In view of the stepwise integration methods to be used, the exponent designated the profile shape can be varied along the surface between the integral fraction limits 1/5 and 1/11 through interpolation. Agreement obtained between theoretical and experimental boundary-layer development in a supersonic nozzle at a nominal Mach number of 2 indicates the general validity of the approximations used in the analysis - in particular, the method of extrapolating low-speed skin-friction relations to high Mach number flows. The extrapolation method used assumes that the skin-friction coefficient depend primarily on Reynolds number, provided that the density and the kinematic viscosity are evaluated at surface conditions.

Description: The condensation pressure of air was determined over the range of temperature from 60 to 85 K. The experimental results were slightly higher than the calculated values based on the ideal solution law. Heat of vaporization of oxygen was determined at four temperatures ranging from about 68 to 91 K and of nitrogen similarly at four temperatures ranging from 62 to 78 K.

Description: The heat requirements for the icing protection of two radome configurations have been studied over a range of design icing conditions. Both the protection limits of a typical thermal protection system and the relative effects of the various icing variables have been determined. For full evaporation of all impinging water, an effective heat density of 14 watts per square inch was required. When a combination of the evaporation and running wet surface systems was employed, a heat requirement of 5 watts per square inch provided protection at severe icing and operating conditions.

Description: The Navier-Stokes equations of motion and the equation of continuity are transformed so as to apply to an orthogonal curvilinear coordinate system rotating with a uniform angular velocity about an arbitrary axis in space. A usual simplification of these equations as consistent with the accepted boundary-layer theory and an integration of these equations through the boundary layer result in boundary-layer momentum-integral equations for three-dimensional flows that are applicable to either rotating or nonrotating fluid boundaries. These equations are simplified and an approximate solution in closed integral form is obtained for a generalized boundary-layer momentum-loss thickness and flow deflection at the wall in the turbulent case. A numerical evaluation of this solution carried out for data obtained in a curving nonrotating duct shows a fair quantitative agreement with the measures values. The form in which the equations are presented is readily adaptable to cases of steady, three-dimensional, incompressible boundary-layer flow like that over curved ducts or yawed wings; and it also may be used to describe the boundary-layer flow over various rotating surfaces, thus applying to turbomachinery, propellers, and helicopter blades.

Description: An experiment was performed on oscillatory thermocapillary flow in the Glovebox aboard the USML-1 Spacelab which was launched in July, 1992. Cylindrical containers of 1 and 3 em in diameter were used. Silicone oils of 2 and 5 cSt viscosity were the test fluids. The fluid was heated by a cylindrical heater placed along the centerline of the container. The diameter of the heater was 10% of the container diameter. The fluid motion was studied by flow visualization. Although oscillations were observed briefly, bubbles generated in the fluid during the experiment disturbed the flow substantially so that the critical temperature differences could not be determined.

Description: The trajectories of droplets in the air flowing past NACA 65(1)-208 airfoil and an NACA 65(1)-212 airfoil, both at an angle of attack of 4 degrees, were determined. The amount of water in droplet form impinging on the airfoils, the area of droplet impingement, and the rate of droplet impingement per unit area on the airfoil surface affected were calculated from the trajectories and are presented. The amount, extent, and rate of impingement of the NACA 65(1)-208 airfoil are compared with the results for the NACA 65(1)1-212 airfoil. Under similar conditions of operation, the NACA 65(1)-208 airfoil collects less water than the NACA 65(1)-212 airfoil. The extent of impingement on the upper surface of the NACA 65(1)-208 airfoil is much less than on the upper surface of the NACA 65(1)-212 airfoil, but on the lower surface the extents of impingement are about the same.

Description: The field of Computational Fluid Dynamics (CFD) has advanced to the point where it can now be used for many applications in fluid mechanics research and aerospace vehicle design. A few applications being explored at NASA Ames Research Center will be presented and discussed. The examples presented will range in speed from hypersonic to low speed incompressible flow applications. Most of the results will be from numerical solutions of the Navier-Stokes or Euler equations in three space dimensions for general geometry applications. Computational results will be used to highlight the presentation as appropriate. Advances in computational facilities including those associated with NASA's CAS (Computational Aerosciences) Project of the Federal HPCC (High Performance Computing and Communications) Program will be discussed. Finally, opportunities for future research will be presented and discussed. All material will be taken from non-sensitive, previously-published and widely-disseminated work.

Description: Line Integral Convolution (LIC) is a powerful technique for imaging and animating vector fields. We extend the LIC paradigm in three ways: (1) The existing technique is limited to vector fields over a regular Cartesian grid. We extend it to vector fields over parametric surfaces, such as those found in curvilinear grids, used in computational fluid dynamics simulations; (2) Periodic motion filters can be used to animate the flow visualization. When the flow lies on a parametric surface, however, the motion appears misleading. We explain why this problem arises and show how to adjust the LIC algorithm to handle it; (3) We introduce a technique to visualize vector magnitudes as well as vector direction. Cabral and Leedom have suggested a method for variable-speed animation, which is based on varying the frequency of the filter function. We develop a different technique based on kernel phase shifts which we have found to show substantially better results. Our implementation of these algorithms utilizes texture-mapping hardware to run in real time, which allows them to be included in interactive applications.

Description: The proposed paper will present a numerical investigation of the flow characteristics and boundary layer development in the nozzles of high enthalpy shock tunnel facilities used for hypersonic propulsion testing. The computed flow will be validated against existing experimental data. Pitot pressure data obtained at the entrance of the test cabin will be used to validate the numerical simulations. It is necessary to accurately model the facility nozzles in order to characterize the test article flow conditions. Initially the axisymmetric nozzle flow will be computed using a Navier Stokes solver for a range of reservoir conditions. The calculated solutions will be compared and calibrated against available experimental data from the DLR HEG piston-driven shock tunnel and the 16-inch shock tunnel at NASA Ames Research Center. The Reynolds number is assumed to be high enough at the throat that the boundary layer flow is assumed turbulent at this point downstream. The real gas affects will be examined. In high Mach number facilities the boundary layer is thick. Attempts will be made to correlate the boundary layer displacement thickness. The displacement thickness correlation will be used to calibrate the quasi-1D codes NENZF and LSENS in order to provide fast and efficient tools of characterizing the facility nozzles. The calibrated quasi-1D codes will be implemented to study the effects of chemistry and the flow condition variations at the test section due to small variations in the driver gas conditions.

Description: In a numerical flow simulation. it is common to generate several thousand time steps of unsteady (time-dependent) flow data. Each time step may require tens to hundreds of megabytes for disk storage, and the total disk requirement for storing the unsteady flow data may be hundreds of gigabytes. Interactive visualization of unsteady flow data of this magnitude is presently impossible with the current hardware technology. This chapter describes the current approaches for unsteady flow visualization. An effective particle tracing technique for unsteady flow is also described. First, the life cycle of a typical numerical flow simulation is outlined. Several unsteady flow data sets from real-world problems are then given. The current approaches for visualizing unsteady flow are then described. There are many existing systems for flow visualization, and some of them are discussed. Streaklines depict time-varying phenomena that are sometimes difficult or impossible to see with other visualization techniques. The algorithms for computing streaklines are described. Several unsteady flow data sets have been visualized using streaklines, and the results are presented. Finally, some current issues in unsteady flow visualization are discussed.

Description: We have implemented a three-dimensional compressible Navier-Stokes code on the Connection Machine CM-5. The code is set up for implicit time-stepping on single or multiple structured grids. For multiple grids and geometrically complex problems, we follow the 'chimera' approach, where flow data on one zone is interpolated onto another in the region of overlap. We will describe our design philosophy and give some timing results for the current code. A parallel machine like the CM-5 is well-suited for finite-difference methods on structured grids. The regular pattern of connections of a structured mesh maps well onto the architecture of the machine. So the first design choice, finite differences on a structured mesh, is natural. We use centered differences in space, with added artificial dissipation terms. When numerically solving the Navier-Stokes equations, there are liable to be some mesh cells near a solid body that are small in at least one direction. This mesh cell geometry can impose a very severe CFL (Courant-Friedrichs-Lewy) condition on the time step for explicit time-stepping methods. Thus, though explicit time-stepping is well-suited to the architecture of the machine, we have adopted implicit time-stepping. We have further taken the approximate factorization approach. This creates the need to solve large banded linear systems and creates the first possible barrier to an efficient algorithm. To overcome this first possible barrier we have considered two options. The first is just to solve the banded linear systems with data spread over the whole machine, using whatever fast method is available. This option is adequate for solving scalar tridiagonal systems, but for scalar pentadiagonal or block tridiagonal systems it is somewhat slower than desired. The second option is to 'transpose' the flow and geometry variables as part of the time-stepping process: Start with x-lines of data in-processor. Form explicit terms in x, then transpose so y-lines of data are in-processor. Form explicit terms in y, then transpose so z-lines are in processor. Form explicit terms in z, then solve linear systems in the z-direction. Transpose to the y-direction, then solve linear systems in the y-direction. Finally transpose to the x direction and solve linear systems in the x-direction. This strategy avoids inter-processor communication when differencing and solving linear systems, but requires a large amount of communication when doing the transposes. The transpose method is more efficient than the non-transpose strategy when dealing with scalar pentadiagonal or block tridiagonal systems. For handling geometrically complex problems the chimera strategy was adopted. For multiple zone cases we compute on each zone sequentially (using the whole parallel machine), then send the chimera interpolation data to a distributed data structure (array) laid out over the whole machine. This information transfer implies an irregular communication pattern, and is the second possible barrier to an efficient algorithm. We have implemented these ideas on the CM-5 using CMF (Connection Machine Fortran), a data parallel language which combines elements of Fortran 90 and certain extensions, and which bears a strong similarity to High Performance Fortran. We make use of the Connection Machine Scientific Software Library (CMSSL) for the linear solver and array transpose operations.

Description: Steady and unsteady flows for propulsion systems are efficiently simulated by solving the incompressible Navier-Stokes equations. The solution method is based on the pseudo compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Current computations use one equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard kappa - epsilon model equations. The resulting computer code is applied to the flow analysis inside an advanced rocket pump impeller in steadily rotating reference frames. Numerical results are compared with experimental measurements. The effects of exit and shroud cavities with the leak-age flow are investigated. Time-accurate incompressible Navier-Stokes formulation with the overlapped grid scheme capability was evaluated by using MIT flapping foil experiment. The grid dependency, turbulence model effects, and the effect of order of differencing were investigated. Numerical results were compared against experimental data. The resulting procedure were applied to unsteady flapping foil calculations. Two upstream NACA 0025 foils perform high frequency synchronized motion and generate unsteady flow conditions to the downstream larger stationary foil. Comparison between unsteady experimental data and numerical results from two different moving boundary procedures will be presented.

Description: The diffusion controlled binary reaction between initially segregated reactants in a two-dimensional low Mach number mixing layers is studied via numerical simulation. The stoichiometric ratio of the reactants is chosen to be much larger than one, as is typical of hydrocarbon flames in air. This results is a flame that is offset from the main vortical region of the mixing layer. In agreement with experimental observations, the flame remains surprisingly uncontorted during the flow evolution and is not entrained into the mixing layer. The effect of the heat release of the flame on the evolution of the mixing layer is thus similar to the effect of a difference in free-stream density between the two sides of the layer. The resulting baroclinic torque inhibits the familiar rolup and pairing of mixing layer vortices common in constant density flows. This also contributes to the layers inability to entrain the flame. The increase in viscosity caused by the heating of the flame reduces the effective Reynolds number of the flow. But, contrary to what has commonly been suggested, this is not the major reason for the inhibition of the usual large-scale mixing layer structures.

Description: Computational fluid dynamic (CFD) analysis is performed on the Lockheed Lifting Body Single-Stage-to-Orbit vehicle to determine the heat transfer to the vehicle during its descent trajectory. Seven species, chemical nonequilibriurn computations using the GASP code will be completed at several trajectory points to assess the thermal protection requirements of the vehicle. Sophisticated surface boundary conditions including in-depth conduction, catalycity, and a variable temperature wall have been incorporated into the flow solver.

Description: MIT flapping foil experiment was used as a validation case to evaluate the current incompressible Navier-Stokes approach with overlapped grid schemes. Steady-state calculations were carried out for overlapped and patched grids. The grid dependency, turbulence model effects, and the effect of order of differencing were investigated. Numerical results were compared against experimental data. The resulting procedure were applied to unsteady flapping foil calculations. Two upstream NACA 0025 foils perform high-frequency synchronized motion and generate unsteady flow conditions to the downstream larger stationary foil. Comparison between unsteady experimental data and numerical results from two different moving boundary procedures will be presented.

Description: An adiabatic demagnetization refrigerator (ADR) is under development at NASA-Ames Research Center that will operate between 2 K and 10 K and will provide 50 mW of cool ng at 2 K. Gadolinium Gallium Garnet (GGG) is selected as the refrigerant for the ADR, To minimize temperature gradients in the GGG, thick slices of GGG are sandwiched together with strips of high-purity copper in between them. The copper strips are used to exchange heat between the GGG and the 2 K and the 10 K heat switches. The heat transfer across the Cu-GGG interfaces is improved by placing thin foils of' high-purity indium at the interfaces. The heat switches employed in the ADR have no moving parts. The 10 K heat switch is a helium gas-gap heat switch; while, the 2 K heat switch is a He ll-gap heat switch. A switch is on when its gap Is filled with helium and is off' when the gap is emptied. This is accomplished with an activated carbon pump (ACP). The ACP adsorbs helium when cooled and desorbs it when heated. A superconducting magnet capable of providing 9 T at 2 K is used for the ADR cycle. A prototype of this refrigerator has been built and is currently under test. A detailed design of the ADR and preliminary test results performed on the prototype ADR will be presented.

Description: Fluid dynamics of turbomachines are complicated by inherently three dimensional structures such as endwall boundary layers, hub corner separation bubbles and tip-leakage flows. In addition, the relative motion between rotors and stators causes unsteady aerodynamic interactions to occur between blade rows. It is necessary to understand the aerodynamics associated with these interactions in order to design turbomachines that are both light and compact as well as reliable and efficient. An unsteady, three-dimensional, thin-layer, Navier-Stokes zonal algorithm is used to investigate the unsteady aerodynamics of multi-stage turbines and compressors. Relative motion between rotors and stators is made possible by the use of systems of patched and overlaid grids. Time-averaged pressures and pressure envelopes have been computed for several two- and three-dimensional single- and multi-stage configurations. Flow visualizations and computed results are in good agreement with experimental data.

Description: Based on the geometry of Mars Environment Survey (MESUR) Pathfinder aeroshell and an estimated Mars entry trajectory, two-dimensional axisymmetric time dependent calculations have been obtained using GIANTS (Gauss-Siedel Implicit Aerothermodynamic Navier-Stokes code with Thermochemical Surface Conditions) code and CMA (Charring Material Thermal Response and Ablation) Program for heating analysis and heat shield material sizing. These two codes are interfaced using a loosely coupled technique. The flowfield and convective heat transfer coefficients are computed by the GIANTS code with a species balance condition for an ablating surface, and the time dependent in-depth conduction with surface blowing is simulated by the CMA code with a complete surface energy balance condition. In this study, SLA-561V has been selected as heat shield material. The solutions, including the minimum heat shield thicknesses over aeroshell forebody, pyrolysis gas blowing rates, surface heat fluxes and temperature distributions, flowfield, and in-depth temperature history of SLA-561V, are presented and discussed in detail.

Description: Up to today, preconditioning methods on massively parallel systems have faced a major difficulty. The most successful preconditioning methods in terms of accelerating the convergence of the iterative solver such as incomplete LU factorizations are notoriously difficult to implement on parallel machines for two reasons: (1) the actual computation of the preconditioner is not very floating-point intensive, but requires a large amount of unstructured communication, and (2) the application of the preconditioning matrix in the iteration phase (i.e. triangular solves) are difficult to parallelize because of the recursive nature of the computation. Here we present a new approach to preconditioning for very large, sparse, unsymmetric, linear systems, which avoids both difficulties. We explicitly compute an approximate inverse to our original matrix. This new preconditioning matrix can be applied most efficiently for iterative methods on massively parallel machines, since the preconditioning phase involves only a matrix-vector multiplication, with possibly a dense matrix. Furthermore the actual computation of the preconditioning matrix has natural parallelism. For a problem of size n, the preconditioning matrix can be computed by solving n independent small least squares problems. The algorithm and its implementation on the Connection Machine CM-5 are discussed in detail and supported by extensive timings obtained from real problem data.

Description: A nonequilibrium, axisymmetric, Navier-Stokes flow solver with coupled radiation has been developed to use in the design of thermal protection systems for vehicles where radiation effects are important. The present method has been compared with an existing flow and radiation solver and with the Project Fire II experimental data. Very good agreement has been obtained over the entire Fire II trajectory with the experimentally determined values of the stagnation radiation intensity in the .2 to 6.2 eV range and with the total stagnation heating. The agreement was significantly better than previous numerical predictions. The effects of a number of flow models are examined to determine which combination of physical models produces the best agreement with the experimental data. These models include radiation coupling, multi-temperature thermal models, finite-rate chemistry, and a quasi-steady-state or Boltzmann assumption for the calculation of the excited electronic states. Finally, the computational efficiency of the present model is evaluated. The radiation properties model developed for this study is shown to offer significant computational savings compared to existing codes.

Description: Rotational temperatures have been measured in rarefied, nonequilibrium, heated freejet expansions of nitrogen using the electron beam fluorescence technique at the University of California at Berkeley Low Density Wind Tunnel facility. Spectroscopic measurements of the (0,0) band of the first negative system of nitrogen reveal the nonequilibrium behavior in the flowfield upstream of, and through the Mach disk, which forms as the freejet expands into a region of finite back pressure. Results compare well with previous freejet expansion data and computations regarding location of the Mach disk and terminal rotational temperature in the expansion. Measurements are also presented for shock thickness based on the rotational temperature changes in the flow. Thickening shock layers, departures of rotational temperature from equilibrium in the expansion region, and downstream rotational temperature recovery much below that of an isentropic normal shock provide indications of the rarefied, nonequilibrium flow behavior. The data are analyzed to infer constant values of the rotational-relaxation collision number from 2.2 to 6.5 for the various flow conditions. Collision numbers are also calculated in a consistent manner for data from other investigations for which is seen a qualitative increase with increasing temperature. Rotational-relaxation collision numbers are seen as not fully descriptive of the rarefied freejet flows. This may be due to the high degree of nonequilibrium in the flowfields, and/or to the use of a temperature-insensitive rotational-relaxation collision number model in the data analyses.

Description: An earlier proposed constitutive relation for normal stresses originated by random particle fluctuations is used to describe a joint effect of thermal and shear-induced fluctuations on concentrational distributions in suspension flow. Averaged products of components of the fluctuation velocity are evaluated on a basis of the rational mechanics approach combined with a simple kinematic consideration. The equation of momentum conservation of the dispersed phase of a suspension closed with this constitutive relation is applied to unidirectional shear flow in the gravity field and to rotational Couette flow. Coupling of the thermal and shear-induced fluctuations results in that the ability of shear flow to suspend particles has a minimum at a certain particle size, all other things being equal. The developed model provides also for a reasonable explanation of particle distributions observed in Couette flow. The approach based on the consideration of momentum balance for the dispersed phase is proved to lead to an effective equation of convective diffusion of the suspended particles. Coefficients of mutual diffusion due to both thermal and shear-induced fluctuations are drastically different from corresponding self-diffusivities as regards both their scaling and their concentrational dependence.

Description: The dynamic regime of gas injection through a circular plate orifice into an ideally wetting liquid is considered, when successively detached bubbles may be regarded as separate identities. In normal gravity and at relatively low gas flow rates, a growing bubble is modeled as a spherical segment touching the orifice perimeter during the whole time of its evolution. If the flow rate exceeds a certain threshold value, another stage of the detachment process takes place in which an almost spherical gas envelope is connected with the orifice by a nearly cylindrical stem that lengthens as the bubble rises above the plate. The bubble shape resembles then that of a mushroom and the upper envelope continues to grow until the gas supply through the stem is completely cut off. Such a stage is always present under conditions of sufficiently low gravity, irrespective of the flow rate. Two major reasons make for bubble detachment: the buoyancy force and the force due to the momentum inflow into the bubble with the injected gas. The former force dominates the process at normal gravity whereas the second one plays a key role under negligible gravity conditions. It is precisely this fundamental factor that conditions the drastic influence on bubble growth and detachment that changes in gravity are able to cause. The frequency of bubble formation is proportional to and the volume of detached bubbles is independent of the gas flow rate in sufficiently low gravity, while at normal and moderately reduced gravity conditions the first variable slightly decreases and the second one almost linearly increases as the flow rate grows. Effects of other parameters, such as the orifice radius, gas and liquid densities, and surface tension are discussed.

Description: The occurrence of large-scale coherent structures in turbulent free shear flows (especially the planar mixing layer) has been recognized for some time. Indeed, the observation of such structures in mixing layers did much to promote interest in the study of coherent structures in turbulence. It has been widely assumed that the large-scale structures in these flows are responsible for the entrainment of free-stream fluid and the overall growth of the layer, while the small-scale structures provide mixing and dissipation. A model of scalar mixing based on these ideas was proposed for these flows. However, recent experimental and computational evidence suggests that the dominance of the large-scale structures in turbulent mixing layers is not universal. In addition, there is a substantial variation among experiments in several statistical measures of self-similar mixing layers, for example growth rate and velocity variances. To investigate the importance of large-scale structures, several free shear flows (mixing layers and wakes) have been simulated via direct numerical simulation. The simulations are designed to mimic experimental mixing layers in which the splitter plate boundary layers are turbulent. Different levels of two-dimensional forcing are included resulting in large-scale structures of differing strength and importance. These simulations are used to investigate the role of large-scale coherent structures in free shear layers and the effect of these structures on relevant turbulence statistics and scalar mixing. It is found that the statistics and structures in several experiments involving turbulent mixing layers are in better agreement with simulations that do not exhibit dominant large-scale structures than those in which the common mixing layer structures do dominate. It is also found that the level of forcing can have a profound effect on the qualitative and quantitative features of these shear layer, even when they are nominally self-similar.

Description: Computational fluid dynamics (CFD) is beginning to play a major role in the aircraft industry of the United States because of the realization that CFD can be a new and effective design tool and thus could provide a company with a competitive advantage. It is also playing a significant role in research institutions, both governmental and academic, as a tool for researching new fluid physics, as well as supplementing and complementing experimental testing. In this presentation, some of the progress made to date in CFD at NASA Ames will be reviewed. The presentation addresses the status of CFD in terms of methods, examples of CFD solutions, and computer technology. In addition, the role CFD will play in supporting the revolutionary goals set forth by the Aeronautical Policy Review Committee established by the Office of Science and Technology Policy is noted. The need for validated CFD tools is also briefly discussed.

Description: Large Dewars often use aluminum radiation shields and stainless steel vent lines. A simple, low cost method of making thermal contact between the shield and the line is to deform the shield around the line. A knowledge of the thermal conductance of such a joint is needed to thermally analyze the system. The thermal conductance of pressed metal contacts consisting of one aluminum and one stainless steel contact has been measured at 77 K, with applied forces from 8.9 N to 267 N. Both 5052 or 5083 aluminum were used as the upper contact. The lower contact was 304L stainless steel. The thermal conductance was found to be linear in temperature over the narrow temperature range of measurement. As the force was increased, the thermal conductance ranged from roughly 9 to 21 mW/K within a range of errors from 3% to 8%. Within the range of error no difference could be found between the using either of the aluminum alloys as the upper contact. Extrapolating the data to zero applied force does not result in zero thermal conductance. Possible causes of this anomalous effect are discussed.

Description: Compressibility effects on turbulent transport of a passive scalar are studied within homogeneous turbulence using a kinematic decomposition of the velocity field into solenoidal and dilatational parts. It is found that the dilatational velocity does not produce a passive scalar flux, and that all of the passive scalar flux is due to the solenoidal velocity.

Description: A conservative zoning technique, wherein the flow field for a finite-difference calculation is divided into several regions to simplify grid generation, is discussed and is applied in the solution of a two-dimensional problem of complex topology. Calculations are performed on two zonal, or patched, grid systems for the supersonic flow over a double-airfoil configuration. The solution is smooth and continuous across the zonal interfaces, and shock waves pass through the boundaries without distortion. In addition, the time accuracy of the zonal boundary method is verified by a two-zone cylinder calculation with a stationary inner and a rotating outer mesh. The feasibility of the zonal approach for use in the solution of geometrically complex and unsteady problems is thus demonstrated.

Description: The index of refraction can considerably influence the temperature distribution and radiative heat flow in semitransparent materials such as some ceramics. For external radiant heating, the refractive index influences the amount of energy transmitted into the interior of the material. Emission within a material depends on the square of its refractive index, and hence this emission can be many times that for a biackbody radiating into a vacuum. Since radiation exiting through an interface into a vacuum cannot exceed that of a blackbody, there is extensive reflection at the internal surface of an interface, mostly by total internal reflection. This redistributes energy within the layer and tends to make its temperature distribution more uniform. The purpose of the present analysis is to show that, for radiative equilibrium in a gray layer with diffuse interfaces, the temperature distribution and radiative heat flux for any index of refraction can be obtained very simply from the results for an index of refraction of unity. For the situation studied here, the layer is subjected to external radiative heating incident on each of its surfaces. The material emits, absorbs, and isotropically scatters radiation. For simplicity the index of refraction is unity in the medium surrounding the layer. The surfaces of the layer are assumed diffuse. This is probably a reasonable approximation for a ceramic layer that has not been polished. When transmitted radiation or radiation emitted from the interior reaches the inner surface of an interface, the radiation is diffused and some of it thereby placed into angular directions for which there is total internal reflection. This provides a trapping effect for retaining energy within the layer and tends to equalize its temperature distribution. An analysis of temperature distributions in absorbing-emitting layers, including index of refraction effects, was developed by Gardon (1958) to predict cooling and heat treating of glass plates. The interfaces were optically smooth; the resulting specular reflections were computed from the Fresnel reflection laws. This provides a somewhat different behavior than for diffuse interfaces. A similar application was for heating that occurs in a window of a re-entry vehicle (Fowle et al., 1969). A number of recent papers (Rokhsaz and Dougherty, 1989; Ping and Lallemand, 1989; Crosbie and Shieh, 1990) further examined the effects of Fresnel boundary reflections and nonunity refractive index. Other examples of analyses of both steady and transient heat transfer to single or multiple plane layers (Amlin and Korpela, 1979; Tarshis et al., 1969) have used diffuse assumptions at the interfaces as in the present study

Description: An investigation of forced-convection heat transfer and associated pressure drops was conducted with air flowing through electrically heated Inconel tubes having various degrees of square-thread-type roughness, an inside diameter of 1/2 inch, and a length of 24 inches. were obtained for tubes having conventional roughness ratios (height of thread/radius of tube) of 0 (smooth tube), 0.016, 0.025, and 0.037 over ranges of bulk Reynolds numbers up to 350,000, average inside-tube-wall temperatures up to 1950deg R, and heat-flux densities up to 115,000 Btu per hour per square foot. Data The experimental data showed that both heat transfer and friction increased with increase in surface roughness, becoming more pronounced with increase in Reynolds number; for a given roughness, both heat transfer and friction were also influenced by the tube wall-to-bulk temperature ratio. Good correlation of the heat-transfer data for all the tubes investigated was obtained by use of a modification of the conventional Nusselt correlation parameters wherein the mass velocity in the Reynolds number was replaced by the product of air density evaluated at the average film temperature and the so-called friction velocity; in addition, the physical properties of air were evaluated at the average film temperature. The isothermal friction data for the rough tubes, when plotted in the conventional manner, resulted in curves similar to those obtained by other investigators; that is, the curve for a given roughness breaks away from the Blasius line (representing turbulent flow in smooth tubes) at some value of Reynolds number, which decreases with increase in surface roughness, and then becomes a horizontal line (friction coefficient independent of Reynolds number). A comparison of the friction data for the rough tubes used herein indicated that the conventional roughness ratio is not an adequate measure of relative roughness for tubes having a square-thread-type element. The present data, as well as those of other investigators, were used to isolate the influence of ratios of thread height to width, thread spacing to width, and the conventional roughness ratio on the friction coefficient. A fair correlation of the friction data was obtained for each tube with heat addition when the friction coefficient and Reynolds number were defined on the basis of film properties; however, the data for each tube retained the curve characteristic of that particular roughness. The friction data for all the rough tubes could be represented by a single line for the complete turbulence region by incorporating a roughness parameter in the film correlation. No correlation was obtained for the region of incomplete turbulence.

Description: Research was conducted to determine the effect of the electrode parameters of spacing, configuration, and material' on the energy required for ignition of a flowing propane-air mixture. In addition, the data were used to indicate the energy distribution along the spark length and to confirm previous observations concerning the effect of spark duration on ignition energy requirements. The data were obtained with a mixture at a fuel-air ratio of 0.0835 (by weight), a pressure of 3 inches of mercury absolute, a temperature of 80 F, and a mixture velocity of 5 feet per second. Results showed that the energy required for ignition decreased as the electrode spacing was increased; a minimum energy occurred at. a spacing of 0.65 inch for large electrodes. For small electrodes, the spacing for minimum energy was not sharply defined. Small-diameter electrodes required less energy than large-diameter electrodes if the spacing was less than the optimum distance of 0.65 inch; at a spacing equal to the optimum distance, no difference was noted. Significant effects of electrode material on ignition energy were ascribed to differences in the type of spark discharges produced; glow discharges required higher energy than the arc-glow discharges. With pure glow discharges, the ignition energy was substantially constant for lead, cadmium, brass, aluminum, and tungsten electrodes. A method is described for determining the energy distribution along a glow discharge. It was found that one-third to one-half of the energy in the spark was concentrated in a small region near the cathode electrode, and the remainder was uniformly distributed across the spark gap. It was impossible to ascertain the dependence of ignition on. this distribution. It was also observed that long-duration (600 microsec) sparks required much less energy for ignition than did short-duration (1 microsec) sparks.

Description: No abstract available

Description: No abstract available

Description: This lecture attempts to illustrate the basic ideas of how the recent advances in nonlinear dynamical systems theory (dynamics) can provide new insights into the understanding of numerical algorithms used in solving nonlinear differential equations (DEs). Examples will be given of the use of dynamics to explain unusual phenomena that occur in numerics. The inadequacy of the use of linearized analysis for the understanding of long time behavior of nonlinear problems will be illustrated, and the role of dynamics in studying the nonlinear stability, accuracy, convergence property and efficiency of using time- dependent approaches to obtaining steady-state numerical solutions in computational fluid dynamics (CFD) will briefly be explained.

Description: A collaborative team of researchers from fields of Computational Fluid Dynamics (CFD), fluid physics, computer architectures, and computer science and knowledge engineering have begun work on a prototype system that addresses several of industry's concerns in using NASA-developed CFD codes as part of the design cycle. A major problem exists in the application of CFD technologies within the aeronautics design cycle due primarily to misunderstandings in the ranges of applicability of the various solver codes or turbulence models. Features that arise during the CFD solution process need to be discriminated and recognized as actual flow features with physical support in the geometry and flow conditions of the problem being solved, or as numerical or non-physical errors arising from mis-application of solver code and its parameters, gridding strategies, or discretization. interpolations. The fundamental concept is to develop an intelligent computational system that can accept the engineer's definition of the problem and construct an optimal CFD solution. To do this requires capturing both the knowledge of how to apply the various CFD tools and how to adapt the application of those tools to flow structures as they evolve during the flow simulation. Embedded within this adaptive system approach is the additional desire to automatically identify and quantify the quality of resolution of the pertinent flow structures, be they genuine or error-induced, and then to adjust the solution strategy accordingly. This paper discusses the status of that prototyping effort.

Description: We present views and analysis of the execution of several PVM codes for Computational Fluid Dynamics on a network of Sparcstations, including (a) NAS Parallel benchmarks CG and MG (White, Alund and Sunderam 1993); (b) a multi-partitioning algorithm for NAS Parallel Benchmark SP (Wijngaart 1993); and (c) an overset grid flowsolver (Smith 1993). These views and analysis were obtained using our Automated Instrumentation and Monitoring System (AIMS) version 3.0, a toolkit for debugging the performance of PVM programs. We will describe the architecture, operation and application of AIMS. The AIMS toolkit contains (a) Xinstrument, which can automatically instrument various computational and communication constructs in message-passing parallel programs; (b) Monitor, a library of run-time trace-collection routines; (c) VK (Visual Kernel), an execution-animation tool with source-code clickback; and (d) Tally, a tool for statistical analysis of execution profiles. Currently, Xinstrument can handle C and Fortran77 programs using PVM 3.2.x; Monitor has been implemented and tested on Sun 4 systems running SunOS 4.1.2; and VK uses X11R5 and Motif 1.2. Data and views obtained using AIMS clearly illustrate several characteristic features of executing parallel programs on networked workstations: (a) the impact of long message latencies; (b) the impact of multiprogramming overheads and associated load imbalance; (c) cache and virtual-memory effects; and (4significant skews between workstation clocks. Interestingly, AIMS can compensate for constant skew (zero drift) by calibrating the skew between a parent and its spawned children. In addition, AIMS' skew-compensation algorithm can adjust timestamps in a way that eliminates physically impossible communications (e.g., messages going backwards in time). Our current efforts are directed toward creating new views to explain the observed performance of PVM programs. Some of the features planned for the near future include: (a) ConfigView, showing the physical topology of the virtual machine, inferred using specially formatted IP (Internet Protocol) packets; and (b) LoadView, synchronous animation of PVM-program execution and resource-utilization patterns.

Description: It is well known that slender bodies of revolution will develop an asymmetric, unsteady flow pattern in experimental tests, if the angle of incidence to the oncoming stream is above a critical value. It has been suggested that the origin of these asymmetric flows may stem from geometric imperfections of the model being tested, or from disturbances in the oncoming stream. In numerical simulations, it is possible to generate bodies of revolution which are perfectly symmetric about their longitudinal axis, and to impose uniform flow conditions which are free from disturbances. The current work presents numerical simulations of the flow about an ogive-cylinder configuration at 40 and 60 degree angle of incidence. These simulations. were performed using numerical algorithms which are also symmetric about the lateral plane of the cylinder body. The flowfields at 40 degree angle of attack were seen to remain symmetric to the round-off accuracy of the computer. At 60 degree angle of attack, a lateral force coefficient developed of O(1) which progressed to an alternate vortex shedding in time. The nature of this lateral force generation and vortex shedding was dependent on the choice of numerical algorithm. The origin of the asymmetries observed in the 60 degree angle of attack computations were traced to round-off errors in the implicit block-matrix inverter. A means of inverting the implicit operator matrices, which maintains the symmetry of the overall numerical algorithm was implemented.

Description: The goal of this work is to add insight about the flow within expansion tubes by using computational fluid dynamics. This is accomplished by comparing the results of axisymmetric numerical simulations with finite-rate chemistry to data from the HYPULSE expansion tube facility which was previously the NASA Langley expansion tube. The numerical simulations begin at the opening of the primary diaphragm and compute the flow throughout the whole facility and, thus, are able to follow and assess the effect of many of the flow features created during operation of the facility. One particular issue that will be investigated is the effect of boundary layer formation in the acceleration tube on the test gas volume and test gas conditions. Both laminar and turbulent boundary layers will be implemented. The effect of momentary shock reflection off the secondary diaphragm will also be investigated. There is concern that such a reflection will stagnate the test gas and create high levels of dissociated molecules. This is particularly important in propulsion experiments where a freestream composition different from flight conditions may influence ignition and burning data. Several different models of diaphragm rupture will be implemented in order to help understand the importance of this issue.

Description: Detailed experimental data have been obtained on several advanced TUFI systems during their exposure to high temperature high shear supersonic turbulent flow. Dimensional stability of these systems were determined at surface temperatures above 3000 F. Effect of step gap configuration on the thermal performance of the TUFI systems were also evaluated.

Description: An alternative theoretical model of joint filtration flow of immiscible incompressible fluids is presented. The model takes into account relaxation processes due to the interchange of the fluids between pores of difference sizes which is driven by capillary forces. The fluids occupy connected regions in a four-dimensional space formed by three coordinates and the pore length scale. When the fluid exchange between pores of given sizes is effected by way of successive flow through pores of all the intermediate sizes, the pressure within each region is governed by a hyperbolic equation, the role of time being played by the pore linear scale. Pressure jumps across hypersurfaces separating the regions equal corresponding values of the capillary pressure. A supplementary condition at any such hypersurface requires the speed of its displacement in the four-dimensional space to coincide with the normal velocity components of both the adjoining fluids. As a result, a principally new statement of multiphase filtration flow problems is gained with allowance for capillary relaxation in the porous space.

Description: The existing pseudocompressibility method for the system of incompressible Navier-Stokes equations is extended to heat transfer problems by including the energy equation. The solution method is based on the pseudo compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Current computations use one-equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard k-epsilon model equations. Both forced and natural convection problems are examined. Numerical results from turbulent reattaching flow behind a backward-facing step will be compared against experimental measurements for the forced convection case. The validity of Boussinesq approximation to simplify the buoyancy force term will be investigated. The natural convective flow structure generated by heat transfer in a vertical rectangular cavity will be studied. The numerical results will be compared by experimental measurements by Morrison and Tran.

Description: We present views and analysis of the execution of several PVM (Parallel Virtual Machine) codes for Computational Fluid Dynamics on a networks of Sparcstations, including: (1) NAS Parallel Benchmarks CG and MG; (2) a multi-partitioning algorithm for NAS Parallel Benchmark SP; and (3) an overset grid flowsolver. These views and analysis were obtained using our Automated Instrumentation and Monitoring System (AIMS) version 3.0, a toolkit for debugging the performance of PVM programs. We will describe the architecture, operation and application of AIMS. The AIMS toolkit contains: (1) Xinstrument, which can automatically instrument various computational and communication constructs in message-passing parallel programs; (2) Monitor, a library of runtime trace-collection routines; (3) VK (Visual Kernel), an execution-animation tool with source-code clickback; and (4) Tally, a tool for statistical analysis of execution profiles. Currently, Xinstrument can handle C and Fortran 77 programs using PVM 3.2.x; Monitor has been implemented and tested on Sun 4 systems running SunOS 4.1.2; and VK uses XIIR5 and Motif 1.2. Data and views obtained using AIMS clearly illustrate several characteristic features of executing parallel programs on networked workstations: (1) the impact of long message latencies; (2) the impact of multiprogramming overheads and associated load imbalance; (3) cache and virtual-memory effects; and (4) significant skews between workstation clocks. Interestingly, AIMS can compensate for constant skew (zero drift) by calibrating the skew between a parent and its spawned children. In addition, AIMS' skew-compensation algorithm can adjust timestamps in a way that eliminates physically impossible communications (e.g., messages going backwards in time). Our current efforts are directed toward creating new views to explain the observed performance of PVM programs. Some of the features planned for the near future include: (1) ConfigView, showing the physical topology of the virtual machine, inferred using specially formatted IP (Internet Protocol) packets: and (2) LoadView, synchronous animation of PVM-program execution and resource-utilization patterns.

Description: FlowViz is a flow visualization application that uses Line Integral Convolution and the texture mapping capabilities of a graphics workstation to create an animation of flow over a curvilinear grid surface.

Description: This paper discusses the development of a general method for the determination of very low leak rates from limiting enclosures. There are many methods that can be used to detect and repair leaks from enclosures. Many methods have also been proposed that allow the estimation of actual leak rates, usually expressed as enclosure volume turnover. The proposed method combines measurements of the state variables (pressure, temperature, and volume) as well as the change in the concentration of a tracer gas to estimate the leak rate. The method was applied to the containment enclosure of the Engineering Development Unit of the CELSS Test Facility, currently undergoing testing at the NASA Ames Research Center.

Description: The case of isotropic compressible turbulence subjected to rapid isotropic compression is studied using inviscid rapid distortion theory and direct numerical simulation. An exact solution to the rapid distortion problem is given, and results are compared to those of direct numerical simulation. Implications for modelling turbulent flows are discussed.

Description: Combining multiple engineering workstations into a network-based heterogeneous parallel computer allows application of aerodynamic optimization with advance computational fluid dynamics codes, which is computationally expensive in mainframe supercomputer. This paper introduces a nonlinear quasi-Newton optimizer designed for this network-based heterogeneous parallel computer on a software called Parallel Virtual Machine. This paper will introduce the methodology behind coupling a Parabolized Navier-Stokes flow solver to the nonlinear optimizer. This parallel optimization package has been applied to reduce the wave drag of a body of revolution and a wing/body configuration with results of 5% to 6% drag reduction.

Description: In order to gain insights into the strong dependence of numerical solutions on initial data for finite time steps, a set of nonlinear test problems rich enough to capture the behavior of difference schemes were recently identified and the numerical basins of attraction for these problems were computed using commonly used time discretizations in CFD. Our study revealed a wealth of surprisingly nonlinear behavior of numerical schemes that were not observed before, in particular for the implicit time discretizations that are commonly used in CFD. The goal of this work is to apply these tools to study a practical model from non-equilibrium flowfield relaxation. This type of problem arises in chemically nonequilibrium hypersonic flows such as in a shock tube experiment or an expanding nozzle. Here we consider a reacting mixture of (N2, N) for an inviscid one-dimensional steady model. Preliminary numerical results indicate that, aside from the possibility of spurious numerical solutions being introduced by the time discretizations, limitations on the model for physical or accurate solutions may also play a part in the dynamics observed.

Description: The present paper addresses some topical issues in modeling compressible turbulent shear flows. The work is based on direct numerical simulation of two supersonic fully developed channel flows between very cold isothermal walls. Detailed decomposition and analysis of terms appearing in the momentum and energy equations are presented. The simulation results are used to provide insights into differences between conventional time-and Favre-averaging of the mean-flow and turbulent quantities. Study of the turbulence energy budget for the two cases shows that the compressibility effects due to turbulent density and pressure fluctuations are insignificant. In particular, the dilatational dissipation and the mean product of the pressure and dilatation fluctuations are very small, contrary to the results of simulations for sheared homogeneous compressible turbulence and to recent proposals for models for general compressible turbulent flows. This provides a possible explanation of why the Van Driest density-weighted transformation is so successful in correlating compressible boundary layer data. Finally, it is found that the DNS data do not support the strong Reynolds analogy. A more general representation of the analogy is analysed and shown to match the DNS data very well.

Description: A shear-stress--sensitive liquid crystal coating (LCC) was used to visualize the surface shear stress distribution on the flat sidewall of the MSU quiet supersonic tunnel as a function of flow stagnation pressure. Under conditions of quiet operation, the LCC color-change response indicated the existence of a nonuniform surface shear stress distribution. This shear pattern was characterized by an elongated, down st ream-point ing triangular region of relatively low shear, with its apex on and its axis coincident with, the sidewall centerline. This low-shear zone was bounded symmetrically by two regions of relatively higher shear; these high-shear zones originated within the concave-curvature portion of the nozzle, in the corners between the flat sidewall and the contoured upper and lower nozzle surfaces. A 3-D Navier-Stokes code was used to compute the pressure and surface shear distributions on the sidewall. Flow-expansion-induced transverse pressure gradients on the nozzle sidewall generated symmetric inflows from the corners towards the sidewall centerline; these inflows caused a thickening of the sidewall boundary layer along the centerline, resulting in lower shear stresses consistent with the liquid crystal results. Peripherally nonuniform laminar boundary layer development, and the associated stability of such complex 3-D flows, must be considered in quiet-tunnel applications using rectangular nozzles. A color video will be shown.

Description: This study was conducted to experimentally characterize the flow field created by the interaction of a single-expansion-ramp-nozzle (SERN) flow with a hypersonic external stream Data were obtained from a generic nozzle/afterbody model in the 3.5-Foot Hypersonic Wind Tunnel of the NASA Ames Research Center in a cooperative experimental program involving Ames and the McDonnell Douglas Aerospace. The model design and test planning were performed in close cooperation with members of the Ames computational fluid dynamics (CFD) team for the National Aero-Space Plane (NASP) program. This paper presents experimental results consisting of oil-flow and shadowgraph flow-visualization photographs, afterbody surface-pressure distributions, boundary-layer rake measurements, and Preston-tube skin-friction measurements.

Description: Experiments were performed to study the evolution of the heat transfer structure in a separated free shear layer region of an incompressible separated turbulent boundary layer flow behind a backward-facing step. While there is an abundance of velocity field measurements of separated flows, heat transfer measurements are rather scarce, thus limiting assessment of the heat transfer physics and its accurate modeling. The purpose of the paper is twofold: to improve an understanding of effects of flow separation on heat transfer characteristics, and to provide data for turbulence modeling and computation. The boundary layer upstream of the step was turbulent and fully developed. A constant temperature surface boundary condition was imposed upstream and downstream of the step for the heat transfer study. An internal mixing-layer like flow forms and grows from the step lip within the original boundary layer. The turbulent structure of the flow evolving downstream, however, does not switch immediately to that of a mixing layer over the entire shear layer thickness. Measurements of mean and fluctuating velocity and temperature fields indicate that the internal layer spreads gradually in the transverse direction while the outer part of the original boundary layer is effectively unperturbed. The results in this paper have not been previously reported.

Description: Direct simulation Monte Carlo (DSMC) calculations of rarefied flows about entry bodies typically employ a fixed surface temperature or a radiative-equilibrium energy balance to compute that temperature. Such boundary conditions neglect any effects of heat capacitance and heat conduction in the spacecraft heat shield and, therefore, provide an upper bound for the surface temperature. Such calculations also neglect pyrolysis from the heat shield which can be significant for a high-energy incident flow at very low densities. Accurate prediction of both heating and aerodynamic forces requires including pyrolysis and surface heat transfer in the models for gas-surface interaction employed in DSMC methods. Although these physical models have long appeared in various continuum flow calculation codes, they have only recently appeared in DSMC codes which are required to simulate rarefied flows during entry at high altitudes. In the current implementation, routines from the widely distributed Charring Material Thermal Response and Ablation (CMA) program are coupled into a DSMC code to calculate the one-dimensional heat transfer into the carbon phenolic heat shield at each point on a vehicle surface. Temperature-dependent material properties, surface re-radiation, and in-depth pyrolysis were included in the calculation, but surface ablation was neglected. Sample calculations for entry of the Galileo probe into the atmosphere of Jupiter demonstrate that including pyrolysis in the model leads to significant differences in predicted aerodynamics. Granted, the drag coefficient does not depend strongly on the surface temperature which can itself be significantly below the radiative equilibrium value during entry. However, the surface mass flux due to pyrolysis of the material is significant once the probe drops to altitudes characterized by transition flow. This leads to a noticeable increase in drag and a decrease in heating compared to a body without pyrolysis.

Description: A critical assessment and compilation of data are presented on attached hypersonic turbulent boundary layers in pressure gradients and compressible turbulent mixing layers. Extensive searches were conducted to identify candidate experiments, which were subjected to a rigorous set of acceptance criteria. Accepted datasets are both tabulated and provided in machine-readable form. The purpose of this database effort is to make existing high quality data available in detailed form for the turbulence-modeling and computational fluid dynamics communities. While significant recent data were found on the subject of compressible turbulent mixing, the available boundary-layer/pressure-gradient experiments are all older ones of which no acceptable data were found at hypersonic Mach numbers.

Description: Hydrogen-air mixing and combustion experiments at air velocities approaching 17,000 ft/s have been conducted in the expansion tube in CASL's HYPULSE facility. The data were obtained in the facility's sub-millisecond steady flow test period. The results proved to be repeatable and indicated increasing extent of reaction as the fuel-air ratio was increased from stoichiometric to three times stoichiometric. Comparisons were made with the results from a one-dimensional, finite-rate chemistry code and from one-dimensional cycle codes which assume equilibrium chemistry. These comparisons indicate that the magnitudes and trends of the measured longitudinal pressure distributions are predicted by the one-dimensional codes provided that the relevant physical phenomena are accounted for in the computations.

Description: Comparisons between scramjet combustor data and a three-dimensional full Navier-Stokes calculation have been made to verify and substantiate computational fluid dynamics (CFD) codes and application procedures. High Mach number scramjet combustor development will rely heavily on CFD applications to provide wind tunnel-equivalent data of quality sufficient to design, build and fly hypersonic aircraft. Therefore. detailed comparisons between CFD results and test data are imperative. An experimental case is presented, for which combustor wall static pressures were measured and flow-fieid interferograms were obtained. A computer model was done of the experiment, and counterpart parameters are compared with experiment. The experiment involved a subscale combustor designed and fabricated for the National Aero-Space Plane Program, and tested in the Calspan Corporation 96" hypersonic shock tunnel. The combustor inlet ramp was inclined at a 20 angle to the shock tunnel nozzle axis, and resulting combustor entrance flow conditions simulated freestream M=10. The combustor body and cowl walls were instrumented with static pressure transducers, and the combustor lateral walls contained windows through which flowfield holographic interferograms were obtained. The CFD calculation involved a three-dimensional time-averaged full Navier-Stokes code applied to the axial flow segment containing fuel injection and combustion. The full Navier-Stokes approach allowed for mixed supersonic and subsonic flow, downstream-upstream communication in subsonic flow regions, and effects of adverse pressure gradients. The code included hydrogen-air chemistry in the combustor segment which begins near fuel injection and continues through combustor exhaust. Combustor ramp and inlet segments on the combustor lateral centerline were modelled as two dimensional. Comparisons to be shown include calculated versus measured wall static pressures as functions of axial flow coordinate, and calculated path-averaged density contours versus an holographic Interferogram.

Description: A vapor block resistant liquid artery structure for heat pipes. A solid tube artery with openings is encased in the sintered material of a heat pipe wick. The openings are limited to that side of the artery which is most remote from the heat source. The liquid in the artery can thus exit the artery through the openings and wet the sintered sheath, but vapor generated at the heat source is unlikely to move around the solid wall of the artery and reverse its direction in order to penetrate the artery through the openings. An alternate embodiment uses finer pore size wick material to resist vapor entry.

Description: A thermal microstructure measurement system (TMMS) operates autonomously h its own internal power supply and telemeters data to a platform. A thermal array is mounted on a cross-braced frame designed to orient itself normal to existing currents with fixed sensor positioning bars protruding from the cross bars. A plurality of matched thermistors, conductivity probes and inclinometers are mounted on the frame. A compass and pressure transducer are contained in an electronics package suspended below the array. The array is deployed on a taut mooring below a subsurface float. Data are digitized, transmitted via cable to a surface buoy and then telemetered to the platform where the data is processed via a computer, recorded and/or displayed. The platform computer also sends commands to the array via telemetry.

Description: The invention concerns the mounting of propeller blades to a ring-shaped rotor. The blades are of the variable pitch type, and the shank of each blade extends through a respective hole in the rotor. Each hole contains an annular shelf which is fastened to the wall of the hole and surrounds each shank. Each shank bears a pair of bearing races which sandwich the annular shelf in order to connect the blade to the rotor. Bearing rollers are positioned between the annular shelf and the bearing races.

Description: A fluid channeling system includes a fluid ejector, a heat exchanger, and a fluid pump disposed in series flow communication The ejector includes a primary inlet for receiving a primary fluid, and a secondary inlet for receiving a secondary fluid which is mixed with the primary fluid and discharged therefrom as ejector discharge. Heat is removed from the ejector discharge in the heat exchanger, and the heat exchanger discharge is compressed in the fluid pump and channeled to the ejector secondary inlet as the secondary fluid In an exemplary embodiment, the temperature of the primary fluid is greater than the maximum operating temperature of a fluid motor powering the fluid pump using a portion of the ejector discharge, with the secondary fluid being mixed with the primary fluid so that the ejector discharge temperature is equal to about the maximum operating temperature of the fluid motor.

Description: This paper describes a parallel implementation of the direct simulation Monte Carlo (DSMC) method. Runtime library support is used for scheduling and execution of communication between nodes, and domain decomposition is performed dynamically to maintain a good load balance. Performance tests are conducted using the code to evaluate various remapping and remapping-interval policies, and it is shown that a one-dimensional chain-partitioning method works best for the problems considered. The parallel code is then used to simulate the Mach 20 nitrogen flow over a finite-thickness flat plate. It is shown that the parallel algorithm produces results which compare well with experimental data. Moreover, it yields significantly faster execution times than the scalar code, as well as very good load-balance characteristics.

Description: Reynolds number and cowl position effects on the internal shock structure and the resulting performance of a generic three-dimensional sidewall compression scramjet inlet with a leading edge sweep of 45 degrees at Mach 10 have been examined both computationally and experimentally. Prior to the experiment, a three-dimensional Navier-Stokes code was adapted to perform preliminary parametric studies leading to the design of the present configuration. Following this design phase, the code was then utilized as an analysis tool to provide a better understanding of the flow field and the experimental static pressure data for the final experimental configuration. The wind tunnel model possessed 240 static pressure orifices distributed on the forebody plane, sidewalls, and cowl and was tested in the NASA Langley 31 Inch Mach 10 Tunnel.

Description: An apparatus for the levitation of a liquid drop by a fluid flow comprising a profile generator, a fluid flow supply means operatively connected to the profile generator. The profile generator includes an elongate cylindrical shell in which is contained a profiling means for configuring the velocity profile of the fluid flow exiting the profile generator.

Description: A clearance control apparatus has a positioning mechanism for controlling clearance between rotor blade tips and a shroud segment of a gas turbine engine casing. The positioning mechanism is supported by the casing, coupled to the shroud segment, and actuatable for moving the shroud segment toward and away from the rotor blade tips to reach a position at which a desired clearance is established. The mechanism includes a shaft and a bellcrank actuating assembly. The shaft is mounted through a passage defined by a mounting structure on the casing for radial movement along a longitudinal axis of the shaft toward and away from the rotor axis and is coupled at its inner end to the shroud segment. The bellcrank actuating assembly is coupled to an outer end of the shaft and mounted to the casing adjacent the shaft for pivotal movement about an axis spaced from and extending transverse to the longitudinal axis of the shaft. The pivotal movement of the bellcrank actuating assembly produces radial movement of the shaft and shroud segment therewith toward and away from the rotor blade tips.

Description: Four turbulence models are described and evaluated for transonic flows using the upwind code CFL3D and the central-difference code TLNS3D. In particular, the effects of recent modifications to the half-equation model of Johnson-King are explored in detail, and different versions of the model are compared. This model can obtain good results for both two-dimensional (2D) and three-dimensional (3D) separated flows. The one-equation models of Baldwin-Barth and Spalart-Allmaras perform well for separated airfoil flows, but can predict the shock too far forward at the outboard stations of a separated wing. The equilibrium model of Baldwin-Lomax predicts the shock location too far aft for both 2D and 3D separated flows, as expected. In general, all models perform well for attached or mildly separated flows.

Description: A capability to calculate surface heating rates has been incorporated in an approximate three-dimensional inviscid technique. Surface streamlines are calculated from the inviscid solution, and the axisymmetric analog is then used along with a set of approximate convective-heating equations to compute the surface heat transfer. The method is applied to blunted axisymmetric and three-dimensional ellipsoidal cones at angle of attack for the laminar flow of a perfect gas. The method is also applicable to turbulent and equilibrium-air conditions. The present technique predicts surface heating rates that compare favorably with experimental (ground-test and flight) data and numerical solutions of the Navier-Stokes (NS) and viscous shock-layer (VSL) equations. The new technique represents a significant improvement over current engineering aerothermal methods with only a modest increase in computational effort.

Description: For aerodynamic and hydrodynamic vehicles, it is highly desirable to reduce drag and noise levels. A reduction in drag leads to fuel savings. In particular for submersible vehicles, a decrease in noise levels inhibits detection. A suggested means to obtain these reduction goals is by delaying the transition from laminar to turbulent flow in external boundary layers. For hydrodynamic applications, a passive device which shows promise for transition delays is the compliant coating. In previous studies with a simple mechanical model representing the compliant wall, coatings were found that provided transition delays as predicted from the semi-empirical e(sup n) method. Those studies were concerned with the linear stage of transition where the instability of concern is referred to as the primary instability. For the flat-plate boundary layer, the Tollmien-Schlichting (TS) wave is the primary instability. In one of those studies, it was shown that three-dimensional (3-D) primary instabilities, or oblique waves, could dominate transition over the coatings considered. From the primary instability, the stretching and tilting of vorticity in the shear flow leads to a secondary instability mechanism. This has been theoretical described by Herbert based on Floquet theory. In the present study, Herbert's theory is used to predict the development of secondary instabilities over isotropic and non-isotropic compliant walls. Since oblique waves may be dominant over compliant walls, a secondary theory extention is made to allow for these 3-D primary instabilities. The effect of variations in primary amplitude, spanwise wavenumber, and Reynolds number on the secondary instabilities are examined. As in the rigid wall case, over compliant walls the subharmonic mode of secondary instability dominates for low-amplitude primary disturbances. Both isotropic and non-isotropic compliant walls lead to reduced secondary growth rates compared to the rigid wall results. For high frequencies, the non-isotropic wall suppresses the amplification of the secondary instabilities, while instabilities over the isotropic wall may grow with an explosive rate similar to the rigid wall results. For the more important lower frequencies, both isotropic and non-isotropic compliant walls suppress the amplification of secondary instabilities compared to the rigid wall results. The twofold major discovery and demonstration of the present investigation are: (1) the use of passive devices, such as compliant walls, can lead to significant reductions in the secondary instability growth rates and amplification; (2) suppressing the primary growth rates and subsequent amplification enable delays in the growth of the explosive secondary instability mechanism.

Description: The effect of far-field boundary conditions on the evolution of a finite-amplitude two-dimensional wave in the Blasius boundary layer is assessed. With the use of the parabolized stability equations (PSE) theory for the numerical computations, either asymptotic, Dirichlet, Neumann or mixed boundary conditions are imposed at various distances from the wall. The results indicate that asymptotic and mixed boundary conditions yield the most accurate mean-flow distortion and unsteady instability modes in comparison with the results obtained with either Dirichlet or Neumann conditions.

Description: The Interface Configuration Experiment (ICE) was carried out on USML-1 to investigate liquid-gas interfaces in certain rotationally-symmetric containers having prescribed, mathematically derived shapes. These containers have the property that they admit an entire continuum of distinct equilibrium rotationally-symmetric interfaces for a given liquid volume and contact angle. Furthermore, it can be shown that none of these interfaces can be stable. It was found, after the containers were filled in orbit, that an initial equilibrium interface from the symmetric continuum re-oriented, when perturbed, to a stable interface that was not rotationally symmetric, in accordance with the mathematical theory.

Description: Tests of a 1/5 scale model of a proposed 153-foot high-speed submarine have been conducted in the Langley full-scale tunnel at the request of the Bureau of Ships, Department of the Navy. The test program included: (1) force tests to determine the drag, control effectiveness, and static stability characteristics for a number of model configurations, both in pitch and in yaw, (2) pressure measurements to determine the boundary-layer conditions and flow characteristics in the region of the propeller, and (3) an investigation of the effects of propeller operation on the model aerodynamic characteristics. In response to oral requests from the Bureau of Ships representatives t hat the basic data obtained in these tests be made available to them as rapidly as possible, this data report has been prepared to present some of the more pertinent results. All test results given in the present paper are for the propeller-removed condition and were obtained at a Reynolds number of approximately 22,300,000 based on model length.

Description: Transfer functions descriptive of the response of most engine variables were determined from transient data that were obtained from approximate step inputs in fuel flow and in exhaust-nozzle area. The speed responses of both spools to fuel flow and to turbine-inlet temperature appeared as identical first-order lags. Response to exhaust-nozzle area was characterized by a first-order lag response of the outer-spool speed, accompanied by virtually no change in inner-spool speed.

Description: The static lateral- and directional-stability characteristics of a high-speed fighter-type airplane, obtained from wind-tunnel tests of a model, are presented. The model consisted of a thin, unswept wing of aspect ratio 2.3 and taper ratio 0.385, a body, and a horizontal tail mounted in a high position on a vertical tail. Rolling-moment, yawing moment, and cross-wind-force coefficients are presented for a range of sideslip angles of -5 deg. to +5 deg, for Mach numbers of 0.90, 1.45, and 1.90. Data are presented which show the effects on the lateral and directional stability of: (1) component parts of the complete model, (2) modification of the empennage so as to provide different heights of the horizontal tail above the wing plane, (3) angle of attack, and (4) dihedral of the wing.

Description: An improved efficiency arcjet thruster has a constrictor and electrically-conductive nozzle anode defining an arc chamber, and an electrically-conductive rod having a tip spaced upstream from the constrictor and defining a cathode spaced from the anode by a gap generally coextensive with the arc chamber. An electrical potential is applied to the anode and cathode to generate an electrical arc in the arc chamber from the cathode to anode. Catalytically decomposed hydrazine is supplied to the arc chamber with generation of the arc so as to produce thermal heating and expansion thereof through the nozzle. The constrictor can have a electrically insulative portion disposed between the cathode tip and the nozzle anode, and an electrically-conductive anode extension disposed along the insulative portion so as to define an auxiliary gap with the cathode tip substantially smaller than the gap defined between the cathode and nozzle anode for facilitating startup of arc generation. The constrictor can also include an electrically-conductive electrode with a variable electrical potential to vary the shape of the arc generated in the arc chamber. Also, the cathode is mounted for axial movement such that the gap between its tip and the nozzle anode can be varied to facilitate a generally nonerosive generation of the electrical arc at startup and reliable steady state operation. Further, the arc chamber can have a nonparallel subsonic-to-supersonic transition configuration, or alternatively solely a nonparallel supersonic configuration, for improved arc attachment.

Description: A lift producing device is disclosed which is adapted to be connected to a vehicle to provide lift to the vehicle when the vehicle is moved relative to a first fluid medium having a first density and viscosity and being in contact with a second fluid medium adjacent the vehicle. The second fluid medium has a second fluid density which is different from the first fluid density. The lift producing device comprises opposed first and second major surfaces joined at a longitudinally extending leading edge and at a longitudinally extending trailing edge, with at least a portion of the longitudinally extending leading edge being spaced from the longitudinally extending trailing edge by a predetermined mean chord length. When the vehicle is moved relative to the first fluid medium at a velocity within a range of predetermined velocities, with each of the velocities having a direction inclined from a plane extending through the leading edge and the trailing edge within a predetermined angular range, a region of high pressure is generated in the first fluid medium adjacent the first major surface and a region of low pressure is generated in the first fluid medium adjacent the second major surface. The lift producing device has a cross-sectional shape which will generate a pressure distribution around the device when the vehicle is moved relative to the first fluid medium at a velocity within the range of predetermined velocities such that the first fluid medium exhibits attached laminar flow along the device for a portion of the predetermined mean chord length from the leading edge to the trailing edge and will neither form a laminar separation bubble adjacent the second major surface of the device, nor exhibit turbulent separation adjacent the second major surface for substantially all of the predetermined mean chord length from the leading edge to the trailing edge. The portion along which attached laminar flow is maintained is the longest portion which will still fulfill the flow separation requirements. A method for producing the foil is also disclosed.

Description: A poppet is modulated between closed and full open positions by a brushless DC motor operating magnetically through a housing to drive a permanent magnet rotor which carries the poppet. The rotor is supported on several parallel cables which are stationarily fixed at one end and attached to the rotor at the other end, whereby rotation of the rotor twists the cables, causing axial foreshortening and axial translation of rotor and poppet. Axial translation is enhanced by placing a spacer between the cables, intermediate their ends. A permanent magnet ring is disposed around the valve seat directly axially attracting the rotor to a valve closed position.

Description: An optically indicating de-icing solution for surfaces comprising a freezing point depressant liquid and a compound which exhibits a visually observable change as solid phase domains become present in the solution is disclosed. When applied to a surface, particularly the surface(s) of aircraft, the formation of solid phase ice domains in the liquid provides a distinct and visible change in the appearance of the solution. This allows a determination of ice formation upon the aircraft as well as a determination of the effectiveness of the de-icing solution.

Description: A system that incorporates inertial sensor information into optical flow computations to detect obstacles and to provide alternative navigational paths free from obstacles. The system is a maximally passive obstacle detection system that makes selective use of an active sensor. The active detection typically utilizes a laser. Passive sensor suite includes binocular stereo, motion stereo and variable fields-of-view. Optical flow computations involve extraction, derotation and matching of interest points from sequential frames of imagery, for range interpolation of the sensed scene, which in turn provides obstacle information for purposes of safe navigation.

Description: A system and method for non-intrusively obtaining the thrust value of combustion by-products of a jet engine is disclosed herein. The system includes laser elements for inducing absorption for use in determining the axial velocity and density of the jet flow stream and elements for calculating the thrust value therefrom.

Description: Apparatus for decontaminating a contaminated fluid by using photocatalytic particles. The apparatus includes a reactor tank for holding a slurry of the contaminated fluid and the photocatalytic particles ultraviolet light irradiates the surface of the slurry, thereby activating the photocatalytic properties of the particles. Stirring blades for continuously agitate the irradiated fluid surface maintaining the particles in a suspended state within the fluid. A cross flow filter is used for separating the fluid from the semiconductor powder after the decomposition reaction is ended. The cross flow filter is occasionally back flushed to remove any caked semiconductor powder. The semiconductor powder may be recirculated back to the tank for reuse, or may be stored for future use. A series of reactor tanks may be used to gradually decompose a chemical in the fluid. The fluid may be pretreated to remove certain metal ions which interfere with the photocatalytic process. Such pretreatment may be accomplished by dispersing semiconductor particles within the fluid, which particles adsorb ions or photodeposit the metal as the free metal or its insoluble oxide or hydroxide, and then removing the semiconductor particles together with the adsorbed metal ions/oxides/hydroxide/free metal from the fluid.

Description: A method of constructing a nozzle having cooling channels comprises a shell and a liner which are formed into a body of revolution having an axis of revolution. Helical welds are formed to hold the liner and shell to each other with a channel position being defined between each pair of helical welds. Pressurized fluid which may be a gas or a liquid, is introduced between the weld pairs to outwardly bulge the material of at least one of the liner and shell to define the channels.

Description: A numerical study was performed to investigate the shock-wave/boundary-layer interactions on a flat plate with bleed through one or more circular holes that vent into a plenum. The bleed-hole patterns considered for the study include in-line multiple holes and staggered multiple-row holes that are configured to simulate the patterns used in inlet bleed systems of high performance aircraft. The focus of the study was to examine how the bleed through multiple holes affect bleed rate and the pressure and Mach number distributions. Since the bleed performance was found sensitive to the change in bleed conditions, a computational procedure was developed to give a good turnaround computational time for parametric studies involving changes in bleed hole geometry and the structure of shock-wave/boundary-layer flowfield. The procedure includes the grid-generation methodology and the flow simulation with solutions from the Navier-Stokes equations. The computational techniques permit analysis of complex bleed systems and make possible the investigation of a broader range of design variables associated with inlet bleed operation.

Description: A gas turbine engine hot section combustor liner is provided a non-film cooled portion of a heat transfer wall having a hot surface and a plurality of longitudinally extending micro-grooves disposed in the portion of the wall along the hot surface in a direction parallel to the direction of the hot gas flow. The depth of the micro-grooves is very small and on the order of magnitude of a predetermined laminar sublayer of a turbulent boundary layer. The micro-grooves are sized so as to inhibit heat transfer from the hot gas flow to the hot surface of the wall while reducing NOx emissions of the combustor relative to an otherwise similar combustor having a liner wall portion including film cooling apertures. In one embodiment the micro-grooves are about 0.001 inches deep and have a preferred depth range of from about 0.001 inches to 0.005 inches and which are square, rectangular, or triangular in cross-section and the micro-grooves are spaced about one width apart.

Description: Heat gauges are used to measure heat flow in industrial activities. They must periodically be certified by instruments designed to provide a heat flux measurement standard. CSTAR, a NASA CCDS, and REMTECH have developed a portable heat flux checker/calibrator. The Q-CHEC can be carried to the heat gauge for certification, reducing out of service time for the gauge and eliminating the need for a replacement gauge during certification. It can provide an "end-to-end" check of the instrumentation measurement system or be used as a standalone calibrator. Because Q-CHEC offers on-site capability to detect and eliminate measurement errors, measurements do not have to be repeated, and money is saved.