ALBERT

Description: The importance of understanding and modeling the unsteady flow phenomena in turbomachinery is discussed. Historical events in the application and development of gas turbines for aircraft propulsion are traced. Technology advancements over the years are highlighted with focus on the compression system components. Trends in compressor research within the National Advisory Committee for Aeronautics (NACA)/National Aeronautics and Space Administration (NASA) are noted. The impact of technology advancements on the increased occurrences of unsteady flow related problems in advanced engine development programs is discussed. The impact of the new and more demanding requirements being imposed on the propulsion system to meet advanced aircraft mission needs are also noted. Brief discussions on the present day understanding and modeling capability of the unsteady flow phenomena are presented to include discussions on rotating stall, surge, flutter, forced response and noise generation.

Description: The objective is to verify the capability of a cascade variable conductance heat pipe (CVCHP) system to provide precise temperature control of long life spacecraft without the need for a feedback heater or other power sources for temperature adjustment under conditions of widely varying power input and ambient environment. Solar energy is the heat source and space the heat sink for thermally loading two series connected variable conductance heat pipes. Electronics and power supply equipment requirements are minimal. A 7.5 V lithium battery supplies the power for thermistor type temperature sensors for monitoring system performance, and a 28 V lithium battery supplies power for valve actuation.

Description: The objective of this experiment is to evaluate the zero-g performance of a number of transverse flat plate heat pipe modules. Performance will include the transport capability of the pipes, the temperature drop, and the ability to maintain temperature over varying duty cycles and environments. Additionally, performance degradation, if any, will be monitored over the length of the Long Duration Exposure Facility (LDEF) mission. This information is necessary if heat pipes are to be considered for system designs where they offer benefits not available with other thermal control techniques.

Description: The principal objectives of the experiment are to determine zero-g start-up performance for conventional and diode low temperature heat pipes, to evaluate heat pipe performance in zero-g for an extended period of time, to determine zero-g transport capability of each heat pipe, and to determine diode operation, including forward conductance, turndown ratio, and transient behavior. Two heat pipes, a fixed conductance transporter heat pipe and a thermal diode heat pipe, are coupled with a radiant cooler system. Both pipes are charged with ethane. Also integrated with the radiator is a phase change material (PCM) canister which provides temperature stability during transport tests. N-heptane, which has a melting/freezing point of 182 K, is used as the PCM. The high heat capacity (28 W-hr of latent heat) provided by the canister permits high power heat pipe testing at constant temperature.

Description: Numerical experiments are used to study thermally driven flows which occur during vertical Bridgman crystal growth of a single component fluid. The solid-liquid interface was specified as parabolic and flow patterns were calculated for various insulation thicknesses, Grashof, Prandtl, and Biot numbers. When the melt is on top and the gravity vector is axially downward it is shown that flow persists as long as a radial temperature gradient is present. If the interface is convex, as viewed from the liquid, a single cell is observed. A concave interface exhibits multiple counterrotating cells. The insulation thickness and Grashof, Prandtl, and Biot numbers influence the flow in a quantitative manner.

Description: The stability of the parallel flow between a vertical crystal-melt interface and a vertical wall held at a temperature above the melting point of the crystal is analyzed for Prandtl numbers, P, ranging from 0.01 to 100. Three modes of instability occur: (1) a buoyant mode, (2) a shear mode, and (3) a coupled crystal-melt mode. The buoyant and shear modes are similar to those that occur for flow between two vertical rigid walls held at different temperatures. For Prandtl numbers greater than approximately two, the coupled crystal-melt mode occurs at a lower Grashof number than the other two modes. Specific results are given for succinonitrile (P = 22.8) and lead (P = 0.0225). These calculations and similar calculations for a cylindrical geometry were motivated by and are in general agreement with recent experiments on succinonitrile.

Description: The solid-liquid interface position and the temperature gradients in both the solid and liquid at the interface have been studied in a modified Bridgman-Stockbarger crystal growth furnace. These crystal growth factors have been studied as a function of ampoule translation rate, materials properties, and the size and temperature of a small auxiliary heater placed at the edge of the furnace hot zone. It has been found that the interface position with respect to a furnace reference point is essentially constant during a run for a low thermal conductivity material whereas the interface position changes continuously during a run with high thermal conductivity material. However, the ampoule translation rate and auxiliary heater conditions produce interface position changes in both high and low thermal conductivity materials.

Description: The motion of two and four rectilinear vortices inside a cylindrical pipe is studied under the restriction that the total circulation be zero. In the two-vortex case, it is shown that the motion is always periodic and an expression for the period is derived. In the four-vortex case, the motion is determined not to be periodic in general. However, a class of solutions where the motion is periodic is found. Several sample calculations of the vortex motion are included.

Description: Some turbulent solutions of the unaveraged Navier-Stokes equations (equations of fluid motion) are reviewed. Those equations are solved numerically in order to study the nonlinear physics of incompressible turbulent flow. Initial three-dimensional cosine velocity fluctuations and periodic boundary conditions are used in most of the work considered. The three components of the mean-square velocity fluctuations are initially equal for the conditions chosen. The resulting solutions show characteristics of turbulence such as the linear and nonlinear excitation of small-scale fluctuations. For the stronger fluctuations, the initially nonrandom flow develops into an apparently random turbulence. Thus randomness or turbulence can arise as a consequence of the structure of the Navier-Stokes equations. The cases considered include turbulence which is statistically homogeneous or inhomogeneous and isotropic or anisotropic. A mean shear is present in some cases. A statistically steady-state turbulence is obtained by using a spatially periodic body force. Various turbulence processes, including the transfer of energy between eddy sizes and between directional components, and the production, dissipation, and spatial diffusion of turbulence, are considered. It is concluded that the physical processes occurring in turbulence can be profitably studied numerically.

Description: An experiment was conducted to measure the heat transfer from a heated cylinder in crossflow in an array of circular cylinders. All cylinders had a length-to-diameter ratio of 3.0. Both in-line and staggered array patterns were studied. The cylinders were spaced 2.67 diameters apart center-to-center in both the axial and transverse directions to the flow. The row containing the heated cylinder remained in a fixed position in the channel and the relative location of this row within the array was changed by adding up to five upstream rows. The working fluid was nitrogen gas at pressures from 100 to 600 kPa. The Reynolds number range based on cylinder diameter and average unobstructed channel velocity was from 5,000 to 125,000. Turbulence intensity profiles were measured for each case at a point one half space upstream of the row containing the heated cylinder. The basis of comparison for all the heat transfer data was the single row with the heated cylinder. For the in-line cases the addition of a single row of cylinders upstream of the row containing the heated cylinder increased the heat transfer by an average of 50 percent above the base case. Adding up to five more rows caused no increase or decrease in heat transfer. Adding rows in the staggered array cases resulted in average increases in heat transfer of 21, 64, 58, 46, and 46 percent for one to five upstream rows, respectively. Previously announced in STAR as N82-19493

Description: The effects of mass injection and pressure gradients on the drag of surfaces were studied theoretically with the aid of boundary-layer and Navier-Stokes codes. The present investigation is concerned with the effects of spatially varying the injection in the case of flat-plate drag. Effects of suction and injection on wavy wall surfaces are also explored. Calculations were performed for 1.2 m long surfaces, one flat and the other sinusoidal with a wavelength of 30.5 cm. Attention is given to the study of the effect of various spatial blowing variations on flat-plate skin friction reduction, local skin friction coefficient calculated by finite difference boundary-layer code and Navier-Stokes code, and the effect of phase-shifting sinusoidal mass transfer on the drag of a sinusoidal surface.

Description: The numerical aspects of simulation unsteady flows which arise in turbomachinery are addressed. In particular the simulation of rotating stall and surge is discussed.

Description: A numerical algorithm that is second-order accurate in time has been developed for the conjugated problem of a separated, compressible flow field and a conductive solid body. The full two-dimensional time-dependent Navier-Stokes equations are coupled with the time-dependent energy equation for the solid body and are solved simultaneously. using implicit algorithms. The energy equation for the solid body may include arbitrarily distributed heat sources. The algorithm has been exmined for the case of two-dimensional supersonic compression-corner interaction, with a heat source embedded in the wall in the vicinity of the separation bubble and the attached boundary layer. The effect of the heat source on the flow field is studied for steady and transient cases.

Description: A power-law relationship between the average erosion rate and cumulative erosion is presented. Data analyses from Venturi, magnetostriction, and liquid-impingement devices conform to this unified relation. A normalization technique is also suggested for prediction purposes.

Description: The coil planet centrifuge designed by Ito employs flow of a single liquid phase, through a rotating coiled tube in a centrifugal force field, to provide a separation of particles based on sedimentation rates. Mathematical solutions are derived for the linear differential equations governing particle behavior in the coil planet centrifuge device. These solutions are then applied as the basis of a model for optimizing particle separations.

Description: For the past 25 years, there has existed in the Thermosciences Laboratory of the Mechanical Engineering Department of Stanford University a research program, primarily experimental, concerned with heat transfer through turbulent boundary layers. In the early phases of the program, the topics considered were the simple zero-pressure-gradient turbulent boundary layer with constant and with varying surface temperature, and the accelerated boundary layer. Later equilibrium boundary layers were considered along with factors affecting the boundary layer, taking into account transpired flows, flows with axial pressure gradients, transpiration, acceleration, deceleration, roughness, full-coverage film cooling, surface curvature, free convection, and mixed convection. A description is provided of the apparatus and techniques used, giving attention to the smooth plate rig, the rough plate rig, the full-coverage film cooling rig, the curvature rig, the concave wall rig, the mixed convection tunnel, and aspects of data reduction and uncertainty analysis.

Description: The current investigation has the objective to provide data which will make it possible to obtain a better estimate regarding the roughness drag for surface waviness. The data employed for this investigation were acquired in connection with a wavy wall study which was conducted as part of an overall program to reduce the skin friction of turbulent boundary layers in external flows. The results of the present investigation show that the low-speed roughness drag of small-amplitude sinusoidal wave trains having wavelengths of the order of the boundary-layer thickness is not only a function of h/lambda (h = total wave height, lambda = wavelength), but, in addition, is also a function of the roughness Reynolds number.

Description: Galerkin finite-element approximations are combined with computer-implemented perturbation methods for tracking families of solutions to calculate the steady axisymmetric flows in a differentially rotated cylindrical drop as a function of Reynolds number Re, drop aspect ratio and the rotation ratio between the two end disks. The flows for Reynolds numbers below 100 are primarily viscous and reasonably described by an asymptotic analysis. When the disks are exactly counter-rotated, multiple steady flows are calculated that bifurcate to higher values of Re from the expected solution with two identical secondary cells stacked symmetrically about the axial midplane. The new flows have two cells of different size and are stable beyond the critical value Re sub c. The slope of the locus of Re sub c for drops with aspect ratio up to 3 disagrees with the result for two disks of infinite radius computed assuming the similarity form of the velocity field. Changing the rotation ratio for exact counter-rotation ruptures the junction of the multiple flow fields into two separated flow families.

Description: A new system of approximation equations is derived for three-dimensional steady viscous compressible flows in which a primary-flow direction is present, but in which both transverse velocity components can be large. Previous approaches which address simplification of the steady Navier-Stokes equations are discussed, and a new approach is proposed. The transverse velocity vector which corrects a given potential flow has been decomposed into potential and rotational components. It is found that the potential-velocity vector may be assumed small, whereas the rotational-velocity vector may be assumed small, whereas the rotational velocity vector and hence the composite secondary flow can be of order unity. This assumption leads to a system of governing equations whose characteristic polynomial has a non-elliptic form for arbitrary Mach numbers. The resulting non-elliptic approximation equations can be solved as an initial/boundary-value problem. Computed results confirm the small scalar-potential approximation.

Description: A cooled porous insert in a curved wall has a specified spatially varying heat flux applied to one side. It is desired to control the distribution of coolant flow out through this curved surface so that the surface will be kept at a desired uniform temperature. The flow regulation is accomplished by shaping the surface through which the coolant enters the region to obtain the required variation of flow resistance within the region. The proper surface shape is found by solving a Cauchy boundary value problem. Analytical solutions are given in two dimensions for various shapes of the heated boundary subjected to different heating distributions.

Description: Galerkin finite-element approximations and Newton's method for solving free boundary problems are combined with computer-implemented techniques from nonlinear perturbation analysis to study solidification problems with natural convection in the melt. The Newton method gives rapid convergence to steady state velocity, temperature and pressure fields and melt-solid interface shapes, and forms the basis for algebraic methods for detecting multiple steady flows and assessing their stability. The power of this combination is demonstrated for a two-phase Rayleigh-Benard problem composed of melt and solid in a veritical cylinder with the thermal boundary conditions arranged so that a static melt with a flat melt-solid interface is always a solution. Multiple cellular flows bifurcating from the static state are detected and followed as Rayleigh number is varied. Changing the boundary conditions to approach those appropriate for the vertical Bridgman solidification system causes imperfections that eliminate the static state. The flow structure in the Bridgman system is related to those for the Rayleigh-Benard system by a continuous evolution of the boundary conditions.

Description: Computational models of turbulence in incompressible Newtonian fluids governed by the Navier-Stokes equations are reviewed. The governing equations are presented, and both direct and large-eddy-simulations are examined. Resolution requirements and numerical techniques of spatial representation, definition of initial and boundary conditions, and time advancement are considered. Results of simulations of homogeneous turbulence in uniform shear, the evolution of a turbulent mixing layer, and turbulent channel flow are presented graphically and discussed.

Description: Some physical, analytical, and computational aspects of viscous flow are examined with reference to examples of computed flows. The discussion of the physical aspects covers the development of important scales used to reference flow phenomena in laminar and turbulent shear layers; the usefulness of the concepts of circulation and vorticity; and some relatively large-scale organized structures that have recently been identified in transitional and fully developed turbulent flows. Among the analytical aspects discussed are a compact presentation of the compressible Navier-Stokes equations, the Reynolds-averaged form of these equations, and a simplified description of some forms of turbulent models. Finally, results of a simulation of the onset of transition, direct turbulence simulations, and large-eddy simulations are given.

Description: It is pointed out that most practical power generation and propulsion systems involve the burning of different types of fuel sprays, taking into account aircraft propulsion, industrial furnaces, boilers, gas turbines, and diesel engines. There has been a lack of data which can serve as a basis for spray model development and validation. A major aim of the present investigation is to fill this gap. Experimental apparatus and techniques for studying the characteristics of fuel sprays are discussed, taking into account two-dimensional still photography, cinematography, holography, a laser diffraction particle sizer, and a laser anemometer. The considered instruments were used in a number of experiments, taking into account three different types of fuel spray. Attention is given to liquid fuel sprays, high pressure pulsed diesel sprays, and coal-water slurry sprays.

Description: The solution by multigrid techniques of the steady inviscid compressible equations of gas dynamics, the Euler equations is investigated. Steady two dimensional transonic flow over an airfoil section is studied intensively. Most of the material is applicable to three dimensional flow problems of aerodynamic interest.

Description: Rotating baroclinic flow for Richardson number lower than about 1 is studied by means of a finite difference Navier-Stokes model assuming no variations except in the vertical plane that completely contains the density gradient vector. The horizontally infinite channel to which attention is given further assumes periodic boundary conditions at the vertical computational boundaries and no-slip conducting horizontal boundaries. Two configurations are considered. Symmetric baroclinic waves developed in the flows in a manner consistent with linear theory, and it is noted that the structures and energetics of the fully developed waves were dependent on the Prandtl number Pr. For Pr greater than 1, the conversion from potential energy to wave kinetic energy was direct, via temperature and vertical motion correlation, while for Pr of less than 1, conversion proceeded from potential energy to average kinetic energy by means of an induced meridional flow, and then to wave kinetic energy.

Description: Linear stability of the one-dimensional flow between infinite vertical coaxial cylinders induced by heating the inner cylinder is considered for various ratios (kappa) of the inner radius to the outer radius, and for Prandtl numbers (P) appropriate to air and water. For air with P = 0.71 the least stable disturbance is nonaxisymmetric for kappa less than 0.44 and is axisymmetric for kappa greater than 0.44, and in either case the instability is due to the action of the shear forces. For P = 3.5, the situation is similar, except that the asymmetric shear mode is superseded by an axisymmetric instability driven by buoyancy forces for kappa = 0.03-0.16. Wave speeds, wavenumbers, and critical Grashof numbers for these cases and for the case of zero Prandtl number are given.

Description: The extent of convective and radiative heating for a Saturn entry probe is investigated in the absence and presence of ablation mass injection. The flow in the shock layer is assumed to be axisymmetric, viscous and in local thermodynamic equilibrium. The importance of chemical nonequilibrium effects for both the radiative and convective nonblowing surface heating rates is demonstrated for prescribed entry conditions. Results indicate that the nonequilibrium chemistry can significantly influence the rate of radiative heating to the entry probes. With coupled carbon-phenolic ablation injection, the convective heating rates are reduced substantially. Turbulence has little effect on radiative heating but it increases the convective heating considerably.

Description: Ames Research Center has the lead role among NASA centers to conduct research in computational fluid dynamics. The past, the present, and the future prospects in this field are reviewed. Past accomplishments include pioneering computer simulations of fluid dynamics problems that have made computers valuable in complementing wind tunnels for aerodynamic research. The present facilities include the most powerful computers built in the United States. Three examples of viscous flow simulations are presented: an afterbody with an exhaust plume, a blunt fin mounted on a flat plate, and the space shuttle. The future prospects include implementation of the Numerical Aerodynamic Simulation Processing System that will provide the capability for solving the viscous flow field around an aircraft in a matter of minutes.

Description: The influence of high frequency excitations (HFE) on a fluid is investigated. The response to these excitations is decomposed in two parts: 'slow' motion, which practically remains unchanged during the vanishingly small period tau, and 'fast' motion whose value during this period is negligible in terms of displacements, but is essential in terms of the kinetic energy. After such a decomposition the 'slow' and 'fast' motions become nonlinearly coupled by the corresponding governing equations. This coupling leads to an 'effective' potential energy which imparts some 'elastic' properties to the fluid and stabilizes laminar flows.

Description: A solution method for finding the unknown solidification interface in manufacturing slab ingots as a continuous casting is presented, which involves a product solution in the potential plane and the use of conjugate harmonic functions. It is argued that the method may be more direct for some geometries than the Cauchy boundary value method. Moreover, the usefulness of the Cauchy boundary value method is demonstrated through the example of a nonsymmetric horizontal mold where the walls are offset to support the lower ingot boundary.

Description: A method for approximately analyzing the feedback between downstream and upstream edges in incompressible shear flow is described. The shear flow is modeled by a vortex sheet. Equations for resonance eigenvalues are derived. After the reduction of growth rate by finite shear layer thickness is allowed for, agreement is found between calculated resonances and those that have been observed experimentally.

Description: An implicit finite-difference method is presented for obtaining steady-state solutions to the time-dependent, conservative Euler equations for flows containing shocks. The method uses a two-point central-difference scheme for the flux derivatives with dissipation added at supersonic points via the retarded density concept. Application of the method to 1-dimensional nozzle flow equations for various combinations of subsonic and supersonic boundary conditions show the method to be very efficient. Residuals are typically reduced to machine zero in approximately 35 time steps for 50 mesh points. For 1-dimensional Euler calculations, it is shown that the scheme offers two advantages over the more widely-used three-point schemes. The first is in regard to application of boundary conditions, and the second relates to the fact that the two-point algorithm is well-conditioned for large time steps.

Description: This paper is concerned with small-amplitude, unsteady, vortical and entropic motion imposed on steady potential flows. It is restricted to the case where the spatial scale of the unsteady motion is small compared to that of the mean flow. Under such conditions, the unsteady motion may be influenced by viscosity even if the mean flow is not. An exact high-frequency (small-wavelength) solution is obtained for the small-amplitude viscous motion imposed on a steady potential flow. It generalizes the one obtained by Pearson (1959) for the homogeneous-strain case to the case of quasi-homogeneous strain. This result is used to study the effect of viscosity on rapidly distorted turbulent flows. Specific numerical results are given for a turbulent flow near a two-dimensional stagnation point.

Description: Past work on the influence of Mach number on the viscous and inviscid instability of flat-plate boundary layers is reviewed, and new spatial calculations are presented. These calculations support the previous view that viscosity is only stabilizing for both two- and three-dimensional first-mode waves above M1 = 3.0, and for second-mode waves at all Mach numbers. It is concluded that the calculations of Wazzan, Taghavi, and Keltner that show viscous instability at M1 = 6.0 for first-mode 50 deg waves, and at M1 = 3.0 for two-dimensional second-mode waves, are not correct.

Description: Previously cited in issue 05, p. 635, Accession no. A83-16649

Description: A theory is developed for the stagnation point boundary layer with injection under the hypothesis that turbulence is produced at the wall by injection. From the existing experimental heat transfer rate data obtained in wind tunnels, the wall mixing length is deduced to be a product of a time constant and an injection velocity. The theory reproduces the observed increase in heat transfer rates at high injection rates. For graphite and carbon-carbon composite, the time constant is determined to be 0.0002 sec from the existing ablation data taken in an arc-jet tunnel and a balistic range.

Description: Many of the gas turbine combustors in operation use multiple rows of dilution jets, and some have geometries that are different from circular holes. The data base available in literature is generally applicable to a single row of circular holes. Tests were performed with uniform mainstream conditions for several orifice plate configurations. Temperature and pressure measurements were made in the test section at 4 axial and 11 transverse stations. These measurements were made with a 60-element rake probe. Test results for some of these cases are discussed.

Description: Improved methods of predicting airfoil local metal temperatures require advances in the understanding of the physics and methods of analytically predicting the following four aerothermal loads: hot gas flow over airfoils, heat transfer rates on the gas-side of airfoils, cooling air flow inside airfoils, and heat transfer rates on the coolant-side of airfoils. A systematic building block research approach is being pursued to investigate these four areas of concern from both the experimental and analytical sides. Experimental approaches being pursued start with fundamental experiments using simple shapes and flat plates in wind tunnels, progress to more realistic cold and hot cascade tests using airfoils, continue to progress in large low-speed rigs and turbines and warm turbines, and finally, combine all the interactive effects in tests using real engines or real engine type turbine rigs. Analytical approaches being pursued also build from relatively simple steady two dimensional inviscid flow and boundary layer heat transfer codes to more advanced steady two and three dimensional viscous flow and heat transfer codes. These advanced codes provide more physics to model better the interactive effects and the true real-engine environment.

Description: The objective of this 36-month experimental and analytical program is to develop a heat transfer and pressure drop database, computational fluid dynamic techniques, and correlations for multipass rotating coolant passages with and without flow turbulators. The experimental effort will be focused on the simulation of configurations and conditions expected in the blades of advanced aircraft high pressure turbines so that the effects of Coriolis and buoyancy forces on the coolant side flow can be rationally included in the design of turbine blades.

Description: The effects of high-intensity, large-scale turbulence on turbulent boundary-layer heat transfer are studied. Flow fields were produced with turbulence intensities up to 40% and length scales up to several times the boundary layer thickness. In addition, three different types of turbulence will be compared to see whether they have the same effect on the boundary layer. The three are: the far field of a free jet, flow downstream of a grid, and flow downstream of a simulated gas turbine combustor. Each turbulence field will be characterized by several measures: intensity (by component), scale, and spectrum. Heat transfer will be measured on a 2.5 m long, 0.5 m wide flat plate using the energy-balance technique. The same plate will be used in each of the four flow fields; a low-turbulence tunnel for baseline data, and the three flow situations mentioned.

Description: The primary basis for heat transfer analysis of turbine airfoils is experimental data obtained in linear cascades. A detailed set of heat transfer coefficients was obtained along the midspan of a stator and a rotor in a rotating turbine stage. The data are to be compared to standard analyses of blade boundary layer heat transfer. A detailed set of heat transfer coefficients was obtained along the midspan of a stator located in the wake of a full upstream turbine stage. Two levels of inlet turbulence (1 and 10 percent) were used. The analytical capability will be examined to improve prediction of the experimental data.

Description: Advanced 3-D inelastic structural/stress analysis methods and solution strategies for more accurate and yet more cost-effective analysis of combustors, turbine blades, and vanes are being developed. The approach is to develop four different theories, one linear and three higher order with increasing complexities including embedded singularities. Progress in each area is reported.

Description: The objectives of this work are: (1) to extend the technique of direct numerical simulations to turbulent, chemically reacting flows, (2) to test the validity of the method by comparing computational results with laboratory data, and (3) to use the simulations to gain a better understanding of the effects of turbulence on chemical reactions. The effects of both the large scale structure and the smaller scale turbulence on the overall reaction rates are addressed. The relationship between infinite reaction rate and finite reaction rate chemistry is compared with some of the results of calculations with existing theories and laboratory data. The direct numerical simulation method involves the numerical solution of the detailed evolution of the complex turbulent velocity and concentration fields. Using very efficient numerical methods (e.g., pseudospectral methods), the fully nonlinear (possibly low pass filtered) equations of motion are solved and no closure assumptions or turbulence models are used. Statistical data are obtained by performing spatial, temporal, and/or ensemble averages over the computed flow fields.

Description: The present paper will compare temperature field measurements from selected cases in these investigations with distributions calculated with an empirical model based on assumed vertical profile similarity and superposition and with a 3-D elliptic code using a standard K-E turbulence model. The results will show the capability (or lack thereof) of the models to predict the effects of the principle flow and geometric variables.

Description: Fuel spray analyses which are a necessary input to the analytical modeling of the complex mixing and combustion processes which occur in advanced combustor systems are discussed. It is anticipated that by controlling fuel air reaction conditions, combustor temperatures can be better controlled, leading to improved combustion system durability. The capability to measure liquid droplet size, velocity, and number density throughout a fuel spray and to utilize this measurement technique in laboratory benchmark experiments was demonstrated. The experiment to characterize fuel sprays is described. The experiments and data are useful for application to and validation of turbulent flow modeling to improve the design systems of future advanced technology engines.

Description: Most fluid flows are turbulent rather than laminar and the reason for this was studied. One of the earliest explanations was that laminar flow is unstable, and the linear instability theory was first developed to explore this possibility. A series of early papers by Rayleigh produced many notable results concerning the instability of inviscid flows, such as the discovery of inflectional instability. Viscosity was commonly thought to act only to stabilize the flow, and flows with convex velocity profiles appeared to be stable. The investigations that led to a viscous theory of boundary layer instability was reported. The earliest application of linear stability theory to transition prediction calculated the amplitude ratio of the most amplified frequency as a function of Reynolds number for a Blasius boundary layer, and found that this quantity had values between five and nine at the observed Ret. The experiment of Schubauer and Skramstad (1947) completely reversed the prevailing option and fully vindicated the Gottingen proponents of the theory. This experiment demonstrated the existence of instability waves in a boundary layer, their connection with transition, and the quantitative description of their behavior by the theory of Tollmien and Schlichting. It is generally accepted that flow parameters such as pressure gradient, suction and heat transfer qualitatively affect transition in the manner predicted by the linear theory, and in particular that a flow predicted to be stable by the theory should remain laminar. The linear theory, in the form of the e9, or N-factor is today in routine use in engineering studies of laminar flow. The stability theory to boundary layers with pressure gradients and suction was applied. The only large body of numerical results for exact boundary layer solutions before the advent of the computer age by calculating the stability characteristics of the Falkner-Skan family of velocity profiles are given. When the digital computer reached a stage of development which permit the direct solution of the primary differential equations, numerical results were obtained from the linear theory during the next 10 years for many different boundary layer flows: three dimensional boundary layers; free convention boundary layers; compressible boundary layers; boundary layers on compliant walls; a recomputation of Falkner-Skan flows; unsteady boundary layers; and heated wall boundary layers.

Description: An iterative collocation technique is described for modeling implicit viscosity in three-dimensional incompressible wall bounded shear flow. The viscosity can vary temporally and in the vertical direction. Channel flow is modeled with a Fourier-Legendre approximation and the mean streamwise advection is treated implicitly. Explicit terms are handled with an Adams-Bashforth method to increase the allowable time-step for calculation of the implicit terms. The algorithm is applied to low amplitude unstable waves in a plane Poiseuille flow at an Re of 7500. Comparisons are made between results using the Legendre method and with Chebyshev polynomials. Comparable accuracy is obtained for the perturbation kinetic energy predicted using both discretizations.

Description: The effect of shoulder radiusing and grooving (longitudinally or circumferentially) the afterbodies of bluff bodies of reduce the base drag at low speeds is discussed. Shoulder radii as large as 2.75 body diameters are examined. Reynolds number based on body diameter varied from 20,000 to 200,000. Results indicate that increasing the shoulder radius to 2.75 body diameters can reduce the drag levels to those of a streamline body having 67 percent greater fineness ratio. For the relatively sharp shoulder case, body drag reductions as large as 50% are obtained using circumferential or longitudinal grooves.

Description: A detailed computerized sensitivity analysis of the triple hot-wire equations has been performed in order to delineate the uncertainties associated with measurements of the velocity components. Absolute and relative uncertainties for the instantaneous hot-wire outputs are calculated as functions of roll and pitch angles, based on a constant probability combination of the uncertainties in the inputs. From the results, it is concluded that the small inherent difficulties associated with the triple hot-wire data do not reflect artifacts introduced by the data processing. Fixed errors present in the V and W channels of the output are due to the nonzero measuring volume of the triple wire probe, and are entirely predictable.

Description: A computationally efficient method is proposed for obtaining fine-mesh solutions of the Navier-Stokes equations for high-Re flow, including separation. The method involves implementing a multi-grid solution procedure with suitably chosen elements, such that the resulting overall solution technique is efficient as well as robust. The robustness and efficiency of the solution technique are demonstrated by applying it to three model problems: flow in a driven cavity, downstream asymptotic flow in curved ducts of square and polar sections, and Newmann boundary-value problem in clustered curvilinear orthogonal coordinates.

Description: The differential formulations and computational techniques currently used for the incompressible Navier-Stokes (NS) and parabolic Navier-Stokes (PNS) equations are reviewed. In particular, attention is given to problems associated with the choice of difference equations, the method of solution and the choice of algorithm, the coupling of dependent variables and discretized equations, the application of boundary conditions, and grid generation. A new composite velocity NS and PNS formulation in (u,v,p) variables is presented, and the applicability of a 'forward' difference global pressure iteration for the (u,v,p) PNS system is demonstrated.

Description: An attempt was made to model the so called leading edge vortex which forms over the leading edge of delta wings at high angles of attack. A simplified model was considered, namely that of a two-dimensional, inviscid, incompressible steady flow around a flat plate at an angle of attack with a stationary vortex detached on top, as well as a sink to simulate the strong spanwise flow. The results appear to agree qualitatively with experiments. A comparison was also made between the lift and the drag of this model and the corresponding results for two classical solutions: (1) that of totally attached flow over the plate with the Kutta condition satisfied at the trailing edge only: and (2) the Helmholtz solution of totally separated flow over the plate.

Description: Zonal concepts are utilized to delineate regions of application of three-dimensional boundary layer (DBL) theory. The zonal approach requires three distinct analyses. A modified version of the 3-DBL code named TABLET is used to analyze the boundary layer flow. This modified code solves the finite difference form of the compressible 3-DBL equations in a nonorthogonal surface coordinate system which includes coriolis forces produced by coordinate rotation. These equations are solved using an efficient, implicit, fully coupled finite difference procedure. The nonorthogonal surface coordinate system is calculated using a general analysis based on the transfinite mapping of Gordon which is valid for any arbitrary surface. Experimental data is used to determine the boundary layer edge conditions. The boundary layer edge conditions are determined by integrating the boundary layer edge equations, which are the Euler equations at the edge of the boundary layer, using the known experimental wall pressure distribution. Starting solutions along the inflow boundaries are estimated by solving the appropriate limiting form of the 3-DBL equations.

Description: The experimental contract objective is to provide a complete set of benchmark quality data for the flow within a large rectangular turning duct. The data are to be used to evaluate and verify three-dimensional internal viscous flow models and computational codes. The analytical contract objective is to select such a computational code and define the capabilities of this code to predict the experimental results. Details of the proper code operation will be defined and improvements to the code modeling capabilities will be formulated.

Description: Work is currently underway to develop and characterize an analytical approach, based on boundary layer theory, for predicting the effects of leading edge (showerhead) film cooling on downstream gas side heat transfer rates. Parallel to this work, experiments are being conducted to build a relevant data base for present and future methods verification.

Description: A Navier-Stokes analysis employing the time-dependent Linearized Block Implicit scheme (LBI) was applied to two-dimensional and three-dimensional transonic turbulent cascade flows. In general, the geometrical configuration of the turbine blade impacts both the grid construction procedure and the implementation of the numerical algorithm. Since modern turbine blades of interest are characterized by very blunt leading edges, rounded trailing edges and high stacking angles, a robust grid construction procedure is required that can accommodate the severe body shape while resolving regions of large flow gradients. A constructive O-type grid generation technique, suitable for cascades with rounded trailing edges, was developed and used to construct the C3X turbine cascade coordinate grid. Two-dimensional calculations were performed employing the Navier-Stokes procedure for the C3X turbine cascade, and the predicted pressure coefficients and heat transfer rates were compared with the experimental data. Three-dimensional Navier-Stokes calculations were also performed.

Description: An experimental study of mixing downstream of axial and swirling coaxial jets is being conducted to obtain data for the evaluation and improvement of turbulent transport models currently employed in a variety of computational procedures used throughout the propulsion community. Effort was directed toward the acquisition of length scale and dissipation rate data that will provide more accurate inlet boundary conditions for the computational procedures and a data base to evaluate the turbulent transport models in the near jet region where recirculation does not occur. Mass and momentum turbulent transport data with a blunt inner-jet inlet configuration will also be acquired.

Description: Combustor models for the aircraft gas turbine industry have been obtained because of the need to reduce the costs of developing improved performance and more durable engines. A few years ago, it became apparent that the mass concentration and velocity predictions provided by the computer codes were not representing the data measured in some confined recirculating flows. It is pointed out that errors in the mass concentration distribution are an especially serious problem because of their influence on the heat release, temperature, and reactant distributions. Combined mass and momentum turbulent transport experiments with swirling and nonswirling flow have been conducted with the objective to obtain an experimental data base which can be used to evaluate and improve the turbulent transport submodes employed in the aerothermal models. The present paper is mainly concerned with the overall characteristics of the mass turbulent transport processes in complex flows with recirculation and the deficiencies of the conventional models.

Description: Experimental studies have been conducted with the aim to obtain a better understanding of the fluid dynamics of mixing in gas turbine combustors, and solid fuel ramjet combustors subject to spin. The present investigation represents a continuation of studies conducted by So et al. (1984). It is also concerned with the verification of some conclusions reported by Ahmed et al. (1984). Attention is given to the experimental facility and instrumentation, centerline concentration measurements, mean concentration profiles, and a comparison of concentration and axial velocity results in the case of swirling flow.

Description: The acoustical excitation of shear layers is investigated. Acoustical excitation causes the so-called orderly structures in shear layers and jets. Also, the deviations in the spreading rate between different shear layer experiments are due to the same excitation mechanism. Measurements in the linear interaction region close to the edge from which the shear layer is shed are examined. Two sets of experiments (Houston 1981 and Berlin 1983/84) are discussed. The measurements were carried out with shear layers in air using hot wire anemometers and microphones. The agreement between these measurements and the theory is good. Even details of the fluctuating flow field correspond to theoretical predictions, such as the local occurrence of negative phase speeds.

Description: The thermal control requirements of a large space station are considered. Motivations for advanced thermal technology are discussed. Two test programs, designed to evaluate the analytical and theoretical basis from which thermal technology directions are determined, are described.

Description: An experimental study of mixing downstream of axial and swirling coaxial jets is being conducted to obtain data for the evaluation and improvement of turbulent transport models currently employed in a variety of computational procedures used throughout the propulsion community. The axial coaxial jet study was completed under Phase 1. The swirling coaxial jet study, which is the subject of this paper, was conducted under Phase 2 of the contract. A TEACH code was acquired, checked out for several test cases, and is reported. A study to measure length scales and to obtain a limited number of measurements with a blunt trailing edge inlet is being conducted under Phase 3 of the contract.

Description: In recent years, there has been a growing interest in the development and use of mathematical models for the simulation of fluid flow, heat transfer and combustion processes in engineering equipment. The equations representing the multi-dimensional transport of mass, momenta and species are numerically solved by finite-difference or finite-element techniques. However despite the multiude of differencing schemes and solution algorithms, and the advancement of computing power, the calculation of multi-dimensional flows, especially three-dimensional flows, remains a mammoth task. The following discussion is concerned with the author's recent work on the construction of accurate discretization schemes for the partial derivatives, and the efficient solution of the set of nonlinear algebraic equations resulting after discretization. The present work has been jointly supported by the Ramjet Engine Division of the Wright Patterson Air Force Base, Ohio, and the NASA Lewis Research Center.

Description: Several laser measurement methods are being studied to provide the capability to make droplet size and velocity measurements under a variety of spray conditions. The droplet sizing interferometer (DSI) promises to be a successful technique because of its capability for rapid data acquisition, compilation and analysis. Its main advantage is the ability to obtain size and velocity measurements in air-fuel mixing studies and hot flows. The existing DSI at NASA Lewis is a two-color, two-component system. Two independent orthogonal measurements of size and velocity components can be made simultaneously. It also uses an off-axis large-angle light scatter detection. The fundamental features of the system are optics, signal processing and data management system. The major component includes a transmitter unit, two receiver units, two signal processors, two data management systems, two Bragg cell systems, two printer/plotters, a laser, power supply and color monitor.

Description: An experimental program to study the characteristics of inhomogeneous jets in confined swirling flows to obtain detailed and accurate data for the evaluation and improvement of turbulent transport modeling for combustor flows is discussed. The work was also motivated by the need to investigate and quantify the influence of confinement and swirl on the characteristics of inhomogeneous jets. The flow facility was constructed in a simple way which allows easy interchange of different swirlers and the freedom to vary the jet Reynolds number. The velocity measurements were taken with a one color, one component DISA Model 55L laser-Doppler anemometer employing the forward scatter mode. Standard statistical methods are used to evaluate the various moments of the signals to give the flow characteristics. The present work was directed at the understanding of the velocity field. Therefore, only velocity and turbulence data of the axial and circumferential components are reported for inhomogeneous jets in confined swirling air flows.

Description: Turbulent sudden expansion flows are of significant theoretical and practical importance. Such flows have been the subject of extensive analytical and experimental study for decades, but many issues are still unresolved. Detailed information on reacting sudden expansion flows is very limited, since suitable measurement techniques have only been available in recent years. The present study of reacting flow in an axisymmetric sudden expansion was initiated under NASA support in December 1983. It is an extension of a reacting flow program which has been carried out with Air Force support under Contract F33615-81-K-2003. Since the present effort has just begun, results are not yet available. Therefore a brief overview of results from the Air Force program will be presented to indicate the basis for the work to be carried out.

Description: Techniques yielding simultaneous, multiple-point measurements of velocity in reacting or nonreacting flow fields have the potential to significantly impact basic and applied studies of fluid mechanics. This research program is aimed at investigating several candidate schemes which could provide such measurement capability. The concepts under study have in common the use of a laser source (to illuminate a column, a grid, a plane or a volume in the flow) and the collection of light at right angles (from Mie scattering, fluorescence, phosphorescence or chemiluminescence) using a multi-element solid-state camera (100 x 100 array of photodiodes). The work will include an overview and a status report of work in progress with particular emphasis on the method of Doppler-modulated absorption.

Description: Studies are concerned with experimental and theoretical research on 2-D axisymmetric geometries under low speed, nonreacting, turbulent, swirling flow conditions. The flow enters the test section and proceeds into a larger chamber (the linear expansion ratio D/d = 2, 1.5 and 1) via a sudden or gradual expansion (side wall angle alpha = 90 and 45 degrees). A weak or strong nozzle (of area ratio A/a = 2 and 4) may be positioned downstream at x/D = 2 to form a contraction exit to the test section. Inlet swirl vanes are adjustable to a variety of vane angles with values of theta = 0, 38, 45, 60 and 70 degrees being emphasized. The objective is to determine the effect of these parameters on isothermal flow field patterns, time mean velocities and turbulence quantities, and to establish an improved simulation in the form of a computer prediction code equipped with a suitable turbulence model. The goal of the on going research is to perform experiments and complementary computations with the idea of doing the necessary type of research that will yield improved calculation capability. This involves performing experiments where time mean turbulence quantities are measured and taking input conditions and running an existing prediction code for a variety of test cases so as to compare predictions against experiment.

Description: Three dimensional combustor calculations are currently stretching the computer hardware capabilities and the computing budgets of gas turbine manufacturers. One of the main reasons for this relates to the large number of complex physical processes occurring in the combustor. Airflow, fuel spray, reaction kinetics, flame radiation, and not the least of which, turbulence must be modeled and the related differential equations solved. Discussions in this conference will address methods to improve the accuracy of combustor flow field calculations and methods to speed the convergence of the modeled equations. This report will focus on aspects of merging these two new technologies. The improved accuracy discretization schemes have a negative impact on the speed of convergence of the modeled equations that the improved solution algorithms may not overcome. A description of the causes of this problem and potential solutions will be examined.

Description: A proposed two equation turbulence model for incompressible dilute two phase flows was validated and extended for steady incompressible two phase flow including phase change. The model was tested for the flow of a turbulent axisymmetric gaseous jet laden with multisize evaporating liquid droplets. Predicted results include distributions of the mean velocity; volume fractions of different phases concentration of the evaporated material in the carrier phase; turbulence intensity and shear stress of the carrier phase; droplet diameter distribution; and the jet spreading rate. Results are analyzed based on a qualitative comparison with the corresponding single phase jet flow.

Description: Techniques for studying the events directly preceding combustion in the liquid fuel sprays are being used to provide information as a function of space and time on droplet size, shape, number density, position, angle of flight and velocity. Spray chambers were designed and constructed for: (1) air-assist liquid fuel research sprays; (2) high pressure and temperature chamber for pulsed diesel fuel sprays; and (3) coal-water slurry sprays. Recent results utilizing photography, cinematography, and calibration of the Malvern particle sizer are reported. Systems for simultaneous measurement of velocity and particle size distributions using laser Doppler anemometry interferometry and the application of holography in liquid fuel sprays are being calibrated.

Description: Equations governing the incident heat flux distribution beneath a lamp-reflector system were developed. Analysis of a particular radiant heating facility showed good agreement between theory and experiment when a lamp power loss correction was used. In addition, the theory was employed to estimate thermal disruption in the radiation field caused by a protruding probe.

Description: External heat-transfer rates were measured on a flat plate in the junction region of a circular cylinder mounted normal to the plate at a location 90 deg from the stagnation point. This configuration simulates the junction of the shroud with gas-turbine vanes. Heat-transfer results are presented for laminar, transition, and turbulent boundary layers for a Mach number of 0.14 with gas temperatures of approximately 750 deg R over a flat plate at room temperature. The measurements were made in air for a unit Reynolds number, Re/cm, range of 11,000 to 58,000. Heat-transfer measurements were conducted in the 70-ft long, 4-in. diameter shock tube. A shock-wave reflection technique was used to produce a flow Mach number of 0.14. Thin-film platinum heat gages were mounted on the flat plate and along the line of the stagnation point of the cylinder to measure the local heat flux in the junction region. The experimental heat-transfer data were correlated with the laminar and turbulent boundary-layer theories for the flat plate. With the cylinder the heat fluxes on the flat plate were greatly increased in the junction region compared to the heat flux for the plate alone.

Description: The results of laser anemometer and flow visualization based fluid mechanics studies of porous media are used to obtain heat transfer coefficients for porous materials. Average pore flow Re ranging from 0.16-700 were examined. Darcy, inertial steady laminar, unsteady laminar and turbulent flow regimes were detected. A passage length model was devised to derive the heat transfer coefficient. Sample data from flows through porous metals composed of powders and fibers validated the passage length for Darcy and inertial flow regimes. Unsteady laminar and turbulent flow coefficients require the identification of new parameters.

Description: The basic heat pipe principle is employed to provide a self-contained passively cooled probe that may be placed into a high temperature environment. The probe consists of an evaporator region of a heat pipe and a sensing instrument. Heat is absorbed as the working fluid evaporates in the probe. The vapor is transported to the vapor space of the condenser region. Heat is dissipated from the condenser region and fins causing condensation of the working fluid, which returns to the probe by gravity and the capillary action of the wick. Working fluid, wick and condenser configurations and structure materials can be selected to maintain the probe within an acceptable temperature range.

Description: For the one-dimensional nozzle flow equations, the number of analytical boundary conditions that may be applied corresponds to the number of exterior characteristics reaching the boundary points. Since the slopes of the characteristics are the characteristic speeds u, u+c, and u-c, one concludes that two analytical boundary conditions are required at a subsonic inflow boundary and one at a subsonic outflow boundary. However, using this guideline alone can lead to nonunique solutions in subsonic cases. This paper presents a rule for selecting analytical boundary conditions for the dependent variables which will yield unique solutions if they exist.

Description: Phase averaged measurements of the axial mass velocity, radial velocity, and the product of the two have been made in an artificially excited cold air jet. The radial velocity fluctuations had a lower amplitude for both the full wave and phase averaged fluctuations compared to the axial mass velocity fluctuations. The Reynolds stress tended to decay faster in amplitude than the axial mass velocity and radial velocity fluctuations for both the full wave and phase averaged components. The radial velocity fluctuations possessed higher levels of coherence than the axial mass velocity fluctuation and the Reynolds stress. When the jet was excited at the frequency of St = 0.474, a high amplitude component at that frequency dominated the entire phase averaged flowfield from X/D = 1 to 10. This fundamental frequency was accompanied by several relatively high amplitude harmonics. These harmonics were more pronounced in the radial velocity fluctuations than for the axial fluctuations or the Reynolds stress.

Description: Experimental measurements of turbulent mass transit in a bluff-body diffusion-flame combustor (using CO2 instead of fuel) are analyzed to evaluate the accuracy of physical models used in flow computations. The data of Lightman and Magill (1981), Lightman et al. (1983), and Roquemore et al. (1983) are used to calculate apparent turbulent Schmidt numbers (TSN) for a series of flow conditions by a modified TEACH viscous-flow code. The modeling principles, calculation grid, and boundary conditions are discussed, and the results are presented in graphs comparing calculated and measured values. It is shown that models such as the two-equation (K - epsilon) model which use a single value of the TSN are inappropriate for conditions typical of gas-turbine combustors: local TSN variations due to turbulence behavior must be taken into account in improved models.

Description: Structure measurements were completed within the dilute portion of axisymmetric nonevaporating sprays (SMD of 30 and 87 microns) injected into a still air environment, including: mean and fluctuating gas velocities and Reynolds stress using laser-Doppler anemometry; mean liquid fluxes using isokinetic sampling; drop sizes using slide impaction; and drop sizes and velocities using multiflash photography. The new measurements were used to evaluate three representative models of sprays: (1) a locally homogeneous flow (LHF) model, where slip between the phases was neglected; (2) a deterministic separated flow (DSF) model, where slip was considered but effects of drop interaction with turbulent fluctuations were ignored; and (3) a stochastic separated flow (SSF) model, where effects of both interphase slip and turbulent fluctuations were considered using random sampling for turbulence properties in conjunction with random-walk computations for drop motion. The LHF and DSF models were unsatisfactory for present test conditions-both underestimating flow widths and the rate of spread of drops. In contrast, the SSF model provided reasonably accurate predictions, including effects of enhanced spreading rates of sprays due to drop dispersion by turbulence, with all empirical parameters fixed from earlier work.

Description: Previously cited in issue 08, p. 1213, Accession no. A82-22085

Description: Estimation of the integral scale and intensity of a generic turbulence record is treated as a statistical problem of parameter estimation. Properties of parameter estimators and the method of maximum likelihood are reviewed. Likelihood equations are derived for estimation of the integral scale and intensity applicable to a general class of turbulence spectra that includes the von Karman and Dryden transverse and longitudinal spectra as special cases. The method is extended to include the Bullen transverse and longitudinal spectra. Coefficients of variation are given for maximum likelihood estimates of the integral scale and intensity of the von Karman spectra. Application of the method is illustrated by estimating the integral scale and intensity of an atmospheric turbulence vertical velocity record assumed to be governed by the von Karman transverse spectrum.

Description: A technique for fabricating uniform thin-wall metallic heat-transfer models and which simulates a Shuttle thermal protection system tile is described. Two 6- by 6- by 2.5-in. tiles were fabricated to obtain local heat transfer rates. The fabrication process is not limited to any particular geometry and results in a seamless thin-wall heat-transfer model which uses a one-wire thermocouple to obtain local cold-wall heat-transfer rates. The tile is relatively fragile because of the brittle nature of the material and the structural weakness of the flat-sided configuration; however, a method was developed and used for repairing a cracked tile.

Description: A two-point finite difference unsteady laminar and turbulent boundary layer computational method was used to investigate the properties of the flow around an airfoil (NACA 0012) oscillating through angles of attack up to 18 degrees, for reduced frequencies of 0.01 and 0.20. The unsteady potential flow was determined using the unsteady potential flow method of Geissler. The influence of transition location on stal behavior was investigated, using both experimentally determined transition information, and transition located at the pressure peak; the results show the need for viscous-inviscid interaction in future computation of such flows.

Description: The status of a study to investigate the convective structures found in Doppler lidar windfield data is reported. Projected research is mentioned briefly.

Description: The development of a space-compatible He(3) refrigerator would provide a significant improvement in several areas of research in the 0.3 to 1 K temperature range. There are several methods of achieving these temperatures on Earth: He(3) refrigeration, dilution refrigeration, and adiabatic demagnetization refrigeration. The progress of adapting He(3) refrigeration for use in space is described. Various cycles and possible embodiments of He(3) refrigerators are described. Also included is an analysis of the liquid confinement and liquid-vapor phase-separation system. A possible configuration is then analyzed. Finally, the results of ground-based experiments will be discussed.

Description: A procedure is presented, as well as some results, to calculate the flow over the winged orbiter. This necessitates the use of two computer codes. A parabolized marching Navier-Stokes code is used to obtain the solution up to the bow shock-wing shock interaction region and for the region after the interaction. An unsteady Navier-Stokes code is to be used in the region of the shock interaction. Only resuls for the marching code are presented. For the flow conditions calculated, M infinity = 7.9, alpha = 25 deg, T(wall) = 540 R, Re(L) = 60728 per inch, laminar or turbulent, the PNS code was marched up to an X/L = 0.7 which is where the bow shock-wing shock interaction region occurs.

Description: A 6-ft by 6-ft wind tunnel section was modification to accommodate the 7-ft wide NASA dual-jet flate model in an effort to determine the effects of nonuniform and/or noncircular jet exhaust profiles on the pressure field induced on a nearby surface. Tests completed yield surface pressure measurements for a 90 deg circular injector producing exit profiles representative of turbofan nozzles (such as the TF-34 nozzle). The measurements were obtained for both tandem and side-by-side jet configurations, jet spacing of S/D =2, and velocity ratios of R=2.2 and 4.0. Control tests at the same mass flow rate but with uniform exit velocity profiles were also conducted, for comparison purposes. Plots for 90 deg injection and R=2.2 show that the effects of exit velocity profile nonuniformity are quite significant.

Description: A Fourier-Chebyshev spectral method for the incompressible Navier-Stokes equations is described. It is applicable to a variety of problems including some with fluid properties which vary strongly both in the normal direction and in time. In this fully spectral algorithm, a preconditioned iterative technique is used for solving the implicit equations arising from semi-implicit treatment of pressure, mean advection and vertical diffusion terms. The algorithm is tested by applying it to hydrodynamic stability problems in channel flow and in external boundary layers with both constant and variable viscosity.

Description: The two dimensional boundary layer on a concave wall is centrifugally unstable with respect to vortices aligned with the basic flow for sufficiently high values of the Goertler number. However, in most situations of practical interest the basic flow is three dimensional and previous theoretical investigations do not apply. The linear stability of the flow over an infinitely long swept wall of variable curvature is considered. If there is no pressure gradient in the boundary layer the instability problem can always be related to an equivalent two dimensional calculation. However, in general, this is not the case and even for small values of the crossflow velocity field dramatic differences between the two and three dimensional problems emerge. When the size of the crossflow is further increased, the vortices in the neutral location have their axes locally perpendicular to the vortex lines of the basic flow.

Description: An efficient and general method of analyzing a coupled dynamic system of fluid flow and elastic structures is investigated. The improvement of Structural Performance Analysis and Redesign (SPAR) code is summarized. All error codes are documented and the SPAR processor/subroutine cross reference is included.

Description: A three-dimensional solution procedure for the approximate form of the Navier-Stokes equation was exercised in the two- and three-dimensional modes to compute the unsteady turbulent boundary layer on a flat plate corresponding to the data of Karlsson. The procedure is based on the use of a consistently split Linearized Block Implicit technique in conjunction with a QR operator scheme. New time-dependent upstream boundary conditions were developed that yielded realistic solutions for the interior in the vicinity of the upstream boundary. Comparisons of the computation employing these boundary conditions with the data indicate that both qualitative and quantitative agreement was obtained for the mean velocity and the in phase and out of phase components of the first harmonic of the velocity. In addition, the calculation gave results for the skin friction phase angle that had expected physical behavior for large distances downstream of the inflow boundary. For the three-dimensional case, the two-dimensional data of Karlsson was considered, but in a coordinate system skewed at 45 deg to the free stream direction. The results of the calculations were in excellent agreement with the data and the two-dimensional computations.

Description: MacCormack's implicit finite-difference scheme was used to solve the two-dimensional parabolized Navier-Stokes (PNS) equations. This method for solving the PNS equations does not require the inversion of block tridiagonal systems of algebraic equations and permits the original explicit MacCormack scheme to be employed in those regions where implicit treatment is not needed. The advantages and disadvantages of the present adaptation are discussed in relation to those of the conventional Beam-Warming scheme for a flat plate boundary layer test case. Comparisons are made for accuracy, stability, computer time, computer storage, and ease of implementation. The present method was also applied to a second test case of hypersonic laminar flow over a 15% compression corner. The computed results compare favorably with experiment and a numerical solution of the complete Navier-Stokes equations.

Description: A survey of measuring devices or probes, which were used to investigate gas side fouling, was carried out. Five different types of measuring devices are identified and discussed including: heat flux meters, mass accumulation probes, optical devices, deposition probes, and acid condensation probes. A total of 32 different probes are described in detail and summarized in matrix or tabular form. The important considerations of combustion gas characterization and deposit analysis are also given a significant amount of attention. The results show that considerable work was done in the development of gas side fouling probes. However, it is clear that the design, construction, and testing of a durable, versatile probe - capable of monitoring on-line fouling resistances - remains a formidable task.

Description: The low temperature provided by He-3 refrigerators (0.3 to 3 K) have useful space applications. However, the low temperatures and the low surface tension of He-3 require special design considerations. The considerations include the need for small pores to contain the liquid in a matrix; the effects of bubble nucleation and growth; and the effects of the thermal conductivity within the matrix. These design considerations are discussed and a possible confinement system is analyzed.

Description: A closed loop thermoelectric temperature control system is developed for stabilizing sensitive RF integrated circuits within a microwave radiometer to an accuracy of + or - 0.1 C over a range of ambient conditions from -20 C to +45 C. The dual mode (heating and cooling) control concept utilizes partial thermal isolation of the RF units from an instrument deck which is thermally controlled by thermoelectric coolers and thin film heaters. The temperature control concept is simulated with a thermal analyzer program (MITAS) which consists of 37 nodes and 61 conductors. A full scale thermal mockup is tested in the laboratory at temperatures of 0 C, 21 C, and 45 C to confirm the validity of the control concept. A flight radiometer and temperature control system is successfully flight tested on the NASA Skyvan aircraft.

Description: Heat transfer coefficients were measured using both dry and humid air in the same forced convection cooling scheme and were compared using appropriate nondimensional parameters (Nusselt, Prandtl and Reynolds numbers). A forced convection scheme with a complex flow field, two dimensional arrays of circular jets with crossflow, was utilized with humidity ratios (mass ratio of water vapor to air) up to 0.23. The dynamic viscosity, thermal conductivity and specific heat of air, steam and air/steam mixtures are examined. Methods for determining gaseous mixture properties from the properties of their pure components are reviewed as well as methods for determining these properties with good confidence. The need for more experimentally determined property data for humid air is discussed. It is concluded that dimensionless forms of forced convection heat transfer data and empirical correlations based on measurements with dry air may be applied to conditions involving humid air with the same confidence as for the dry air case itself, provided that the thermophysical properties of the humid air mixtures are known with the same confidence as their dry air counterparts.