ALBERT

Description: The low-g fluids management group with the Center for Space Construction is engaged in active research on the following topics: gauging; venting; controlling contamination; sloshing; transfer; acquisition; and two-phase flow. Our basic understanding of each of these topics at present is inadequate to design space structures optimally. A brief report is presented on each topic showing the present status, recent accomplishings by our group and our plans for future research. Reports are presented in graphic and outline form.

Description: Microgravity offers a unique environment for studying polymer diffusion and polymer polymerization reactions. The absence of convection currents, which are the major mode of mixing at the molecular level on Earth, are eliminated or reduced in the microgravity environment. More importantly, the prediction of unique copolymer composition development in microgravity allows controlled formation of new compositions of matter. The absence of mixing at the molecular level should produce unique short block copolymers available for the first time for comonomer compositions which normally lead to random or long block copolymer under good mixing. The investigation of fundamental polymer diffusion and polymer polymerization processes in microgravity is proposed. This effort will involve fundamental studies of monomer and polymer diffusion; their effects on initiation, propagation, and especially termination kinetics rate constant; and the accurate evaluation of copolymerization reactivity ratios in microgravity. The experimental design is presented for these studies along with an evaluation technique for in situ monitoring of polymer diffusion and polymerization kinetics.

Description: The analysis of light scattered from an ensemble of particles has long been a preferred method for characterizing their physical properties. Instrumentation to perform the measurements which forms the basis for such analysis is available in many forms based upon a variety of different experimental techniques. A system is presented which is singularly applicable for making many types of measurements in a microgravity environment. The commercial version of this device, the DAWN-F, has been used in many labs throughout the world to perform analyses of particular importance for both research and production. Light scattering theory is reviewed and the structure and function of the system is described.

Description: Avalanche photodiodes (APD) are excellent small area, solid state detectors for photon counting. Performance possibilities include: photon detection efficiency in excess of 50 percent; wavelength response from 400 to 1000 nm; count rate to 10 (exp 7) counts per sec; afterpulsing at negligible levels; timing resolution better than 1 ns. Unfortunately, these performance levels are not simultaneously available in a single detector amplifier configuration. By considering theoretical performance predictions and previous and new measurements of APD performance, the anticipated performance of a range of proposed APD-based photon counting modules is derived.

Description: Traditional optical systems for photon correlation spectroscopy and laser anemometry have relied upon physically large and fairly expensive lasers, bulk-optics such as lenses of a few inches diameter, large mechanical mounts and carefully selected, fragile and bulky photon counting photomultiplier detectors. In some cases, experimental fluid dynamics at a desired position in a flow, perhaps deep inside complex machinery, is physically impossible or very difficult. Similar problems exist with photon correlation spectroscopy, e.g., remote and heterodyne experiments. Various optical and electro optical components were investigated and characterized with the aim of replacing existing photon correlation laser spectroscopy and anemometry techniques in miniaturized form, and with significant cost reduction. Very recently, a range of miniature, modular light scattering systems were constructed from little solid state optical and electro optical components, and experimentally verified measurement performance comparable to standard lab photon correlation spectroscopy and laser anemometry equipment.

Description: Experimental and theoretical studies have been conducted to determine critical parameters at the onset of nonlinear counterflow in He II below the lambda point of He-4. Critical temperature differences have been measured in porous media for zero net mass flow and for Darcy permeabilities in the order of magnitude range from 10 to the -10th to 10 to the -8th sq cm. The normalized critical temperature gradients, which covered the liquid temperature range of 1.5 K to the lambda temperature, are found to vary with T proportional to the ratio of the superfluid density to the normal fluid density. This liquid temperature dependence appears to be consistent with duct data which are limited at low temperature by a Reynolds number criterion.

Description: Images from an airborne, scanning radiometer operating at a frequency of 98 GHz have been analyzed. The millimeter-wave images were obtained in 1985-1986 using the JPL millimeter-wave imaging sensor. The goal of this study was to enhance the information content of these images and make their interpretation easier. A visual-interpretative approach was used for information extraction from the images. This included application of nonlinear transform techniques for noise reduction and for color, contrast, and edge enhancement. Results of using the techniques on selected millimeter-wave images are discussed.

Description: An analysis is made based upon the concept that the velocity fluctuations, and therefore, the Reynolds stresses, driven by the instability of the original flow grow until a new stable state is approached. The Reynolds stresses incorporated into the Orr-Sommerfeld equation are coupled with the main flow such that all the imaginary parts of the complex eigenvalues vanish, i.e., the original instability is eliminated. Using this stabilization principle, it is possible to find the Reynolds stresses as well as the mean velocity for plane Poiseuille flow with the Reynolds number slightly higher than the critical.

Description: The behavior of the reverse flow ceiling jet against the ventilation flow from 0.58 to 0.87 m/s was investigated in a 1/3 scale model of a wide body aircraft interior. For all tests, strong reverse-flow ceiling jets of hot gases were detected well upstream of the fire. Both thicknesses of the reverse-flow ceiling jet and the smoke layer increased with the fire-crossflow parameter. The thickness of the smoke layer where the smoke flows along the main flow below the reverse-flow ceiling jet was almost twice that of the reverse-flow ceiling jet. Detailed spatial and time-varying temperatures of the gas in the test section were measured, and velocity profiles were also measured using a temperature compensated hot film.

Description: Certain theoretical studies of boundary-layer transition are described, based on high Reynolds numbers and with attention drawn to the various nonlinear interactions and scales present. The article concentrates in particular on theories for which the mean-flow profile is completely altered from its original state. Two- and three-dimensional flow theory and conjectures on turbulent-boundary-layer structures are included. Specific recent findings noted, and in qualitative agreement with experiments, are: nonlinear finite-time break-ups in unsteady interactive boundary layers; strong vortex/wave interactions; and prediction of turbulent boundary-layer displacement- and stress sublayer-thicknesses.

Description: This paper presents the application of a class of multi-grid methods to the solution of the Navier-Stokes equations for two-dimensional laminar flow problems. The methods consists of combining the full approximation scheme-full multi-grid technique (FAS-FMG) with point-, line- or plane-relaxation routines for solving the Navier-Stokes equations in primitive variables. The performance of the multi-grid methods is compared to those of several single-grid methods. The results show that much faster convergence can be procured through the use of the multi-grid approach than through the various suggestions for improving single-grid methods. The importance of the choice of relaxation scheme for the multi-grid method is illustrated.

Description: It is shown how a commercial time interval counter can be used to measure the relative stability of two signals that are offset in frequency and mixed down to a beat note of about 1 Hz. To avoid the dead-time problem, the counter is set up to read the time interval between each beat note upcrossing and the next pulse of a 10 Hz reference pulse train. The actual upcrossing times are recovered by a simple algorithm whose outputs can be used for computing residuals and Allan variance. A noise floor-test yielded a delta f-f Allan deviation of 1.3 times 10 to the minus 9th power/tau relative to the beat frequency.

Description: Rapid distortion theory is applied to study distortion of homogeneous turbulence subject to two different axisymmetric strain modes: the axisymmetric contraction (AC, nozzle-type flow), and the axisymmetric expansion (AE, diffuser-type flow). The paper explores the differences in effects of the two axisymmetric strain modes on the anisotropy of correlations and structures of turbulence; examines the effect of dilatation on the distortion of turbulence; and provides a theoretical background for turbulence model development. It is found that velocity and vorticity fluctuations are enhanced more efficiently by contraction than by expansion; contraction produces much higher anisotropy in velocity and vorticity than expansion; root-mean-square pressure is slightly reduced during contraction, whereas it increases rapidly during expansion; and vortical structures of rodlike shape develop in a contraction flow, while disklike structures develop in an expansion flow. A simple model that reflects the dependence of turbulence evolution on structural parameters such as the Reynolds-stress anisotropy and total strain is proposed, and is shown to outperform all other models for all cases examined, regardless of the mean strain rate.

Description: A correction factor to the number density measured by the Forward Scattering Spectrometer Probe (FSSP) which compensates for dead time and coincidence errors was determined by calculating the probabilities of, and the average number of particles in, the six possible types of dead time and coincidence events. These probabilities and averages were calculated by means of a probabilistic model based on Poisson statistics. A Monte Carlo computer simulation of the FSSP operation was also carried out and the number density correction factor was compared with the Monte Carlo data. For an actual number density of 2000/cu cm, it was found that the measured number density was of the order of 300/cu cm.

Description: Measurement of times of flight of sound waves can be used to determine temperatures in a gas. This paper describes a system, based on this principle, that is capable of giving the temperature profile in a nonisothermal gas volume, for example, prevalent in a large furnace. The apparatus is simple, rugged, accurate, and capable of being automated for process control applications. It is basically an acoustic waveguide where the outside temperature profile is transferred to a chosen gas contained inside the guide.

Description: The development of a unique noncontact temperature measurement device utilizing rotating analyzer ellipsometry is described. The technique circumvents the necessity of spectral emissivity estimation by direct measurement concomitant with radiance brightness. Simultaneous determinations of dielectric constants and refractive indices allow changes in the physical and chemical state of a heated surface to be monitored. The results of optical property measurements at 633 nm as functions of temperature between 1000 and 2500 K for eight transition metals including Hf, Ir, Mo, Nb, Pd, Pt, Ta, and V are presented together with preliminary results of oxidation studies on iridium.

Description: Numerical studies of turbulent flow in an axisymmetric 45-deg-expansion combustor and bifurcated diffuser are presented. The Navier-Stokes equations incorporating a k-epsilon model were solved in a nonorthogonal curvilinear coordinate system. A zonal-grid method, where the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. Boundary overlap and interpolating techniques were used, and an adjustment of the flow variables was required to assure conservation of mass flux. Three finite-differencing methods (hybrid, quadratic upwind, and skew upwind) were used to represent the convection terms. Results were compared with existing experimental data. In general, good agreement between predicted and measured values was obtained.

Description: A design concept developed for a polarimeter on the vector magnetograph of the SAMEX satellite that would be very sensitive to solar vector magnetic fields is described. A description of the Poincare sphere is presented, along with the instrument scientific requirements, to provide an understanding of how the polarimeter design has been selected. It is shown that the design goal of a polarimetric sensitivity of 0.0001 can be achieved in the linear measurements using a hybrid analyzer. It is also noted that the systematic errors that produce linear crosstalk into the circular measurement will require the use of the redundant polarization measurements for in-flight calibrations and corrections of the data.

Description: The present analysis method for hot-wire data in supersonic turbulence takes sound field effects into account and yields greater accuracy in its treatment of flow variable fluctuations than existing methods despite requiring only a moderately accurate estimate of static pressure fluctuations. The method demonstrates the way in which neglecting pressure fluctuations will affect hot-wire data analysis, as well as indicating the probable direction the errors will take.

Description: Displays of multifrequency passive microwave data from the Special Sensor Microwave/Imager (SSM/I) flying on the Defense Meteorological Satellite Program (DMSP) spacecraft are presented. Observed brightness temperatures at 85.5 GHz (vertical and horizontal polarizations) and 37 GHz (vertical polarization) are respectively used to 'drive' the red, green, and blue 'guns' of a color monitor. The resultant false-color images can be used to distinguish land from water, highlight precipitation processes and structure over both land and water, and detail variations in other surfaces such as deserts, snow cover, and sea ice. The observations at 85.5 GHz also add a previously unavailable frequency to the problem of rainfall estimation from space. Examples of mesoscale squall lines, tropical and extra-tropical storms, and larger-scale land and atmospheric features as 'viewed' by the SSM/I are shown.

Description: A postdispersion system for astronomical observations with Fourier transform spectrometers in the thermal infrared has been developed which improves the sensitivity of radiation noise limited observations by reducing the spectral range incident on the detector. Special attention is given to the first-generation blocked impurity band detector. Planetary, solar, and stellar observations are reported.

Description: The measurement by a three-dimensional laser Doppler velocimeter of a turbulent flow has been numerically simulated. Errors associated with the probe volume geometry and the coincidence time window concept are revealed. One type of error occurs for high system data rates when multiple particles lead to system realizations. Another error occurs associated with a geometric bias discovered in the present study. This three-dimensional ldv geometric bias exists even for single-particle realizations and regardless of the system data rate. A technique for the elimination of the geometric bias is presented.

Description: A computer-controlled tunable diode laser spectrometer and spectral analysis software are described. The three-channel system records simultaneously the transmission of a subject gas, a temperature-stabilized etalon, and a calibration gas. The software routines are applied to diode laser spectra of HNO3 and NO2 to illustrate the procedures adopted for conversion of raw spectral data to useful transmission and harmonic spectra. Extraction of line positions, absorption intensities, collisional broadening coefficients, and gas concentrations from recorded spectra is also described.

Description: Optical interconnects are being considered for control signal distribution in phased array antennas. A packaged hybrid GaAs optical controller with a 1:16 demultiplexed output that is suitable for this application is described. The controller, which was fabricated using enhancement/depletion mode MESFET technology, operates at demultiplexer-limited input data rates up to 305 Mb/s and requires less than 200 microW optical input power.

Description: A technique has been developed to support the study of the effects of cosmic rays on integrated circuits. The system is designed to determine the particle distribution across the surface of an integrated circuit accurately while the circuit is bombarded by a particle beam. The system uses photomultiplier tubes, an octal discriminator, a computer-controlled NIM quad counter, and an IBM PC. It provides real-time operator feedback for fast beam tuning and monitors momentary fluctuations in the particle beam. The hardware, software, and system performance are described.

Description: A three-channel transportable radiometer that operates at 20.6, 31.65, and 90.0 GHz has been developed. The design of the radiometer is described and the results of various experimental applications of the radiometer are presented. The experiments include a comparison of brightness temperature measurements at various frequencies with the water vapor absorption models of Waters (1976) and Lieve (1989). Measurements of absorption ratios for cloud attenuation are given and the radiometer measurements are compared with oxygen absorption calculated using the Rasenkranz (1988) model.

Description: This paper describes a program of airborne radiometric imaging at 90 GHz and 140 GHz. Using high sensitivity (below 1 K) and high angular resolution (0.5-1.0 degree), high quality images have been made. The following measurements are discussed: cloud and fog penetration at 90 GHz, discrimination between agricultural and urban areas, discrimination between different vegetation types, detection of vehicles on roads, detection and classification of airports and airplanes, ship detection and quantitative oil spill sensing. The application of information enhancement techniques with automatic and real time application aspects is also described, and results of applied techniques for contrast and contour enhancement are shown.

Description: A numerical procedure in which the Navier-Stokes equations are discretized using tightly coupled discretizations of pressure derivatives and continuity equations is used here to extend the range of known terminal velocities of gaseous bubbles in liquids well beyond that in previous investigations. Computations performed for Reynolds numbers up to 2000 and Marangoni numbers up to 1000 show only a modest variation of the scaled bubble velocity between 0.16 and 0.5. The bubble velocity is influenced more by the Marangoni number than by the Reynolds number.

Description: A numerical analysis is performed on thermocapillary buoyancy convection induced by phase change in a liquid droplet. A finite-difference code is developed using an alternating-direction implicit (ADI) scheme. The intercoupling relation between thermocapillary force, buoyancy force, fluid property, heat transfer, and phase change, along with their effects on the induced flow patterns, are disclosed. The flow is classified into three types: thermocapillary, buoyancy, and combined convection. Among the three mechanisms, the combined convection simulates the experimental observations quite well, and the basic mechanism of the observed convection inside evaporating sessile drops is thus identified. It is disclosed that evaporation initiates unstable convection, while condensation always brings about a stable density distribution which eventually damps out all fluid disturbances. Another numerical model is presented to study the effect of boundary recession due to evaporation, and the 'peeling-off' effect (the removal of the surface layer of fluid by evaporation) is shown to be relevant.

Description: Two-dimensional solidification influenced by anisotropic heat conduction has been considered. The interfacial energy balance was derived to account for the heat transfer in one direction (x or y) depending on the temperature gradient in both the x and y directions. A parametric study was made to determine the effects of the Stefan number, aspect ratio, initial superheat, and thermal conductivity ratios on the solidification rate. Because of the imposed boundary conditions, the interface became skewed and sometimes was not a straight line between the interface position at the upper and lower adiabatic walls (spatially nonlinear along the height). This skewness depends on the thermal conductivity ratio k(yy)/k(yx). The nonlinearity of the interface is influenced by the solidification rate, aspect ratio, and k(yy/k(yx).

Description: Mercuric iodide X-ray detectors have been undergoing tests in a prototype scanning electron microscope system being developed for unmanned space flight. The detector program addresses the issues of geometric configuration in the SEM, compact packaging that includes separate thermoelectric coolers for the detector and FET, X-ray transparent hermetic encapsulation and electrical contacts, and a clean vacuum environment.

Description: A multiple-time-scale turbulence model of a single point closure and a simplified split-spectrum method is presented. In the model, the effect of the ratio of the production rate to the dissipation rate on eddy viscosity is modeled by use of the multiple-time-scales and a variable partitioning of the turbulent kinetic energy spectrum. The concept of a variable partitioning of the turbulent kinetic energy spectrum and the rest of the model details are based on the previously reported algebraic stress turbulence model. Example problems considered include: a fully developed channel flow, a plane jet exhausting into a moving stream, a wall jet flow, and a weakly coupled wake-boundary layer interaction flow. The computational results compared favorably with those obtained by using the algebraic stress turbulence model as well as experimental data. The present turbulence model, as well as the algebraic stress turbulence model, yielded significantly improved computational results for the complex turbulent boundary layer flows, such as the wall jet flow and the wake boundary layer interaction flow, compared with available computational results obtained by using the standard kappa-epsilon turbulence model.

Description: The emission from a gray radiating medium is analyzed for transient cooling in surroundings at a low temperature. The medium is rectangular with no variations in the direction normal to the cross section. The integral equation for the transient temperature distribution is solved numerically using a two-dimensional Gaussian integration subroutine. The emissive ability for a rectangle at uniform temperature is compared with that for transient cooling where the temperature distribution of the region has reached a fully developed shape, as shown by a separation of variables solution. The two solutions provide the upper and lower bounds for the emittance of a rectangle during transient cooling. The emittances for various aspect ratios are presented as a function of the mean length of the rectangle and are compared with results for a plane layer and a cylinder.

Description: Calculations for n-decane drops evaporating in a spherical cluster surrounded by unvitiated ambient air at atmospheric pressure were performed using two previously proposed cluster models. Both cluster models predict that turbulent transport effects are more important in the case of small clusters. This is due to the smaller volume to surface ratio and thus to the greater transport of hot unvitiated gas to the drops in order to promote evaporation. The results obtained are compared with those of two turbulent models for each one of the 'trapping factors' and similarity models.

Description: The effects of critical layer nonlinearity are considered on spatially growing oblique instability waves on nominally two-dimensional shear layers between parallel streams. The analysis shows that three-dimensional effects cause nonlinearity to occur at much smaller amplitudes than it does in two-dimensional flows. The nonlinear instability wave amplitude is determined by an integro-differential equation with cubic type nonlinearity. The numerical solutions to this equation are worked out and discussed in some detail. The numerical solutions always end in a singularity at a finite downstream distance.

Description: The effects of critical-layer nonlinearity on spatially growing oblique instability waves on compressible shear layers between two parallel streams are considered. The analysis shows that mean temperature nonuniformities cause nonlinearity to occur at much smaller amplitudes than it does when the flow is isothermal. The nonlinear instability wave growth rate effects are described by an integrodifferential equation which bears some resemblance to the Landau equation, in that it involves a cubic-type nonlinearity. The numerical solutions to this equation are worked out and discussed in some detail. Inviscid solutions always end in a singularity at a finite downstream distance, but viscosity can eliminate this singularity for certain parameter ranges.

Description: A simulation is performed of a passive scalar field convected by a rapidly fluctuating velocity field whose correlation time approaches zero. By using a code proposed in a previous study (Chasnov et al., 1988), the turbulence spectrum of the passive temperature field in the conductive subrange is determined. A theoretical model is proposed which explains the result obtained by representing the transfer of scalar variance by an eddy conductivity, whose correlation time is limited by the correlation time of the velocity field.

Description: The recent development of Scanning Tunneling Microscopy technology allows the application of electron tunneling to position detectors for the first time. The vacuum tunnel junction is one of the most sensitive position detection mechanisms available. It is also compact, simple, and requires little power. A prototype accelerometer based on electron tunneling, and other sensor applications of this promising new technology are described.

Description: Long's self-similar vortex is known to have two solutions for each supercritical value of the flow force. Each of those solutions is shown to have a double structure if the flow force is large. The inertial instabilities of one of those large-flow-force limit solutions are investigated, showing that they are related to the instabilities of the Bickley jet in one regime. However, the swirl in the vortex becomes important for long waves, very strongly modifying the sinuous and varicose, Bickley modes. The asymptotic results obtained agree well with the numerical solutions for the sinuous mode, but not for the varicose mode, the difficulty in the latter case being apparently due to mode jumping.

Description: Numerical techniques are developed to solve the Navier-Stokes equations for unsteady incompressible flow. The extension of the finite-difference Galerkin (FDG) method of Stephens et al. (1984) to the continuous-time case in two or three space dimensions is explained, and the numerical implementation of the method is discussed with particular attention to the staggered-MAC-grid primitive-variable discretization, the application of discrete mass balance to avoid problems inherent in FDG schemes, the direct interpretation of the FDG expansion variables as a discrete streamfunction, and a mass-balance approach to two-dimensional problems with throughflow or obstacles. Numerical results are presented graphically for the evolution of asymptotic steady flow in a driven cavity at Reynolds number 400, 1000, or 3200; good agreement with published experimental data is demonstrated, with accurate predictions of secondary-vortex formation from wall bubble recirculations at Reynolds number 1000.

Description: Flow-field measurements of unsteady turbulent flow downstream of a rotating spoked-wheel wake generator were performed in a short-duration light-piston tunnel, and the instantaneous-velocity data were phase averaged based on a signal synchronized with the bar-passing frequency. Mean axial velocities were found to agree well with those obtained from measurements behind a stationary cylinder and to be independent of both Reynolds and bar-passing Strouhal numbers. Reynolds stresses were found to be consistent with related cylinder-wake measurements, but were significantly higher than corresponding measurements obtained in large-scale research turbomachines. Phase-averaged triple velocity correlations were calculated from the digital velocity records, revealing the sign and the magnitude of skewness in the velocity probability density distributions for the two components.

Description: The bifurcation diagram corresponding to the Eckhaus stability curve has been constructed for the one-dimensional Swift-Hohenberg equation in a finite domain. Finite-amplitude solutions with particular spatial wavelength recover linear stability, as predicted by the Eckhaus curve, after a sequence of secondary bifurcations from the branch of solutions with this wavelength. No connectivity between the primary-solution branches is admissible if the stability predicted by this bifurcation diagram is to correspond to the prediction of the Eckhaus analysis. The Eckhaus curve does not exist if nonlinear couplings destroy this pattern. This is demonstrated by analysis of a coupled pair of Swift-Hohenberg equations.

Description: The effects of asymmetry in furnace temperature profile and pulling velocity on the crystal interface shape are demonstrated while neglecting the latent heat of solidification. It is seen that the furnace temperature profile may be varied in order to influence the shape of the melt-crystal interface. An exact thermal analysis is then performed on the Bridgman technique by including the latent heat of solidification as a source term. The exact temperature field required for yielding a flat melt-crystal interface is obtained. The earlier observation regarding the influence of furnace temperature profile on the interface shape is confirmed and a criterion for achieving a flat interface is obtained. Various furnace temperature profiles are selected and their corresponding melt-crystal interface results are presented.

Description: New methods have been developed to implement position sensors based on electron tunneling. The electron tunneling methods enable position to be detected with sub-Angstrom resolution using a compact mechanical structure and simple electronic control elements. A prototype accelerometer is the first sensor based on these principles; it shows reliable operation with a noise-limited sensitivity and a bandwidth of approximately 3 kHz. Based on these results, it is expected that tunnel sensors optimized for many applications will enable an entirely new class of sensors to be developed.

Description: Transient cooling by radiation is analyzed for a cylindrical region filled with axially flowing streams of drops that are becoming solidified. This is of interest for the dissipation of waste heat from orbiting power system in space. The drops absorb, emit, and scatter radiation, and the surroundings are at a lower uniform temperature. The radiative properties are assumed gray, and the scattering is isotropic. The radiating region is a two-phase mixture that remains at the melting temperature of the drops. Its temperature uniformity maintains a high emissive power as energy is lost. This is an advantage over a sensible heat radiator in which the temperature decreases, thereby reducing the emissive power. The results provide the axial length that remains two-phase and the fraction of energy dissipated within this length in which the emissive power has not decreased because of sensible cooling. It is also shown how the radial distribution of the axial velocity of the drops can be modified to increase this energy fraction.

Description: A number of successful applications of a spectral collocation method extended by a multi-domain patching technique are shown. The multi-domain technique can be used to improve resolution for problems with widely disparate scales, and to reduce the ill-conditioning of the spectral operators for problems in which a large number of points are required for distributed resolution. A new nonreflecting outflow boundary treatment for unsteady transition-to-turbulence simulations is also presented, which relies on the multi-domain technique. The role of multi-domain in improving the efficiency of such calculations is discussed.

Description: Spectral element methods are high-order weighted residual techniques based on spectral expansions of variables and geometry for the Navier-Stokes (NS) and transport equations. Here, practical aspects of these methods and their efficient implementation are examined, and several examples of flows in truly complex geometries are presented. The spectral element discretization for NS equations is introduced, and the convergence of the method is addressed. An efficient data management scheme is discussed in the context of parallel processing computations. The method is validated by comparing the spectral element solutions with the exact eigensolutions for the Orr-Sommerfeld equations in two and three dimensions. Computer-aided flow visualizations are presented for an impulsive flow past a sharp edge wedge. Three-dimensional states of channel flow disrupted by an array of cylindrical eddy promoters are studied, and the results of a direct simulation of the turbulent flow in a plane channel are presented.

Description: The Particle Measuring Systems (PMS) Forwared Scattering Spectrometer Probe (FSSP-100) which was flown on the ER-2 during the Airborne Antarctic Ozone Experiment for the measurement of particles in polar stratospheric clouds has been evaluated and calibrated. The sample volume of the probe per 10-s sampling period increases from 257 cu cm for 1-micron particles to 412 for 15-micron particles, but there is substantial uncertainty in this value. Limitations in the measurements from this instrument and possible corrections are discussed. The uncertainty in the total particle mass measured by the probe may be as large as + or - 100 percent. Recommendations are given for the processing of data from the FSSP used in this project.

Description: A description of the Axial Flow Turbine Research Facility (AFTRF) being built at the Turbomachinery Laboratory of the Pennsylvania State University is presented. The purpose of the research to be performed in this facility is to obtain a better understanding of the rotor/stator interaction, three dimensional viscous flow field in nozzle and rotor blade passages, spanwise mixing and losses in these blade rows, transport of wake through rotor passage, and unsteady aerodynamics and heat transfer of rotor blade row. The experimental results will directly feed and support the analytical and the computational tool development. This large scale low speed facility is heavily instrumented with pressure and temperature probes and has provision for flow visualization and laser Doppler anemometer measurement. The facility design permits extensive use of the high frequency response instrumentation on the stationary vanes and more importantly on the rotating blades. Furthermore it facilitates detailed nozzle wake, rotor wake, and boundary layer surveys. The large size of the rig also has the advantage of operating at Reynolds numbers representative of the engine environment.

Description: A quasi-three-dimensional analysis has been developed for unsteady rotor-stator interaction in turbomachinery. The analysis solves the unsteady Euler or thin-layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stress-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations are integrated in time using an explicit four-stage Runge-Kutta scheme with a constant time step. Implicit residual smoothing is used to increase the stability limit of the time-accurate computations. The scheme is described, and stability and accuracy analyses are given.

Description: A 3-D model was developed for simulating multistage turbomachinery flows using supercomputers. This average passage flow model described the time averaged flow field within a typical passage of a bladed wheel within a multistage configuration. To date, a number of inviscid simulations were executed to assess the resolution capabilities of the model. Recently, the viscous terms associated with the average passage model were incorporated into the inviscid computer code along with an algebraic turbulence model. A simulation of a stage-and-one-half, low speed turbine was executed. The results of this simulation, including a comparison with experimental data, is discussed.

Description: A recent study by Rocketdyne for NASA identified laser anemometry, using a compact optical head, as a feasible diagnostic instrument for the Space Shuttle Main Engine (SSME) Model Verification experiments. Physical Research, Inc. (PRI) is presently under contract from NASA Lewis to develop and deliver such a laser anemometer system. For this application, it is desired to place the laser at a remote distance from the engine, and use single mode polarization preserving fiber optics for the transmission of the laser light to and from the measurement head. Other requirements are given. Analytical and experimental tools are being used to develop the technologies required for the laser anemometer. These include finite element analysis of the optical head and vibration tests for various optical and mechanical components. Design of the optical head and the fiber optic connectors are driven by the temperature and vibration requirements for the measurement environment. Results of the finite element analysis and the vibration tests of the components are included. Conceptual design of the fiber optic launcher and the optical probe has also been complete. Detailed design of the probe as well as the fabrication and assembly of the components is in progress.

Description: Spacecraft range measurements have provided the most accurate tests, to date, of some relativistic gravitational parameters, even though the measurements were made with ranging systems having error budgets of about 10 meters. Technology is now available to allow an improvement of two orders of magnitude in the accuracy of spacecraft ranging. The largest gains in accuracy result from the replacement of unstable analog components with high speed digital circuits having precisely known delays and phase shifts.

Description: EDO Corporation, Barnes Engineering Division designed and constructed a high resolution thermal imaging system on contract to Lockheed for use in the SDI Star Lab. This employs a Pt Si CCD array which is sensitive in the spectral range of 3 to 5 microns. Star Lab will be flown in the Shuttle bay and consists basically of a large, reflecting, tracking telescope with associated sensors and electronics. The thermal imaging system is designed to operate in the focal plane of this telescope. The configuration of the system is illustrated. The telescope provides a collimated beam output which is focussed onto the detector array by a silicon objective lens. The detector array subtends a field of view of 1.6 degrees x 1.22 degrees. A beam switching mirror permits bypassing the large telescope to give a field of 4 degrees x 3 degrees. Two 8 position filter wheels are provided, and background radiation is minimized by Narcissus mirrors. The detector is cooled with a Joule-Thompson cryostat fed from a high pressure supply tank. This was selected instead of a more convenient closed-cycle system because of concern with vibration. The latter may couple into the extremely critical Starlab tracking telescope. The electronics produce a digitized video signal for recording. Offset and responsivity correction factors are stored for all pixels and these corrections are made to the digitized output in real time.

Description: A preliminary evaluation was made by ORNL of a two-color ratio pyrometer (TCRP) for temperature control in the Modular Electromagnetic Levitation (MEL) experiment. A discussion was presented by Eric Spjut at the 1987 NASA Non-Contact Temperature Measurement Workshop (NASA Conf. Publ. 2503, pp. 182-213) in which he described the non-linear characteristics of the time response of TCPs. Researchers replicated his model and results and note that the non-linear response behavior is minimized for small temperature steps at high temperatures. They then used the predicted response in a model for a proportional or integral feedback controller and predicted the control characteristics for heating and cooling a 5-mm diameter sphere of niobium at high (1500 to 2750 K) temperatures. The analysis shows that for a slow (25-ms) time response for a commercial RCRP, overshoots of several hundred kelvins will result from a 100-K decrease in the setpoint, and temperature tracking errors of 14 to 45 K will occur for control temperature ramps of 1000K/s. For a fast (greater than 0.1 ms) time response, the overshoot and ramp response errors are largely eliminated.

Description: The development of a noncontact temperature measurement device utilizing rotating analyzer ellipsometry is described. The technique circumvents the necessity of spectral emissivity estimation by direct measurement concomittant with radiance brightness. Using this approach, the optical properties of electromagnetically levitated liquid metals Cu, Ag, Au, Ni, Pd, Pt, and Zr were measured in situ at four wavelengths and up to 600 K superheat in the liquid. The data suggest an increase in the emissivity of the liquid compared with the incandescent solid. The data also show moderate temperature dependence of the spectral emissivity. A few measurements of the optical properties of undercooled liquid metals were also conducted. The data for both solids and liquids show excellent agreement with available values in the literature for the spectral emissivities as well as the optical constants.

Description: Commercially available types of infrared thermal imaging instruments, both viewers (qualitative) and imagers (quantitative) are discussed. The various scanning methods by which thermal images (thermograms) are generated will be reviewed. The performance parameters (figures of merit) that define the quality of performance of infrared radiation thermometers will be introduced. A discussion of how these parameters are extended and adapted to define the performance of thermal imaging instruments will be provided. Finally, the significance of each of the key performance parameters of thermal imaging instruments will be reviewed and procedures currently used for testing to verify performance will be outlined.

Description: A multiphase turbulence closure model is presented which employs one transport equation, namely the turbulence kinetic energy equation. The proposed form of this equation is different from the earlier formulations in some aspects. The power spectrum of the carrier fluid is divided into two regions, which interact in different ways and at different rates with the suspended particles as a function of the particle-eddy size ratio and density ratio. The length scale is described algebraically. A mass/time averaging procedure for the momentum and kinetic energy equations is adopted. The resulting turbulence correlations are modeled under less retrictive assumptions comparative to previous work. The closures for the momentum and kinetic energy equations are given. Comparisons of the predictions with experimental results on liquid-solid jet and gas-solid pipe flow show satisfactory agreement.

Description: Generally, two types of theory are used to describe the field equations for suspensions. The so-called postulated equations are based on the kinetic theory of mixtures, which logically should give reasonable equations for solutions. The basis for the use of such theory for suspensions is tenuous, though it at least gives a logical path for mathematical arguments. It has the disadvantage that it leads to a system of equations which is underdetermined, in a sense that can be made precise. On the other hand, the so-called averaging theory starts with a determined system, but the very process of averaging renders the resulting system underdetermined. A third type of theory is proposed in which the kinetic theory of gases is used to motivate continuum equations for the suspended particles. This entails an interpretation of the stress in the particles that is different from the usual one. Classical theory is used to describe the motion of the suspending medium. The result is a determined system for a dilute suspension. Extension of the theory to more concentrated systems is discussed.

Description: The optimal placement of discrete actuators and sensors is posed as a combinatorial optimization problem. Two examples for truss structures were used for illustration; the first dealt with the optimal placement of passive dampers along existing truss members, and the second dealt with the optimal placement of a combination of a set of actuators and a set of sensors. Except for the simplest problems, an exact solution by enumeration involves a very large number of function evaluations, and is therefore computationally intractable. By contrast, the simulated annealing heuristic involves far fewer evaluations and is best suited for the class of problems considered. As an optimization tool, the effectiveness of the algorithm is enhanced by introducing a number of rules that incorporate knowledge about the physical behavior of the problem. Some of the suggested rules are necessarily problem dependent.

Description: A sensor has been designed and tested for precise pointing applications. The device is able to sense extremely small rotary motion and is immune to cross-axis forces. The hardware and design characteristics of the torque sensor are presented. Test data, integrated control methodology, and future applications are included.

Description: The development of a prototype Control Moment Gyroscope (CMG) is discussed. Physical characteristics and the results of functional testing are presented to demonstrate the level of system performance obtained. Particular attention is given to how the man-rated mission requirement influenced the choice of the materials, fabrication, and design details employed. Comparisons are made of the measured system responses against the prediction generated by computer simulation.

Description: The effect of transient heat flux on heat flux sensor response and calibration is analyzed. A one dimensional case was studied in order to elucidate the key parameters and trends for the problem. It has the added advantage that the solutions to the governing equations can be obtained by analytic means. The analytical results obtained to date indicate that the transient response of a heat flux sensor depends on the thermal boundary conditions, the geometry and the thermal properties of the sensor. In particular it was shown that if the thermal diffusivity of the sensor is small, then the transient behavior must be taken into account.

Description: Measurements of heat flux to space shuttle main engine (SSME) turbine blade surfaces are being made in the Lewis heat flux calibration facility. Surface heat flux information is obtained from transient temperature measurements taken at points within the gauge. A 100-kW Vortek arc lamp is used as a source of thermal radiant energy. Thermoplugs, with diameters of about 0.190 cm and lengths varying from about 0.190 to 0.320 cm, are being investigated. The thermoplug is surrounded on all surfaces except the active surface by a pocket of air located in the circular annulus and under the back cover. Since the thermoplug is insulated, it is assumed that heat is conducted in a one-dimensional manner from the hot active surface to the cooler back side of the thermoplug. It is concluded that the miniature plug-type gauge concept is feasible for measurement of blade surface heat flux. It is suggested that it is important to measure heat flux near the hub on the suction surface and at the throat on SSME blades rotating in engines because stress and heat transfer coefficients are high in this region.

Description: A classification of infrared sensing instruments by type and application, listing commercially available instruments, from single point thermal probes to on-line control sensors, to high speed, high resolution imaging systems is given. A review of performance specifications follows, along with a discussion of typical thermographic display approaches utilized by various imager manufacturers. An update report on new instruments, new display techniques and newly introduced features of existing instruments is given.

Description: The design of a two color pyrometer with infrared optical fiber bundles for collection of the infrared radiation is described. The pyrometer design is engineered to facilitate its use for measurement of the temperature of small, falling samples in a microgravity materials processing experiment using a 100 meter long drop tube. Because the samples are small and move rapidly through the field of view of the pyrometer, the optical power budget of the detection system is severly limited. Strategies for overcoming this limitation are discussed.

Description: A multicolor imaging pyrometer was designed for accurately and precisely measuring the temperature distribution histories of small moving samples. The device projects six different color images of the sample onto a single charge coupled device array that provides an RS-170 video signal to a computerized frame grabber. The computer automatically selects which one of the six images provides useful data, and converts that information to a temperature map. By measuring the temperature of molten aluminum heated in a kiln, a breadboard version of the device was shown to provide high accuracy in difficult measurement situations. It is expected that this pyrometer will ultimately find application in measuring the temperature of materials undergoing radiant heating in a microgravity acoustic levitation furnace.

Description: The vaporization of a droplet, interacting with its neighbors in a non-dilute spray environment is examined as well as a vaporization scaling law established on the basis of a recently developed theory of renormalized droplet. The interacting droplet consists of a centrally located droplet and its vapor bubble which is surrounded by a cloud of droplets. The distribution of the droplets and the size of the cloud are characterized by a pair-distribution function. The vaporization of a droplet is retarded by the collective thermal quenching, the vapor concentration accumulated in the outer sphere, and by the limited percolative passages for mass, momentum and energy fluxes. The retardation is scaled by the local collective interaction parameters (group combustion number of renormalized droplet, droplet spacing, renormalization number and local ambient conditions). The numerical results of a selected case study reveal that the vaporization correction factor falls from unity monotonically as the group combustion number increases, and saturation is likely to occur when the group combustion number reaches 35 to 40 with interdroplet spacing of 7.5 diameters and an environment temperature of 500 K. The scaling law suggests that dense sprays can be classified into: (1) a diffusively dense cloud characterized by uniform thermal quenching in the cloud; (2) a stratified dense cloud characterized by a radial stratification in temperature by the differential thermal quenching of the cloud; or (3) a sharply dense cloud marked by fine structure in the quasi-droplet cloud and the corresponding variation in the correction factor due to the variation in the topological structure of the cloud characterized by a pair-distribution function of quasi-droplets.

Description: The use of a gyrotron for space power beaming, especially in the form of a lunar orbiting power station is discussed. The advantages of phased array power beaming, output power, and the design of a quasi-optical gyrotron are discussed.

Description: The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Description: The improvement in the sphere data processing, concerning the signal to noise ratio, is discussed. Frequency analysis of the radar data is effectuated. It reveals a specific frequency component in the radar angle error, which may originate from the tracking radar mechanism itself. An optimal (Wiener) filter is applied to the radar data in order to suppress the systematic angular error components selectively. Using this technique, a significant improvement in the signal to noise ratio is achieved. The resolution of sphere measurements, previously limited by the length of the polynomial filter in the sphere data processing algorithm, is improved.

Description: An explicit multistage Runge-Kutta type of time-stepping scheme is used for solving transonic flow past a transport type wing/fuselage configuration. Solutions for both Euler and Navier-Stokes equations are obtained for quantitative assessment of boundary layer interaction effects. The viscous solutions are obtained on both a medium resolution grid of approximately 270,000 points and a find grid of 460,000 points to assess the effects of grid density on the solution. Computed pressure distributions are compared with the experimental data.

Description: Some considerations toward developing numerical procedures for simulating viscous compressible flows are discussed. Both Navier-Stokes and boundary layer field methods are considered. Because efficient viscous-inviscid interaction methods have been difficult to extend to complex 3-D flow simulations, Navier-Stokes procedures are more frequently being utilized even though they require considerably more work per grid point. It would seem a mistake, however, not to make use of the more efficient approximate methods in those regions in which they are clearly valid. Ideally, a general purpose compressible flow solver that can optionally take advantage of approximate solution methods would suffice, both to improve accuracy and efficiency. Some potentially useful steps toward this goal are described: a generalized 3-D boundary layer formulation and the fortified Navier-Stokes procedure.

Description: The construction and development of the multi-component traversing system and associated control hardware and software are presented. A hydrogen bubble/laser sheet flow visualization technique was developed to visually study the characteristics of the mixing layers. With this technique large-scale rollers arising from the Taylor-Gortler instability and its interaction with the primary Kelvin-Helmholtz structures can be studied.

Description: The sizes and arrangement of the wind tunnel used for the experimentation are described. The specifications for the cold-wire anemometers, hot-wire anemometers, cold-wire rakes, and miniature 3-wire probe are proveded. The results of the experiment are briefly discussed.

Description: A very low Reynolds number turbulent boundary layer subject to an adverse pressure gradient is studied. The aim is to obtain highly accurate mean-flow and turbulence measurements under conditions that can be closely related to the numerical simulations of Philippe Spalart for the purposes of CFD validation. Much of the Boundary Layer Wind Tunnel was completely rebuilt with a new wider contraction and working section which will improve compatibility with the simulations. A unique sophisticated high-speed computer controlled 3-D probe traversing mechanism was integrated into the test section. Construction of the tunnel and traverse is discussed in some detail. The hardware is now complete, and measurements are in progress. The mean-flow data indicate that a suitably two-dimensional base flow was established. Automation of the probe positioning and data acquistion have led to a decreased running time for total pressure measurements. However, the most significant benefits are expected to occur when using hot-wire probes. Calibrations can be performed automatically and there is no need to handle fragile probes when moving between measuring stations. Techniques are being developed which require sampling of the signals from moving hot-wire probes on the basis of their position in the flow. Measurements can be made in high intensity turbulence by flying probes upstream at high speed so that the relative magnitude of the turbulent velocity fluctuations are reduced. In regions, where the turbulence intensity is not too large, the probe can also be repetitively scanned across very dense spatial grids in other directions. With this technique, a complete profile can be measured in about 1/3 the time and with a spatial density about 50 times that obtainable using a stationary probe.

Description: Thermal convection was proposed as a possible mechanism for generation and maintenance of turbulence in the inner accretion disk regime of the primordial solar nebula. It is of fundamental interest to design experiments with the basic physical features of the solar nebula conditions cannot be produced in the laboratory, numerical simulations of hydrodynamic flows, which have been very successful in describing aerodynamic flows, can be suitable modified to provide experimental data for solar nebula modelling. The goals are to modify an extant, proven hydrodynamics code with the most important features of the solar nebula and other thin accretion disks: bouyancy terms to generate convection, internal heating representing the release of gravitational potential energy, a variable gravity linearly proportional the the distance from the vertical midplane due to centrifugal balance, rapid rotation with axis aligned with gravity, and Keplerian rotational shear; to determine the effect that these features have on the turbulent convection by introducing them individually and to determine the cumulative nature of the turbulent convection for accretion disk conditions; and to model the convection and the turbulence. In this manner, prior solar nebula models can be tested and their deficiencies rectified.

Description: It is likely that turbulence played a major role in the evolution of the solar nebula, which is the flattened disk of dust and gas out of which our solar system formed. Relevant turbulent processes include the transport of angular momentum, mass, and heat, which were critically important to the formation of the solar system. This research will break ground in the modeling of compressible turbulence and its effects on the evolution of the solar nebula. The computational techniques which were developed should be of interest to researchers studying other astrophysical disk systems (e.g., active galactic nuclei), as well as turbulence modelers outside the astrophysics community.

Description: A research program for direct numerical simulations of compressible reacting flows is described. Two main research subjects are proposed: the effect of pressure waves on turbulent combustion and the use of direct simulation methods to validate flamelet models for turbulent combustion. The interest of a compressible code to study turbulent combustion is emphasized through examples of reacting shear layer and combustion instabilities studies. The choice of experimental data to compare with direct simulation results is discussed. A tentative program is given and the computation cases to use are described as well as the code validation runs.

Description: Outline of the research program and a recent progress in the studies of sheared turbulence are described. The research program reported is directed at two goals: (1) understanding of fundamental physics of organized structures in turbulent shear flows; and (2) development of phenomenological models of turbulence based on physical arguments. Three projects that were carried out are: (1) structure of sheared turbulence near a plane boundary; (2) distortion of turbulence by axisymmetric strain and dilation; and (3) study of energy transfer in turbulent shear flow.

Description: Vortex interactions and their role in turbulent flow are examined. The objectives are twofold. First, to use the existing axisymmetric code to study the annihilation process of colliding vortex rings and determine the relevance of this problem to similar 3-D phenomena. The second objective is to extend the code to three dimensions. The code under development is unique in that it can compute flows in a truly infinite domain (i.e., without periodic boundary conditions or approximations from truncating the domain). Because of this, the far field sound can be computed, and therefore, contribute to improved models of turbulence generated noise for this class of flows. Issues which can be addressed by the code include: effects of viscosity on mode selection in azimuthal breakdown of vortex rings (i.e., the Widnall instability); reconnection, the associated production of small scales, and the time scale of the process.

Description: The motivation for studying close vortex interactions is briefly discussed in the light of turbulence and coherent structures. Particular attention is given to the interaction known as reconnection. Two reconnection mechanisms are discussed. One is annihilation of vorticity by cross-diffusion, the other is an inviscid head-tail formation. At intermediate Reynolds numbers both mechanisms are operating.

Description: The objective of the present research was to extend the Direct Numerical Simulation (DNS) approach to particle-laden turbulent flows using a simple model of particle/flow interaction. The program addressed the simplest type of flow, homogeneous, isotropic turbulence, and examined interactions between the particles and gas phase turbulence. The specific range of problems examined include those in which the particle is much smaller than the smallest length scales of the turbulence yet heavy enough to slip relative to the flow. The particle mass loading is large enough to have a significant impact on the turbulence, while the volume loading was small enough such that particle-particle interactions could be neglected. Therefore, these simulations are relevant to practical problems involving small, dense particles conveyed by turbulent gas flows at moderate loadings. A sample of the results illustrating modifications of the particle concentration field caused by the turbulence structure is presented and attenuation of turbulence by the particle cloud is also illustrated.

Description: Direct numerical simulations are being performed on two different fluid flows in an attempt to discover the mechanism underlying the transition to turbulence in each. The first system is Taylor-Couette flow; the second, two-dimensional flow over an airfoil. Both flows exhibit a gradual transition to high-dimensional turbulence through low-dimensional chaos. The hope is that the instabilities leading to chaos will be easier to relate to physical processes in this case, and that the understanding of these mechanisms can then be applied to a wider array of turbulent systems.

Description: Cocke (1969) proved that in incompressible, isotropic turbulence the average material line (material surface) elements increase in comparison with their initial values. Good estimates of how much they increase in terms of the eigenvalues of the Green deformation tensor were rigorously obtained.

Description: The objective is to understand and extend a recent theory of turbulence based on dynamic renormalization group (RNG) techniques. The application of RNG methods to hydrodynamic turbulence was explored most extensively by Yakhot and Orszag (1986). An eddy viscosity was calculated which was consistent with the Kolmogorov inertial range by systematic elimination of the small scales in the flow. Further, assumed smallness of the nonlinear terms in the redefined equations for the large scales results in predictions for important flow constants such as the Kolmogorov constant. It is emphasized that no adjustable parameters are needed. The parameterization of the small scales in a self-consistent manner has important implications for sub-grid modeling.

Description: Use of the Smagorinsky eddy-viscosity formulation and related schemes for subgrid-scale parameterization of large eddy simulation models requires specification of a single length scale, earlier related by Lilly to the scale of filtering and/or numerical resolution. An anisotropic integration of the Kolmogoroff enstrophy spectrum allows generalization of that relationship to anisotropic resolution. It is found that the Deardorff assumption is reasonably accurate for small anisotropies and can be simply improved for larger values.

Description: In recent years codes that use the Navier-Stokes equations to compute aerodynamic flows have evolved from computing two-dimensional flows around simple airfoils to computing flows around full scale aircraft configurations. Most flows of engineering interest are turbulent and turbulence models are needed for their prediction. Yet, it is known that present turbulence models are adequate only for simple flows and do poorly in complicated flows such as three-dimensional separation, or large-scale unsteadiness. The same progress that allowed the development of these aerodynamic codes, namely the introduction of supercomputers, has allowed us to compute directly turbulent flows, albeit only for simple flows at moderate Reynolds numbers. These direct turbulence simulations provide us with detailed data that experimentalists were not able to measure. This work is motivated by the fact that data exists for developing better turbulence models and by the need for better models to compute flows of engineering interest. The objective is to develop turbulence models for engineering applications. The model categories that show promise for immediate use are on the two-equation level and the Reynolds-stress level.

Description: Proposed systems for Deep Space Network (DSN) microwave antenna holography are analyzed. Microwave holography, as applied to antennas, is a technique which utilizes the Fourier Transform relation between the complex far-field radiation pattern of an antenna and the complex aperture field distribution to provide a methodology for the analysis and evaluation of antenna performance. Resulting aperture phase and amplitude distribution data are used to precisely characterize various crucial performance parameters, including panel alignment, subreflector position, antenna aperture illumination, directivity at various frequencies, and gravity deformation. Microwave holographic analysis provides diagnostic capacity as well as being a powerful tool for evaluating antenna design specifications and their corresponding theoretical models.

Description: In the last two decades there have been extensive developments in computational unsteady transonic aerodynamics. Such developments are essential since the transonic regime plays an important role in the design of modern aircraft. Therefore, there has been a large effort to develop computational tools with which to accurately perform flutter analysis at transonic speeds. In the area of Computational Fluid Dynamics (CFD), unsteady transonic aerodynamics are characterized by the feature of modeling the motion of shock waves over aerodynamic bodies, such as wings. This modeling requires the solution of nonlinear partial differential equations. Most advanced codes such as XTRAN3S use the transonic small perturbation equation. Currently, XTRAN3S is being used for generic research in unsteady aerodynamics and aeroelasticity of almost full aircraft configurations. Use of Euler/Navier Stokes equations for simple typical sections has just begun. A brief history of the development of CFD for aeroelastic applications is summarized. The development of unsteady transonic aerodynamics and aeroelasticity are also summarized.

Description: A finite difference code was implemented for the compressible Navier-Stokes equations on the Connection Machine, a massively parallel computer. The code is based on the ARC2D/ARC3D program and uses the implicit factored algorithm of Beam and Warming. The codes uses odd-even elimination to solve linear systems. Timings and computation rates are given for the code, and a comparison is made with a Cray XMP.

Description: The hardware and software currently used for visualization of fluid dynamics at NASA Ames is described. The software includes programs to create scenes (for example particle traces representing the flow over an aircraft), programs to interactively view the scenes, and programs to control the creation of video tapes and 16mm movies. The hardware includes high performance graphics workstations, a high speed network, digital video equipment, and film recorders.

Description: The present investigation has focused on a computational methodology for the fundamental case of transition in channel flow, in which recently published experimental data are utilized both as a stimulus and as a measure of merit of the method. The research has proceeded along three avenues in parallel. The first task has consisted of the development and verification of a computer code which calculates the mean evolution of flow in a channel similar to the one employed experimentally by Blair and Anderson. An analytical test case was created for the dual purposes of code verification and of highlighting the interactions between the Reynolds stress and the mean velocity profile. This test case generated a Reynolds stress by the residue in the momentum equation which is produced by a typical analytical velocity profile. By a substitution of this Reynolds stress into the appropriate code module, the correctness of the code may be verified, along with the accuracy of the computational method. The second task pursued has involved the development of a triple layer model for the Reynolds stress profile, which was suggested and derived from experimental velocity profiles. It is demonstrated that the innermost length scale is based on the local friction velocity, the intermediate layer corresponds to the usual logarithmic law of the wall region in which the normalized Reynolds stress is approximately unity, and the outermost layer is represented by a closed mathematical form depending explicitly on the velocity profile in the wake region. The third task was comprised of scrutiny of the excellent databases developed by Blair and others, and the planning of its incorporation into the transition analysis. These extensive measurements indicate that turbulent statistics in the transition regime may be considered to alternate between laminar and fully turbulent types, the proportions of which are quantified by a measured intermittency function.