ALBERT

Description: We study the onset of a pure Marangoni convection in a liquid layer with two deformable interfaces in the no-gravity environment. Both oscillatory and stationary instabilities are considered for a wide range of parameters. It is shown that only stationary instability is possible when surface tension at the colder interface is lower than that at the hotter one. Oscillatory instability tends to disappear and to be replaced by the stationary instability with increase of the Prandtl number and decrease of surface tension at the colder interface.

Description: A consistent solution of the radiative transfer equation characterizing photon transport in a semi-infinite medium of refractive index greater than or equal to one is obtained following the method of Sobolev. Fresnel specular reflection, Snell's law and isotropic scattering are assumed. An algorithm is developed and its accuracy is demonstrated. A numerical Laplace transform inversion leads to an efficient evaluation for the interior flux and source function distributions.

Description: Results are reported of the Surface Tension Driven Convection Experiment (STDCE) aboard USML-1 Spacelab. Steady and transient thermocapillary flows were investigated in a 10 cm dia. circular container filled with 10 Cs silicone oil. The velocity and temperature fields were studied in detail under various conditions. It is shown in this paper how the Marangoni number affects the velocity field. A numerical analysis of the flows was also conducted and its results were compared to the experimental data. Good agreement is shown.

Description: Anchorage dependent cell cultures in fluidized beds are tested. Feasibility calculations indicate the allowed parameters and estimate the shear stresses therein. In addition, the diffusion equation with first order reaction is solved for the spherical shell (double bubble) reactor with various constraints.

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: Technical improvements of a long life heat rejection system, suitable for long duration high power missions, that can be constructed and deployed in orbit is discussed. A mathematical model is formulated and a computer program developed which describes the transient priming characteristics of a dual passage heat pipe. An experimental test package is described for flight in the KC-135 Zero-g Aircraft, to be used to verify the modeling predictions.

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: Spacelab experiment to investigate two-phase flow patterns under gravity uses a water-air mixture experiment. Air and water are circulated through the system. The quality or the mixture or air-water is controlled. Photographs of the test section are made and at the same time pressure drop across the test section is measured. The data establishes a flow regime map under reduced gravity conditions with corresponding pressure drop correlations. The test section is also equipped with an electrical resistance heater in order to allow a flow boiling experiment to be carried out using Freon II. High-speed photographs of the test section are used to determine flow patterns. The temperature gradient and pressure drop along the duct can be measured. Thus, quality change can be measured, and heat transfer calculated.

Description: The equations of motion governing an incompressible fluid contained in an orbiting laboratory were examined to isolate various fictitious forces and their relative influence on the fluid. The forces are divided into those arising from the orbital motions and those arising from small local motions of the spacecraft about its center of mass. The latter dominate the nonrotating experiments. Both are important for rotating experiments. A brief discussion of the onset of time-dependence and violent instability in earth-based rotating and processing systems is given.

Description: Natural convection is not always harmful and, therefore, to be avoided. In some situations it may be desirable to have fluid flows in space processes, e.g., to stir the fluid phase for mixing and cooling or to help maintain concentration gradients. In may event, it is important to know the extent and nature of convection in space and the factors on which it depends, in order either to minimize the effects to convection, or to utilize the convection to advantage. The information needed to assess both conventional and unstable convection includes: (1) the magnitude and direction of accelerations; (2) geometric configuration; (3) imposed boundary conditions; and (4) material properties.

Description: Systematic scaling or dimensional analysis reveals that certain scales of geophysical fluid flows (such as stellar, ocean, and planetary atmosphere circulations) can be accurately modeled in the laboratory using a procedure which differs from conventional engineering modeling. Rather than building a model to obtain numbers for a specific design problem, the relative effects of the significant forces are systematically varied in an attempt to deepen understanding of the effects of these forces. Topics covered include: (1) modeling a large-scale planetary atmospheric flow in a rotating cylindrical annulus; (2) achieving a radial dielectric body force; (3) spherical geophysical fluid dynamics experiments for Spacelab flights; (4) measuring flow and temperature; and (5) the possible effect of rotational or precessional disturbances on the flow in the rotating spherical containers.

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: After the external tank separates from the Orbiter about 2000 pounds of residual liquid oxygen remain in the main propulsion system lines. The pressurization of liquid oxygen from a subcritical to a supercritical state by the use of the heaters of the PRSA tanks while in a low-g environment is investigated. The performance of the heaters while bringing the state of the substance from the subcritical state to the supercritical one is studied, with particular emphasis on the time the pressurization process takes, and the temperature of the heater as the process proceeds.

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: This is a review of the influence of convection on the growth of crystals from solution. The growth rate is increased by convection up to the point where interface kinetics becomes rate controlling. Compositional inhomogeneity and morphological instability (inclusion formation) are probably worse for gentle convection than for either no convection or for vigorous stirring. Stirring, particularly of crystal suspensions, can cause an orders of magnitude increase in the rate of formation of new crystals. This is called 'secondary nucleation'.

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: Attention is given to the way in which external turbulence affects an initially turbulence-free region in which there is a mean velocity gradient. External turbulence induces irrotational fluctuations in the sheared region which interact with the shear to produce rotational velocity fluctuations and mean Reynolds stresses. Since the actual front between the initial external turbulence and the shear flow is a randomly contorted surface, the turbulence near the front is intermittent, and is presently included in the form of a simple statistical model. In wind tunnel tests, turbulent shear stress was found to grow from zero to significant values in the interaction region. Observed stress magnitude and extent agrees with predictions, and it is concluded that turbulent stresses can be produced by irrotational fluctuations in a region of mean shear.

Description: A liquid, contained in a quarter plane, undergoes steady motion due to thermocapillary forcing on its upper boundary, a free surface separating the liquid from a passive gas. The rigid vertical sidewall has a strip whose temperature is elevated compared with the liquid at infinity. A boudnary-layer analysis is performed that is valid for large Marangoni numbers M and Prandtl numbers P. It is found that the Nusselt number N for the horizontal heat transport satisfies N proportional to min (M to the 1 2/7/power, M to the 1 1/5/power, M to the 1 1/10/power) Generalizations are discussed.

Description: Six flat-disk models made of carbon-carbon and carbon-phenolic materials were launched in an argon-filled track-range facility to test ablation characteristics in a radiation-dominated, massive-blowing environment. The shock standoff distances deduced from the shadowgraphs agree with theoretical predictions during the earlier portion of the flight, while the wall temperatures determined by the image-converter photographs agree with predictions during the later portion. The measured surface recessions exceed the calculated values by about 60 percent for carbon-phenolic and 30 percent for carbon-carbon. The discrepancies are attributed to spallation. The measured char thicknesses agree with theoretical predictions.

Description: The results of modifications in continuation methods applied to obtain solutions to the Navier-Stokes systems of equations for incompressible, two-dimensional, steady flows are reported. It is shown that parameter continuation permits prediction of accurate, initial estimates for iterative processing of nonlinear finite difference and finite element equations of motions. The new parameter steps are derived from values of the preceding parameter steps. The accuracy of the estimates is ensured with appropriate choices of the step size. The continuation predictor/iterative corrector is demonstrated to trace the branches of parameter space along which steady flow states are found, and techniques are available for tracing multiply branching paths. The techniques are applied to solving the Navier-Stokes equations for flow through a rotating square channel, the formation of a falling liquid curtain, and gyrostatic equilibria of rotating cylindrical drops.

Description: The four-point, centered implicit scheme that is extensively used in open channel flow simulation is shown to be applicable to rapid and slow pressure transient problems in conduits with nearly single phase and two-phase flows. It is only necessary to choose the proper weighting factor value, theta, of the Courant number. For rapid pressure transients such as waterhammer, the implicit method can yield reasonable results with limited numerical dispersion and attenuation if theta is only slightly greater than the critical value of 0.5. For slower pressure gradients in single and two-phase flows, reasonable numerical solutions may be achieved for Courant number values as high as 20.

Description: This paper examines, both theoretically and experimentally, the effect produced by irrotational fluctuations, associated with a nearby turbulent field, in a region where the turbulence is initially very low but where there is a mean shear. Calculations are based on rapid distortion theory and experiments use linearized hot wire anemometers in an open circuit wind tunnel. Turbulent shear stress is observed to grow from zero to significant values in the interaction region. The magnitude and extent of this observed shear stress agree reasonably well with predictions of the analysis, when intermittency effects are included. It is concluded that turbulent stresses can be produced by irrotational fluctuations in a region of mean shear and that this effect can be estimated using rapid distortion theory if the overall strain ratio is not large.

Description: A revised version of Dodge's split-velocity method for numerical calculation of compressible duct flow was developed. The revision incorporates balancing of mass flow rates on each marching step in order to maintain front-to-back continuity during the calculation. The (checkerboard) zebra algorithm is applied to solution of the three dimensional continuity equation in conservative form. A second-order A-stable linear multistep method is employed in effecting a marching solution of the parabolized momentum equations. A checkerboard iteration is used to solve the resulting implicit nonlinear systems of finite-difference equations which govern stepwise transition. Qualitative agreement with analytical predictions and experimental results was obtained for some flows with well-known solutions. Previously announced in STAR as N82-16363

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 systems of truncated differential equations that have been proposed to reduce the complexity and large computational costs of solutions to the full Navier-Stokes equations are considered. These systems are computationally efficient and capture all the physically relevant behavior. The systems follow a certain hierarchy: (1) the classical boundary-layer equations with specified edge properties (usually the streamwise pressure distribution); (2) the coupled boundary-layer/inviscid equations; (3) the so-called thin-layer equations that discard streamwise diffusion; and (4) the Navier-Stokes equations. Consideration is given to each of these approximations applied to an incompressible, laminar-separating flow at low and moderate Reynolds numbers. It is pointed out that for any flow or region of flow for which viscous-inviscid interaction effects are small, classical boundary-layer equations will provide a satisfactory description of the viscous flow at a fraction of the computational cost of any higher approximations.

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: An investigation of the growth of the three-dimensional, counter-rotating, longitudinal type vortices is considered in two-dimensional laminar compressible boundary-layer flow. The basic approximation of the disturbance equations that includes the terms due to boundary layer growth is considered and solved numerically. These terms are shown to have large local effects near the neutral stability region. The study shows that the instability of the boundary layer with respect to the three-dimensional vortices sets in at higher Goertler number as Mach number increases. Also the maximum amplitude ratio of the vortices is reduced by about 20 percent as Mach number increases from 0 to 5.

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: Eleven tests were carried out in DLR's high enthalpy tunnel in Goettingen (HEG), Germany, for reservoir conditions ranging from 10 to 23 MJ/kg and a free stream Mach number of approximately 10. A blunted cone model with a cylindrical afterbody (sting) was investigated. To obtain information on the influence of defined parameters on the body back flow, especially of the reacting gas, the heat transfer rate along the body contour was measured with fast response surface thermocouples on the forebody and sensitive thin film heat transfer gauges on the base and sting of the model. Flow visualization with a holographic interferometry system was provided. Tests are described, and a preliminary interpretation of the observed flow effects are given.

Description: Get Away Special (GAS) payload G-093, also called VORTEX (VOrtex Ring Transit EXperiment), is an investigation of the propagation of a vortex ring through a liquid-gas interface in microgravity. This process results in the formation of one or more liquid droplets similar to earth based liquid atomization systems. In the absence of gravity, surface tension effects dominate the drop formation process. The Shuttle's microgravity environment allows the study of the same fluid atomization processes as using a larger drop size than is possible on Earth. This enables detailed experimental studies of the complex flow processes encountered in liquid atomization systems. With VORTEX, deformations in both the vortex ring and the fluid surface will be measured closely for the first time in a parameters range that accurately resembles liquid atomization. The experimental apparatus will record images of the interactions for analysis after the payload has been returned to earth. The current design of the VORTEX payload consists of a fluid test cell with a vortex ring generator, digital imaging system, laser illumination system, computer based controller, batteries for payload power, and an array of housekeeping and payload monitoring sensors. It is a self-contained experiment and will be flown on board the Space Shuttle in a 5 cubic feet GAS canister. The VORTEX Project is entirely run by students at the University of Michigan but is overseen by a faculty advisor acting as the payload customer and the contact person with NASA. This paper summarizes both the technical and programmatic aspects of the VORTEX Project.

Description: This paper focuses on the flight results of the Cryogenic Two-Phase Flight Experiment (CRYOTP), which was a Hitchhiker based experiment that flew on the space shuttle Columbia in March of 1994 (STS-62). CRYOTP tested two new technologies for advanced cryogenic thermal control; the Space Heat Pipe (SHP), which was a constant conductance cryogenic heat pipe, and the Brilliant Eyes Thermal Storage Unit (BETSU), which was a cryogenic phase-change thermal storage device. These two devices were tested independently during the mission. Analysis of the flight data indicated that the SHP was unable to start in either of two attempts, for reasons related to the fluid charge, parasitic heat leaks, and cryocooler capacity. The BETSU test article was successfully operated with more than 250 hours of on-orbit testing including several cooldown cycles and 56 freeze/thaw cycles. Some degradation was observed with the five tactical cryocoolers used as thermal sinks, and one of the cryocoolers failed completely after 331 hours of operation. Post-flight analysis indicated that this problem was most likely due to failure of an electrical controller internal to the unit.

Description: The Capillary Pumped Loop Flight Experiment (CAPL) employs a passive two-phase thermal control system that uses the latent heat of vaporization of ammonia to transfer heat over long distances. CAPL was designed as a prototype of the Earth Observing System (EOS) instrument thermal control systems. The purpose of the mission was to provide validation of the system performance in micro-gravity, prior to implementation on EOS. CAPL was flown on STS-60 in February, 1994, with some unexpected results related to gravitational effects on two-phase systems. Flight test results and post flight investigations will be addressed, along with a brief description of the experiment design.

Description: This paper presents an overview of the Cryogenic Test Bed (CTB) experiments including experiment results, integration techniques used, and lessons learned during integration, test and flight phases of the Cryogenic Heat Pipe Flight Experiment (STS-53) and the Cryogenic Two Phase Flight Experiment (OAST-2, STS-62). We will also discuss the Cryogenic Flexible Diode Heat Pipe (CRYOFD) experiment which will fly in the 1996/97 time frame and the fourth flight of the CTB which will fly in the 1997/98 time frame. The two missions tested two oxygen axially grooved heat pipes, a nitrogen fibrous wick heat pipe and a 2-methylpentane phase change material thermal storage unit. Techniques were found for solving problems with vibration from the cryo-collers transmitted through the compressors and the cold heads, and mounting the heat pipe without introducing parasitic heat leaks. A thermally conductive interface material was selected that would meet the requirements and perform over the temperature range of 55 to 300 K. Problems are discussed with the bi-metallic thermostats used for heater circuit protection and the S-Glass suspension straps originally used to secure the BETSU PCM in the CRYOTP mission. Flight results will be compared to 1-g test results and differences will be discussed.

Description: With the recent focus on the needs of design and applications CFD, research groups have begun to address the traditional bottlenecks of grid generation and surface modeling. Now, a host of emerging technologies promise to shortcut or dramatically simplify the simulation process. This paper discusses the current status of these emerging technologies. It will argue that some tools are already available which can have positive impact on portions of the design cycle. However, in most cases, these tools need to be integrated into specific engineering systems and process cycles to be used effectively. The rapidly maturing status of unstructured and Cartesian approaches for inviscid simulations makes suggests the possibility of highly automated Euler-boundary layer simulations with application to loads estimation and even preliminary design. Similarly, technology is available to link block structured mesh generation algorithms with topology libraries to avoid tedious re-meshing of topologically similar configurations. Work in algorithmic based auto-blocking suggests that domain decomposition and point placement operations in multi-block mesh generation may be properly posed as problems in Computational Geometry, and following this approach may lead to robust algorithmic processes for automatic mesh generation.

Description: The last decade has witnessed a vigorous and sustained research effort on unstructured methods for computational fluid dynamics. Unstructured mesh generators and flow solvers have evolved to the point where they are now in use for design purposes throughout the aerospace industry. In this paper we survey the various mesh types, structured as well as unstructured, and examine their relative strengths and weaknesses. We argue that unstructured methodology does offer the best prospect for the next generation of computational fluid dynamics algorithms.

Description: A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree data structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: A gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment and other accepted computational results for a series of low and moderate Reynolds number flows.

Description: Grid related issues of the Chimera overset grid method are discussed in the context of a method of solution and analysis of unsteady three-dimensional viscous flows. The state of maturity of the various pieces of support software required to use the approach is considered. Current limitations of the approach are identified.

Description: This paper presents 'A view from the trenches' on CFD grid generation from a Pratt & Whitney perspective. We anticipate that other organizations have similar views. We focus on customer expectations and the consequent requirements. We enunciate a vision for grid generation, discuss issues that developers must recognize.

Description: This paper presents a perspective on the requirements that Computational Fluid Dynamics (CFD) technology must meet for its effective use in aerospace design. General observations are made on current aerospace design practices and deficiencies are noted that must be rectified for the U.S. aerospace industry to maintain its leadership position in the global marketplace. In order to rectify deficiencies, industry is transitioning to an integrated product and process development (IPPD) environment and design processes are undergoing radical changes. The role of CFD in producing data that design teams need to support flight vehicle development is briefly discussed. An overview of the current state of the art in CFD is given to provide an assessment of strengths and weaknesses of the variety of methods currently available, or under development, to produce aerodynamic data. Effectiveness requirements are examined from a customer/supplier view point with design team as customer and CFD practitioner as supplier. Partnership between the design team and CFD team is identified as an essential requirement for effective use of CFD. Rapid turnaround, reliable accuracy, and affordability are offered as three key requirements that CFD community must address if CFD is to play its rightful role in supporting the IPPD design environment needed to produce high quality yet affordable designs.

Description: The efforts in geometry modeling and grid generation at the NASA Lewis Research Center, as applied to the computational fluid dynamic (CFD) analysis of aeropropulsion systems, are presented. The efforts are mainly characterized by a focus on the analysis of components of an aeropropulsion system, which involve turbulent viscous flow with heat transfer and chemistry. Thus, this discussion will follow that characterization and will sequence through the components of typical propulsion systems consisting of inlets, compressors, combustors, turbines, and nozzles. For each component, some applications of CFD analysis will be presented to show how CFD is used to compute the desired performance information, how geometry modeling and grid generation are performed, and what issues have developed related to geometry modeling and grid generation. The discussion will illustrate the following needs related to geometry modeling and grid generation as observed in aeropropulsion analysis: (1) accurate and efficient resolution of turbulent viscous and chemically-reacting flowfields; (2) easy-to-use interfaces with CAD data for automated grid generation about complex geometries; and (3) automated batch grid generation software for use with design and optimization software.

Description: This viewgraph presentation discusses key features of current combustion system CFD modeling capabilities at GE Aircraft Engines provided by the CONCERT code; CONCERT development history; modeling applied for designing engine combustion systems; modeling applied to improve fundamental understanding; CONCERT3D results for current production combustors; CONCERT3D model of NASA/GE E3 combustor; HYBRID CONCERT CFD/Monte-Carlo modeling approach; and future modeling directions.

Description: A procedure was developed to improve the turn-around time for computational fluid dynamics (CFD) simulations of an inlet-bleed problem involving oblique shock-wave/boundary-layer interactions on a flat plate with bleed into a plenum through one or more circular holes. This procedure is embodied in a preprocessor called AUTOMAT. With AUTOMAT, once data for the geometry and flow conditions have been specified (either interactively or via a namelist), it will automatically generate all input files needed to perform a three-dimensional Navier-Stokes simulation of the prescribed inlet-bleed problem by using the PEGASUS and OVERFLOW codes. The input files automatically generated by AUTOMAT include those for the grid system and those for the initial and boundary conditions. The grid systems automatically generated by AUTOMAT are multi-block structured grids of the overlapping type. Results obtained by using AUTOMAT are presented to illustrate its capability.

Description: The role of Computational Fluid Dynamics (CFD) at Ames Research Center has expanded to address a broad range of aeronautical problems, including wind tunnel support, flight test support, design, and analysis. Balancing the requirements of each new problem against the available resources - software, hardware, time, and expertise - is critical to the effective use of CFD. Several case studies of recent applications highlight the depth of CFD capability at Ames, the tradeoffs involved in various approaches, and lessons learned in the use of CFD as an engineering tool.

Description: When computing the flow around complex three dimensional configurations, the generation of the mesh is the most time consuming part of any calculation. With some meshing technologies this can take of the order of a man month or more. The requirement for a number of design iterations coupled with ever decreasing time allocated for design leads to the need for a significant acceleration of this process. Of the two competing approaches, block-structured and unstructured, only the unstructured approach will allow fully automatic mesh generation directly from a CAD model. Using this approach coupled with the techniques described in this paper, it is possible to reduce the mesh generation time from man months to a few hours on a workstation. The desire to closely couple a CFD code with a design or optimization algorithm requires that the changes to the geometry be performed quickly and in a smooth manner. This need for smoothness necessitates the use of Bezier polynomials in place of the more usual NURBS or cubic splines. A two dimensional Bezier polynomial based design system is described.

Description: The problem of forced fluid vibrations in a partially filled spinning spherical tank is solved numerically by using the finite element method. The governing equations include Coriolis acceleration and spatially homogeneous vorticity. An exponential instability is detected in the present simulation for fill ratios below 0.5 and centrifugal acceleration to thrust ratios less than 1.7. This fictitious instability appears in the model as a result of the homogeneous vortex assumption since the free slosh equations are neutrally stable in the Liapunov sense.

Description: In the present paper, a simple numerical model is used to study the warming of the mixed layer during the early summer. It is shown that the springtime temperature increase in the layer below the mixed layer (for example, in the cold pool on a continental shelf) has a maximum value which occurs for a limiting value of the surface heat flux. This is a result of the positive feedback at large Richardson numbers between stability and vertical diffusion of heat. The springtime temperature increase in the mixed layer increases nonlinearly with surface heat flux, because of the same positive feedback. The effects of interseasonal fluctuations of the surface heat flux on the spring and summer mixed layer and deeper temperature increases can be as great as the effect of interseasonal fluctuations of the average heat flux.

Description: Shapes and stability of surface-tension-endowed drops rotating rigidly at fixed angular momentum are calculated by finite-element analysis. A new family of asymmetric two-lobed drop shapes is discovered that branches from, and rejoins, the Pik-Pichak family of symmetric two-lobed shapes. The computations are verified for axisymmetric and symmetric two-lobed drop shape by comparison with previous approximations.

Description: The influences of memory effects, coupling between velocity and temperature fluctuations and tensorial transport properties on momentum and heat transfers in turbulent flow which do not follow the Boussinesq relation are examined. It is shown that the memory effect, represented by the Lagrangian of the velocity gradients, can account for the decoupling between the flux and the gradient, while the tensorial properties of the transport coefficients allow a normal Boussinesq-type transfer with memory and anomalous counter-gradient or gradient-less transport.

Description: An experimental work is discussed whose objective was to obtain data that show the effect of temperature and temperature fluctuations on surface noise. This was accomplished experimentally by immersing a small chord airfoil in the turbulent airstream of a hot jet. The theory and experiment reported by Olsen (1976) provided a guide for designing and validating the hot jet experiment and for interpreting the data. It is shown that increased temperature causes a small decrease in the sound levels; at the same time it causes a shift in the spectra that is smaller but similar to the shift observed with subsonic hot jet noise.

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: Heat transfer characteristics are analyzed for a cooled two-dimensional porous medium having a curved boundary. A general analytical procedure is given in combination with a numerical conformal mapping method used to transform the porous region into an upper half plane. To illustrate the method, results are evaluated for a cosine shaped boundary subjected to uniform external heating. The results show the effects of coolant starvation in the thick regions of the medium, and the extent that internal heat conduction causes the heated surface to have a more uniform temperature.

Description: An interactive method is proposed for the solution of two-dimensional, laminar flow fields with identifiable regions of recirculation, such as the shear-layer-driven cavity flow. The method treats the flow field as composed of two regions, with an appropriate mathematical model adopted for each region. The shear layer is computed by the compressible boundary layer equations, and the slowly recirculating flow by the incompressible Navier-Stokes equations. The flow field is solved iteratively by matching the local solutions in the two regions. For this purpose a new matching method utilizing an overlap between the two computational regions is developed, and shown to be most satisfactory. Matching of the two velocity components, as well as the change in velocity with respect to depth is amply accomplished using the present approach, and the stagnation points corresponding to separation and reattachment of the dividing streamline are computed as part of the interactive solution. The interactive method is applied to the test problem of a shear layer driven cavity. The computational results are used to show the validity and applicability of the present approach.

Description: The shear-free turbulent boundary layer is calculated by the large-eddy simulation technique. The filtered Navier-Stokes equations are used; the method of integration employs Fourier expansions in the homogeneous directions and finite differences in the cross-stream direction. Results indicate that the simulation is capable of predicting the primary Reynolds-number effects.

Description: A model is proposed for the study of the growth and shrinkage of gas bubbles in systems containing many gas bubbles. The key feature of this model is the replacement of the bubbles by point sources of gas concentration. Calculations are performed in the simple case of an initial uniform array of bubbles of equal radii.

Description: Recent progress in the development of vortex methods and their applications to the numerical simulation of incompressible fluid flows are reviewed. Emphasis is on recent results concerning the accuracy of these methods, improvements in computational efficiency, and the development of three-dimensional methods. Simulations of several example flows which display some of the strengths and weaknesses of vortex methods are presented.

Description: The quantum mechanical technique is used to study ionic, configurational, and impurity defects in the ice surface. In addition to static calculations of the energetics of the water monomer-ice surface interactions, molecular dynamics studies were initiated. The calculations of the monomer-ice surface interaction, molecular dynamics studies were initiated. The calculations of monomer-ice surface interactions indicate that many adsorption sites exist on the ice surfaces and that the barriers between bonding sites are relatively low. Bonding on the prism face of ice is preferentially above lattice sites.

Description: The central field empirical pair potential model is applied to studying the effects of kinks, ledges, and vacancies on the absorption of water molecules from the vapor. Molecular dynamics simulations indicate that cluster and/or surface modes play a primary role in the absorption process, the flexibility of the hydrogen bond serves to decrease the energy required for structural interconversion, and the rapid distribution of added energy in a hydrogen bonded system lead to aggregate stability which greatly exceeds that predicted by static energy calculations.

Description: Areas of investigation in fluid dynamics, recommended experiments, and use of the facility for theory evaluation are discussed. Tunnel flow quality and calibration of the NTF are considered. Recent technological advances affecting tunnel design are surveyed.

Description: The conservation-law form of the inviscid gasdynamic equations has the remarkable property that the nonlinear flux vectors are homogeneous functions of degree one. This property readily permits the splitting of flux vectors into subvectors by similarity transformations so that each subvector has associated with it a specified eigenvalue spectrum. As a consequence of flux vector splitting, new explicit and implicit dissipative finite-difference schemes are developed for first-order hyperbolic systems of equations. Appropriate one-sided spatial differences for each split flux vector are used throughout the computational field even if the flow is locally subsonic. The results of some preliminary numerical computations are included.

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: Two iterative schemes based on the mixed finite element method are developed for analyzing steady natural convection in a melt adjacent to its solid phase. The simplest method decouples the calculation of the field variables and the shape of the melt/solid interface into two interlocked iterations that are performed successively. The second method uses Newton's iteration to solve simultaneously for both types of unknowns and has a quadratic convergence rate. Results for a model problem of melt and solid in a cylindrical ampoule show the Newton algorithm to be a factor of three more efficient.

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 multiple-grid algorithm for use in efficiently obtaining steady solution to the Euler and Navier-Stokes equations is presented. The convergence of a simple, explicit fine-grid solution procedure is accelerated on a sequence of successively coarser grids by a coarse-grid information propagation method which rapidly eliminates transients from the computational domain. This use of multiple-gridding to increase the convergence rate results is substantially reduced work requirements for the numerical solution of a wide range of flow problems. Computational results are presented for subsonic and transonic inviscid flows and for laminar and turbulent, attached and separated, subsonic viscous flows. Work reduction factors as large as eight, in comparison to the basic fine-grid algorithm, were obtained. Possibilities for further performance improvement are discussed. Previously announced in STAR as N83-21847

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: It is noted that several terms in the two-point spectral equation for homogeneous turbulence can be interpreted as spectral-transfer terms; that is, they represent the net rate of energy transfer into a wavenumber region from all other wavenumbers. This holds for terms associated with both turbulence and self-interaction and interaction between turbulence and mean gradients. It is not seen as obvious, however, that similar interpretations apply when the turbulence is not homogeneous. In particular, one might question the interpretation for the terms associated with turbulence self-interaction because the condition of homegeneity is generally used in making the interpretation. It is the purpose here to consider whether terms interpretable as transfer terms exist in the equations for inhomogeneous turbulence. It is found that certain terms in the two-point spectral equation can be interpreted as transfer terms.

Description: A solution to the rapid-distortion theory for small-scale turbulence in flow round an axisymmetric obstacle is derived. General formulae for velocity covariances and Eulerian time scales are obtained and are evaluated for the particular case of flow round a sphere. The large-scale limit for this flow is also discussed.

Description: A change in thermal conductivity associated with melting or solidification can have a profound influence on the isotherms near the solidification interface if the material is being directionally solidified in an ampoule whose walls carry a substantial portion of the heat. This analysis was prompted by a recent discovery that the thermal conductivity of Hg(1-x)CD(x)Te increased dramatically as the material is heated above the solidus curve. An illustrative example is shown in which the sample is approximated as an infinite cylinder with constant but diffferent thermal properties in the solid and melt. The boundary conditions are fixed on the surface by a conductive ampoule in a two-zone Bridgman furnace with an adiabatic region separating the two zones. The effect of the adiabatic zone in this case is to intensify the curvature of the interface rather than to lessen it.

Description: A cooled porous region has a plane surface exposed to a specified spatially varying heat flux. The coolant leaves the region through this surface, and it is desired to control the flow distribution to maintain a specified uniform surface temperature. This is accomplished by having the coolant entrance surface shaped to provide in the region the necessary variation of path length and, hence, flow resistance. The surface shape at the coolant entrance is found by solving a Cauchy boundary value problem. An exact solution is obtained that will deal with a wide variety of heating distributions for both two- and three-dimensional shapes.

Description: A set of three-dimensional flow-field data for the region around a cylinder impulsively spun-up from rest was derived with a numerical model based on the Navier-Stokes equations. Laser-Doppler anemometer data in the azimuthal direction was employed to test the model predictions, and data was developed for a flowfield with Ekman numbers from 9.18/1,000,000 to 9.18/10,000. The contributions of inviscid and viscous terms were determined as functions of radius and time. It was found that immediately after start-up viscous diffusion is the dominant factor, which is replaced by nonlinear radial advection. The Coriolis force dominates in the later stages of spin-up. The inward radial flow is a maximum near the front, where the vertical velocity is small, but features strong radial gradients, as it does at the edge of the Ekman layer.

Description: Ground-wind environments for Monte Carlo simulations of Space Shuttle liftoff at KSC are developed. Input parameters include randomly selected 18.3-m-altitude peak wind speed (from which mean wind profile and turbulence intensity are calculated), randomly selected mean wind direction, and longitudinal and lateral turbulence components obtained from the Shuttle-simulation turbulence tapes (SSTT: Tatom et al., 1982). The steps in the simulation of turbulence time histories and horizontal wind fields are listed. It is found that separate statistical analysis of each hour-season pair, applying data on the time fraction of occurrence of peak winds and wind directions at KSC, will be necessary to interpret simulation results consistently.

Description: The method of matched asymptotic expansions is used to study the generation of Tollmien-Schlichting waves by free-stream disturbances incident on a flat-plate boundary layer. Near the leading edge, the motion is governed by the unsteady boundary-layer equation, while farther downstream it is governed (to lowest order) by the Orr-Sommerfeld equation with slowly varying coefficients. It is shown that there is an overlap domain where the Tollmien-Schlichting wave solutions to the Orr-Sommerfeld equation and appropriate asymptotic solutions of the unsteady boundary-layer equation match, in the matched-asymptotic-expansion sense. The analysis explains how long-wavelength free-stream disturbances can generate Tollmien-Schlichting waves of much shorter wavelength. It also leads to a set of scaling laws for the asymptotic structure of the unsteady boundary layer.

Description: The results of a first order perfect gas correction for the effects of the boundary layer formation within expansion tubes with nozzles are presented. The analytical model developed to describe the boundary layer formation within the expansion tube and an expansion nozzle located at the end of the acceleration tube is based on the Karman integral equations. The results of this analytical model are compared with experimental data from an expansion diffuser. The model provides a useful tool for the preliminary design of nozzles for such facilities.