ALBERT

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: 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: 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: 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 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: 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: 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: 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: 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 subsonic wind tunnel investigation of pneumatic vortex flow control on a chined forebody using slots was accomplished at a dynamic pressure of 50 psf resulting in a R(n)/ft of 1.3 x 10(exp 6). Data were acquired from angles of attack ranging from -4deg to +34deg at side slips of +0.4deg and +10.4deg. The test article used in this study was the 10% scale Fighter Lift and Control (FLAC) advanced diamond winged, vee-tailed fighter configuration. Three different slot blowing concepts were evaluated; outward, downward, and tangential with ail blowing accomplished asymmetrically. The results of three different mass flows (0.067, 0.13, and 0.26 lbm/s; C(sub mu)'s of less than or equal to 0.006, 0.011. and 0.022 respectively) were analyzed and reported. Test data are presented on the effects of mass flows, slot lengths and positions and blowing concepts on yawing moment and side force generation. Results from this study indicate that the outward and downward blowing slots developed yawing moment and side force increments in the direction opposite of the blowing side while the tangential blowing slots generated yawing moment and side force increments in the direction towards the blowing side. The outward and downward blowing slots typically produced positive pitching moment increments while the tangential blowing slots typically generated negative pitching moment increments. The slot blowing nearest the forebody apex was most effective at generating the largest increments and as the slot was moved aft or increased in length, its effectiveness at generating forces and moments diminished.

Description: The phenomena of rotational relaxation of nitrogen has been examined by numerous investigators over many years. One of the experiments which has been performed examines nonequilibrium flow in low-density free jet expansions. Data have been taken in such flows using a variety of techniques, including time-of-flight methods and electron beam fluorescence spectroscopy. The direct flow properties measured in these different investigations, such as density and translational, rotational and vibrational temperatures generally show reasonable agreement. However, this kind of correlation from experiment to experiment tends to be lost when these data are analyzed to obtain rotational relaxation time or collision number. The goal of such data analyses is to generate a succinct model for rotational relaxation in nitrogen which is essential for the computation of nonequilibrium rarefied flows. The objective of the present work is to process a large body of experimental data in a consistent manner to yield relaxation model parameters of the greatest utility for flow computations.

Description: This study focuses on the suppression of instability growth using an automated active-control technique. The evolution of 2D disturbances that are spatially growing in a flat-plate boundary layer are computed with a spatial DNS code. A controller receives wall sensor information (pressure or shear) as input and provides a signal that controls an actuator response as output. The control law assumes that wave cancellation is valid. The results indicate that a measure of wave cancellation can be obtained for small- and large-amplitude instabilities without feedback; however, feedback is required to optimize the control amplitude and phase for exact wave cancellation.

Description: An apparatus comprising a rotatable mass of structured packing for mass or heat transfer between two contacting fluids of different densities wherein the packing mass is made up of corrugated sheets of involute shape relative to the axis of the packing mass and form a logarithmic spiral curved counter to the direction of rotation.

Description: Study of turbulent flows in rotating reference frames has proven to be one of the more challenging areas of turbulence research. The large number of theoretical, experimental, and computational studies performed over the years have demonstrated that the effect of solid-body rotation on turbulent flows is subtle and remains exceedingly difficult to predict. Because of the complexities associated with non-homogeneous turbulence, it is worthwhile to examine the effect of steady system rotation on the evolution of an initially isotropic turbulent flow. The assumption of statistical homogeneity considerably simplifies analysis and computation; calculation of homogeneous turbulence is further motivated since it possesses the essential physics found in more complex rotating flows. The principal objectives of the present study have therefore been to increase our fundamental understanding of turbulent flows in rotating reference frames through an examination of the asymptotic state of homogeneous rotating turbulence; particularly as to the existence of an asymptotic state which is self similar. Knowledge of an asymptotic similarity state permits prediction of the ultimate statistical evolution of the flow without requiring detailed knowledge of the complex, and not well understood, non-linear transfer processes. Aside from examination of possible similarity states in rotating turbulence, of further interest in this study has been an examination of the degree to which solid-body rotation induces a two-dimensional state in an initially isotropic flow.

Description: The results from a concentrated computational effort are presented with the primary objective being directed at evaluating the vortical-flow-prediction capability of an unstructured-grid Euler solver. Both viscous and inviscid solutions, obtained from an established structured-grid method, along with an experimental wind-tunnel data are used as bench-mark measures to assess the validity of the unstructured-grid Euler results. Viscous effects on vortical flows are first identified by comparing the viscous and inviscid solutions obtained form the structured-grid method. Computational data analysis are then presented which reveal excellent correlations between the inviscid structured and unstructured-grid results in terms of off-surface flow structures, surface pressure distribution and the predicted longitudinal aerodynamic characteristics. The sensitivity of the unstructured-grid inviscid solutions to grid refinement is also discussed along with an analysis of the convergence and performance characteristics for each method.

Description: The urgent need for dramatic reductions in aircraft design cycle time is focusing scrutiny upon all aspects of computational fluid dynamics (CFD). These reductions will most likely come not from increased reliance upon user-interactive (and therefore time-expensive) methods, but instead from methods that can be fully automated and incorporated into 'black box' solutions. In comparison with tetrahedral methods, three-dimensional Cartesian grid approaches are in relative infancy, but initial experiences with automated Cartesian techniques are quite promising. Our research is targeted at furthering the development of Cartesian methods so that they can become key elements of a completely automatic grid generation/flow solution procedure applicable to the Euler analysis of complex aircraft geometries.

Description: The proposed paper presents flow visualization performed during experiments conducted on a full-scale F/A-18 aircraft in the 80- by 120-Foot Wind-Tunnel at NASA Ames Research Center. The purpose of the flow-visualization experiments was to document the forebody and leading edge extension (LEX) vortex interaction along with the wing flow patterns at high angles of attack and low speed high Reynolds number conditions. This investigation used surface pressures in addition to both surface and off-surface flow visualization techniques to examine the flow field on the forebody, canopy, LEXS, and wings. The various techniques used to visualize the flow field were fluorescent tufts, flow cones treated with reflective material, smoke in combination with a laser light sheet, and a video imaging system for three-dimension vortex tracking. The flow visualization experiments were conducted over an angle of attack range from 20 deg to 45 deg and over a sideslip range from -10 deg to 10 deg. The various visualization techniques as well as the pressure distributions were used to understand the flow field structure. The results show regions of attached and separated flow on the forebody, canopy, and wings as well as the vortical flow over the leading-edge extensions. This paper will also present flow visualization comparisons with the F-18 HARV flight vehicle and small-scale oil flows on the F-18.

Description: A formalism will be presented that allows the transformation of two-equation eddy viscosity turbulence models into one-equation models. The transformation is based on an assumption that is widely accepted over a large range of boundary layer flows and that has been shown to actually improve predictions when incorporated into two-equation models of turbulence. Based on that assumption, a new one-equation turbulence model will be derived. The new model will be tested in great detail against a previously introduced one-equation model and against its parent two-equation model.

Description: This viewgraph presentation gives an overview of the visual environments for computational fluid dynamics (CFD) research. It includes details on critical needs from the future computer environment, features needed to attain this environment, prospects for changes in and the impact of the visualization revolution on the human-computer interface, human processing capabilities, limits of personal environment and the extension of that environment with computers. Information is given on the need for more 'visual' thinking (including instances of visual thinking), an evaluation of the alternate approaches for and levels of interactive computer graphics, a visual analysis of computational fluid dynamics, and an analysis of visualization software.

Description: The Web provides new Methods for accessing Information world-wide, but the current text-and-pictures approach neither utilizes all the Web's possibilities not provides for its limitations. While the inclusion of pictures and animations in a paper communicates more effectively than text alone, It Is essentially an extension of the concept of "publication." Also, as use of the Web increases putting images and animations online will quickly load even the "Information Superhighway." We need to find forms of communication that take advantage of the special nature of the Web. This paper presents one approach: the use of the Internet and the Mosaic interface for data sharing and collaborative analysis. We will describe (and In the presentation, demonstrate) our approach: using FAST (Flow Analysis Software Toolkit), a scientific visualization package, as a data viewer and interactive tool called from MOSAIC. Our intent is to stimulate the development of other tools that utilize the unique nature of electronic communication.

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 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: 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: 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: 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 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: 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: 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: 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: Compressibility effects within decaying isotropic turbulence and homogeneous turbulent shear flow have been studied using direct numerical simulation. The objective of this work is to increase our understanding of compressible turbulence and to aid the development of turbulence models for compressible flows. The numerical simulations of compressible isotropic turbulence show that compressibility effects are highly dependent on the initial conditions. The shear flow simulations, on the other hand, show that measures of compressibility evolve to become independent of their initial values and are parameterized by the root mean square Mach number. The growth rate of the turbulence in compressible homogeneous shear flow is reduced compared to that in the incompressible case. The reduced growth rate is the result of an increase in the dissipation rate and energy transfer to internal energy by the pressure-dilatation correlation. Examination of the structure of compressible homogeneous shear flow reveals the presence of eddy shocklets, which are important for the increased dissipation rate of compressible turbulence.

Description: A simple set of equations is derived for predicting the temperature distribution and radiative energy flow in a semitransparent layer consisting of an arbitrary number of laminated sublayers that absorb, emit, and scatter radiation. Each sublayer can have a different refractive index and optical thickness. The plane composite region is heated on each exterior side by a different amount of incident radiation. The results are for the limiting case where heat conduction within the layers is very small relative to radiative transfer, and is neglected. The interfaces are assumed diffuse, and all interface reflections are included in the analysis. The thermal behavior is readily calculated from the analytical expressions that are obtained. By using many sublayers, the analytical expressions provide the temperature distribution and heat flow for a diffusing medium with a continuously varying refractive index, including internal reflection effects caused by refractive index gradients. Temperature and heat flux results are given to show the effect of variations in refractive index and optical thickness through the multilayer laminate.

Description: An interrupted monogroove slot in a heat pipe facilitates priming of the heat pipe under zero gravity conditions by preventing the monogroove slot from completely priming before the liquid channel is primed.

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: 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: 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: 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: 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: 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 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: A new measurement technique, Doppler global velocimetry, has been described along with results from an experimental investigation of a vortical flow field. The results of this investigation, when compared to conventional laser velocimeter measurements, indicate that this technique is capable of describing the global velocity field within a measurement plane in real time. Further investigations are being conducted to correct the deficiencies described and to refine the measurement technique for flight applications.

Description: PHOENICS was used to simulate turbulent flow of supercritical hydrogen through asymmetrically heated rectangular ducts with one side rib-roughened, in order to provide design data for active cooling systems of hypersonic aircraft where the hydrogen fuel is used as coolant. Variable property and conjugate heat transfer effects were considered. The computations using PHOENICS required provision of property sub-routines for supercritical hydrogen, a low-Reynolds number form of the k-epsilon turbulence model, a rough-wall wall function and a procedure for calculating proper diffusion coefficients at the fluid-solid interface. The effect of each design variable such as width and height of duct and wall thickness, on pressure drop and maximum wall temperature was investigated. PHOENICS version 1.4 was used on an IBM 3090 computer.

Description: A life support back pack 14 for use during extravehicular activity in space incorporates a cooling apparatus 20 comprised of five panels 22 each of which include in layered fashion a LCG coolant heat exchange coil 32, a heat distribution plate 42, and a heat dissipation module 50A or 50B having an outer radiator surface 52. Each module 50A houses a first phase change material 55A, for example hexadecane paraffin, and each module 50B houses a second phase change material 55B, for example tetradecane paraffin, which has a phase change temperature which is less than the phase change temperature of the first phase change material 55A. The cooling apparatus 20 is equipped with a coolant heat exchange circuit provided with mode selection valves 84 and 86 which are operated by a controller 88 to selectively direct the LCG coolant to be cooled through the cooling apparatus in one of three operating modes.

Description: A moment method for computing 3-D radiative transport is applied to axisymmetric flows in thermochemical nonequilibrium. Such flows are representative of proposed aerobrake missions. The method uses the P-1 approximation to reduce the governing system of integro-di erential equations to a coupled set of partial di erential equations. A numerical solution method for these equations given actual variations of the radiation properties in thermochemical nonequilibrium blunt body flows is developed. Initial results from the method are shown and compared to tangent slab calculations. The agreement between the transport methods is found to be about 10 percent in the stagnation region, with the difference increasing along the flank of the vehicle.

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: 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: 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: 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: 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 method for predicting radiative heating and coupling effects in nonequilibrium flow-fields has been developed. The method resolves atomic lines with a minimum number of spectral points, and treats molecular radiation using the smeared band approximation. To further minimize computational time, the calculation is performed on an optimized spectrum, which is computed for each flow condition to enhance spectral resolution. Additional time savings are obtained by performing the radiation calculation on a subgrid optimally selected for accuracy. Representative results from the new method are compared to previous work to demonstrate that the speedup does not cause a loss of accuracy and is sufficient to make coupled solutions practical. The method is found to be a useful tool for studies of nonequilibrium flows.