ALBERT

Description: A modified strip analysis has been developed for rapidly predicting flutter of finite-span, swept or unswept wings at subsonic to hypersonic speeds. The method employs distributions of aerodynamic parameters which may be evaluated from any suitable linear or nonlinear steady-flow theory or from measured steady-flow load distributions for the underformed wing. The method has been shown to give good flutter results for a broad range of wings at Mach number from 0 to as high as 15.3. The principles of the modified strip analysis may be summarized as follows: Variable section lift-curve slope and aerodynamic center are substituted respectively, for the two-dimensional incompressible-flow values of 2 pi and quarter chord which were employed by Barmby, Cunningham, and Garrick. Spanwise distributions of these steady-flow section aerodynamic parameters, which are pertinent to the desired planform and Mach number, are used. Appropriate values of Mach number-dependent circulation functions are obtained from two-dimensional unsteady compressible-flow theory. Use of the modified strip analysis avoids the necessity of reevaluating a number of loading parameters for each value of reduced frequency, since only the modified circulation functions, and of course the reduced frequency itself, vary with frequency. It is therefore practical to include in the digital computing program a very brief logical subroutine, which automatically selects reduced-frequency values that converge on a flutter solution. The problem of guessing suitable reduced-frequency values is thus eliminated, so that a large number of flutter points can be completely determined in a single brief run on the computing machine. If necessary, it is also practical to perform the calculations manually. Flutter characteristics have been calculated by the modified strip analysis and compared with results of other calculations and with experiments for Mach numbers up to 15.3 and for wings with sweep angles from 0 degrees to 52.5 degrees, aspect ratios from 2.0 to 7.4, taper ratios from 0.2 to 1.0, and center-of-gravity positions between 34% chord and 59% chord. These ranges probably cover the great majority of wings that are of practical interest with the exception of very low-aspect-ratio surfaces such as delta wings and missile fins. This program has been implemented on the IBM 7094.

Description: The GTRAN program was developed to solve transient, as well as steady state, problems for gas piping systems. GTRAN capabilities allow for the analysis of a variety of system configurations and components. These include: multiple pipe junctions; valves that change position with time; fixed restrictions (orifices, manual valves, filters, etc.); relief valves; constant pressure sources; and heat transfer for insulated piping and piping subjected to free or forced convection. In addition, boundary conditions can be incorporated to simulate specific components. The governing equations of GTRAN are the one-dimensional transient gas dynamic equations. The three equations for pressure, velocity, and density are reduced to numerical equations using an implicit Crank-Nicholson finite difference technique. Input to GTRAN includes a description of the piping network, the initial conditions, and any events (e.g. valve closings) occuring during the period of analysis. Output includes pressure, velocity, and density versus time. GTRAN is written in FORTRAN 77 for batch execution and has been implemented on a DEC VAX series computer. GTRAN was developed in 1983.

Description: The SUBAERF2 program was developed to provide for the aerodynamic analysis and design of low speed wing flap systems. SUBAERF2 is based on a linearized theory lifting surface solution. It is particularly well suited to configurations which, because of high speed flight requirements, must employ thin wings with highly swept leading edges. The program is applicable to wings with either sharp or rounded leading edges. This program is a new and improved version of LAR-13116 and LAR-12987, which it replaces. The low speed aerodynamic analysis method used in SUBAERF2 provides estimates of wing performance which include the effects of attainable leading-edge thrust and vortex lift. This basic aerodynamic analysis method has been improved to provide for the convenient, efficient and accurate treatment of simple leading-edge and trailing-edge flap systems. The user inputs flap geometry directly. Solutions can be found for various combinations of leading and trailing edge flap deflections. The program provides for the simultaneous analysis of up to 25 pairs of leading-edge and trailing-edge flap deflection schedules. A revised attainable thrust algorithm improves accuracy at the low Mach numbers sometimes encountered in wind tunnel testing. Also added is a means of estimating the distribution of leading edge separation vortex forces. The revised program has been particularly useful in the subsonic analysis of vehicles designed for supersonic cruise. The SUBAERF2 program is written in FORTRAN V for batch execution and has been implemented on a CDC 175 computer operating under NOS 2.4 with a central memory requirement of approximately 115K (octal) of 60 bit words. This program was originally developed in 1983 and later revised in 1988.

Description: WAVDRAG calculates the supersonic zero-lift wave drag of complex aircraft configurations. The numerical model of an aircraft is used throughout the design process from concept to manufacturing. WAVDRAG incorporates extended geometric input capabilities to permit use of a more accurate mathematical model. With WAVDRAG, the engineer can define aircraft components as fusiform or nonfusiform in terms of non-intersecting contours in any direction or more traditional parallel contours. In addition, laterally asymmetric configurations can be simulated. The calculations in WAVDRAG are based on Whitcomb's area-rule computation of equivalent-bodies, with modifications for supersonic speed. Instead of using a single equivalent-body, WAVDRAG calculates a series of equivalent-bodies, one for each roll angle. The total aircraft configuration wave drag is the integrated average of the equivalent-body wave drags through the full roll range of 360 degrees. WAVDRAG currently accepts up to 30 user-defined components containing a maximum of 50 contours as geometric input. Each contour contains a maximum of 50 points. The Mach number, angle-of-attack, and coordinates of angle-of-attack rotation are also input. The program warns of any fusiform-body line segments having a slope larger than the Mach angle. WAVDRAG calculates total drag and the wave-drag coefficient of the specified aircraft configuration. WAVDRAG is written in FORTRAN 77 for batch execution and has been implemented on a CDC CYBER 170 series computer with a central memory requirement of approximately 63K (octal) of 60 bit words. This program was developed in 1983.

Description: An approximate inverse solution is presented for the nonequilibrium flow in the inviscid shock layer about a vehicle in hypersonic flight. The method is based upon a thin-shock-layer approximation and has the advantage of being applicable to both subsonic and supersonic regions of the shock layer. The relative simplicity of the method makes it ideally suited for programming on a digital computer with a significant reduction in storage capacity and computing time required by other more exact methods. Comparison of nonequilibrium solutions for an air mixture obtained by the present method is made with solutions obtained by two other methods. Additional cases are presented for entry of spherical nose cones into representative Venusian and Matrian atmospheres. A digital computer program written in FORTRAN language is presented that permits an arbitrary gas mixture to be employed in the solution. The effects of vibration, dissociation, recombination, electronic excitation, and ionization are included in the program.

Description: Small areas of high heat transfer and pressure can occur on a vehicle surface due to the influence of an impinging shock on the local flow. A method was needed to determine peak pressure and heating of these areas. This package is a system of computer programs designed to calculate two-dimensional shock interference patterns for six types of interference flows. Results also include properties of the inviscid flow field and the inviscid-viscous interaction at the surface along with peak pressure and peak heating at the impingement point. The six types of interference flow patterns considered are: 1) Type I interference patterns, occurring when two weak shocks of opposite families, BS (bow shock) and IS (impingment shock), intersect when the flow upstream of the impingement point is supersonic, or in the case of a blunt body, takes place well below the sonic point. 2) Type II interference pattern occurs when two shocks of opposite families (bow shock and impinging shock) intersect. Both shocks are weak as in type I, but are of such strength that in order to turn the flow, a Mach reflection must exist in the center of the flow field with an embedded subsonic region occurring between the intersection points (A & B) and the accompanying shear layers. Type II interference occurs on a blunt body when the impinging shock intersects the bow shock near the sonic point. 3) Type III shock interference pattern occurs when a weak impinging shock intersects a strong detached bow shock. On a blunt body the shock intersection occurs near or above the lower sonic point. 4) Type IV interference can occur when the impinging shock intersects a strong bow shock ahead of a subsonic flow region. On a blunt body this shock intersection is located between the lower sonic point and just above the body axis. The impinging shock causes a displacement of the bow shock and the formation of a supersonic jet that is embedded in the subsonic region. A jet bow shock is produced when the jet impinges on the surface, creating a small region with high stagnation heating. 5) Type V interference involves the interaction of two weak shocks of the same family. The interaction produces a shear layer, a supersonic jet, and a transmitted impinging shock. On a blunt body the shock interaction occurs near the upper sonic point. 6) Type VI interference involves the intersection of two weak shocks of the same family, which leads to an entirely supersonic flow field. This type of interference is important because it provides a means for predicting the onset of type V. Peak-heating correlations for laminar and turbulent shock-boundary-layer interactions are included in the programs for types I, II, V, and VI interference patterns. Heating correlations for laminar and turbulent reattaching shear layers obtained from separation studies are included in the program for type III interference. This program is written in FORTRAN IV for batch execution and has been implemented on a CDC 6000 Series computer. This program was developed in 1973.

Description: A computer program has been written to obtain the wave and friction drag of configurations with bodies of revolution and fins. These inviscid flow fields are superimposed and the wave drag of the configuration is obtained by integration of the surface pressures. The friction drag is obtained from the viscous flow field of the body and a flat-plate friction analysis of the fins. The numerical solution of these flow fields, superposition, and integration to obtain total drag have been programmed for high-speed digital computation. A large portion of the input required by the program is involved with the description of the configuration geometry and the specific surface positions where pressures are to be evaluated. In addition to drag forces, an output is available whereby the pressure distributions on the body and fins can be obtained.

Description: This paper is on the control of nonlinear-nonstationary vibration of a frame-stringer structure resulting from high levels of excitaation from a nearby supersonic jet exhaust. The structure exhibits periodic, chaotic, or random behaviors when forced by high-intensity sound from a supersonic jet exhaust with shock loading superimposed on the broadband response. The time history of the pressure, showing the rotation and flapping of the shock structure in the jet column due to large-scale instabilities, indicates that the response is not only nonlinear but also nonstationary. The acoustic pressure radiated by the structure also contains shocks and the formation of harmonics with distance. Control of the structural response is achieved by actively forcing the structure with an actuator at the shock oscillation frequency whose amplitude is locked into a self-control cycle. Results show that the peak power level is reduced by a factor of 63, or 18 dB. As a result, new broadband components emerge with at least four harmonics. At accelerating and decelerating supersonic speeds, the exhaust from the jet induces higher transient loading on the nearby flexible structure due to the occurence of multiple shocks from the jet.

Description: Results are presented for the numerical simulation of unsteady viscous incompressible flow past thick airfoils. Specifically, flow past a NACA 4424 at an angle of attack of 2.5 deg and Reynolds numbers in the range of 1700-4000 has been simulated using the spectral element method. At these conditions the flow is separatedd and an unsteady wake is formed. Application of the method of empirical eigenfunction reveals the structure of the most energetic components of the flow. These are found to occur in pairs that, through phase exchange, are responsible for the vortex shedding. A set of ordinary differential equations is obtained for the amplitudes of these eigenfunctions by a Galerkin projection of the Navier-Stokes equations. The solutions of the model system are compared with the full simulation. The work is of relevance to the transition process and observed routes to chaos in airfoil wakes.

Description: The objective of this study is to develop a reduced mechanism for ethylene oxidation. The authors are interested in a model with a minimum number of species and reactions that still models the chemistry with reasonable accuracy for the expected combustor conditions. The model will be validated by comparing the results to those calculated with a detailed kinetic model that has been validated against the experimental data.

Description: Artificial viscosity is added either implicity or explicitly in practically every numerical scheme for suppressing spurious oscillations in the solution of fluid-dynamics equations. In the present central-difference scheme, artificial viscosity is added explicitly for suppressing high-frequency oscillations and achieving good convergence properties. The amount of artificial viscosity added is controlled through the use of preselected coefficients. In the standard scheme, scalar coefficients based on the spectral radii of the Jacobian of the convective fluxes are used. However, this can add too much viscosity to the slower waves. Hence, the use of matrix-valued coefficients, which give appropriate viscosity for each wave component, is suggested. With the matrix-valued coefficients, the central-difference scheme produces more accurate solutions on a given grid, particularly in the vicinity of shocks and boundary layers, while still maintaining good convergence properties.

Description: A potential flow based three-dimensional panel method was modified to treat time-dependent conditions in which several submerged bodies can move within the fluid along different trajectories. This modification was accomplished by formulating the momentary solution in an inertial frame of reference, attached to the undisturbed stationary fluid. Consequently, the numerical interpretation of the multiple-body, solid-surface boundary condition and the viscous wake rollup was considerably simplified. The usteady capability of this code was calibrated and validated by comparing computed results with closed-form analytical results available for an airfoil, which was impulsively set into a constant speed forward motion. To demonstrate the multicomponent capability, computations were made for two wings following closely intersecting paths (i.e., simulations aimed at avoiding mid-air collisions) and for a flowfield with relative rotation (i.e., the case of a helicopter rotor rotating relative to the fuselage). Computed results for the cases were compared to experimental data, when such data was available.

Description: We report the discovery of a series of infrared absorption bands between 3600 and 3100/cm (2.8-3.2 micrometers) in the spectrum of Io. Individual narrow bands are detected at 3553, 3514.5, 3438, 3423, 3411.5, and 3401/cm (2.815, 2.845, 2.909, 2.921, 2.931, and 2.940 micrometers, respectively). The positions and relative strengths of these bands, and the difference of their absolute strengths between the leading and trailing faces of Io, indicate that they are due to SO2. The band at 3438/cm (2.909 micrometers) could potentially have a contribution from an additional molecular species. The existence of these bands in the spectrum of Io indicates that a substantial fraction of the SO2 on Io must reside in transparent ices having relatively large crystal sizes. The decrease in the continuum observed at the high frequency ends of the spectra is probably due to the low frequency side of the recently detected, strong 3590/cm (2.79 micrometer) feature. This band is likely due to the combination of a moderately strong SO2 band and an additional absorption from another molecular species, perhaps H2O isolated in SO2 at low concentrations. A broad (FWHM approximately = 40-60/cm), weak band is seen near 3160/cm (3.16 micrometers) and is consistent with the presence of small quantities of H2O isolated in SO2-rich ices. There is no evidence in the spectra for the presence of H2O vapor on Io. Thus, the spectra presented here neither provide unequivocal evidence for the presence of H2O on Io nor preclude it at the low concentrations suggested by past studies.

Description: An approximate method for calculating heating rates on three-dimensional vehicles at angle of attack is presented. The method is based on the axisymmetric analog for three-dimensional boundary layers and uses a generalized body-fitted coordinate system. Edge conditions for the boundary-layer solution are obtained from an inviscid flowfield solution, and because of the coordinate system used, the method is applicable to any blunt body geometry for which an inviscid flowfield solution can be obtained. The method is validated by comparing with experimental heating data and with thin-layer Navier-Stokes calculations on the shuttle orbiter at both wind-tunnel and flight conditions and with thin-layer Navier-Stokes calculations on the HL-20 at wind-tunnel conditions.

Description: The magnetic field configuration-states of the magnetotails of the planets Uranus and Neptune are compared. Earth's case is also briefly treated, as well as some related aspects of the other three magnetic planets. In Uranus' case, due to the large tilt (59 deg) of the planet's magnetic dipole with respect to its spin axis and the unusual obliquity of that axis, the angle of attack (alpha) of the solar wind with respect to dipole alignment goes through all possible angles, 0 deg to 180 deg, yielding a very broad spectrum of configuration-states of its tail. Cases are discussed where the planetary magnetic dipole is either aligned with the Sun-planet-line ('pole-on' state) or perpendicular to it and some intermediate states, for both Uranus and Neptune. Only Uranus experiences the pole-on state, which next occurs in November 1999 (+/- 2 months); last year (1993.2) it had the first 'perpendicular' state since Voyager encounter which resembles Earth's case. Neptune never has a pole-on configuration, but it gets as close as alpha = 14 deg from it; the next occurrence is early in 2003. At Voyager encounter Neptune's magnetotail apparently rapidly migrated through a broad spectrum of field structures with near extreme states resembling an Earth-like case on the one hand and a cylindrically symmetric one on the other. Magnetopause 'openness' should dramatically change in terms of the rapidly changing angle of attack throughout a planetary day for these two planets, and this has important implications for their magnetotails. Any future manetospheric mission plans for Uranus or Neptune should take in to consideration the allowed range of values for alpha for the epoch of interest; this is especially of concern for Uranus which has a pole-on state, and all possible alphas, around the middle of 2014, 20 years from now.

Description: We studied five new Antartic achondrites, MacAlpine Hills (MAC) 88177, Yamato (Y)74357, Y75274, Y791491 and Elephant Moraine (EET)84302 by mineralogical techniques to gain a better understanding of the mineral assemblages of a group of meteorites with an affinity to Lodran (stony-iron meteorite) and their formation processes. This group is being called lodranites. These meteorites contain major coarse-grained orthopyroxene (Opx) and olivine as in Lodran and variable amounts of FeNi metal and troilite etc. MAC88177 has more augite and less FeNi than Lodran; Y74357 has more olivine and contains minor augite; Y791491 contains in addition plagioclase. EET84302 has an Acapulco-like chondritic mineral assembladge and is enriched in FeNi metal and plagioclase, but one part is enriched in Opx and chromite. The EET84302 and MAC88177 Opx crystals have dusty cores as in Acapulco. EET84302 and Y75274 are more Mg-rich than other members of the lodranite group, and Y74357 is intermediate. Since these meteorites all have coarse-grained textures, similar major mineral assemblages, variable amounts of augite, plagioclase, FeNi metal, chromite and olivine, we suggest that they are related and are linked to a parent body with modified chondritic compositions. The variability of the abundances of these minerals are in line with a proposed model of the surface mineral assemblages of the S asteroids. The mineral assemblages can best be explained by differing degrees of loss or movements of lower temperature partial melts and recrystallization, and reduction. A portion of EET84302 rich in metal and plagioclase may represent a type of component removed from the lodranite group meteorites. Y791058 and Caddo County, which were studied for comparison, are plagioclase-rich silicate inclusions in IAB iron meteorites and may have been derived by similar process but in a different body.

Description: We have reduced high-titanium lunar mare soil and iron-rich lunar volcanic glass with hydrogen at temperatures of 900-1100 C. Ilmenite is the most reactive phase in the soil, exhibiting rapid and complete reduction at all temperatures. Ferrous iron in the glass is extensively reduced concurrent with partial crystallization. In both samples pyroxene and olivine undergo partial reduction along with chemical and mineralogical modifications. High-temperature reduction provides insight into the optical and chemical effects of lunar soil maturation, and places constraints on models of that process. Mare soil and volcanic glass are attractive feedstocks for lunar oxygen production, with achievable yields of 2-5 wt%.

Description: Recently Brown et al. (1991) showed that Triton's internal heat source could amount to 5-20% of the absorbed insolation on Triton, thus significantly affecting volatile transport and atmospheric pressure. Subsequently, Kirk and Brown (1991a) used simple analytical models of the effect of internal heat on the distribution of volatiles on Triton's surface, confirming the speculation of Brown et al. that Triton's internal heat flow could strongly couple to the surface volatile distribution. To further explore this idea, we present numerical models of the permanent distribution of nitrogen ice on Triton that include the effects of sunlight, the two-dimensional distribution of internal heat flow, the coupling of internal heat flow to the surface distribution of nitrogen ice, and the finite viscosity of nitrogen ice. From these models we conclude that: (1) The strong vertical thermal gradient induced in Triton's polar caps by internal heat-flow facilitates viscous spreading to lower latitudes, thus opposing the poleward transport of volatiles by sunlight, and, for plausible viscosities and nitrogen inventories, producing permanent caps of considerable latitudinal extent; (2) It is probable that there is a strong coupling between the surface distribution of nitrogen ice on Triton and internal heat flow; (3) Asymmetries in the spatial distribution of Triton's heat flow, possibly driven by large-scale, volcanic activity or convection in Triton's interior, can result in permanent polar caps of unequal latitudinal extent, including the case of only one permanent polar cap; (4) Melting at the base of a permanent polar cap on Triton caused by internal heat flow can significantly enhance viscous spreading, and, as an alternative to the solid-state greenhouse mechanism proposed by Brown et al. (1990), could provide the necessary energy, fluids, and/or gases to drive Triton's geyser-like plumes; (5) The atmospheric collapse predicted to occur on Triton in the next 20 years (Spencer, 1990) may be plausibly avoided because of the large latitudinal extent expected for permanent polar caps on Triton.

Description: In early 1994, the US will once again place a spacecraft in lunar orbit. Popularly called Clementine, this spacecraft will spend about two months mapping the Moon and then will travel on to encounter the near-Earth asteroid, 1620 Geographos. A brief description of the historical development of the mission; the lunar survey; Clementine's payload, equipment, and capabilties; and the encounter with Geographos is presented.

Description: This paper presents the results of a study on the variation of the critical Marangoni number (Ma(sub c)) for the onset of Benard convection in a finite liquid layer bounded horizontally as well as from below. A direct-numerical-simulation procedure is devised to determine the Ma(sub c) for aspect ratios (Ar) ranging from 0.8 to 10. The results predict a strong increase of Ma(sub c) as Ar decreases to below 2. A dip of Ma(sub c) occurs between Ar = 1.45 and 1.3, which is accompanied by a pattern transition from a two-cell convection to a unicellular flow. For Ar above 4, the calculated Ma(sub c) shows little change and asymptotically approach a value of 116.15, with Biot number (Bi) equal to 1.

Description: An experimental investigation carried out to determine aerodynamic and acoustic characteristics of a low area ratio rectangular jet ejector is reported. A supersonic primary jet issuing from a rectangular convergent-divergent nozzle of aspect ratio 4, into a rectangular duct of area ratio 3, was used. Improved performance was found when the ejector screech tone is most intense and appears to match the most unstable Strouhal number of the free rectangular jet. When the primary jet was operating at over and ideally expanded conditions, significant noise reduction was obtained with the ejector as compared to a corresponding free jet. Application of particle image velocimetry to high speed ejector flows was demonstrated through the measurement of instantaneous two dimensional velocity fields.

Description: Particle imaging velocimetry, which is rapidly becoming an essential tool for the synoptic measurement of two dimensional velocity fields, is addressed. This rapid development is carried out by a growing number of fluid mechanics experimentalists who recognize the unique capabilities of the technique to measure velocity fields in both space and time, and is supported by the concomitant development and availability of microcomputer hardware and software. Various applications of the technique to map different flow regimes are described and illustrated. Hardware implementations of the technique which utilize both conventional photography for image acquisition as well as digital 'online' methods for integrated image acquisition and processing are discussed. Recommendations to further enhance the technique and make it possible to map the three velocity components in three dimensional flow regions are given.

Description: The transitional flight characteristics of a geometrically simplified Short Take-Off Vertical Landing (STOVL) aircraft configuration have been measured in the NASA Ames 7- by 10-Foot Wind Tunnel. The experiment is the first in a sequence of tests designed to provide detailed data for evaluating the capability of computational fluid dynamics methods to predict the important flow parameters for powered lift. The model consists of a 60 deg cropped delta wing platform, blended fuselage and two circular in-line jets that exit perpendicularly from the flat lower surface. The measured flows have a maximum freestream Mach number of 0.2. Model angle of attack is varied between -10 and +20 deg. The flow is ambient temperature in both jet exits and the nozzle pressure ratios are varied between 1 and 3. The data presented includes forces and moments, pressures measured at 281 surface pressure ports and the pressures of the jets. Measurements of the flow are also made in the tunnel test section upstream and downstream of the model and at the jet exits to guide boundary condition selection for the planned computations. Flow visualization and total pressure measurements in the jet plumes provide a description of the three-dimensional jet efflux flowfield.

Description: This test was originally conducted to determine the effects of several empennage and afterbody parameters on the aft-end aerodynamic characteristics of a twin-engine fighter-type configuration. Model variables were as follows: horizontal tail axial location and incidence, vertical tail axial location and configuration (twin-vs single-tail arrangements), tail booms, and nozzle power setting. Jet propulsion was simulated by exhausting high-pressure, cold-flow air from the nozzles. Following a successful test conducted on a single engine nacelle model to validate a CFD code, this model was chosen to be instrumented with pressure taps on the afterbody and nozzles and used as a follow-on test, providing a more complex geometry for the CFD code validation. A more limited test matrix was run to collect the pressure data, employing only the twin-tail configuration and varying only the horizontal and vertical tail locations. Mach number was varied from 0.6 to 1.2. Nozzle pressure ratio was varied from jet-off to 8. Angle-of-attack varied from 0 to 8 deg.

Description: This test was initiated to provide validation data on low aspect ratio wings at transonic speeds. The test was conducted so that the data obtained would be useful in the validation of codes, and all boundary condition data required would be measured as part of the test. During the conduct of the test, the measured quantities were checked for repeatability, and when the data would not repeat, the cause was tracked down and either eliminated or included in the measurement uncertainty. The accuracy of the data was in the end limited by wall imperfections of the wind tunnel in which the test was run.

Description: The data presented in this contribution were obtained in the NASA Langley 16-Foot Transonic Tunnel. Multiple test entries were completed and the results have been completely reported in five NASA reports. The objective of the initial investigation was to determine the effect of empennage (tail) interference on the drag characteristics of an axisymmetric model with a single engine fighter aft-end with convergent divergent nozzles. Two nozzle power settings, dry and maximum afterburning, were investigated. Several empennage arrangements and afterbody modifications were investigated during the initial investigation. Subsequent investigations were used to determine the effects of other model variables including tail incidence, tail span, and nozzle shape. For the final investigation, extensive surface pressure instrumentation was added to the model in order to develop and understanding of the flow interactions associated with afterbody/empennage integration and also to provide data for code validation. Extensive computational analysis has been conducted on the staggered empennage configuration at a Mach number of 0.6 utilizing a three-dimensional Navier Stokes code. Most of the investigations were conducted at Mach numbers from 0.60 to 1.20 and at ratios of jet total pressure to free stream static pressure (nozzle pressure ratio) from 0.1 (jet off) to 8.0. Some angle of attack variation was obtained at jet off conditions.

Description: The purpose of this investigation is to provide a comprehensive data base for the validation of numerical simulations. The initial results of the study (single angle of attack) were presented in ref. 1, where the effects of various parameters and the adequacies of selected turbulence models were discussed. The objective of the present paper is to provide a tabulation of the experimental data. The data were obtained in the two-dimensional, transonic flowfield surrounding a supercritical airfoil. A variety of flows were studied in which the boundary layer at the trailing edge of the model was either attached or separated. Unsteady flows were avoided by controlling the Mach number and angle of attack. Surface pressures were measured on both the model and wind tunnel walls, and the flowfield surrounding the model was documented using a laser Doppler velocimeter (LDV). Although wall interference could not be completely eliminated, its effect was minimized by employing the following techniques. Sidewall boundary layers were reduced by aspiration, and upper and lower walls were contoured to accommodate the flow around the model and the boundary-layer growth on the tunnel walls. A data base with minimal interference from a tunnel with solid walls provides an ideal basis for evaluating the development of codes for the transonic speed range because the codes can include the wall boundary conditions more precisely than interference corrections can be made to the data sets.

Description: The computation of unsteady shock waves, which contribute significantly to noise generation in supersonic jet flows, is investigated. This paper focuses on the difficulties of computing slowly moving shock waves. Numerical error is found to manifest itself principally as a spurious entropy wave. Calculations presented are performed using a third order essentially nonoscillatory scheme. The effect of stencil biasing parameters and of two versions of numerical flux formulas on the magnitude of spurious entropy are investigated. The level of numerical error introduced in the calculation in quantified as a function of shock pressure ratio, shock speed, Courant number, and mesh density. The spurious entropy relative to the entropy jump across a static shock decreases with increasing shock strength and shock velocity relative to the grid, but is insensitive to Courant number. The structure of the spurious entropy wave is affected by the choice of flux formulas and algorithm biasing parameters. The effect of the spurious numerical waves on the calculation of sound amplification by a shock wave is investigated. For this class of problem, the acoustic pressure waves are relatively unaffected by the spurious numerical phenomena.

Description: One of the key technical elements in NASA's high speed research program is reducing the noise level to meet the federal noise regulation. The dominant noise source is associated with the supersonic jet discharged from the engine exhaust system. Whereas the turbulence mixing is largely responsible for the generation of the jet noise, a broadband shock-associated noise is also generated when the nozzle operates at conditions other than its design. For both mixing and shock noise components, because the source of the noise is embedded in the jet plume, one can expect that jet noise can be predicted from the jet flowfield computation. Mani et al. developed a unified aerodynamic/acoustic prediction scheme by applying an extension of Reichardt's aerodynamic model to compute turbulent shear stresses which are utilized in estimating the strength of the noise source. Although this method produces a fast and practical estimate of the jet noise, a modification by Khavaran et al. has led to an improvement in aerodynamic solution. The most notable feature in this work is that Reichardt's model is replaced with the computational fluid dynamics (CFD) solution of Reynolds-averaged Navier-Stokes equations. The major advantage of this work is that the essential, noise-related flow quantities such as turbulence intensity and shock strength can be better predicted. The predictions were limited to a shock-free design condition and the effect of shock structure on the jet mixing noise was not addressed. The present work is aimed at investigating this issue. Under imperfectly expanded conditions the existence of the shock cell structure and its interaction with the convecting turbulence structure may not only generate a broadband shock-associated noise but also change the turbulence structure, and thus the strength of the mixing noise source. Failure in capturing shock structures properly could lead to incorrect aeroacoustic predictions.

Description: A dimensionless group, called a pressure loss modulus (N(sub PL)), is introduced that, in conjunction with an appropriately defined Reynolds number, is of considerable engineering utility in correlating steady-state Delta p vs flow calibration data and subsequently as a predictor, using the same or a different fluid, in uniformly distributed pressure loss devices. It is particularly useful under operation in the transition regime. Applications of this simple bivariate correlation to three diverse devices of particular interest for small liquid rocket engine fluid systems are discussed: large L/D capillary tube restrictors, packed granular catalyst beds, and stacked vortex-loss disk restrictors.

Description: A numerical study of axisymmetric overexpanded nozzle is presented. The flow structure of the startup and throttle-down processes are examined. During the impulsive startup process, observed flow features include the Mach disk, separation shock, Mach stem, vortex core, contact surface, slip stream, initial shock front, and shocklet. Also the movement of the Mach disk is not monotonical in the downstream direction. For a range of pressure ratios, hysteresis phenomenon occurs; different solutions were obtained depending on different processes. Three types of flow structures were observed. The location of separation point and the lower end turning point of hysteresis are closely predicted. A high peak of pressure is associated with the nozzle flow reattachment. The reversed vortical structure and affects engine performance.

Description: Incompresible, thin sheet flows have been of research interest for many years. Those studies were mainly concerned with the stability of the flow in a surrounding gas. Squire was the first to carry out a linear, invicid stability analysis of sheet flow in air and compare the results with experiment. Dombrowski and Fraser did an experimental study of the disintegration of sheet flows using several viscous liquids. They also detected the formulation of holes in their sheet flows. Hagerty and Shea carried out an inviscid stability analysis and calculated growth rates with experimental values. They compared their calculated growth rates with experimental values. Taylor studied extensively the stability of thin liquid sheets both theoretically and experimentally. He showed that thin sheets in a vacuum are stable. Brown experimentally investigated thin liquid sheet flows as a method of application of thin films. Clark and Dumbrowski carried out second-order stability analysis for invicid sheet flows. Lin introduced viscosity into the linear stability analysis of thin sheet flows in a vacuum. Mansour and Chigier conducted an experimental study of the breakup of a sheet flow surrounded by high-speed air. Lin et al. did a linear stability analysis that included viscosity and a surrounding gas. Rangel and Sirignano carried out both a linear and nonlinear invisid stability analysis that applies for any density ratio between the sheet liquid and the surrounding gas. Now there is renewed interest in sheet flows because of their possible application as low mass radiating surfaces. The objective of this study is to investigate the fluid dynamics of sheet flows that are of interest for a space radiator system. Analytical expressions that govern the sheet geometry are compared with experimental results. Since a space radiator will operate in a vacuum, the analysis does not include any drag force on the sheet flow.

Description: Jovis Tholus, Ulysses Patera, and Biblis Patera, three small volcanoes in the Tharsis area of Mars, provide important insight into the evolution of volcanism on Mars. All three are interpreted to be shield volcanoes, indicating that shield volcansim was present from the outset in Tharsis. Jovis Tholus is the least complex with simple repeated outpouring of lavas and caldera-forming events. Ulysses Patera is dominated by a giant caldera within which is a line of cinder cones or domes suggesting terminal stages of volcanism in which the magma had either significant volatiles or increased viscosity. Biblis Patera is characterized by nested calderas which have expanded by block faulting of the flank; it also exhibits lava flows erupted onto the flanks from events along concentric fractures. These shields are different from the younger Tharsis Montes shields only in terms of the volume of erupted material. The limited shield volume suggests that the magma source which fed the shields was rapidly depleted. The relatively large size ofthe calderas probably results from relatively large, shallow magma bodies rather than significant burial of the flanks by younger lavas. Eruption rates consistent with typical terrestrial basaltic eruptiuon rates suggest that these volcanoes were probably built over time spans of 10(exp 4) to 10(exp 5) years. Stratigraphic ages range from Early to Upper Hesperian; absolute ages range from 1.9 to 3.4 Ga.

Description: The transition of an incompressible three-dimensional boundary layer with strong cross-flow is considered theoretically and computationally in the context of vortex/wave interactions. Specifically the work centers on two lower-branch Tollmien-Schlichting waves which mutually interact nonlinearly to induce a longitudinal vortex flow. The vortex motion in turn gives rise to significant wave modulation via wall-shear forcing. The characteristic Reynolds number is large and, as a consequence, the waves' and the vortex motion are governed primarily by triple deck theory. The nonlinear interaction is captured by a viscous partial-differential system for the vortex coupled with a pair of amplitude equations for each wave pressure. Following analysis and computation over a wide range of parameters, three distinct responses are found to emerge in the nonlinear behavior of the flow solution downstream: an algebraic finite-distance singularity, far-downstream saturation or far-downstream wave decay leaving pure vortex flow. These depend on the input conditions, the wave angles and the size of the cross flow.

Description: Carbonate samples from the 8.9-Mt nuclear (near-surface explosion) crater, OAK, and a terrestrial impact crater, Meteor Crater, were analyzed for shock damage using electron paramagnetic resonance (EPR). Samples from below the OAK apparent crater floor were obtained from six boreholes, as well as ejecta recovered from the crater floor. The degree of shock damage in the carbonate material was assessed by comparing the sample spectra to the spectra of Solenhofen and Kaibab limestone, which had been skocked to known pressures. Analysis of the OAK Crater borehole samples has identified a thin zone of allocthonous highly shocked (10-13 GPa) carbonate material underneath the apparent crater floor. This approx. 5- to 15-m-thick zone occurs at a maximum depth of approx. 125 m below current seafloor at the borehole, sited at the initial position of the OAK explosive, and decreases in depth towards the apparent crater edge. Because this zone of allocthonous shocked rock delineates deformed rock below, and a breccia of mobilized sand and collapse debris above, it appears to outline the transient crater. The transient crater volume inferred in this way is found to by 3.2 +/- 0.2 times 10(exp 6)cu m, which is in good agreement with a volume of 5.3 times 10(exp 6)cu m inferred from gravity scaling of laboratory experiments. A layer of highly shocked material is also found near the surface outside the crater. The latter material could represent a fallout ejecta layer. The ejecta boulders recovered from the present crater floor experienced a range of shock pressures from approx. 0 to 15 GPa with the more heavily shocked samples all occurring between radii of 360 and approx. 600 m. Moreover, the fossil content, lithology and Sr isotopic composition all demonstrate that the initial position of the bulk of the heavily shocked rock ejecta sampled was originally near surface rock at initial depths in the 32 to 45-m depth (below sea level) range. The EPR technique is also sensitive to prehistoric shock damage. This is demonstrated by our study of shocked Kaibab limestone from the 49,000-year-old Meteor (Barringer) Crater Arizona.

Description: A search through cycle 1, 2, and 3 Magellan radar data covering 98% of the surface of Venus revealed very few dunes. Only two possible dune fields and several areas that may contain microdunes smaller than the resolution of the images (75 m) were identified. The Aglaonice dune field was identified in the cycle 1 images by the specular returns characteristic of dune faces oriented perpendicular to the radar illumination. Cycle 1 and 2 data of the Fortuna-Meshkenet dune field indicate that there has been no noticeable movement of the dunes over an 8-month period. The dunes, which are oriented both parallel and perpendicular to the radar illumination, appear to be dark features on a brighter substrate. Bright and dark patches that were visible in either cycle 1 or 2 data, but not both, allowed identification of several regions in the southern part of Venus that may contain microdunes. The microdunes are associated with several parabolic crater deposits in the region and are probably similar to those formed in wind tunnel experiments under Venus-like conditions. Bragg scattering and/or subpixel relfections from the near-normal face on asymmetric microdunes may account for these bright and dark patches. Look-angle effects and the lack of sufficient sand-size particles seem to be most likely reasons so few dunes were identified in Magellan data. Insufficient wind speeds, thinness of sand cover, and difficulty in identifying isolated dunes may also be contributors to the scarcity of dunes.

Description: The Magellan and Pioneer Venus Orbiter radiometric tracking data sets have been combined to produce a 60th degree and order spherical harmonic gravity field. The Magellan data include the high-precision X-band gravity tracking from September 1992 to May 1993 and post-aerobraking data up to January 5, 1994. Gravity models are presented from the application of Kaula's power rule for Venus and an alternative a priori method using surface accelerations. Results are given as vertical gravity acceleration at the reference surface, geoid, vertical Bouguer, and vertical isostatic maps with errors for the vertical gravity and geoid maps included. Correlation of the gravity with topography for the different models is also discussed.

Description: The joule dissipation inside flux ropes in Venus' ionosphere is so great that they must be formed near, and maintained at, the place where they are observed. Thus ropes are not formed by a Kelvin-Helmholtz instability of the ionopause. The hypothesis that ropes may be formed by the dynamo action of internal gravity waves in Venus' thermosphere (Luhmann and Elphic, 1985; Cole, 1993) is strengthened by discussion of a magnetic evolution equation which includes neutral air motion. However, the dynamo process would work only at altitudes at which v(sub in) is greater than or equal to omega(sub i). At altitudes or parts of a rope where v(sub in) is much less than omega(sub i), the process does not work. A solar wind dynamo is therefore examined to account for the ropes. Thereby a major new heat source for ions of the Venus ionosphere associated with the ropes is uncovered.

Description: Analysis of CCD images of Triton obtained with the 1.5-m telescope on Palomar Mountain shows that in the time period surrounding the Voyager 2 encounter with the satellite (1985-1990), no changes in the satellite's visual albedo or color occurred. The published observations of Triton in the 0.35- to 0.60-micrometer spectral region obtained between 1950 and 1990 were reanalyzed to detect historical variability in both its albedo and visual color. Analysis of the photometry indicates that there is little, if any, change in Triton's visual geometric albedo. This result is consistent with the albedo pattern observed by Voyager and the change in sub-Earth latitude. Two distinct types of color changes are evident: a significant secular increase in the blue region of the visual spectrum since at least the 1950s, and the reported dramatic reddening of Triton's spectrum in the late 1970s. The latter change can be explained only by a short-lived geological phenomenon. Triton's changing pole orientation with respect to a terrestrial observer cannot explain the secular color changes. These changes imply volatile transport on a global scale on Triton's surface during the past 4 decades. We present two models which show that either removal of a red volatile from Triton's polar cap or deposition of a blue volatile in the equatorial regions can explain the secular color changes. A third possibility is that the changes are the result of the alpha-beta phase transition of nitrogen and subsequent fracturing of the polar cap region (N. S. Duxbury and R. H. Brown (1993).

Description: A computational method for predicting the two-phase transient fluid flow and heat transfer characteristics within a reservoir of the capillary-pumped-loop, intended to be used for spacecraft thermal management, has been developed. The model is based on the enthalpy formulation in an axisymmetric configuration. The reservoir operates under a constant thermodynamic pressure by allowing mass exchange between the reservoir and the outside loop. Both 1 g and 0 g environments have been considered to assess the effects of gravity on the reservoir performance. Depending on the gravity level, the power input and the reservoir orientation, three different convection modes have been identified, namely, the thermocapillary mode, the buoyancy mode, and the rapid-expansion mode (caused by interface movement). The impact of these modes on the performance of the reservoir and the associated physical phenomena have been discussed.

Description: The structure of the vortical flowfield over delta wings at high angles of attack was investigated. Three-dimensional Navier-Stokes numerical simulations were carried out to predict the complex leeward-side flowfield characteristics, including leading-edge separation, secondary separation, and vortex breakdown. Flows over a 75- and a 63-deg sweep delta wing with sharp leading edges were investigated and compared with available experimental data. The effect of variation of circumferential grid resolution grid resolution in the vicinity of the wing leading edge on the accuracy of the solutions was addressed. Furthermore, the effect of turbulence modeling on the solutions was investigated. The effects of variation of angle of attack on the computed vortical flow structure for the 75-deg sweep delta wing were examined. At moderate angles of attack no vortex breakdown was observed. When a critical angle of attack was reached, bubble-type vortex breakdown was found. With further increase in angle of attack, a change from bubble-type breakdown to spiral-type vortex breakdown was predicted by the numerical solution. The effects of variation of sweep angle and freestream Mach number were addressed with the solutions on a 63-deg sweep delta wing.

Description: Results are presented from an investigation of the structure of energized wakes in the presence of ground or ceiling planes when powered-lift devices are tested at zero and low wind tunnel velocities. Solutions are presented to indicate how the presence of a nearby ground plane causes the energized wake to constrict less than when in free space. In contrast, another set of solutions indicate that the presence of a ceiling plane enhances wake construction so that the wake width is even smaller than when in free space. Hence, a ground plane tends to reduce the chance for interaction between the wake and the energizing elements, while a nearby ceiling plane tends to increase the likelihood for interference.

Description: Despite the extensive experimental and computational data base in the literature on passive porosity, no clear explanation of the governing flow physics exists. It is theorized that the positive porosity concept modifies the external pressure loading by allowing communication between high- and low-pressure regions on the external surface. This study determines the dominant flow phenomena that govern the effectiveness of passive porosity. It aims to assess the contribution of each phenomenon as related to a porous slender axisymmetric forebody. To assess the influence of the mass transfer and pressure equalization phenomena on the effectiveness of passive porosity on slender axisymmetric forebodies, strakes were attached to the 5.0-caliber solid and porous forebodies to force crossflow separation. Longitudinal force and moment data were obtained at a Mach number of 0.1 over an angle-of-attack range of 0 to 55 deg.

Description: Neodymium, stontium, and chromium isotopic studies of the LEW86010 angrite established its absolute age and the formation interval between its crystallization and condensation of Allende CAIs from the solar nebula. Pyroxene and phosphate were found to contain approximately 8% of its Sm and Nd inventory. A conventional Sm-147-Nd-143 isochron yielded an age of 4.53 +/- 0.04 Ga (2 sigma and Epsilon(sub Nd sup 143)) = 0.45 +/- 1.1. An Sm-146-Nd-142 isochron gives initial Sm-146/Sm-144 = 0.0076 +/- 0.0009 and Epsilon (sub Nd sup 142) = -2.5 +/- 0.4. The Rb-Sr analyses give initial Sr-87/Sr-86 Iota(sub Sr sup 87) = 0.698972 +/- 8 and 0.698970 +/- 18 for LEW and ADOR, respectively, relative to Sr-87/Sr-86 = 0.71025 for NBS987. The difference, Delta Iota(sub Sr Sup 87), between Iota (sub sr sup 87) for the angrites and literature values for Allende CAIs, corresponds to approximately Ma of growth in a solar nebula with a CI chondrite value of Rb-87/Sr-86 = 0.91, or approximately 5 Ma in a nebula with solar photospheric Rb-87/Sr-86 = 1.51. Excess Cr-53 from extinct Mn-53(t(sub 1/2) = 3.7 Ma)in LEW86010 corresponds to initial Mn-53/Mn-55 = 4.4 +/- 1.0 x 10(exp -5) for the inclusions as previously reported by the Paris group (Birck and Allegre, 1988). The Sm-146/Sm-144 value found for LEW86010 corresponds to solar system initial (Sm-146/Sm-144) = 0.0080 +/- 0.0009 for crystallization 8 Ma after Allende, the difference between Pb-Pb ages of angrites and Allende, or 0.0086 +/- 0.0009 for crystallation 18 Ma after Allende, using the Mn-Cr formation interval. The isotopic data are discussed in the context of a model in which an undifferentiated 'chondritic' parent body formed from the solar nebula approximately Ma after Allende CAIs and subsequently underwent differentiation accompanied by loss of volatiles. Parent bodies with Rb/Sr similar to that of CI, CM, or CO chondrites could satisfy the Cr and Sr isotopic systematics. If the angrite parent body had Rb/Sr similar to that of CV meteorites, it would have to form slightly later, approximately 2.6 Ma after the CAIs, to satisfy the Sr and CR isotopic systematics.

Description: The diagnostic analysis of numerical simulations of the Venus/Titan wind regime reveals an overlooked constraint upon the latitudinal structure of their zonal-mean angular momentum. The numerical experiments, as well as the limited planetary observations, are approximately consistent with the hypothesis that within the latitudes bounded by the wind maxima the total Ertel potential vorticity associated with the zonal-mean motion is approximately well mixed with respect to the neutral equatorial value for a stable circulation. The implied latitudinal profile of angular momentum is of the form M equal to or less than M(sub e)(cos lambda)(exp 2/Ri), where lambda is the latitude and Ri the local Richardson number, generally intermediate between the two extremes of uniform angular momentum (Ri approaches infinity) and uniform angular velocity (Ri = 1). The full range of angular momentum profile variation appears to be realized within the observed meridional - vertical structure of the Venus atmosphere, at least crudely approaching the implied relationship between stratification and zonal velocity there. While not itself indicative of a particular eddy mechanism or specific to atmospheric superrotation, the zero potential vorticity (ZPV) constraint represents a limiting bound for the eddy - mean flow adjustment of a neutrally stable baroclinic circulation and may be usefully applied to the diagnostic analysis of future remote sounding and in situ measurements from planetary spacecraft.

Description: We analyzed ropy glasses from Apollo 12 soils 12032 and 12033 by a variety of techniques including SEM/EDX, electron microprobe analysis, INAA, and Ar-39-Ar-40 age dating. The ropy glasses have potassium rare earth elements phosphorous (KREEP)-like compositions different from those of local Apollo 12 mare soils; it is likely that the ropy glasses are of exotic origin. Mixing calculations indicate that the ropy glasses formed from a liquid enriched in KREEP and that the ropy glass liquid also contained a significant amount of mare material. The presence of solar Ar and a trace of regolith-derived glass within the ropy glasses are evidence that the ropy glasses contain a small regolith component. Anorthosite and crystalline breccia (KREEP) clasts occur in some ropy glasses. We also found within these glasses clasts of felsite (fine-grained granitic fragments) very similar in texture and composition to the larger Apollo 12 felsites, which have a Ar-39-Ar-40 degassing age of 800 +/- 15 Ma. Measurements of 39-Ar-40-Ar in 12032 ropy glass indicate that it was degassed at the same time as the large felsite although the ropy glass was not completely degassed. The ropy glasses and felsites, therefore, probably came from the same source. Most early investigators suggested that the Apollo 12 ropy glasses were part of the ejecta deposited at the Apollo 12 site from the Copernicus impact. Our new data reinforce this model. If these ropy glasses are from Copernicus, they provide new clues to the nature of the target material at the Copernicus site, a part of the Moon that has not been sampled directly.

Description: A preconditioned Krylov subspace method (GMRES) is used to solve the linear systems of equations formed at each time-integration step of the unsteady, two-dimensional, compressible Navier-Stokes equations of fluid flow. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux-split formulation. Several preconditioning techniques are investigated to enhance the efficiency and convergence rate of the implicit solver based on the GMRES algorithm. The superiority of the new solver is established by comparisons with a conventional implicit solver, namely line Gauss-Seidel relaxation (LGSR). Computational test results for low-speed (incompressible flow over a backward-facing step at Mach 0.1), transonic flow (trailing edge flow in a transonic turbine cascade), and hypersonic flow (shock-on-shock interactions on a cylindrical leading edge at Mach 6.0) are presented. For the Mach 0.1 case, overall speedup factors of up to 17 (in terms of time-steps) and 15 (in terms of CPU time on a CRAY-YMP/8) are found in favor of the preconditioned GMRES solver, when compared with the LGSR solver. The corresponding speedup factors for the transonic flow case are 17 and 23, respectively. The hypersonic flow case shows slightly lower speedup factors of 9 and 13, respectively. The study of preconditioners conducted in this research reveals that a new LUSGS-type preconditioner is much more efficient than a conventional incomplete LU-type preconditioner.

Description: A grid generation and flow solution algorithm for the Euler equations on unstructured grids is presented. The grid generation scheme utilizes Delaunay triangulation and self-generates the field points for the mesh based on cell aspect ratios and allows for clustering near solid surfaces. The flow solution method is an implicit algorithm in which the linear set of equations arising at each time step is solved using a Gauss Seidel procedure which is completely vectorizable. In addition, a study is conducted to examine the number of subiterations required for good convergence of the overall algorithm. Grid generation results are shown in two dimensions for a National Advisory Committee for Aeronautics (NACA) 0012 airfoil as well as a two-element configuration. Flow solution results are shown for two-dimensional flow over the NACA 0012 airfoil and for a two-element configuration in which the solution has been obtained through an adaptation procedure and compared to an exact solution. Preliminary three-dimensional results are also shown in which subsonic flow over a business jet is computed.

Description: On 1992 August 14 at 12:40 UTC, an ordinary chondrite of type L5/6 entered the atmosphere over Mbale, Uganda, broke up, and caused a strewn field of size 3 x 7 km. Shortly after the fall, an expedition gathered eye witness accounts and located the position of 48 impacts of masses between 0.19 and 27.4 kg. Short-lived radionuclide data were measured for two specimens, one of which was only 12 days after the fall. Subsequent recoveries of fragements has resulted in a total of 863 mass estimates by 1993 October. The surfaces of all fragments contain fusion crust. The meteorite shower caused some minor inconveniences. Most remarkably, a young boy was hit on the head by a small specimen. The data interpreted as to indicate that the meteorite had an initial mass between 400-1000 kg (most likely approximately 1000 kg) and approached Mbale from AZ = 185 +/- 15, H = 55 +/- 15, and V(sub infinity) = 13.5 +/- 1.5/s. Orbital elements are given. Fragmentation of the initial mass started probably above 25 km altitude, but the final catastrophic breakup occurred at an altitude of 10-14 km. An estimated 190 +/- 40 kg reached the Earth's surface minutes after the final breakup of which 150 kg of material has been recovered.

Description: ELM is a simple computational tool for modeling the steady-state thermal-hydraulics of propellant flow through fuel element coolant channels in nuclear thermal rockets. Written for the nuclear propulsion project of the Space Exploration Initiative, ELM evaluates the various heat transfer coefficient and friction factor correlations available for turbulent pipe flow with heat addition. In the past, these correlations were found in different reactor analysis codes, but now comparisons are possible within one program. The logic of ELM is based on the one-dimensional conservation of energy in combination with Newton's Law of Cooling to determine the bulk flow temperature and the wall temperature across a control volume. Since the control volume is an incremental length of tube, the corresponding pressure drop is determined by application of the Law of Conservation of Momentum. The size, speed, and accuracy of ELM make it a simple tool for use in fuel element parametric studies. ELM is a machine independent program written in FORTRAN 77. It has been successfully compiled on an IBM PC compatible running MS-DOS using Lahey FORTRAN 77, a DEC VAX series computer running VMS, and a Sun4 series computer running SunOS UNIX. ELM requires 565K of RAM under SunOS 4.1, 360K of RAM under VMS 5.4, and 406K of RAM under MS-DOS. Because this program is machine independent, no executable is provided on the distribution media. The standard distribution medium for ELM is one 5.25 inch 360K MS-DOS format diskette. ELM was developed in 1991. DEC, VAX, and VMS are trademarks of Digital Equipment Corporation. Sun4 and SunOS are trademarks of Sun Microsystems, Inc. IBM PC is a registered trademark of International Business Machines. MS-DOS is a registered trademark of Microsoft Corporation.

Description: Observations of molten mid-ocean ridge basalt (MORB)-molybdenum (Mo) interactions produced by shock experiments provide insight into impact and differentiation processes involving metal-silicate partitioning. Analysis of fragments recovered from experiments (achieving MORB liquid shock pressures from 0.8 to 6 GPa) revealed significant changes in the composition of the MORB and Mo due to reaction of the silicate and metal liquids on a short time scale (less than 13 s). The FeO concentration of the shocked liquid decreases systematically with increasing pressure. In fact, the most highly shocked liquid (6 GPa) contains only 0.1 wt% FeO compared to an initial concentration of 9 wt% in the MORB. We infer from the presence of micrometer-sized Fe-, Si- and Mo-rich metallic spheres in the shocked glass that the Fe and Si oxides in the MORB were reduced in an estimated oxygen fugacity of 10(exp -17) bar and subsequently alloyed with the Mo. The in-situ reduction of FeO in the shocked molten basalt implies that shock-induced reduction of impact melt should be considered a viable mechanism for the formation of metallic phases. Similar metallic phases may form during impact accretion of planets and in impacted material found on the lunar surface and near terrestrial impact craters. In particular, the minute, isolated Fe particles found in lunar soils may have formed by such a process. Furthermore, the metallic spheres within the shocked glass have a globular texture similar to the textures of metallic spheroids from lunar samples and the estimated, slow cooling rate of less than or equal to 140 C/s for our spheres is consistent with the interpretation that the lunar spheroids formed by slow cooling within a melted target.

Description: We propose a new type of telescope designed specifically for the lunar environment of high vacuum and low temperature. Large area UV-Visible-IR telescope arrays can be built with ultra-light-weight replica optics. High T(sub c) superconductors provide support, steering, and positioning. Advantages of this approach are light-weight payload compatible with existing launch vehicles, configurable large area optical arrays, no excavation or heavy construction, and frictionless electronically controlled mechanisms. We have built a prototype and will be demonstarting some of its working characteristics.

Description: The ultimate imaging resolution in the UV and photometric precision achievable with a small (less than 1-meter) telescope located on the Moon is considered. The imaging resolution and photometric precision that might be practically achieved when the effects of the Lunar environment and equipment limitations are accounted for is then suggested. Finally, the practicality of soft landing such a telescope on the moon is considered, along with suggestions of how it might be directly controlled by using astronomers without any significant permanent staff.

Description: The 21st century is likely to see the start of the manned exploration and settlement of the inner solar system. NASA's plans for this endeavor are focused upon the Space Exploration Initiative which calls for a return to the Moon, to stay, followed by manned missions to Mars. To execute these missions safely provides solar physics with both a challenge and an opportunity. As the past solar maximum has clearly demonstrated, the Sun, through the solar flare process, is capable of generating and accelerating to high energies large fluxes of protons whose cumulative dose to unprotected astronauts can be fatal. It will be the responsibility of solar physicists to develop an accurate physical description of the mechanisms of flare energy storage and release, and of particle acceleration and propagation through interplanetary space upon which to base a sound method of flare and energetic particle prediction.

Description: The Jet Propulsion Laboratory Contamination Analysis Program (CAP) is a generalized transient executive analysis computer code which solves realistic mass transport problems in the free molecular flow environment. These transport problems involve mass flux from surface source emission and re-emission, venting, and engine emission. CAP solution capability allows for one-bounce mass reflections if required. CAP was developed to solve thin-film contamination problems in the free molecular flow environment, the intent being to provide a powerful analytic tool for evaluating spacecraft contamination problems. The solution procedure uses an enclosure method based on a lumped-parameter multinodal approach with mass exchange between nodes. Transient solutions are computed by the finite difference Euler method. First-order rate theory is used to represent surface emission and reemission (user care must be taken to insure the problem is appropriate for such behavior), and all surface emission and reflections are assumed diffuse. CAP does not include the effects of post-deposition chemistry or interaction with the ambient atmosphere. CAP reads in a model represented by a multiple-block data stream. CAP allows the user to edit the input data stream and stack sequential editing operations (or cases) in order to make complex changes in behavior (surface temperatures, engine start-up and shut-down, etc.) in a single run if desired. The eight data blocks which make up the input data stream consist of problem control parameters, nodal data (area, temperature, mass, etc.), engine or vent distribution factors (based upon plume definitions), geometric configuration factors (diffuse surface emission), surface capture coefficient tables, source emission rate constant tables, reemission rate constant tables, and partial node to body collapse capability (for deposition rates only). The user must generate this data stream, since neither the problem-specific geometric relationships, the constituents involved, nor plume distribution functions are a part of CAP. Instead, these are used to generate the data stream model CAP solves. Outputs vary from individual deposition rates of exchange, on an internodal basis and on a constituent basis as a function of time, to deposition on each surface on a constituent basis as a function of time. The type of outputs may be user-specified by control parameters. CAP allows the user to select output intervals within the solution interval and to generate restart nodal data blocks. CAP is composed of several FORTRAN subroutines which serve specific functions and can be easily edited. The code is relatively small (2152 statements), and contains comment statements for all operations. It is written in relatively generic FORTRAN to be adaptable to a variety of computers. CAP was implemented on a DEC VAX 11/780 computer, and is distributed on a 9-track DEC VAX BACKUP format magnetic tape. Virtual memory required is 4.6 MB, which corresponds to a 900 node model capacity. CAP was originally developed under contract for NASA/Goddard Space Flight Center in 1979 by JPL, and was subsequently modified as required for project support at JPL. CAP is a copyrighted work with all copyright vested in NASA.

Description: A robust and efficient technique is discussed for surface-grid generation on a general curvilinear surface. This technique is based on a nonuniform parameter space and allows for the generation of surface grids on highly skewed and nonuniform spaced background surface-grids. This method has been successfully integrated into the GRIDGEN software system.

Description: The Space Exploration Initiative presents an opportunity to construct astronomical telescopes on the Moon using the infrastructure provided by the lunar outpost. Small automatically deployed telescopes can be carried on the survey missions, be deployed on the lunar surface and be operated remotely from the Earth. Possibilities for early, small optical telescopes are a zenith pointed transit telescope, a student telescope, and a 0.5 to 1 meter automatic, fully steerable telescope. After the lunar outpost is established the lunar interferometers will be constructed in an evolutionary fashion. There are three lunar interferometers which have been studied. The most ambitious is the optical interferometer with a 1 to 2 -km baseline and seven 1.5 aperture elements arranged in a 'Y' configuration with a central beam combiner. The Submillimeter interferometer would use seven, 5-m reflectors in a 'Y' or circular configuration with a 1-km baseline. The Very Low Frequency (VLF) array would operate below 30 mHz with as many as 100 elements and a 200-km baseline.

Description: An observational approach to Planetary Sciences and exploration from Earth applies to a quite limited number of targets, but most of these are spatially complex, and exhibit variability and evolution on a number of temporal scales which lie within the scope of possible observations. Advancing our understanding of the underlying physics requires the study of interactions between the various elements of such systems, and also requires study of the comparative response of both a given object to various conditions and of comparable objects to similar conditions. These studies are best conducted in 'campaigns', i.e. comprehensive programs combining simultaneous coherent observations of every interacting piece of the puzzle. The requirements include both imaging and spectroscopy over a wide spectral range, from UV to IR. While temporal simultaneity of operation in various modes is a key feature, these observations are also conducted over extended periods of time. The moon is a prime site offering long unbroken observation times and high positional stability, observations at small angular separation from the sun, comparative studies of planet Earth, and valuable technical advantages. A lunar observatory should become a central piece of any coherent set of planetary missions, supplying in-situ explorations with the synoptic and comparative data necessary for proper advance planning, correlative observations during the active exploratory phase, and follow-up studies of the target body or of related objects.

Description: JTMIX was written to allow the prediction of both ideal and realistic properties of mixed gases in the 65-80K temperature range. It allows mixed gas J-T analysis for any fluid combination of neon, nitrogen, various hydrocarbons, argon, oxygen, carbon monoxide, carbon dioxide, and hydrogen sulfide. When used in conjunction with the NIST computer program DDMIX, JTMIX has accurately predicted order-of-magnitude increases in J-T cooling capacities when various hydrocarbons are added to nitrogen, and it predicts nitrogen normal boiling point depressions to as low as 60K when neon is added. JTMIX searches for heat exchanger "pinch points" that can result from insolubility of various components in each other. These points result in numerical solutions that cannot exist. The length of the heat exchanger is searched for such points and, if they exist, the user is warned and the temperatures and heat exchanger effectiveness are corrected to provide a real solution. JTMIX gives very good correlation (within data accuracy) to mixed gas data published by the USSR and data taken by APD for the U.S. Naval Weapons Lab. Data taken at JPL also confirms JTMIX for all cases tested. JTMIX is written in Turbo C for IBM PC compatible computers running MS-DOS. The National Institute of Standards and Technology's (NIST, Gaithersburg, MD, 301-975-2208) computer code DDMIX is required to provide mixed-fluid enthalpy data which is input into JTMIX. The standard distribution medium for this program is a 5.25 inch 360K MS-DOS format diskette. JTMIX was developed in 1991 and is a copyrighted work with all copyright vested in NASA.

Description: Knowledge of the gravitation field, in combination with surface topography, provides one of the principal means of inferring the internal structure of a planetary body. Previous analyses of the lunar gravitational field have been based on data from the Lunar Orbiters, the Apollo subsatellites, and the low altitude passes of the Apollo spacecraft. Recently, Konopliv et al. have reanalyzed all available Lunar Orbiter and Apollo subsatellite tracking data, producing a 60th degree and order solution. In preparation for the Clementine Mission to the Moon, we have also initiated a reanalysis of the Lunar Orbiter and Apollo subsatellite data. Our reanalysis takes advantage of advanced force and measurement modeling techniques as well as modern computational facilities. We applied the least squares collocation technique which stabilizes the behavior of the solution and high degree and order. The extension of the size of the field reduces the aliasing coming from the omitted portion of the gravitational field. This is especially important for the analysis of the tracking data from the Lunar Orbiters, as the periapse heights frequently ranged from 50 to 100 km.

Description: The 108 photons of the Martian He 584 A airglow detected by the Extreme Ultraviolet Explorer satellite during a two-day exposure (22-23 Jan. 1993) correspond to the effective disk average intensity of 43 (+/-) 10 Rayleigh. Radiative transfer calculations, using a model atmosphere appropriate to the conditions of the observation and having an exospheric temperature of 210 (+/-) 20 K, result in an He mixing ratio of 1.1 (+/-) 0.4 ppm in the lower atmosphere. Nonthermal escape of helium is due to the following: electron impact ionization and pickup of He(+) by the solar wind; collisions with hot oxygen atoms; and charge exchange with molecular species with corresponding column loss rates of 1.4 x 10(exp 5), 3 x 10(exp 4), and 7 x 10(exp 3) cm(exp -2)s(exp -1), respectively. The lifetime of helium on Mars is 5 x 10(exp 4) yr. The He outgassing rate, coupled with the Ar-40 atmospheric abundance and with the K:U:Th ratio measured in the surface rocks, is used as input to a simple two-reservoir degassing model which presumes the loss of all argon accumulated in the atmosphere during the first Byr by large-scale impacts. The model results in total planet mass ratios of 10(exp -5) g/g for K, 2.3 x 10(exp -9) g/g for U, 8.5 x 10(exp -9) g/g for Th, 4 x 10(exp -10) g/g for He, and 1.5 x 10(exp -9) g/g for Ar-40. The predicted radiogenic heat flux is 2 erg cm(exp -2)s(exp -1). Similar modeling for Venus results in total plant mass ratios of 4.7 x 10(exp -5) g/g for K, 6.7 x 10(exp -9) g/g for U, 2.2 x 10(exp -8) g/g for Th, 1.3 x 10(exp -9) for He, 6.7 x 10(exp -9) g/g for Ar-40, and a radiogenic heat flux of 15 erg cm(exp -2)s(exp -1). The implications of these results are discussed. The modeling shows that the radioactive elements were not distributed uniformly in the protoplanetary nebula, and their relative abundances differ very much in the terrestrial planets.

Description: Very small variations in Nd-142 abundance in SNC meteorites lunar basalts, and a terrestrial supracrustal rock, have been attributed to the decay of 103 Ma Sm-146 initially present in basalt source regions in varying abundances as a result of planetary differentiation. We previously interpreted variations in Nd-142 abundances in two Apollo 17 high-Ti basalts, three Apollo 12 low-Ti basalts, and two KREEP basalts as defining an isochron giving a formation interval of approximately 94 Ma for the lunar mantle. Here we report new data for a third Apollo 17 high-Ti basalt, two Apollo 15 low-Ti basalts, the VLT basaltic lunar meteorite A881757 (formerly Asuka 31), basalt-like KREEP impact melt rocks 14310 and 14078, and three terrestrial rock standards. Those lunar samples which were not exposed to large lunar surface thermal neutron fluences yield a revised mantle formation interval of 237 +/- 64 Ma.

Description: A full-scale F/A-18 was tested in the 80 by 120-Foot Wind Tunnel at NASA Ames Research Center to measure the effectiveness of a tangentially blowing slot in generating significant yawing moments while minimizing coupling in the pitch and roll axes. Various slot configurations were tested to determine the optimum configuration. The test was conducted for angles of attack from 25 to 50 deg, angles of sideslip from -15 to +15 deg, and freestream velocities from 67 ft/sec to 168 ft/sec. By altering the forebody vortex flow, yaw control was maintained for angles of attack up to 50 deg. Of particular interest was the result that blowing very close to the radome apex was not as effective as blowing slightly farther aft on the radome, that a 16-inch slot was more efficient, and that yawing moments were generated without inducing significant rolling or pitching moments.

Description: Agglutinitic glass contains much of the reduced Fe in lunar soils, and it contributes to the modification of reflectance spectra from lunar soils. Previous work has shown that agglutinitic glass can be compositionally heterogeneous, but the scale of these heterogeneities is not well known. In addition, few data are available on the characteristics of the inclusions in agglutinitic glass. Here we report on our preliminary transmission electron microscope (TEM) examination of agglutinitic glass fragments from the Apollo 11 soil 10084.

Description: Periodic Comet Shoemaker/Levy 9 will impact Jupiter in late July 1994. The comet, which broke into more than 20 telescopically detectable fragments when it passed with the Roche lobe of Jupiter on July 8, 1992, is captured in a highly eccentric orbit about Jupiter. The 21 recognized nuclei will be spread out in a train of the order 7 x 10(exp 6) km long at the time of impact, and the impacts will be spread in time over about 5 1/2 days centered on about July 21.2 UT. In addition to the train of recognized bright nuclei, the comet consists of 'wings' of unresolved bodies that are the source of a very broad composite dust tail. The linear extent of the wings is about an order of magnitude greater than that of the train of recognized discrete nuclei. Collision of the wings will be spread in time over several months. Thus the impact of P/S-L 9 with Jupiter will be an event of appreciable duration.

Description: Approximately 950 impact craters have been identified on the surface of Venus, mainly in Magellan radar images. From a combination of Earth-based Arecibo, Venera 15/1, and Magellan radar images, we have interpreted 72 as unequivocal peak-ring craters and four as multiringed basins. The morphological and structural preservation of these craters is high owing to the low level of geologic activity on the venusian surface (which is in some ways similar to the terrestrial benthic environment). Thus these craters should prove crucial to understanding the mechanics of ringed crater formation. They are also the most direct analogs for craters formed on the Earth in Phanerozoic time, such as Chicxulub. We summarize our findings to date concerning these structures.

Description: A ray theory approach is used to examine the propagation of sonic booms through a turbulent ground layer, and to make predictions about the received waveform. The rays are not propagated one at a time, as is typical in ray theory; instead, sufficient rays to represent a continuous wave front are propagated together. New rays are interpolated as needed to maintain the continuity of the wave front. In order to predict the received boom signature, the wave front is searched for eigenrays after it has propagated to the receiver. The Comte-Bellot turbulence model is used to generate an instantaneous 'snapshot' of the turbulent field. The transient acoustic wave is assumed to be sufficiently short in duration such that the time-dependacde of the turbulent field may be neglected.

Description: We propose sending a balloon-borne UV photometer sensor package to measure atmospheric ozone on Mars, and this package could be a Discovery Program sensor candidate. Past measurements of ozone on Mars are highly uncertain, perhaps a factor of 3 or so uncertain, due primarily to interference and masking by cloud and dust. In-situ balloon measurements would avoid these problems, and would provide 'ground truth' which remote sensing techniques cannot. We have explored this approach to measure ozone abundance in the terrestrial stratosphere with a balloon-borne UV absorption photometer. Atmospheric pressures and temperatures and ozone concentrations near the surface of Mars are similar to those in the terrestrial stratosphere.

Description: The Hydroblast Research Cell at Marshall Space Flight Center is used to investigate the use of high pressure waterjets to strip paint, grease, adhesive and thermal spray coatings from various substrates. Current methods of cleaning often use ozone depleting chemicals (ODC) such as chlorinated solvents. High pressure waterjet cleaning has proven to be a viable alternative to the use of solvents. A popular method of waterjet cleaning involves the use of a rotating, multijet, high pressure water nozzle which is robotically controlled. This method enables rapid cleaning of a large area, but problems such as incomplete coverage and damage to the substrate from the waterjet have been observed. This report summarizes research consisting of identifying and investigating the basic properties of rotating, multijet, high pressure water nozzles, and how particular designs and modes of operation affect such things as stripping rate, standoff distance and completeness of coverage. The study involved computer simulations, an extensive literature review, and experimental studies of different nozzle designs.

Description: An integrated moon program has often been proposed as a logical next step for today's space efforts. In the context of preparing for the possibility of launching a moon program, the European Space Agency is currently conducting an internal study effort which is focusing on the assessment of key technologies. Current thinking has this moon program organized into four phases. Phase 1 will deal with lunar resource exploration. The goal would be to produce a complete chemical inventory of the moon, including oxygen, water, other volatiles, carbon, silicon, and other resources. Phase 2 will establish a permanent robotic presence on the moon via a number of landers and surface rovers. Phase 3 will extend the second phase and concentrate on the use and exploitation of local lunar resources. Phase 4 will be the establishment of a first human outpost. Some preliminary work such as the building of the outpost and the installation of scientific equipment will be done by unmanned systems before a human crew is sent to the moon.

Description: A new model for the two-point vector stream function correlation has been developed using tensor invariant arguments and evaluated by the comparison of model predictions with DNS data for incompressible homogeneous turbulent shear flow. This two-point vector stream function model correlation can then be used to calculate the two-point velocity correlation function and other quantities useful in turbulence modeling. The model assumes that the two-point vector stream function correlation can be written in terms of the separation vector and a new tensor function that depends only on the magnitude of the separation vector. The model has a single free model coefficient, which has been chosen by comparison with the DNS data. The relative error of the model predictions of the two-point vector stream function correlation is only a few percent for a broad range of the model coefficient. Predictions of the derivatives of this correlation, which are of interest in turbulence modeling, may not be this accurate.

Description: A non-local model is presented for approximating the pressure diffusion in calculations of turbulent free shear and boundary layer flows. It is based on the solution of an elliptic relaxation equation which enables local diffusion sources to be distributed over lengths of the order of the integral scale. The pressure diffusion model was implemented in a boundary layer code within the framework of turbulence models based on both the kappa-epsilon-(bar)upsilon(exp 2) system of equations and the full Reynolds stress equations. Model computations were performed for mixing layers and boundary layer flows. In each case, the pressure diffusion model enabled the well-known free-stream edge singularity problem to be eliminated. There was little effect on near-wall properties. Computed results agreed very well with experimental and DNS data for the mean flow velocity, the turbulent kinetic energy, and the skin-friction coefficient.

Description: The structure of the separated flow region over a 20 degree compression corner at a free-stream Mach number of 2.84 is investigated computationally using a Reynolds averaged Navier Stokes (R.A.N.S.) solver and kappa-epsilon model. At this Mach number and ramp angle, a steady-state recirculation region of order delta(sub o) is observed, with onset of a 'plateau' in the wall pressure distribution near the corner. At lower ramp angles, separation is negligible, while at an angle of 24 degrees, separation regions of length 2 delta(sub o) are expected. Of interest here is the response of the mathematical model to inclusion of the pressure dilatation term for turbulent kinetic energy. Compared with the experimental data of Smits and Muck (1987), steady-state computations show improvement when the pressure dilatation term is included. Unsteady computations, using both unforced and then forced inlet conditions, did not predict the oscillatory motion of the separation bubble as observed in laboratory experiments. An analysis of the separation bubble oscillation and the turbulent boundary layer (T.B.L.) frequencies for this flow suggests that the bubble oscillations are of nearly the same order as the turbulent frequencies, and therefore difficult for the model to separate and resolve.

Description: A direct numerical approach has been developed to simulate supersonic turbulent boundary layers. The mean flow quantities are obtained by solving the parabolized Reynolds-averaged Navier-Stokes equations (globally). Fluctuating quantities are computed locally with a temporal direct numerical simulation approach, in which nonparallel effects of boundary layers are partially modeled. Preliminary numerical results obtained at the free-stream Mach numbers 3, 4.5, and 6 with hot-wall conditions are presented. Approximately 5 million grid points are used in all three cases. The numerical results indicate that compressibility effects on turbulent kinetic energy, in terms of dilatational dissipation and pressure-dilatation correlation, are small. Due to the hot-wall conditions the results show significant low Reynolds number effects and large streamwise streaks. Further simulations with a bigger computational box or a cold-wall condition are desirable.

Description: A major problem in modeling of turbulent supersonic flows is the correct assessment of viscous-inviscid interaction problems. Of particular interest is the interaction of boundary layers with shocks. Present turbulence models give in most cases unsatisfactory results in the region of rapid distortion and in the separation region (if one is present) in particular with regard to mean flow profiles and turbulence quantities. The objective of the present work is the direct numerical simulation of shock boundary layer interaction. This report summarizes the first phase during which a numerical method suitable for this problem has been developed and a computer code has been written and tested.

Description: Since high resolution numerical simulations of turbulent flows, or experiments at high Reynolds numbers, represent a substantial investment in resources, the estimation of the minimum resolution required for the study of a given property has been the subject of continued interest. Different properties require, in general, different resolutions, and the present paper is dedicated to the requirements for the measurement of the probability distribution functions of the velocity gradients and, in particular, of their low order moments. The deviation of these quantities from the values corresponding to a Gaussian distribution was one of the first indications of the presence of Reynolds number-dependent intermittency and has been the object of recent interest as numerical simulations have become able to explore the distribution of gradients in the low Reynolds number range, while new experiments have extended the range to increasingly high Reynolds numbers. We will use progressive filtering of the results of numerical simulations of isotropic turbulence as a model for the effect of a sensor of finite size. The numerical issues will be addressed first to ensure that the simulations are fully resolved from the point of view of the velocity gradients. This will also give us an estimate for the numerical resolution required for the different quantities.

Description: The MSU compact four-wire transverse vorticity probe permits omega(sub z)(t) measurements in a nominally 1 sq mm domain. Note that a conventional coordinate system is used with x and y in the streamwise and normal directions respectively. The purpose of this investigation was to acquire time series data in the same access port at the ceiling of the 80 ft x 120 ft wind tunnel (NASA Ames Research Center) as earlier used by the Wallace group from the University of Maryland and to compare the present results with those of the three-component vorticity probe used in that earlier study.

Description: Boundary layer flow transition has long been suggested as a potential noise source in both marine (sonar-dome self noise) and aeronautical (aircraft cabin noise) applications, owing to the highly transient nature of process. The design of effective noise control strategies relies upon a clear understanding of the source mechanisms associated with the unsteady flow dynamics during transition. Due to formidable mathematical difficulties, theoretical predictions either are limited to early linear and weakly nonlinear stages of transition, or employ acoustic analogy theories based on approximate source field data, often in the form of empirical correlation. In the present work, an approach which combines direct numerical simulation of the source field with the Lighthill acoustic analogy is utilized. This approach takes advantage of the recent advancement in computational capabilities to obtain detailed information about the flow-induced acoustic sources. The transitional boundary layer flow is computed by solving the incompressible Navier-Stokes equations without model assumptions, thus allowing a direct evaluation of the pseudosound as well as source functions, including the Lighthill stress tensor and the wall shear stress. The latter are used for calculating the radiated pressure field based on the Curle-Powell solution of the Lighthill equation. This procedure allows a quantitative assessment of noise source mechanisms and the associated radiation characteristics during transition from primary instability up to the laminar breakdown stage. In particular, one is interested in comparing the roles played by the fluctuating volume Reynolds stress and the wall-shear-stresses, and in identifying specific flow processes and structures that are effective noise generators.

Description: Turbulence modeling continues to be one of the most difficult problems in fluid mechanics. Existing prediction methods are well developed for certain classes of simple equilibrium flows, but are still not entirely satisfactory for a large category of complex non-equilibrium flows found in engineering practice. Direct and large-eddy simulation (LES) approaches have long been believed to have great potential for the accurate prediction of difficult turbulent flows, but the associated computational cost has been prohibitive for practical problems. This remains true for direct simulation but is no longer clear for large-eddy simulation. Advances in computer hardware, numerical methods, and subgrid-scale modeling have made it possible to conduct LES for flows or practical interest at Reynolds numbers in the range of laboratory experiments. The objective of this work is to apply ES and the dynamic subgrid-scale model to the flow of a boundary layer over a concave surface.

Description: Previously the theoretical development leading to the dynamic localization model (DLM) for large-eddy simulation (LES) was presented. The method has been successfully applied to isotropic turbulence, channel flow, and the flow over a backward-facing step. Here we apply the model to the computation of the temporally developing place wake. The two main objectives of this project are: (1) Use the model to perform an LES of a time developing plane wake and compare the results with direction numerical simulation (DNS) data to see if important statistical measures can be readily predicted, and to provide a relative evaluation of the several versions of the model in terms of predictive capability and cost; and (2) If the tests in (1) show that the model generates reliable predictions, then use the LES to study various aspects of the physics of turbulent wakes and mixing layers.

Description: The modeling of subgrid scales in large-eddy simulation (LES) has been rationalized by the introduction of the dynamic localization procedure. This method allows one to compute rather than prescribe the unknown coefficients in the subgrid-scale model. Formally, the LES equations are supposed to be obtained by applying to the Navier-Stokes equations a 'grid filter' operation. Though the subgrid stress itself is unknown, an identity between subgrid stresses generated by different filters has been derived. Although preliminary tests of the Dynamic Localization Model (DLM) with k-equation have been satisfactory, the use of a negative eddy viscosity to describe backscatter is probably a crude representation of the physics of reverse transfer of energy. Indeed, the model is fully deterministic. Knowing the filtered velocity field and the subgrid-scale energy, the subgrid stress is automatically determined. We know that the LES equations cannot be fully deterministic since the small scales are not resolved. This stems from an important distinction between equilibrium hydrodynamics and turbulence. In equilibrium hydrodynamics, the molecular motions are also not resolved. However, there is a clear separation of scale between these unresolved motions and the relevant hydrodynamic scales. The result of molecular motions can then be separated into an average effect (the molecular viscosity) and some fluctuations. Due to the large number of molecules present in a box with size of the order of the hydrodynamic scale, the ratio between fluctuations and the average effect should be very small (as a result of the 'law of large numbers'). For that reason, the hydrodynamic balance equations are usually purely deterministic. In turbulence, however, there is no clear separation of scale between small and large eddies. In that case, the fluctuations around a deterministic eddy viscosity term could be significant. An eddy noise would then appear through a stochastic term in the subgrid-scale model and could be the source of backscatter.

Description: The elliptic relaxation model of Durbin (1993) for modeling near-wall turbulence using second moment closures (SMC) is compared to DNS data for a channel flow at Re(sub t) = 395. The agreement for second order statistics and even the terms in their balance equation is quite satisfactory, confirming that very little viscous effects (via Kolmogoroff scales) need to be added to the high Reynolds versions of SMC for near-wall-turbulence. The essential near-wall feature is thus the kinematic blocking effect that a solid wall exerts on the turbulence through the fluctuating pressure, which is best modeled by an elliptic operator. Above the transition layer, the effect of the original elliptic operator decays rapidly, and it is suggested that the log-layer is better reproduced by adding a non-homogeneous reduction of the return to isotropy, the gradient of the turbulent length scale being used as a measure of the inhomogeneity of the log-layer. The elliptic operator was quite easily applied to the non-linear Craft & Launder pressure-strain model yielding an improved distinction between the spanwise and wall normal stresses, although at higher Reynolds number (Re) and away from the wall, the streamwise component is severely underpredicted, as well as the transition in the mean velocity from the log to the wake profiles. In this area a significant change of behavior was observed in the DNS pressure-strain term, entirely ignored in the models.

Description: This work carries earlier finite-difference calculations of the Widnall instability of vortex rings into the late non-linear stage. Plots of energy in azimuthal Fourier modes indicate that low-order modes dominate at large times; their structure and dynamics remain unexplored, however. An attempt was made to calculate the acoustic signal using the theory of Mohring (1978), valid for unbounded flow. This theory shows that only low-order azimuthal modes contribute to the sound. As a check on the effects of axial periodicity and a slip wall at large radius imposed by the numerical scheme, the acoustic integrals were also computed in a truncated region. Half of the terms contributing to the sound have large differences between the two regions, and the results are therefore unreliable. The error is less severe for a contribution involving only the m = 2 mode, and its low frequency is consistent with a free elliptic bending wave on a thin ring.

Description: The spectra of the pressure, and other higher-order quantities including the dissipation, the enstrophy, and the square of the longitudinal velocity derivative are computed using data obtained from direct numerical simulation of homogeneous isotropic turbulence at Taylor-Reynolds numbers R(sub lambda) in the range 38 - 170. For the pressure spectra we find reasonable collapse in the dissipation range (of the velocity spectrum) when scaled in Kolmogorov variables and some evidence, which is not conclusive, for the existence of a k(exp -7/3) inertial range where k = absolute value of K, is the modulus of the wavenumber. The power spectra of the dissipation, the enstrophy, and the square of the longitudinal velocity derivative separate in the dissipation range, but appear to converge together in the short inertial range of the simulations. A least-squares curve-fit in the dissipation range for one value of R(sub lambda) = 96 gives a form for the spectrum of the dissipation as k(exp 0)exp(-Ck eta), for k(eta) greater than 0.2, where eta is the Kolmogorov length and C is approximately equal to 2.5.

Description: Panel methods are commonly used in computational fluid dynamics for the solution of potential flow problems. The methods are a numerical technique based on the surface distribution of singularity elements. The solution is the process of finding the strength of the singularity elements distributed over the body's surface. This process involves the solution of the matrix problem Pq = p' for a set of unknowns q. The Fast Multipole Method is used to directly compute q without using matrix solvers. The algorithm works in O(N) time for N points, a great improvement over standard matrix solvers. In panel methods, the surface of a body is divided into a series of quadrilateral panels. The methods involve the computation of the influence of all other panels on each individual panel. The influence is based on the surface distribution, though this can be approximated by the area for distant panels. An alternative approximation, though with arbitrary accuracy, is to develop a multipole expansion about the center of the panel to describe the effect of a given panel on distant points in space. The expansion is based on the moments of the panel, thus allow the use of various surface distributions without changing the basic algorithm, just the computation of the various moments. The expansions are then manipulated in a tree walk to develop Taylor series expansions about a point in space which describe the effect of all distant panels on any point within a volume of convergence. The effect of near panels then needs to be computed directly, but the effect of all distant panels can be computed by simply evaluating the resulting expansion. The Fast Multipole Method has been applied to panel methods for the solution of source and doublet distributions. A major feature of the algorithm is that the algorithm does not change to derive the potential and velocity for sources and doublets. The same expansions can be used for both sources and doublets. Since the velocity is related to the potential, and the doublet potential is related to the z-component of the source velocity, all values can be derived from the same expansion by taking a series of partial derivatives. This requires more expansion terms to be kept since terms are lost in the process of taking partial derivatives. Thus to maintain accuracy for the doublet computation, more terms are required than if just evaluating for sources. The resulting Fast Multipole code should then parallelize better than classical panel methods due to the locality of data dependencies found in the Fast Multipole Method. Theoretically the parallelized code should execute in O(log N) time with O(N) processors, though this is not practical. Ongoing work includes implementing the parallel accelerated panel method, including methods to improve the load balancing of the problem by taking advantage of the known geometry of panels, and to encorporate sensitivity analysis into the algorithm.

Description: Asymptotic similarity states at large Reynolds numbers and small Rossby numbers in rotating homogeneous turbulence are investigated using the database obtained from large-eddy simulations of the incompressible Navier-Stokes equations. Previous work has shown that the turbulence kinetic energy and integral length scales are accurately described by simple scaling laws based on the low wavenumbers part of the three-dimensional energy spectrum. The primary interest of the present study is to search for spectrum similarity in the asymptotic state. Four independent energy spectra are defined. It is shown that rescaling of these energy spectra in the asymptotic regime will collapse three out of the four spectra. The spectrum which does not collapse is a function only of the vertical wavenumber and corresponds to two-component motions in the plane normal to the rotation axis. Detailed investigation of the cause of this anomalous behavior reveals the existence of a strong reverse cascade of energy from small-to-large scales of the two-dimensional, two-component motions. This feature of the rotating flow is presumably linked to the lack of a complete similarity state, though further study of this issue is required.

Description: This paper describes a path planning method for planetary rovers to search for paths on planetary surfaces. The planetary rover is required to travel safely over a long distance for many days over unfamiliar terrain. Hence it is very important how planetary rovers process sensory information in order to understand the planetary environment and to make decisions based on that information. As a new data structure for informational mapping, an extended elevation map (EEM) has been introduced, which includes the effect of the size of the rover. The proposed path planning can be conducted in such a way as if the rover were a point while the size of the rover is automatically taken into account. The validity of the proposed methods is verified by computer simulations.

Description: The pressure distribution on a swept wing causes the streamlines at the edge of the boundary layer to be curved. This pressure gradient normal to the external streamline creates a velocity component normal to the external streamline within the boundary layer which is referred to as the crossflow velocity. Because the crossflow velocity profile perpendicular to the wing surface has an inflection point, the profile is unstable. The stationary instability mode takes the form of crossflow vortices. Under these conditions, the boundary layer on the wing is extremely unstable and transition to turbulent flow takes place much closer to the leading edge of the wing than it would on an unswept wing. Higher skin friction drag is associated with turbulent flow, and so better aircraft performance could be obtained if the crossflow could be eliminated One method of controlling crossflow that is being investigated is boundary-layer suction. An extensive airfoil suction experiment in the 8 feet Transonic Pressure Tunnel (TPT) at NASA Langley Research Center will begin late in 1994. Because of the size, complexity, and expense associated with this test, a number of 'risk-reduction' tests are currently being conducted. The 20 x 28 in. Shear Flow Control Tunnel at NASA Langley is being used for some of these tests. Prior to the summer of 1994, a flat plate with a swept leading edge was installed in the 20 x 28 in. tunnel, with a displacement body mounted on the tunnel ceiling that created a pressure distribution on the plate similar to the pressure distribution on a swept wing. The flow over the plate was investigated during the summer of 1994 using a laser Doppler velocimeter (LDV) system. The LDV measurements indicated the possible presence of multiple disturbance modes, a rarely-seen phenomena since, in most tests, one disturbance mode dominates. The possible existence of multiple disturbance modes in the flat plate boundary layer, however, means that the flow in the 20 x 28 in. tunnel is of interest itself, and will be investigated more thoroughly in the future. With a view to these investigations, the boundary layer traverse mechanism in the 20 x 28 in. tunnel was modified to improve its performance, and strain gauges were mounted on the traverse in order to monitor its deflection during a test. Other preliminary work conducted in the 20 x 28 in. tunnel included the use of an infrared camera system. Previous work with this system showed that transition indeed could be detected, but the signal produced by the crossflow vortices was too weak to be detected. It was hoped that spraying the flat plate with naphthalene would augment the heat transfer associated with the crossflow vortices so that they would show up in a IR image; however, experiments showed that this would not work. Another set of tests was conducted in the 20 x 28 in. tunnel to determine the tripping requirements for a set of airfoil-shaped struts that will be used in the 8 feet TPT experiment. Since the Reynolds number associated with these struts is small, a laminar boundary layer would separate early, causing large fluctuations in the flow field. A turbulent boundary layer would remain attached further back, but tripping from laminar to turbulent flow at low Reynolds number is very difficult. However, trip strip configurations were found that should effectively trip the boundary layer at the required conditions. Currently underway is an investigation of the data acquisition requirements for the 8 feet TPT experiment, with the purpose of the finding the minimum amount of data needed to characterize sufficiently the swept-wing boundary layer. This study is being conducted using a numerically generated data set.

Description: The long-term goal of the current research is the study of Large-Eddy Simulation (LES) as a tool for aeroacoustics. New algorithms and developments in computer hardware are making possible a new generation of tools for aeroacoustic predictions, which rely on the physics of the flow rather than empirical knowledge. LES, in conjunction with an acoustic analogy, holds the promise of predicting the statistics of noise radiated to the far-field of a turbulent flow. LES's predictive ability will be tested through extensive comparison of acoustic predictions based on a Direct Numerical Simulation (DNS) and LES of the same flow, as well as a priori testing of DNS results. The method presented here is aimed at allowing simulation of a turbulent flow field that is both simple and amenable to acoustic predictions. A free shear flow is homogeneous in both the streamwise and spanwise directions and which is statistically stationary will be simulated using equations based on the Navier-Stokes equations with a small number of added terms. Studying a free shear flow eliminates the need to consider flow-surface interactions as an acoustic source. The homogeneous directions and the flow's statistically stationary nature greatly simplify the application of an acoustic analogy.

Description: Detailed three-dimensional three-component phase averaged measurements of the spanwise and streamwise vorticity formation and evolution in acoustically forced plane free-shear flows have been obtained. For the first time, phase-averaged measurements of all three velocity components have been obtained in both a mixing layer and a wake on three-dimensional grids, yielding the spanwise and streamwise vorticity distributions without invoking Taylor's hypothesis. Initially, two-frequency forcing was used to phase-lock the roll-up and first pairing of the spanwise vortical structures in a plane mixing layer. The objective of this study was to measure the near-field vortical structure morphology in a mixing layer with 'natural' laminar initial boundary layers. For the second experiment the second and third subharmonics of the fundamental roll-up frequency were added to the previous two-frequency forcing in order to phase-lock the roll-up and first three pairings of the spanwise rollers in the mixing layer. The objective of this study was to determine the details of spanwise scale changes observed in previous time-averaged measurements and flow visualization of unforced mixing layers. For the final experiment, single-frequency forcing was used to phase-lock the Karman vortex street in a plane wake developing from nominally two-dimensional laminar initial boundary layers. The objective of this study was to compare measurements of the three-dimensional structure in a wake developing from 'natural' initial boundary layers to existing models of wake vortical structure.

Description: It was previously shown that the strong vorticity in isotropic turbulence is organized into tubular vortices ('worms') whose properties were characterized through the use of full numerical simulations at several Reynolds numbers. At the time most of the observations were kinematic, and several scaling laws were discovered for which there was no theoretical explanation. In the meantime, further analysis of the same fields yielded new information on the generation of the vortices, and it was realized that even if they had to be formed by stretching, they were at any given moment actually compressed at many points of their axes. This apparent contradiction was partially explained by postulating axial inertial waves induced by the nonuniformity of the vortex cores, which helped to 'spread' the axial strain and allowed the vortices to remain compact even if not uniformly stretched. The existence of such solutions was recently proved numerically. The present report discusses a set of new numerical simulations of isotropic turbulence, and a reanalysis of the old ones, in an effort to prove or disprove the presence of these waves in actual turbulent flows and to understand the dynamics, as opposed to the kinematics, of the vortices.

Description: A method for active fluid flow control based on control theory is discussed. Dynamic programming and fixed point successive approximations are used to accommodate the nonlinear control problem. The long-term goal of this project is to establish an effective method applicable to complex flows such as turbulence and jets. However, in this report, the method is applied to stochastic Burgers equation as an intermediate step towards this goal. Numerical results are compared with those obtained by gradient search methods.

Description: The objective of the present research is the development and application of efficient adaptive numerical algorithms for the study, via direct numerical simulations, of active vortex generators. We are using innovative computational schemes to investigate flows past complex configurations undergoing arbitrary motions. Some of the questions we try to answer are: Can and how may we control the dynamics of the wake? What is the importance of body shape and motion in the active control of the flow? What is the effect of three-dimensionality in laboratory experiments? We are interested not only in coupling our results to ongoing, related experimental work, but furthermore to develop an extensive database relating the above mechanisms to the vortical wake structures with the long-range objective of developing feedback control mechanisms. This technology is very important to aircraft, ship, automotive, and other industries that require predictive capability for fluid mechanical problems. The results would have an impact in high angle of attack aerodynamics and help design ways to improve the efficiency of ships and submarines (maneuverability, vortex induced vibration, and noise).

Description: A challenge for traditional turbulence modeling, based on the Reynolds averaged Navier-Stokes equations, remains the accurate prediction of 'mild', adverse pressure-gradient driven separation from a smooth surface. With this study we want to explore the capability of large-eddy simulation to predict the separation which occurs on the deflected wall of an asymmetric, plane diffuser with opening angle of 10 deg. The flow through the plane diffuser exhibits some additional interesting physical phenomena which make it a challenging test case. In addition to 'mild' separation about halfway down the deflected ramp, the flow is characterized by a small backflow zone with stalled fluid in the rear part of the expanding section. The turbulent flow entering the diffuser is subject to combined adverse and radial pressure gradients stemming from the convex curvature. Finally the flow recovers into a developed, turbulent channel flow in the outlet section. Obi et al. provide measurements of mean flow, Reynolds stresses, and pressure recovery, which were obtained by means of LDV in a wind tunnel. The objective of this study is to investigate whether LES with the standard dynamic model is able to accurately predict the flow in the one-sided diffuser and to explore the resolution requirements and associated costs.

Description: The dynamic subgrid-scale (SGS) model has given good results in the large-eddy simulation (LES) of homogeneous isotropic or shear flow, and in the LES of channel flow, using averaging in two or three homogeneous directions (the DA model). In order to simulate flows in general, complex geometries (with few or no homogeneous directions), the dynamic SGS model needs to be applied at a local level in a numerically stable way. Channel flow, which is inhomogeneous and wall-bounded flow in only one direction, provides a good initial test for local SGS models. Tests of the dynamic localization model were performed previously in channel flow using a pseudospectral code and good results were obtained. Numerical instability due to persistently negative eddy viscosity was avoided by either constraining the eddy viscosity to be positive or by limiting the time that eddy viscosities could remain negative by co-evolving the SGS kinetic energy (the DLk model). The DLk model, however, was too expensive to run in the pseudospectral code due to a large near-wall term in the auxiliary SGS kinetic energy (k) equation. One objective was then to implement the DLk model in a second-order central finite difference channel code, in which the auxiliary k equation could be integrated implicitly in time at great reduction in cost, and to assess its performance in comparison with the plane-averaged dynamic model or with no model at all, and with direct numerical simulation (DNS) and/or experimental data. Other local dynamic SGS models have been proposed recently, e.g., constrained dynamic models with random backscatter, and with eddy viscosity terms that are averaged in time over material path lines rather than in space. Another objective was to incorporate and test these models in channel flow.

Description: Two projects are described in this report. The first involves assessing turbulence models in separated flow. The second addresses the anomalous behavior of certain turbulence models in stagnation point flow. The primary motivation for developing turbulent transport models is to provide tools for computing non-equilibrium, or complex, turbulent flows. Simple flows can be analyzed using data correlations or algebraic eddy viscosities, but in more complicated flows such as a massively separated boundary layer, a more elaborate level of modeling is required. It is widely believed that at least a two-equation transport model is required in such cases. The transport equations determine the evolution of suitable velocity and time-scales of the turbulence. The present study included assessment of second-moment closures in several separated flows, including sharp edge separation; smooth wall, pressure driven separation; and unsteady vortex shedding. Flows with mean swirl are of interest for their role in enhancing mixing both by turbulent and mean motion. The swirl can have a stabilizing effect on the turbulence. An axi-symmetric extension to the INS-2D computer program was written adding the capability of computing swirling flow. High swirl can produce vortex breakdown on the centerline of the jet and it occurs in various combustors.

Description: A model for the effects of rotation on turbulence is proposed and tested. These effects which influence mainly the rate of turbulence decay are modeled in a modified turbulent energy dissipation rate equation that has explicit dependence on the mean rotation rate. An appropriate definition of the rotation rate derived from critical point theory and based on the invariants of the deformation tensor is proposed. The modeled dissipation rate equation is numerically well behaved and can be used in conjunction with any level of turbulence closure. The model is applied to the two-equation kappa-epsilon turbulence model and is used to compute separated flows in a backward-facing step and an axisymmetric swirling coaxial jets into a sudden expansion. In general, the rotation modified dissipation rate model shows some improvements over the standard kappa-epsilon model.

Description: The long-time, asymptotic state of rotating homogeneous shallow-water equations is investigated. Our analysis is based on long-time averaged rotating shallow-water equations describing interactions of large-scale, horizontal, two-dimensional motions with surface inertial-gravity waves field for a shallow, uniformly rotating fluid layer. These equations are obtained in two steps: first by introducing a Poincare/Kelvin linear propagator directly into classical shallow-water equations, then by averaging. The averaged equations describe interaction of wave fields with large-scale motions on time scales long compared to the time scale 1/f(sub o) introduced by rotation (f(sub o)/2-angular velocity of background rotation). The present analysis is similar to the one presented by Waleffe (1991) for 3D Euler equations in a rotating frame. However, since three-wave interactions in rotating shallow-water equations are forbidden, the final equations describing the asymptotic state are simplified considerably. Special emphasis is given to a new conservation law found in the asymptotic state and decoupling of the dynamics of the divergence free part of the velocity field. The possible rising of a decoupled dynamics in the asymptotic state is also investigated for homogeneous turbulence subjected to a background rotation. In our analysis we use long-time expansion, where the velocity field is decomposed into the 'slow manifold' part (the manifold which is unaffected by the linear 'rapid' effects of rotation or the inertial waves) and a formal 3D disturbance. We derive the physical space version of the long-time averaged equations and consider an invariant, basis-free derivation. This formulation can be used to generalize Waleffe's (1991) helical decomposition to viscous inhomogeneous flows (e.g. problems in cylindrical geometry with no-slip boundary conditions on the cylinder surface and homogeneous in the vertical direction).

Description: As part of the subsonic transport high-lift program, flight experiments are being conducted using NASA Langley's B737-100 to measure the flow characteristics of the multi-element high-lift system at full-scale high-Reynolds-number conditions. The instrumentation consists of hot-film anemometers to measure boundary-layer states, an infra-red camera to detect transition from laminar to turbulent flow, Preston tubes to measure wall shear stress, boundary-layer rakes to measure off-surface velocity profiles, and pressure orifices to measure surface pressure distributions. The initial phase of this research project was recently concluded with two flights on July 14. This phase consisted of a total of twenty flights over a period of about ten weeks. In the coming months the data obtained in this initial set of flight experiments will be analyzed and the results will be used to finalize the instrumentation layout for the next set of flight experiments scheduled for Winter and Spring of 1995. The main goal of these upcoming flights will be: (1) to measure more detailed surface pressure distributions across the wing for a range of flight conditions and flap settings; (2) to visualize the surface flows across the multi-element wing at high-lift conditions using fluorescent mini tufts; and (3) to measure in more detail the changes in boundary-layer state on the various flap elements as a result of changes in flight condition and flap deflection. These flight measured results are being correlated with experimental data measured in ground-based facilities as well as with computational data calculated with methods based on the Navier-Stokes equations or a reduced set of these equations. Also these results provide insight into the extent of laminar flow that exists on actual multi-element lifting surfaces at full-scale high-life conditions. Preliminary results indicate that depending on the deflection angle, the slat and flap elements have significant regions of laminar flow over a wide range of angles of attack. Boundary-layer transition mechanisms that were observed include attachment-line contamination on the slat and inflectional instability on the slat and fore flap. Also, the results agree fairly well with the predictions reported in a paper presented at last year's AIAA Fluid Dynamics Conference. The fact that extended regions of laminar flow are shown to exist on the various elements of the high-lift system raises the question what the effect is of loss of laminar flow as a result of insect contamiantion, rain or ice accumulation on high-life performance.

Description: Existing aerodynamic design methods have generally concentrated on the optimization of airfoil or wing shapes to produce a minimum drag while satisfying some basic constraints such as lift, pitching moment, or thickness. Since the minimization of drag almost always precludes the existence of separated flow, the evaluation and validation of these design methods for their robustness and accuracy when separated flow is present has not been aggressively pursued. However, two new applications for these design tools may be expected to include separated flow and the issues of aerodynamic design with this feature must be addressed. The first application of the aerodynamic design tools is the design of airfoils or wings to provide an optimal performance over a wide range of flight conditions (multipoint design). While the definition of 'optimal performance' in the multipoint setting is currently being hashed out, it is recognized that given a wide enough range of flight conditions, it will not be possible to ensure a minimum drag constraint at all conditions, and in fact some amount of separated flow (presumably small) may have to be allowed at the more demanding flight conditions. Thus a multipoint design method must be tolerant of the existence of separated flow and may include some controls upon its extent. The second application is in the design of wings with extended high speed buffet boundaries of their flight envelopes. Buffet occurs on a wing when regions of flow separation have grown to the extent that their time varying pressures induce possible destructive effects upon the wing structure or adversely effect either the aircraft controllability or the passenger comfort. A conservative approach to the expansion of the buffet flight boundary is to simply expand the flight envelope of nonseparated flow under the assumption that buffet will also thus be alleviated. However, having the ability to design a wing with separated flow and thus to control the location, extent, and severity of the separated flow regions may allow aircraft manufacturers to gain an advantage in the early design stages of an aircraft, when configuration changes are relatively inexpensive to make. Continuing the work begun last year, an airfoil design package has been modified to provide some control over the existence and extent of flow separation. This package consists of a 2-D Navier-Stokes flow solver which is coupled to the CDISC (constrained direct/iterative surface curvature) design method. The first modification is a prediction method for determining whether separation is likely based solely upon a given pressure distribution. If separation is predicted but is undesirable, the new routines will modify the pressure distribution to alleviate the problem. This new pressure distribution is then used in the design method to generate a new aerodynamic shape. Since separation may be acceptable in some cases, particularly if the separation does not extend to the trailing edge, another added logic estimates the extent of separation based upon a correlation with calculated separated flow cases. If the flow behind a shock induced separation is not predicted to reattach before the trailing edge, the logic weakens the shock strength and otherwise alters the pressure distribution in order to promote reattachment. This later addition is as yet unreliable due to secondary separation effects, but additional work is being pursued to improve the method.