ALBERT

Description: Helmet streamers on the sun were observed to be the site of coronal mass ejections, dynamic events that eject coronal plasma and magnetic fields into the solar wind. A two dimensional (azimuthally symmetric) helmet streamer configuration was developed by computing solutions of the time dependent magnetohydrodynamic (MHD) equations, for a specified magnetic flux distribution on the sun. The helmet streamer is not symmetric about the equator. The evolution of the configuration, when differential rotation is applied, was investigated. It was found that after many rotations the configuration does not reach a steady state, but disrupts recurrently with the ejection of a plasmoid. These results suggest that differential rotation may be one of the mechanisms by which mass ejections are initiated.

Description: The ground-based observing facilities of the National Solar Observatory (NSO) are reviewed from the perspective of joint observations with SOHO. A specific proposal is presented for observations of the HE-I 1083.0 nm line with the NASA/NSO spectromagnetograph and He 10830 video filtergraph/magnetograph in coordination with ultraviolet sensitive instruments on SOHO. The first task will be to look for associations of low-temperature transition-region lines with He 1083 nm absorption to investigate Andretta's conjecture, i.e. that the He 1083 nm line is formed in two layers where extreme ultraviolet radiation produced both in the low-temperature transition region (the upper layer) and in the surrounding corona products - a lower layer of absorption in the upper chromosphere.

Description: The white light coronograph (WLC) on Skylab provided an opportunity to study the corona at high spatial and temporal resolution. The spatial resolution of the instrument was approximately 25 cm with images taken approximately one per min. One set of images taken over a 10 min period was digitized, providing ten high spatial resolution images for analysis. The progress in data processing techniques available at the time was not sufficient to permit a reliable study of the fine structure in these images. Using current techniques an investigation of the sizes and lifetimes of the smallest scale features in the data was carried out. A preliminary analysis of an area between 2 and 3 Ro was completed. The results show that very narrow rays extend from at least 2 to 3 Ro. The narrowest of these rays has a thickness of approximately 75 cm. The contrast is so low that they are very close to the noise limit of the data. Most of the rays observed become unrecognizable after 10 min, although some remain visible over the entire time. Some notion seems to be detectable in the fine structure rays, but analysis of more frames will be needed to quantify these results.

Description: Filaments, flare sprays, prominences and 'post-flare' loops are familiar to H alpha observers in their frequent appearances 'in absorption', dark against the chromospheric background or plages. Observations of the X-ray corona are generally interpreted as due to emission via optically thin thermal bremsstrahlung. Several cases of X-ray coronal structures in Yohkoh images, due to high opacity, absorbing matter in coronograph loops, are presented. The presence of the absorbing matter, mixed with emitting matter, complicates inference of physical parameters such as emission measures in X-ray sources. In the case of well defined features, absorption provides an opportunity to infer density. Quantitative estimates of the attenuation due to the absorption in example features are presented.

Description: The causal association of major solar particle events seen at earth with coronal mass ejections (CME's), and not with solar flares, is discussed. Evidence that led to the demise of the flare dominated paradigm for major solar energetic particle events are described. The possibility of distinguishing particles from impulsive and gradual events using only observations is described. Particle acceleration at the CME level is discussed. Multi-spacecraft observations of CME events are described. Concerning the interplanetary CME, bidirectional proton events are discussed. Conclusions from progress in understanding the characteristics of solar energetic particles and their relation to the physical mechanisms of acceleration are given.

Description: Two dimensional magnetohydrodynamic simulations of the distortion of a magnetic flux tube, accelerated through ambient solar wind plasma, are presented. Vortices form on the trailing edge of the flux tube, and couple strongly to its interior. If the flux tube azimuthal field is weak, it deforms into an elongated banana-like shape after a few Alfven transit times. A significant azimuthal field component inhibits this distortion. In the case of magnetic clouds in the solar wind, it is suggested that the shape observed at 1 AU was determined by distortion of the cloud in the inner heliosphere. Distortion of the cloud beyond 1 AU takes many days. It is estimated that effective drag coefficients slightly greater than unity are appropriate for modeling flux tube propagation. Synthetic magnetic field profiles as would be seen by a spacecraft traversing the cloud are presented.

Description: Since radio propagation measurements using either natural or spacecraft radio signals are used for probing the solar wind in the vicinity of the sun, they represent a key tool for studying the interplanetary consequences of solar structure and dynamic phenomena. New information on the near sun consequences was obtained from radio scintillation observations of coherent spacecraft signals. The results covering density fluctuations, fractional density fluctuations, coronal streamers, heliospheric current sheets, coronal mass ejections and interplanetary shocks are reviewed. A joint ICE S-band (13 cm wavelength) Doppler scintillation measurement with the SOHO white-light coronograph (LASCO) is described.

Description: The test capabilities of the Stability Wind Tunnel of the Virginia Polytechnic Institute and State University are described, and calibrations for curved and rolling flow techniques are given. Oscillatory snaking tests to determine pure yawing derivatives are considered. Representative aerodynamic data obtained for a current fighter configuration using the curved and rolling flow techniques are presented. The application of dynamic derivatives obtained in such tests to the analysis of airplane motions in general, and to high angle of attack flight conditions in particular, is discussed.

Description: The scaling properties of a time series of Doppler images obtained in good visibility conditions are studied. A 28 cm vacuum telescope and a vacuum spectroheliograph in video spectra-spectroheliograph mode, are used. Sixty line-of-sight Doppler images of an area of the quiet sun are investigated. They were taken at 60 sec intervals over a one hour span and have a 2 arcsec resolution. After the removal of the five-minute oscillations, the time-spatial spectrum is calculated. To study the turbulence of photospheric flows, two scaling parameters in the spectra, are estimated: the exponent of the spatial part of the power spectrum, and the exponent governing the scaling of time correlations. The implied diffusive behavior is discussed. This includes the estimation of a diffusion coefficient and the type of diffusion involved.

Description: The nonlocal non-diffusive transport of passive scalars in turbulent magnetohydrodynamic (MHD) convection is investigated using transilient matrices. These matrices describe the probability that a tracer particle beginning at one position in a flow will be advected to another position after some time. A method for the calculation of these matrices from simulation data which involves following the trajectories of passive tracer particles and calculating their transport statistics, is presented. The method is applied to study the transport in several simulations of turbulent, rotating, three dimensional compressible, penetrative MDH convection. Transport coefficients and other diagnostics are used to quantify the transport, which is found to resemble advection more closely than diffusion. Some of the results are found to have direct relevance to other physical problems, such as the light element depletion in sun-type stars. The large kurtosis found for downward moving particles at the base of the convection zone implies several extreme events.

Description: Long uninterrupted sequences of solar magnetograms from the global oscillations network group (GONG) network and from the solar and heliospheric observatory (SOHO) satellite will provide the opportunity to study the proper motions of magnetic features. The possible use of multiscale regularization, a scale-recursive estimation technique which begins with a prior model of how state variables and their statistical properties propagate over scale. Short magnetogram sequences are analyzed with the multiscale regularization algorithm as applied to optical flow. This algorithm is found to be efficient, provides results for all the spatial scales spanned by the data and provides error estimates for the solutions. It is found that the algorithm is less sensitive to evolutionary changes than correlation tracking.

Description: Ulysses has collected data between 1 and 5 AU during, and just following solar maximum, when the heliospheric current sheet (HCS) can be thought of as reaching its maximum tilt and being subject to the maximum amount of turbulence in the solar wind. The Ulysses solar wind plasma instrument measures the vector velocity and can be used to estimate the flow speed and direction in turbulent 'eddies' in the solar wind that are a fraction of an astronomical unit in size and last (have either a turnover or dynamical interaction time of) several hours to more than a day. Here, in a simple exercise, these solar wind eddies at the HCS are characterized using Ulysses data. This character is then used to define a model flow field with eddies that is imposed on an ideal HCS to estimate how the HCS will be deformed by the flow. This model inherently results in the complexity of the HCS increasing with heliocentric distance, but the result is a measure of the degree to which the observed change in complexity is a measure of the importance of solar wind flows in deforming the HCS. By comparison with randomly selected intervals not located on the HCS, it appears that eddies on the HCS are similar to those elsewhere at this time during the solar cycle, as is the resultant deformation of the interplanetary magnetic field (IMF). The IMF deformation is analogous to what is often termed the 'random walk' of interplanetary magnetic field lines.

Description: The Galileo Ultravilet Spectrometer Experiment (UVS) obtained a partial celestial sphere map of interplanetary Lyman-alpha (IP L alpha) on 13-14 December 1990 during the first Earth encounter. The Galileo spacecraft was near the downwind axis of the local interstellar medium flow. These UVS measurements sampled the downwind, anti-sunward hemisphere. The data were modeled using a hot model of the interplanetary hydrogen density distribution with the goal of studying multiple scattering effects in the inner solar system. The derived ratio in the downwind direction of the observed brightness and a single scattering model brightness, both normalized to unity in the upwind direction, is 1.82 +/- 0.2. This brightness ratio requires a multiple scattering correction which is 36% larger than can be accounted for by theoretical calculations. The hot model may require: (1) a temperature perturbation of the interstellar wind velocity distribution or (2) an additional downstream source of interplanetary hydrogen. However, a more likely exlanation which affects the hot model is the latitude dependence of the radiation pressure. This dependence, based on the known solar L alpha flux latitude variation at solar maximum, causes a downwind brightness enhancement by preferential focusing of H-atoms with trajectory planes containing the solar poles. This result implies that radiation pressure near the solar poles is nearly independent of solar cycle and is insufficient to lead to a net repulsion of hydrogen atoms by the sun, as can occur near the ecliptic plane during the solar maximum. In addition, the UVS performed 13 observations of IP L alpha while in cruise between Venus and the Earth in 3 directions fixed in ecliptic coordinates.

Description: In this paper we generalize earlier gasdynamic analyses of the motion of the heliospheric termination shock in response to upstream disturbances (Barnes, 1993, 1994; Naidu and Barnes, 1994), to include magnetohydrodynamic (MHD) phenomena. We assume that the termination shock is a strong, perpendicular shock and that the initial upstream disturbance is a tangential discontinuity. The resulting configuration after the interaction is very similar to that in the gasdynamic models after an interaction with a contact discontinuity or interplanetary shock, and for an increase (decrease) in dynamic pressure consists of an outward (inward) propagating termination shock and an outward propagating shock (MHD rarefraction wave) that carries the signal of the disturbance into the far downstream plasma. The plasma immediately behind the new termination shock is separated from the downstream signal by a tangential discontinuity. The results of the model show that the speed of the new termination shock depends mainly on the magnitude of the change in dynamic pressure and are typically of order approximately 100 km/s, comparable to the results of the gasdynamic models.

Description: Several physical and observational effects contribute to the significant imbalances of magnetic flux that are often observed in active regions. We consider an effect not previously treated: the influence of electric currents in the photosphere. Electric currents can cause a line-of-sight flux imbalance because of the directionality of the magnetic field they produce. Currents associated with magnetic flux tubes produce larger imbalances than do smoothly-varying distributions of flux and current. We estimate the magnitude of this effect for current densities, total currents, and magnetic geometry consistent with observations. The expected imbalances lie approximately in the range 0-15%, depending on the character of the current-carying fields and the angle from which they are viewed. Observationally, current-induced flux imbalances could be indicated by a statistical dependence of the imbalance on angular distance from disk center. A general study of magnetic flux balance in active regions is needed to determine the relative importance of other- probably larger- effects such as dilute flux (too weak to measure or rendered invisible by radiative transfer effects), merging with weak background fields, and long-range connections between active regions.

Description: Measurements of wing buffeting, using root strain gages, were made in the NASA Langley 0.3 m cryogenic wind tunnel to refine techniques which will be used in larger cryogenic facilities such as the United States National Transonic Facility (NTF) and the European Transonic Wind Tunnel (ETW). The questions addressed included the relative importance variations in frequency parameter and Reynolds number, the choice of model material (considering both stiffness and damping) and the effects of static aeroelastic distortion. The main series of tests was made on three half models of slender 65 deg delta wings with a sharp leading edge. The three delta wings had the same planform but widely differing bending stiffnesses and frequencies (obtained by varying both the material and the thickness of the wings). It was known that the steady flow on this configuration would be insensitive to variations in Reynolds number. On this wing at vortex breakdown the spectrum of the unsteady excitation is unusual, having a sharp peak at particular frequency parameter. Additional tests were made on one unswept half-wing of aspect ratio 1.5 with an NPL 9510 aerofoil section, known to be sensitive to variations in Reynolds number at transonic speeds. The test Mach numbers were M = 0.21 and 0.35 for the delta wings and to M = 0.30 for the unswept wing. On this wing the unsteady excitation spectrum is fairly flat (as on most wings). Hence correct representation of the frequency parameter is not particularly important.

Description: The compressible dynamic stall flowfield over a NACA 0012 airfoil transiently pitching from 0 to 60 deg at a constant rate under compressible flow conditions has been studied using real-time interferometry. A quantitative description of the overall flowfield, including the finer details of dynamic stall vortex formation, growth, and the concomitant changes in the airfoil pressure distribution, has been provided by analyzing the interferograms. For Mach numbers above 0.4, small multiple shocks appear near the leading edge and are present through the initial stages of dynamic stall. Dynamic stall was found to occur coincidentally with the bursting of the separation bubble over the airfoil. Compressibility was found to confine the dynamic stall vortical structure closer to the airfoil surface. The measurements show that the peak suction pressure coefficient drops with increasing freestream Mach number, and also it lags the steady flow values at any given angle of attack. As the dynamic stall vortex is shed, an anti-clockwise vortex is induced near the trailing edge, which actively interacts with the post-stall flow.

Description: The effect of the porous leading edge of an airfoil on the blade-vortex interaction noise, which dominates the far-field acoustic spectrum of the helicopter, is investigated. The thin-layer Navier-Stokes equations are solved with a high-order upwind-biased scheme and a multizonal grid system. The Baldwin-Lomax turbulence model is modified for considering transpiration on the surface. The amplitudes of the propagating acoustic wave in the near field are calculated directly from the computation. The porosity effect on the surface is modeled in two ways: (1) imposition of prescribed transpiration velocity distribution and (2) calculation of transpiration velocity distribution by Darcy's law. Results show leading-edge transpiration can suppress pressure fluctuations at the leading edge during blade-vortex interaction and consequently reduce the amplitude of propagating noise by 30% at a maximum in the near field.

Description: A method has been developed for calculating the viscous flow about airfoils with and without deflected flaps at -90 deg incidence. This method provides for the solution of the unsteady incompressible Navier-Stokes equations by means of an implicit technique. The solution is calculated on a body-fitted computational mesh using a staggered-grid method. The vorticity is defined at the node points, and the velocity components are defined at the mesh-cell sides. The staggered-grid orientation provides for accurate representation of vorticity at the node points and the continuity equation at the mesh-cell centers. The method provides for the noniterative solution of the flowfield and satisfies the continuity equation to machine zero at each time step. The method is evaluated in terms of its stability to predict two-dimensional flow about an airfoil at -90-deg incidence for varying Reynolds number and laminar/turbulent models. The variations of the average loading and surface pressure distribution due to flap deflection, Reynolds number, and laminar or turbulent flow are presented and compared with experimental results. The comparisom indicate that the calculated drag and drag reduction caused by flap deflection and the calculated average surface pressure are in excellent agreement with the measured results at a similar Reynolds number.

Description: A semianalytic method is derived for dealing simultaneously with large numbers of linear stellar oscillation modes trapped in a cavity (a shell) of fluid which is rotating and convecting. A simple generalization of mixing-length theory shows how convection is modulated by weak rotational effects and by the horizontal wind fields of linear r-mode oscillations. The modulated convection is then used to compute the energy lost to turbulent viscosity by a family of nondegenerate oscillations. Viscosity terms of fourth degree in the wind shear can be included if they are a perturbation affecting only a small portion of the r-mode. Viscous energy loss strenghthens convection in a narrow layer near the base of the H and He ionization zone. In the Sun, this layer is about 7 Mm thick and centered at 0.932 of a solar radius where convection cells have a typical size of about 20 Mm and a lifetime of 0.3 Ms, both similar to what is observed in supergranules. If the rms velocity of r-modes at the surface exceeds 5 m/s, then energy is deposited inside the Sun at a sufficient rate to power the supergranulation and impose on it a weak latitude dependence.

Description: We present observational evidence that eruptions of quiescent filaments and associated coronal mass ejections (CMEs) occur as a consequence of the destabilization of large-scale coronal arcades due to interactions between these structures and new and growing active regions. Both statistical and case studies have been carried out. In a case study of a 'bulge' observed by the High-Altitude Observatory Solar Maximum Mission coronagraph, the high-resolution magnetograms from the Big Bear Solar Observatory show newly emerging and rapidly changing flux in the magnetic fields that apparently underlie the bugle. For other case studies and in the statistical work the eruption of major quiescent filaments was taken as a proxy for CME eruption. We have found that two thirds of the quiescent-filament-associated CMEs occurred after substantial amounts of new magnetic flux emerged in the vicinity of the filament. In addition, in a study of all major quiescent filaments and active regions appearing in a 2-month period we found that 17 of the 22 filaments that were associated with new active regions erupted and 26 of the 31 filaments that were not associated with new flux did not erupt. In all cases in which the new flux was oriented favorably for reconnection with the preexisting large-scale coronal arcades; the filament was observed to erupt. The appearance of the new flux in the form of new active regions begins a few days before the eruption and typically is still occurring at the time of the eruption. A CME initiation scenario taking account of these observational results is proposed.

Description: Total solar irradiance measurements from the 1984-1993 Earth Radiation Budget Satellite (ERBS) active cavity radiometer and 1978-1993 Nimbus 7 transfer cavity radiometer spacecraft experiments are analyzed to detect the presence of 11-, 22-, and 80-year irradiance variability components. The analyses confirmed the existence of a significant 11-year irradiance variability component, associated with solar magnetic activity and the sunspot cycle. The analyses also suggest the presence of a 22- or 80-year variability component. The earlier Nimbus 7 and Solar Maximum Mission (SMM) spacecraft irradiance measurements decreased approximately 1.2 and 1.3 W/sq m, respectively, between 1980 and 1986. The Nimbus 7 values increased 1.2 W/sq m between 1986 and 1989. The ERBS irradiance measurements increased 1.3 W/sq m during 1986-1989, and then decreased 0.4 W/sq m (at an annual rate of 0.14 W/sq. m/yr) during 1990-1993. Considering the correlations between ERBS, Nimbus 7, and SMM irradiance trends and solar magnetic activity, the total solar irradiance should decrease to minimum levels by 1997 as solar activity decreases to minimum levels, and then increase to maximum levels by the year 2000 as solar activity rises. The ERBS measurements yielded 165.4 +/- 0.7 W/sq m as the mean irradiance value with measurement accuracies and precisions of 0.2% and 0.02%, respectively. The ERBS mean irradiance value is within 0.2% of the 1367.4, 1365.9, and 1366.9 W/sq m mean values for the SMM, Upper Atmosphere Research Satellite (UARS), and Space Shuttle Atmospheric Laboratory for Applications and Science (ATLAS 1) Solar Constant (SOLCON) active cavity radiometer spacecraft experiments, respectively. The Nimbus 7 measurements yielded 1372.1 W/sq m as the mean value with a measurement accuracy of 0.5%. Empirical irradiance model fits, based upon 10.7 -cm solar radio flux (F10) and photometric sunspot index (PSI), were used to assess the quality of the ERBS, Numbus 7, SMM, and the UARS irradiance data sets and to identify irradiance variability trends which may be caused by drifts or shifts in the spacecraft sensor responses. Comparisons among the fits and measured irradiances indicate that the Nimbus 7 radiometer response shifted by a total of 0.8 W/sq m between September 1989 and April 1990 and that the ERBS and UARS radiometers each drifted approximately 0.5 W/sq m during the first 5 months in orbit.

Description: The term 'magnetic hole' has been used to denote isolated intervals when the magnitude of the interplanetary magnetic field drops to a few tenths, or less, of its ambient value for a time that corresponds to a linear dimension of tens to a few hundreds of proton gyro-radii. Data obtained by the Ulysses magnetometer and solar wind anlayzer have been combined to study the properties of such magnetic holes in the solar wind between 1 AU and 5.4 AU and to 23 deg south latitude. In order to avoid confusion with decreases in field strength at interplanetary discontinuities, the study has focused on linear holes across which the field direction changed by less than 5 deg. The holes occurred preferentially, but not without exception, in the interaction regions on the leading edges of high-speed solar wind streams. Although the plasma surrounding the holes was generally stable against the mirror instability, there are indications that the holes may have been remnants of mirror-mode structures created upstream of the points of observation. Those indications include the following: (1) For the few holes for which proton of alpha-particle pressure could be measured inside the hole, the ion thermal pressure was always greater than in the plasma adjacent to the holes. (2) The plasma surrounding many of the holes was marginally stable for the mirror mode, while the plasma environment of all holes was significantly closer to mirror instability than was the average solar wind. (3) The plasma containing trains of closely spaced holes was closer to mirror instability than was the plasma containing isolated holes. (4) The near-hole plasma had much higher ion beta (ratio of thermal to magnetic pressure) than did the average solar wind. (5) Near the holes, T(sub perp)/T(sub parallel) tended to be either greater than 1 or larger than in the average wind. (6) The proton and alpha-particle distribution functions measured inside the holes occasionally exhibited the flattened phase-space-density contoures in nu(sub perp)/nu(sub parallel) space found in some numerical simulations of the mirror instability.

Description: Rotor noise prediction codes predict the thickness and loading noise produced by a helicopter rotor, given the blade motion, rotor operating conditions, and fluctuating force distribution over the blade surface. However, the criticality of these various inputs, and their respective effects on the predicted acoustic field, have never been fully addressed. This paper examines the importance of these inputs, and the sensitivity of the acoustic predicitions to a variation of each parameter. The effects of collective and cyclic pitch, as well as coning and cyclic flapping, are presented. Blade loading inputs are examined to determine the necessary spatial and temporal resolution, as well as the importance of the chordwise distribution. The acoustic predictions show regions in the acoustic field where significant errors occur when simplified blade motions or blade loadings are used. An assessment of the variation in the predicted acoustic field is balanced by a consideration of Central Processing Unit (CPU) time necessary for the various approximations.

Description: The U.S. National Aeronautics and Space Administration (NASA) Balloon Program has been highly successful since recovering from the catastrophic balloon failure problems of the early to mid 1980s. Balloons have continued to perform at unprecedented success rates. The comprehensive research and development (R&D) effort has continued with advances being made across the spectrum of balloon related disciplines. The long duration balloon project will be transitioning from a development effort to an operational capability this year. Recently, emphasis has been placed on the development and implementation of new support systems and facilities. A new permanent launch facility at Fort Sumner, New Mexico has been established. New ground station support equipment is being implemented, and a new heavy load launch vehicle is scheduled to be implemented in 1992. The progress, status and future plans for these and other aspects of the NASA program will be presented.

Description: The catastrophic balloon failure during the first half of the 1980's identified the need for a comprehensive and continuing balloon research and development (R&D) commitment by NASA. Technical understanding was lacking in many of the disciplines and processes associated with scientific ballooning. A comprehensive balloon R&D plan was developed in 1986 and implemented in 1987. The objectives were to develop the understanding of balloon system performance, limitations, and failure mechanisms. The program consisted of five major technical areas: structures, performance and analysis, materials, chemistry and processing, and quality control. Research activitites have been conducted at NASA/Goddard Space Flight Center (GSFC)-Wallops Flight Facility (WFF), other NASA centers and government facilities, universities, and the balloon manufacturers. Several new and increased capabilities and resources have resulted from this activity. The findings, capabilities, and plan of the balloon R&D program are presented.

Description: Caps have been used to structurally reinforce scientific research balloons since the late 1950's. The scientific research balloons used by the National Aeronautics and Space Administration (NASA) use internal caps. A NASA cap placement specification does not exist since no empirical information exisits concerning cap placement. To develop a cap placement specification, NASA has completed two in-hangar inflation tests comparing the structural contributions of internal caps and external caps. The tests used small scale test balloons designed to develop the highest possible stresses within the constraints of the hangar and balloon materials. An externally capped test balloon and an internally capped test balloon were designed, built, inflated and simulated to determine the structural contributions and benefits of each. The results of the tests and simulations are presented.

Description: During the international campaign of June 1991, the active region AR 6659 produced six very large, long-duration flares (X10/12) during its passage across the solar disk. We present the characteristics of four of them (June 4, 6, 9, 15). Precise measurements of the spot motions from Debrecen and Tokyo white-light pictures are used to understand the fragmentation of the main sunspot group with time. This fragmentation leads to a continuous restructuring of the magnetic field pattern while rapid changes are evidenced due to fast new flux emergence (magnetograms of Marshall Space Flight Center (MSFC), Huairou). The first process leads to a shearing of the field lines along which there is energy storage; the second one is the trigger which causes the release of energy by creating a complex topology. We conjecture that these two processes with different time scales are relevant to the production of flares.

Description: Interplanetary disturbances characterized by plasma that is more turbulence and/or moves faster than the background solar wind are readily defected as transients in Doppler scintillation measurements of the near-Sun solar wind. Systematic analysis of over 23,000 hours of Pioneer Venus Orbiter Doppler measurements obtained inside 0.5 AU during 1979-1987 have made it possible for the first time to investigate the frequency of occurrence of Doppler scintillation transients under solar minimum conditions and to determine its dependence on solar cycle. On the basis of a total of 142 transients, Doppler scintillation transient rates vary from a high of 0.22 in 1979 (one every 4.6 days) to a low of 0.077 transients/d in 1986 (one every 13 days), a decrease by almost a factor of 3 from solar maximum to solar minimum. This solar cycle variation, the strongest yet of any solar wind Doppler scintillation property, is highly correlated with both solar activity characterized by sunspot number and the coronal mass ejection rates deduced from Solswind and Solar Maximum Mission (SMM) coronagraph observations. These results indicate that coronal mass ejections and Doppler scintillation transients are closely related not just during solar maximum, as occasional individual comparisons have shown in the past, but throughout the entire solar cycle, and strengthen the notation that the Doppler scintillation and optical transients are different manifestations of the same physical phenomenon. The magnitudes of the transients, as described by the ratio of peak to pretransient scintillation levels (EF for enhancement factor), and their distribution iwth heliocentric distance also vary with solar cycle. While EF tends to diminish with increasing heliocentric distance during high solar activity, it is more evenly distributed during low solar activity. EF is also lower during solar minimum, as 13% of the transients during solar maximum have values exceeding 23, the highest EF observed during solar minimum. These results are consistent with the fact that occasional major fast-moving interplanetary shocks that are observed during solar maximum are very rate during solar minimum.

Description: It has recently been reported that the total radiative emission variations from solar type stars exceeds the currently solar constant variations (from spacecraft over the last decade) by a factor near 4. Aside from other remote alternatives, this suggests three clear possibilities: (1) the Sun may undergo irradiance variations several times larger than any we have seen; (2) our Sun is highly unusual with regard to its radiative output; or (3) our terrestrial position in the heliosphere provides a special vantage point which reduces the observed solar irradiance variations. We investigate the last possibility by considering the influence of observer latitude upon calculated irradiance variations using a simple model for emission from solar contrast features. We consider modeled sunspots, faculae, and network structures. As the latitude angle of the observer rises relative to the heliographic equator, sunspot deficit contributions diminish and facular plus network contributions escalate. We find that the observing latitude can influence the irradiance variations by a factor near 6. When we integrate the irradiance variations, over the celestial sphere, they average to 3 times the terrestrial effect, suggesting that the solar cycle luminosity variations are proportionally, 3 times larger than the solar constant variations. Thus we suggest the Sun's luminosity output varies even more strongly with the solar cycle than is apparent in the solar constant variations. The influence of the observer viewing angle relative to stellar spin axis, studied here, may be possible to investigate with a thorough statistical examination of other type stars. Additionally, the rotational modulation due to active regions (as a function of observer viewing angle) may also be a valuable are for future investigation.

Description: For more than 90 years, solar extreme ultraviolet (EUV) irradiance modeling has progressed from empirical blackbody radiation formulations, through fudge factors, to typically measured irradiances and reference spectra was well as time-dependent empirical models representing continua and line emissions. A summary of recent EUV measurements by five rockets and three satellites during the 1980s is presented along with the major modeling efforts. The most significant reference spectra are reviewed and threee independently derived empirical models are described. These include Hinteregger's 1981 SERF1, Nusinov's 1984 two-component, and Tobiska's 1990/1991/SERF2/EUV91 flux models. They each provide daily full-disk broad spectrum flux values from 2 to 105 nm at 1 AU. All the models depend to one degree or another on the long time series of the Atmosphere Explorer E (AE-E) EUV database. Each model uses ground- and/or space-based proxies to create emissions from solar atmospheric regions. Future challenges in EUV modeling are summarized including the basic requirements of models, the task of incorporating new observations and theory into the models, the task of comparing models with solar-terrestrial data sets, and long-term goals and modeling objectives. By the late 1990s, empirical models will potentially be improved through the use of proposed solar EUV irradiance measurements and images at selected wavelengths that will greatly enhance modeling and predictive capabilities.

Description: The purpose of this Note is to present results from an analytic/experimental study that investigated the potential for passively changing blade twist through the use of extension-twist coupling. A set of composite model rotor blades was manufactured from existing blade molds for a low-twist metal helicopter rotor blade, with a view toward establishing a preliminary proof concept for extension-twist-coupled rotor blades. Data were obtained in hover for both a ballasted and unballasted blade configuration in sea-level atmospheric conditions. Test data were compared with results obtained from a geometrically nonlinear analysis of a detailed finite element model of the rotor blade developed in MSC/NASTRAN.

Description: The paper considers the compressible Rayleigh equation as a model for the Mach wave emission mechanism associated with high-temperature supersonic jets. Solutions to the compressible Rayleigh equation reveal the existence of several families of supersonically convecting instability waves. These waves directly radiate noise to the jet far field. The predicted noise characteristics are compared to previously acquired experimental data for an axisymmetric Mach 2 fully pressure balanced jet operating over a range of jet total temperatures from ambient to 1370 K. The results of this comparison show that the first-order supersonic instability wave and the Kelvin-Hemlhlotz first-, second-, and third-order modes have directional radiation characteristics that are in agreement with observed data. The assumption of equal initial amplitudes for all of the waves leads to the conclusion that the flapping mode of instability dominates the noise radiatio process of supersonic jets. At a jet temperature of 1370 K, supersonic instability waves are predicted to dominate the noise radiated at high frequency at narrow angles to the jet axis.

Description: The objective of the present work is to study the mixing characteristics of a linear array of supersonic rectangular jets under conditions of screech synchronization. The screech synchronization at a fully expanded jet Mach number of 1.61 is achieved by a precise adjustment of the internozzle spacing. To our knowledge, such an experiment on the resonant mixing of screech synchronized multiple rectangular jets has not been reported before. The results are compared with the case where the screech was suppressed in the multijet configuration.

Description: The objective of the present investigation is to assess the effect of the spatial order of accuracy used for the evaluation of the inviscid fluxes on the resolution of higher order quantitites, such as velocity gradients. The viscous terms are computed as second-order accurate with central difference formulas, even though for the explicit part of the algorithm higher order approximations may be used. A viscous/inviscid method is used, and the outer part of the flowfield is computed with the inviscid flow equations. The viscous boundary-layer type flow region close to the body surface is computed with an algebraic eddy viscosity model. Results obtained with the conservative and nonconservative formulations and the viscous/inviscid approach are compared with available experimental data. The effect of grid refinement on the accuracy of the solution is also presented.

Description: Enhanced fluxes of suprathermal electrons are commonly observed upstream of corotating forward and reverse shocks in the solar wind at heliocentric distances beyond approximately 2 AU by the Los Alamos plasma experiment on Ulysses. The average duration of these events, which are most intense immediately upstream from the shocks and which fade with increasing distance from them, is approximately 2.4 days near 5 AU. These events are caused by the leakage of shock-heated electrons into the upstream region. The upstream regions of these shocks face back toward the Sun along the interplanetary magnetic field, so these leaked electrons commonly counterstream relative to the normal solar wind electron heat flux. The observations suggest that conservation of magnetic moment and scattering typically limit the sunward propagation of these electrons as beams to field-aligned distances of approximately 15 AU. Although it seems unlikely that these shock-associated events are an important source of counterstreaming events near 1 AU, remnants of the backstreaming beams may contribute importantly to the diffuse solar wind halo electron population there.

Description: In May, 1993, the heliospheric current sheet (HCS) ceased to be seen by the Ulysses spacecraft at a heliocentric latitude of approximately 30 deg S and distance of 4.7 AU. The disappearance of the HCS coincided with the solar wind speed remaining greater than 560 km/s and with the disappearance of one of four interaction regions previously seen on each solar rotation. The heliographic latitude of the disappearance of the HCS at Ulysses was 11 deg equatorward of the latitude of the magnetic neutral sheet computed at the source surface at 2.5 solar radii, and it occurred a half year earlier than predicted on the basis of the persistance of the time profile of the neutral sheet tilt from one solar cycle to the next.

Description: The unique vantage point of the Ulysses spacecraft throughout 1992 and the beginning of 1993, at a close to constant heliocentric distance of about 5 AU and a slowly varying heliographic latitude from 5 deg to 30 deg south is used to describe and discuss the evolution of the sector structure of the interplanetary magnetic field during the declining phase of the solar cycle. From the end of 1990 to the beginning of 1992 the sector structure changed from a four sector to a two sector structure, but remained constant in solar longitude. From about June-July 1992, the structure, matching the evolution in the computed coronal magnetic fields, drifted eastwards, with a recurrence period of about 28 days. This result may indicate a slower rotation rate for the dipolar component of the solar magnetic field which becomes dominant about this time in the solar cycle.

Description: Near-ecliptic solar wind observations by Ulysses on its way to the polar regions of the Sun, compared with those from IMP 8 at 1 AU, showed that high-speed streams decay and broaden with heliocentric distance from IMP 8 to Ulysses, as expected. In July 1992 while traveling south at approximately 13 deg S and 5.3 AU, Ulysses encountered a recurrent high-speed stream, that may also have been observed at IMP 8. The stream has been observed a total of 14 times, once in each solar rotation through June 1993 at approximately 34 deg S. The source of the high-speed stream is an equatorward extension of the south polar coronal hole. From July 1992 through June 1993, averages of solar wind peak speed increased while density decreased with heliographic latitude. Both the stream and a low-speed, high-density flow, presumably associated with the heliomagnetic (coronal) streamer belt encircling the heliomagnetic equator, crossed Ulysses with the solar rotation period until April 1993 when the spacecraft was at approximately 29 deg S heliographic latitude. After this time, as the spacecraft climbed to higher latitudes, the central portion of the streamer belt with lowest speed and highest density disappeared. Therefore, at its maximum inclination, the belt was tilted at approximately 29 deg to the heliographic equator at this point in the solar cycle.

Description: Compressible MHD simulations in one dimension with three-dimensional vectors are used to investigate a number of processes relevant to problems in interplanetary physics. The simulations indicate that a large-amplitude nonequilibrium (e.g., linearly polarized) Alfvenic wave, which always starts with small relative fluctuations in the magnitude B of the magnetic field, typically evolves to flatten the magnetic profile in most regions. Under a wide variety of conditions B and the density rho become anticorrelated on average. If the mean magnetic field is allowed to decrease in time, the point where the transverse magnetic fluctuation amplitude delta B(sub T) is greater than the mean field B(sub 0) is not special, and large values of delta B(sub T)/B(sub 0) do not cause the compressive thermal energy to increase remarkably or the wave energy to dissipate at an unusually high rate. Nor does the 'backscatter' of the waves that occurs when the sound speed is less than the Alfven speed result, in itself, in substantial energy dissipation, but rather primarily in a phase change between the magnetic and velocity fields. For isolated wave packets the backscatter does not occur for any of the parameters examined; an initial radiation of acoustic waves away from the packet establishes a stable traveling structure. Thus these simulations, although greatly idealized compared to reality, suggest a picture in which the interplanetary fluctuations should have small deltaB and increasingly quasi-pressure balanced compressive fluctuations, as observed, and in which the dissipation and 'saturation' at delta B(sub T)/B(sub 0) approximately = 1 required by some theories of wave acceleration of the solar wind do not occur. The simulations also provide simple ways to understand the processes of nonlinear steepening and backscattering of Alfven waves and demonstrate the existence of previously unreported types of quasi-steady MHD states.

Description: The benefits of using a hypersonic waverider for spacecraft trajectory modification are presented. A waverider is a hypersonic vehicle specifically designed so that the undersurface bow shock is attached to the leading edge, which provides for the highest known lift-to-drag ratios achievable at high Mach number flight. Several viable space missions are suggested which could use such configurations for low-drag aero-assisted maneuvers in planetary atmospheres. It is shown that large changes in the spacecraft velocity vector can be accomplished with acceptably small losses in energy due to drag using a waverider aeroshell. The primary advantage of an aero-assist maneuver is suggested by comparison to a traditional gravity-assist trajectory. Some scaling laws are presented for comparing waveriders designed for different planetary atmospheres, and it is shown that the compositional differences between the terrestrial planets has a minimal impact on waverider design.

Description: During the late declining phase of the solar cycle, the tilt of the solar magnetic dipole with respect to the Sun's rotation axis leads to large-scale organization of the solar wind, such that alternating regions of high- and low-speed solar wind are observed in the ecliptic plane. In this paper, we use Doppler scintillation measurements to investigate mass flux of these two types of solar wind in the ecliptic plane and inside 0.3 AU, where in situ measurements have not been possible. To the extent that Doppler scintillation reflects mass flux, we find that mass flux in high-speed streams: (1) is lower (by a factor of approximately 2.2) than the mass flux of the average solar wind in the heliocentric distance range of 0.3-0.5 AU; (2) is lower still (by as much as a factor of about 4) than the mass flux of the slow solar wind associated with the streamer belt; and (3) appears to grow with heliocentric distance. These Doppler scintillation results are consistent with the equator to pole decrease in mass flux observed in earlier spectral broadening measurements, and with trends and differences between high- and low-speed solar wind observed by in situ measurements in the range of 0.3-0.1 AU. The mass flux results suggest that the solar wind flow in high-speed streams is convergent towards the ecliptic near the Sun, becoming less convergent and approaching radial with increasing heliocentric distance beyond 0.3 AU. The variability of mass flux observed within equatorial and polar high-speed streams close to the Sun is strikingly low. This low variability implies that, as Ulysses currently ascends to higher latitudes and spends more time in the south polar high-speed stream after crossing the heliocentric current sheet, it can expect to observe a marked decrease in variations of both mass flux and solar wind speed, a trend that appears to have started already.

Description: Highly structured Langmuir waves, also known as electron plasma oscillations, have been observed in the foreshock of Venus using the plasma wave experiment on the Galileo spacecraft during the gravity assist flyby on February 10, 1990. The Galileo wideband sampling system provides digital electric field waveform measurements at sampling rates up to 201,600 samples per second, much higher than any previous instrument of this type. The main Langmuir wave emission band occurs near the local electron plasma frequency, which was approximately 43 kHz. The Langmuir waves are observed to shift above and below the plasma frequency, sometimes by as much as 20 kHz. The shifts in frequency are closely correlated with the downstream distance from the tangent field line, implying that the shifts are controlled by the electron beam velocity. Considerable fine structure is also evident, with timescales as short as 0.15 ms, corresponding to spatial scales of a few tens of Debye lengths. The frequency spectrum often consists of beat-type waveforms, with beat frequencies ranging from 0.2 to 7 kHz, and in a few cases, isolated wave packets. The peak electric field strengths are approximately 1 mV/m. These field strengths are too small for strongly nonlinear processes to be important. The beat-type waveforms are suggestive of a parametric decay process.

Description: The irregular polygonal pattern of solar granulation is analyzed for size-shape relations using statistical crystallography. In contrast to previous work which has assumed perfectly hexagonal patterns for granulation, more realistic accounting of cell (granule) shapes reveals a broader basis for quantitative analysis. Several features emerge as noteworthy: (1) a linear correlation between number of cell-sides and neighboring shapes (called Aboav-Weaire's law); (2) a linear correlation between both average cell area and perimeter and the number of cell-sides (called Lewis's law and a perimeter law, respectively) and (3) a linear correlation between cell area and squared perimeter (called convolution index). This statistical picture of granulation is consistent with a finding of no correlation in cell shapes beyond nearest neighbors. A comparative calculation between existing model predictions taken from luminosity data and the present analysis shows substantial agreements for cell-size distributions. A model for understanding grain lifetimes is proposed which links convective times to cell shape using crystallographic results.

Description: Dual-frequency ranging and Doppler measurements were conducted in support of the Ulysses Solar Corona Experiment (SCE) at and around the spacecraft's first solar conjunction in 1991 August. The differential group delay time between range codes on the two downlink carrier signals at the wavelengths 13.1 and 3.6 cm, a direct measure of the total electron content between spacecraft and ground station, was used to derive the electron density distribution in the solar corona. Linear power-law representations of the coronal electron density were derived for the range of solar distances from 4 solar radii to 40 solar radii on both sides of the Sun. The corona was found to be very nearly symmetric; the radial falloff exponent being 2.54 +/- 0.05 for occultation ingress (east solar limb) and 2.42 +/- 0.05 for egress (west limb), respectively. The departure of these exponents from the inverse equare relation implies that significant solar wind acceleration is occurring within the radial range of the observations. The electron density level was found to be considerably lower than that observed during the 1988 December solar occultation of Voyager 2. Although the smoothed sunspot number R(sub z) (a standard indicator of solar activity) was almost the same in 1988 December and 1991 August, the mean electron density at 20 solar radii was found to be 1.7 +/- 0.1 x 10(exp 3)/cu cm during the Ulysses conjunction, a decline by almost a factor of 4 from the value obtained during the Voyager conjunction.

Description: One of the long-standing uncertainties in the wave-resonance theory of coronal heating is the stability of the resonance layer. The wave motions in the resonance layer produce highly localized shear flows which vary sinusoidally in time with the resonance period. This configuration is potentially susceptible to the Kelvin-Helmholtz instability (KHI), which can enhance small-scale structure and turbulent broadening of shear layers on relatively rapid ideal timescales. We have investigated numerically the response of a characteristic velocity profile, derived from resonance absorption models, to finite fluid perturbations comparable to photospheric fluctuations. We find that the KHI primarily should affect long (approximately greater than 6 x 10(exp 4) km) loops where higher velocity flows (M approximately greater than 0.2) exist in resonance layers of order 100 km wide. There, the Kelvin-Helmholtz growth time is comparable to or less than the resonance quarter-period, and the potentially stabilizing magnetic effects are not felt until the instability is well past the linear growth stage. Not only is the resonance layer broadened by the KHI, but also the convective energy transport out of the resonance layer is increased, thus adding to the efficiency of the wave-resonance heating process. In shorter loops, e.g., those in bright points and compact flares, the stabilization due to the magnetic field and the high resonance frequency inhibit the growth of the Kelvin-Helmholtz instability beyond a minimal level.

Description: The solar flare and coronal mass ejection (CME) events associated with the large and complex March 1989 active region are discussed. This active region gave us a chance to study the relation of CME with truly major solar flares. The work concentrates on questions of the relation of CMEs and flares to one another and to other types of activity on the Sun. As expected, some major (X-3B class) flares had associated CMEs. However, an unexpected finding is that others did not. In fact, there is strong evidence that the X4-4B flare of March 9th had no CME. This lack of a CME for such an outstanding flare event has important implications to theories of CME causation.Apparently, not all major flares cause CMEs or are caused by CMEs. The relations between CMEs and other types of solar activity are also discussed. No filament disappearances are reported for major CMEs studied here. Comparing these results with other studies, CMEs occur in association with flares and with erupting prominences, but neither are required for a CME. The relation between solar structures showing flaring without filament eruptions and structures showing filament eruptions without flares becomes important. The evolutionary relation between an active flaring sunspot region and extensive filaments without sunspots is reviewed, and the concept of an 'evolving magnetic structure' (EMS) is introduced. It is suggested that all CMEs arise in EMSs and that CMEs provide a major path through which azimuthal magnetic fields escape form the Sun during the solar cycle.

Description: The large gradual solar-energetic-particle (SEP) events, where abundances are commonly measured, are produced when coronal mass ejections (CMEs) drive shock waves through the corona and the interplanetary medium. The shock accelerates particles from the highly-ionized, approximately 1.5 MK, plasma in a manner that depends only weakly upon the Q/A of the ion, except at very high energies. Averaging the approximately 1 MeV/amu abundances over many events compensates for the acceleration effects to produce abundances that appear to correspond directly to those in the coronal source for all observed elements, including H. The resulting abundances reflect the 4 x enhancement of ions with low values of first ionization potential (FIP) arising from ion-neutral fractionation that occurs as the atoms are transported up from the photosphere. A different pattern of fractionation is found for ions that are shock-accelerated from the high speed solar wind emerging from coronal holes.

Description: ARC2D is a computational fluid dynamics program developed at the NASA Ames Research Center specifically for airfoil computations. The program uses implicit finite-difference techniques to solve two-dimensional Euler equations and thin layer Navier-Stokes equations. It is based on the Beam and Warming implicit approximate factorization algorithm in generalized coordinates. The methods are either time accurate or accelerated non-time accurate steady state schemes. The evolution of the solution through time is physically realistic; good solution accuracy is dependent on mesh spacing and boundary conditions. The mathematical development of ARC2D begins with the strong conservation law form of the two-dimensional Navier-Stokes equations in Cartesian coordinates, which admits shock capturing. The Navier-Stokes equations can be transformed from Cartesian coordinates to generalized curvilinear coordinates in a manner that permits one computational code to serve a wide variety of physical geometries and grid systems. ARC2D includes an algebraic mixing length model to approximate the effect of turbulence. In cases of high Reynolds number viscous flows, thin layer approximation can be applied. ARC2D allows for a variety of solutions to stability boundaries, such as those encountered in flows with shocks. The user has considerable flexibility in assigning geometry and developing grid patterns, as well as in assigning boundary conditions. However, the ARC2D model is most appropriate for attached and mildly separated boundary layers; no attempt is made to model wake regions and widely separated flows. The techniques have been successfully used for a variety of inviscid and viscous flowfield calculations. The Cray version of ARC2D is written in FORTRAN 77 for use on Cray series computers and requires approximately 5Mb memory. The program is fully vectorized. The tape includes variations for the COS and UNICOS operating systems. Also included is a sample routine for CONVEX computers to emulate Cray system time calls, which should be easy to modify for other machines as well. The standard distribution media for this version is a 9-track 1600 BPI ASCII Card Image format magnetic tape. The Cray version was developed in 1987. The IBM ES/3090 version is an IBM port of the Cray version. It is written in IBM VS FORTRAN and has the capability of executing in both vector and parallel modes on the MVS/XA operating system and in vector mode on the VM/XA operating system. Various options of the IBM VS FORTRAN compiler provide new features for the ES/3090 version, including 64-bit arithmetic and up to 2 GB of virtual addressability. The IBM ES/3090 version is available only as a 9-track, 1600 BPI IBM IEBCOPY format magnetic tape. The IBM ES/3090 version was developed in 1989. The DEC RISC ULTRIX version is a DEC port of the Cray version. It is written in FORTRAN 77 for RISC-based Digital Equipment platforms. The memory requirement is approximately 7Mb of main memory. It is available in UNIX tar format on TK50 tape cartridge. The port to DEC RISC ULTRIX was done in 1990. COS and UNICOS are trademarks and Cray is a registered trademark of Cray Research, Inc. IBM, ES/3090, VS FORTRAN, MVS/XA, and VM/XA are registered trademarks of International Business Machines. DEC and ULTRIX are registered trademarks of Digital Equipment Corporation.

Description: Panel method computer programs are software tools of moderate cost used for solving a wide range of engineering problems. The panel code PMARC_12 (Panel Method Ames Research Center, version 12) can compute the potential flow field around complex three-dimensional bodies such as complete aircraft models. PMARC_12 is a well-documented, highly structured code with an open architecture that facilitates modifications and the addition of new features. Adjustable arrays are used throughout the code, with dimensioning controlled by a set of parameter statements contained in an include file; thus, the size of the code (i.e. the number of panels that it can handle) can be changed very quickly. This allows the user to tailor PMARC_12 to specific problems and computer hardware constraints. In addition, PMARC_12 can be configured (through one of the parameter statements in the include file) so that the code's iterative matrix solver is run entirely in RAM, rather than reading a large matrix from disk at each iteration. This significantly increases the execution speed of the code, but it requires a large amount of RAM memory. PMARC_12 contains several advanced features, including internal flow modeling, a time-stepping wake model for simulating either steady or unsteady (including oscillatory) motions, a Trefftz plane induced drag computation, off-body and on-body streamline computations, and computation of boundary layer parameters using a two-dimensional integral boundary layer method along surface streamlines. In a panel method, the surface of the body over which the flow field is to be computed is represented by a set of panels. Singularities are distributed on the panels to perturb the flow field around the body surfaces. PMARC_12 uses constant strength source and doublet distributions over each panel, thus making it a low order panel method. Higher order panel methods allow the singularity strength to vary linearly or quadratically across each panel. Experience has shown that low order panel methods can provide nearly the same accuracy as higher order methods over a wide range of cases with significantly reduced computation times; hence, the low order formulation was adopted for PMARC_12. The flow problem is solved by modeling the body as a closed surface dividing space into two regions: the region external to the surface in which an unknown velocity potential exists representing the flow field of interest, and the region internal to the surface in which a known velocity potential (representing a fictitious flow) is prescribed as a boundary condition. Both velocity potentials are required to satisfy Laplace's equation. A surface integral equation for the unknown potential external to the surface can be written by applying Green's Theorem to the external region. Using the internal potential and zero flow through the surface as boundary conditions, the unknown potential external to the surface can be solved for. When the internal flow option, which allows the analysis of closed ducts, wind tunnels, and similar internal flow problems, is selected, the geometry is modeled such that the flow field of interest is inside the geometry and the fictitious flow is outside the geometry. Items such as wings, struts, or aircraft models can be included in the internal flow problem. The time-stepping wake model gives PMARC_12 the ability to model both steady and unsteady flow problems. The wake is convected downstream from the wake-separation line by the local velocity field. With each time step, a new row of wake panels is added to the wake at the wake-separation line. Time stepping can start from time t=0 (no initial wake) or from time t=t0 (an initial wake is specified). A wide range of motions can be prescribed, including constant rates of translation, constant rate of rotation about an arbitrary axis, oscillatory translation, and oscillatory rotation about any of the three coordinate axes. Investigators interested in a visual representation of the phenomenon they are studying with PMARC_12 may want to consider obtaining the program GVS (ARC-13361), the General Visualization System. GVS is a Silicon Graphics IRIS program which was created for the purpose of supporting the scientific visualization needs of PMARC_12. GVS is available separately from COSMIC. PMARC_12 is written in standard FORTRAN 77, with the exception of the NAMELIST extension used for input. This makes the code fairly machine independent. A compiler which supports the NAMELIST extension is required. The amount of free disk space and RAM memory required for PMARC_12 will vary depending on how the code is dimensioned using the parameter statements in the include file. The recommended minimum requirements are 20Mb of free disk space and 4Mb of RAM. PMARC_12 has been successfully implemented on a Macintosh II running System 6.0.7 or 7.0 (using MPW/Language Systems Fortran 3.0), a Sun SLC running SunOS 4.1.1, an HP 720 running HP-UX 8.07, an SGI IRIS running IRIX 4.0 (it will not run under IRIX 3.x.x without modifications), an IBM RS/6000 running AIX, a DECstation 3100 running ULTRIX, and a CRAY-YMP running UNICOS 6.0 or later. Due to its memory requirements, this program does not readily lend itself to implementation on MS-DOS based machines. The standard distribution medium for PMARC_12 is a set of three 3.5 inch 800K Macintosh format diskettes and one 3.5 inch 1.44Mb Macintosh format diskette which contains an electronic copy of the documentation in MS Word 5.0 format for the Macintosh. Alternate distribution media and formats are available upon request, but these will not include the electronic version of the document. No executables are included on the distribution media. This program is an update to PMARC version 11, which was released in 1989. PMARC_12 was released in 1993. It is available only for use by United States citizens.

Description: This program determines the supersonic flowfield surrounding three-dimensional wing-body configurations of a delta wing. It was designed to provide the numerical computation of three dimensional inviscid, flowfields of either perfect or real gases about supersonic or hypersonic airplanes. The governing equations in conservation law form are solved by a finite difference method using a second order noncentered algorithm between the body and the outermost shock wave, which is treated as a sharp discontinuity. Secondary shocks which form between these boundaries are captured automatically. The flowfield between the body and outermost shock is treated in a shock capturing fashion and therefore allows for the correct formation of secondary internal shocks . The program operates in batch mode, is in CDC update format, has been implemented on the CDC 7600, and requires more than 140K (octal) word locations.

Description: This theoretical aerodynamics program, TAD, was developed to predict the aerodynamic characteristics of vehicles with sounding rocket configurations. These slender, axisymmetric finned vehicle configurations have a wide range of aeronautical applications from rockets to high speed armament. Over a given range of Mach numbers, TAD will compute the normal force coefficient derivative, the center-of-pressure, the roll forcing moment coefficient derivative, the roll damping moment coefficient derivative, and the pitch damping moment coefficient derivative of a sounding rocket configured vehicle. The vehicle may consist of a sharp pointed nose of cone or tangent ogive shape, up to nine other body divisions of conical shoulder, conical boattail, or circular cylinder shape, and fins of trapezoid planform shape with constant cross section and either three or four fins per fin set. The characteristics computed by TAD have been shown to be accurate to within ten percent of experimental data in the supersonic region. The TAD program calculates the characteristics of separate portions of the vehicle, calculates the interference between separate portions of the vehicle, and then combines the results to form a total vehicle solution. Also, TAD can be used to calculate the characteristics of the body or fins separately as an aid in the design process. Input to the TAD program consists of simple descriptions of the body and fin geometries and the Mach range of interest. Output includes the aerodynamic characteristics of the total vehicle, or user-selected portions, at specified points over the mach range. The TAD program is written in FORTRAN IV for batch execution and has been implemented on an IBM 360 computer with a central memory requirement of approximately 123K of 8 bit bytes. The TAD program was originally developed in 1967 and last updated in 1972.

Description: This program, which is called 'AOFA', determines the complete viscous and inviscid flow around a body of revolution at a given angle of attack and traveling at supersonic speeds. The viscous calculations from this program agree with experimental values for surface and pitot pressures and with surface heating rates. At high speeds, lee-side flows are important because the local heating is difficult to correlate and because the shed vortices can interact with vehicle components such as a canopy or a vertical tail. This program should find application in the design analysis of any high speed vehicle. Lee-side flows are difficult to calculate because thin-boundary-layer theory is not applicable and the concept of matching inviscid and viscous flow is questionable. This program uses the parabolic approximation to the compressible Navier-Stokes equations and solves for the complete inviscid and viscous regions of flow, including the pressure. The parabolic approximation results from the assumption that the stress derivatives in the streamwise direction are small in comparison with derivatives in the normal and circumferential directions. This assumption permits the equation to be solved by an implicit finite difference marching technique which proceeds downstream from the initial data point, provided the inviscid portion of flow is supersonic. The viscous cross-flow separation is also determined as part of the solution. To use this method it is necessary to first determine an initial data point in a region where the inviscid portion of the flow is supersonic. Input to this program consists of two parts. Problem description is conveyed to the program by namelist input. Initial data is acquired by the program as formatted data. Because of the large amount of run time this program can consume the program includes a restart capability. Output is in printed format and magnetic tape for further processing. This program is written in FORTRAN IV and has been implemented on a CDC 7600 with a central memory requirement of approximately 35K (octal) of 60 bit words.

Description: The Comprehensive Analytical Model of Rotorcraft Aerodynamics, CAMRAD, program is designed to calculate rotor performance, loads, and noise; helicopter vibration and gust response; flight dynamics and handling qualities; and system aeroelastic stability. The analysis is a consistent combination of structural, inertial, and aerodynamic models applicable to a wide range of problems and a wide class of vehicles. The CAMRAD analysis can be applied to articulated, hingeless, gimballed, and teetering rotors with an arbitrary number of blades. The rotor degrees of freedom included are blade/flap bending, rigid pitch and elastic torsion, and optionally gimbal or teeter motion. General two-rotor aircrafts can be modeled. Single main-rotor and tandem helicopter and sideby-side tilting proprotor aircraft configurations can be considered. The case of a rotor or helicopter in a wind tunnel can also be modeled. The aircraft degrees of freedom included are the six rigid body motion, elastic airframe motions, and the rotor/engine speed perturbations. CAMRAD calculates the load and motion of helicopters and airframes in two stages. First the trim solution is obtained; then the flutter, flight dynamics, and/or transient behavior can be calculated. The trim operating conditions considered include level flight, steady climb or descent, and steady turns. The analysis of the rotor includes nonlinear inertial and aerodynamic models, applicable to large blade angles and a high inflow ratio, The rotor aerodynamic model is based on two-dimensional steady airfoil characteristics with corrections for three-dimensional and unsteady flow effects, including a dynamic stall model. In the flutter analysis, the matrices are constructed that describe the linear differential equations of motion, and the equations are analyzed. In the flight dynamics analysis, the stability derivatives are calculated and the matrices are constructed that describe the linear differential equations of motion. These equations are analyzed. In the transient analysis, the rigid body equations of motion are numerically integrated, for a prescribed transient gust or control input. The CAMRAD program product is available by license for a period of ten years to domestic U.S. licensees. The licensed program product includes the CAMRAD source code, command procedures, sample applications, and one set of supporting documentation. Copies of the documentation may be purchased separately at the price indicated below. CAMRAD is written in FORTRAN 77 for the DEC VAX under VMS 4.6 with a recommended core memory of 4.04 megabytes. The DISSPLA package is necessary for graphical output. CAMRAD was developed in 1980.

Description: ARC2D is a computational fluid dynamics program developed at the NASA Ames Research Center specifically for airfoil computations. The program uses implicit finite-difference techniques to solve two-dimensional Euler equations and thin layer Navier-Stokes equations. It is based on the Beam and Warming implicit approximate factorization algorithm in generalized coordinates. The methods are either time accurate or accelerated non-time accurate steady state schemes. The evolution of the solution through time is physically realistic; good solution accuracy is dependent on mesh spacing and boundary conditions. The mathematical development of ARC2D begins with the strong conservation law form of the two-dimensional Navier-Stokes equations in Cartesian coordinates, which admits shock capturing. The Navier-Stokes equations can be transformed from Cartesian coordinates to generalized curvilinear coordinates in a manner that permits one computational code to serve a wide variety of physical geometries and grid systems. ARC2D includes an algebraic mixing length model to approximate the effect of turbulence. In cases of high Reynolds number viscous flows, thin layer approximation can be applied. ARC2D allows for a variety of solutions to stability boundaries, such as those encountered in flows with shocks. The user has considerable flexibility in assigning geometry and developing grid patterns, as well as in assigning boundary conditions. However, the ARC2D model is most appropriate for attached and mildly separated boundary layers; no attempt is made to model wake regions and widely separated flows. The techniques have been successfully used for a variety of inviscid and viscous flowfield calculations. The Cray version of ARC2D is written in FORTRAN 77 for use on Cray series computers and requires approximately 5Mb memory. The program is fully vectorized. The tape includes variations for the COS and UNICOS operating systems. Also included is a sample routine for CONVEX computers to emulate Cray system time calls, which should be easy to modify for other machines as well. The standard distribution media for this version is a 9-track 1600 BPI ASCII Card Image format magnetic tape. The Cray version was developed in 1987. The IBM ES/3090 version is an IBM port of the Cray version. It is written in IBM VS FORTRAN and has the capability of executing in both vector and parallel modes on the MVS/XA operating system and in vector mode on the VM/XA operating system. Various options of the IBM VS FORTRAN compiler provide new features for the ES/3090 version, including 64-bit arithmetic and up to 2 GB of virtual addressability. The IBM ES/3090 version is available only as a 9-track, 1600 BPI IBM IEBCOPY format magnetic tape. The IBM ES/3090 version was developed in 1989. The DEC RISC ULTRIX version is a DEC port of the Cray version. It is written in FORTRAN 77 for RISC-based Digital Equipment platforms. The memory requirement is approximately 7Mb of main memory. It is available in UNIX tar format on TK50 tape cartridge. The port to DEC RISC ULTRIX was done in 1990. COS and UNICOS are trademarks and Cray is a registered trademark of Cray Research, Inc. IBM, ES/3090, VS FORTRAN, MVS/XA, and VM/XA are registered trademarks of International Business Machines. DEC and ULTRIX are registered trademarks of Digital Equipment Corporation.

Description: Panel methods are moderate cost tools for solving a wide range of engineering problems. PMARC (Panel Method Ames Research Center) is a potential flow panel code that numerically predicts flow fields around complex three-dimensional geometries. PMARC's predecessor was a panel code named VSAERO which was developed for NASA by Analytical Methods, Inc. PMARC is a new program with many additional subroutines and a well-documented code suitable for powered-lift aerodynamic predictions. The program's open architecture facilitates modifications or additions of new features. Another improvement is the adjustable size code which allows for an optimum match between the computer hardware available to the user and the size of the problem being solved. PMARC can be resized (the maximum number of panels can be changed) in a matter of minutes. Several other state-of-the-art PMARC features include internal flow modeling for ducts and wind tunnel test sections, simple jet plume modeling essential for the analysis and design of powered-lift aircraft, and a time-stepping wake model which allows the study of both steady and unsteady motions. PMARC is a low-order panel method, which means the singularities are distributed with constant strength over each panel. In many cases low-order methods can provide nearly the same accuracy as higher order methods (where the singularities are allowed to vary linearly or quadratically over each panel). Low-order methods have the advantage of a shorter computation time and do not require exact matching between panels. The flow problem is solved by assuming that the body is at rest in a moving flow field. The body is modeled as a closed surface which divides space into two regions -- one region contains the flow field of interest and the other contains a fictitious flow. External flow problems, such as a wing in a uniform stream, have the external region as the flow field of interest and the internal flow as the fictitious flow. This arrangement is reversed for internal flow problems where the internal region contains the flow field of interest and the external flow field is fictitious. In either case it is assumed that the velocity potentials in both regions satisfy Laplace's equation. PMARC has extensive geometry modeling capabilities for handling complex, three-dimensional surfaces. As with all panel methods, the geometry must be modeled by a set of panels. For convenience, the geometry is usually subdivided into several pieces and modeled with sets of panels called patches. A patch may be folded over on itself so that opposing sides of the patch form a common line. For example, wings are normally modeled with a folded patch to form the trailing edge of the wing. PMARC also has the capability to automatically generate a closing tip patch. In the case of a wing, a tip patch could be generated to close off the wing's third side. PMARC has a simple jet model for simulating a jet plume in a crossflow. The jet plume shape, trajectory, and entrainment velocities are computed using the Adler/Baron jet in crossflow code. This information is then passed back to PMARC. The wake model in PMARC is a time-stepping wake model. The wake is convected downstream from the wake separation line by the local velocity flowfield. With each time step, a new row of wake panels is added to the wake at the wake separation line. PMARC also allows an initial wake to be specified if desired, or, as a third option, no wakes need be modeled. The effective presentation of results for aerodynamics problems requires the generation of report-quality graphics. PMAPP (ARC-12751), the Panel Method Aerodynamic Plotting Program, (Sterling Software), was written for scientists at NASA's Ames Research Center to plot the aerodynamic analysis results (flow data) from PMARC. PMAPP is an interactive, color-capable graphics program for the DEC VAX or MicroVAX running VMS. It was designed to work with a variety of terminal types and hardcopy devices. PMAPP is available separately from COSMIC. PMARC was written in standard FORTRAN77 using adjustable size arrays throughout the code. Redimensioning PMARC will change the amount of disk space and memory the code requires to be able to run; however, due to its memory requirements, this program does not readily lend itself to implementation on MS-DOS based machines. The program was implemented on an Apple Macintosh (using 2.5 MB of memory) and tested on a VAX/VMS computer. The program is available on a 3.5 inch Macintosh format diskette (standard media) or in VAX BACKUP format on TK50 tape cartridge or 9-track magnetic tape. PMARC was developed in 1989.

Description: The Transonic Airfoil analysis computer code, TAIR, was developed to employ a fast, fully implicit algorithm to solve the conservative full-potential equation for the steady transonic flow field about an arbitrary airfoil immersed in a subsonic free stream. The full-potential formulation is considered exact under the assumptions of irrotational, isentropic, and inviscid flow. These assumptions are valid for a wide range of practical transonic flows typical of modern aircraft cruise conditions. The primary features of TAIR include: a new fully implicit iteration scheme which is typically many times faster than classical successive line overrelaxation algorithms; a new, reliable artifical density spatial differencing scheme treating the conservative form of the full-potential equation; and a numerical mapping procedure capable of generating curvilinear, body-fitted finite-difference grids about arbitrary airfoil geometries. Three aspects emphasized during the development of the TAIR code were reliability, simplicity, and speed. The reliability of TAIR comes from two sources: the new algorithm employed and the implementation of effective convergence monitoring logic. TAIR achieves ease of use by employing a "default mode" that greatly simplifies code operation, especially by inexperienced users, and many useful options including: several airfoil-geometry input options, flexible user controls over program output, and a multiple solution capability. The speed of the TAIR code is attributed to the new algorithm and the manner in which it has been implemented. Input to the TAIR program consists of airfoil coordinates, aerodynamic and flow-field convergence parameters, and geometric and grid convergence parameters. The airfoil coordinates for many airfoil shapes can be generated in TAIR from just a few input parameters. Most of the other input parameters have default values which allow the user to run an analysis in the default mode by specifing only a few input parameters. Output from TAIR may include aerodynamic coefficients, the airfoil surface solution, convergence histories, and printer plots of Mach number and density contour maps. The TAIR program is written in FORTRAN IV for batch execution and has been implemented on a CDC 7600 computer with a central memory requirement of approximately 155K (octal) of 60 bit words. The TAIR program was developed in 1981.

Description: The TAWFIVE program calculates transonic flow over a transport-type wing and fuselage. Although more complex Euler and Navier-Stokes methods are available, TAWFIVE combines a multi-grid acceleration technique in the iterative solution of the potential equation with the use of integral-form boundary-layer equations to provide a computationally efficient and sufficiently accurate design tool. TAWFIVE simplifies the solution process by breaking the problem into a loosely coupled set of modified equations. The inviscid method, using standard inviscid equations (nonlinear full potential), is valid in the "outer" region away from the wing, whereas the boundary-layer equations are valid in the thin region near the solid surface of the wing. The two types of equations are coupled by a technique of modifying surface boundary conditions for the inviscid equations. This interaction process starts with a solution of the outer flow field. Pressures are computed at the wing surface and are used to calculate the boundary layer. The boundary-layer and wake properties are then computed using a three-dimensional integral method, and the computed displacement thickness is added to the surface of the "hard" geometry. This new displaced wing surface is then regridded and the inviscid flowfield is recomputed. New values of the inviscid pressures are then used by the boundary-layer method to predict a new displacement thickness distribution. An under-relaxed update of the previously predicted displacement thickness is then made to obtain a new displacement thickness correction that is added to the "hard" geometry. These global iterations are continued until suitable convergence is obtained. Input to TAWFIVE is limited to geometric definition of the configuration, free-stream flow quantities, and iteration control parameters. The geometric input consists of the definition of a series of airfoil sections to define the wing and a series of fuselage cross sections to model the fuselage. High-aspect-ratio wings are modeled more accurately than low-aspect-ratio wings since no special provisions are made to accurately model the wing-fuselage juncture or the wingtip region. The user can specify the solution either in terms of lift or in terms of angle of attack. TAWFIVE can produce tabular output and input files for PLOT3D (COSMIC program number ARC-12779). TAWFIVE is written in FORTRAN 77 for CRAY series computers running UNICOS. The main memory requirement is 2.7Mb for execution. This program is available on a 9-track 1600 BPI UNIX tar format magnetic tape. TAWFIVE was under development from 1979 to 1989 and first released by COSMIC in 1991. CRAY and UNICOS are registered trademarks of Cray Research, Inc.

Description: A summary is presented of vortex control applications and current techniques for the control of longitudinal vortices produced by bodies, leading edges, tips and intersections. Vortex control has up till now been performed by many approaches in an empirical fashion, assisted by the essentially inviscid nature of much of longitudinal vortex behavior. Attention is given to Reynolds number sensitivities, vortex breakdown and interactions, vortex control on highly swept wings, and vortex control in juncture flows.

Description: The flow-field within an axial flow turbomachine, such as a turbine or compressor, is extremely complex because of three-dimensional features such as hub-corner stall, tip-leakage flows, and airfoil wakes. These flow features interact with each other and with rotor and stator airfoils inducing time-varying forces on the airfoils. These complicated rotor-stator interactions must be understood in order to design turbomachines that are light and compact as well as reliable and efficient. Two codes, STAGE-2 and STAGE-3, have been developed to compute these unsteady rotor-stator interaction flows in multistage turbomachines. An implicit, thin-layer Euler/Navier-Stokes zonal algorithm is used to compute the unsteady flow-field within both turbine and compressor configurations. Results include surface pressures and wake profiles for two-dimensional turbine and compressor configurations and surface pressures for a three-dimensional single-stage turbine configuration. The results compare well with experimental data and other unsteady computations.

Description: Jet noise and jet-induced structural loads have become key issues in the design of commercial and military aircraft. Computational Fluid Dynamics (CFD) can be of use in predicting the underlying jet shear-layer instabilities and, in conjunction with classical acoustic theory, jet noise. The computational issues involved in the resolution of high Reynolds number unsteady jet flows are addressed in this paper. Once these jet flows can be accurately resolved, it should be possible to use acoustic theory to extract, for example, the far-field jet noise. An assessment of future work and computational resources required for directly computing far-field jet noise is also presented.

Description: The signal from a stable periodicity can seem to be intermittent when it is partially masked by an unmodelled window function or when the data set is too short to resolve closely spaced periodicities. By taking this into account, short-lived periodicities in solar data can be reinterpreted as evidence for continuously periodic behavior. The periodic sources are located in the solar interior and caused by global oscillation modes. The convective envelope acts as the window for these sources. Recent reports of seven periodicities from 100 to 1000 days are compared with this model. Precise long-term values for the periodicities are predicted and they agree closely with observations. Some elements are suggested that might explain the well-documented 155-day periodicity. Conventional filtering methods to suppress effects of the 11-year cycle are criticized as inadequate.

Description: Data on solar wind electrons at the ISEE 3 spacecraft located 0.01 AU upstream from the earth (McComas et al., 1989) showed periods of time when the flux of antisunward suprathermal electrons would decrease suddenly, leading to heat flux dropouts (HFDs). This paper examines data from ISEE 1 at the 1.5 x 10 exp 6 km downstream location to determine whether HFDs identified at ISEE 3 by McComas et al. can be detected at this location and whether the ISEE 1 observations can provide information to one or the other possible interpretations of HFDs: that HFDs are due to enhanced Coulomb scattering, or to disconnection from the sun of the magnetic flux tube. The results of the examination identified the presence of HFD events in the ISEE 1 data, and the findings indicate that Coulomb scattering plays a substantial role in at least some HFD events.

Description: This computer program is designed to calculate the flow fields in two-dimensional and three-dimensional axisymmetric supersonic inlets. The method of characteristics is used to compute arrays of points in the flow field. At each point the total pressure, local Mach number, local flow angle, and static pressure are calculated. This program can be used to design and analyze supersonic inlets by determining the surface compression rates and throat flow properties. The program employs the method of characteristics for a perfect gas. The basic equation used in the program is the compatibility equation which relates the change in stream angle to the change in entropy and the change in velocity. In order to facilitate the computation, the flow field behind the bow shock wave is broken into regions bounded by shock waves. In each region successive rays are computed from a surface to a shock wave until the shock wave intersects a surface or falls outside the cowl lip. As soon as the intersection occurs a new region is started and the previous region continued only in the area in which it is needed, thus eliminating unnecessary calculations. The maximum number of regions possible in the program is ten, which allows for the simultaneous calculations of up to nine shock waves. Input to this program consists of surface contours, free-stream Mach number, and various calculation control parameters. Output consists of printed and/or plotted results. For plotted results an SC-4020 or similar plotting device is required. This program is written in FORTRAN IV to be executed in the batch mode and has been implemented on a CDC 7600 with a central memory requirement of approximately 27k (octal) of 60 bit words.

Description: This program was developed to predict turbine stage performance taking into account the effects of complex passage geometries. The method uses a quasi-3D inviscid-flow analysis iteratively coupled to calculated losses so that changes in losses result in changes in the flow distribution. In this manner the effects of both the geometry on the flow distribution and the flow distribution on losses are accounted for. The flow may be subsonic or shock-free transonic. The blade row may be fixed or rotating, and the blades may be twisted and leaned. This program has been applied to axial and radial turbines, and is helpful in the analysis of mixed flow machines. This program is a combination of the flow analysis programs MERIDL and TSONIC coupled to the boundary layer program BLAYER. The subsonic flow solution is obtained by a finite difference, stream function analysis. Transonic blade-to-blade solutions are obtained using information from the finite difference, stream function solution with a reduced flow factor. Upstream and downstream flow variables may vary from hub to shroud and provision is made to correct for loss of stagnation pressure. Boundary layer analyses are made to determine profile and end-wall friction losses. Empirical loss models are used to account for incidence, secondary flow, disc windage, and clearance losses. The total losses are then used to calculate stator, rotor, and stage efficiency. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370/3033 under TSS with a central memory requirement of approximately 4.5 Megs of 8 bit bytes. This program was developed in 1985.

Description: Turbomachinery components are often connected by ducts, which are usually annular. The configurations and aerodynamic characteristics of these ducts are crucial to the optimum performance of the turbomachinery blade rows. The ANDUCT computer program was developed to calculate the velocity distribution along an arbitrary line between the inner and outer walls of an annular duct with axisymmetric swirling flow. Although other programs are available for duct analysis, the use of the velocity gradient method makes the ANDUCT program fast and convenient while requiring only modest computer resources. A fast and easy method of analyzing the flow through a duct with axisymmetric flow is the velocity gradient method, also known as the stream filament or streamline curvature method. This method has been used extensively for blade passages but has not been widely used for ducts, except for the radial equilibrium equation. In ANDUCT, a velocity gradient equation derived from the momentum equation is used to determine the velocity variation along an arbitrary straight line between the inner and outer wall of an annular duct. The velocity gradient equation is used with an assumed variation of meridional streamline curvature. Upstream flow conditions may vary between the inner and outer walls, and an assumed total pressure distribution may be specified. ANDUCT works best for well-guided passages and where the curvature of the walls is small as compared to the width of the passage. The ANDUCT program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 60K of 8 bit bytes. The ANDUCT program was developed in 1982.

Description: The Panel Code for Planar Cascades was developed as an aid for the designer of turbomachinery blade rows. The effective design of turbomachinery blade rows relies on the use of computer codes to model the flow on blade-to-blade surfaces. Most of the currently used codes model the flow as inviscid, irrotational, and compressible with solutions being obtained by finite difference or finite element numerical techniques. While these codes can yield very accurate solutions, they usually require an experienced user to manipulate input data and control parameters. Also, they often limit a designer in the types of blade geometries, cascade configurations, and flow conditions that can be considered. The Panel Code for Planar Cascades accelerates the design process and gives the designer more freedom in developing blade shapes by offering a simple blade-to-blade flow code. Panel, or integral equation, solution techniques have been used for several years by external aerodynamicists who have developed and refined them into a primary design tool of the aircraft industry. The Panel Code for Planar Cascades adapts these same techniques to provide a versatile, stable, and efficient calculation scheme for internal flow. The code calculates the compressible, inviscid, irrotational flow through a planar cascade of arbitrary blade shapes. Since the panel solution technique is for incompressible flow, a compressibility correction is introduced to account for compressible flow effects. The analysis is limited to flow conditions in the subsonic and shock-free transonic range. Input to the code consists of inlet flow conditions, blade geometry data, and simple control parameters. Output includes flow parameters at selected control points. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 590K of 8 bit bytes. This program was developed in 1982.

Description: An exact, full-potential-equation model for the steady, irrotational, homoentropic, and homoenergetic flow of a compressible, inviscid fluid through a two-dimensional planar cascade together with its appropriate boundary conditions has been derived. The CAS2D computer program numerically solves an artificially time-dependent form of the actual full-potential-equation, providing a nonrotating blade-to-blade, steady, potential transonic cascade flow analysis code. Comparisons of results with test data and theoretical solutions indicate very good agreement. In CAS2D, the governing equation is discretized by using type-dependent, rotated finite differencing and the finite area technique. The flow field is discretized by providing a boundary-fitted, nonuniform computational mesh. This mesh is generated by using a sequence of conformal mapping, nonorthogonal coordinate stretching, and local, isoparametric, bilinear mapping functions. The discretized form of the full-potential equation is solved iteratively by using successive line over relaxation. Possible isentropic shocks are captured by the explicit addition of an artificial viscosity in a conservative form. In addition, a four-level, consecutive, mesh refinement feature makes CAS2D a reliable and fast algorithm for the analysis of transonic, two-dimensional cascade flows. The results from CAS2D are not directly applicable to three-dimensional, potential, rotating flows through a cascade of blades because CAS2D does not consider the effects of the Coriolis force that would be present in the three-dimensional case. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 200K of 8 bit bytes. The CAS2D program was developed in 1980.

Description: A computer program, QSONIC, has been developed for calculating the full potential, transonic quasi-three-dimensional flow through a rotating turbomachinery blade row. The need for lighter, more efficient turbomachinery components has led to the consideration of machines with fewer stages, each with blades capable of higher speeds and higher loading. As speeds increase, the numerical problems inherent in the transonic regime have to be resolved. These problems include the calculation of imbedded shock discontinuities and the dual nature of the governing equations, which are elliptic in the subcritical flow regions but become hyperbolic for supersonic zones. QSONIC provides the flow analyst with a fast and reliable means of obtaining the transonic potential flow distribution on a blade-to-blade stream surface of a stationary or rotating turbomachine blade row. QSONIC combines several promising transonic analysis techniques. The full potential equation in conservative form is discretized at each point on a body-fitted period mesh. A mass balance is calculated through the finite volume surrounding each point. Each local volume is corrected in the third dimension for any change in stream-tube thickness along the stream tube. The nonlinear equations for all volumes are of mixed type (elliptic or hyperbolic) depending on the local Mach number. The final result is a block-tridiagonal matrix formulation involving potential corrections at each grid point as the unknowns. The residual of each system of equations is solved along each grid line. At points where the Mach number exceeds unity, the density at the forward (sweeping) edge of the volume is replaced by an artificial density. This method calculates the flow field about a cascade of arbitrary two-dimensional airfoils. Three-dimensional flow is approximated in a turbomachinery blade row by correcting for stream-tube convergence and radius change in the through flow direction. Several significant assumptions were made in developing the QSONIC program, including: (1) the flow is inviscid and adiabatic, (2) the flow relative to the blade is steady, (3) the fluid is a perfect gas with constant specific heat, (4) the flow is isentropic and any discontinuities (shocks) are weak enough to be approximated as isentropic jumps, (5) there is no velocity component normal to the stream surface, and (6) the flow relative to a fixed frame in space (absolute velocity) is completely irrotational. These assumptions place some limitations on the application of QSONIC. Sharp leading edges at high incidence and high-Mach-number turbine blade trailing edges with substantial deviation will both cause large velocity peaks on the blade. In addition, the program may have difficulty converging if the passage is nearly choked. Input to QSONIC consists of case control parameters, a geometry description, upstream boundary conditions, and a rotor description. Output includes solution scheme parameters and flow field parameters. A data file is also output which contains data on the solution mesh, surface Mach numbers, surface static pressures, isomachs, and the velocity vector field. This data may be used for further processing or for plotting. The QSONIC is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 500K of 8 bit bytes. QSONIC was developed in 1982.

Description: This computer program, WIND, was developed to numerically solve the exact, full-potential equation for three-dimensional, steady, inviscid flow through an isolated wind turbine rotor. The program automatically generates a three-dimensional, boundary-conforming grid and iteratively solves the full-potential equation while fully accounting for both the rotating and Coriolis effects. WIND is capable of numerically analyzing the flow field about a given blade shape of the horizontal-axis type wind turbine. The rotor hub is assumed representable by a doubly infinite circular cylinder. An arbitrary number of blades may be attached to the hub and these blades may have arbitrary spanwise distributions of taper and of the twist, sweep, and dihedral angles. An arbitrary number of different airfoil section shapes may be used along the span as long as the spanwise variation of all the geometeric parameters is reasonably smooth. The numerical techniques employed in WIND involve rotated, type-dependent finite differencing, a finite volume method, artificial viscosity in conservative form, and a successive overrelaxation combined with the sequential grid refinement procedure to accelerate the iterative convergence rate. Consequently, WIND is cabable of accurately analyzing incompressible and compressible flows, including those that are locally transonic and terminated by weak shocks. Along with the three-dimensional results, WIND provides the results of the two-dimensional calculations to aid the user in locating areas of possible improvement in the aerodynamic design of the blade. Output from WIND includes the chordwise distribution of the coefficient of pressure, the Mach number, the density, and the relative velocity components at spanwise stations along the blade. In addition, the results specify local values of the lift coefficient and the tangent and axial aerodynamic force components. These are also given in integrated form expressing the total torque and the total axial force acting on the shaft. WIND can also be used to analyze the flow around isolated aircraft propellers and helicopter rotors in hover as long as the relative oncoming flow is subsonic. The WIND program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 253K of 8 bit bytes. WIND was developed in 1980.

Description: This computer program calculates the flow field in the supersonic portion of a mixed-compression aircraft inlet at non-zero angle of attack. This approach is based on the method of characteristics for steady three-dimensional flow. The results of this program agree with those produced by the two-dimensional method of characteristics when axisymmetric flow fields are calculated. Except in regions of high viscous interaction and boundary layer removal, the results agree well with experimental data obtained for threedimensional flow fields. The flow field in a variety of axisymmetric mixed compression inlets can be calculated using this program. The bow shock wave and the internal shock wave system are calculated using a discrete shock wave fitting procedure. The internal flow field can be calculated either with or without the discrete fitting of the internal shock wave system. The influence of molecular transport can be included in the calculation of the external flow about the forebody and in the calculation of the internal flow when internal shock waves are not discretely fitted. The viscous and thermal diffussion effects are included by treating them as correction terms in the method of characteristics procedure. Dynamic viscosity is represented by Sutherland's law and thermal conductivity is represented as a quadratic function of temperature. The thermodynamic model used is that of a thermally and calorically perfect gas. The program assumes that the cowl lip is contained in a constant plane and that the centerbody contour and cowl contour are smooth and have continuous first partial derivatives. This program cannot calculate subsonic flow, the external flow field if the bow shock wave does not exist entirely around the forebody, or the internal flow field if the bow flow field is injected into the annulus. Input to the program consists of parameters to control execution, to define the geometry, and the vehicle orientation. Output consists of a list of parameters used, solution planes, and a description of the shock waves. This program is written in FORTRAN IV for batch execution and has been implemented on a CDC 6000 series machine with a central memory requirement of 110K (octal) of 60 bit words when it is overlayed. This flow analysis program was developed in 1978.

Description: This computer program was developed for calculating the subsonic or transonic flow on the hub-shroud mid-channel stream surface of a single blade row of a turbomachine. The design and analysis of blades for compressors and turbines ideally requires methods for analyzing unsteady, three-dimensional, turbulent viscous flow through a turbomachine. Since an exact solution is impossible at present, solutions on two-dimensional surfaces are calculated to obtain a quasi-three dimensional solution. When three-dimensional effects are important, significant information can be obtained from a solution on a cross-sectional surface of the passage normal to the flow. With this program, a solution to the equations of flow on the meridional surface can be carried out. This solution is chosen when the turbomachine under consideration has significant variation in flow properties in the hubshroud direction, especially when input is needed for use in blade-to-blade calculations. The program can also perform flow calculations for annular ducts without blades. This program should prove very useful in the design and analysis of any turbomachine. This program calculates a solution for two-dimensional, adiabatic shockfree flow. The flow must be essentially subsonic, but there may be local areas of supersonic flow. To obtain the solution, this program uses both the finite difference and the quasi-orthogonal (velocity gradient) methods combined in a way that takes maximum advantage of both. The finite-difference method solves a finite-difference equation along the meridional stream surface in a very efficient manner but is limited to subsonic velocities. This approach must be used in cases where the blade aspect ratios are above one, cases where the passage is curved, and cases with low hub-tip-ratio blades. The quasi-orthogonal method solves the velocity gradient equation on the meridional surface and is used if it is necessary to extend the range of solutions into the transonic regime. In general the blade row may be fixed or rotating and the blades may be twisted and leaned. The flow may be axial, radial, or mixed. The upstream and downstream flow conditions can vary from hub to shroud with provisions made for an approximate correction for loss of stagnation pressure. Also, viscous forces are neglected along solution mesh lines running from hub to tip. The capabilities of this program include handling of nonaxial flows without restriction, annular ducts without blades, and specified streamwise loss distributions. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 360 computer with a central memory requirement of approximately 700K of 8 bit bytes. This core requirement can be reduced depending on the size of the problem and the desired solution accuracy. This program was developed in 1977.

Description: A computer program has been developed for the design of supersonic rotor blades where losses are accounted for by correcting the ideal blade geometry for boundary layer displacement thickness. The ideal blade passage is designed by the method of characteristics and is based on establishing vortex flow within the passage. Boundary-layer parameters (displacement and momentum thicknesses) are calculated for the ideal passage, and the final blade geometry is obtained by adding the displacement thicknesses to the ideal nozzle coordinates. The boundary-layer parameters are also used to calculate the aftermixing conditions downstream of the rotor blades assuming the flow mixes to a uniform state. The computer program input consists essentially of the rotor inlet and outlet Mach numbers, upper- and lower-surface Mach numbers, inlet flow angle, specific heat ratio, and total flow conditions. The program gas properties are set up for air. Additional gases require changes to be made to the program. The computer output consists of the corrected rotor blade coordinates, the principal boundary-layer parameters, and the aftermixing conditions. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 7094. This program was developed in 1971.

Description: This program obtains a transonic flow solution on a blade-to-blade surface between blades of a turbomachine. The flow must be essentially subsonic, but there may be locally supersonic flow. The solution is two-dimensional, isentropic, and shock free. The blades may be fixed or rotating. The flow may be axial, radial, or mixed, and there may be a change in stream-channel thickness in the through-flow direction. A loss in relative stagnation pressure may be accounted for. The program input consists of blade and stream-channel geometry, stagnation flow conditions, inlet and outlet flow angles, and blade-to-blade stream-channel weight flow. The output includes blade surface velocities, velocity magnitude and direction at all interior mesh points in the blade-to-blade passage, and streamline coordinates throughout the passage. The transonic solution is obtained by a combination of a finite-difference, stream-function solution and a velocity-gradient solution. The finite-difference solution at a reduced weight flow provides information needed to obtain a velocity-gradient solution. This program is written in FORTRAN IV for batch execution and has been implemented on the IBM 360 computer with a central memory requirement of approximately 36K of 8 bit bytes. This program was developed in 1969 and last updated in 1979.

Description: This program is a revision of an existing program for blade-to-blade aerodynamic analysis of turbomachine blades and it is a simpler program while consistent with related programs. The analysis is for two-dimensional, subsonic, compressible (or incompressible), nonviscous flow in a circular or straight infinite cascade of blades, which may be fixed or rotating. The flow may be axial, radial, or mixed, and the stream channel thickness may change in the through-flow direction. The program input consists of blade and stream channel geometry, total flow conditions, inlet and outlet flow angles, and blade-to-blade stream channel weight flow. The output includes blade surface velocities, velocity magnitude and direction at all interior mesh points in the blade-to-blade passage, and streamline coordinates throughout the passage. This program was developed on an IBM 7094/7044 DCS.

Description: This computer program gives the blade-to-blade solution of the two-dimensional, subsonic, compressible (or incompressible), nonviscous flow problem for a circular or straight infinite cascade of tandem or slotted turbomachine blades. The blades may be fixed or rotating. The flow may be axial, radial , or mixed. The method of solution is based on the stream function using an iterative solution of nonlinear finite-difference equations. These equations are solved using two major levels of iteration. The inner iteration consists of the solution of simultaneous linear equations by successive over-relaxation, using an estimated optimum over-relaxation factor. The outer iteration then changes the coefficients of the simultaneous equations to correct for compressibility. The program input consists of the basic blade geometry, the meridional stream channel coordinates, fluid stagnation conditions, weight flow and flow split through the slot, and inlet and outlet flow angles. The output includes blade surface velocities, velocity magnitude and direction throughout the passage, and the streamline coordinates.

Description: This FORTRAN IV computer program which incorporates the method of characteristics was written to assist in the design of supersonic inlets. There were two objectives: (1) to study a greater variety of supersonic inlet configurations and (2) to reduce the time required for trial-and-error procedures to arrive at optimum inlet design. The computer program was written with the intention of being able to construct a variety of inlet configurations by interchanging specific subroutines. In this manner, greater flexibility of choice was attained, and the time required to program a specific inlet configuration was greatly reduced. The second objective was accomplished by a reformulation of the boundary value problem for hyperbolic equations. By this reformulation of the boundary data, the engineering design quantities, throat Mach number and flow angle, were introduced as direct input quantities to the computer program. As a consequence of introducing the engineering parameters as input, the computer program will calculate the surface contours required to satisfy the specific throat conditions. Inviscid flow is assumed and the method used to calculate the inlet contour results in minimum distortion to the flow in the throat. This program was developed on an IBM 7094.

Description: This program represents a subsonic aerodynamic method for determining the mean camber surface of trimmed noncoplaner planforms with minimum vortex drag. With this program, multiple surfaces can be designed together to yield a trimmed configuration with minimum induced drag at some specified lift coefficient. The method uses a vortex-lattice and overcomes previous difficulties with chord loading specification. A Trefftz plane analysis is used to determine the optimum span loading for minimum drag. The program then solves for the mean camber surface of the wing associated with this loading. Pitching-moment or root-bending-moment constraints can be employed at the design lift coefficient. Sensitivity studies of vortex-lattice arrangements have been made with this program and comparisons with other theories show generally good agreement. The program is very versatile and has been applied to isolated wings, wing-canard configurations, a tandem wing, and a wing-winglet configuration. The design problem solved with this code is essentially an optimization one. A subsonic vortex-lattice is used to determine the span load distribution(s) on bent lifting line(s) in the Trefftz plane. A Lagrange multiplier technique determines the required loading which is used to calculate the mean camber slopes, which are then integrated to yield the local elevation surface. The problem of determining the necessary circulation matrix is simplified by having the chordwise shape of the bound circulation remain unchanged across each span, though the chordwise shape may vary from one planform to another. The circulation matrix is obtained by calculating the spanwise scaling of the chordwise shapes. A chordwise summation of the lift and pitching-moment is utilized in the Trefftz plane solution on the assumption that the trailing wake does not roll up and that the general configuration has specifiable chord loading shapes. VLMD is written in FORTRAN for IBM PC series and compatible computers running MS-DOS. This program requires 360K of RAM for execution. The Ryan McFarland FORTRAN compiler and PLINK86 are required to recompile the source code; however, a sample executable is provided on the diskette. The standard distribution medium for VLMD is a 5.25 inch 360K MS-DOS format diskette. VLMD was originally developed for use on CDC 6000 series computers in 1976. It was originally ported to the IBM PC in 1986, and, after minor modifications, the IBM PC port was released in 1993.

Description: This code was developed to aid design engineers in the selection and evaluation of aerodynamically efficient wing-canard and wing-horizontal-tail configurations that may employ simple hinged-flap systems. Rapid estimates of the longitudinal aerodynamic characteristics of conceptual airplane lifting surface arrangements are provided. The method is particularly well suited to configurations which, because of high speed flight requirements, must employ thin wings with highly swept leading edges. The code is applicable to wings with either sharp or rounded leading edges. The code provides theoretical pressure distributions over the wing, the canard or horizontal tail, and the deflected flap surfaces as well as estimates of the wing lift, drag, and pitching moments which account for attainable leading edge thrust and leading edge separation vortex forces. The wing planform information is specified by a series of leading edge and trailing edge breakpoints for a right hand wing panel. Up to 21 pairs of coordinates may be used to describe both the leading edge and the trailing edge. The code has been written to accommodate 2000 right hand panel elements, but can easily be modified to accommodate a larger or smaller number of elements depending on the capacity of the target computer platform. The code provides solutions for wing surfaces composed of all possible combinations of leading edge and trailing edge flap settings provided by the original deflection multipliers and by the flap deflection multipliers. Up to 25 pairs of leading edge and trailing edge flap deflection schedules may thus be treated simultaneously. The code also provides for an improved accounting of hinge-line singularities in determination of wing forces and moments. To determine lifting surface perturbation velocity distributions, the code provides for a maximum of 70 iterations. The program is constructed so that successive runs may be made with a given code entry. To make additional runs, it is necessary only to add an identification record and the namelist data that are to be changed from the previous run. This code was originally developed in 1989 in FORTRAN V on a CDC 6000 computer system, and was later ported to an MS-DOS environment. Both versions are available from COSMIC. There are only a few differences between the PC version (LAR-14458) and CDC version (LAR-14178) of AERO2S distributed by COSMIC. The CDC version has one main source code file while the PC version has two files which are easier to edit and compile on a PC. The PC version does not require a FORTRAN compiler which supports NAMELIST because a special INPUT subroutine has been added. The CDC version includes two MODIFY decks which can be used to improve the code and prevent the possibility of some infrequently occurring errors while PC-version users will have to make these code changes manually. The PC version includes an executable which was generated with the Ryan McFarland/FORTRAN compiler and requires 253K RAM and an 80x87 math co-processor. Using this executable, the sample case requires about four hours to execute on an 8MHz AT-class microcomputer with a co-processor. The source code conforms to the FORTRAN 77 standard except that it uses variables longer than six characters. With two minor modifications, the PC version should be portable to any computer with a FORTRAN compiler and sufficient memory. The CDC version of AERO2S is available in CDC NOS Internal format on a 9-track 1600 BPI magnetic tape. The PC version is available on a set of two 5.25 inch 360K MS-DOS format diskettes. IBM AT is a registered trademark of International Business Machines. MS-DOS is a registered trademark of Microsoft Corporation. CDC is a registered trademark of Control Data Corporation. NOS is a trademark of Control Data Corporation.

Description: This program provides a wing design algorithm based on modified linear theory which takes into account the effects of attainable leading-edge thrust. A primary objective of the WINGDES2 approach is the generation of a camber surface as mild as possible to produce drag levels comparable to those attainable with full theoretical leading-edge thrust. WINGDES2 provides both an analysis and a design capability and is applicable to both subsonic and supersonic flow. The optimization can be carried out for designated wing portions such as leading and trailing edge areas for the design of mission-adaptive surfaces, or for an entire planform such as a supersonic transport wing. This program replaces an earlier wing design code, LAR-13315, designated WINGDES. WINGDES2 incorporates modifications to improve numerical accuracy and provides additional capabilities. A means of accounting for the presence of interference pressure fields from airplane components other than the wing and a direct process for selection of flap surfaces to approach the performance levels of the optimized wing surfaces are included. An increased storage capacity allows better numerical representation of those configurations that have small chord leading-edge or trailing-edge design areas. WINGDES2 determines an optimum combination of a series of candidate surfaces rather than the more commonly used candidate loadings. The objective of the design is the recovery of unrealized theoretical leading-edge thrust of the input flat surface by shaping of the design surface to create a distributed thrust and thus minimize drag. The input consists of airfoil section thickness data, leading and trailing edge planform geometry, and operational parameters such as Mach number, Reynolds number, and design lift coefficient. Output includes optimized camber surface ordinates, pressure coefficient distributions, and theoretical aerodynamic characteristics. WINGDES2 is written in FORTRAN V for batch execution and has been implemented on a CDC CYBER computer operating under NOS 2.7.1 with a central memory requirement of approximately 344K (octal) of 60 bit words. This program was developed in 1984, and last updated in 1990. CDC and CYBER are trademarks of Control Data Corporation.

Description: Unsteady flow behavior and load characteristics of a 2D VR-7 airfoil with and without a leading-edge slat were studied in the water tunnel of the Aeroflightdynamics Directorate, NASA Ames Research Center. Both airfoils were oscillated sinusoidally between 5 and 25 deg at Re = 200,000 to obtain the unsteady lift, drag, and pitching moment data. A fluorescent dye was released from an orifice located at the leading edge of the airfoil for the purpose of visualizing the boundary layer and wake flow. The flowfield and load predictions of an incompressible Navier-Stokes code based on a velocity-vorticity formulation were compared with the test data. The test and predictions both confirm that the slatted VR-7 airfoil delays both static and dynamic stall as compared to the VR-7 airfoil alone.

Description: The dynamics of radially expanding magnetic clouds is rigorously analyzed within the framework of ideal MHD. The cloud is modelled as a cylindrically symmetric magnetic flux rope. In the force balance we include the gas pressure gradient and the Lorentz force. Interaction with the ambient solar wind due to expansion of the magnetic cloud is represented by a drag force proportional to the bulk velocity. We consider the self-similar expansion of a polytrope, and reduce the problem to an ordinary nonlinear differential equation for the evolution function. Analyzing the asymptotic behavior of the evolution function, we formulate theoretical expectations for the long-term behavior of cloud parameters. We focus on the temporal evolution of (1) the magnetic field strength; (2) the twist of the field lines; (3) the asymmetry of the total field profile; and (4) the bulk flow speed. We present data from two magnetic clouds observed at 1 AU and 2 AU, respectively, and find good agreement with theoretical expectations. For a peak magnetic field strength at 1 AU of 25 nT and a polytropic index of 0.5, we find that a magnetic cloud can be distinguished from the background interplanetary field up to a distance of about 5 AU. Taking larger magnetic fields and bigger polytropic indices this distance can double.

Description: Using unique 3-d velocity space measurements by the Ulysses solar wind plasma experiment from 1.15 to 5.34 AU, we assess the radial gradient in thermal electron temperature. Until 3.8 AU, the gradient was steeper than previously reported but flatter than adiabatic; after 3.8 AU the gradient flattened. Trends in the observed electron distribution shapes qualitatively support predictions for regulation by Coulomb collisions and by expansion in a spiral IMF.

Description: We use the horizontal momentum balance equation to infer the strength of the meridional circulation (MC) and Reynolds stresses (RS) at the sun's surface from the observed properties of the differential rotation (DR). Both MC and RS are important for maintaining the equatorial acceleration. The results indicate that the average value of MC is about 1.1 m/s, with circulation directed towards the poles in both the hemispheres, and the average value of RS is 3.6 10 exp 7 sq cm2/s, with transport of angular momentum directed towards the equator in both hemispheres, this latter in good agreement with observations. With the above values of MC and RS, we integrate the momentum equation in time, starting from a state of rigid rotation, to investigate the competitive role of MC and RS in producing the presently observed average DR. Our results show that DR is consistent with observations only if both MC and RS have opposite effects, with the strengths given above, in order to balance the viscous torque.

Description: The problem of the hypersonic double ellipse in rarefied flow is treated by a particle method using the collision model first described by McDonald (1988). In the approach used here, the computational overhead is reduced by using simple cubic cells. The problem of the definition of complex geometries is addressed by developing an algorithm to define the relation of a body surface to the network of cells.

Description: Program LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is an upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3D viscous hypersonic flows in chemical and thermal nonequilibrium. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with a modified Roe's averaging and with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The code has been applied to Problem 8.2 of this workshop for the case of thermochemical nonequilibrium flow through a nozzle. Chemical reaction rates are defined with the model of Park (1987). Thermal nonequilibrium is modeled using a two-temperature approximation in which the vibrational energies of all molecules are assumed to be in equilibrium at a single temperature which is generally different from the translational-rotational temperature. Two grids were used to define the flow for the original problem, with a stagnation temperature of 6500 K. A third case with a stagnation temperature of 10,000 K is also presented. The solution domain includes the converging nozzle, subsonic flow domain in which the gas is substantially in thermochemical equilibrium and the diverging nozzle, hypersonic flow domain in which the gas is substantially in thermochemical nonequilibrium.

Description: Solutions have been computed and results are presented for Problem 1, the case of Mach 9 transitional flow past a 7 deg half-angle cone at zero incidence. The solutions were computed using a code developed for the integration of the parabolized Navier-Stokes equations. The algorithm employed in the code is based on a Roe-type flux-difference-splitting scheme applied following a finite-volume approach. The basic algorithm has been modified to make it implicit and second-order accurate in the crossflow directions. Results are presented in terms of surface pressure and heat transfer as well as boundary layer profiles of pitot pressure, Mach number, and tangential velocity. The case was recalculated several times in an effort to determine sensitivities to such parameters as grid density, wall temperature, turbulence model parameters, as well as freestream expansion. Comparisons with the experimental data are presented and discussed.

Description: A development status evaluation is presented for the theoretical understanding and design conceptualization of boundary layer control (BLC) systems applicable to supersonic transports, such as the currently envisioned NASA High Speed Civil Transport. By reducing fuel burned, supersonic BLC techniques could expand ranges to Pacific-crossing scales, while lowering sonic boom effects and upper-atmosphere pollution and even reducing skin friction temperature. The critical consideration for supersonic BLC is the presence of wave effects.

Description: Attention is drawn to the influence of preexisting finite-amplitude instabilities on the growth of other disturbances; current design tools for LFC take no notice of this kind of interaction. When a rational accounting is accomplished for the evolution of incoming disturbances in finite-amplitude solutions of the equations of motion, future transition-prediction methods will need to take these wave interactions into account. Attention is given here to interactions in the presence of crossflow vortices and interactions involving Goertler vortices.

Description: A development history and a development-trends evaluation are presented for laminar flow controlled airfoil technologies and design concepts, including the search for 'natural' laminar flow and actively controlled flow via suction through small pores on the airfoil surface. While most NASA activities in this field have been concerned with subsonic aircraft, it has been projected that the control of boundary layer turbulence may be even more critical to the aerodynamic efficiency of supersonic aircraft. Developmental programs for these techniques have been conducted with several modified conventional aircraft.

Description: An account is given of the development history of natural laminar-flow (NLF) airfoil profiles under guidance of an experimentally well-verified theoretical method for the design of airfoils suited to virtually all subcritical applications. This method, the Eppler Airfoil Design and Analysis Program, contains a conformal-mapping method for airfoils having prescribed velocity-distribution characteristics, as well as a panel method for the analysis of potential flow about given airfoils and a boundary-layer method. Several of the NLF airfoils thus obtained are discussed.

Description: We correlate EUV spectroheliograms of the quiet sun taken by the Harvard EUV Spectroheliometer on Skylab with a magnetogram taken within 15-20 minutes by the 40 channel magnetograph at Kitt Peak. Focus is on the O VI and Ne VII ions which are formed in the hotter transition region and the Mg x ion which is formed in the corona. Since no imaged observations in lines from T = 2 x 10 exp 5 - 10 exp 6 K temperature region have been made since Skylab, these data are unique in providing opportunities to investigate the hotter transition region and its relationship to the corona. The signature of an EUV intranetwork loop is a small EUV brightening located on the neutral line of a small magnetic bipole within the supergranular network. The data clearly show such features. Some of these features are visible at coronal temperatures, indicating that network loops can reach 10 exp 6 K. Others are bright in the lines formed below 10 exp 6 K but occur in regions of depleted coronal emission, which cannot be explained by models based on back-heating from a large-scale corona. We conclude that some fraction of the quiet solar EUV output is generated within network loops that are effectively insulated from the corona.

Description: We consider the gyroresonant interaction of protons with parallel electromagnetic plasma waves. These waves have either right- or left-hand circular polarization and include as a subset Alfven and whistler waves. We identify three comoving gyroresonances, which can lead to divergences in the Fokker-Planck coefficients. Taking into account thermal damping, we calculate the Fokker-Planck coefficient along with momentum diffusion coefficient D(p) and the mean-free path. Resulting acceleration time scales are compared with solar flare observations.

Description: We consider ion energization by Alfven waves, and suggest that a combination of nonlinear Landau damping and gyroresonant acceleration could account for proton acceleration in impulsive solar flares. Based on simulations which include both processes in a H plasma, we find that 10 exp -5 of the ambient protons can be accelerated above 1 MeV and that these particles contain about 10 percent of the total available wave energy. From a simulation in a H-He plasma, we find a large ratio of accelerated alpha particles to protons. Furthermore, arguments based upon the nonlinear Landau damping rate indicate heavy element enhancements. These results are consistent with accelerated particle observations from impulsive flares.