ALBERT

Description: The exchange of models is one of the most serious problems currently encountered in the practice of spacecraft thermal analysis. Essentially, the problem originates in the diversity of computing environments that are used across different sites, and the consequent proliferation of native tool formats. Furthermore, increasing pressure to reduce the development's life cycle time has originated a growing interest in the so-called spacecraft concurrent engineering. In this context, the realization of the interdependencies between different disciplines and the proper communication between them become critical issues. The use of a neutral format represents a step forward in addressing these problems. Such a means of communication is adopted by consensus. A neutral format is not directly tied to any specific tool and it is kept under stringent change control. Currently, most of the groups promoting exchange formats are contributing with their experience to STEP, the Standard for Exchange of Product Model Data, which is being developed under the auspices of the International Standards Organization (ISO 10303). This paper presents the different efforts made in Europe to provide the spacecraft thermal analysis community with a Thermal Neutral Format (TNF) based on STEP. Following an introduction with some background information, the paper presents the characteristics of the STEP standard. Later, the first efforts to produce a STEP Spacecraft Thermal Application Protocol are described. Finally, the paper presents the currently harmonized European activities that follow up and extend earlier work on the area.

Description: Using global interpolation functions (GIF's) boundary element solutions are obtained for two-dimensional laminar flows. Two schemes are proposed for handling the convective terms. The first treats convection as a forcing function, and converts the flow equations to pseudo-Poisson equations. In the second scheme, some convective effect is incorporated into the fundamental solution used in constructing the pertinent integral equations. The lid-driven cavity flow is selected as the benchmark problem.

Description: A two dimensional finite volume method is used to predict the film coefficients in the transitional flow region (laminar or turbulent) for the radiator panel tubes. The code used to perform this analysis is CAST (Computer Aided Simulation of Turbulent Flows). The information gathered from this code is then used to augment a Sinda85 model that predicts overall performance of the radiator. A final comparison is drawn between the results generated with a Sinda85 model using the Sinda85 provided transition region heat transfer correlations and the Sinda85 model using the CAST generated data.

Description: This paper presents a probabilistic one-dimensional finite element model for heat transfer processes in porous heat exchangers. The Galerkin approach is used to develop the finite element matrices. Some of the submatrices are asymmetric due to the presence of the flow term. The Neumann expansion is used to write the temperature distribution as a series of random variables, and the expectation operator is applied to obtain the mean and deviation statistics. To demonstrate the feasibility of the formulation, a one-dimensional model of heat transfer phenomenon in superfluid flow through a porous media is considered. Results of this formulation agree well with the Monte-Carlo simulations and the analytical solutions. Although the numerical experiments are confined to parametric random variables, a formulation is presented to account for the random spatial variations.

Description: This paper presents a numerical scheme, based on the finite element method, to solve strongly coupled fluid flow and heat transfer problems. The surface radiation effect for gray, diffuse and isothermal surfaces is considered. A procedure for obtaining the view factors between the radiating surfaces is discussed. The overall solution strategy is verified by comparing the available results with those obtained using this approach. An analysis of a thermosyphon is undertaken and the effect of considering the surface radiation is clearly explained.

Description: An algorithm has been developed for the forced convective diffusion-reaction problem for convection inside and outside a droplet by a recirculating flow field hydrodynamically coupled at the droplet interface with an external flow field that at infinity becomes a uniform streaming flow. The concentration field inside the droplet is likewise coupled with that outside by boundary conditions at the interface. A chemical reaction can take place either inside or outside the droplet or reactions can take place in both phases. The algorithm has been implemented and results are shown here for the case of no reaction and for the case of an external first order reaction, both for unsteady behavior. For pure interphase mass transfer, concentration isocontours, local and average Sherwood numbers, and average droplet concentrations have been obtained as a function of the physical properties and external flow field. For mass transfer enhanced by an external reaction, in addition to the above forms of results, we present the enhancement factor, with the results now also depending upon the (dimensionless) rate of reaction.

Description: The critical state of vortex cores downstream of vortex breakdown has been studied. Base vortical flows were computed using the Reynolds-averaged, axisymmetric Navier-Stokes equations. Standard K - epsilon, RNG and second-order Reynolds stress models were employed. Results indicate that the return to supercriticality is highly dependent on the turbulence model. The K - epsilon model predicted a rapid return of the vortex to supercritical conditions, the location of which showed little sensitivity to changes in the swirl ratio. The Reynolds stress model predicted that the vortex remains subcritical to the end of the domain for each of the swirl ratios employed, and provided results in qualitative agreement with experimental work. The RNG model produced intermediate results, with a downstream movement in the critical location with increasing swirl. Calculations for which area reductions were introduced at the exit in a subcritical flow were also performed using the Reynolds stress model. The structure of the resulting recirculation zone was altered significantly. However, when area reductions were employed within supercritical flows as predicted using the two-equation models, no significant influence on the recirculation zone was noted.

Description: The view factors which are used in diffuse-gray radiation enclosure calculations are often computed by approximate numerical integrations. These approximately calculated view factors will usually not satisfy the important physical constraints of reciprocity and closure. In this paper several view-factor rectification algorithms are reviewed and a rectification algorithm based on a least-squares numerical filtering scheme is proposed with both weighted and unweighted classes. A Monte-Carlo investigation is undertaken to study the propagation of view-factor and surface-area uncertainties into the heat transfer results of the diffuse-gray enclosure calculations. It is found that the weighted least-squares algorithm is vastly superior to the other rectification schemes for the reduction of the heat-flux sensitivities to view-factor uncertainties. In a sample problem, which has proven to be very sensitive to uncertainties in view factor, the heat transfer calculations with weighted least-squares rectified view factors are very good with an original view-factor matrix computed to only one-digit accuracy. All of the algorithms had roughly equivalent effects on the reduction in sensitivity to area uncertainty in this case study.

Description: A transient, one-dimensional numerical code is developed to model the liquid motion in an axial groove with square cross section. Axial variation in liquid level, shear stress and heat transfer between the groove wall and the liquid, evaporation and transient body forces are accounted for in the model. Dryout and rewet of the groove are allowed; the front location is determined numerically using conservation of mass and linear extrapolation. Several numerical test results are presented and discussed.

Description: A higher-order finite-difference technique is developed to calculate the developing-flow field of steady incompressible laminar flows in the entrance regions of circular pipes. Navier-Stokes equations governing the motion of such a flow field are solved by using this new finite-difference scheme. This new technique can increase the accuracy of the finite-difference approximation, while also providing the option of using unevenly spaced clustered nodes for computation such that relatively fine grids can be adopted for regions with large velocity gradients. The velocity profile at the entrance of the pipe is assumed to be uniform for the computation. The velocity distribution and the surface pressure drop of the developing flow then are calculated and compared to existing experimental measurements reported in the literature. Computational results obtained are found to be in good agreement with existing experimental correlations and therefore, the reliability of the new technique has been successfully tested.

Description: The Systems Improved Numerical Differencing Analyzer and Fluid Integrator (SINDA/FLUINT) code has often been used to determine the transient and steady-state response of various thermal and fluid flow networks. While this code is an often used design and analysis tool, the validation of this program has been limited to a few simple studies. For the current study, the SINDA/FLUINT code was compared to four different analytical solutions. The thermal analyzer portion of the code (conduction and radiative heat transfer, SINDA portion) was first compared to two separate solutions. The first comparison examined a semi-infinite slab with a periodic surface temperature boundary condition. Next, a small, uniform temperature object (lumped capacitance) was allowed to radiate to a fixed temperature sink. The fluid portion of the code (FLUINT) was also compared to two different analytical solutions. The first study examined a tank filling process by an ideal gas in which there is both control volume work and heat transfer. The final comparison considered the flow in a pipe joining two infinite reservoirs of pressure. The results of all these studies showed that for the situations examined here, the SINDA/FLUINT code was able to match the results of the analytical solutions.

Description: An IBM Personal Computer (PC) version of the Groove Analysis program (GAP) was developed to predict the steady state heat transport capability of an axially grooved heat pipe for a specified groove geometry and working fluid. In the model, the capillary limit is determined by the numerical solution of the differential equation for momentum conservation with the appropriate boundary conditions. This governing equation accounts for the hydrodynamic losses due to friction in liquid and vapor flows and due to liquid/vapor shear interaction. Back-pumping in both 0-g and 1-g is accounted for in the boundary condition at the condenser end. Slug formation in 0-g and puddle flow in 1-g are also considered in the model. At the user's discretion, the code will perform the analysis for various fluid inventories (undercharge, nominal charge, overcharge, or a fixed fluid charge) and heat pipe elevations. GAP will also calculate the minimum required heat pipe wall thickness for pressure containment at design temperatures that are greater than or lower than the critical temperature of the working fluid. This paper discusses the theory behind the development of the GAP model. It also presents the many useful and powerful capabilities of the model. Furthermore, a correlation of flight test performance data and the predictions using GAP are presented and discussed.

Description: The paper presents the numerical solution of heat and mass transfer during cross-flow (orthogonal) mixed convection. In this class of flow, a buoyancy-driven transport in the vertical direction and a forced convective flow in the horizontal direction results in a three-dimensional boundary layer structure adjacent to the plate. The rates of heat and mass transfer are determined by a combined influence of the two transport processes. The equations for the conservation of mass, momentum, energy, and species concentration were solved along with appropriate boundary conditions to determine the distributions of velocity components, temperature, and concentration across the thickness of the boundary layer at different locations on the plate. Results were expressed in dimensionless form using Reynolds number, Richardson number for heat transfer, Richardson number for mass transfer, Prandtl number, and Schmidt number as parameters. It was found that the transport is dominated by buoyancy at smaller vertical locations and at larger distances away from the forced convection leading edge. Effects of forced convection appeared to be very strong at smaller horizontal distances from the leading edge. The cross stream forced convection enhanced the rate of heat and mass transfer by a very significant amount.

Description: The paper presents a detailed theoretical analysis of the process of gas absorption to a thin liquid film adjacent to a horizontal rotating disk. The film is formed by the impingement of a controlled liquid jet at the center of the disk and subsequent radial spreading of liquid along the disk. The chemical reaction between the gas and the liquid film can be expressed as a zero-order homogeneous reaction. The process was modeled by establishing equations for the conservation of mass, momentum, and species concentration and solving them analytically. A scaling analysis was used to determine dominant transport processes. Appropriate boundary conditions were used to solve these equations to develop expressions for the local concentration of gas across the thickness of the film and distributions of film height, bulk concentration, and Sherwood number along the radius of the disk. The partial differential equation for species concentration was solved using the separation of variables technique along with the Duhamel's theorem and the final analytical solution was expressed using confluent hypergeometric functions. Tables for eigenvalues and eigenfunctions are presented for a number of reaction rate constants. A parametric study was performed using Reynolds number, Ekman number, and dimensionless reaction rate as parameters. At all radial locations, Sherwood number increased with Reynolds number (flow rate) as well as Ekman number (rate of rotation). The enhancement of mass transfer due to chemical reaction was found to be small when compared to the case of no reaction (pure absorption), but the enhancement factor was very significant when compared to pure absorption in a stagnant liquid film. The zero-order reaction processes considered in the present investigation included the absorption of oxygen in aqueous alkaline solutions of sodiumdithionite and rhodium complex catalyzed carbonylation of methanol. Present analytical results were compared to previous theoretical results for limiting conditions, and were found to have very good agreement.

Description: Forced convective diffusion-reaction is considered for viscous axisymmetric extensional convecting velocity in the neighborhood of a sphere. For Peclet numbers in the range 0.1 less than or equal to Pe less than or equal to 500 and for Damkohler numbers increasing with increasing Pe but in the overall range 0.02 less than or equal to Da less than or equal to 10, average and local Sherwood numbers have been computed. By introducing the eigenfunction expansion c(r, Theta) = Sum of c(n)(r)P(n)(cos Theta) into the forced convective diffusion equation for the concentration of a chemical species undergoing a first order homogeneous reaction and by using properties of the Legendre functions Pn(cos Theta), the variable coefficient PDE can be reduced to a system of N + 1 second order ODEs for the radial functions c(sub n)(r), n = 0, 1, 2,..., N. The adaptive grid algorithm of Pereyra and Lentini can be used to solve the corresponding 2(N + 1) first order differential equations as a two-point boundary value problem on 1 less than or equal to r less than or equal to r(sub infinity). Convergence of the expansion for a specific value of N can thus be established and provides 'spectral' behavior as well as the full concentration field c(r, Theta).

Description: Annotation is a key activity of data analysis. However, current systems for data analysis focus almost exclusively on visualization. We propose a system which integrates annotations into a visualization system. Annotations are embedded in 3D data space, using the Post-it metaphor. This embedding allows contextual-based information storage and retrieval, and facilitates information sharing in collaborative environments. We provide a traditional database filter and a Magic Lens filter to create specialized views of the data. The system has been customized for fluid flow applications, with features which allow users to store parameters of visualization tools and sketch 3D volumes.

Description: Fluctuating surface pressure measurements have been made to investigate the effectiveness of boundary layer separators (BLS's) in reducing the fluctuating pressure loads produced by separated shock wave turbulent boundary layer interactions. Measurements have been made under unswept and swept compression corner interactions in a Mach 5 flow. BLS's fix the separation location and eliminate the large-amplitude, low-frequency fluctuating pressure loads upstream of the compression corners. The loads on the unswept compression corner face are reduced by as much as 59%. The BLS's also shift the mean pressure distribution on the unswept corner face in the streamwise direction. Results show that the loads on the corner face vary with the BLS height and the distance between the BLS and the compression corner. Suggestions for the optimum placement and the use of the BLS's are also made.

Description: Fluctuating wall pressure measurements have been made in a separated shock wave/turbulent boundary layer interaction produced by an unswept compression corner in a Mach 5 flow. Wheeler doublet vortex generators were placed 15.8 boundary layer thicknesses upstream of the corner to study their effect on the fluctuating pressure loads produced by the translating separation shock. The vortex generators produced significant three-dimensionality in an otherwise two-dimensional interaction. They reduced the upstream influence and the length of the region of shock motion by 60% and 64%, respectively, decreased the maximum wall pressure rms by 23%, and shifted the fluctuations to a higher frequency band. The maximum fraction of energy in the 100-500 Hz frequency band is decreased by 11%. These changes are due to a fuller boundary layer profile, a weaker separation shock, and increased boundary layer turbulence causing increased separation shock jitter.

Description: A semianalytical solution is presented for the unsteady pressure field of a vortical gust interacting with a flat-plate airfoil in subsonic flow. The solution will serve as a benchmark for evaluating the accuracy and efficiency of time dependent numerical schemes. The specific case considered corresponds to the ICASE benchmark problem number 6. The results are compared with those of asymptotic theories for high frequency and show excellent agreement.

Description: A numerical method has been developed in order to address aeroacoustic problems modeled by the linearized Euler equations. A weak formulation of the equations leads to a time-dependent equation for the test functions. The basic solver being one dimensional, two dimensional problems are handled by directional splitting. This method shows low dissipation and dispersion errors.

Description: The underlying bases and developments in two techniques of detailed turbulence modeling are described where the flow is treated in the Eulerian sense, and one technique where the Lagrangian motions of vortices are followed. First, a technique is described for solving the single-point statistically averaged conservation equations. The Reynolds stresses that appear in these equations are evaluated by solving supplemental differential equations which contain terms that are modeled. A sequence of increasingly complex, but also increasingly general, modeling equations is described and computations based on these equations are compared with experimental data. The hierarchy of models described terminates with equations for the individual components of the Reynolds stress tensor. The second Eulerian technique approach to turbulence modeling is the direct numerical simulation of turbulent fields. In this approach, all three dimensional eddies between a predetermined range of sizes are computed in time within a specified volume of flow. Present day computers require a tradeoff between the size of the volume that can be considered and the degree of resolution of the turbulent eddies. Techniques of modeling the smallest eddies are described that permit enlarging the volume, or Reynolds number, that can be considered.

Description: The paper describes a theoretical investigation of the stability under gravitational and surface forces of a liquid in a circular cylindrical container with a concave spheroidal bottom for the case in which the volume of the liquid is sufficiently small so that the bottom is not covered completely. The gravitational field is assumed to be directed along the symmetry axis of the container, and for a specific container shape the critical Bond number is calculated as a function of liquid volume for contact angles of gamma = 0, 1, 2, and 4 deg. For gamma = 0 deg, some critical equilibrium configurations and corresponding perturbation modes are presented.

Description: The presented review is concerned with the problem of calculating compressible viscous flows. Basic numerical considerations and problems associated with calculating viscous flows are examined and current numerical approaches toward the solution of the Navier-Stokes equations are discussed. It is pointed out that the numerical solution of the full time-dependent equations for turbulent flow is not practical with present computers. Therefore, turbulence effects must be accounted for by modeling. Developments related to turbulence modeling are described. In connection with a discussion of numerical methods for solving viscous flow equations, attention is given to numerical domains of dependence of typical explicit and implicit methods, the diffusion problem, the convection-diffusion problem, and the split-hybrid method.

Description: Steady state solutions to two time dependent partial differential systems have been obtained by the Method of Lines (MOL) and compared to those obtained by efficient standard finite difference methods: (1) Burger's equation over a finite space domain by a forward time central space explicit method, and (2) the stream function - vorticity form of viscous incompressible fluid flow in a square cavity by an alternating direction implicit (ADI) method. The standard techniques were far more computationally efficient when applicable. In the second example, converged solutions at very high Reynolds numbers were obtained by MOL, whereas solution by ADI was either unattainable or impractical. With regard to 'set up' time, solution by MOL is an attractive alternative to techniques with complicated algorithms, as much of the programming difficulty is eliminated.

Description: Experimental results show conclusively that the presence of a small quantity of a noncondensable gas (NCG) mixed with the working fluid has a considerable effect on the condensation process in a rotating heat pipe. The temperature distribution in the condenser shows the blanketing effect of the NCG and the ratio of the molecular weight of the working fluid to that of the NCG has a very definite effect on the shape of this distribution. Some of the effects are quite similar to the well-established data on stationary heat pipes.

Description: Experimental evidence shows the importance of external boundary conditions on the overall performance of a rotating heat pipe condenser. Data are presented for the boundary conditions of constant heat flux and constant wall temperature for rotating heat pipes containing either pure vapor or a mixture of vapor and noncondensable gas as working fluid.

Description: The turbulence generated by random entropy fluctuations in an accelerating stream is analyzed. The results are obtained by using rapid distortion theory together with a high frequency solution of a previously developed wave equation that governs the small-amplitude unsteady vortical and entropic motion on steady potential flows (Goldstein, 1978). Simple results are obtained for the case of symmetric contraction, expansion or combination of the two. It is shown that the energy of the entropy-generated turbulence increases more rapidly with the contraction ratio of a subsonic flow than that of any imposed upstream turbulence. This result indicates that the entropy-generated turbulence may be more significant than the hydrodynamically generated turbulence in the turbine stages of aircraft engines.

Description: On the basis of this investigation of the high-temperature behavior of polytetrafluoroethylene (PTFE), the transient one-dimensional ablation of PTFE has been developed by taking into account the optical transmittance of both the amorphous zone and the crystalline zone of PTFE layer. Results show that although the exposed surface receded at an apparently steady state, both the internal temperature and the thickness of the gel layer increase continuously due to the internal absorption of radiation.

Description: It is shown that the pressure and velocity fluctuations of the unsteady motion on a transversely sheared mean flow can be expressed entirely in terms of the derivatives of two potential functions. One of these is a convected quantity that can be specified as a boundary condition and is related to a transverse component of the upstream velocity field. The other can be determined by solving an inhomogeneous wave equation whose source term is also a convected quantity that can be specified as a boundary condition in any given problem. The general theory is used to study the interaction of an unsteady flow with a semi-infinite plate embedded in a shear layer. The acoustic field produced by this interaction is calculated in the limits of low and high frequency. The results are compared with experimental one-third octave sound pressure level radiation patterns. The agreement is found to be excellent, especially in the low frequency range, where the mean-flow and convective effects are shown to have a strong influence on the directivity of the sound.

Description: Surface water waves generated by surface and near surface point explosions are calculated. Taking the impulse distribution imparted at the water surface by the explosion as the overriding mechanism for transferring energy of the explosive to surface wave motion, the linearized theory of Kranzer and Keller is used to obtain the wave displacement in the far field. The impulse distribution is obtained by integrating the pressure wave over an appropriate time interval on a horizontal surface just beneath the undisturbed water surface. For surface explosions, a modified form of the similarity method first used by Collins and Holt is used to obtain the flow field. In the case of submerged explosions, the flow field is estimated by making necessary modifications to Sedov's similarity solution to account for the venting that accompanies the interaction of the leading (blast) wave with the ocean surface. Surface waves generated by a charge at six depths of placement (0.15 m, 0.30 m, 0.61 m, 0.91 m, 1.37 m, 3.05 m) are considered in addition to surface explosions. The results seem to support the existence of an upper critical depth phenomenon (of the type already established for chemical explosions) for point (nuclear) explosions.

Description: Compressible Navier-Stokes equations for quasi-one-dimensional flow in a converging-diverging nozzle have been solved using Stetter's three-step predictor-corrector technique. Particular emphasis is given to the minimum iterative step feature for steady-state solutions. It is found that for the nonviscous-dominated case, Stetter's method attains the steady-state solution in the fewest steps when compared to four other currently used techniques.

Description: Predictions of mixing length by Pletcher's (1976) method (using a two-layer eddy viscosity model for a turbulent boundary layer at low Reynolds numbers) are compared to some values derived from turbulent boundary-layer profiles by other authors. The model is incorporated into a finite-difference scheme, to accurately predict low Reynolds number skin friction in supersonic flow.

Description: The use of hot-wire anemometry for obtaining fluctuating data in transonic flows has been evaluated. From hot-wire heat loss correlations based on previous transonic data, the sensitivity coefficients for velocity, density, and total temperature fluctuations have been calculated for a wide range of test conditions and sensor parameters. For sensor Reynolds number greater than 20 and high sensor overheat ratios, the velocity sensitivity remains independent of Mach number and equal to the density sensitivity. These conditions were verified by comparisons of predicted sensitivities with those from recent direct calibrations in transonic flows. Based on these results, techniques are presented to obtain meaningful measurements of fluctuating velocity, density, and Reynolds shear stress using hot-wire and hot-film anemometers. Example of these measurements are presented for two transonic boundary layers.

Description: The accretion of rotating gas on to a gravitating object is investigated by means of a perturbation to the spherically symmetrical flow. An expression is found for the correction to the accretion rate due to rotation of the gas in terms of the boundary conditions far from the object. In the case of accretion from a cloud with uniform angular velocity, the gas is accreted preferentially from a 'jet' near the axis of rotation. The angular-momentum distribution in the cloud can be altered by the propagation of inertia waves.

Description: An iterative method for numerically solving the time independent Navier-Stokes equations for viscous compressible flows is presented. The method is based upon partial application of the Gauss-Seidel principle in block form to the systems of the nonlinear algebraic equations which arise in construction of finite element (Galerkin) models approximating solutions of fluid dynamic problems. The continuous cubic element on triangles is employed for function approximation. Computational results for a free shear flow at Re = 1000 indicate significant achievement of economy in iterative convergence rate over finite element and finite difference models which employ the customary time dependent equations and symptotic time marching procedure to steady solution. Numerical results are in excellent agreement with those obtained for the same test problem employing time marching finite element and finite difference solution techniques.

Description: An analysis of condensation problems in rotating heat pipes containing vapors with different concentrations of non-condensable gases is given. In situations such as this, temperature and concentration gradients are set up in the vapor-gas mixture. There is a transport of mass due to temperature gradients accompanied by an energy transport phenomena due to a concentration gradient. A Nusselt type analysis is not suited to this type of problem; however, a boundary layer type approach has successfully been used to analyze stationary condensation systems with non-condensable gases present. The present boundary layer analysis is presented for condensation processes on the inside of a rotating heat pipe in the presence of non-condensable gases.

Description: The aims of the experiment are outlined. Flight experiments included in this program were provided by NASA, Goddard Space Flight Center, ESA (European Space Agency), the German Ministry of Technology, Hughes Aircraft Company and NASA, Ames Research Center.

Description: A theoretical model was derived for laminar film condensation on the inside of a rotating, truncated cone which includes the effects of vapor shear and vapor pressure drop. Results are compared to those of previous investigations. Experimental data are presented for rotational speeds of 700, 1,400, 2,100, and 2,800 rpm using water, ethyl alcohol, and Freon 113 as working fluids. Agreement between theory and experiment is within + or - 20 percent.

Description: Preliminary compliant wall skin friction test results obtained in a low-turbulence pressure tunnel are reported. Compliant surface skins consisted of 0.0025 cm thick mylar, stretched under tension and area-bonded or longitudinally strip-bonded with silicone rubber adhesive to polyurethane foam. Mean velocity and fluctuating survey data were obtained with a single slanted hot wire. Mean velocity profiles for the area-bonded mylar surface skins indicated up to a 20% reduction in boundary-layer thickness (and lower momentum thickness) over rigid surfaces. This reduction in boundary-layer thickness indicates that a drag reduction occurred. In addition, a 16% reduction in wall shear was evident for the mylar/compressor foam compliant surface.

Description: Experiments investigating the stability characteristics of a single-phase free convection loop are reported. Results of the study confirm the contention made by previous workers that instabilities near the thermodynamic critical point can occur for ordinary fluids as well as those with unusual behavior in the near-critical region. Such a claim runs counter to traditional beliefs, but it is supported by the observation of such instabilities for water at atmospheric pressure and moderate temperatures in the present work.

Description: The essential ingredients and the effectiveness of several levels of turbulent-flow partial differential equation models are considered. Zero-equation models use only the partial differential equation (pde) for the mean velocity field and do not employ turbulence pde's. One-equation models make use of an additional pde relating to the turbulence velocity scale. Attention is also given to two-equation models, stress-equation models, and large eddy simulations. Large-eddy simulations are concerned with a three-dimensional time-dependent numerical computation of the large-scale turbulence.

Description: The paper determines the effect of various available drag coefficient equations on particle velocity calculations for typical two phase flows encountered in supersonic and turbulent laser velocimeter applications. The predictions of the particle drag coefficient equations are compared with experimental sphere drag data. For the laser velocimeter applications, the relative Mach number less than 2 and the relative Reynolds number less than 200 are of particular importance.

Description: The paper studies the effect of low Reynolds number in high-speed turbulent boundary layers on variations of mixing length. Boundary layers downstream of natural transition on plates, cones and cylinders, and boundary layers on nozzle walls without laminarization-retransition are considered. The problem of whether low Reynolds number amplification of shear stress is a result of transitional flow structure is considered. It is concluded that a knowledge of low Reynolds number boundary layer transition may be relevant to the design of high-speed vehicles.

Description: A previous analysis of the acoustic radiation from multipole sources is extended to include additional components of the dipole and quadrupole sources. It is found that, unlike the components of the sources considered in the previous paper, the exponent of the Doppler factor now depends on the location of the sources within the jet.

Description: A possible alternative explanation is proposed for compliant wall drag reductions measured in previous investigations. Standing waves were observed to form on the surfaces of compliant wall models in air with water substrates as the freestream velocity was increased from 15 to 30 m/s. These waves resembled sine waves with half of the wave protruding over the upstream portion of the model and the other half being recessed over the downstream end of the model. These data coupled with results of recent drag reduction experiments suggest that standing waves could have caused a shift in the model center of gravity creating a bending moment that was interpreted as a reduction in the skin friction drag.

Description: A general analysis is presented of the steady nondissipative constant-property laminar boundary layer flow over a two-dimensional body of uniform surface heat flux situated in an infinite ambient fluid of undisturbed temperature. The analysis is then applied to a long horizontal circular cylinder. Numerical solutions to the universal functions associated with the first two terms in the derived series are given for Prandtl numbers 0.7, 1, 3, 5, 7, and 10. The results are compared with those obtained by Koh (1964) whose method is patterned after the Blasius-Frossling procedure for forced convection flow. The study reveals that Wilks' (1972) analysis concerning the external natural convection about two-dimensional bodies with constant heat flux is in error.

Description: Accurate capturing of discontinuities within compressible flow computations is achieved by coupling a suitable solver with an automatic adaptive mesh algorithm for unstructured triangular meshes. The mesh adaptation procedures developed rely on non-hierarchical dynamical local refinement/derefinement techniques, which hence enable structural optimization as well as geometrical optimization. The methods described are applied for a number of the ICASE test cases are particularly interesting for unsteady flow simulations.

Description: This paper discusses the issues which arise when combining multigrid strategies with adaptive meshing techniques for solving steady-state problems on unstructured meshes. A basic strategy is described, and demonstrated by solving several inviscid and viscous flow cases. Potential inefficiencies in this basic strategy are exposed, and various alternate approaches are discussed, some of which are demonstrated with an example. Although each particular approach exhibits certain advantages, all methods have particular drawbacks, and the formulation of a completely optimal strategy is considered to be an open problem.

Description: The dynamic solution adaptive grid algorithm, DSAGA3D, is extended to automatically adapt 2-D structured multi-block grids, including adaption of the block boundaries. The extension is general, requiring only input data concerning block structure, connectivity, and boundary conditions. Imbedded grid singular points are permitted, but must be prevented from moving in space. Solutions for workshop cases 1 and 2 are obtained on multi-block grids and illustrate both increased resolution of and alignment with the solution. A mesh quality assessment criteria is proposed to determine how well a given mesh resolves and aligns with the solution obtained upon it. The criteria is used to evaluate the grid quality for solutions of workshop case 6 obtained on both static and dynamically adapted grids. The results indicate that this criteria shows promise as a means of evaluating resolution.

Description: This paper describes adaptive grid methods developed specifically for compressible flow computations. The basic flow solver is a finite-volume implementation of Roe's flux difference splitting scheme or arbitrarily moving unstructured triangular meshes. The grid adaptation is performed according to geometric and flow requirements. Some results are included to illustrate the potential of the methodology.

Description: Solution-adaptive methods based on cutting bodies out of Cartesian grids are gaining popularity now that the ways of circumventing the accuracy problems associated with small cut cells have been developed. Researchers are applying Cartesian-based schemes to a broad class of problems now, and, although there is still development work to be done, it is becoming clearer which problems are best suited to the approach (and which are not). The purpose of this paper is to give a candid assessment, based on applying Cartesian schemes to a variety of problems, of the strengths and weaknesses of the approach as it is currently implemented.

Description: A three dimensional finite volume scheme based on hybrid grids containing both tetrahedral and hexahedral cells is presented. The application to hybrid grids offers the possibility to combine the flexibility of tetrahedral meshes with the accuracy of hexahedral grids. An algorithm to compute a dual mesh for the entire computational domain was developed. The dual mesh technique guarantees conservation in the whole flow field even at interfaces between hexahedral and tetrahedral domains and enables the employment of an accurate upwind flow solver. The hybrid mesh can be adapted to the solution by dividing cells in areas of insufficient resolution. The method is tested on different viscous and inviscid cases for hypersonic, transonic and subsonic flows.

Description: A method for generating high quality unstructured triangular grids for high Reynolds number Navier-Stokes calculations about complex geometries is described. Careful attention is paid in the mesh generation process to resolving efficiently the disparate length scales which arise in these flows. First the surface mesh is constructed in a way which ensures that the geometry is faithfully represented. The volume mesh generation then proceeds in two phases thus allowing the viscous and inviscid regions of the flow to be meshed optimally. A solution-adaptive remeshing procedure which allows the mesh to adapt itself to flow features is also described. The procedure for tracking wakes and refinement criteria appropriate for shock detection are described. Although at present it has only been implemented in two dimensions, the grid generation process has been designed with the extension to three dimensions in mind. An implicit, higher-order, upwind method is also presented for computing compressible turbulent flows on these meshes. Two recently developed one-equation turbulence models have been implemented to simulate the effects of the fluid turbulence. Results for flow about a RAE 2822 airfoil and a Douglas three-element airfoil are presented which clearly show the improved resolution obtainable.

Description: This paper presents an hp-adaptive discontinuous Galerkin method for linear hyperbolic conservation laws. A priori and a posteriori error estimates are derived in mesh-dependent norms which reflect the dependence of the approximate solution on the element size (h) and the degree (p) of the local polynomial approximation. The a posteriori error estimate, based on the element residual method, provides bounds on the actual global error in the approximate solution. The adaptive strategy is designed to deliver an approximate solution with the specified level of error in three steps. The a posteriori estimate is used to assess the accuracy of a given approximate solution and the a priori estimate is used to predict the mesh refinements and polynomial enrichment needed to deliver the desired solution. Numerical examples demonstrate the reliability of the a posteriori error estimates and the effectiveness of the hp-adaptive strategy.

Description: Fluid dynamic research with the objective of developing new and improved technology in both test facility concepts and test techniques is being reported. A summary of efforts and results thus far obtained in four areas is presented. The four area are: (1) the use of heavy gases to obtain high Reynolds numbers at transonic speeds: (2) high Reynolds number tests of the C-141A wing configuration; (3) performance and flow quality of the pilot injector driven wind tunnel; and (4) integration time required to extract accurate static and dynamic data from tests in transonic wind tunnels. Some of the principal conclusions relative to each of the four areas are: (1) Initial attempts to apply analytical corrections to test results using gases with gamma other than 1.4 to simulate conditions in air show promise but need significant improvement; (2) for the C-141A configuration, no Reynolds number less than the full scale flight value provides an accurate simulation of the full scale flow; (3) high ratios of tunnel mass flow rate to injection mass flow rate and high flow quality can be obtained in an injector driven transonic wind tunnel; and (4) integration times of 0.5 to 1.0 sec may be required for static force and pressure tests, respectively, at some transonic test conditions in order to obtain the required data accuracy.

Description: As part of a special international effort, three nozzles were designed and tested on single nacelle models in wind tunnels of several nations belonging to the North Atlantic Treaty Organization. All three of these nozzles were investigated in the Langley 16-foot transonic wind tunnel at the National Aeronautics and Space Administration's Langley Research Center. Langley Research Center also contributed theoretical calculations of the jet plume boundary and afterbody pressures. The calculations were obtained using an iterative solution which combined the inviscid Douglas Neumann method for the external flow with the method of characteristics for the flow in the jet plume. For the investigation, the nozzles were mounted on a single nacelle model 15.24 centimeters in diameter and 162.56 centimeters long. Tests were made at free stream Mach number from 0.4 to 1.2, and at Reynolds numbers per meter from 7.38 million to 13.78 million depending on the Mach number. Four types of data were recorded: afterbody pressure data, afterbody force data, model boundary layer data, and tunnel wall pressure data. The ratio of jet total pressure to free stream static pressure ranged up to 8.5. A description of the wind tunnel, model, and test procedure is included.

Description: The time-splitting explicit numerical method of MacCormack is applied to separated turbulent boundary layer flow problems. Modifications of this basic method are developed to counter difficulties associated with complicated geometry and severe numerical resolution requirements of turbulence model equations. The accuracy of solutions is investigated by comparison with exact solutions for several simple cases. Procedures are developed for modifying the basic method to improve the accuracy. Numerical solutions of high-Reynolds-number separated flows over an airfoil and shock-separated flows over a flat plate are obtained. A simple mixing length model of turbulence is used for the transonic flow past an airfoil. A nonorthogonal mesh of arbitrary configuration facilitates the description of the flow field. For the simpler geometry associated with the flat plate, a rectangular mesh is used, and solutions are obtained based on a two-equation differential model of turbulence.

Description: Various types of series solutions for predicting laminar, free-convection boundary-layer heat transfer over both isothermal and nonisothermal boundaries are reviewed. The methods include finite difference, Merk series, Blasius series, and Goertler series. Comparative results are presented for heat transfer over an isothermal, horizontal, elliptical cylinder in both slender and blunt configurations.

Description: Although the Navier-Stokes equations describe most flows of interest in aerodynamics, the inviscid conservation law equations may be used for small regions with viscous forces. Thus, Euler equations and several time-accurate finite difference procedures, explicit and implicit, are discussed. Although implicit techniques require more computational work, they permit larger time steps to be taken without instability. It is noted that the Jacobian matrices for Euler equations in conservation-law form have certain eigenvalue-eigenvector properties which may be used to construct conservative-form coefficient matrices. This reduces the computation time of several implicit and semiimplicit schemes. Extensions of the basic approach to other areas are suggested.

Description: The problem of decaying isotropic turbulence has been studied using a Wiener-Hermite expansion with a renormalized time-dependent base. The theory is largely deductive and uses no modeling approximations. It has been found that many properties of large-Reynolds-number turbulence can be calculated (at least for moderate time) using the moving-base expansion alone. Such properties found are the spectrum shape in the dissipation range, the Kolmogorov constant, and the energy cascade in the inertial subrange. Furthermore, by using a renormalization scheme, it is possible to extend the calculation to larger times and to initial conditions significantly different from the equilibrium form. If the initial spectrum is the Kolmogorov spectrum perturbed with a spike or dip in the inertial subrange, the process proceeds to eliminate the perturbation and relax to the preferred spectrum shape. The turbulence decays with the proper dissipation rate, and several other properties are found to agree with measured data. The theory is also used to calculate the energy transfer and the flatness factor of turbulence.

Description: A steady-state analysis is conducted to examine the basic flow structure of a non-Newtonian fluid in a domain including an inflow region, a contraction region, and an outflow region. A Cartesian grid system is used throughout the entire flow domain, including the contraction region, thus creating an irregular grid cell structure adjacent to the curved boundary. At node points adjacent to the curved boundary symmetry conditions are derived for the different flow variables in order to solve the governing difference equations. Attention is given to the motion and non-Newtonian constitutive equations, the boundary conditions, the numerical modeling of the non-Newtonian equations, the stream function contour lines for the non-Newtonian fluid, the vorticity contour lines for the non-Newtonian fluid, the velocity profile across the contraction, and the shear stress contour lines for the non-Newtonian fluid.

Description: A three-dimensional turbulent boundary layer relaxing behind a transverse hump (shaped as a 30-deg swept 5-ft chord wing-type model) was studied in a low-speed wind tunnel. Data obtained with hot-wire probes showed that the apparent dimensionless velocity profiles in the viscous sublayer prevail universally; evidence for wall similarity in the relaxing flow field was confirmed in the form of a log law. An unusual region of slightly decreasing cross-flow angle was found in collateral regions, and a near-wall noncollateral flow was posited. Streamwise relaxation of the mean flow field was also investigated.

Description: The paper describes turbulence simulation experiments based on the principles of control system theory, that is, the construction of a system characterized by a system function such that upon exciting the system with prescribed noise processes the output of the system is a realization of a random processing the desired statistical attributes of turbulence. An experimental autocorrelation of Jimsphere measurements of wind velocity was approximated to simulate turbulent wind. From the approximate autocorrelation function, the required system function is obtained, and a discrete time system is designed. Another method of simulation is to solve the convolution integral by filter techniques. Other methods include discrete Fourier simulation and self-similar simulation.

Description: The problem of closure in turbulence in the case of two-point correlations resides in the existence of two unknowns E and W, the energy spectrum function and the transfer function, respectively, in the spectrum equation. In the case of weak turbulence, W is negligible. In case of higher correlations, closure can be effective by neglecting the inertia term in the highest order term used. Specifying a certain number of spectra at an initial time is also a way of getting around the closure problem. A simple case of turbulent shear flow is then considered, where two-point correlation equations are used and the velocity is broken into mean and fluctuating components. This yields a differential equation for the energy spectrum, the three terms of which are the energy spectrum, production term and dissipation term. They are plotted for a particular time. Similar analyses and comparisons with experiment are made for pipe and boundary layer flows.

Description: A solution for the two-dimensional and axisymmetric laminar boundary-layer momentum equation of power-law non-Newtonian fluid is presented. The analysis makes use of the Merk-Chao series solution method originally devised for the flow of Newtonian fluid. The universal functions for the leading term in the series are tabulated for n from 0.2 to 2. Equations governing the universal functions associated with the second and the third terms are provided. The solution together with either Lighthill's formula or Chao's formula constitutes a simple yet general procedure for the calculation of wall shear and surface heat transfer rate. The theory was applied to flows over a circular cylinder and a sphere and the results compared with published data.

Description: Evolution of a rotating flow in a body of fluid bounded by a stationary flat surface is discussed. The calculated results show that the radial pressure gradient is substantially reduced in the region close to the surface, so that letting that gradient be independent of distance from the surface would be expected to give only rough or qualitative estimates. However, the reduced rotation near the stationary surface is still large enough to cause an inflow near the surface and to set up a recirculation pattern. The concentration of vorticity by the radial inflow is not great enough to increase the tangential velocities near the center of rotation.

Description: An analysis is described of long, finite-amplitude internal waves in a stratified shear flow. Both regular and singular modes are considered with a nonlinear critical layer employed in the latter case. A three-layer model is used to develop the theory and closed-form expressions are obtained relating the phase speed to the Richardson number, the latter quantity being taken as O(1). The amplitude evolution equation is found to be either the Korteweg-de Vries equation or the Benjamin-Davis-Ono equation depending upon the distance of the more remote boundary from the edge of the shear layer.

Description: Explicit, implicit, and characteristic finite-difference methods are applied to solve model equations representative of the compressible Navier-Stokes equations. An approach is then formulated for solving the Navier-Stokes equation at high Reynolds numbers. The approach has drastically reduced the computation time required to obtain viscous flow solutions. Computational results for shock wave separated flows are presented.

Description: Formal solutions of the static equilibrium equations for the form of the outer surface of a pendent liquid drop are studied. An approach is adopted in which only the one-parameter family determined by vertex height (u sub 0) need be described. Attention is restricted to rotationally symmetric configurations, and all symmetric solutions are characterized for the case where the Lagrange parameter lambda is equal to zero. It is shown that for any u sub 0 the function u(r; u sub 0) can be extended as a parametric solution of a system of equations for all arc lengths, yielding a curve without limit sets or double points, and that the resulting capillary rotation surface spreads out indefinitely away from the axis r = 0. The asymptotic form of the surface in the case of large absolute values of u sub 0 is characterized quantitatively, along with the global structure of all such surfaces.

Description: A numerical study of transient thermal response of a blunt-nosed axisymmetric body made of Teflon is presented using a two-layer thermal model. It is shown that phase change and transverse heat conduction have a considerable effect on the internal temperature field. Comparison of the numerical results with experimental data shows that the single-layer thermal model does not predict the real feature of the thermal field, whereas the results of the two-layer thermal model agree reasonably well with the experiment.

Description: Turbulence decay is calculated by using experimental initial conditions and discarding quadruple-correlation terms in the correlation equations. Agreement with experiment is good only for moderately small times, but there are no perceptible negative spectral energies even at large times.

Description: The theoretical study of flow in biological vascular systems is made very difficult in connection with local and temporal changes of the cross section. Experimental investigations with the aid of model tests are, however, not enough for a solution of the problems, and numerical solutions are more and more employed for a description of such flows. A description is presented of a difference procedure for the solution of the complete Navier-Stokes equations with curvilinear coordinates for three-dimensional flows in containers. The integration of the Navier-Stokes equations for flows in containers with rigid and moving walls is considered.

Description: This paper presents an approximate, but general, analysis for the thermal response behavior of incompressible, constant property, laminar boundary layer flow over a smooth object of arbitrary shape. It encompasses the classic time-dependent Leveque problem as its special case. Comparisons with available data show that, for fluids with Prandtl number of the order of unity or larger, this analytical solution is able to provide reasonably accurate results for most engineering applications. Under certain restrictive conditions, it can also be used to predict the thermal response of a compressible boundary layer flow.

Description: A finite-element thermal analysis procedure for elements with several temperature-dependent thermal parameters is presented. The procedure, based on an application of the Newton-Raphson iteration technique, is formulated by resolving element matrices into component matrices, one component for each thermal parameter. Component conductance matrices are evaluated by assuming constant thermal parameters within an element and are computed once per unit thermal parameter. Significant savings in computer time result from the unit thermal parameter concept. The solution procedure applied to a convectively cooled structure with significantly varying thermal parameters converged in four iterations.

Description: The reported investigation shows that the hot-wire probe induces stable upstream oscillations in a free shear layer, similar to the jet edge tone mechanism. This effect can be significant also in measurements involving large-scale organized structure, conditional sampling, space-time correlation, and convection velocity, when a reference or indicator probe may be used near the origin of the free shear layer. It appears that even in a free shear layer without any wedge, an object in the flow sufficiently downstream can also provide feedback to the flow upstream. A description is given of the edge-tone phenomenon which is observed when a thin slit jet impinges on a plane wedge. Attention is given to the free shear layer tone induced by a hot-wire probe, the free shear layer tone phenomenon, and shear layer tone eigenvalues and eigenfunctions.

Description: A procedure is described for computation of incompressible, steady, two-dimensional flows in fully stalled diffusers with plenum exit. The procedure is successful in predicting pressure distributions and patterns to the accuracy of the data. The procedure employs a zonal model; this maintains close connections between the modeling and the physics thereby providing insight into critical aspects of modeling separated flows. The procedure presented is also convenient for computing unstalled flows in passages with turbulent boundary layers for either direct or indirect design problems. Computing times are well within engineering feasibility. The concepts developed can be extended to other classes of separated flows; some of these extensions have already been completed and are referenced.

Description: A technique is described for the numerical solution of non-axisymmetric flow problems posed in cylindrical coordinates when the z-axis is included in the flowfield. The highlight of the technique is the manner in which the singularities at the centerline are handled. Specifically, the governing flowfield equations at r = 0 are put in a special form by applying L'Hospital's Rule. The required radial derivatives are evaluated using a one-sided, second-order accurate, first-difference. This leads to a smooth, convergent calculation of the flowfield at the centerline. This appears to be the first generally applicable numerical method for avoiding coordinate system singularities in the context of a finite-difference scheme, and could have application to many nonaxisymmetric flows. The technique is illustrated by specific results for the time-dependent flowfield inside an internal combustion engine.

Description: Stereoscopic photography was used to record the three-dimensional motions of fluid elements along a flat-plate turbulent boundary layer flow. The outer region of the boundary layer was dominated by the formation and convection of transverse vortices, formed at the high-speed front between low- and high-speed fluid elements. Vortical motions observed in the wall region consisted of (1) motions of single particles around part of a circle, (2) streamwise vortices, and (3) transverse vortices.

Description: When a point source explosion is initiated at the ocean surface, the shock propagated into the water is reflected at the surface as a centered expansion wave. The solution in the neighborhood of the interaction point is obtained by writing the equations of motion in the appropriate similarity variables and then changing the independent variables to polar coordinates based at the interaction point. From the zero-order solution of the resulting equations the slopes of boundaries at the interaction point are obtained. A first-order perturbation of this solution provides more accurate representation of the flow variables and the curvature of the shock surface near the interaction point.

Description: A statistical description of turbulence in an incompressible fluid obeying the Navier-Stokes equations is proposed, where pressure is regarded as a potential for the interaction between fluid elements. A scaling procedure divides a fluctuation into three ranks representing the three transport processes of macroscopic evolution, transport property, and relaxation. Closure is obtained by relaxation, and a kinetic equation is obtained for the fluctuation of the macroscopic rank of the distribution function. The solution gives the transfer function and eddy viscosity. When applied to the inertia subrange of the energy spectrum the analysis recovers the Kolmogorov law and its numerical coefficient.

Description: The steady normal shock wave solutions of parahydrogen at various total pressures and total temperatures were numerically determined by iterating the upstream Mach number and by using a modified interval halving technique. The results obtained are compared with the ideal diatomic gas values and are presented in tabulated form.

Description: A mathematical model for the heat transfer within the electronics package of a Chaparral missile was performed. The Grashof number for this configuration was less than 2000 which indicated that the primary mode of heat transfer was conduction. The Vodicka theory for heat conduction in laminated composite media was utilized to obtain the solution for the model.

Description: A unique method was developed for the determination of heat transfer coefficients for water flowing through capillary tubes using a rastered electron beam heater. Heat flux levels of 150 and 500 watts/sq cm were provided on the top surface of four square tubes. Temperature gradient along the tube length and mass flow rates versus pressure drop were measured.

Description: Surface temperature gradients were measured with miniature thermocouples installed in a 58.5 cm (23-inch) square window. Test measurements at 25 locations were made under vacuum and with the window operating in radiant heat transfer mode. The analysis of thermocouple design and installation is presented along with a lead wire routing scheme to allow for both differential and absolute temperature measurements while using a minimum number of signal feedthru paths through the test chamber wall. Typical test data and operational precautions are presented along with the accuracy analysis for installation effects and measurement effects to support differential temperature measurement precision values of + or - 0.06 C RMS + or - 0.1 F RMS).

Description: Lock-hopper systems are the most common means for feeding solids to and from coal conversion reactor vessels. The rate at which crushed solids flow by gravity through the vertical pipes and valves in lock-hopper systems affects the size of pipes and valves needed to meet the solids-handling requirements of the coal conversion process. Methods used to predict flow rates are described and compared with experimental data. Preliminary indications are that solids-handling systems for coal conversion processes are over-designed by a factor of 2 or 3.

Description: Techniques for calculating high Reynolds number flow around an airfoil undergoing dynamic stall are reviewed. Emphasis is placed on predicting the values of lift, drag, and pitching moments. Methods discussed include: the discrete potential vortex method; thin boundary layer method; strong interaction between inviscid and viscous flows; and solutions to the Navier-Stokes equations. Empirical methods for estimating unsteady airloads on oscillating airfoils are also described. These methods correlate force and moment data from wind tunnel tests to indicate the effects of various parameters, such as airfoil shape, Mach number, amplitude and frequency of sinosoidal oscillations, mean angle, and type of motion.

Description: Unsteady flow phenomena are reviewed with emphasis on separated flow in the subsonic and transonic regimes. Specific topics discussed include external flows past bluff bodies, unsteady separation on slender bodies, and internal flows.

Description: Quasi-steady three dimensional separated flows about bodies of large fineness ratio operating at large angles of incidence or yaw are discussed. The general character of the three dimensional attached boundary layer, the concept of limiting streamlines, and the physics of three dimensional separation and reattachment are among the factors considered. Specific examples are given. The advantages of swept, sharp edges that generate controlled (or fixed) three dimensional flow separations on a vehicle, due to the qualitatively unchanging flow field developed throughout the range of flight conditions, are emphasized.

Description: The mathematical model and results of numerical solutions are given for the one dimensional problem when the linear equations are written in a rectangular coordinate system. All the computations are easily realizable for two and three dimensional problems when the equations are written in any coordinate system. Explicit and implicit schemes are shown in tabular form for stability and oscillations criteria; the initial temperature distribution is considered uniform.

Description: The calculation of flow fields past aircraft configuration at flight Reynolds numbers is considered. Progress in devising accurate and efficient numerical methods, in understanding and modeling the physics of turbulence, and in developing reliable and powerful computer hardware is discussed. Emphasis is placed on efficient solutions to the Navier-Stokes equations.

Description: The unsteady Kutta-Joukowski condition, dynamic stall on oscillating airfoils, and unsteady shock wave-boundary layer interaction are discussed. Emphasis is placed on developing reliable prediction techniques and suppression of unsteady separation on oscillating control surfaces, wings, and rotating blades to improve aerodynamic stability.

Description: Various techniques to control and reduce radiated noise and the application of these techniques to a 1/2-water Mach 5 quiet tunnel are reviewed. Measurements in a small scale nozzle have shown that the upstream part of the supersonic wall boundary layer could be maintained laminar up to Reynolds numbers of nearly 4 x 1 million based on the test region length upstream of the nozzle exit. Turbulent noise levels in this test region were then reduced by an order of magnitude. To maintain low noise levels at higher Reynolds numbers, laminar flow noise shields are required. Data are presented for shields that consist of small diameter rods alined nearly parallel to the entrance flow with small gaps between the rods for boundary layer suction. Analysis and data presented on the noise shielding and reflection characteristics of flat plates and a rod-wall test panel indicate that freestream turbulent noise can be reduced by 70 to 90 deg at high Reynolds numbers. Performance estimates for the 1/2-meter tunnel are based on these results.

Description: Navier-Stokes equation as discretized by new flux conserving method proposed by Chang and Scott results in the system: vector F(vector x) = 0, where F is a vector valued function. The Optimization method we use is based on Quasi-Newton methods: given a nonlinear function vector F(vector x) = 0, we solve, Delta(vector x) = -BF(vector x), where Delta(vector x) is the correction term and B is the inverse Jacobian of F(x). Then, iteratively, vector(x(sub (i+1))) = vector(x (sub i)) + alpha.Delta(vector x(sub i)), where alpha is a line search correction term determined by a line search routine. We use the BFCG's update the Jacobian matrix B(sub k) at each iteration. It is well known that B(sub k) approaches B(*) at the solution X(*). This algorithm has several advantages over the Newton-Raphson method. For example, we do not need to calculate the Jacobian matrix at each iteration which is computationally very expensive.

Description: Advanced airbreathing propulsion systems used in Mach 4-6 mission scenarios, usually consist of a single integrated turboramjet or as in this study, a turbojet housed in an upper bay with a separate ramjet housed in a lower bay. As the engines transition from turbojet to ramjet, there is an operational envelope where both engines operate simultaneously. One nozzle concept under consideration has a common nozzle, where the plumes from the turbojet and ramjet interact with one another as they expand to ambient conditions. In this paper, the two plumes interact at the end of a common 2-D cowl, when they both reach an approximate Mach 3.0 condition and then jointly expand to Mach 3.6 at the common nozzle exit plane. At this condition, the turbojet engine operated at a higher NPR than the ramjet, where the turbojet overpowers the ramjet plume, deflecting it approximately 12 degrees downward and in turn the turbojet plume is deflected 6 degrees upward. In the process, shocks were formed at the deflections and a shear layer formed at the confluence of the two jets. This particular case was experimentally tested and the data used to compare with the PARC3D code with k-kl two equation turbulence model. The 2-D and 3-D centerline CFD solutions are in good agreement, but as the CFD solutions approach the outer sidewall, a slight variance occurs. The outer wall boundary layers are thin and do not present much of an interaction, however, where the confluence interaction shocks interact with the thin boundary layer on the outer wall, strong vortices run down each shock causing substantial disturbances in the boundary layer. These disturbances amplify somewhat as they propagate downstream axially from the confluence point. The nozzle coefficient (CFG) is reduced 1/2 percent as a result of this sidewall interaction, from 0.9850 to 0.9807. This three-dimensional reduction is in better agreement with the experimental value of 0.9790.

Description: Fluid behavior in a low-g environment is controlled primarily by surface tension forces. Certain fluid and system characteristics determine the magnitude of these forces for both a free liquid surface and liquid in contact with a solid. These characteristics, including surface tension, wettability or contact angle, system geometry, and the relationships governing their interaction, are discussed. Various aspects of fluid behavior in a low-g environment are then presented. This includes the formation of static interface shapes, oscillation and rotation of drops, coalescence, the formation of foams, tendency for cavitation, and diffusion in liquids which were observed during the Skylab fluid mechanics science demonstrations. Liquid reorientation and capillary pumping to establish equilibrium configurations for various system geometries, observed during various free-fall (drop-tower) low-g tests, are also presented. Several passive low-g fluid storage and transfer systems are discussed. These systems use surface tension forces to control the liquid/vapor interface and provide gas-free liquid transfer and liquid-free vapor venting.

Description: The effect of roughness on the heat transfer distributions and the transition criteria for the windward pitch plane of the shuttle orbiter at an angle of attack of 30 deg was studied using data obtained in hypersonic wind tunnels. The heat transfer distributions and the transition locations for the roughened models were compared with the corresponding data for smooth models. The data were correlated using theoretical solutions for a nonsimilar, laminar boundary layer subject to two different flow field models for the orbiter.

Description: The stability of an electrically conducting fluid subjected to two dimensional disturbance was investigated. A physical system consisting of two parallel infinite vertical plates which are thermally insulated was studied. An external magnetic field of constant strength was applied to normal plates. The fluid was heated from below so that a steady temperature gradient was maintained in the fluid. The governing equations were derived by perturbation technique, and solutions were obtained by a modified Galerkin method. It was found that the presence of the magnetic field increases the stability of the physical system and instability can occur in the form of neutral or oscillatory instability.

Description: The stability problem of thermoviscoelastic fluid flow between rotating coaxial cylinders is investigated using nonlinear thermoviscoelastic constitutive equations due to Eringen and Koh. The velocity field is found to be identical with that of the classical viscous case and the case of the viscoelastic fluid, but the temperature and pressure fields are found to be different. By imposing some physically reasonable mechanical and geometrical restrictions on the flow, and by a suitable mathematical analysis, the problem is reduced to a characteristic value problem. The resulting problem is solved and stability criteria are obtained in terms of critical Taylor numbers. In general, it is found that thermoviscoelastic fluids are more stable than classical viscous fluids and viscoinelastic fluids under similar conditions.