ALBERT

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: 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: 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 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 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: 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: The rate of heat transfer between a fluid stream in turbulent flow and a smooth, solid wall is largely controlled by the relatively high resistance of the laminar sublayer next to the wall. Although this laminar layer ii extremely thin, heat can be transferred through it only by molecular diffusion. Hence the resistance of this layer is very much greater than for a layer the same thickness farther out in the stream where turbulent exchange is the controlling factor. The thickness of the laminar layer is difficult to define precisely, since there is a gradual transition to the turbulent flow outside, but for the usual scale of many engineering applications almost half the temperature difference between the fluid and the wall occurs in a layer of a few thousands of an inch in thickness. When the wall is made of porous material and a coolant gas is forced through the wall into the stream, it has been found that a very small flow rate of the coolant is remarkably effective in keeping the wall at a low temperature. The coolant flow rate required is such as to give an average velocity normal cooling wall of the order of 1 per cent of the main stream velocity. This flow rate is so low that clearly the injected gas must act as an insulator rather than as a normal coolant. Because of its relatively low velocity, the injected gas can have very little influence on heat convection or momentum transfer in the turbulent stream, and its effect must be confined to the laminar sublayer. The possible influence of the coolant flow on the thickness of the laminar layer will be discussed in Section V.

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: 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: A large number of papers have been devoted to the problem of integration of equations of two-dimensional steady nonvertical adiabatic motion of a gas. Most of these papers are based on the application of the hodograph method of S. A. Chaplygin in which the plane of the hodograph of the velocity is taken as the region of variation of the independent variables in the equations of motion; the equations become linear in this plane. The exact integration of these equations is, however, obtained in the form of infinite series containing hypergeometric functions. The obtaining of such solutions and their investigation involves extensive computations. As a result, methods have been developed for the approximate integration of the equations of motion first transformed to a linear form. S. A. Chaplygin first pointed out such an approximate method applicable to flows in which the Mach number does not exceed 0.4.

Description: A study is made herein of the irrotational adiabatic motion of a gas in the transition from subsonic to supersonic velocities. A shape of the de Laval nozzle is given, which transforms a homogeneous plane-parallel flow at large subsonic velocity into a supersonic flow without any shockwaves beyond the transition line from the subsonic to the supersonic regions of flow. The method of solution is based on integration near the transition line of the gas equations of motion in the form investigated by S. A. Christianovich.

Description: By means of characteristics theory, formulas for the numerical treatment of stationary compressible supersonic flows for the two-dimensional and rotationally symmetrical cases have been obtained from their differential equations.

Description: The turbulent flow in a conical diffuser represents the type of turbulent boundary layer with positive longitudinal pressure gradient. In contrast to the boundary layer problem, however, it is not necessary that the pressure distribution along the limits of the boundary layer(along the axis of the diffuser) be given, since this distribution can be obtained from the computation. This circumstance, together with the greater simplicity of the problem as a whole, provides a useful basis for the study of the extension of the results of semiempirical theories to the case of motion with a positive pressure gradient. In the first part of the paper,formulas are derived for the computation of the velocity and.pressure distributions in the turbulent flow along, and at right angles to, the axis of a diffuser of small cone angle. The problem is solved.

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: Strain gages were used to measure blade vibrations possibly causing failure in the 10-stage compressor of the 19XB jet-propulsion engine. The seventh and tenth stages were of great concern as a result of failures experienced by the manufacturer. Strain-gage records were obtained from all stages during acceleration, deceleration, and constant speed runs. Curves are presented herein showing the maximum allowable vibratory stress for a given speed, the change of the damping coefficient with the mounting of a strain gage at the base of the blade, the effect of rotor speed, on blade natural frequency, and the effect of the order of first bending-mode vibration on stress. It was found that for all stages the lower the order of vibration the higher the stress but no destructive vibrations were detected.

Description: An investigation has been conducted on a one-sixth segment of an annular turbojet combustor to determine the effects of modification in air-flow distribution and total-pressure loss on the performance of the segment. The performance features investigated during this series of determinations were the altitude operational limits and the temperature-rise efficiency. Altitude operational limits of the combustor segment, for the 19XB engine using the original combustor-basket design were approximately 38,000 feet at 17,000 rpm and 26,000 feet at 10,000 rpm. The altitude operational limits were approximately 50,000 feet at 17,000 rpm and 38,000 feet at 10,000 rpm for a combustor-basket design in which the air-passage area in the basket was redistributed so as to admit gradually no more than 20 percent of the air along the first half of the basket. In this case the total pressure loss through the combustor segment was not appreciably changed from the total-pressure loss for the original combustor basket design. Altitude operational limits of the combustor segment for the 19XB engine were above 52,000 feet at 17,000 rpm and were approximately 23,000 feet at 10,000 rpm for a combustor-basket design in which the distribution of the air-passage area in the basket was that of the original design but where the total-pressure loss was increased to 19 times the inlet reference kinetic pressure at an inlet-to-outlet density ratio of 2.4. The total-pressure loss for the original design was 14 times the inlet kinetic reference pressure at an inlet-to-outlet density ratio of 2.4.

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.

Description: In one type of aircraft propulsion system, propeller blades are mounted on a ring which surrounds a turbine. An annular space exists between the turbine and the ring. If a propeller blade should break free, the unbalanced centrifugal load tends to deform the ring. The invention reduces the deformation, as by locating spacers between the turbine and the ring.

Description: The goal of the symposium was to examine the state of technology of all areas of magnetic suspension and to review related recent developments in sensors and controls approaches, superconducting magnet technology, and design/implementation practices. The symposium included 17 technical sessions in which 55 papers were presented. The technical session covered the areas of bearings, sensors and controls, microgravity and vibration isolation, superconductivity, manufacturing applications, wind tunnel magnetic suspension systems, magnetically levitated trains (MAGLEV), space applications, and large gap magnetic suspension systems.

Description: A heat-transfer investigation was conducted with air flowing through an electrically heated silicon carbide tube with a rounded entrance, an inside diameter of 3/4 inch, and effective heat-transfer length of 12 inches over a range of Reynolds numbers up to 300,000 and a range of average inside-tube-wall temperatures up to 2500 R. The highest corresponding local outside-tube-wall temperature was 3010 R.

Description: A method and apparatus for sensing wave flow across a surface wherein at least two pressure levels are sensed and combined to provide a representation of waves within the flow. In the preferred embodiment holes bored through the aircraft surface at an interval of one-half the wavelength of the flow being measured introduce pressure perturbations into a cavity so they may acoustically interfere. The interfering waveform is sensed by at least one microphone disposed in the cavity.

Description: The calculation of the phenomena within the boundary layer of bodies immersed in a flow underwent a decisive development on the basis of L. Prandtl's trains of thought, stated more than forth years ago, and by numerous later treatises again and again touching upon them. The requirements of the steadily improving aerodynamics of airplanes have greatly increased with the passing of time and recently research became particularly interested in such phenomena in the boundary layer as are caused by small external disturbances. Experimental results suggest that, for instance, slight fluctuations in the free stream velocities as they occur in wind tunnels or slight wavelike deviations of outer wing contours from the prescribed smooth course as they originate due to construction inaccuracies may exert strong effects on the extent of the laminar boundary layer on the body and thus on the drag. The development of turbulence in the last part of the laminar portion of the boundary layer is, therefore, the main problem, the solution of which explains the behavior of the transition point of the boundary layer. A number of reports in literature deal with this problem,for instance, those of Tollmien, Schlichting, Dryden, and Pretsch. The following discussion of the behavior of the laminar boundary layer for periodically oscillating pressure variation also purports to make a contribution to that subject.

Description: Some aerodynamic relations are derived which exist between two infinitely long airfoils if one is in a straight flow and the other in oblique flow, and both present the same profile in the direction of flow.

Description: At the request of the Junkers Aircraft and Engine Construction Company, Engine Division, Dessau Main Plant, an investigation was made using the interferometer method on the two turbine-blade profiles submitted. The interferometer method enables making visible the differences in density and consequently the boundary layers that develop when a flow is directed on the profile. Recognition of the points on the profile at which separation of flow occurs is thus possible. By means of the interference photographs the extent of the dead-water region may be ascertained. The size of the dead-water region provides evidence as to the quality of the flow and allows a qualitative estimate of the amount of the flow losses. Interference photographs thus provide means of judging the utility of profiles under specific operating conditions and provide suggestions for possible changes of profile contours that might help to improve flow relations. Conclusions may be drawn concerning the influence of the blade-spacing ratio, the inlet-air angle, and the connection between the curvature of the profile contour and the point of separation of the flow from the profile surface.

Description: The cavitation in nozzles on airfoils of various shape and on a sphere are experimentally investigated. The limits of cavitation and the extension of the zone of the bubbles in different stages of cavitation are photographically established. The pressure in the bubble area is constant and very low, jumping to high values at the end of the area. The analogy with the gas compression shock is adduced and discussed. The collapse of the bubbles under compression shock produces very high pressures internally, which must be contributory factors to corrosion. The pressure required for purely mechanical corrosion is also discussed.

Description: This paper includes the following topics: 1) Characteristic differential equations; 2) Treatment of practical examples; 3) First example: Diffuser; and 4) Second Example: Nozzle.

Description: This paper contains a tabulation of functions of the Mach number which are frequently used in high-speed aerodynamics. The tables extend from M = 0 to M = 10.0 in increments of 0.01 and are based on the assumption that air is a perfect gas having a specific heat ratio of 1.400.

Description: The weakly nonlinear evolution of an inviscid marginally unstable wave growing on a boundary layer supporting a streamwise vortex structure is investigated. The nonlinear growth of the wave is found to be controlled by the diffusion layer located at the edge of the critical layer associated with the wave. The evolution equation is found to depend on the upstream history of the wave and the solution of the equation suggests that the wave either restructures the mean state so as to make it stable or develops a singularity at a finite distance downstream of the point of neutral stability.

Description: One or the key objectives of the Applied Research Branch in the Numerical Aerodynamic Simulation (NAS) Systems Division at NASA Allies Research Center is the accelerated introduction of highly parallel machines into a full operational environment. In this report we discuss the performance results obtained from the implementation of some computational fluid dynamics (CFD) applications on the Connection Machine CM-2 and the Intel iPSC/860. We summarize some of the experiences made so far with the parallel testbed machines at the NAS Applied Research Branch. Then we discuss the long term computational requirements for accomplishing some of the grand challenge problems in computational aerosciences. We argue that only massively parallel machines will be able to meet these grand challenge requirements, and we outline the computer science and algorithm research challenges ahead.

Description: A steam reformer is disclosed having an annular steam reforming catalyst bed formed by concentric cylinders and having a catalytic combustor located at the center of the innermost cylinder. Fuel is fed into the interior of the catalytic combustor and air is directed at the top of the combustor, creating a catalytic reaction which provides sufficient heat so as to maintain the catalytic reaction in the steam reforming catalyst bed. Alternatively, air is fed into the interior of the catalytic combustor and a fuel mixture is directed at the top. The catalytic combustor provides enhanced radiant and convective heat transfer to the reformer catalyst bed.

Description: The conference on Turbojet-Engine Thrust-Augmentation Research was organized by the NACA to present in summarized form the results of the latest experimental and analytical investigations conducted at the Lewis Flight Propulsion Laboratory on methods of augmenting the thrust of turbojet engines. The technical discussions are reproduced herewith in the same form in which they were presented. The original presentation in this record are considered as complementary to, rather than substitutes for, the committee's system of complete and formal reports.

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: A two-dimensional unstructured Navier-Stokes code is utilized for computing the flow around multielement airfoil configurations. Comparisons are shown for a landing configuration with an advanced-technology flap. Grid convergence studies are conducted to assess inaccuracies caused by inadequate grid resolution. Although adequate resolution is obtained for determining the pressure distributions, further refinement is needed to sufficiently resolve the velocity profiles at high angles of attack. For the advanced flap configuration, comparisons of pressure distributions and lift are made with experimental data. Here, two flap riggings and two Reynolds numbers are considered. In general, the trends caused by variations in these quantities are well predicted by the computations, although the angle of attack for maximum lift is overpredicted.

Description: The condensation of water vapor in an air consequences: acquisition of heat (liberated heat vaporization; loss of mass on the part of the flowing gas (water vapor is converted to liquid); change in the specific gas constants and of the ratio k of the specific heats (caused by change of gas composition). A discontinuous change of state is therefore connected with the condensation; schlieren photographs of supersonic flows in two-dimensional Laval nozzles show two intersecting oblique shock fronts that in the case of high humidities may merge near the point of intersection into one normal shock front.