ALBERT

Description: Heat pipes may be employed as temperature regulators, heat diodes, transformers, storage batteries, or utilized for transforming thermal energy into mechanical, electric, or other forms of energy. General concepts were established for the analysis of the transfer process in heat pipes. A system of equations was developed to describe the thermodynamics of steam passage through a cross section of a heat pipe.

Description: A description is presented of the split finite-volume method which is a viable numerical procedure for performing with the aid of a modern special purpose vector computer numerical simulation studies of complicated flow fields, including chemical reactions, about geometrically complex bodies. Such numerical studies are needed for the development of atmospheric entry vehicles such as the space shuttle. The equations which are approximated are quite general and can be used in studies of combustion, pollution, and other chemically reacting flow phenomena, where convective transport effects dominate the influence of radiative, viscous, and other transport mechanisms. The shock perturbed flow about a shuttle orbiter flying at a large angle of attack during atmospheric entry is illustrated. The method uses a time splitting of the convection differencing operator to achieve efficient data management.

Description: The possibilities of producing heat pipes and, especially, the necessary capillary structures are discussed. Several types of heat pipes were made from stainless steel and tested at temperatures between 400 and 1055 deg C. The thermal power was determined by a calorimeter. Results indicate: bubble-free evaporation of sodium from rectangular open chennels is possible with a heat flux of more than 1,940 W/sq cm at 1055 C. The temperature drop along the tube could be measured only at low temperatures. A subdivided heat pipe worked against the gravitational field. A heat pipe with a capillary structure made of a rolled screen was supported by rings and bars operated at 250 W/sq cm heat flux in the evaporating region.

Description: The reflection factor at a tube which ends at a plate over which a flow is forming was determined as a function of the Strouhal number, formed from the flow velocity, the aperture radius, and the acoustic frequency. Several adjacent openings were investigated to determine the interactions between several openings.

Description: The heat transfer performance of baffled cooling fins on cylinder heads of small, air-cooled, general-aviation aircraft engines was analyzed to determine the potential for improving cooling fin design. Flow baffles were assumed to be installed tightly against the fin end edges, an ideal baffle configuration for guiding all flow between the fins. A rectangular flow passage is thereby formed between each set of two adjacent fins, the fin base surface, and the baffle. These passages extend around each side of the cylinder head, and the cooling air absorbs heat as it flows within them. For each flow passage length, the analysis was concerned with optimizing fin spacing and thickness to achieve the best heat transfer for each fin width. Previous literature has been concerned mainly with maximizing the local fin conductance and has not considered the heating of the gas in the flow direction, which leads to higher wall temperatures at the fin passage exits. If the fins are close together, there is a large surface area, but the airflow is restricted.

Description: The objective of this investigation was to study the effects of small amplitude sinusoidal pulsations on fully developed turbulent flow in a tube from both experimental and theoretical viewpoints. Theoretical models for the macroscopic behavior of pulsating turbulent tube flow were developed for the two cases of very low and very high pulsation frequencies. The models are based on assumptions of quasi-steady and frozen eddy viscosity flow behavior, respectively. The models successfully predict unsteady velocity profiles, thereby supporting the currently proposed definitions of frequency regimes in pulsating turbulent flow. Experimental measurements were made of the time-dependent pressure drop and velocity profiles over the range of frequency-to-Reynolds number ratios from 0.0095 to 0.24. The two macroscopic models developed in this study predict unsteady velocity profiles which are in moderately good agreement with the experiments in their respective frequency regimes, and a previously developed quasi-steady model is found to predict experimental velocity profiles well in both the quasisteady and the frozen eddy viscosity frequency regimes. The effect of flow pulsations on the dissipation of turbulence energy in the vicinity of the wall was measured in the lower transition frequency regime. The long-time averaged dissipation was observed to be unchanged from the steady flow dissipation, within the accuracy of the experiment. A theoretical model of the periodic viscous sublayer was also developed and applied to pulsating flow in a tube, in order to investigate the effects of flow pulsations on the rate of production of turbulence in the region of the wall. The periodic viscous sublayer model predicts sublayer growth periods in steady flow which agree with the published experimental data. When the model is applied to pulsating flow, the response of the sublayer growth period falls into three frequency regimes, the parameters of which are in approximate agreement with the frequency regimes which are defined on the basis of macroscopic flow behavior. The sublayer renewal cycle exhibits quasi-steady flow behavior when the sublayer growth period is much less than the pulsation period, transition behavior when these two periods are approximately equal, and frozen eddy viscosity behavior when the sublayer period is much longer than the pulsation period. The effect of the sublayer growth and renewal cycle on the level of turbulence was investigated by two methods. The velocity fluctuations seen by a point velocity probe located close to the wall were predicted from the model in one method and the rate of turbulence production was estimated from the frequency of sublayer renewal events in the other.

Description: This paper describes the development of a probe to study turbulent reactive mixing of two hypergolic streams of gas. The streams contain CO and O2 respectively and are mixed at 1400 K. The probe is a water cooled sonic orifice probe followed by a Fanno duct. Tests showed that it rapidly quenched the reaction, and that time-average, ingested gas composition accurately reflected the time-average composition at the sampling point. A new way to measure mixedness to the molecular level in a reactive flow situation is described, based on time-average compositions at the probe. This definition is compared to the definition of intensity of segregation which can be used for unreactive mixing situations. Finally, the probe is used to examine centerline composition in a confined coaxial mixing situation.

Description: Radiation intensity profiles behind incident normal shock waves in pure CO2 have been measured spectroscopically in the Langley Arc-Driven Shock Tube. These profiles, which were obtained for shock velocities between 9 and 13 km/sec and ambient densities corresponding to Venus altitudes between 100 and 80 km, were measured in the vacuum ultraviolet regime. Wavelengths of 127.7, 158.0, 177.5, and 195.0 nm were monitored simultaneously using a four-channel vacuum spectrograph equipped with sodium-salicylate-coated photomultipliers, thereby including the CO(4+) band system which is the most prominent radiator. Measured nonequilibrium overshoots are modeled to provide a means of estimating the effect of nonequilibrium radiation heating to the stagnation region of proposed aero-shells for Venusian entry. These results indicate a significant increase in radiative heating due to nonequilibrium effects. The measurements are believed to represent the most accurate data available on the effect of nonequilibrium radiative heat transfer for Venus entry. This accuracy is primarily due to improved spectrographic instrumentation, which is discussed in some detail regarding its application in related studies.

Description: The effects of vibration, flow transients, and warm gas pressurization on capillary acquisition system performance were evaluated. The degradation observed under wide band random and high frequency sinusoidal vibration was of a substantially different nature from that obtained under low frequency sinusoidal vibration. With the former, ingestion of small gas bubbles into the liquid region was correlated by a hydrostatic model, while the capillary stability was destroyed and liquid was lost from the liquid region with the latter. No degradation was observed as a result of flow transients in a flight-type multichannel screen device, but it was observed in a transparent laboratory device. Liquid hydrogen outflow tests were conducted with a multilayer dual-screen-liner system with both helium and hydrogen pressurant gases. The tendency towards dryout of the device with hydrogen pressurant was found to increase with increasing pressurant temperature and length of prepressurization period. Dryout did not occur with helium pressurant.

Description: A mean turbulent closure method with a length scale function is developed for predicting the reduced spreading rate in high Mach number, high Reynolds number free shear layers. The new feature of the method is the inclusion of effects of the pressure-dilation term which appears in the compressible form of the turbulence energy equation. The model indicates that the pressure-dilation correlation may act as a turbulence energy sink or source depending on whether the transverse gradients of mean Mach number and mean density have the same or opposite signs.

Description: Mean velocity and stagnation pressures of a coaxial air jet injected into an ambient atmosphere were measured and local free-stream static pressure profiles were derived from these measurements. The mean velocity profiles show a dramatic decrease in the velocity of the central jet immediately downstream of the injection station. It is found that the high velocity outer jet creates a reduction in static pressure through an aspiration effect, which causes the central jet to spread, thereby reducing its velocity. Results demonstrate that even for the simple coaxial injector geometry often used in H2/O2 rockets, conventional modeling techniques do not apply.

Description: Unmixedness in flames refers to the condition in which eddies of fuel and oxidant have engulfed each other but where the molecular mixing necessary for complete reaction has not taken place. This phenomenon is studied experimentally using a flow geometry consisting of a coaxial Mach 2 hydrogen jet surrounded by a parallel cold Mach 2 air stream. Gas sample measurements taken downstream in the region near the flame revealed the presence of unreacted hydrogen and oxygen in the same sample. The data were analyzed by means of the CHARNAL program which uses a step function variation with time for instantaneous concentration.

Description: A tentative scheme is proposed for constructing second-order closure models of turbulent reacting flows, with particular reference to mixing processes in turbulent diffusion flows. The model is based on the concepts of typical structure and typical eddy, and treats various hypothetical turbulent flame situations, including the time dependence of species mass fractions passing a point in the flame, the most typical occurrence of species mass fractions, and the case of a rather bizarre typical eddy which shows the effects of nonmixedness.

Description: LOX and LN2 data for two-phase critical flow through nozzles have been acquired with precision control. The principal measured parameters were inlet conditions, critical flow rate and critical flow pressure ratio. The data conclusively demonstrate that the principle of corresponding states can be applied to two-phase choked flow through nozzles. These data also demonstrate that the proper normalizing parameters have been developed and current theories can provide an adequate means for extrapolating to other fluids.

Description: Low speed wind tunnel tests were conducted to determine the aeroacoustic performance of several model sonic inlets. The results were analyzed to indicate how inlet aeroacoustic characteristics were affected by inlet design and operating conditions. A system for regulating sonic inlet noise reduction was developed and tested. Results indicate that pressure losses at forward velocity may be substantially less than those at static conditions. This is particularly true for translating centerbody inlets with the centerbody extended in the approach and landing position. Operation to simulated takeoff incidence angles of 50 degrees was demonstrated with good inlet performance. Inlet sound pressure level reduction was regulation was regulated to within approximately + or - 1 dB by controlling inlet surface static pressure measured at the diffuser exit.

Description: Choked flow rates and axial pressure distributions were measured for subcooled nitrogen in a converging-diverging nozzle with a constant area section in the throat region. Stagnation pressures ranged from slightly above saturation to twice the thermodynamic critical pressure. Stagnation temperatures ranged from 0.75 to 1.03 times the thermodynamic critical temperature. The choking plane is at the divergence end of the constant area throat section. At high stagnation pressures the fluid stays liquid well into the constant area throat region; at near saturation stagnation pressures it appears that vaporization occurs at or before the entrance to the constant area throat region. The throat-to-stagnation pressure ratio data exhibits an anomalous flat region, and this anomaly is related to the two-phase process. The fluid is metastably all liquid below the saturation pressure.

Description: Approaches for developing an analytical model capable of determining the effects of rotor flow and blade parameters and turbulence properties (i.e. energy, velocity correlations, and length scale) on the rotor wake characteristics and its diffusion properties are discussed. The three-dimensional model will employ experimental measurements, instantaneous velocities, and turbulence properties at various stations downstream from a rotor. A triaxial probe and a rotating conventional probe, which is mounted on a traverse gear operated by two step motors, are to be used for these measurements. The final rotor wake model will be capable of predicting the discrete and broadband noise generated in a fan rotor and of evaluating the aerodynamic losses, efficiency and optimum spacing between a rotor and stator in turbomachinery.

Description: Axial and swirling airflows were used to break up water jets and sheets into sprays of droplets to determine the overall effects of orifice diameter, weight flow of air, and the use of an air swirler on fineness of atomization as characterized by mean drop size. A scanning radiometer was used to determine the mean drop diameter of each spray. Swirling airflows were produced with an axial combustor, 70 deg blake angle, air swirling. Water jets were injected axially upstream, axially downstream and cross stream into the airflow. In addition, pressure atomizing fuel nozzles which produced a sheet and ligament type of breakup were investigated. Increasing the weight flow rate of air or the use of an air swirling markedly reduced the spray mean drop size.

Description: The development of the pressure spectrum has been followed from its separation into turbulent-turbulent and turbulent-mean square contributions to a unified spectral model. The model is obtained by Fourier transforming the integral solution to Poisson's equation for an isotropic, homogeneous, constant mean shear flow. The variations in the asymptotic form for the mean-shear spectrum is discussed. It is noted that the spectra are associated with the particular velocity field selected, and that by an appropriate choice of spectral forms, the mean-square pressure of the flow may be calculated.

Description: Paper describes a numerical calculation scheme for tangential slot injection (wall-wake) flows; application of the scheme over a wide range of flow conditions indicates increased accuracy compared to previous work. Predictions from the numerical code were in good agreement with experiment (velocity profile, skin-friction, and effectiveness data) for low- and high-speed flows. To achieve improved accuracy, modifications in the turbulence modeling, compared to previous research, were necessary for the imbedded shear layer region in the near field and for the wall region near shear layer impingement. Anomalous behavior was noted for far field experimental velocity profiles in low-speed flow when the slot-to-free stream velocity ratio was near one

Description: This paper is concerned with the two-dimensional supersonic flow of a thick turbulent boundary layer over a train of relatively small wave-like protuberances. The flow conditions and the geometry are such that there exists a strong interaction between the viscous and inviscid flow. The problem cannot be solved without inclusion of interaction effects due to the occurrence of the separation singularity in classical boundary layer methods. Here the interacting boundary layer equations are solved numerically using a time-like relaxation method with turbulence effects represented by the inclusion of the eddy viscosity model of Cebeci and Smith. Results are presented for flow over a train of up to six waves for Mach numbers of 2.5 and 3.5, Reynolds numbers of 10,000,000/m and 32,000,000/m, and wall temperature ratios of 0.4 and 0.8. Limited comparisons with independent experimental and analytical results are also given.

Description: Data for the energy transfer function are used to estimate the degree of localness of energy transfer in homogeneous turbulence. It is found that in regions where the energy which enters a wavenumber band is greater than the energy leaving, much of the energy entering the band is produced by wavenumbers an order of magnitude smaller. Thus for both low and high Reynolds numbers, spectral energy transfer is nonlocal. The tendency of the energy to jump between separated wavenumber regions agrees with the theory that turbulence forms concentrated regions of large velocity gradients. It is also felt that the universal equilibrium theory may be applicable if the Reynolds number of the turbulence is very high.

Description: Forced convection and subcooled boiling heat transfer data for liquid nitrogen and liquid neon were obtained in support of a design study for a 30 tesla cryomagnet cooled by forced convection of liquid neon. This design precludes nucleate boiling in the flow channels as they are too small to handle vapor flow. Consequently, it was necessary to determine boiling incipience under the operating conditions of the magnet system. The cryogen data obtained over a range of system pressures, fluid flow rates, and applied heat fluxes were used to develop correlations for predicting boiling incipience and convective boiling heat transfer coefficients in uniformly heated flow channels. The accuracy of the correlating equations was then evaluated. A technique was also developed to calculate the position of boiling incipience in a uniformly heated flow channel. Comparisons made with the experimental data showed a prediction accuracy of plus or minus 15 percent

Description: Two phase and gaseous choked flow data for fluid nitrogen were obtained for a test section which was a long constant area duct of 16 200 L/D with a diverging diffuser attached to the exit. Flow rate data were taken along five isotherms (reduced temperature of 0.81, 0.96, 1.06, 1.12, and 2.34) for reduced pressures to 3. The flow rate data were mapped in the usual manner using stagnation conditions at the inlet mixing chamber upstream of the entrance length. The results are predictable by a two-phase homogeneous equilibrium choking flow model which includes wall fraction. A simplified theory which in essence decouples the long tube region from the high acceleration choking region also appears to predict the data reasonably well, but about 15 percent low.

Description: Previously reported vaporization time data of liquid nitrogen drops in film boiling on a flat plate are about 30 percent shorter than predicted from standard laminar film boiling theory. This theory, however, had been found to successfully correlate the data for conventional fluids such as water, ethanol, benzene, or carbon tetrachloride. This paper presents experimental evidence that some of the discrepancy for cryogenic fluids results from ice contamination due to condensation. The data indicate a fairly linear decrease in droplet evaporation time with the diameter of the ice crystal residue. After correcting the raw data for ice contamination along with convection, a comparison of theory with experiment shows good agreement.

Description: A scanning radiometer was used to determine the effect of airstream velocity on the mean drop diameter of water sprays produced by pressure atomizing and air atomizing fuel nozzles used in previous combustion studies. Increasing airstream velocity from 23 to 53.4 meters per second reduced the Sauter mean diameter by approximately 50 percent with both types of fuel nozzles. The use of a sonic cup attached to the tip of an air assist nozzle reduced the Sauter mean diameter by approximately 40 percent. Test conditions included airstream velocities of 23 to 53.4 meters per second at 293 K and atmospheric pressure.

Description: In the state of Leidenfrost boiling, liquid drops are observed to vibrate in a variety of modal patterns. Theories are presented which predict the frequency of oscillation and show that the observed model patterns of drops correspond to the minimum energy oscillatory excitation state. High-speed photographic techniques were used to record these motions and substantiate the theories. An incipient temperature was also found for water drops in film boiling below which free oscillations do not exist. In addition to these oscillations, photographic sequences are presented which show that wave motion can exist along the circumference of the drop. Following the study of free oscillations, the system was mounted on a shaker table and the drop subjected to a range of forced frequencies and accelerations.

Description: During a test program to investigate low-cycle thermal fatigue, 21 of 22 cylindrical test sections of a cylindrical rocket thrust chamber were thermally cycled to failure. Cylinder liners were fabricated from OFHC copper, Amzirc, and NARloy-Z. The cylinders were fabricated by milling cooling channels into the liner and closing out the backside with electrodeposited copper. The tests were conducted at a chamber pressure of 4.14 MN/sq m (600 psia) and an oxidant-fuel ratio of 6.0 using hydrogen-oxygen as propellants. The average throat heat flux was 54 MW/sq m (33 Btu/sq in./sec). All of the failures were characterized by a thinning of the cooling channel wall and eventual failure by tensile rupture. The 1/2-hard Amzirc material showed little improvement in cyclic life when compared with OFHC copper; while the NARloy-Z and aged Amzirc materials had the best cyclic life characteristics. One OFHC copper cylinder was thermall cycled 2044 times at a steady-state hot-gas-side wall temperature of 514 K (925 R) without failing.

Description: The sensitivity of silica heat shield requirements to gap width, tile edge radius, and heat transfer distribution within tile gaps was investigated. A two-dimensional thermal model was modified and used to determine the effect of two dimensional heat transfer distributions at high temperature reusable surface insulation edges on shuttle thermal protection system (TPS) requirements. The sensitivity of TPS requirements to coating thickness, emissivity, substructure thickness, and changes in gap heating for several locations on shuttle was also studied. An inverse solution technique was applied to temperature data obtained in the Ames 20 MW turbulent duct in order to examine the effect of tile edge radius on TPS requirements. The derived heating values were then used to predict TPS requirements. Results show that increasing tile radius reduces TPS requirements.

Description: An analytical study of nonsimilar jet mixing is made for compressible, nonisoenergetic flows. The conservation equations are solved for each of the streams above and below the dividing streamline by using Meksyn's asymptotic method of integration for solving boundary-layer problems. The problem of laminar mixing between two parallel streams is investigated for the case of a constant pressure gradient. It is found that the velocity and temperature profiles from the exact solution to the nonsimilar governing equations can be well approximated by the locally similar solution.

Description: The unsteady, three-dimensional flowfield resulting from the interaction of a plane shock with a cone-shaped vehicle traveling supersonically is determined, using a second-order, shock-capturing, finite-difference approach. The time-dependent, inviscid gasdynamic equations are transformed to include the self-similar property of the flow, to align various coordinate surfaces with known shock waves, and to cluster points in the vicinity of the intersection of the transmitted incident shock and the surface of the vehicle. The governing partial differential equations in conservation-law form are then solved iteratively using MacCormack's (1969) algorithm.

Description: A recuperator consisting of two fluidized bed regenerators with circulating solid particles is considered for use in a Brayton cycle. These fluidized beds offer the possibility of high temperature operation if ceramic particles are used. Calculations of the efficiency and size of fluidized bed regenerators for typical values of operating parameters have been made and compared to a shell and tube recuperator. The calculations indicate that the fluidized beds will be more compact than the shell and tube as well as offering a high temperature operating capability.

Description: A study was performed to determine the desirability of replacing the hydrogen peroxide settling system on the Centaur D-1S with a capillary acquisition system. A comprehensive screening was performed to select the most promising capillary device fluid acquisition, thermal conditioning, and fabrication techniques. Refillable start baskets and bypass feed start tanks were selected for detailed design. Critical analysis areas were settling and refilling, start sequence development with an initially dry boost pump, and cooling the fluid delivered to the boost pump in order to provide necessary net position suction head (NPSH). Design drawings were prepared for the start basket and start tank concepts for both LO2 and LH2 tanks. System comparisons indicated that the start baskets using wicking for thermal conditioning, and thermal subcooling for boost pump NPSH, are the most desirable systems for future development.

Description: Data are presented of an experiment in which subcooled liquid nitrogen was discharged through a sharp-edged orifice at flow rates near the maximum. The data covered a range of inlet stagnation pressures from slightly above saturation to twice the thermodynamic critical pressure. The data were taken along five separate inlet stagnation isotherms ranging from 0.75 to 1.035 times the thermodynamic critical temperature. The results indicate that subcooled liquids do not choke or approach maximum flow in an asymptotic manner even though the back pressure is well below saturation; and orifice flow coefficients are not constant as is frequently assumed. A metastable jet appears to exist which breaks down if the difference between back pressure and saturation pressure is large enough.

Description: A two-phase mathematical model is proposed for calculating the induced turbulent vertical liquid flow. Bubbles provide a large buoyancy force and the associated drag on the liquid moves the liquid upward. The liquid pumped upward consists of the bubble wakes and the liquid brought into the jet region by turbulent entrainment. The expansion of the gas bubbles as they rise through the liquid is taken into account. The continuity and momentum equations are solved numerically for an axisymmetric air jet submerged in water. Water pumping rates are obtained as a function of air flow rate and depth of submergence. Comparisons are made with limited experimental information in the literature.

Description: Analytical investigations of fluid dynamics problems of relevance to the gaseous core nuclear reactor program are presented. The vortex type flow which appears in the nuclear light bulb concept is analyzed along with the fluid flow in the fuel inlet region for the coaxial flow gaseous core nuclear reactor concept. The development of numerical methods for the solution of the Navier-Stokes equations for appropriate geometries is extended to the case of rotating flows and almost completes the gas core program requirements in this area. The investigations demonstrate that the conceptual design of the coaxial flow reactor needs further development.

Description: The effect of noncondensable gases on high performance arterial heat pipes has been investigated both analytically and experimentally. Models have been generated which characterize the dissolution of gases in condensate and the diffusional loss of dissolved gases from condensate in arterial flow. These processes, and others, have been used to postulate stability criteria for arterial heat pipes. Experimental observations of gas occlusions were made using a stainless steel heat pipe equipped with viewing ports, and the working fluids methanol and ammonia with the gas additives helium, argon, and xenon. Observations were related to gas transport models.

Description: The laminar boundary layer has been theoretically studied for six gases for flows over cold walls with zero pressure gradient at Mach numbers between 5.5 and 12.5 to correlate boundary layer quantities for the various gases. The flow conditions considered correspond to those that can be generated in test facilities such as the shock tunnel and the expansion tube. Computed results obtained using real gas properties indicate that the Eckert number based on edge conditions serves to correlate the results in terms of the wall shear stress and enthalpy gradient, the Stanton number, and the momentum thickness for the various gases within plus or minus 10 per cent for Te = Tw and Te approximately 3Tw. Computed Reynolds analogy factors exhibit very good agreement with those predicted by the Colburn analogy. Velocity and displacement thicknesses correlate well with Eckert number for Te = Tw, but fail to correlate for Te approximately 3Tw. Differences in results are traced to property variations. Results show that the Eckert number is a significant correlating variable for the flows considered.

Description: Liquid zones of water, soap solution and soap foam were deployed between two aligned circular disks which were free to rotate about the zone axis in the microgravity environment of Skylab IV. Such a configuration is of interest in the containerless handling of melts for possible future space processing crystal growth experiments. Three basic types of zone surface deformation and instability were observed for these rotational conditions; axisymmetric shape changes under single disk rotation, nonaxisymmetric, whirling, C-modes for long zones with equal rotation of both disks, and capillary wave phenomena for short zones with equal rotation of both disks. The sources of these instabilities and the conditions promoting them are analyzed in detail from video tape recordings of the Skylab experiments.

Description: On October 4, 1974 the International Heat Pipe Experiment was launched aboard a Black Brant Sounding Rocket from White Sands, New Mexico. The flight provided six minutes of near zero gravity during which a total of ten separate heat pipe experiments were performed. The fifteen heat pipes which were tested represent some of the latest American and European technology. This flight provided the first reported zero gravity data on cryogenic and flat plate vapor chamber heat pipes. Additionally, valuable design and engineering data was obtained on several other heat pipe configurations. This paper will discuss the payload and four of the individual experiments.

Description: In the Skylab experiment a drop of water was caused to oscillate at its natural oscillation frequency. The drop was observed until its oscillation began to decay due to its internal damping. A determination of the change in amplitude with time made it possible to verify the applicability of a theoretical model. The fundamental concept of the model is that fluid surfaces tend to an equilibrium shape produced by the balance of the forces of fluid pressure and surface tension. The theoretical model is based on the assumption that the wave amplitude is small compared with the wavelength.

Description: Unsteady three-dimensional fluid flows which are characterized by low viscosity and the presence of distinct regions of high vorticity embedded in an otherwise irrotational flow are considered, taking into account the case in which an intermingling of vortex filaments appears. Incompressible flows of the considered type are simulated with the aid of an approach in which the vorticity distribution is modeled in terms of continuous closed filaments. It is attempted to track these filaments in a Lagrangian reference frame.

Description: The reported investigation represents an extension of the time-dependent solution of separated laminar flows based on the complete Navier-Stokes equations reported by MacCormack (1971) and Carter (1973). The current study includes turbulence models in conjunction with the compressible flow equations. The calculations start with a uniform flow except for values imposed along the upstream and outer boundaries. The basic numerical method is discussed along with questions concerning the exponential accuracy and the resolution of the viscous sublayer in a compressed region.

Description: A novel method for treating certain troublesome boundary conditions in the numerical solution of time-dependent incompressible viscous flow problems is presented. This new method is developed on the basis of an integral representation for the velocity vector which contains the entire kinematics of the problem, including the boundary conditions of concern. It is shown that for the exterior flow problem the free-stream condition is satisfied at infinity exactly and the need to treat a far-field condition is removed by the use of the integral representation. The distribution of a nonvelocity variable on the solid boundary - i.e., the 'extraneous' boundary condition needed for both the exterior and the interior flows - is shown to be governed by the kinematics of the problem.

Description: The two-phase flows studied were liquid and gas mixtures in a straight flow channel of circular cross-section. Boundaries between flow regimes have been defined for normogravity on coordinates of gas quality and total mass velocity; and, when combined with boundary expressions having a Froude number term, an analytical model was derived predicting boundary shifts with changes in gravity level. Experiments with air and water were performed, first in the normogravity environment of a ground laboratory and then in 'zero gravity' aboard a KC-135 aircraft flying parabolic trajectories. Data reduction confirmed regime boundary shifts in the direction predicted, although the magnitude was a little less than predicted. Pressure drop measurements showed significant increases for the low gravity condition.

Description: A major source of the noise generated by high bypass ratio aircraft fan engines is the unsteady pressures (forces and moments) acting on the fan blades. These unsteady pressures are the result of the interaction of the blades with the wakes of upstream blade rows, inlet flow distortions, and inlet turbulence. An experimental system was developed which will permit the measurement of the unsteady pressure distributions; a series of experiments were conducted to investigate these effects, and the results are compared with predictions from existing theoretical analyses.

Description: Research accomplishments in the following areas are discussed: (1) the calibration of the gas chromatograph for acetaldehyde and ethanol; (2) the development of data reduction and analysis methods; (3) the generation and analysis of experimental data for the transport of 100 ppm acetaldehyde through a cylindrical bed packed with activated carbon granules; (4) the generation and analysis of experimental data for the transport of 100 ppm ethanol through a cylindrical bed packed with activated carbon granules; and (5) a comparison of the volume adsorption capacity of activated carbon for 100 ppm concentrations of acetaldehyde, ethanol, and acetone. Activities in progress and planned activities are reviewed.

Description: A finite element analysis was used to study dynamic instability in ducts conveying high speed fluids. Ducts examined include cantilevered curved, flexibly supported, arbitrarily shaped, and composite duct systems. Partial differential equations were used to study the duct systems.

Description: A 30-tesla magnet design is studied which calls for forced convection liquid neon heat transfer in small coolant channels. The design also requires suppressing boiling by subjecting the fluid to high pressures through use of magnet coils enclosed in a pressure vessel which is maintained at the critical pressure of liquid neon. This high pressure reduces the possibility of the system flow instabilities which may occur at low pressures. The forced convection heat transfer data presented were obtained by using a blowdown technique to force the fluid to flow vertically through a resistance heated, instrumented tube.

Description: A recuperator consisting of two fluidized bed regenerators with circulating solid particles is considered for use in a Brayton cycle. These fluidized beds offer the possibility of high temperature operation if ceramic particles are used. Calculations of the efficiency and size of fluidized bed regenerators for typical values of operating parameters were made and compared to a shell and tube recuperator. The calculations indicate that the fluidized beds will be more compact than the shell and tube as well as offering a high temperature operating capability.

Description: The results of the data reduction and analysis of the Skylab fluid mechanics demonstrations are presented. All the fluid mechanics data available from the Skylab missions were identified and surveyed. The significant fluid mechanics phenomena were identified and reduced to measurable quantities wherever possible. Data correlations were performed using existing theories. Among the phenomena analyzed were: static low-g interface shapes, oscillation frequency and damping of a liquid drop, coalescence, rotating drop, liquid films and low-g ice melting. A survey of the possible applications of the results was made and future experiments are recommended.

Description: LOX and LN2 data for two-phase critical flow through nozzles were acquired with precision control. The principal measured parameters were inlet conditions, critical flow rate and critical flow pressure ratio. It is conclusively demonstrated that the principle of corresponding states can be applied to two-phase choked flow through nozzles. It is also shown that the proper normalizing parameters are developed, and that current theories can provide an adequate means for extrapolating data to other fluids.

Description: The closure problem in turbulence is reviewed, and some simple closure schemes are introduced. Processes occurring in turbulent flow are discussed on the basis of solutions for some elementary flows such as homogeneous turbulence with and without uniform shear. Closure by specification of initial conditions, and finally practical closure schemes for more complicated flows are considered. The latter include Reynolds stress, eddy viscosity, and mixing-length closures.

Description: The operating characteristics of a porous plug which has liquid helium on one side and which is pumped on under vacuum on the other side are discussed. The system investigated consists of a container of liquid helium which is well isolated, and the only means for mass flow out of the container is through a plug mode of porous material. The plug was assumed to have liquid helium on the container side while the other side of the plug is evacuated. Three cases were considered: (1) perfect evacuation with zero pressure, (2) evacuation through a chocked orifice, and (3) evacuation through a long, small diameter pipe with heating. Mass flow rates were determined along with mass flow at temperature above the lambda point temperature. Solutions were obtained for normal and superfluid velocity.

Description: A finite difference formulation is presented for thermal wave propagation resulting from periodic heat sources. The numerical technique can handle complex problems that might result from variable thermal diffusivity, such as heat flow in the earth with ice and snow layers. In the numerical analysis, the continuous temperature field is represented by a series of grid points at which the temperature is separated into real and imaginary terms. Computer routines previously developed for acoustic wave propagation are utilized in the solution for the temperatures. The calculation procedure is illustrated for the case of thermal wave propagation in a uniform property semi-infinite medium.

Description: A finite difference formulation is presented for thermal wave propagation resulting from periodic heat sources. The numerical technique can handle complex problems that might result from variable thermal diffusivity, such as heat flow in the earth with ice and snow layers. In the numerical analysis, the continuous temperature field is represented by a series of grid points at which the temperature is separated into real and imaginary terms. Next, computer routines previously developed for acoustic wave propagation are utilized in the solution for the temperatures. The calculation procedure is illustrated for the case of thermal wave propagation in a uniform property semi-infinite medium.

Description: Vaporization times of mercury droplets in Leidenfrost film boiling on a flat horizontal plate are measured in an air atmosphere. Extreme care was used to prevent large amplitude droplet vibrations and surface wetting; therefore, these data can be compared to film boiling theory. For these data, diffusion from the upper surface of the drop is a dominant mode of mass transfer from the drop. A closed-form analytical film boiling theory is developed to account for the diffusive evaporation. Reasonable agreement between data and theory is seen.

Description: A procedure for numerical solution of the time-dependent, two-dimensional incompressible Navier-Stokes equations that can treat the unsteady laminar flow about bodies of arbitrary shape, such as two-dimensional airfoils, multiple airfoils, and submerged hydrofoils, as naturally as it can deal with the flow about simple bodies. The solution is based on a method of automatic numerical generation of a general curvilinear coordinate system with coordinate lines coincident with all boundaries of a general multiconnected region containing any number of arbitrarily shaped bodies. The curvilinear coordinates are generated as the solution of two elliptical partial differential equations with Dirichlet boundary conditions, one coordinate being specified to be constant on each of the boundaries, and a distribution of the other being specified along the boundaries. The solution compares excellently with the Blasius boundary layer solution for the flow past a semiinfinite flat plate.

Description: A numerical method for solving the compressible form of the unsteady Navier-Stokes equations is described. This method was originally presented in 1970 and has since been modified during the development of computer programs at Ames for implementing models that account for the effects of turbulence in shock-induced separated flows. Although this paper does not describe the turbulence models themselves, a complete description of the basic numerical method is given with emphasis on the choice of a computational mesh for high Reynolds number flows, finite-difference approximations for mixed partial derivatives, extension of the Courant-Friedrichs-Lewy stability condition for viscous flows, mesh boundary conditions, and numerical smoothing for strong shock-wave calculations.

Description: The design, fabrication, and evaluation of a full scale shuttle-type condensing heat exchanger constructed of aluminum and utilizing aluminum clad titanium parting sheets is described. A long term salt spray test of candidate parting sheet specimens is described. The results of an investigation into an alternate method of making composite sheet material are discussed.

Description: The purpose of the present study is the numerical simulation of transition to turbulence in a boundary layer. Starting from a three-dimensional disturbance, the evolution of the Navier-Stokes solution until breakdown is followed. Qualitative comparison with available experimental observations is carried out.

Description: A simplified expression to estimate surface temperatures in forced convection boiling was developed using a liquid nitrogen data base. Using the principal of corresponding states and the Kutateladze relation for maximum pool boiling heat flux, the expression was normalized for use with other fluids. The expression was applied also to neon and water. For the neon data base, the agreement was acceptable with the exclusion of one set suspected to be in the transition boiling regime. For the water data base at reduced pressure greater than 0.05 the agreement is generally good. At lower reduced pressures, the water data scatter and the calculated temperature becomes a function of flow rate.

Description: Wavy wall experiments using solid waves and progressive waves have been reported. For this paper, the major effects of waviness of the wall on the flow are identified as due to oscillatory curvature (convex-concavity) and oscillatory acceleration/deceleration of the flow, which imposes a highly nonequilibrium influence upon the turbulence structure. The theoretical analysis in this presentation takes into account proper turbulence modeling (including the nonequilibrium effects) for the wavy wall problem. The analysis proceeds in three stages: (1) inviscid solution for induced pressure due to the physical wall, (2) solution of a turbulent boundary layer with pressure gradients and curvature effects in the modeling from which the profile correction is computed, and (3) induced pressure computations for the corrected profile. The phase shift of pressure perturbations with respect to the physical wall can be predicted, and pressure drag and skin friction drag can be estimated, with nonlinear viscous effects included. Comparison of the theoretical estimates with experimental data are also presented.

Description: Thermal convection driven by uniform volumetric energy sources was studied in a horizontal fluid layer bounded from above by a rigid, isothermal surface and from below by a rigid, zero heat-flux surface. The side walls of the fluid domain were assumed to be rigid and perfectly insulating. The computations were formally restricted to two-dimensional laminar convection but were carried out for a range of Rayleigh numbers which spans the regimes of laminar and turbulent flow. The results of the computations consists of streamline and isotherm patterns, horizontally averaged temperature distributions, and horizontally averaged Nusselt numbers at the upper surface. Flow and temperature fields do not exhibit a steady state, but horizontally averaged Nusselt numbers reach limiting, quasi-steady values for all Rayleigh numbers considered. Correlations of the Nusselt number in terms of the Rayleigh and Prandtl numbers were determined.

Description: An analytical and experimental investigation of the effect of vibration on the retention characteristics of screen acquisition systems was performed. The functioning of surface tension devices using fine-mesh screens requires that the pressure differential acting on the screen be less than its pressure retention capability. When exceeded, screen breakdown will occur and gas-free expulsion of propellant will no longer be possible. An analytical approach to predicting the effect of vibration was developed. This approach considers the transmission of the vibration to the screens of the device and the coupling of the liquid and the screen in establishing the screen response. A method of evaluating the transient response of the gas/liquid interface within the screen was also developed.

Description: The plane motion of an articulated pipe made of two segments is examined and the flow velocity at which flutter manifests itself is sought. The pressure in the reservoir feeding the pipe is kept constant. In contrast to previous works, the flow velocity is not taken as a prescribed parameter of the system but is left to follow the laws of motion. This approach requires a nonlinear formulation of the problem and the equations of motion are solved using Krylov-Bogoliubov's method. A graph of the amplitude of the limit cycles, as a function of the fluid-system mass ratio, is presented and conclusions are drawn as to the necessity of considering nonlinearities in the analysis.

Description: Hydrodynamic effects are analyzed for a stepped piston moving within a tight clearance tube filled with an incompressible fluid. Together with the hydrostatic effects that were analyzed in an earlier paper, a complete solution is obtained and an optimum step design for centering of the piston is suggested. The axial speed resulting from an axial driving force is calculated, and some experimental results for pistons falling in a water filled tube are presented.

Description: The status of the continuing compliant wall drag reduction research at NASA-Langley Research Center is discussed. Preliminary surface motion calculations are reported along with compliant surface design concepts and their numerical models. A compliant drag reduction theory based on stabilizing the turbulent substructure is proposed and previous experiments have been examined relative to that theory. Results of recent low speed compliant surface experiments have been reported which include initial attempts to measure local compliant surface motion.

Description: A code has been written to describe the behavior of an ablating heat shield throughout the entire maneuver into a planetary atmosphere. The code includes a trajectory computation for a variable mass, a computation of the heating history and distribution over the probe, ablation effects, indepth material response, and shape change. The code draws the initial and final shapes of the heat shield and computes its initial and final mass. Results are compared with those of other investigators, and solutions are obtained for a silica heat shield entering three model atmospheres of Jupiter. The effects of varying half cone angle and entry angle are examined as well as shape change effects. The so-called 'cusping effect' in the stagnation region is examined.

Description: Test data obtained for a low temperature phase change material (PCM) canisters are presented. The canister was designed to provide up to 30 w-hrs of storage capacity at approximately -90 C with an overall thermal conductance which is greater than 8 w/deg C. N-heptane which is an n-paraffin and has a -90.6 C freezing point was used as the working fluid. The canister was fabricated from aluminum and has an aluminum honeycomb core. Its void volume permits service temperatures up to 70 C. Results obtained from component and system's tests indicate well defined melting and freezing points which are repeatable and within 1 C of each other. Subcooling effects are less than 0.5 C and are essentially negligible. Measured storage capacities are within 94 to 88% the theoretical.

Description: This paper summarizes the 'state-of-the-art' of axially grooved heat pipes. Applications are identified and the related heat pipe design and performance are defined. Recent developments in the analysis, design and fabrication of axially grooved hardware are also discussed. A mathematical model which predicts the hydrodynamic behavior and accounts for liquid recession, liquid/vapor shear interaction and 1-g puddle flow is also presented. Performance data for various fluids in the 100-500 K range is compared to predictions from the Groove Analysis Program (GAP). Finally, a simplified closed form solution which accounts for gravity effects, self-priming and composite pumping by the grooves as well as all of the hydrodynamic losses is also discussed.

Description: The GASJET nose tip is proposed to reduce the erosion at the apex of a missile nose flying at hypersonic speeds through a rain storm as it reenters the atmosphere. A forward facing sonic jet is directed through the tip introducing a secondary counter-flow which displaces the bow shock and blankets the tip with a protective layer of relatively cool gas. Wind tunnel experiments are described which proved the validity of measuring nose recession in flight by recording the pressure in the blast tube supplying the GASJET.

Description: An evaporative surface is described for heat pipes and other two-phase heat transfer applications that consists of a hybrid composition of V-grooves and capillary wicking. Characteristics of the surface include both a high heat transfer coefficient and high heat flux capability relative to conventional open faced screw thread surfaces. With a groove density of 12.6/cm and ammonia working fluid, heat transfer coefficients in the range of 1 to 2 W/sq cm K have been measured, along with maximum heat flux densities in excess of 20 W/sq cm. A peak heat transfer coefficient in excess of 2.3 W/sq cm K at 20 W/sq cm was measured with a 37.8/cm hybrid surface.

Description: This paper describes theoretical and experimentally verified heat pipe characteristics of an axially grooved aluminum extrusion with a re-entrant groove profile. The extrusion is 13 mm diameter with 20 axial grooves, each groove consisting of a nominal .8 mm diameter channel with a .2 mm wide passageway connecting the channel to the hollow core. A computer program was written to compute the zero gravity heat transport capability of the extrusion. A heat pipe was fabricated and its performance characteristics measured. The characteristics of the pipe with ammonia at 20 C are: zero gravity pumping limit 143 w-meters; static wicking height 21.5 mm; evaporator and condenser coefficients 7300 and 20,500 watt/sq m C, respectively.

Description: An investigation was conducted to study the concept of thermal control by means of physical or chemical reaction heats for applications involving the storage of cryogens during long-term space voyages. The investigation included some preliminary experimental tests of energy storage material (ESM) effectiveness. The materials considered can store and liberate large amounts of thermal energy by means of mechanisms such as sensible heat, heat of fusion, and physical or chemical reaction heat. A differential thermal analysis was utilized in the laboratory tests. Attention is given to the evaluation of cryogenic ESM thermal control concepts, the experimental determination of phase change materials characteristics, and adsorption ESMs. It is found that an ESM shield surrounded by multiple layer insulation provides the best protection for a cryogen store.

Description: A transient thermal analysis and test of a thermal control canister is described. The 1 x 1 x 3 m canister provides a uniform thermal environment for shuttle instrument payloads requiring fine temperature control, the design goal being operation between 0 C and 20 C with a range of plus or minus 1 C at any selected set-point temperature. The canister side walls are isothermalized by a system of longitudinal and circumferential heat pipes rejecting heat through feedback controlled, variable conductance heat pipes to side mounted radiators. A breadboard model of two side walls and two radiators was tested in a thermal vacuum chamber. The breadboard was stable over a wide range of effective environments, experiment dissipations, and control point temperature levels.

Description: The decomposition and corrosion of two-phase heat transfer liquids and metal envelopes have been investigated on the basis of molecular, bond strengths and chemical thermodynamics. Potentially stable heat transfer fluids for the temperature range 100 to 350 C have been identified, and reflux heat pipe tests initiated with 10 fluids and carbon steel and aluminum envelopes to experimentally establish corrosion behavior and noncondensable gas generation rates.

Description: Mathematical models were developed for one-, two-, and three-stage radiator systems to determine optimum stage areas and system performance as a function of such parameters as insulation effectiveness, cold stage temperature, and heat load to the cold and intermediate stages. This study shows that multistage radiator systems can be optimized on the basis of weight or projected area, and that cold stage temperatures as low as 15 K are theoretically possible with present technology levels for insulation emittance. For the baseline design, analyses were performed to determine optimum radiator fin geometry and heat pipe spacing as a function of temperature, material properties, and heat pipe weight. In addition, a ground test system was designed for the baseline design with heat rejection requirements of 10 MW at 35 K on the cold stage and 100 MW at the second stage.

Description: A mathematical model is developed of excess liquid in heat pipes that is used to calculate the parameters governing the axial flow of liquid in fillets and puddles that form in vapor spaces. In an acceleration field, the hydrostatic pressure variation is taken into account, which results in noncircular meniscus shapes. The two specific vapor-space geometries considered are circular and the 'Dee-shape' that is formed by a slab wick in a circular tube. Also presented are theoretical and experimental results for the conditions under which liquid slugs form at the ends of the vapor spaces. These results also apply to the priming of arteries.

Description: The entrainment-shear performance limit in heat pipes is investigated. The entrainment heat flux limit is defined as the condition where the Weber number is greater than or equal to one. In this analysis, the critical value for the entrainment Weber number is between 2pi and 3pi. Performance degradation due to vapor-liquid shearing stress is predicted. Preliminary qualitative experiments were conducted to observe the shear stress wave formation phenomena. The equations presented may be used to predict and minimize the vapor-liquid shear stress performance effects in axial groove and puddle flow artery heat pipes.

Description: Much attention is currently being focused on turbulence modeling by separation of large-scale and small-scale turbulence. One major difficulty with this method is the problem of separating and defining the large-scale structure. In this paper, a method is presented which does this separation by analyzing time-averaged, two-point double velocity correlation measurements. The basis of the method is a maximization of the inner product of a candidate velocity field with the measured turbulent velocity field. An eigenvalue problem results, whose solution provides a set of orthogonal eigenfunctions associated with specific modes of turbulent motion. The large-scale eddy, containing most of the turbulent energy, is associated with the first (largest) mode. The turbulent velocity field is then represented as a sum of the eigenfunctions. A short description of the equations is given, the methods of solution are discussed, and some results for a two-dimensional fully developed turbulent boundary layer are presented.

Description: A fourth-order box method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method is the natural extension of the second-order Keller Box Scheme to fourth order and is demonstrated with application to the incompressible, laminar and turbulent boundary-layer equations. The efficiency of the present method is compared with other two-point and three-point higher-order methods; namely, the Keller Box Scheme with Richardson extrapolation, the method of deferred corrections, and the three-point spline methods. For equivalent accuracy, numerical results show the present method to be more efficient than the other higher-order methods for both laminar and turbulent flows.

Description: An implicit finite difference scheme is developed for the numerical solution of the compressible Navier-Stokes equations in conservation-law form. The algorithm is second-order-time accurate, noniterative, and spatially factored. In order to obtain an efficient factored algorithm, the spatial cross-derivatives are evaluated explicitly. However, the algorithm is unconditionally stable and, although a three-time-level scheme, requires only two-time-levels of data storage. The algorithm is constructed in a 'delta' form (i.e., increments of the conserved variables and fluxes) that provides a direct derivation of the scheme and leads to an efficient computational algorithm. In addition, the delta form has the advantageous property of a steady-state (if one exists) independent of the size of the time step. Numerical results are presented for a two-dimensional shock boundary-layer interaction problem.

Description: This paper presents results from a comprehensive investigation of large-eddy simulations of homogeneous isotropic turbulence. Calculations have been made using grid meshes, a number of flowfield filters, and two simple subgrid scale turbulence models. Particular attention has been paid to the degree of isotropy and the ability of the approach to predict higher-order statistical quantities (skewness and flatness). Directions for needed improvement in the simulation approach are indicated, and a plea is made to experimenters to process their data in a way which will facilitate comparison with large eddy simulations.

Description: By using Fourier transforms for evaluating spatial derivatives, the accuracy of the large-eddy simulation of homogeneous isotropic turbulence is improved. Especially, the treatment of certain terms that arise in filtering the equations is considerably improved in both speed and accuracy. Use of vorticity as the principal variable is shown to be a viable and potentially useful alternative to the primitive variables. A method of deriving conservation properties of numerical schemes is given which is much simpler than previous methods and is widely applicable. The methods are applied to the computation of homogeneous isotropic turbulence, and it is found that the subgrid scale model is improved by using finite differences in place of 'exact' derivatives.

Description: Shock-wave-induced patterns and turbulent supersonic shear flows over a stagnant bubble were determined for cases where recirculation and final mixing occurred within a very short time and within a distance of about 1 cm. The effect of a shock wave upon stratified laminas of gases in a shock tube was studied for a series of gases of different densities and molecular weights. Steady-state mixing of gas streams were investigated in wind tunnel experiments by means of a chemiluminescent reaction. Shock-wave-induced mixing which involves the introduction of a liquid jet into a supersonic gas stream is also described.

Description: The effect, of the laminar/turbulent boundary layer state on the mean and rms velocities of a developing plane mixing layer, has been investigated. The maximum Ux/nu allowed by the facility was 6.7 x 100.000; the maximum x station surveyed corresponds to approximately 2000 and 700 initial momentum thicknesses for the laminar and turbulent cases respectively. The use of commonly accepted non-dimensional representations of the data confirm (at least) an approximately self-preserving condition. They also suggest that the effects of the laminar/turbulent initial condition persist in the self-preserving region. A direct comparison of the data reveals that the persistence so observed is illusory. An interpretation of the reason for this misunderstanding is advanced.

Description: The interaction between the large-scale structure and the disparately fine-grained turbulence in a free turbulent shear flow, in particular, the mixing region is considered. New results pertaining to the spatial problem are presented. An incompressible fluid is considered. An energy integral description is discussed, giving attention to shape assumptions, the interaction problem, energy transfer mechanisms, the spectral dependence of the interaction, and the control of the large-scale structure.

Description: Numerical solutions of the Navier-Stokes equations for shock separated turbulent boundary-layer flows are presented. Several turbulence models are investigated and assessed by their ability to predict the physical phenomena associated with two extensively documented experiments. The experimental flows consist of shock-wave boundary-layer interactions in axisymmetric internal and external geometries at Mach numbers of 1.5 and 7, respectively. Algebraic and one-equation eddy viscosity models are used to describe the Reynolds shear stress. Calculated values of skin friction, wall pressure distribution, kinetic energy of turbulence, and heat transfer are compared with measurements.

Description: A finite-element method for steady-state thermal analysis of convectively cooled structures is presented. The method is based on representing the coolant passages by finite elements with fluid bulk temperature nodes and fluid/structure interface modes. Four finite elements are described: two basic elements (mass transport and surface convection) and two integrated elements for applications to discrete tube and plate-fin convectively cooled structural configurations. Comparative finite-element and lumped-parameter thermal analyses of several convectively cooled structures demonstrate the practicality of utilizing finite-element methodology for thermal analysis of realistic convectively cooled structures.

Description: A complete second-order closure program for the investigation of chemically reacting turbulent flowfields is being developed. Many of the turbulence models being used have been developed and tested during previous studies of a variety of incompressible and compressible shear flows. A number of additional scalar correlations have to be modeled for chemically reacting flows, and a 'typical eddy' model has been developed for the joint probability density function for all the scalars. A simpler model in which higher-order scalar correlations are set zero is also being used. The program has been recently extended to multi-species, multi-step chemical reactions. Computations have been performed using the second-order closure program to study the turbulence-chemistry interaction in DF chemical laser flows and in a hydrogen-air diffusion flame. The results demonstrate the importance of the inclusion of the concentration and temperature fluctuation correlations in turbulent reacting flows.

Description: Data were obtained for two phase and gaseous choked flow nitrogen in a long constant area duct of 16200 L/D with a diverging diffuser attached to the exit. Flow rate data were taken along five isotherms (reduced temperature of 0.81, 0.96, 1.06, 1.12, and 2.34) for reduced pressures to 3. The flow rate data were mapped in the usual manner using stagnation conditions at the inlet mixing chamber upstream of the entrance length. The results are predictable by a two phase homogeneous equilibrium choking flow model which includes wall friction. A simplified theory which in essence decouples the long tube region from the high acceleration choking region also appears to predict the data resonably well, but about 15 percent low.

Description: A design analysis, is developed based on experimental data, to predict the effects of transient flow and pressure surges (caused either by valve or pump operation, or by boiling of liquids in warm lines) on the retention performance of screen acquisition systems. A survey of screen liquid acquisition system applications was performed to determine appropriate system environment and classification. A screen model was developed which assumed that the screen device was a uniformly distributed composite orthotropic structure, and which accounted for liquid inflow/outflow, gas ingestion quality, screen stress, and liquid spill. A series of 177 tests using 13 specimens (5 screen meshes, 4 screen device construction/backup methods, and 2 orientations) with three test fluids (isopropyl alcohol, Freon 114, and LH2) provided data which verified important features of the screen model and resulted in a design tool which could accurately predict the transient startup performance acquisition devices.

Description: Semidirect methods are discussed, their present role, as well as some developments for their application in computational fluid dynamics. A semidirect method is a computational scheme that uses a fast, direct, elliptic solver as the driving algorithm for the iterative solution of finite difference equations. Specific subtopics include: (1) direct Cauchy Riemann solvers for first order elliptic equations; (2) application of the semidirect method to the mixed elliptic hyperbolic problem of steady, inviscid transonic flow; and (3) the treatment of interior conditions, such as those on an airfoil or wing, in semidirect methods.

Description: The effects of: (1) a reactive environment on film cooling effectiveness, and (2) film cooling on rocket engine performance were determined experimentally in a rocket thrust chamber assembly operating with hydrogen and oxygen propellants at 300 psi chamber pressure. Tests were conducted using hydrogen, helium, and nitrogen film coolants in an instrumented, thin walled, steel thrust chamber. The film cooling, performance loss, and heat transfer coefficient data were correlated with the ALRC entrainment film cooling model which relates film coolant effectiveness and mixture ratio at the wall to the amount of mainstream gases entrained with the film coolant in a mixing layer. In addition, a comprehensive thermal analysis computer program, HOCOOL, was prepared from previously existing ALRC computer programs and analytical techniques.

Description: An experimental study of two coaxial gas streams, which react very rapidly, was performed to investigate the mixing characteristics of turbulent flow fields. The center stream consisted of a CO-N2 mixture and the outer annular stream consisted of air vitiated by H2 combustion. The streams were at equal velocity (50 m/sec) and temperature (1280 K). Turbulence measurements were obtained using hot film anemometry. A sampling probe was used to obtain time averaged gas compositions. Six different turbulence generators were placed in the annular passage to alter the flow field mixing characteristics. The turbulence generators affected the bulk mixing of the streams and the extent of CO conversion to different degrees. The effects can be related to the average eddy size (integral scale) and the bulk mixing. Higher extents of conversion of CO to CO2 were found be increasing the bulk mixing and decreasing the average eddy size.

Description: The paper summarizes results of theoretical and experimental investigations of wall temperature change and heat transfer in the laminar boundary layer formed within noncentered unsteady expansion waves in shock tubes. The study is restricted to a class of noncentered plane waves which have one key feature in common with centered waves, i.e., the first derivatives of the inviscid flow quantities have values or are discontinuous at the wavehead. The wall temperature is calculated implicitly by matching the local heat transfer rates of gas and wall thermal boundary layers. Experimentally, the wall temperature was measured using a thin-film resistance thermometer and heat transfer was then deduced using a one-dimensional unsteady heat conduction equation.

Description: A one-dimensional nonequilibrium flow analysis has been formulated for a one component two phase flow. The flow is considered homogeneous and essentially isothermal. Phase change is assumed to occur at heterogeneous nucleation sites and the growth of the vapor bubbles is governed by heat conduction from the liquid to the bubble. The analysis adjusted for friction is applied to liquid nitrogen flow in a venturi and comparison is made with the NASA experimental results of Simoneau. Good agreement with the experiments is obtained when one assumes the effective activation energy for nucleus formation to be small but nonzero. The computed pressure distributions deviate from the experimental results in the throat region of the venturi in a manner consistent with centrifugal effects not accounted for in the one-dimensional theory. The results are shown to depend not only on cavitation number but on additional dimensionless parameters governing the nonequilibrium production and subsequent growth of nuclei.