ALBERT

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Beschreibung: The analysis of balloon envelopes by the finite element (FE) method is plagued by convergence problems. A pratical FE analysis approach is based on the fact that in thin shells with non-zero Gaussian curvature the membrane solution component is essentially decoupled from the bending solution component. A proxy-problem is solved by using a small artificial bending stiffness that assures convergence without significantly affecting the membrane solution component. This approach has been previously validated on slightly overpressurized balloon envelopes. Extensions of this approach to more difficult problems in the structural analysis of balloon envelopes are presented. The convergence forcing modelling measures are discussed. Implications of the findings of the analysis results to future balloon designs are also discussed.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Beschreibung: The Composite Plate Buckling Analysis Program (COMPPAP) was written to help engineers determine buckling loads of orthotropic (or isotropic) irregularly shaped plates without requiring hand calculations from design curves or extensive finite element modeling. COMPPAP is a one element finite element program that utilizes high-order displacement functions. The high order of the displacement functions enables the user to produce results more accurate than traditional h-finite elements. This program uses these high-order displacement functions to perform a plane stress analysis of a general plate followed by a buckling calculation based on the stresses found in the plane stress solution. The current version assumes a flat plate (constant thickness) subject to a constant edge load (normal or shear) on one or more edges. COMPPAP uses the power method to find the eigenvalues of the buckling problem. The power method provides an efficient solution when only one eigenvalue is desired. Once the eigenvalue is found, the eigenvector, which corresponds to the plate buckling mode shape, results as a by-product. A positive feature of the power method is that the dominant eigenvalue is the first found, which is this case is the plate buckling load. The reported eigenvalue expresses a load factor to induce plate buckling. COMPPAP is written in ANSI FORTRAN 77. Two machine versions are available from COSMIC: a PC version (MSC-22428), which is for IBM PC 386 series and higher computers and compatibles running MS-DOS; and a UNIX version (MSC-22286). The distribution medium for both machine versions includes source code for both single and double precision versions of COMPPAP. The PC version includes source code which has been optimized for implementation within DOS memory constraints as well as sample executables for both the single and double precision versions of COMPPAP. The double precision versions of COMPPAP have been successfully implemented on an IBM PC 386 compatible running MS-DOS, a Sun4 series computer running SunOS, an HP-9000 series computer running HP-UX, and a CRAY X-MP series computer running UNICOS. COMPPAP requires 1Mb of RAM and the BLAS and LINPACK math libraries, which are included on the distribution medium. The COMPPAP documentation provides instructions for using the commercial post-processing package PATRAN for graphical interpretation of COMPPAP output. The UNIX version includes two electronic versions of the documentation: one in LaTex format and one in PostScript format. The standard distribution medium for the PC version (MSC-22428) is a 5.25 inch 1.2Mb MS-DOS format diskette. The standard distribution medium for the UNIX version (MSC-22286) is a .25 inch streaming magnetic tape cartridge (Sun QIC-24) in UNIX tar format. For the UNIX version, alternate distribution media and formats are available upon request. COMPPAP was developed in 1992.

Beschreibung: The Composite Plate Buckling Analysis Program (COMPPAP) was written to help engineers determine buckling loads of orthotropic (or isotropic) irregularly shaped plates without requiring hand calculations from design curves or extensive finite element modeling. COMPPAP is a one element finite element program that utilizes high-order displacement functions. The high order of the displacement functions enables the user to produce results more accurate than traditional h-finite elements. This program uses these high-order displacement functions to perform a plane stress analysis of a general plate followed by a buckling calculation based on the stresses found in the plane stress solution. The current version assumes a flat plate (constant thickness) subject to a constant edge load (normal or shear) on one or more edges. COMPPAP uses the power method to find the eigenvalues of the buckling problem. The power method provides an efficient solution when only one eigenvalue is desired. Once the eigenvalue is found, the eigenvector, which corresponds to the plate buckling mode shape, results as a by-product. A positive feature of the power method is that the dominant eigenvalue is the first found, which is this case is the plate buckling load. The reported eigenvalue expresses a load factor to induce plate buckling. COMPPAP is written in ANSI FORTRAN 77. Two machine versions are available from COSMIC: a PC version (MSC-22428), which is for IBM PC 386 series and higher computers and compatibles running MS-DOS; and a UNIX version (MSC-22286). The distribution medium for both machine versions includes source code for both single and double precision versions of COMPPAP. The PC version includes source code which has been optimized for implementation within DOS memory constraints as well as sample executables for both the single and double precision versions of COMPPAP. The double precision versions of COMPPAP have been successfully implemented on an IBM PC 386 compatible running MS-DOS, a Sun4 series computer running SunOS, an HP-9000 series computer running HP-UX, and a CRAY X-MP series computer running UNICOS. COMPPAP requires 1Mb of RAM and the BLAS and LINPACK math libraries, which are included on the distribution medium. The COMPPAP documentation provides instructions for using the commercial post-processing package PATRAN for graphical interpretation of COMPPAP output. The UNIX version includes two electronic versions of the documentation: one in LaTex format and one in PostScript format. The standard distribution medium for the PC version (MSC-22428) is a 5.25 inch 1.2Mb MS-DOS format diskette. The standard distribution medium for the UNIX version (MSC-22286) is a .25 inch streaming magnetic tape cartridge (Sun QIC-24) in UNIX tar format. For the UNIX version, alternate distribution media and formats are available upon request. COMPPAP was developed in 1992.

Beschreibung: Structural flaws and cracks may grow under fatigue inducing loads and, upon reaching a critical size, cause structural failure to occur. The growth of these flaws and cracks may occur at load levels well below the ultimate load bearing capability of the structure. The Fatigue Crack Growth Computer Program, NASA/FLAGRO, was developed as an aid in predicting the growth of pre-existing flaws and cracks in structural components of space systems. The earlier version of the program, FLAGRO4, was the primary analysis tool used by Rockwell International and the Shuttle subcontractors for fracture control analysis on the Space Shuttle. NASA/FLAGRO is an enhanced version of the program and incorporates state-of-the-art improvements in both fracture mechanics and computer technology. NASA/FLAGRO provides the fracture mechanics analyst with a computerized method of evaluating the "safe crack growth life" capabilities of structural components. NASA/FLAGRO could also be used to evaluate the damage tolerance aspects of a given structural design. The propagation of an existing crack is governed by the stress field in the vicinity of the crack tip. The stress intensity factor is defined in terms of the relationship between the stress field magnitude and the crack size. The propagation of the crack becomes catastrophic when the local stress intensity factor reaches the fracture toughness of the material. NASA/FLAGRO predicts crack growth using a two-dimensional model which predicts growth independently in two directions based on the calculation of stress intensity factors. The analyst can choose to use either a crack growth rate equation or a nonlinear interpolation routine based on tabular data. The growth rate equation is a modified Forman equation which can be converted to a Paris or Walker equation by substituting different values into the exponent. This equation provides accuracy and versatility and can be fit to data using standard least squares methods. Stress-intensity factor numerical values can be computed for making comparisons or checks of solutions. NASA/FLAGRO can check for failure of a part-through crack in the mode of a through crack when net ligament yielding occurs. NASA/FLAGRO has a number of special subroutines and files which provide enhanced capabilities and easy entry of data. These include crack case solutions, cyclic load spectrums, nondestructive examination initial flaw sizes, table interpolation, and material properties. The materials properties files are divided into two types, a user defined file and a fixed file. Data is entered and stored in the user defined file during program execution, while the fixed file contains already coded-in property value data for many different materials. Prompted input from CRT terminals consists of initial crack definition (which can be defined automatically), rate solution type, flaw type and geometry, material properties (if they are not in the built-in tables of material data), load spectrum data (if not included in the loads spectrum file), and design limit stress levels. NASA/FLAGRO output includes an echo of the input with any error or warning messages, the final crack size, whether or not critical crack size has been reached for the specified stress level, and a life history profile of the crack propagation. NASA/FLAGRO is modularly designed to facilitate revisions and operation on minicomputers. The program was implemented on a DEC VAX 11/780 with the VMS operating system. NASA/FLAGRO is written in FORTRAN77 and has a memory requirement of 1.4 MB. The program was developed in 1986.

Beschreibung: Delaminations near the outer surface of a laminate are susceptible to local buckling and buckling-induced delamination propagation when the laminate is subjected to transverse impact loading. This results in a loss of stiffness and strength. TRBUCKL is an unique dynamic delamination buckling and delamination propagation analysis capability that can be incorporated into the structural analysis program, NASTRAN. This capability will aid engineers in the design of structures incorporating composite laminates. The capability consists of: (1) a modification of the direct time integration solution sequence which provides a new analysis algorithm that can be used to predict delamination buckling in a laminate subjected to dynamic loading; and (2) a new method of modeling the composite laminate using plate bending elements and multipoint constraints. The capability now exists to predict the time at which the onset of dynamic delamination buckling occurs, the dynamic buckling mode shape, and the dynamic delamination strain energy release rate. A procedure file for NASTRAN, TRBUCKL predicts both impact induced buckling in composite laminates with initial delaminations and the strain energy release rate due to extension of the delamination. In addition, the file is useful in calculating the dynamic delamination strain energy release rate for a composite laminate under impact loading. This procedure simplifies the simulation of progressive crack extension. TRBUCKL has been incorporated into COSMIC NASTRAN. TRBUCKL is a DMAP Alter for NASTRAN. It is intended for use only with the COSMIC NASTRAN Direct Transient Analysis (RF 9) solution sequence. The program is available as a listing only. TRBUCKL was developed in 1987.

Beschreibung: Predictions of fatigue crack growth behavior can be made with the Fatigue Crack Growth Structural Analysis (FASTRAN II) computer program. As cyclic loads are applied to a selected crack configuration with an initial crack size, FASTRAN II predicts crack growth as a function of cyclic load history until either a desired crack size is reached or failure occurs. FASTRAN II is based on plasticity-induced crack-closure behavior of cracks in metallic materials and accounts for load-interaction effects, such as retardation and acceleration, under variable-amplitude loading. The closure model is based on the Dugdale model with modifications to allow plastically deformed material to be left along the crack surfaces as the crack grows. Plane stress and plane strain conditions, as well as conditions between these two, can be simulated in FASTRAN II by using a constraint factor on tensile yielding at the crack front to approximately account for three-dimensional stress states. FASTRAN II contains seventeen predefined crack configurations (standard laboratory fatigue crack growth rate specimens and many common crack configurations found in structures); and the user can define one additional crack configuration. The baseline crack growth rate properties (effective stress-intensity factor against crack growth rate) may be given in either equation or tabular form. For three-dimensional crack configurations, such as surface cracks or corner cracks at holes or notches, the fatigue crack growth rate properties may be different in the crack depth and crack length directions. Final failure of the cracked structure can be modelled with fracture toughness properties using either linear-elastic fracture mechanics (brittle materials), a two-parameter fracture criterion (brittle to ductile materials), or plastic collapse (extremely ductile materials). The crack configurations in FASTRAN II can be subjected to either constant-amplitude, variable-amplitude or spectrum loading. The applied loads may be either tensile or compressive. Several standardized aircraft flight-load histories, such as TWIST, Mini-TWIST, FALSTAFF, Inverted FALSTAFF, Felix and Gaussian, are included as options. FASTRAN II also includes two other methods that will help the user input spectrum load histories. The two methods are: (1) a list of stress points, and (2) a flight-by-flight history of stress points. Examples are provided in the user manual. Developed as a research program, FASTRAN II has successfully predicted crack growth in many metallic materials under various aircraft spectrum loading. A computer program DKEFF which is a part of the FASTRAN II package was also developed to analyze crack growth rate data from laboratory specimens to obtain the effective stress-intensity factor against crack growth rate relations used in FASTRAN II. FASTRAN II is written in standard FORTRAN 77. It has been successfully compiled and implemented on Sun4 series computers running SunOS and on IBM PC compatibles running MS-DOS using the Lahey F77L FORTRAN compiler. Sample input and output data are included with the FASTRAN II package. The UNIX version requires 660K of RAM for execution. The standard distribution medium for the UNIX version (LAR-14865) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. It is also available on a 3.5 inch diskette in UNIX tar format. The standard distribution medium for the MS-DOS version (LAR-14944) is a 5.25 inch 360K MS-DOS format diskette. The contents of the diskette are compressed using the PKWARE archiving tools. The utility to unarchive the files, PKUNZIP.EXE, is included. The program was developed in 1984 and revised in 1992. Sun4 and SunOS are trademarks of Sun Microsystems, Inc. IBM PC is a trademark of International Business Machines Corp. MS-DOS is a trademark of Microsoft, Inc. F77L is a trademark of the Lahey Computer Systems, Inc. UNIX is a registered trademark of AT&T Bell Laboratories. PKWARE and PKUNZIP are trademarks of PKWare, Inc.

Beschreibung: The Panel Analysis and Sizing Code (PASCO) was developed for the buckling and vibration analysis and sizing of prismatic structures having an arbitrary cross section. PASCO is primarily intended for analyzing and sizing stiffened panels made of laminated orthotropic materials and is of particular value in analyzing and sizing filamentary composite structures. When used in the analysis mode, PASCO calculates laminate stiffnesses, lamina stress and strains (including the effects of temperature and panel bending), buckling loads, vibration frequencies, and overall panel stiffness. When used in the sizing mode, PASCO adjusts sizing variables to provide a low-mass panel design that carries a set of specified loadings without exceeding buckling or material strength allowables and that meets other design requirements such as upper and lower bounds on sizing variables, upper and lower bounds on overall bending, extensional and shear stiffnesses, and lower bounds on vibration frequencies. Although emphasis in PASCO is placed on flat panels having several identical bays, the only restriction on configuration modeling is that the structure is assumed to be prismatic. In addition, it is assumed that loads and temperatures do not vary along the length of a panel. Because of their wide application in aerospace structures, stiffened panels are readily handled by PASCO. The panel cross section may be composed of an arbitrary assemblage of thin, flat, rectangular plate elements that are connected together along their longitudinal edges. Each plate element consists of a balanced symmetric laminate of any number of layers of orthotropic material. Any group of element widths, layer thicknesses, and layer orientation angles may be selected as sizing variables. Substructuring is available to increase the efficiency of the analysis and to simplify the modeling of complicated structures. The Macintosh version of PASCO includes an interactive, graphic preprocessor called MacPASCO. The main objective of MacPASCO is to make the use of PASCO faster, simpler, and less error-prone. By using a graphical user interface (GUI), MacPASCO simplifies the specification of panel geometry and reduces user input errors, thus making the modeling and analysis of panel designs more efficient. The user draws the initial structural geometry on the computer screen, then uses a combination of graphic and text inputs to: refine the structural geometry, specify information required for analysis such as panel load conditions, and define design variables and constraints for minimum-mass optimization. Composite panel design is an ideal application because the graphical user interface can: serve as a visual aid, eliminate the tedious aspects of text-based input, and eliminate many sources of input errors. The current version of MacPASCO does not implement all the modeling features of PASCO, but has been found to be sufficient for many users. Many difficulties common to text-based inputs are avoided because MacPASCO uses a GUI. First, the graphic displays eliminate syntax errors, like misplaced commas and incorrect command names, because there is no text-based syntax. Second, graphic displays allow the user to see the geometry as it is created and immediately detect and correct any errors. Third, MacPASCO's drawing tools have been designed to avoid modeling errors. Fourth, the graphic displays make revisions to existing structural designs much easier and less error-prone by eliminating the need for the user to conceptualize the text input as geometry. The user can work directly with the geometry displayed on the screen. Finally, MacPASCO automatically generates the correct PASCO input from the geometry displayed on the screen. This input file can be used with any machine version of PASCO to actually perform the analysis and sizing and to output results. The DEC VAX version of PASCO is written in FORTRAN IV for batch execution and has been implemented on a DEC VAX series computer. The Macintosh version of PASCO was developed for Macintosh II series computers with at least 2Mb of RAM running MPW Pascal 3.0 and Language Systems FORTRAN 2.0 under the MPW programming environment. It includes MPW compatible makefiles for compiling the source code. The Macintosh version uses input files compatible with versions of PASCO running on different platforms. MacPASCO is written in Macintosh Programmers Workbench 3.0, MPW Pascal 3.0, and MacAPP 2.0. The Pascal source code is included on the distribution diskette. MacAPP is a development library which is not included. MacPASCO requires a Mac Plus, SE/30, or MacII, IIx, IIcx, IIci, or IIfx running System 6.0 or greater. MacPASCO is System 7.0 compatible. A minimum of 2Mb of RAM is required for execution. The Macintosh version of PASCO is distributed on four 3.5 inch 800K Macintosh format diskettes. The DEC VAX version is distributed on a 9-track 1600 BPI magnetic tape. The PASCO program was developed in 1981, adapted to the DEC VAX in 1983 and to the Macintosh in 1991. MacPASCO was released in 1992.

Beschreibung: PLAN2D is a FORTRAN computer program for the plastic analysis of planar rigid frame structures. Given a structure and loading pattern as input, PLAN2D calculates the ultimate load that the structure can sustain before collapse. Element moments and plastic hinge rotations are calculated for the ultimate load. The location of hinges required for a collapse mechanism to form are also determined. The program proceeds in an iterative series of linear elastic analyses. After each iteration the resulting elastic moments in each member are compared to the reserve plastic moment capacity of that member. The member or members that have moments closest to their reserve capacity will determine the minimum load factor and the site where the next hinge is to be inserted. Next, hinges are inserted and the structural stiffness matrix is reformulated. This cycle is repeated until the structure becomes unstable. At this point the ultimate collapse load is calculated by accumulating the minimum load factor from each previous iteration and multiplying them by the original input loads. PLAN2D is based on the program STAN, originally written by Dr. E.L. Wilson at U.C. Berkeley. PLAN2D has several limitations: 1) Although PLAN2D will detect unloading of hinges it does not contain the capability to remove hinges; 2) PLAN2D does not allow the user to input different positive and negative moment capacities and 3) PLAN2D does not consider the interaction between axial and plastic moment capacity. Axial yielding and buckling is ignored as is the reduction in moment capacity due to axial load. PLAN2D is written in FORTRAN and is machine independent. It has been tested on an IBM PC and a DEC MicroVAX. The program was developed in 1988.

Beschreibung: Predictions of fatigue crack growth behavior can be made with the Fatigue Crack Growth Structural Analysis (FASTRAN II) computer program. As cyclic loads are applied to a selected crack configuration with an initial crack size, FASTRAN II predicts crack growth as a function of cyclic load history until either a desired crack size is reached or failure occurs. FASTRAN II is based on plasticity-induced crack-closure behavior of cracks in metallic materials and accounts for load-interaction effects, such as retardation and acceleration, under variable-amplitude loading. The closure model is based on the Dugdale model with modifications to allow plastically deformed material to be left along the crack surfaces as the crack grows. Plane stress and plane strain conditions, as well as conditions between these two, can be simulated in FASTRAN II by using a constraint factor on tensile yielding at the crack front to approximately account for three-dimensional stress states. FASTRAN II contains seventeen predefined crack configurations (standard laboratory fatigue crack growth rate specimens and many common crack configurations found in structures); and the user can define one additional crack configuration. The baseline crack growth rate properties (effective stress-intensity factor against crack growth rate) may be given in either equation or tabular form. For three-dimensional crack configurations, such as surface cracks or corner cracks at holes or notches, the fatigue crack growth rate properties may be different in the crack depth and crack length directions. Final failure of the cracked structure can be modelled with fracture toughness properties using either linear-elastic fracture mechanics (brittle materials), a two-parameter fracture criterion (brittle to ductile materials), or plastic collapse (extremely ductile materials). The crack configurations in FASTRAN II can be subjected to either constant-amplitude, variable-amplitude or spectrum loading. The applied loads may be either tensile or compressive. Several standardized aircraft flight-load histories, such as TWIST, Mini-TWIST, FALSTAFF, Inverted FALSTAFF, Felix and Gaussian, are included as options. FASTRAN II also includes two other methods that will help the user input spectrum load histories. The two methods are: (1) a list of stress points, and (2) a flight-by-flight history of stress points. Examples are provided in the user manual. Developed as a research program, FASTRAN II has successfully predicted crack growth in many metallic materials under various aircraft spectrum loading. A computer program DKEFF which is a part of the FASTRAN II package was also developed to analyze crack growth rate data from laboratory specimens to obtain the effective stress-intensity factor against crack growth rate relations used in FASTRAN II. FASTRAN II is written in standard FORTRAN 77. It has been successfully compiled and implemented on Sun4 series computers running SunOS and on IBM PC compatibles running MS-DOS using the Lahey F77L FORTRAN compiler. Sample input and output data are included with the FASTRAN II package. The UNIX version requires 660K of RAM for execution. The standard distribution medium for the UNIX version (LAR-14865) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. It is also available on a 3.5 inch diskette in UNIX tar format. The standard distribution medium for the MS-DOS version (LAR-14944) is a 5.25 inch 360K MS-DOS format diskette. The contents of the diskette are compressed using the PKWARE archiving tools. The utility to unarchive the files, PKUNZIP.EXE, is included. The program was developed in 1984 and revised in 1992. Sun4 and SunOS are trademarks of Sun Microsystems, Inc. IBM PC is a trademark of International Business Machines Corp. MS-DOS is a trademark of Microsoft, Inc. F77L is a trademark of the Lahey Computer Systems, Inc. UNIX is a registered trademark of AT&T Bell Laboratories. PKWARE and PKUNZIP are trademarks of PKWare, Inc.

Beschreibung: The SNIP program is a FORTRAN computer code that generates NASTRAN structural model thermal loads when given SINDA (or similar thermal model) temperature results. SNIP correlates thermal nodes to structural elements to interface SINDA finite difference thermal models with NASTRAN finite element structural models. Node-to-element correlation includes determining which SINDA nodes should be related to each NASTRAN element and calculating a weighing factor for temperatures associated with each element-related thermal node. SNIP provides structural model thermal loads that accurately reflect thermal model results while reducing the time required to interface thermal and structural models as compared to other methods. SNIP uses thermal model geometry to search the three-dimensional space around each structural element for the nearest thermal nodes. Thermal model geometry is the combination of standard thermal model temperature results from SINDA and structural model geometry from NASTRAN. Thermal and structural models must both be defined in the same, single Cartesian coordinate system. The thermal nodes located nearest each element are used to determine element temperature for thermal distortion and stress analysis. The program shapes the three-dimensional search region while the user controls the size. With these region specifications, the numerical coding of thermal nodes, and the structural element numbers; the code can provide for the separation of substructures during correlation. The input to SNIP contains a file of thermal model temperature results and a physical location of each thermal node in three-dimensional space, combined in a SNIP-unique format. The input also contains a standard NASTRAN input deck for a model made up of plate, shell, beam, and bar elements. SNIP supports the CTRIA, CQUAD, CBAR, and CBEAM elements of NASTRAN. The user adjusts the input parameters in the source code which control the node-to-element correlation. The program outputs NASTRAN element temperature load cards for each element and NASTRAN case control cards for each temperature load set. SNIP also outputs a list of elements that contains the numbers of the SINDA nodes related to each NASTRAN element and the weight that is given to each node in temperature calculations. SNIP is written in ANSI standard FORTRAN 77. The PC version requires a PC FORTRAN compiler and has compiled successfully using Lahey FORTRAN v. 3.0. A core memory of 300k is recommended. The program was developed in 1987.

Beschreibung: NPLOT is an interactive computer graphics program for plotting undeformed and deformed NASTRAN finite element models (FEMs). Although there are many commercial codes already available for plotting FEMs, these have limited use due to their cost, speed, and lack of features to view BAR elements. NPLOT was specifically developed to overcome these limitations. On a vector type graphics device the two best ways to show depth are by hidden line plotting or haloed line plotting. A hidden line algorithm generates views of models with all hidden lines removed, and a haloed line algorithm displays views with aft lines broken in order to show depth while keeping the entire model visible. A haloed line algorithm is especially useful for plotting models composed of many line elements and few surface elements. The most important feature of NPLOT is its ability to create both hidden line and haloed line views accurately and much more quickly than with any other existing hidden or haloed line algorithms. NPLOT is also capable of plotting a normal wire frame view to display all lines of a model. NPLOT is able to aid in viewing all elements, but it has special features not generally available for plotting BAR elements. These features include plotting of TRUE LENGTH and NORMALIZED offset vectors and orientation vectors. Standard display operations such as rotation and perspective are possible, but different view planes such as X-Y, Y-Z, and X-Z may also be selected. Another display option is the Z-axis cut which allows a portion of the fore part of the model to be cut away to reveal details of the inside of the model. A zoom function is available to terminals with a locator (graphics cursor, joystick, etc.). The user interface of NPLOT is designed to make the program quick and easy to use. A combination of menus and commands with help menus for detailed information about each command allows experienced users greater speed and efficiency. Once a plot is on the screen the interface becomes command driven, enabling the user to manipulate the display or execute a command without having to return to the menu. NPLOT is also able to plot deformed shapes allowing it to perform post-processing. The program can read displacements, either static displacements or eigenvectors, from a MSC/NASTRAN F06 file or a UAI/NASTRAN PRT file. The displacements are written into a unformatted scratch file where they are available for rapid access when the user wishes to display a deformed shape. All subcases or mode shapes can be read in at once. Then it is easy to enable the deformed shape, to change subcases or mode shapes and to change the scale factor for subsequent plots. NPLOT is written in VAX FORTRAN for DEC VAX series computers running VMS. As distributed, the NPLOT source code makes calls to the DI3000 graphics package from Precision Visuals; however, a set of interface routines is provided to translate the DI3000 calls into Tektronix PLOT10/TCS graphics library calls so that NPLOT can use the standard Tektronix 4010 which many PC terminal emulation software programs support. NPLOT is available in VAX BACKUP format on a 9-track 1600 BPI DEC VAX BACKUP format magnetic tape (standard media) or a TK50 tape cartridge. This program was developed in 1991. DEC, VAX, VMS, and TK50 are trademarks of Digital Equipment Corporation. Tektronix, PLOT10, and TCS are trademarks of Tektronix, Inc. DI3000 is a registered trademark of Precision Visuals, Inc. NASTRAN is a registered trademark of the National Aeronautics and Space Administration. MSC/ is a trademark of MacNeal-Schwendler Corporation. UAI is a trademark of Universal Analytics, Inc.

Beschreibung: LATDYN is a computer code for modeling the Large Angle Transient DYNamics of flexible articulating structures and mechanisms involving joints about which members rotate through large angles. LATDYN extends and brings together some of the aspects of Finite Element Structural Analysis, Multi-Body Dynamics, and Control System Analysis; three disciplines that have been historically separate. It combines significant portions of their distinct capabilities into one single analysis tool. The finite element formulation for flexible bodies in LATDYN extends the conventional finite element formulation by using a convected coordinate system for constructing the equation of motion. LATDYN's formulation allows for large displacements and rotations of finite elements subject to the restriction that deformations within each are small. Also, the finite element approach implemented in LATDYN provides a convergent path for checking solutions simply by increasing mesh density. For rigid bodies and joints LATDYN borrows extensively from methodology used in multi-body dynamics where rigid bodies may be defined and connected together through joints (hinges, ball, universal, sliders, etc.). Joints may be modeled either by constraints or by adding joint degrees of freedom. To eliminate error brought about by the separation of structural analysis and control analysis, LATDYN provides symbolic capabilities for modeling control systems which are integrated with the structural dynamic analysis itself. Its command language contains syntactical structures which perform symbolic operations which are also interfaced directly with the finite element structural model, bypassing the modal approximation. Thus, when the dynamic equations representing the structural model are integrated, the equations representing the control system are integrated along with them as a coupled system. This procedure also has the side benefit of enabling a dramatic simplification of the user interface for modeling control systems. Three FORTRAN computer programs, the LATDYN Program, the Preprocessor, and the Postprocessor, make up the collective LATDYN System. The Preprocessor translates user commands into a form which can be used while the LATDYN program provides the computational core. The Postprocessor allows the user to interactively plot and manage a database of LATDYN transient analysis results. It also includes special facilities for modeling control systems and for programming changes to the model which take place during analysis sequence. The documentation includes a Demonstration Problem Manual for the evaluation and verification of results and a Postprocessor guide. Because the program should be viewed as a byproduct of research on technology development, LATDYN's scope is limited. It does not have a wide library of finite elements, and 3-D Graphics are not available. Nevertheless, it does have a measure of "user friendliness". The LATDYN program was developed over a period of several years and was implemented on a CDC NOS/VE & Convex Unix computer. It is written in FORTRAN 77 and has a virtual memory requirement of 1.46 MB. The program was validated on a DEC MICROVAX operating under VMS 5.2.

Beschreibung: The increasing number of applications of fiber-reinforced composites in industry demands a detailed understanding of their material properties and behavior. A three-dimensional finite-element computer program called PAFAC (Plastic and Failure Analysis of Composites) has been developed for the elastic-plastic analysis of fiber-reinforced composite materials and structures. The evaluation of stresses and deformations at edges, cut-outs, and joints is essential in understanding the strength and failure for metal-matrix composites since the onset of plastic yielding starts very early in the loading process as compared to the composite's ultimate strength. Such comprehensive analysis can only be achieved by a finite-element program like PAFAC. PAFAC is particularly suited for the analysis of laminated metal-matrix composites. It can model the elastic-plastic behavior of the matrix phase while the fibers remain elastic. Since the PAFAC program uses a three-dimensional element, the program can also model the individual layers of the laminate to account for thickness effects. In PAFAC, the composite is modeled as a continuum reinforced by cylindrical fibers of vanishingly small diameter which occupy a finite volume fraction of the composite. In this way, the essential axial constraint of the phases is retained. Furthermore, the local stress and strain fields are uniform. The PAFAC finite-element solution is obtained using the displacement method. Solution of the nonlinear equilibrium equations is obtained with a Newton-Raphson iteration technique. The elastic-plastic behavior of composites consisting of aligned, continuous elastic filaments and an elastic-plastic matrix is described in terms of the constituent properties, their volume fractions, and mutual constraints between phases indicated by the geometry of the microstructure. The program uses an iterative procedure to determine the overall response of the laminate, then from the overall response determines the stress state in each phase of the composite material. Failure of the fibers or matrix within an element can also be modeled by PAFAC. PAFAC is written in FORTRAN IV for batch execution and has been implemented on a CDC CYBER 170 series computer with a segmented memory requirement of approximately 66K (octal) of 60 bit words. PAFAC was developed in 1982.

Beschreibung: The Panel Analysis and Sizing Code (PASCO) was developed for the buckling and vibration analysis and sizing of prismatic structures having an arbitrary cross section. PASCO is primarily intended for analyzing and sizing stiffened panels made of laminated orthotropic materials and is of particular value in analyzing and sizing filamentary composite structures. When used in the analysis mode, PASCO calculates laminate stiffnesses, lamina stress and strains (including the effects of temperature and panel bending), buckling loads, vibration frequencies, and overall panel stiffness. When used in the sizing mode, PASCO adjusts sizing variables to provide a low-mass panel design that carries a set of specified loadings without exceeding buckling or material strength allowables and that meets other design requirements such as upper and lower bounds on sizing variables, upper and lower bounds on overall bending, extensional and shear stiffnesses, and lower bounds on vibration frequencies. Although emphasis in PASCO is placed on flat panels having several identical bays, the only restriction on configuration modeling is that the structure is assumed to be prismatic. In addition, it is assumed that loads and temperatures do not vary along the length of a panel. Because of their wide application in aerospace structures, stiffened panels are readily handled by PASCO. The panel cross section may be composed of an arbitrary assemblage of thin, flat, rectangular plate elements that are connected together along their longitudinal edges. Each plate element consists of a balanced symmetric laminate of any number of layers of orthotropic material. Any group of element widths, layer thicknesses, and layer orientation angles may be selected as sizing variables. Substructuring is available to increase the efficiency of the analysis and to simplify the modeling of complicated structures. The Macintosh version of PASCO includes an interactive, graphic preprocessor called MacPASCO. The main objective of MacPASCO is to make the use of PASCO faster, simpler, and less error-prone. By using a graphical user interface (GUI), MacPASCO simplifies the specification of panel geometry and reduces user input errors, thus making the modeling and analysis of panel designs more efficient. The user draws the initial structural geometry on the computer screen, then uses a combination of graphic and text inputs to: refine the structural geometry, specify information required for analysis such as panel load conditions, and define design variables and constraints for minimum-mass optimization. Composite panel design is an ideal application because the graphical user interface can: serve as a visual aid, eliminate the tedious aspects of text-based input, and eliminate many sources of input errors. The current version of MacPASCO does not implement all the modeling features of PASCO, but has been found to be sufficient for many users. Many difficulties common to text-based inputs are avoided because MacPASCO uses a GUI. First, the graphic displays eliminate syntax errors, like misplaced commas and incorrect command names, because there is no text-based syntax. Second, graphic displays allow the user to see the geometry as it is created and immediately detect and correct any errors. Third, MacPASCO's drawing tools have been designed to avoid modeling errors. Fourth, the graphic displays make revisions to existing structural designs much easier and less error-prone by eliminating the need for the user to conceptualize the text input as geometry. The user can work directly with the geometry displayed on the screen. Finally, MacPASCO automatically generates the correct PASCO input from the geometry displayed on the screen. This input file can be used with any machine version of PASCO to actually perform the analysis and sizing and to output results. The DEC VAX version of PASCO is written in FORTRAN IV for batch execution and has been implemented on a DEC VAX series computer. The Macintosh version of PASCO was developed for Macintosh II series computers with at least 2Mb of RAM running MPW Pascal 3.0 and Language Systems FORTRAN 2.0 under the MPW programming environment. It includes MPW compatible makefiles for compiling the source code. The Macintosh version uses input files compatible with versions of PASCO running on different platforms. MacPASCO is written in Macintosh Programmers Workbench 3.0, MPW Pascal 3.0, and MacAPP 2.0. The Pascal source code is included on the distribution diskette. MacAPP is a development library which is not included. MacPASCO requires a Mac Plus, SE/30, or MacII, IIx, IIcx, IIci, or IIfx running System 6.0 or greater. MacPASCO is System 7.0 compatible. A minimum of 2Mb of RAM is required for execution. The Macintosh version of PASCO is distributed on four 3.5 inch 800K Macintosh format diskettes. The DEC VAX version is distributed on a 9-track 1600 BPI magnetic tape. The PASCO program was developed in 1981, adapted to the DEC VAX in 1983 and to the Macintosh in 1991. MacPASCO was released in 1992.

Beschreibung: BUCKO is a computer program developed to predict the buckling load of a rectangular compression-loaded orthotropic plate with a centrally located cutout. The plate is assumed to be a balanced, symmetric laminate of uniform thickness. The cutout shape can be elliptical, circular, rectangular, or square. The BUCKO package includes sample data that demonstrates the essence of the program and its ease of usage. BUCKO uses an approximate one-dimensional formulation of the classical two-dimensional buckling problem following the Kantorovich method. The boundary conditions are considered to be simply supported unloaded edges and either clamped or simply supported loaded edges. The plate is loaded in uniaxial compression by either uniformly displacing or uniformly stressing two opposite edges of the plate. The BUCKO analysis consists of two parts: calculation of the inplane stress distribution prior to buckling, and calculation of the plate axial load and displacement at buckling. User input includes plate planform and cutout geometry, plate membrane and bending stiffnesses, finite difference parameters, boundary condition data, and loading data. Results generated by BUCKO are the prebuckling strain energy, inplane stress resultants, buckling mode shape, critical end shortening, and average axial and transverse strains at buckling. BUCKO is written in FORTRAN V for batch execution and has been implemented on a CDC CYBER 170 series computer operating under NOS with a central memory requirement of approximately 343K of 60 bit words. This program was developed in 1984 and was last updated in 1990.

Beschreibung: Beta Lyr is an enigmatic object. Despite several decades of concentrated observational and theoretical effort, our understanding of the system is still rather poor. We discuss two existing structural models of the system, the massive torus model of Wilson (1981), and the low mass accretion disk model of Hubeny & Plavec (1991). In particular, we answer recent criticism of the low mass disk model expressed by the proponents of the massive disk hypothesis. We show that although both theories have a large degree of internal consistency, there are several serious physical objections against the massive disk model. In particular, such a model requires unrealistically low viscosity (large Reynolds number). Moreover, such massive disk is likely to be dynamically unstable. Finally, we propose several observational and theoretical approaches that could lead to construction of a more physically realistic model of the Beta Lyr system.

Beschreibung: Charge transfer rate coefficients are estimated using Landau-Zener and modified Demkov approximations. The coefficients, augmented by those available from the literature, are used in statistical equilibrium equations describing the state of helium-rich supernova plasma. Such a plasma may describe both Type Ib and Type Ic supernova ejecta. The hypothesis that extensive mixing of metals with helium in Type Ic supernovae may provide a catalyst for rapid charge transfer that weakens the helium line emission by altering the excitation balance is tested. It is shown that charge transfer as a mechanism for suppressing helium line emission is ineffective unless the metal abundance is comparable to or larger than the helium abundance. This result supports an earlier conclusion that Type Ic supernovae must be helium poor relative to Type Ib events.

Beschreibung: The instability of the line-driven winds of hot stars leads to the formation of strong shocks. These shocks not only emit thermal X-rays, but also accelerate a small fraction of the thermal electrons and ions to relativistic energies. Synchrotron radiation from these energetic particles can account for the non-thermal radio emission observed from some hot stars, and can also explain the hard X-rays detected in the Einstein X-ray spectra. Our calculations indicate that the gamma-ray emission from non-thermal particles should be detectable by Gamma Ray Observatory (GRO). The detection (or non-detection) of these emissions over a wide energy range, from the radio to gamma-rays, should provide a great deal of information on the structure of the unstable winds and the physics of particle acceleration by shocks.

Beschreibung: We present observations obtained with the Goddard High Resolution spectrograph (GHRS) of three O-type stars: xi per (O7 III), lambda Cep (O6 Iaf) and 10 Lac (O9 V). These observations show evidence of instabilities in the winds of all three stars, although not necessarily those expected from current models of structured winds.

Beschreibung: A thermoviscoplastic finite element method employing the Bodner-Partom constitutve model is used to investigate the response of simplified thermal-structural models to intense local heating. The computational method formulates the problem in rate and advances the solution in time by numerical integration. The thermoviscoplastic response of simplified structures with prescribed temperatures is investigated. With rapid rises of temperature, the nickel alloy structures display initially higher yield stresses due to strain rate effects. As temperatures approach elevated values, yield stress and stiffness degrade rapidly and pronounced plastic deformation occurs.

Beschreibung: We present ultraviolet, optical and infrared observations of the RV Tauri star SX Cen. From a fit of model atmospheres to the data we conclude that Z/Z solar = 0.033, that SX Cen has an extended atmosphere, and that the spectral types implied by the ultraviolet data close to deep minimum are consistent with those implied from optical spectra. Unlike AC Her, there seems to be no silicate component in the circumstellar dust shell, and the star pulsates in purely radial modes. The implication is that, despite their many similarities, SX Cen seems to be at a significantly earlier phase of post-AGB evolution than AC Her. Additional data, particularly in the infrared, are necessary to confirm this conclusion.

Beschreibung: The impact of each fragment of comet SL-9 will produce a downward-propagating pressure wave which will travel at the sound speed through the jovian interior. Since the sound speed increases with depth, most of the energy in the pressure pulse will be strongly refracted and return to the surface, as recently computed by Marley (1994). This wave may in principle be observable as it propagates into the stratosphere, using sufficiently sensitive thermal infrared imaging. If so, it will provide a unique opportunity to constrain models of the jovian interior. This paper extends Marley's calculations to include the effect of the limited spatial resolution which will be characteristic of real observations. The wave pattern on the disk will consist of closely spaced regions of alternating temperature increases and decreases. Spatial averaging will significantly reduce the observed amplitude for resolutions attainable using earth-based telescopes, but the waves should remain above the detection limit.

Beschreibung: Three-dimensional numerical simulations of the impact of Comet Shoemaker - Levy 9 on Jupiter and the resulting vapor plume expansion were conducted using the Smoothed Particle Hydrodynamics (SPH) method. An icy body with a diameter of 2 km can penetrate to an altitude of -350 km (0 km = 1 bar) and most of the incident kinetic energy is transferred to the atmosphere between -100 km to -250 km. This energy is converted to potential energy of the resulting gas plume. The unconfined plume expands vertically and has a peak radiative power approximately equal to the total radiation from Jupiter's disk. The plume rises a few tens of atmospheric scale heights in approximately 10(exp 2) seconds. The rising plume reaches the altitude of approximately 3000 km, but no atmospheric gas is accelerated to the escape velocity (approximately 60 km/s).

Beschreibung: We describe a semiempirical methodology-based on measurements of far-infrared (FIR) lines-that yields information on electron densities in regions where various ionic species exist, effective temperatures (T(sub eff)) for stars ionizing H II regions, and gas-phase heavy element abundances. Although this capability has long been available via optical data, the special features of FIR lines-relative insensitivity to extinction and electron temperature variations-extend the analysis ability. Several line ratios serve as diagnostics of electron density, N(sub e), probing different ionization conditions and different density regimes. The more N(sub e)-diagnostic observations made, the more reliable will be the deciphering of the actual variation in density throughout a nebula. A method to estimate T(sub eff) from the FIR (N III)/(N II) line ratio requires that the nebula be ionization bounded and that substantially all of the flux from the revevant lines be observed. However, to estimate T(sub eff) by a second method that uses the ratio of FIR (S III)/(O III) lines, an ionization-bounded nebula is a sufficient, but not necessary, condition. These restrictions are unnecessary for estimating densities and heavy element abundances. We show that a fairly general determination of metallicity, via the S/H ratio, may be made for H II regions with observations of just two lines-(S III) 19 micron and a hydrogen recombination line (or appropriate substitute). These techniques are applied to recent FIR data for the G333.6-0.2 H II region, including application to the recently measured (N II) 122 and 205 micron lines.

Beschreibung: Ultraviolet observations of the gamma-ray balzar 3C 279 were carried out in 1991 July with the International Ultraviolet Explorer (IUE) satellite, 28 days after the outburst of intense gamma-ray emission detected from this source with the high-energy Energetic Gamma Ray Experiment Telescope (EGRET) instrument aboard the Compton Gamma-Ray Observatory. IUE observations were conducted over the wavelength range 1200-3200 A (5-10 eV) and are compared with archival UV measurements spanning the period 1988-1991. This set of observations was analyzed with a uniform, standardized data-reduction procedure. No significant variability in the UV spectrum of 3C 279 is noted over timescales of hours, though variations at the 3 sigma level or higher exist on longer timescales of months to years. The UV observations which most closely bracket the gamma-ray flare detected by EGRET show approximately 3 sigma differences in flux and spectral index, consistent with (but not definitive proof of) variations in the UV flux period of the gamma-ray flare. Statistically marginal evidence for a monotonic correlation is found between UV spectral index and the emission observed from 3C 279, in that a harder spectrum is associated with increased UV flux.

Beschreibung: The cosmic microwave background radiation (CMBR) spectrum measured by the Far-Infrared Absolute Spectrophotometer (FIRAS) instrument on NASA's Cosmic Background Explorer (COBE) is indistinguishable from a blackbody, implying stringent limits on energy release in the early universe later than the time t = 1 yr after the big bang. We compare the FIRAS data to previous precise measurements of the cosmic microwave background spectrum and find a reasonable agreement. We discuss the implications of the absolute value of y is less than 2.5 x 10(exp -5) and the absolute value of mu is less than 3.3 x 10(exp -4) 95% confidence limits found by Mather et al. (1994) on many processes occurring after t = 1 yr, such as explosive structure formation, reionization, and dissipation of small-scale density perturbations. We place limits on models with dust plus Population III stars, or evolving populations of IR galaxies, by directly comparing the Mather et al. spectrum to the model predictions.

Beschreibung: The cosmic microwave background radiation (CMBR) has a blackbody spectrum within 3.4 x 10(exp -8) ergs/sq cm/s/sr cm over the frequency range from 2 to 20/cm (5-0.5 mm). These measurements, derived from the Far-Infrared Absolute Spectrophotomer (FIRAS) instrument on the Cosmic Background Explorer (COBE) satellite, imply stringent limits on energy release in the early universe after t approximately 1 year and redshift z approximately 3 x 10(exp 6). The deviations are less than 0.30% of the peak brightness, with an rms value of 0.01%, and the dimensionless cosmological distortion parameters are limited to the absolute value of y is less than 2.5 x 10(exp -5) and the absolute value of mu is less than 3.3 x 10(exp -4) (95% confidence level). The temperature of the CMBR is 2.726 +/- 0.010 K (95% confidence level systematic).

Beschreibung: The Far-Infrared Absolute Spectrophotometer (FIRAS) instrument on the Cosmic Background Explorer (COBE) satellite was designed to accurately measure the spectrum of the cosmic microwave background radiation (CMBR) in the frequency range 1-95/cm with an angular resolution of 7 deg. We describe the calibration of this instrument, including the method of obtaining calibration data, reduction of data, the instrument model, fitting the model to the calibration data, and application of the resulting model solution to sky observations. The instrument model fits well for calibration data that resemble sky condition. The method of propagating detector noise through the calibration process to yield a covariance matrix of the calibrated sky data is described. The final uncertainties are variable both in frequency and position, but for a typical calibrated sky 2.6 deg square pixel and 0.7/cm spectral element the random detector noise limit is of order of a few times 10(exp -7) ergs/sq cm/s/sr cm for 2-20/cm, and the difference between the sky and the best-fit cosmic blackbody can be measured with a gain uncertainty of less than 3%.

Beschreibung: We present self-consistent solutions of the full set of ideal MHD equations which describe steady-state relativistic cold outflows from thin accretion disks. The magnetic field forms a spiral which is anchored in the disk, rotates with it, and accelerates the flow out of the disk plane. The collimation at large distances depends on the total amount of electric current that flows along the jet. We considered various distributions of electric current and derived the result that in straight jets which extend to infinite distances, a strong electric current flows along their axis of symmetry. The asymptotic flow velocities are of the order of the initial rotational velocity at the base of the flow (a few tenths of the speed of light). The solutions are applied to both galactic (small-scale) and extragalactic (large-scale) jets.

Beschreibung: Improved electric dipole transition matrix elements for rovibrational transitions in the ground state X(1)Sigma(+) of (12)C(6)O and (13)C(16)O are calculated 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. 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. These convenient to use polynomial representations are given in Tables 1-4 for (12)C(16)O and in Tables 5-8 for (13)C(16)O. We observe that there is intensity enhancement due to vibration-rotation interaction for the P-branch transitions at the expense of the R-branch transitions for delta v = +1. This enhancement can be as large as 40% at the highest J. For the delta v = +2 and +3 transitions, the R-branch transitions are enhanced by as much as a factor of 2.75 and 10 at the highest J, respectively. The P-branch transitions exhibit only minor decreases. Comparisons with previous calculations show good agreement for the delta v = +1 transitions. The comparison for delta v = +2 and +3 transitions show differences as large as a factor of 5.

Beschreibung: We examine the importance of general relativistic corrections to the production of gamma rays near the surface of a neutron star. Due to the change in the magnetic dipole field in curved spacetime, the polar cap angle decreases by 30% compared to flat spacetime. However, the curved photon trajectories compensate for the decrease in the polar cap angle, and, as a result, the pulse profile of the photons emitted parallel to the field is expected to be very similar in curved spacetime and in flat spacetime. We find that the curved spacetime metric significantly increases the magnitude of the magnetic field and, therefore, the attenuation coefficients of curvature radiation gamma rays for pair production in a magnetic field can be increased by factors as large as 100. As a resutl the survival distance of 1 GeV photons for pair production is decreased by a factor of 2 for B is approximately 10(exp 12) G.

Beschreibung: We present a V-magnitude light curve of the cataclysmic binary V794 Quilae covering an 800 day time span. The system shows variations of ups to 3 mag. In particular, there are two dips of approximately 1 mag which last for about 50 days, and are followed by a rapid recovery to the original brightness, and a third dip of approximately 3 mag lasting about 100 days. These fluctuations are thought to be caused by the response of the accretion disk in the system to the cessation of mass transfer from the mass-losing secondary star. We present computations using a time-dependent accretion disk code to delineate a range of allowed behavior for th accretion disk. To model the observed light curve, we require that the model parameters alpha(sub cold) and alpha(sub hot), which characterize the degree of coupling of the viscous dissipation to the orbital shear in the low and high states of the accretion disk, be smaller than they have been inferred to be in the dwarf novae, and that the ratio alpha(sub hot)/alpha(sub cold) also be smaller. The fact that alpha seems to vary with the rate of mass transfer may provide support for the model of Vishniac and Diamond, in which the impact of the mass transfer stream onto the outer edge of the accretion disk excites internal waves which transport angular momentum outward and provide the viscous dissipation.

Beschreibung: A three-bay, space, cantilever truss is probabilistically evaluated to describe progressive buckling and truss collapse in view of the numerous uncertainties associated with the structural, material, and load variables that describe the truss. Initially, the truss is deterministically analyzed for member forces, and members in which the axial force exceeds the Euler buckling load are identified. These members are then discretized with several intermediate nodes, and a probabilistic buckling analysis is performed on the truss to obtain its probabilistic buckling loads and the respective mode shapes. Furthermore, sensitivities associated with the uncertainties in the primitive variables are investigated, margin of safety values for the truss are determined, and truss end node displacements are noted. These steps are repeated by sequentially removing buckled members until onset of truss collapse is reached. Results show that this procedure yields an optimum truss configuration for a given loading and for a specified reliability.