ALBERT

Description: Transpiration cooling is treated and then full coverage discrete hole injection for three injection orientations. Spacings with pitch to diameter ratios of 5 and 10 are discussed. The array is staggered, with the transverse pitch and the streamwise pitch the same. Results are presented in terms of the Stanton number using the heat transfer coefficient defined in terms of the difference between the wall temperature and the free stream temperature. Two values of Stanton number are provided for each situation: one with the injectant at wall temperature, and the other with the injectant at free stream temperature. These two values are equivalent to knowing the heat transfer coefficient and the adiabatic effectiveness. The heat transfer coefficient thus defined is used with the actual wall temperature to and the actual gas temperature to calculate the heat load. The principle of superposition thus invoked is valid exactly when the governing equations are linear.

Description: Analysis techniques for three aspects of the performance of the NASA/MSFC 32 meter drop tube are considered. Heat loss through the support wire in a pendant drop sample, temperature history of a drop falling through the drop tube when the tube is filled with helium gas at various pressures, and drag and resulting g-levels experienced by a drop falling through the tube when the tube is filled with helium gas at various pressures are addressed. The developed methods apply to systems with sufficiently small Knudsen numbers for which continuum theory may be applied. Sample results are presented, using niobium drops, to indicate the magnitudes of the effects. Helium gas at one atmosphere pressure can approximately double the amount of possible undercooling but it results in an apparent gravity levels of up to 0.1 g.

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: Several environmental parameters presently acknowledged to affect heat transfer are discussed including: (1) the experimental apparatus used, (2) uniform and variable wall temperatures, (3) acceleration effects, (4) deceleration, (5) free stream turbulence, (6) surface roughness, (7) unsteady effects, and (8) secondary flows. Calculation procedures, and some physically based models that are successful in computing heat transfer rates are discussed.

Description: Results are examined from an experiment conducted to determine quantitatively the secondary factors which affect the response of a turbulent boundary layer to convex curvature and to examine the recovery process after curvature ended. The variation of Stanton number with streamwise distance and with enthalpy thickness Reynolds number for the baseline case is shown. The effect of delta sub .99/R on the velocity of the potential core would have if we extended to the wall with no viscous effects, of free stream acceleration, of an unheated starting length, and of boundary layer maturity are discussed. Mixing length and turbulent Prandtl number models are reviewed.

Description: The Lagrangian dispersion theory of Durbin (1980) is used to analyze experiments by Warhaft and Lumley (1978) and by Sreenivasan et al. (1980) on temperature fluctuations in grid-generated turbulence. Both theory and experiment show that the decay exponent m depends on the ratio of the initial length scales of velocity and temperature, although when this ratio is greater than 2.5 such dependence is negligible. The theory shows that m is not truly constant, but within the range covered by the experiments it is nearly so. The agreement between theory and experiment lends credence to the idea that the decay of fluctuations is controlled largely by turbulent relative dispersion.

Description: The effect of large deformation in the flow between the bellmouth and centerbody is considered analytically for application to studies of vortex breakdown in a pipe. Basic equations are defined for axisymmetric inviscid swirling flows at the inflow and outflow sections. Axial and circumferential velocity component profiles are presented, and comparisons are made with trials involving vane angles of 42 deg and Re of 2300. Axial components of the prediction matched well in the inner half of the pipe radius and not well with the outer, while circumferential predictions were good only at the axis. A lack of viscosity was concluded to result in the inaccuracies near the wall.

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: Incompressible turbulent channel flow is investigated by large eddy simulation using improved numerical methods and boundary conditions. In downstream and spanwise directions, cyclic boundary conditions are imposed for velocity and pressure, and two types of boundary conditions near the wall are used and compared. One type is based on the logarithmic law of the mean velocity near the wall and has a slip boundary condition where the molecular-viscous term is neglected. The other type is based on a no-slip boundary condition, where fine mesh spacing near the wall is used to take account of the molecular viscosity. Although the present study employs a coarse mesh (16 x 16 x 21), its results are in good agreement with those of Moin and Kim (1981).

Description: The velocity characteristics of laminar and turbulent developing flow in an S-duct formed from two 22.5-deg bends of rectangular cross-section have been studied experimentally using laser Doppler velocimetry. It is shown that pressure-driven secondary flows arise in the first bend of the duct and reach maxima of 0.22 and 0.15 of the bulk velocity in the laminar and turbulent flows, respectively. The velocities are greater in the laminar flow, mainly because of the thicker inlet boundary layers. On passing through the second half of the S-duct, a secondary flow is established over most of the section in the direction opposite to that in the first half. Near the outer wall of the second bend, however, the secondary flow generated in the first bend is sustained because of the local sign of radial vorticity. This effect contributes to a redistribution of the streamwise isotachs, by the end of the duct, comparable with that in unidirectional bends.

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: The apparent stability of erythrocyte suspensions layered on stationary and flowing Ficoll solutions was studied considering the effects of particle concentration, type and size, and the different flow rates of the particle suspensions and chamber liquid. The data from the flowing system were empirically fitted and, when extrapolated to zero chamber liquid flow rate, gave values comparable to the data from the stationary system, thus confirming the validity of the data and our approach to obtain that data.

Description: A method is developed to determine the shape of steady state solidification interfaces formed when liquid above its freezing point circulates over a cold surface. The solidification interface, which is at uniform temperature, will form in a shape such that the non-uniform energy convected to it is locally balanced by conduction into the solid. The interface shape is of interest relative to the crystal structure formed during solidification; regulating the crystal structure has application in casting naturally strengthened metallic composites. The results also pertain to phase-change energy storage devices, where the solidified configuration and overall heat transfer are needed. The analysis uses a conformal mapping technique to relate the desired interface coordinates to the components of the temperature gradient at the interface. These components are unknown because the interface shape is unknown. A Cauchy integral formulation provides a second relation involving the components, and a simultaneous solution yields the interface shape.

Description: The accuracy of calculations of the radiation emissions from argon plasmas produced by the shock layers over blunt bodies is assessed. The existing theoretical and experimental spectroscopic data on argon are collated. A set of such data is selected for use in the radiative transfer calculations. Calculations are performed for the stagnation regions of the shock layers over laboratory-sized models using these data, and the results are compared with the existing experimental results obtained in a shock-tube. Through this comparison and a parametric study it is shown that radiative heat fluxes at the stagnation point in an argon environment can be calculated within an uncertainty of about 15%. It is shown also that radiative heat fluxes of the order of 100 kW/sq cm can be produced in the existing laboratory facilities.

Description: This paper reports the discovery of a new resonant entrainment phenomenon associated with a confined, pulsed jet flow. It was found that a confined jet, when pulsed at an organ-pipe resonant frequency of the confinement tube, experiences greatly enhanced entrainment and mixing near the exit end of the confinement tube compared to a steady confined jet. The mixing and entrainment rates for the resonantly pulsed confined jet approach, and in some cases slightly exceed, those for an unconfined pulsed jet. Both visual and quantitative evidence of this phenomenon is presented. The new effect should be of considerable interest in ejector and combustor design, both of which benefit from any enhancement in mixing between a primary and a secondary flow

Description: In accordance with the Marangoni effect, immiscible droplets in a host fluid in which a temperature gradient exists move in the direction of increasing temperature. It is proposed that this thermocapillary migration could be used to construct a 'liquid wick' that would return the condensed vapor at the condenser end of a heat pipe back to the evaporator, thus completing the fluid circuit. The droplets would be formed by capillary pressure forcing the condensate through a perforated diaphragm whose temperature would control the droplet flux, and hence the heat flux between the two ends of the heat pipe, thus making it a controllable heat valve.

Description: The shape of a cooled porous wall section is found that will provide a uniform surface temperature, as dictated by material limitations, when the surface is subjected to spatially nonuniform heating. In the analysis, local temperatures and pressures in the porous material are expressed in terms of a potential function. From the imposed thermal conditions, this potential function is governed by the dual constraints of both its value and its normal derivative being specified along the heated surface. The unknown shape of this surface is obtained by meeting these dual conditions. The analytical method uses a generalized conformal mapping procedure that includes a curved boundary. The coolant flow can be compressible or incompressible, and its viscosity can depend on temperature.

Description: A set of nonlinear partial differential equations suitable for the description of a class of turbulent three-dimensional flow fields in select geometries is identified. On the basis of the concept of enforcing a penalty constraint to ensure accurate accounting of ordering effects, a finite element numerical solution algorithm is established for the equation set and the theoretical aspects of accuracy, convergence and stability are identified and quantized. Hypermatrix constructions are used to formulate the reduction of the computational aspects of the theory to practice. The robustness of the algorithm, and the computer program embodiment, have been verified for pertinent flow configurations.

Description: A relatively simple one-dimensional thermal model of the Bridgman growth process has been developed which is applicable to the growth of small diameter samples with conductivities similar to those of metallic alloys. The heat flow in a translating rod is analyzed in a way that is applicable to Biot numbers less than unity. The model accommodates an adiabatic zone, different heat transfer coefficients in the hot and cold zones, and changes in sample material properties associated with phase change. The analysis is applied to several simplified cases. The effect of the rod's motion is studied in a three-zone furnace for a rod sufficiently long that end effects can be neglected; end effects are then investigated for a motionless rod. Finally, the addition of a fourth zone, an independently controlled booster heater between the main heater and the adiabatic zone, is evaluated for its ability to increase the gradient in the sample at the melt interface and to control the position of the interface.

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: Two orthogonal components of velocity and associated Reynolds stresses are determined in a square-sectioned, 90 degree bend of 2.3 radius ratio by utilizing laser-Doppler velocimetry for Reynolds numbers of 790 and 40,000. Results show that boundary layers at the bend inlet of 0.25 and 0.15 of the hydraulic diameter create secondary velocity maxima of 0.6 and 0.4 of the bulk flow velocity, respectively. It is concluded that the boundary layer thickness is important to the flow development, mainly in the first half of the bend, especially when it is reduced to 0.15 of the hydraulic diameter. Smaller secondary velocities are found for turbulent flow in an identical duct with a radius ratio of 7.0 than in the strongly curved bend, although their effect is more important to the streamwise flow development because of the smaller pressure gradients. In addition, the detail and accuracy of the measurements make them suitable for evaluation of numerical techniques and turbulence models.

Description: The model of Warn-Varnas et al. (1978) is used to numerically examine the spin-up flow of a thermally stratified fluid in a cylinder with an insulating side wall, and comparison of the results with the laser-Doppler measurements of Lee (1975) shows excellent agreement. It is shown that flow gradients are created in the interior of the fluid during the meridional circulation spin-up phase, and that the azimuthal flow decayed faster than has been predicted by Wallin (1969). It is established that viscous diffusion in the interior, arising from the interior-flow gradients, is the cause of the discrepancy with Wallin's theory.

Description: A Legendre formulation for a pattern of streamlines adjacent to a surface considered as trajectories with properties consistent with those of a constant vector field is used to develop a mathematical framework for three-dimensional separated flows. Convergence of skin-friction lines onto a particular skin-friction line originating from a particular saddle point is defined as a necessary condition for flow separation. Steady, three-dimensional flow is considered, and singular points occurring in the skin-friction lines are shown to happen where the skin friction or the surface vorticity become zero, and become nodal or saddle points. The separation initiates and continues only globally, as a mixture of an infinite set of friction lines, or locally, with one line. The topography of streamlines in two-dimensional sections of three-dimensional flows is discussed, and examples are provided of a round-nosed body of revolution at varying angles of attack.

Description: This paper describes an experimental program to determine the heat-transfer characteristics of a combustor and heat-exchanger system in a hybrid solar receiver which utilizes a Stirling engine. The system consists of a swirl combustor with a crossflow heat exchanger composed of a single row of 48 closely spaced curved tubes. In the present study, heat-transfer characteristics of the combustor/heat-exchanger system without a Stirling engine have been studied over a range of operating conditions and output levels using water as the working fluid. Nondimensional heat-transfer coefficients based on total heat transfer have been obtained and are compared with available literature data. The results show significantly enhanced heat transfer for the present geometry and test conditions. Also, heat transfer along the length of the tubes is found to vary, the effect depending upon test condition.

Description: Asymptotic and numerical techniques in bifurcation theory are applied to the Young-Laplace equation governing meniscus shape in order to analyze the dependence of the shape and stability of rigidly rotating drops held captive between corotating solid faces on multiple parameters. Asymptotic analysis of the evolution of drop shape from the cylindrical as a function of distance between the solid faces, drop volume, rotational Bond number and gravitational Bond number shows that some shape bifurcations from cylinders to wavy, axisymmetric menisci are ruptured by small changes in drop volume or gravity. Computer calculations of axisymmetric drop shapes based on a finite element representation of the interface and numerical algorithms for tracking shape families and singular points are then used to map drop stability for the four-dimensional parameter space. The results of the asymptotic and numerical analyses are shown to agree well within the limited range of parameters where the asymptotic analysis is valid.

Description: An alternating direction implicit (ADI) method has been applied to a staggered grid for the computation of convection in a highly stratified fluid. Since artificial viscosity is not needed, subtle effects like the onset of convection can be studied. These computations compare well with the 2-D results by Graham (1975) and also agree with standard Boussinesq results when taken to that limit. Good efficiency has been achieved with a time step hundreds of times larger than the stability limit imposed by the explicit treatment of diffusion and the Courant number is not restricted to be below 1. The Navier-Stokes equation contains cross spatial derivatives which are treated explicitly in most ADI schemes. The destabilizing effect of such a practice on a 2-D model system with second-order spatial derivative terms only was analyzed and found to be not excessive.

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 shock resolution of Harten's (1982) second-order explicit method for one-dimensional hyperbolic conservation laws is investigated for a two-dimensional gas-dynamic problem. The possible extension to a high resolution implicit method for both one- and two-dimensional problems is also investigated. Applications of Harten's method to the quasi-one-dimensional nozzle problem with two nozzle shapes (divergent and convergent-divergent) and the two-dimensional shock-reflection problem resulted in high shock resolution steady-state numerical solutions.

Description: Incompressible viscous flow fields induced by initial vorticity distributions with bounded support or exponential decay in the far field are investigated. A numerical scheme for the solution of the vorticity distribution and the velocity field is presented with special emphasis on the treatment of the boundary data. The efficiency of the scheme is demonstrated. The present method has been applied to the study of the merging and collision of vortex rings.

Description: The theoretical basis for well posed marching of a Parabolic Navier-Stokes (PNS) computational technique for supersonic flow is discussed and examples given to verify the analysis. It is demonstrated that stable computations can be made even with very small steps in the marching direction. The method is applied to cones at large angle of attack in high Reynolds number, supersonic flow. Streamline trajectories generated from the numerical solutions demonstrate the development of vortex structures of the lee side of the cone. Previously announced in STAR as N83-22551

Description: A numerical procedure for the relaxation solution of the full steady Euler equations is described. By embedding the Euler system in a second order surrogate system, central differencing may be used in subsonic regions while retaining matrix forms well suited to iterative solution procedures and convergence acceleration techniques. Hence, this method allows the development of stable, fully conservative differencing schemes for the solution of quite general inviscid flow problems. Results are presented for both subcritical and shocked supercritical internal flows. Comparisons are made with a standard time dependent solution algorithm. Previously announced in STAR as N82-24859

Description: The physical features of steady and unsteady freestream separating turbulent boundary layers that have been determined by pointwise laser anemometer measurements are outlined. It is seen that the large-scale structures control the outer region's backflow behavior. Near the wall, the mean backflow velocity profile for both the steady and unsteady cases is found to scale on the maximum negative mean velocity and its distance from the wall. A description is given of a scanning laser anemometer that produces nearly instantaneous velocity profiles for examing the temporal features of these large-scale structures. Also described is a 'zero-wake' seeder that supplies particles to the outer shear layer and freestream flow with a minimal disturbance.

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: A new numerical technique for simulating three dimensional, unsteady, incompressible pipe flows is presented and its utility and accuracy is shown. Each vector function in the expansion of the velocity field is divergence free and satisfies the boundary conditions for viscous flow. Some of the benefits of the expansion technique are that pressure is eliminated from the dynamics, only two unknowns per mesh point are required, implicit treatment of the viscous terms is provided at no extra computational cost, and no fractional time steps are required. The method uses spectral expansions: Fourier series in the azimuthal and streamwise directions, and Jacobi polynominals in the radial direction. Previously announced in STAR as N82-31644

Description: Developments in three dimensional, time dependent numerical simulation of turbulent flows bounded by a wall are reviewed. Both direct and large eddy simulation techniques are considered within the same computational framework. The computational spatial grid requirements as dictated by the known structure of turbulent boundary layers are presented. The numerical methods currently in use are reviewed and some of the features of these algorithms, including spatial differencing and accuracy, time advancement, and data management are discussed. A selection of the results of the recent calculations of turbulent channel flow, including the effects of system rotation and transpiration on the flow are included. Previously announced in STAR as N82-28577

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 use of superminicomputers for solving a series of increasingly complex thermal analysis problems is investigated. The approach involved (1) installation and verification of the SPAR thermal analyzer software on superminicomputers at Langley Research Center and Goddard Space Flight Center, (2) solution of six increasingly complex thermal problems on this equipment, and (3) comparison of solution (accuracy, CPU time, turnaround time, and cost) with solutions on large mainframe computers.

Description: Thermocapillary stability characteristics of a horizontal liquid layer heated from below rotating about a vertical axis and subjected to a uniform vertical magnetic field are analyzed under a variety of thermal and electromagnetic boundary conditions. Results based on analytical solutions to the pertinent eigenvalue problems are discussed in the light of earlier work on special cases of the more general problem considered here to show in particular the effects of the heat transfer, nonzero curvature and gravity waves at the two-fluid interface. Although the expected stabilizing action of the Coriolis and Lorentz force fields in this configuration are in evidence the optimal choice of an appropriate range for the relevant parameters is shown to be critically dependent on the interfacial effects mentioned above.

Description: Present knowledge of the mechanisms for production and enrichment and film drops by bursting bubbles is summarized, with particular emphasis on the unsolved problems. Sea salt is by far the major constituent cycled through the Earth's atmosphere each year. Bursting bubbles in the oceans appear to be primarily responsible. These salt particles play a role in the formation of maritime clouds, which in turn affect the Earth's radiation budget. Along with the salt are carried various chemical pollutants and potentially pathogenic microorganisms, often in highly enriched form.

Description: A model of the change in shape of a raindrop is presented. Raindrops measured by two orthogonal cameras were classified by shape and orientation to determine the nature of the oscillation. A physical model based on potential energy was then developed to study the amplitude variation of oscillating drops. The model results show that oscillations occur about the equilibrium axis ratio, but the time average axis ratio if significantly more spherical for large amplitudes because of asymmetry in the surface potential energy. A generalization of the model to oscillations produced by turbulence yields average axis ratios that are consistent with the camera measurements. The model results for average axis ratios were applied to rainfall studies with a dual polarized radar.

Description: Results of a wind tunnel experiment in which electrically uncharged water drops of 500 to 3000 microns equivalent radius are freely suspended in the vertical air stream of the UCLA cloud tunnel are presented. During this suspension the drops were exposed to external vertical electric fields of 500 to 8,000 volts/cm. The change in drop shape with drop size and electric field strength was noted and is discussed in the light of theoretical work cited in the literature which unfortunately does not take into account the effects of air flow past the drop. The wind tunnel study is documented by stills from a 16 mm film record that demonstrates the shape of water drops in response to both hydrodynamic and electric forces.

Description: The shape of a gas bubble which rises through a quiescent incompressible, Newtonian fluid at intermediate Reynolds numbers is considered. Exact numerical solutions for the velocity and pressure fields, as well as the bubble shape, are obtained using finite difference techniques and a numerically generated transformation to an orthogonal, boundary-fitted coordinate system. No restriction is placed on the allowable magnitude of deformation.

Description: Holographic studies were performed which examined the fragmentation process during vapor explosion of a water-in-fuel (hexadecane/water) emulsion droplet. Holograms were taken at 700 to 1000 microseconds after the vapor explosion. Photographs of the reconstructed holograms reveal a wide range of fragment droplet sizes created during the explosion process. Fragment droplet diameters range from below 10 microns to over 100 microns. It is estimated that between ten thousand and a million fragment droplets can result from this extremely violent vapor explosion process. This enhanced atomization is thus expected to have a pronounced effect on vaporization processes which are present during combustion of emulsified fuels.

Description: The configuration of liquid hydrogen inside spherical glass shell ICF target was studied both theoretically and experimentally. Because of the zero contact angle between the .D2 liquid and glass substrate and the limited wetting surface that is continuous, the liquid hydrogen completely covers the interior of the glass shell, resulting in the formation of a void at the center. For this reason, the present problem distinguishes itself from that for a sessile drop sitting on a flat surface. A theory was formulated to calculate the liquid hydrogen configuration by including the London-dispersion force between the liquid and the substrate molecules. The net result is an augmented Bashforth-Adams equation appropriate to a spherical substrate, which is considered to be the major contribution of the present work. Preliminary calculations indicate that this equation accurately models the liquid hydrogen behavior inside a spherical microshell.

Description: The NASTRAN thermal analyzer (NTA) which performs large-scale unified thermo-structural analyses with the NASTRAN (NASA structural analysis) computer program is described. The mathematical similitude between these two distinct disciplines of thermal and structure is examined. It serves as the theoretical basis upon which the implementation of the thermal capability in NASTRAN was accomplished. The program structure, the functional flow, the solution algorithms, the organization of an input data deck and the solution capabilities of NTA are summarized. Emphasis is placed on the interface of the unified approach in thermo-structural analyses where stresses, deflections, vibrations and bucklings induced by the effect of temperature change are of concern. Attentions are also directed to the preprocessor and post processors. As a specially designed preprocessor, the VIEW program is capable of generating exchange factors which can be output, at user's option, in formats compatible with that required by NTA. Two post processors that serve specific objectives are included. They are the thermal variance analysis and the graphical displaying capability of temperatures in color or black and white.

Description: Axisymmetric equilibrium shapes of conducting drops and bubbles, (1) pendant or sessile on one face of a circular parallel-plate capacitor or (2) free and surface-charged, are found by solving simultaneously the free boundary problem consisting of the augmented Young-Laplace equation for surface shape and the Laplace equation for electrostatic field, given the surface potential. The problem is nonlinear and the method is a finite element algorithm employing Newton iteration, a modified frontal solver, and triangular as well as quadrilateral tessellations of the domain exterior to the drop in order to facilitate refined analysis of sharply curved drop tips seen in experiments. The stability limit predicted by this computer-aided theoretical analysis agrees well with experiments.

Description: The asymptotic properties for the small Bond number B of the equilibrium capillary interface interior to a circular cylindrical tube vertically dipped in an infinite reservoir of liquid are discussed. (The Bond number B is a dimensionless parameter which is the ratio of gravitational to capillary forces.) The formal expansion in powers of B of the solution to the differential equation describing the equilibrium surface (as can be obtained by standard perturbation methods) is proved to be truly asymptotic -- to all orders and uniformly in the variable and parameter gamma, the contact angle. Sequences of general estimates, in closed form, from above and from below, are also given for the solution and related functions. The M-th term in these sequences are asymptotically exact to order m. An idiosyncrasy of the problem, crucial in obtaining these estimates, is the absolute monotonicity of the structural function of the system in integral form.

Description: On the basis of both a conventional relativistic nuclear fluid dynamic model and a two fluid generalization that takes into account the interpenetration of the target and projectile upon contact, collisions between heavy nuclei moving at relativistic speeds are calculated. This is done by solving the relevant equations of motion numerically in three spatial dimensions by use of particle in cell finite difference computing techniques. The effect of incorporating a density isomer, or quasistable state, in the nuclear equation of state at three times normal nuclear density, and the effect of doubling the nuclear compressibility coefficient are studied. For the reaction 20Ne + 238U at a laboratory bombarding energy per nucleon of 393 MeV, the calculated distributions in energy and angle of outgoing charged particles are compared with recent experimental data both integrated over all impact parameters and for nearly central collisions.

Description: A particular configuration of a vertical capillary tube for which S is the equilibrium interface between two fluids in the presence of a downward pointing gravitational field was investigated. S is the graph a function u whose domain is the (horizontal) cross section gamma of the tube. The mean curvature of S is proportional to its height above a fixed reference plane and lambda is a prescribed constant and may be taken between zero and pi/2. Domains gamma for which us is a bounded function but does not extend continuously to d gamma are sought. Simple domains are found and the behavior of u in those domains is studied. An important comparison principle that has been used in the literature to derive many of the results in capillarity is reviewed. It allows one to deduce the approximate shape of a capillary surface by constructing comparison surfaces with mean curvature and contact angle close to those of the (unknown) solution surface. In the context of nonparametric problems the comparison principle leads to height estimates above and below for the function u. An example from the literature where these height estimates have been used successfully is described. The promised domains for which the bounded u does not extend continuously to the boundary are constructed. The point on the boundary at which u has a jump discontinuity will be the vertext of a re-entrant corner having any interior angle theta pi. Using the comparison principle the behavior of u near this point is studied.

Description: The possibility of rebound for colliding cloud drops was measured by determining the collection efficiency. The collection efficiency for 17 size pairs of relatively uncharged drops in over 500 experimental runs was measured using two techniques. The collection efficiencies fall in a narrow range of 0.60 to 0.70 even though the collection drop was varied between 63 and 326 microns and the size ratio from 0.05 to 0.33. In addition the measured values of collection efficiencies (Epsilon) were below the computed values of collision efficiencies (E) for rigid spheres. Therefore it was concluded that rebound was occurring for these sizes since inferred coalescence (epsilon = Epsilon/E) efficiencies are about 0.6 yo 0.8. At a very small size ratio (r/R = p = 0.05, R = 326 microns) the coalescence efficiency inferred is in good agreement with the experimental findings for a supported collector drop. At somewhat large size ratios the inferred values of epsilon are well above results of supported drop experiments, but show a slight correspondence in collected drop size dependency to two models of drop rebound. At a large size ratio (p = 0.73, R = 275) the inferred coalescence efficiency is significantly different from all previous results.

Description: Acoustic levitation and the response of fluid spheres to spherical harmonic projections of the radiation pressure are described. Simplified discussions of the projections are given. A relationship between the tangential radiation stress and the Konstantinov effect is introduced and fundamental streaming patterns for drops are predicted. Experiments on the forced shape oscillation of drops are described and photographs of drop fission are displayed. Photographs of critical angle and glory scattering by bubbles and rainbow scattering by drops are displayed.

Description: A series of experimental tests was carried out on an 'OGEE' shaped planform, liquid air-shear electrostatic nozzle. Liquid was ejected from the upper surface of the nozzle and was then dispersed and atomized efficiently by a high speed air flow passing over the nozzle and by the effect of two very strong coherent air vortices generated by the 'OGEE' shaped nozzle surface. Initial test results which are presented show the nozzle to perform far superior to a similar delta wing shaped design which is used extensively in various industries applications.

Description: Experiments are presented which were conducted on flow fields produced by a circulation control airfoil utilizing the Coanda effect at the trailing edge. The application of holographic interferometry to obtain both visualization and quantitative data on the flow field about a circulation control airfoil at transonic flow speed is covered. A brief description of the flow model and measurement techniques is given. The data reduction procedure, results, and interpretation are presented. The results have provided a good deal of information on the character of the flow field, particularly in the neighborhood of the trailing edge. As to the airfoil design, it is apparent that improved performance can be achieved if jet detachment is delayed. Another design improvement would involve the development of an optimum trailing-edge shape for the expected operating Mach and Reynolds number ranges.

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.