ALBERT

Description: We study the transfer of momentum from photons to dust grains to (molecular) gas in the outflow around cool giants (carbon-stars, Mira variables and OH/IR stars) beyond the radius where the dust grains condense. The problem is circular: radiation pressure determines the outflow velocity of the dust and thus also the dust density; on the other hand the dust density determines, via radiative transfer effects, the spectrum of the photons and thus the effective radiation pressure. This circular problem is solved by a rapidly converging iterative procedure. We compare our predictions with observed properties of a large sample of OH/IR stars and of Miras and find a good qualitative and quantitative agreement. We confirm a conclusion by Wood et al. (1993) that very luminous OH/IR stars in the Large Magellanic Cloud (LMC) owe their low outflow velocity to the low dust-to-gas ratio, a consequence of the low metallicity of the LMC. Similarly we consider a sample of about 100 OH/IR stars within 200 pc from the galactic center that has an average asymptotic giant branch (AGB) luminosity and an uncommonly high value of v(sub out); we conclude that these stars are probably very metal rich, perhaps even more than the stars in the Baade window studied by Rich (1990).

Description: Magnetic field data measured by the MAGMA instrument in the Martian magnetotail lobes are compared with the ram pressure of the upstream solar wind observed by the TAUS instrument in the circular orbits of the Phobos 2 spacecraft. High correlation was found between the magnetic field intensity in the Martian magnetotail lobes and the solar wind ram pressure. From this relationship the average flaring angle of the Martian magnetotail was determined as approximately 13 deg, and the average magnetosonic Mach number was estimated as approximately 5. The observed relationship between the Martian magnetotail magnetic field intensity and the solar wind magnetic field reflects the correlation of the solar wind magnetic field to the ram pressure providing a value of approximately 7 for the average Alfvenic Mach number. The flaring angle obtained for the Martian magnetotail was found to be an intermediate value between the flaring angle of the magnetotail of the Earth and that of Venus at comparable distances.

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: We present a model of the focused transport of approximately 1 MeV solar energetic protons through interplanetary Alfven waves that the protons themselves amplify or damp. It is based on the quasi-linear theory but with a phenomenological pitch angle diffusion coefficient in the 'resonance gap.' For initial Alfven wave distributions that give mean free paths greater than approximately 0.5 AU for approximately 1 MeV protons in the inner heliosphere, the model predicts greater than roughly an order of magnitude amplification (damping) in the outward (inward) propagating resonant Alfven waves at less than or approximately equal to o.3 AU heliocentric distance. As the strength of proton source is increased, the peak differential proton intensity at approximately 1 MeV at 1 AU increases to a maximum of approximately 250 particles (/(sq cm)(s)(sr)(MeV)) and then decreases slowly. It may be attenuated by a factor of 5 or more relative to the case without wave evolution, provided that the proton source is sufficiently intense that the resulting peak differential intensity of approximately 1 MeV protons at 1 AU exceeds approximately 200 particles (/(sq cm)(s)(sr)(MeV)). Therefore, in large solar proton events, (1) one may have to take into account self-amplified waves in studying solar particle propagation, (2) the number of accelerated protons escaping from a flare or interplanetary shock may have been underestimated in past studies by a significant factor, and (3) accelerated protons escaping from a traveling interplanetary shock at r less than or approximately equal to 0.3 AU should amplify the ambient hydromagnetic waves siginificantly to make the shock an efficient accelerator, even if initially the mean free path is greater than or approximately equal to 1 AU.

Description: Direct one-photon annihilation rate of positrons with a bound atomic electron is evaluated in the nonrelativistic limit. The K- and L-shell contributions are estimated including the screening and effective Coulomb repulsion effects. The annihilation rate of thermal positrons is calculated for various temperatures. The total number of one-photon annihilation events in the interstellar medium is discussed. These results provide the directional and structural information for cosmic gamma-ray sources.

Description: The origin of tektites has been obscure because of the following dilemma. The application of physical principles to the data available on tektites points strongly to origin from one or more lunar volcanoes; but few glasses of tektite composition have hitherto been reported from the lunar samples. Instead, the lunar silicic glasses consist chiefly of a material very rich in K2O and poor in MgO. The ratio of K2O/MgO is higher in these glasses than in any tektites reported. The solution of the dilemma seems to come from the study of some recently discovered terrestrial deposits of tektite glass with high values of K2O/MgO at the Cretaceous Tertiary boundary. These glasses are found to be very vulnerable to crystallization into sandine or to alteration to smectite. These end products are known and are more abundant than any terrestrial deposits of tektite glass. It seems possible that, in fact, the moon produces tektite glass, mostly of the high K2O-low MgO type; but on Earth these deposits are destroyed. The much less abundant deposits with lower K and higher Mg are observed because they survive. Other objections to the lunar origin hypothesis appear to be answerable.

Description: We investigate the radiative shock overstability for finite-sized objects. We follow the analysis of Chevalier & Imamura (1982), but we take into account the transverse flow of material out of the potshock region. The mass loss from the postshock region stabilizes the flow. As a rough estimate, the shock radiative instability takes place when the shock wave position with no radiative cooling (only mass loss present) is larger than the shock position with no mass loss (only radiative cooling present). For typical conditions of planetary nebulae we find that in order for the shock radiative overstability to occur, the nebular radius should be R approximately less than 10(exp 19) n(sub a)(exp -1) cm, where n(sub alpha) is the total number density of the interstellar medium (in units of cm(exp -3). We give several examples of interacting planetary nebulae in light of this condition.

Description: The energy spectra of the black-hole candidate GX 339-4 in the low-intensity state were observed on four occasions through 1989 to 1991 with the Large Area Counter on board the Ginga satellite. The spectra showed significant deviations from a power-law, with an iron K(sub alpha) emission line at approximetaly 6.4 keV and a broad iron K-edge structure above approximately 7 keV. The enrgy spectra above approximately 4 keV were successfully explained with a reflection model, in which part of the incident X-rays with a power-law spectrum is Compton reflected by optically thick matter, resulting in a harder continuum component with iron K-edge absorption and an iron flourescent line. The line equivalent width with respect to the reflection component decreases as the source flux increases. This is consistent with an increase in the ionization state of the material, so that resonant absorption followed by Auger ionization depletes the line. The photon-index of the power-law component was clearly variable, and it correlated with the relative amount of the reflection component. Such a correlation may be explained in the context of the anisotropic Comptonization models of Haardt et al. (1993), or by a variation of the relative geometry of the source and disk.

Description: We suggest that prior to its impact with Jupiter, comet Shoemaker-Levy 9 will behave as an electrical generator in the Jovian magnetosphere, converting planetary rotational energy to electrical energy via a dust/plasma interaction. This electrical energy will then be deposited in the dayside auroral region where it may drive various auroral phenomena including cyclotron radio emission. Such emission could be detected by spacecraft like Ulysses and Galileo many hours prior to the actual comet impact with the upper atmosphere. We apply the theory originally developed to explain the spokes in Saturn's rings. This theory allows us to quantify the driving potential associated with the comet and, consequently, to determine the radio power created in the auroral region. We conclude that if enough fine dust is present in the cometary system, comet-induced auroral radio emissions will reach detectable levels. This emission should be observable in the dayside hemisphere about 12-24 hours prior to each fragment impact.

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 Far-Infrared Absolute Spectrophotometer (FIRAS) instrument on the Cosmic Background Explorer (COBE) has determined the dipole spectrum of the cosmic microwave background radiation (CMBR) from 2 to 20/cm. For each frequency the signal is decomposed by fitting to a monopole, a dipole, and a Galactic template for approximately 60% of the sky. The overall dipole spectrum fits the derivative of a Planck function with an amplitude of 3.343 +/- 0.016 mK (95% confidence level), a temperature of 2.714 +/- 0.022 K (95% confidence level), and an rms deviation of 6 x 10(exp -9) ergs/sq cm/s/sr cm limited by a detector and cosmic-ray noise. The monopole temperature is consistent with that determined by direct measurement in the accompanying article by Mather et al.

Description: Theoretical O I density-sensitive emission-line ratios R = I(2s(sup 2))(2p(sup 4))((sup 3)P(sub 0))-((2s(sup 2))(2p(sup 4))((sup 3)P(sub 1)))/I((2s(sup 2))(2p(sup 4))((sup 3)P(sub 1))-(2s(sup 2))(2p(sup 4))((sup 3)P(sub 2))) = I(146 micrometers)/I(63 micrometers) are presented for a range of temperatures (T = 100-10,000 K), neutral hydrogen densities (N(sub H) = 10(exp -2) to 10(exp 7)/cu cm) and radiation fields (G(sub 0) = 1-10(exp 6)) applicable to both photodissociation regions (PDRs) and H II regions and the diffuse ionized medium (DIM). The observed values of R for several PDRs, measured from far-infrared spectra obtained with the Kuiper Airborne Observatory (KAO), imply hydrogen densities which are in good agreement with those determined using other methods. This provides observational support for the validity of the theoretical O I line ratios, and hence the atomic data used in their derivation.

Description: The 1.6-1.8 micron spectrum of the planetary nebula, IRAS 21282+5050, a strong emitter of the unidentified interstellar bands, contains a 0.02 micron wide eimission feature centered at 1.680 micron, which is well matched by laboratory spectra of the 0-2 CH stretching mode in polycyclic aromatic hydrocarbons (PAHs). We identify the new feature as the overtone of the well-known 3.3 micron band. In view of the high excitation required for emission in this band, the identification indicates that the emission is by free molecules rather than molecular moieties in solid dust grains. Modeling of the intensity ratio of the 2-0 to 1-0 band implied that the PAHs emitting in these bands contain about 60 carbon atoms. It is inferred that the nu = 2-1 hot band of the CH stretching mode occurs at about 3.43 micron and contributes to the long-wavelength shoulder of the 3.40 micron feature. The main 3.40 micron feature probably is due to aliphatic sidegroups on PAH molecules.

Description: Weakly nonlinear Magneto Hydrodynamic (MHD) stability of the Halley cometosheath determined by the balance between the outward ion-neutral drag force and the inward Lorentz force is investigated including the transverse plasma motion as observed in the flanks with the help of the method of multiple scales. The eigenvalues and the eigenfunctions are obtained for the linear problem and the time evolution of the amplitude is obtained using the solvability condition for the solution of the second order problem. The diamagnetic cavity boundary and the adjacent layer of about 100 km thickness is found unstable for the travelling waves of certain wave numbers. Halley ionopause has been observed to have strong ripples with a wavelength of several hundred kilometers. It is found that nonlinear effects have stabilizing effect.

Description: We have developed a radiative transfer model of the dust and gas envelopes around late-type stars. The gas kinetic temperature for each star is calculated by solving equations of motion and the energy balance simultaneously. The main processes include viscous heating and adiabatic and radiative cooling. Heating is dominated by viscosity as the grains stream outward through the gas, with some contribution in oxygen-rich stars by near-IR pumping of H2O followed by collisional de-excitation in the inner envelope. For O-rich stars, rotational H2O cooling is a dominant mechanism in the middle part of the envelope, with CO cooling being less significant. We have applied our model to three well-studied oxygen-rich red giant stars. The three stars cover a wide range of mass-loss rates, and hence they have different temperature structures. The derived temperature structures are used in calculating CO line profiles for these objects. Comparison of the dust and gas mass-loss rates suggests that mass-loss rates are not constant during the asymptotic giant branch phase. In particular, the results show that the low CO 1-0 antenna temperatures of OH/IR stars reflect an earlier phase of much lower mass-loss rate.

Description: We have studied the spectacular 1991 June X-class flares using gamma-ray data from the Charged Particle Detectors (CPDs) of the Burst and Transient Source Experiment (BATSE) on the Compton Gamma Ray Observatory (CGRO) and 80 GHz millimeter data from Nobeyama, Japan. The CPDs were the only CGRO instrument that did not saturate during the extremely intense 1991 June 4 flare. We have shown that for this flare the CPDs respond to MeV photons, most of which are due to bremsstrahlung produced by relativistic electrons at the Sun. We have further shown that the gamma-ray and millimeter observations agree numerically if the 80 GHz radiation is gyrosynchrotron radiation produced by trapped electrons and the gamma rays are thick-target bremsstrahlung due to electrons precipitating out of the trap. The requirement that the trapping time obtained from the numerical comparison be consistent with the observed time profiles implies a magnetic field between about 200 and 300 G and an electron spectral index between about 3 to 5. By comparing the CPD observations with both the 80 GHz data and nuclear line data from the Energetic Gamma Ray Experiment Telescope (EGRET) and the Oriented Scintillation Spectroscopy Experiment (OSSE) on CGRO for the flares of June 4, 6, 9, and 11, we found that the ratio of the CPD counts to both the millimeter flux densities and the nuclear line fluences decreases with decreasing flare heliocentric angle. All of these flares were produced in the same active region. We interpreted this result in terms of a loop model in which the gyrosynchrotron emission is produced in the coronal portion of the loop where the electrons are kept isotropic by pitch angle scattering due to plasma turbulence, while the bremsstrahlung is produced by precipitating electrons that interact anisotropically. We found that the trapping time in the coronal portion is time dependent, reaching a minimum of about 10 s at the peak of the CPD count rate. We suggested the damping of the turbulence as a possible reason for the variation of the trapping time. turbulence as a possible reason for the variation of the trapping time.

Description: Several possible models have been suggested to explain the observed distribution of gamma-ray bursts: heliocentric distributions such as the Oort cloud, large galactic halos, and cosmological models. We report here on an investigation into the implications of the Burst and Transient Source Experiment (BATSE) gamma-ray burst distribution (Meegan et al. 1992a) data on the possible helocentric origin of gamma-ray bursts. We find no statistically significant anisotropy in the angular distribution of the bursts in a Sun-referenced coordinate system; there is no dipole moment in the direction of the Sun, and no quardrupole moment associated with the ecliptic plane. We have employed direct analytic calculations and Monte Carlo simulations of sources in the Oort cloud to constrain possible helicentric burst distributions. These can produce distributions consistent with the observed angular isotropy, the meal value of V/V(sub max), and the observed C/C(sub min) distribution of BATSE, and provide limits to burst energy of a few times approximately 10(exp 27) ergs. However, the agreement of the heliocentric C/C(sub min) distributions with the BATSE data is attributable to the relatively limited sampling of strong, nearby bursts. These bursts are known from observation to be homogeneously distributed, yet the density of sources in the Oort cloud is not constant in this region. Integral number-intensity distributions from the Oort cloud for larger numbers of bursts cannot reproduce the known homogeneity of the strong bursts without modification to the computed cometary number density and are therefore unlikely explanations of the gamma-ray burst distribution.

Description: We report the first X-ray observations of the Seyfert 1 galaxy MCG-6-30-15 obtained at medium spectral resolution. The partially-ionized, 'warm' absorber is resolved and shown to be due to O VII and O VIII. The main absorption edge agrees with that of O VII at the redshift of the galaxy to within 1%. The column density of the absorbing material is greater by a factor of 2 in the first of our two obsevations, which were 3 weeks apart, while the mean flux is slightly lower and the ionization parameter slightly higher. We also discuss the flourescent iron emssion line seen in the source, which is at 6.40 keV. The line is significantly broadened, with a Full Width at Half Maximum (FWHM) of about 0.4 keV.

Description: Near supernova 1987a, the rare honeycomb structure of 20-30 galactic bubbles measures 30 x 90 light years. Its remarkable regularity in bubble size suggests a single-event origin which may correlate with the nearby supernova. To test the honeycomb's regularity in shape and size, the formalism of statistical crystallography is developed here for bubble sideness. The standard size-shape relations (Lewis's law, Desch's law, and Aboav-Weaire's law) govern area, perimeter and nearest neighbor shapes. Taken together, they predict a highly non-equilibrium structure for the galactic honeycomb which evolves as a bimodal shape distribution without dominant bubble perimeter energy.

Description: If, as many believe, Sgr A* is a massive black hole at the Galactic center, one should expect it to be a source of X-ray and gamma-ray activity, behaving basically as a scaled-down active galactic nucleus. An unavoidable source of accretion is the wind from IRS 16, a nearby group of hot, massive stars. Since the density and velocity of the accreting matter are known from observations, the accretion rate is basically a function of the putative black hole mass, M(sub h), only; this value represents a reliable lower limit to a real rate, given the other possible sources of accreting matter. Based on this and on the theories about shock acceleration in active galactic nuclei, we have estimated the expected production of relativistic particles and their hard radiation. These values turn out to be a function of M(sub h) as well. Comparing our results with available X-ray and gamma-ray observations which show Sgr A* to have a relatively low activity level, we conclude tentatively that the putative black hole in the Galactic center cannot have a mass greater than approximately 6 x 10(exp 3) solar mass. This conclusion is consistent with the upper limits to the black hole mass found by different methods earlier, although much more work is needed to make calculations of shock acceleration around black holes more reliable.

Description: The light curve of the Type Ia supernova SN 1937C (in IC 4182) is important because Sandage et al. have measured a distance to the host galaxy by means of Cepheid variables and thus have derived the Hubble constant. However, the peak brightness of SN 1937C has only been derived with the relatively poor original comparison star brightnesses and without regard to a large body of data in the literature. In this paper, I will correct these and other procedural difficulties. I find that the late time photographic light curve appears to have a broken exponential decay with equivalent half-lives of 46 and 58 days with the break near 300 days after maximum. I also find that the peak B-magnitude was 8.71 +/- 0.14 on JD 2428770.0 +/- 1.0 at which time the B-V was -0.03 +/- 0.13. With these improved peak brightnesses, the distance modulus of Sandage et al., and peak absolute magnitudes in the center of the range of modern estimates, I derive the Hubble constant to be 50 km/s Mpc.

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: We have found a new way to make Thorne-Zytkow objects, which are massive stars with degenerate neutron cores. The asymmetric kick given to the neutron star formed when the primary of a massive tight binary system explodes as a supernova sometimes has the appropriate direction and amplitude to place the newly formed neutron star into a bound orbit with a pericenter distance smaller than the radius of the secondary. Consequently, the neutron star becomes embedded in the secondary. Thorne-Zytkow objects are expected to look like extreme M-type supergiants, assuming that they can avoid a runaway neutrino instability. Accretion onto the embedded neutron star will produce either an isolated, spun-up neutron star (possibly a short-period pulsar) or a black hole. Whether neutron star or black hole remnants predominate depends on the lifetime of Thorne-Zytkow objects, the accretion rates involved, and the maximum neutron star mass, none of which are definitively understood.

Description: Observations of the red supergiant (M2 Iab) alpha Ori with the Goddard High Resolution Spectrograph (GHRS) on board the Hubble Space Telescope (HST) have provided an unambiguous detection of a far-ultraviolet (far-UV) chromospheric continuum on which are superposed strong molecular absorption bands. The absorption bands have been identified by Carpenter et al. (1994) with the fourth-positive A-X system of CO and are likely formed in the circumstellar shell. Comparison of these GHRS data with archival International Ultraviolet Explorer (IUE) spectra of alpha Ori indicates that both the continuum and the CO absorption features can be seen with IUE, especially if multiple IUE spectra, reduced with the post-1981 IUESIPS extraction procedure (i.e., with an oversampling slit), are carefully coadded to increase the signal to noise over that obtainable with a single spectrum. We therefore initiated a program, utilizing both new and archival IUE Short Wavelength Prime (SWP) spectra, to survey 15 cool, low-gravity stars, including alpha Ori, for the presence of these two new chromospheric and circumstellar shell diagnostics. We establish positive detections of far-UV stellar continua, well above estimated IUE in-order scattered light levels, in spectra of all of the program stars. However, well-defined CO absorption features are seen only in the alpha Ori spectra, even though spectra of most of the program stars have sufficient signal to noise to allow the dectection of features of comparable magnitude to the absorptions seen in alpha Ori. Clearly if CO is present in the circumstellar environments of any of these stars, it is at much lower column densities.

Description: The giant radio galaxy NGC 6251 is a particularly good object for observational tests of relativistic jet models. Due to its high declination and approximately 0.5 Jy radio nucleus, high-quality Very Long Baseline Interferometry (VLBI) images of the central regions of the source can be made with northern hemisphere arrays. In addition, the large-scale radio morphology strongly suggests that the radio axis lies close to the plane of the sky, so Doppler boosting should be less extreme than in the core-dominated superluminal sources. Earlier 18 cm VLBI observations of NGC 6251 revealed an unexpectedly large jet/counterjet brightness ratio and small transverse motion of a feature in the parsec-scale jet. These early results are difficult to reconcile with the simplest symmetric relativistic jet models. In this paper we present a third-epoch 18 cm VLBI image of the parsec-scale radio jet in NGC 6251, and compare jet morphology over a 5 year time span. The jet shows a minor brightness peak at nearly the same distance from the core as the '25 mas knot' seen in the first- and second-epoch VLBI images. This feature is much less pronounced in the third epoch, and a relatively bright, new knot has appeared approximately 12 mas from the core. If this new component had a constant brightness during the 5 years separating the first and third observing epochs, then it must have moved away from the core with an apparent speed of at least 1.2c (compared with an upper limit of 0.23c for motion of the 25 mas knot). However, we cannot yet rule out a local brightening of the inner jet in favor of a new moving component. We determine a lower limit for the jet/couterjet brightness ratio of 100:1 within 6 mas of the core. We also present a new Very Large Array (VLA) image of the kpc-scale jet with 3 sec resolution, made from data obtained during the VLBI observations. The rate of decrease in jet surface brightness from parsec to kiloparsec scales is similar to jets in known superluminal radio sources.

Description: We compute the three-point temperature correlation function of the Cosmic Background Explorer (COBE) Differential Microwave Radiometer (DMR) first-year sky maps to search for non-Gaussian temperature fluctuations. The level of fluctuations seen in the computed correlation function are too large to be attributable solely to instrument noise. However the fluctuations are consistent with the level expected to result from a superposition of istrument noise and sky signal arising from a Gaussian power-law model of initial fluctuations, with a quadrupole normalized amplitude of 17 micro K and a power-law spectral index n = 1. We place limits on the amplitude of intrinsic three-point correlations with a variety of predicted functional forms.

Description: The Burst and Transient Source Experiment on the Compton Gamma-Ray Observatory detected 260 cosmic gamma-ray bursts during the period 19 Apr 1991 to 5 Mar 1992. This paper presents the occurrence times, locations, peak count rates, peak fluxes, fluences, durations, and plots of time histories for these bursts. The angular distribution is consistent with isotropy. The intensity distribution shows a deficit in the number of weak bursts, which is not consistent with a homogeneous distribution of burst sources in Euclidean space. The duration distribution shows evidence for a separate class of bursts with durations less than about 2 seconds.

Description: A simple method for deriving well-behaved temperature solutions to the equation of hydrostatic equilibrium for intracluster media with X-ray imaging observations is presented and applied to a series of generalized models as well as to observations of the Perseus cluster and Abell 2256. In these applications the allowed range in the ratio of nonbaryons to baryons as a function of radius is derived, taking into account the uncertainties and crude spatial resolution of the X-ray spectra and considering a range of physically reasonable mass models with various scale heights. Particular attention is paid to the central regions of the cluster, and it is found that the dark matter can be sufficiently concentrated to be consistent with the high central mass surface densities for moderate-redshift clusters from their gravitational lensing properties.

Description: Electric dipole transition matrix elements for rovibrational transitions in the X (sup 1)Sigma(sup +) state of the CO minor isotopes (14)C(16)O and (13)C(17)O are calculated for the first time for all the delta v = +1, +2, and +3 transitions for which v less than or equal to 20 and J less than or equal to 150. Improved electric dipole transition matrix elements are also calculated for the minor isotopes (12)C(17)O, (12)C(18)O, (13)C(18)O. We have fitted polynomials to these matrix elements as a function of the parameter m which is defined in terms of the lower state angular momentum quantum number J; the convenient to use polynomial representations are given in tabular form. These results for the minor species of CO complement those previously reported by us for (12)C(16)O and (13)C(16)O.

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: The observational selection bias properties of the large Mathewson-Ford-Buchhorn (MFB) sample of axies are demonstrated by showing that the apparent Hubble constant incorrectly increases outward when determined using Tully-Fisher (TF) photometric distances that are uncorreted for bias. It is further shown that the value of H(sub 0) so determined is also multivlaued at a given redshift when it is calculated by the TF method using galaxies with differenct line widths. The method of removing this unphysical contradiction is developed following the model of the bias set out in Paper II. The model developed further here shows that the appropriate TF magnitude of a galaxy that is drawn from a flux-limited catalog not only is a function of line width but, even in the most idealistic cases, requires a triple-entry correction depending on line width, apparent magnitude, and catalog limit. Using the distance-limited subset of the data, it is shown that the mean intrinsic dispersion of a bias-free TF relation is high. The dispersion depends on line width, decreasing from sigma(M) = 0.7 mag for galaxies with rotational velocities less than 100 km s(exp-1) to sigma(M) = 0.4 mag for galaxies with rotational velocities greater than 250 km s(exp-1). These dispersions are so large that the random errors of the bias-free TF distances are too gross to detect any peculiar motions of individual galaxies, but taken together the data show again the offset of 500 km s(exp-1) fond both by Dressler & Faber and by MFB for galaxies in the direction of the putative Great Attractor but described now in a different way. The maximum amplitude of the bulk streaming motion at the Local Group is approximately 500 km s(exp-1) but the perturbation dies out, approaching the Machian frame defined by the CMB at a distance of approximately 80 Mpc (v is approximately 4000 km s(exp -1)). This decay to zero perturbation at v is approximately 4000 km s(exp -1) argues against existing models with a single attraction at approximately 4500 km s(exp -1) (the Great Attactor model) pulling the local region. Rather, the cause of the perturbation appears to be the well-known clumpy mass distribution within 4000 km s(exp -1) in the busy directions of Hydra, Centaurus, Antila and Dorado, as postulated earlier (Tammann & Sandage 1985).

Description: Nonthermal radio emission has been observed from some of the most luminous hot star winds. It is understood to be synchrotron radiation of the relativistic electrons in the winds. To understand how the electrons are accelerated to such high energies and to correctly explain the observed radio flux and spectra require an exhaustive investigation of all the relevant physical processes involved and possibly point to a complex wind structure. In this paper we discuss the logical path toward a comprehensive model of the nonthermal radio emission from hot star winds. Based on the available observational data and fundamental theoretical considerations, we found that the only physically viable and self-consistent scenario is: the nonthermal radio emission is synchrotron radiation of relativistic electrons the electrons are accelerated by shocks via the first-order Fermi mechanism the acceleration has to be in situ in the radio emitting region and the shocks formed at the base of the winds have to propagate to beyond the radio photosphere.

Description: Extensive observations of the molecular gas in the young, compact planetary nebula M1-16 have been made, using the Swedish-ESO-Submillimeter Telescope. A map of the CO J = 2-1 emission shows that the molecular envelope contains both a slow and a fast outflow with expansion velocities of 19 km/s and greater than 34 km/s, respectively. The slow outflow is mildly elliptical, while the fast molecular outflow is bipolar. This fast outflow is roughly aligned with the very fast outflows recently found in the optical, while the long axis of the slow elliptical outflow is roughly orthogonal to the optical outflow axis. The kinematic timescales for the CO fast outflow and the optical very fast outflow agree closely, supporting the view that the former represents material in the slow outflow accelerated by the very fast outflow. The kinematic signature of a disk expanding with about 15.5 km/s can also be seen in the CO J = 2-1 data. The mass-loss rate (a) for the slow outflow is greater than or equal to 2.8 x 10(exp -5) solar mass/yr and possibly as large as 9 x 10(exp -5) solar mass/yr, (b) for the fast outflow is greater than or equal to 5 x 10(exp -6) solar mass/yr, and (c) for the very fast optically visible outflow is approximately equal 5 x 10(exp -7) solar mass/yr. The disk mass is approximately equal 6 x 10(exp -3) solar mass. Grain photoelectric heating results in temperatures of 20-70 K in molecular gas of the slow outflow. The (13)C/(12)C abundance ratio in M1-16 is found to be 0.33, quite possibly the highest found for any evolved object. Upper limits for the (18)O/(16)O and (17)O/(16)O ratios were found to be consistent with the values found in AGB stars. A search for other molecular species in M1-16 resulted in the detection of the high-excitation species HCN, CN, (13)CN, HCO(+), and H(13)CO(+) and possibly N2H(+). Both the HCO(+)/HCN and CN/HCN line-intensity ratios are enhanced, the former by a very large factor, over the values found in the envelopes of AGB stars, probably as a result of enhancement of the CN and HCO(+) abundances due to photochemistry induced by the stellar UV. The CS J = 2-1, SiO J = 2-1 (v = 0), and SiS J = 6-5 lines were not detected to low levels. For the high-excitation molecules, adequate collisional excitation of rotational levels and survival against photodissociation by the UV radiation requires significant clumping of the molecular gas into clumps with H2 densities approximately 10(exp 5)/cu cm. The IRAS fluxes of M1-16, assuming negligible contribution from line emission, imply the presence of about (1.7-0.4) x 10(exp -3) solar mass of cool dust (temperature around 50 K) and a smaller quantity, (2.7-3.1) x 10(exp -6) solar mass, of warmer dust (temperature around 125 K) for a power-law emissivity index p = 1-2. The evolutionary nature of M1-16 cannot be explained by existing single-star models of post-AGB evolution. The very high (13)C/(12)C abundance ratio in M1-16 suggests a possible evolutionary connection between M1-16 and the rare class of J-type silicate-carbon stars which also have high (13)C/(12)C ratios and are thought to be binary systems with accretion disks.

Description: We have observed the C91 alpha radio recombination line toward the Orion H II region. This narrow (approximately 3-5 km per sec full width at half maximum (FWHM)) line is spatially very extended (approximately 8 arcmin or 1 pc). These charateristics compare well with the observed characteristics of the C II fine structure line at 158 microns. Thus, the C91 alpha line originates in the predominantly neutral photodissociation regions separating the H II region from the molecular cloud. We have developed theoretical models for the C II radio recombination lines from photodissociation regions. The results show that the I(C91 alpha)/I(C158) intensity ratio is a sensitive function of the temperature and density of the emitting gas. We have also extended theoretical models for photodissociation regions to include the C II recombination lines. Comparison with these models show that, in the central portion of the Orion region, the C91 alpha line originates in dense (10(exp 6) per cu cm), warm (500-1000 K) gas. Even at large projected distances (approximately 1 pc), the inferred density is still high (10(exp 5) per cu cm) and implies extremely high thermal pressures. As in the case of the (C II) 158 microns line, the large extent of the C91 alpha line shows that (FUV) photons can penetrate to large distances from the illuminating source. The decline of the intensity of the incident radiation field with distance from Theta(sup 1) C seems to be dominated by geometrical dilution, rather than dust extinction. Finally, we have used our models to calculate the intensity of the 9850 A recombination line of C II. The physical conditions inferred from this line are in good agreement with those determined from the radio recombination and the far-infrared fine-structure lines. We show that the ratio of the 9850 A to the C91 alpha lines is a very good probe of very high density clumps.

Description: Statistical uncertainties in determining the temperatures of hot (0.5-10 keV) coronal plasmas are investigated. The statistical presicion of various spectral temperature diagnostics is established by analyzing synthetic ASCA solid-state imaging spectrometer (SIS) CCD spectra. The diagnostics considered are the ratio of hydrogen-like to helium-like line complexes of Z greater than or = 14 elements, line-free portions of the continuum, and the entire spectrum. While fits to the entire spectrum yield the highest statistical precision, it is argued that fits to the line-free continuum are less susceptible to atomic data uncertainties but lead to a modest increase in statistical uncertainty over full spectral fits. Temperatures deduced from line ratios can have similar accuracy, but only over a narrow range of temperatures. Convenient estimates of statistical accuracies for the various temperature diagnostics are provided which may be used in planning ASCA SIS observations.

Description: The shape of the velocity distribution of water group ions observed by the Giotto ion mass spectrometer on its approach to comet Halley is modeled to derive empirical values for the rates of ionization, energy diffusion, and loss in the midcometosheath. The model includes the effect of rapid pitch angle scattering into a bispherical shell distribution as well as the effect of the magnetization of the plasma on the charge exchange loss rate. It is found that the average rate of ionization of cometary neutrals in this region of the cometosheath appears to be of the order of a factor 3 faster than the `standard' rates approx. 1 x 10(exp -6)/s that are generally assumed to model the observations in most regions of the comet environment. For the region of the coma studied in the present work (approx. 1 - 2 x 10(exp 5) km from the nucleus), the inferred energy diffusion coefficient is D(sub 0) approx. equals 0.0002 to 0.0005 sq km/cu s, which is generally lower than values used in other models. The empirically obtained loss rate appears to be about an order of magnitude greater than can be explained by charge exchange with the `standard' cross section of approx. 2 x 10(exp -15)sq cm. However such cross sections are not well known and for water group ion/water group neutral interactions, rates as high as 8 x 10(exp -15) sq cm have previously been suggested in the literature. Assuming the entire loss rate is due to charge exchange yields a rate of creation of fast neutral atoms of the order of approx. 10(exp -4)/s or higher, depending on the level of velocity diffusion. The fast neutrals may, in turn, be partly responsible for the higher-than-expected ionization rate.

Description: HH 80-81 are two optically visible Herbig-Haro (HH) objects located about 5 minutes south of their exciting source IRAS 18162-2048. Displaced symmetrically to the north of this luminous IRAS source, a possible HH counterpart was recently detected as a radio continuum source with the very large array (VLA). This radio source, HH 80 North, has been proposed to be a member of the Herbig-Haro class since its centimeter flux density, angular size, spectral index, and morphology are all similar to those of HH 80. However, no object has been detected at optical wavelengths at the position of HH 80 North, possibly because of high extinction, and the confirmation of the radio continuum source as an HH object has not been possible. In the prototypical Herbig-Haro objects HH 1 and 2, ammonia emission has been detected downstream of the flow in both objects. This detection has been intepreted as a result of an enhancement in the ammonia emission produced by the radiation field of the shock associated with the HH object. In this Letter we report the detection of the (1,1) and (2,2) inversion transitions of ammonia downstream HH 80 North. This detection gives strong suppport to the interpretation of HH 80 North as a heavily obscured HH object. In addition, we suggest that ammonia emission may be a tracer of embedded Herbig-Haro objects in other regions of star formation. A 60 micrometer IRAS source could be associated with HH 80 North and with the ammonia condensation. A tentative explanation for the far-infrared emission as arising in dust heated by their optical and UV radiation of the HH object is presented.

Description: The Giotto ion mass spectrometer high-intensity spectrometer (IMS-HIS) measured fluxes of ions from about 260,000 km before (1008:37 UT spacecraft time) to about 86,000 km after (1701:33 UT spacecraft time) closest approach to comet P/Grigg-Skjellerup during the encounter on July 10, 1992. Although the HIS sensor was not designed to measure protons, these ions were measured far from the comet. Close in to the comet, the ions observed were probably also protons, although heavier ions cannot be completely ruled out. Considerable temporal structure appears in the data, well-correlated with the data of other instruments onboard, especially those of the magnetometer. In particular, the ion count rate correlates with the direction of the magnetic field. This strong modulation at the water group ion cyclotron period (approx. 90 s) inside the inbound bow wave indicates a very narrow ion pitch angle distribution. Hence at Grigg-Skjellerup the ions appear to experience very little pitch angle scattering. This may result from strong compression in the rapidly increasing magnetic field.

Description: Measurements of the very broadband structure in the visual part of the extinction curve are compared to IUE extinction curve parameters in the scheme of Fitzpatrick and Massa. The correlation of Very-broadband-structure (VBS) depth with FUV rise found by Reimann & Friedemann (1991) is shown to be a correlation with linear rise and not with far UV non-linear rise. The correlation with linear rise suggests that the VBS is due to an onset of extinction at about 1.8/micrometer, rather than what has previously been suggested: due to luminescence or the presence of two extinction components longward and shortward of the VBS. The optical properties inferred for the linear rise carrier are consistent with some amorphous solid with a large optical gap. Small carbonaceous grains may be the carrier of the linear rise in extinction, because the erosion of core-mantle particles is expected to produce many of such small grains and offers a natural mechanism for the existence of the inferred two populations of big (a approximately 0.13 micrometer) and small (a approximately 0.005 micrometer) grains.

Description: The modulated light in the Intermediate Polar FO Aqr at the three periods P(sub spin) = 20.9 min, P(sub orb) = 4.85 hr and P(sub beat) = 22.5 min is studied in different spectral ranges to derive information on their nature. In this system the accretion geometry, with or without an accretion disk, is still a matter of debate (Hellier 1991; Norton et al. 1992). The different orbital behavior of phase coherence between the spin and beat pulses in the X-rays (Norton et al. 1992) and in the optical/IR regions cannot be easily accounted for by only a diskless dominated geometry where the accretion flow is switching from one pole to the other each half of the beat period. We therefore propose an accretion scenario where a non-axisymmetric disk is present. In such a non-standard accretion disk an azimuthal structure provides not only the source of variable mass transfer to the white dwarf, but also a reprocessing site which is mainly viewing the X-ray emission from the lower pole. Our spectral analysis shows that reprocessing is also occurring at the surface of the secondary star. The spin pulsation in the optical and IR continua can be explained by the so-called 'accretion curtain' model (Rosen et al. 1988) though an additional reprocessing component at the spin period cannot be excluded. In contrast to the X-rays, the beat optical/IR modulation is not intrinsic. Reprocessing at the surface of the secondary star and at the thickened part of the disk can also account for the orbital modulation in the UV, optical and IR regions.

Description: We present Vela 5B observations of the Small Magellanic Cloud (SMC) region obtained during the years 1969 - 1979. We detect the 3.89 day orbital modulation of SMC X-1 with a false alarm probability of 4.3 x 10(exp -3) and observe a high state beginning in 1970 September and lasting approximately 100 days. We also detect three outbursts by a transient source consistent with the position of the candidate Be-neutron star system H0107-750 (= 1H 0103-762). These events occur roughly 100 days apart and last for approximately 35 days. No detections of SMC X-2 or SMC X-3 are apparent above a limit of approximately 7 x 10(exp 37) ergs/s.

Description: The hard X-ray transient, GRO J1008-57, was discovered during a bright outburst in 1993 July by the Burst and Transient Source Experiment (BATSE) instrument on board the Compton Gamme-Ray Observatory (CGRO). There are no published reports of previous emission from this 93.6 sec X-ray pulsar. Recent optical results have suggested a Be star as the companion. A search of the EXOSAT archives shows that an ME observation centered on the star HD 88661 includes GRO J1008-57 within the field of view. The characteristics of the medium-energy detection including a hard spectrum and pulsed emission at 91.36 sec (with chance probability of 3 x 10(exp -4)) indicates that EXOSAT detected GRO J1008-57 rather than the field star. The estimated flux is 2.4 x 10(exp -11) ergs/cm sq/sec for a luminosity of 1.1 x 10(exp 34) (D/2 kpc)(exp 2) ergs/sec. The X-ray spectrum is hard (photon index approximately equals 1.2) and highly absorbed (N(sub H) approximately equals 0.7 x 10(exp 22)). The detection of this transient suggests that there is a pulse shape dependence on luminosity, a spectral shape independent of luminosity, a large period derivative of P-dot/P = 2.3 x 10(exp -3)/yr, and a dynamic range of at least 100 in L(sub chi). We discuss GRO J1008-57 data in the context of wind-accretion models for this new member of the Be class of X-ray binaries.

Description: In 1991 and 1992, the dust detector onboard the Ulysses spacecraft detected several dust streams apparently originating from the jovian system. The timing and measured speeds of the final two dust streams are compatible with dust from comet Shoemaker-Levy 9's (SL9) disruption in 1992. Our further investigations of stream characteristics and dust acceleration mechanisms, however, shed some doubt that two of the eleven dust streams are of SL9 origin. In July 1994 when SL9 impacts Jupiter, the Galileo spacecraft will be about 3500 jovian radii away from the planet. Submicronsized dust released into, and accelerated by, the jovian magnetosphere during this event may reach Galileo and impact its dust detector between September and November 1994. We also discuss the possibility of directly sampling dust from SL9 during Galileo's orbital tour.

Description: We have investigated the influence of the color temperature of the illuminating radiation field on the chemical and thermal structure of photon-dominated regions (PDRs). We present the results of a study of the photoelectric efficiency of heating by large molecules such as polycyclic aromatic hydrocarbons (PAHs) and very small grains for radiation fields characterized by different effective temperatures. We show that the efficiency for cooler (T(sub eff) approximately = 6000-10,000 K) stars is at most an order of magnitude smaller than that for hotter (T(sub eff) approximately = 20,000-30,000 K) stars. While cooler radiation fields result in less ultraviolet photons capable of heating, the efficiency per absorbed photon is higher, because the grains become less positively charged. We also present detailed calculations of the chemistry and thermal balance for generic PDRs (n(sub 0) approximately = 10(exp 3), G(sub 0) approximately = 10(exp 3)). For cooler radiation fields, the H/H2 and C(+)/C/CO transition layers shift toward the surface of the PDR, because fewer photons are available to photodissociate H2 and CO and to ionize C. The dominant cooling lines are the (C II) 158 micron and the (O I) 63 micron lines for the hotter radiation fields, but cooling by CO becomes dominant for a color temperature of 6000 K or lower. The (C II)/CO and (O I)/CO ratios are found to be very good diagnostics for the color temperature of the radiation field.

Description: The theory of spectral formation in thermal X-ray sources, where the effects of Comptonization and Klein-Nishina corrections are important, is presented. Analytical expressions are obtained for the produced spectrum as a function of such input parameters as the plasma temperature, the optical depth of the plasma cloud and the injected soft photon spectrum. The analytical theory developed here takes into account the dependence of the scattering opacity on the photon energy. It is shown that the plasma temperature as well as the asymptotic rate of photon escape from the plasma cloud determine the shape of the upscattered hard tail in the emergent spectra, even in the case of very small optical depths. The escape distributions of photons are given for any optical depth of the plasma cloud and their asymptotic dependence for very small and large optical depths are examined. It is shown that this new generalized approach can fit spectra for a large variety of hard X-ray sources and determine the plasma temperature in the region of main energy release in Cyg X-1 and the Seyfert galaxy NGC 4151.

Description: We have measured the spectra of H and He isotopes during the 1987 solar minimum with the cosmic-ray detector system (CRS) on the Voyager 2 spacecraft. By carrying out the measurement near solar minimum and at large heliospheric distances, the effects of solar modulations were reduced. In particular, the adiabatic energy losses were smaller, and these results from 23 AU over the solar minimum period of cycle 21 represent observations at energies not accessible from previous measurements near 1 AU. The modulated spectra with the diffusion coefficient constant k(sub 0) = 3.15 x 10(exp 22) sq cm/s (which corresponds to a solar modulation parameter of 360 MV at 23 AU and 500 MV at 1 AU) agree well with both our data at 23 AU and the previous solar minimum measurements at 1 AU. The measured H-1 and H-2 spectra are both consistent with the calculated spectra, using standard Galactic and heliospheric propagation models without invoking an anomalous hydrogen component. With the fixed modulation parameter of 360 MV, the mean pathlengths, source spectra, and cross sections were varied to study the effects of different input parameters on the spectra and relative abundances. At this stage of our work, we have not found any strong evidence from the low-energy H-2 and He-3 data that H-1 and He-4 should have a different propagation history, or different types of source spectra from the heavier cosmic-ray nuclei.

Description: Significant differences are found in the IRAS color-color diagrams of small regions (2 min x 2 min, or 0.4 x 1.8 kpc) within the disk of M31 compared to Galactic cirrus, most noticeably demonstrated by a trend of low 60 to 100 micrometer surface brightness ratio and high 12 to 25 micrometer ratio. Based on physical arguments, we conclude that these color differences are best explained by assuming that 'very small grains' (VSG; but not polycylic aromatic hydrocarbons) are only half as abundant in M31 as they are in Galactic cirrus. We confirm this conclusion and test its detailed agreement with data by using the phenomenological model of Desert et al. (1990). In particular, we show that the data cannot be explained by postulating weaker UV heating in the disk of M31. We also show that the VSG-deficient model predicts correctly the correspondence between the IRAS colors and the 100 micrometer emissivity per H I atom in the outer disk of M31. 'Very small grains' are a leading candidate for the carrier of the 2175 A bump in the extinction curve. Our suggested VSG deficiency in M31 is thus consistent with recent Hubble Space Telescope (HST) observations which show evidence for a weaker and narrower 2175 A bump on the M31 extinction curve. Some speculation is offered as to possible links between very small grains and the low rate of current star formation in M31.

Description: We present a model for explaining the recent combined X-ray and low-energy gamma-ray observations of the Seyfert galaxy NGC 4151. According to this model, soft photons become Comptonized in a hot spot producing simultaneously the low-energy power law as observed by Ginga and the high-energy cutoff observed by the Oriented Scintillation Spectrometer Experiment (OSSE). Implementing recently developed theoretical calculations toward a generalized theory of Comptonization, we were able to find fits to the observations using only two parameters which characterize the physical quantities of the emission region: the plasma cloud optical depth and its temperature. We find that there is no need for additional nonthermal, reflection, or higher temperature thermal components to fit the aforementioned OSSE and Ginga observations. We derive in addition the size of the photon region and the temperature of the upscattered soft photons. We should emphasize, also, that any attempt at fitting only the high-energy parts of the spectrum (photon energies greater than 60 keV) by the Sunyaev & Titarchuk (1980) nonrelativistic Comptonization model leads to an underestimate of the Comptonization parameter y (or, equivalently, to an overestimation of the X-ray power-law spectral slope) and leads, as a result, to incorrect proportions between the low-energy and high-energy parts of the spectrum.

Description: Observational data for Population I stars have shown that blue loops on the Hertzsprung-Russell (H-R) diagram form for stellar masses as low as approximately 4 solar mass. However, current state-of-the-art stellar models, unlike the older ones that were based on smaller opacities, fail to loop out of the red-giant region during core helium burning for masses less than 7 solar mass. A possible explanation is that the currently used Livermore opacities need to be further increased, by at least 70%, at temperatures characteristic of the base of the outer convection zone, around 1 x 10(exp 6) K. Indeed, no other suggested remedy seems to yield a blue loop at the lowest observed loop luminosities.

Description: The recent development of trapped-ion frequency standards, which offer high stability over long time periods, provides us with a potential new method for detecting unseen matter in the outer solar system. A distribution of matter or a planetary body could produce a measurable gravitational redshift of the radio signal received from a spacecraft equipped with an ultrastable frequency standard. Trapped-ion standards have a potential frequency stability of 1 part in 10(exp 16) or better over long time periods (greater than 10(exp 6) s). We consider the potential improvements this method could yield over conventional dynamical tests for unseen matter in the outer solar system possible now or anticipated in the near future.

Description: The intensities and optical depths of the (1, 1), (2, 2), and (2, 1) inversion transitions of ammonia can be calculated quite accurately without solving the equations of statistical equilibrium. A two-temperature partition function suffices. The excitation of the K-ladders can be approximated by using a temperature obtained from a two-level model with the (2, 1) and (1, 1) levels. Distribution of populations between the ladders is described with the kinetic temperature. This enables one to compute the (1, 1) and (2, 1) inversion transition excitation temperatures and optical depths. To compute the (2, 2) brightness temperatures, the fractional population of the (2, 2) doublet is computed from the population of the (1, 1) doublet using the 'true rotation temperature,' which is calculated using a three-level model with the (2, 1), (2, 2), and (1, 1) levels. In spite of some iterative steps, the calculation is quite fast.

Description: We report a 1.7 day ASCA X-ray observation of the 2.87 day binary Algol (Beta Per), centered on the secondary eclipse. Spectra accumulated for different intensity states show a prominent He-like iron K line at 6.7 keV. A two-temperature variable abundance plasma model applied to the spectra yielded temperatures of approximately 8 and approximately 30 MK. The modeled coronal abundances of Fe, O, Mg, Si, S, Ar, and Ca were a factor of 2-3 below the solar photospheric value, and N less than 0.1. These model abundance anomalies are similar to those found from the ASCA spectra of other late-type stars and may indicate either true deviations from solar abundances or problems with the assumptions and atomic physics of the plasma models. The X-ray light curve shows a factor of 2 increase in flux over 13 hr beginning in the middle of the secondary eclipse, with a shallow eclipse centered on phase 0.45. The orbital light curve is similar to that observed by ROSAT 1 year earlier. The rise in flux is caused by an increase in the emission measure of the higher temperature component. The intensity variation is not associated with any major change in the abundances or temperature of the two components.

Description: We present our analysis of ASCA PV phase observations of the elliptical galaxies NGC 1404 and NGC 4374 (M84). The average metallicities in the hot gas derived from the SIS spectra are exceptionally low, Z approximately 0.15 solar, while the temperatures are 'typical,' kT approximately 0.75 keV. We also place upper limits on intrinsic column densities. The low abundances lend support to the theory of Fe enrichment of intracluster media by protogalactic Type II supernova-driven winds and raise the possibility of a fundamental connection between baryon fraction, dissipation, and abundances in elliptical galaxies.

Description: The rate of occurrence of interplanetary discontinuities (ROID) is examined using Ulysses magnetic field and plasma data from 1 to 5 AU radial distance from the Sun and at high heliographic latitudes. We find two regions where the ROID is high: in stream-stream interaction regions and in Alfven wave trains. This latter feature is particularly obvious at high latitudes when Ulysses enters a high speed stream associated with a polar coronal hole. These streams are characterized by the presence of continuous, large-amplitude (Delta (vector 13)/absolute value of B is about 1-2 Alfven waves and an extraordinarily high ROID value (approximately 150 discontinuities/day). In a number of intervals examined, it is found that (rotational) discontinuities are an integral part of the Alfven waves. The nonlinear Alfven waves are spherically polarized, i.e., the tip of the perturbation vector resides on the surface of a sphere (a consequence of constant absolute value of B). The slowly rotating part of the wave rotates approximately 270 deg in phase. There is a slight arc in the B(sub 1) - B(sub 2) hodogram, suggesting an almost linear polarization. The phase rotation associated with the discontinuity is about 90 deg, lies in the same plane as the slowly rotaing part of the Alfven wave, and therefore completes the 360 deg phase rotation. The best description of the overall Alfven wave plus discontinuity is a spherical, arc-polarized, phase-steepened wave.

Description: Gamma-ray bursts have always been intriguing sources to study in terms of particle acceleration, but not since their discovery two decades ago has the theory of these objects been in such turmoil. Prior to the launch of Compton Gamma-Ray Observatory and observations by Burst and Transient Source Experiment (BATSE), there was strong evidence pointing to magnetized Galactic neutron stars as the sources of gamma-ray bursts. However, since BATSE the observational picture has changed dramatically, requiring much more distant and possibly cosmological sources. I review the history of gamma-ray burst theory from the era of growing consensus for nearby neutron stars to the recent explosion of halo and cosmological models and the impact of the present confusion on the particle acceleration problem.

Description: We discuss the fundamental ideas of particle acceleration in plasma shocks with emphasis on those features that are required to produce the 'universal' power-law spectrum. We compare shock acceleration with the more familiar second-order or stochastic acceleration and see that they are not too different in many respects. We discuss the features of shock acceleration that make it appealing and some of its problems as well.

Description: During the period from 1992 May to early 1992 November, the Energetic Gamma-Ray Experiment Telescope (EGRET) on board the Compton Gamma Ray Observatory obtained high-energy gamma-ray data for most of the sky. A total of 18 active galaxies have been seen with high certainty, and it is expected that more will be found in the data when a more thorough analysis is complete. All of those that have been seen are radio-loud quasars or BL Lacertae objects; most have already been identified as blazars. No Seyfert galaxies have been found thus far. If the spectra are represented as a power law in energy, spectral slopes ranging from approximately -1.7 to -2.4 are found. A wide range of z-values exits in the observed sample, eight having values in excess of 1.0. Time variations have been seen, with the timescale for a significant change being as short as days in at least one case. These results imply the existence of very large numbers of relativistic particles, probably close to the central object. Although a large extrapolation is required, their existence also suggests that these active galactic nuclei may be the source of the extragalactic cosmic rays.

Description: In a sample of impulsive bursts with rise times less than 30 s, a correlation between burst rise times and the frequency of maximum microwave emission has been found. The immplications for source structure and dynamics are discussed. Previously evidence was found that such bursts are caused by some propagating disturbance such as a shock wave or thermal conduction front. Combining that evidence with the microwave and hard X-ray spectral information suggests that the most rapid bursts are emitted from the most compact and intensely magnetized sources. The most rapid bursts also exhibited the hardest X-ray spectra, as published previously. These facts are important clues to understanding the physical process responsible for impulsive bursts. A model for the bursts is suggested, based on the observations and inferences described.

Description: The Fermi shock acceleration mechanism may be responsible for the production of high-energy cosmic rays in a wide variety of environments. Modeling of this phenomenon has largely focused on plane-parallel shocks, and one of the most promising techniques for its study is the Monte Carlo simulation of particle transport in shocked fluid flows. One of the principal problems in shock acceleration theory is the mechanism and efficiency of injection of particles from the thermal gas into the accelerated population. The Monte Carlo technique is ideally suited to addressing the injection problem directly, and previous applications of it to the quasi-parallel Earth bow shock led to very successful modeling of proton and heavy ion spectra, as well as other observed quantities. Recently this technique has been extended to oblique shock geometries, in which the upstream magnetic field makes a significant angle Theta(sub B1) to the shock normal. Spectral resutls from test particle Monte Carlo simulations of cosmic-ray acceleration at oblique, nonrelativistic shocks are presented. The results show that low Mach number shocks have injection efficiencies that are relatively insensitive to (though not independent of) the shock obliquity, but that there is a dramatic drop in efficiency for shocks of Mach number 30 or more as the obliquity increases above 15 deg. Cosmic-ray distributions just upstream of the shock reveal prominent bumps at energies below the thermal peak; these disappear far upstream but might be observable features close to astrophysical shocks.

Description: We report on combined ultraviolet and optical observations of early stages of the outburst of Nova Cygni 1992. We show that the initial optical rise and decline were produced by the dramatic increase in UV line opacity and its gradual subsequent lifting as the ejecta expand. The source of the M(sub v) - t(sub 3) relation is easily understood in light of these low-resolution UV spectra. The multiwavelength data confirm the theoretical prediction that a nova evolves at a constant bolometric luminosity in the early stages of outburst. We verify the prediction for this nova for at least the first month of the outburst. The detection of far-UV radiation with the Voyager Ultraviolet Spectrometer and the detection of an increasing X-ray flux with Roentgen Satellite (ROSAT) imply that this phase lasted for more than 6 months.

Description: The theory of cosmic-ray acceleration by subluminal magnetohydrodynamics (MHD) astrophysical shocks according to the test particle approximation is extended to highly oblique shocks where flow speeds which appear nonrelativistic in the shock rest frame appear on transformation to the de Hoffmann-Teller or E identically equal to 0 frame to have upstream flow speeds approaching c. Monte Carlo simulation shows that relative to the predictions of diffusion theory, as the upstream E identically equal to 0 frame velocity approaches c, flatter spectra and faster acceleration rates occur. These spectral and acceleration time changes are similar to those found for relativistic, parallel MHD shocks and may affect all nonthermal active galactic nuclei emission from relativistic electrons. Also there is an additional means of increasing the upper limit of the cosmic-ray spectrum expected from active galactic nuclei, although the approximation used may not accurately reproduce the spectral shape.

Description: The Advanced Satellite for Cosmology and Astrophysics (ASCA) satellite has obtained a moderate-resolution energy spectrum of E0102-72, the brightest Supernova Remnant (SNR) in the Small Magellanic Cloud (SMC). This paper reports on the first results of the analysis of the high quality spectrum of E0102-72. The spectrum shows resolved emission lines of He-like K alpha, H-like K alpha and K beta from oxygen, neon, and magnesium. The intensity ratios of these lines cannot be explained by a multi-component plasma model with uniform abundances, but requires abundance inhomogeneity in the plasma. We demonstrate how the spectral capabilities of the ASCA SIS make available new diagnostics of X-ray plasmas in a state of non-equilibrium ionization. Some interpretation based on the spectral analysis is also given.

Description: We present the first results from Advanced Satellite for Cosmology and Astrophysics (ASCAS) observations of the Cygnus Loop. A 22 min square of the north-east rim was observed. The imaging capability of ASCA permits clear indentification of the shock front structure. The energy resolving power of ASCA in the energy band between 0.5 keV to 10 keV reveals several emission lines from oxygen to silicon. The intensity ratio between the O VIII and O VII emission lines demonstrates that the plasma has not reached collisional ionization equilibrium. A nonequilibrium ionization model is applied to determine the gradients both in electron temperature T(sub e), and in electron density n(sub e). T(sub e) ranges from 0.22 keV at the shock front to 0.30 keV at 17 min inside the shock front. n(sub e) can be represented as a power law function of radius with an index of about 6 at an ambient gas density of 0.5 cm(exp -3). The plasma parameter n(sub e)t is about 2 x 10(exp 11) cm(exp -3), resulting in an age of about 20000 yr for the Loop.

Description: The spectra and images of the nearby jet galaxy NGC 4258 (M106) obtained with ASCA indicate presence of several distinct X-ray emission components. The emission above 3 keV is pointlike and coincident in position with the optical nucleus, exhibiting a hard (photon index approximately 1.78) and absorbed N(sub H) approximately 1.5 x 10(exp 23) cm(exp -2) spectrum. This provides clear evidence that NGC 4258 hosts an obscured active nucleus of low luminosity, about 4 x 10(exp 40) ergs s(exp -1) in 2-10 keV after removing the absorption. Iron K-line emission with an equivalent width 0.25 +/- 0.10 keV was detected. The emission below 1 keV is dominated by an extended approximately 4 min thin-thermal component with a temperature approximately 0.5 keV exhibiting atomic emission lines, possibly associated with the jet. There exists a third continuum component with an intermediate spectral hardness, which is brightest at approximately 1 min south-east of the nucleus.

Description: Spatially resolved energy spectra in the energy range 0.5-10 keV have been measured for the Centaurus cluster of galaxies with Advanced Satellite for Cosmology and Astrophysics (ASCA). Within 10 min (200 kpc) from the cluster center, the helium-like iron K emission line exhibits a dramatic increase toward the center rising from an equivalent width approximately 500 eV to approximately 1500 eV corresponding to an abundance change from 0.3 to 1.0 solar. The presence of strong iron L lines indicates an additional cool component (kT approximately 1 keV) within 10 min from the center. The cool component requires absorption in excess of the galactic value and this excess absorption increases towards the central region of the cluster. In the surrounding region with radius greater than 10 min, the spectra are well described by a single temperature thermal model with kT approximately 4 keV and spatially uniform abundances at about 0.3-0.4 times solar. The detection of metal-rich hot and cool gas in the cluster center implies a complex nature of the central cluster gas which is likely to be related to the presence of the central cD galaxy NGC 4696.

Description: We report preliminary results of Advanced Satellite for Cosmology and Astrophysics (ASCA) observations of two high-redshift quasars (zeta approximately 3) PKS 0438-436 and PKS 2126-158. We find spectral flattening towards low energies in both sources, first found with Roentgen Satellite (ROSAT) and interpreted as being due to excess X-ray absorption. However, the ROSAT data lacked the bandpass and sensitivity to unambiguously support this interpretation. We find that for PKS 2126-158, absorption in intervening low-redshift material along the line-of-sight provides a statistically better description of the data than absorption in material intrinsic to the quasar. The redshift of the putative absorber is required to be less than a few tenths. Given the redshift constraints and the fact that absorption is not common in lower redshift quasars, this result is not easily explained by a simple model. Either the absorber may be complex or else the intrinsic X-ray spectrum may be complex, or both. Compton-reflection dominated models are ruled out but a broken power-law model provides good fits to the data without requiring excess absorption. The break energy is at approximately 6-8 keV in the quasar frame and the energy index of the soft component is required to be extremely flat. This may mean that the X-ray emission in these objects is dominated by synchrotron losses as in BL Lac objects.

Description: A new computational method and algorithm, based on complex Fourier analysis, is used to derive the spectral density of plane and circularly polarized fluctuation components of the interplanetary magnetic field. Applications of the method have been made using HEOS 2 (1 AU), Pioneer 10 (5 AU), Pioneer 11 (20 AU), and International Cometary Explorer (ICE) (Giocabini-Zinner's comet) data sets. The results show the existence of circularly polarized magnetohydrodynamic (MHD) waves in all cases.

Description: We present the results of spectroscopic observations of 10 nucleated dwarf elliptical galaxies (dE's) in the Fornax cluster. The blue spectra of Fornax dE galaxies indicate a wide range of metallicities at a given luminosity, similar to those of intermediate to metal-rich globular clusters. Metal abundances derived in this paper are well correlated with optical colors and agree with previous spectroscopic results. A discrepancy with metallicities inferred from infrared colors is evident; possible causes include an intermediate age population and dilution of spectral features by a blue light excess. Dwarf ellipticals exhibit a wide variation of hydrogen line strength which points to a complex star formation history. Prominent Balmer absorption lines are the signature of a young stellar population in the nuclei of some (but not all) dE's, while moderately strong Balmer lines in relatively metal-rich dE's are more consistent with an extended main sequence. In a few metal-poor dE galaxies, the hydrogen lines are consisent with, or perhaps weaker than, those found in Galactic globulars of similar metallicity. In the limited magnitude range of this sample, there is no apparent correlation of metallicity either with effective and central surface brightness, or with total and nuclear magnitudes. The velocity distribution of the Fornax dwarfs is flatter than that of brighter galaxies at the 75% confidence level, possibly indicating a difference in the kinematics of the two samples.

Description: This paper presents comprehensive results on the spectra of 30 bright gamma ray bursts (GRBs) as observed by the Spectroscopy Detectors (SDs) of the Burst And Transient Source Experiment (BATSE). The data selection was strict in including only spectra that are of high reliability for continuum shape studies. This BATSE Spectroscopy Catalog presents fluences, model fits (for five spectral models for three energy ranges), and photon spectra in a standard manner for each burst. Complete information is provided to describe the data selection and analysis procedures. The catalog results are also presented in electronic format (from the Compton Observatory Science Support Center) and CD-ROM format (AAS CD-ROM series, Vol. 2). These electronic formats also present the count spectra and detector response matrices so as to allow for independent study and fitting by researchers outside the BATSE Team. This BATSE Spectroscopy Catalog complements the catalog from BATSE Large Area Detector (LAD) data by Fishman et al. (1994).

Description: The implications of recent near-infrared imaging and spectroscopy of the Galactic center stellar cluster are discussed. The central parsec appears to be powered by a cluster of hot, massive stars of which the IRS 16 complex is the central core. In the 1 to 2 micrometer band, the brightest members of this cluster are 10 to 15 HeI/HI emission line stars that can be characterized as approximately 20000 K, helium rich, very luminous supergiants. The He-I/H-I stars can account for a major fraction of the total and Lyman continuum luminosity of the central parsec, but hotter, earlier type stars are probably required in addition to account for the He-continuum. The brightest cool stars in the central parsec are red supergiants, and asymptotic giant branch stars. Two scenarios for the evolution of the central stellar core are presented: one involves a small star formation burst years ago that was the result of substantial prior gas influx into the core. In this scenario the Galactic center is presently in a short-lived, post-main sequence 'wind phase'. The second scenario involves the buildup of massive stars by sequential merging of lower mass stars. The intense mass loss from the hot stars probably affects strongly the gas dynamics in the central 0.1 pc and may prevent gas to accelerate onto the possible central hole.

Description: We present Very Large Array (VLA) radio continuum and Infrared Astronomy Satellite (IRAS) far-infrared (FIR) observations of 16 low luminosity galaxies of mostly low surface brightness. All galaxies had previously claimed single dish radio continuum detections. However, at the frequencies of our observations (1.49 and 8.48 GHz), we find significant radio emission for two objects only. We show that the other previously claimed detections are due to confusion with physically unrelated background sources. This implies a low radio continuum detection rate for these galaxies. Re-reduced IRAS scans yield significant far-infrared flux densities in at least one IRAS band for 6 of the 16 galaxies. These, together with the FIR and radio continuum upper limits, are consistent with the well established radio/FIR relation, where most of our galaxies populate the low-luminosity end. From the radio continuum and FIR flux densities and their upper limits we estimate the current star formation rates and demonstrate that the galaxies are currently passive in forming stars, in agreement with previous optical investigations. There is an indication that the galaxies were forming stars more intensively averaged over their lifetime than they are presently.

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: 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.