ALBERT

Description: Communications Research Laboratory (CRL) carried out a VLBI (Very Long Baseline Interferometer) experiment of the Pulsar 0329 plus 54 on November in 1992 between Kashima 26mPhi and Usuda 64mPhi. We got a correlation for it by using K-3 VLBI correlator which was developed by CRL. Though we observed a slow pulsar this time, we are going to make millisecond pulsar experiments by making good use of this result and by using K-4 correlator which is under developing in CRL. Its result will give us precise positions and proper motions of pulsars which are useful for obtaining precise pulsar timing.

Description: Basic algorithms for unstructured mesh generation and fluid flow calculation are discussed. In particular the following are addressed: preliminaries of graphs and meshes; duality and data structures; basic graph operations important in CFD (Computational Fluid Dynamics); triangulation methods, including Varonoi diagrams and Delaunay triangulation; maximum principle analysis; finite volume schemes for scalar conservation law equations; finite volume schemes for the Euler and Navier-Stokes equations; and convergence acceleration for steady state calculations.

Description: Transient solutions were obtained for a square region of heat conducting semitransparent material cooling by thermal radiation. The region is in a vacuum environment, so energy is dissipated only by radiation from within the medium leaving through its boundaries. The effect of heat conduction during the transient is to partially equalize the internal temperature distribution. As the optical thickness of the region is increased, the temperature gradients increase near the boundaries and corners, unless heat conduction is large. The solution procedure must provide accurate temperature distributions in these regions to prevent error in the calculated radiation losses. Two-dimensional numerical Gaussian integration is used to obtain the local radiative source term. A finite difference procedure with variable space and time increments is used to solve the transient energy equation. Variable spacing was used to concentrate grid points in regions with large temperature gradients.

Description: Very Long Baseline Interferometry (VLBI) observations of extragalactic radio sources provide the basis for defining an accurate non-rotating reference frame in terms of angular positions of the sources. Measurements of the distance from the Earth to the Moon and to the inner planets provide the basis for defining an inertial planetary ephemeris reference frame. The relative orientation, or frame tie, between these two reference frames is of interest for combining Earth orientation measurements, for comparing Earth orientation results with theories referred to the mean equator and equinox, and for determining the positions of the planets with respect to the extragalactic reference frame. This work presents an indirect determination of the extragalactic-planetary frame tie from a combined reduction of VLBI and Lunar Laser Ranging (LLR) observations. For this determination, data acquired by LLR tracking stations since 1969 have been analyzed and combined with 14 years of VLBI data acquired by NASA's Deep Space Network since 1978. The frame tie derived from this joint analysis, with an accuracy of 0.003 sec, is the most accurate determination obtained so far. This result, combined with a determination of the mean ecliptic (defined in the rotating sense), shows that the mean equinox of epoch J2000 is offset from the x-axis of the extragalactic frame adopted by the International Earth Rotation Service for astrometric and geodetic applications by 0.078 sec +/- 0.010 sec along the y-direction and y 0.019 sec +/- 0.001 sec. along the z-direction.

Description: Bolometer observations at 250 GHz of fifteen minor planets have shown that the emissivity of these objects is close to unity. This results in an independent method to determine the absolute calibration scale of radio observations at mm wavelengths: Applying our results to Mars, the prime calibrator at this wavelength, gives a mean absolute disk temperature at mean solar distance of approximately 210 K. Further, the diameters of circularly symmetric asteroids can be determined or the surface area of asteroids can be estimated assuming some geometric constraints on their shape.

Description: Using Ulysses radio wave data taken during the 1992 Jupiter encounter, we conclude that there are significant large and small spatial scale azimuthal asymmetries at high latitudes in the Io plasma torus. During a period of time near perijove when the spacecraft motion was predominantly in the azimuthal direction and was relatively fixed in both latitude and radial distance, inferred electron densities depart significantly from the common assumption of longitudinal symmetry. Specifically, electron plasma concentrations near 0 deg system III longitude (and 0400 LT) are greater than those near 180 deg (and 0000 LT). Superposed on this large-scale variation are regularly spaced density depletions, 30-50% in magnitude, and having a spatial periodicity of about 17 deg. Some of these depletions may drive various known radio and plasma wave sources by means of large B parallel electric potentials. The observations are compared with recent models and with the in-situ Voyager observations.

Description: It is shown that to satisfy the general accepted compressible law of the wall derived from the Van Driest transformation, turbulence modeling coefficients must actually be functions of density gradients. The transformed velocity profiles obtained by using standard turbulence model constants have too small a value of the effective von Karman constant kappa in the log-law region (inner layer). Thus, if the model is otherwise accurate, the wake component is overpredicted and the predicted skin friction is lower than the expected value.

Description: The Be/X-ray binary system A 118-616 has been observed to undergo a major outburst in January 1992 by the Burst and Transient Source Experiment (BATSE) all-sky monitor on the Compton Gamma Ray Observatory. Ground-based optical and IR observations, supported by UV observations obtained under an IUE Target of Opportunity program, have provided us with an excellent multiwavelength study of this system to complement the X-ray data set. The results from this campaign are presented showing the details of the X-ray timing studies, the very strong H-alpha emission and the bright IR excess from the Be star's circumstellar disk (the fuel for the accretion process). Implications for the physics of the system are discussed.

Description: Observations of Einstein Large Magellanic Cloud (LMC) X-ray point sources have been made with ROSAT's High-Resolution Imager to obtain accurate positions from which to search for optical counterparts. This paper is the first in a series reporting results of the ROSAT observations and subsequent optical observations. It includes the X-ray positions and fluxes, information about variability, optical finding charts for each source, a list of identified counterparts, and information about candidates which have been observed spectroscopically in each of the fields. Sixteen point sources were measured at a greater than 3 sigma level, while 15 other sources were either extended or less significant detections. About 50% of the sources are serendipitous detections (not found in previous surveys). More than half of the X-ray sources are variable. Sixteen of the sources have been optically identified or confirmed: six with foreground cool stars, four with Seyfert galaxies, two with signal-to-noise ratio (SNR) in the LMC, and four with peculiar hot LMC stars. Presumably the latter are all binaries, although only one (CAL 83) has been previously studied in detail.

Description: Imaging spectrometer observations were made of the surface of the Moon during the December 1990 flyby of the Earth-Moon system by the Galileo spacecraft. This article documents this data set and presents analyses of some of the data. The near infrared mapping spectrometer (NIMS) investigation obtained 17 separate mosaics of the Moon in 408 spectral channels between about 0.7 and 5.2 micrometers. The instrument was originally designed to operate in orbit about Jupiter and therefore saturates at many spectral channels for most measurement situations at 1 AU. However, sufficient measurements were made of the Moon to verify the proper operation of the instrument and to demonstrate its capabilities. Analysis of these data show that the NIMS worked as expected and produced measurements consistent with previous ground-based telescopic studies. These are the first imaging spectrometer measurements of this type from space for the Moon, and they illustrate several major points concerning this type of observation and about the NIMS capabilities specifically. Of major importance are the difference between framing and scanning instruments and the effects of the spacecraft and the scan platform on the performance of such and experiment. The science return of subsequent NIMS and other investigation measurements will be significantly enhanced by the experience and results gained.

Description: Spectral reflectances and geometic albedos between 2300 and 3250 A are determined for 45 asteroids from data acquired by the International Ultraviolet Explorer satellite. The geometric albedos are consistently low, ranging from approximately 0.02 for C-type asteroids to approximately 0.08 for M-type asteroids. An exception is the single E-type asteroid (44 Nysa) with a geometric albedo of 0.3 at 2950 A. We find that the three major asteroid taxonomic classes persist into the UV. The taxonomic classes are distinguished primarily by their albedos, but S types are generally redder than C or M types. The first ultraviolet phase curved of asteroids are presented.

Description: The multigrid method has been applied to an existing three-dimensional compressible Euler solver to accelerate the convergence of the implicit symmetric relaxation scheme. This lower-upper symmetric Gauss-Seidel implicit scheme is shown to be an effective multigrid driver in three dimensions. A grid refinement study is performed including the effects of large cell aspect ratio meshes. Performance figures of the present multigrid code on Cray computers including the new C90 are presented. A reduction of three orders of magnitude in the residual for a three-dimensional transonic inviscid flow using 920 k grid points is obtained in less than 4 min on a Cray C90.

Description: The present paper explores the use of large-eddy simulations as a tool for predicting noise from first principles. A high-order numerical scheme is used to perform large-eddy simulations of a supersonic jet flow with emphasis on capturing the time-dependent flow structure representating the sound source. The wavelike nature of this structure under random inflow disturbances is demonstrated. This wavelike structure is then enhanced by taking the inflow disturbances to be purely harmonic. Application of Lighthill's theory to calculate the far-field noise, with the sound source obtained from the calculated time-dependent near field, is demonstrated. Alternative approaches to coupling the near-field sound source to the far-field sound are discussed.

Description: The Ultraviolet Imaging Telescope (UIT) cameras include two-stage magnetically focused image intensifiers that introduce small but significant geometrical distortions into the data. These distortions, which create positional offsets as large as 25 arcsec at the field edges, are corrected by the procedure described here to 2-3 arcsec, approximately the resolution of the images. The distortion is measured by comparing and correcting UIT images to digitized Guidestar survey plates of the same fields. Two-dimensional third-order polynomials are used to model the distortion. The models assume that the distortion is an instrumental effect, independent of mission elapsed time and target, and that the effect of distortion is an instrumental effect, independent of mission elapsed time and target, and that the effect of distortion in the center of each field is minimal. The models are used to improve computed astrometric plate solutions and to remove the geometric distortion while transforming the image to a standard north-up, ease-left orientation.

Description: Photoelectric observations of the asteroid 55 Pandora were carried out in Feb. - Mar. 1989 (6 nights) in range of phase angles Delta alpha = (2.5 - 1.4) deg, and in Sep. - Nov. 1991 (15 nights) Delta alpha = (0.5 - 16.3) deg. Average amplitudes of lightcurves in these oppositions are 0.22 m and 0.10 m, respectively. The value of linear phase coefficients and the absence of spike-effect of alpha less than 2 deg indicate that Pandora is a typical M-asteroid and the high albedo measured by IRAS-satellite is not real.

Description: Delta Orionis A was observed with ROSAT in February 1991 at ten different phases of the 5.7 day binary orbit. The soft X-ray emission shows little dependence on binary phase favoring models for the production of the soft X-rays in the outer parts of the stellar wind. The energy spectrum obtained with the Position Sensitive Proportional Counter (PSPC) shows two thermal components with temperatures of 0.1 and 0.2 keV. The temperature and absorption of the hotter component is consistent with the Einstein solid state spectrometers (SSS) spectrum obtained in 1979. The two temperatures suggest either two emission regions with different temperatures, or one emission region with a temperature distribution where the hotter parts are closer to the OB stars than the cooler parts as indicated by the higher absorption of the hotter spectral component. This would be in contrast to stellar wind shock models which predict shock temperatures increasing with distance to the star.

Description: Although Am stars have been assumed to be nonmagnetic stars, we found recently that the hot Am stars omicron Peg has a magnetic field of the order of 2 kG, with a complex structure. Therefore, we reconsidered the question of magnetism among upper main sequence chemically peculiar stars. We report here new observations of a pair of magnetically sensitive Fe II lines in the red spectrum of Am stars and of HgMn stars. Two other hot Am stars, HD 29173 and HD 195479A, are likely to possess also a magnetic field of about the same strength. It is striking that 3 out of 4 hot Am stars observed up to now could be magnetic. On the other hand, this fast method is limited by blends and cannot provide a satisfactory diagnosis for classical Am stars and for HgMn stars. A mutli-line approach must be the next step towards our understanding of magnetism among these chemically peculiar stars.

Description: The steady state solution of the system of equations consisting of the full Navier-Stokes equations and two turbulence equations has been obtained using a multigrid strategy of unstructured meshes. The flow equations and turbulence equations are solved in a loosely coupled manner. The flow equations are advanced in time using a multistage Runge-Kutta time-stepping scheme with a stability-bound local time step, while turbulence equations are advanced in a point-implicit scheme with a time step which guarantees stability and positivity. Low-Reynolds-number modifications to the original two-equation model are incorporated in a manner which results in well-behaved equations for arbitrarily small wall distances. A variety of aerodynamic flows are solved, initializing all quantities with uniform freestream values. Rapid and uniform convergence rates for the flow and turbulence equations are observed.

Description: Analysis of data from the Spectroscopy Detectors (SDs) of the Burst and Transient Source Experiment (BATSE) on the Compton Gamma Ray Observatory (GRO) has found no convincing line features in the spectra of gamma-ray bursts (GRBs) in almost 3 years of operation, in contrast to expectations based on results from other experiments. In this Letter we discuss the visual search for narrow lines in the SD data. The search has examined 192 bursts, of which approximately 18 were intense enough that lines similar to those seen by instruments on the Ginga satellite would have been visible between approximately 20 and approximately 100 keV. A simplified calculation shows that the BATSE and Ginga results are consistent at the 13% level.

Description: We report the first observations of the 13.5 micron fundamental band of SiS in the spectrum of the heavily obscured carbon star IRC +10216. The lines are formed in the inner region of the circumstellar envelope where the gas is accerlerating and where the temperature ranges from 800-500 K. We have carried out a detailed model of the observed line profiles. Our observations are best fit by a gradient in the abundance of SiS. We derive an abundance relative to molecular hydrogen of x(SiS) = 4.3 x 10(exp -6) at a distance of twelve stellar radii from the central star rising to x(SiS) = 4.3 x 10(exp -5) at a few stellar radii from the surface of the star.

Description: I present the ROSAT Position-Sensitive Proportional Counter (PSPC) image of the nearby, nearly face-on, spiral galaxy NGC 6946. The galaxy was observed with Einstein in 1980, uncovering a nuclear source, a source on or near the northern spiral arm, SN 1980K, and evidence for non-point-source (i.e., diffuse) emission. The ROSAT image resolves the nuclear region into approximately three sources, yields an accurate position for the north arm source, and reveals diffuse extended emission across the galaxy face, tracing at least the bright northern spiral arm of the galaxy. The diffuse emission is almost certainly the very hot component of the galaxy's interstellar medium and is probably similar to that found in the Milky Way and the Large Magellanic Cloud.

Description: The first two years of Cosmic Background Explorer (COBE) Differential Microwave Radiometers (DMR) observations of the cosmic microwave background (CMB) anisotropy are analyzed and compared with our previously published first year results. The results are consistent, but the addition of the second year of data increases the precision and accuracy detected CMB temperature fluctuations. The 2 yr 53 GHz data are characterized by rms temperature fluctuations of (delta-T)(sub rms) (7 deg) = 44 +/- 7 micro-K and (delta-T)(sub rms) (10 deg) = 30.5 +/- 2.7 micro-K at 7 deg and 10 deg angular resolution, respectively. The 53 x 90 GHz cross-correlation amplitude at zero lag is C(0)(sup 1/2) = 36 +/- 5 micro-K (68% CL) for the unsmoothed (7 deg resolution) DMR data. We perform a likelihood analysis of the cross-correlation function, with Monte Carlo simulations to infer biases of the method, for a power-law model of initial density fluctuations, P(k) proportional to R(exp n). The Monte Carlo simulations indicate that derived estimates of n are biased by +0.11 +/- 0.01, while the subset of simulations with a low quadrupole (as observed) indicate a bias of +0.31+/- 0.04. Derived values for 68% confidence intervals are given corrected (and not corrected) for our estimated biases. Including the quadrupole anisotropy, the most likely quadrupole-normalized amplitude is Q(sub rms-PS) = 14.3(sup + 5.2 sub -3.3) micro-K (12.8(sup + 5.2 sub -3.3) micro-K0 with a spectral index n = 1.42(sup + 0.49 sub -0.55)(n = 1.53(sup + 0.49 sub -0.55). With n fixed to 1.0 the most likely amplitude is 18.2 +/- 11.5 micro-K (17.4 +/- 1.5 micro-K). The marginal likelihood of n is 1.42 +/- 0.37 (1.53 +/- 0.37). Excluding the quadrupole anisotropy, the most likely quadrupole-normalized amplitude is Q(sub rms-PS) = 17.4(sup + 7.5 sub -5.2) micro-K (15.8(sup + 7.5 sub -5.2) micro-K) with a spectral index n = 1.11(sup + 0.60 sub -0.55) (n = 1.22(sup + 0.60 sub -0.55). With n fixed to 1.0 the most likely amplitude is 18.6 +/- 1.6 micro-K (18.2 +/- 1.6 micro-K). The marginal likelihood of n is 1.11 +/- 0.40 (1.22 +/- 0.40). Our best estimate of the dipole from the 2 yr DMR data is 3.363 +/- 0.024 mK toward Galactic coordinates (l, b) = (264.4 deg +/- 0.2 deg, 48.1 deg +/- 0.4 deg), and our best estimate of the rms quadrupole amplitude in our sky is 6 +/- 3 micro-K (68% CL).

Description: The Global Oscillation Network Group (GONG) Project will place a network of instruments around the world to observe solar oscillations as continuously as possible for three years. The Project has now chosen the six network sites based on analysis of survey data from fifteen sites around the world. The chosen sites are: Big Bear Solar Observatory, California; Mauna Loa Solar Observatory, Hawaii; Learmonth Solar Observatory, Australia; Udaipur Solar Observatory, India; Observatorio del Teide, Tenerife; and Cerro Tololo Interamerican Observatory, Chile. Total solar intensity at each site yields information on local cloud cover, extinction coefficient, and transparency fluctuations. In addition, the performance of 192 reasonable networks assembled from the individual site records is compared using a statistical principal components analysis. An accompanying paper descibes the analysis methods in detail; here we present the results of both the network and individual site analyses. The selected network has a duty cycle of 93.3%, in good agreement with numerical simulations. The power spectrum of the network observing window shows a first diurnal sidelobe height of 3 x 10(exp -4) with respect to the central component, an improvement of a factor of 1300 over a single site. The background level of the network spectrum is lower by a factor of 50 compared to a single-site spectrum.

Description: We report an experimental absolute oscillator strengths for 18 UV lines of Ru II, obtained by combining laser-induced flourescence measurements of radiative lifetimes and branching fractions from line intensities in a calibrated Fourier-transform spectrum Hubble Space Telescope/Goddard High Resolution Spectrograph (HST/GHRS) observations of the spectrum of the sharp-lined B star chi Lupi contain six of these lines, for which 'astrophysical' relative f-values have been determined. The agreement is within 0.10 dex for a Ru abundance of log N(Ru)/N(H) = -7.90, which is 2.3 dex above the solar abundance.

Description: The six-band ultraviolet light curves of beta Lyrae obtained with the Orbiting Astronomical Observatory (OAO) A-2 in 1970 exhibited a very unusual behavior. The secondary minimum deepened at shorter wavelength, indicating that one was not observing light variations caused primarily by the eclipses of two stars having a roughly Planckian energy distribution. It was then suggested that the light variations were caused by a viewing angle effect of an optically thick, ellipsoidal circumbinary gas cloud. Since 1978 beta Lyrae has been observed with the International Ultraviolet Explorer (IUE) satellite. We have constructed ultraviolet light curves from the IUE archival data for comparison with the OAO A-2 results. We find that they are in substantial agreement with each other. The Voyager ultraviolet spectrometer was also used to observe this binary during a period covered by IUE observations. The Voyager results agree with those of the two other satellite observatories at wavelengths longer than about 1350 A. However, in the wavelength region shorter than the Lyman-alpha line at 1216 A, the light curves at 1085 and 965 A show virtually no light variation except an apparent flaring near phase 0.7, which is also in evidence at longer wavelengths. We suggest that the optically thick circumbinary gas cloud, which envelops the two stars completely, assumes a roughly spherical shape when observed at these shorter wavelengths.

Description: From two sets of the spectroscopic observations covering a ten year period we have obtained 59 radial velocities of the chromospherically-active star HD 28591 = V492 Per. It is a G9III single-lined spectroscopic binary with a period of 21.2910 days and a circular orbit. The upsilon sin i of 24.6 km/sec, results in a minimum radius 10.3 solar radii. We estimate a distance of 165 +/- 40 pc and an orbital inclination of 65 +/- 25 degrees. The secondary is probably a mid to late-type K dwarf. The star is brighter than the limiting magnitude of the Bright Star Catalogue. The mean photometric and the orbital periods are identical within their uncertainties. Since the star fills a significant fraction of its Roche lobe, about 62%, the photometric light curve may be the result of starspots and a modest ellipticity effect.

Description: We report preliminary results of an ASCA observation of the Seyfert 2 galaxy Markarian 3 (Mkn 3). Comparison with previous Ginga and Broad-Band X-ray Telescope (BBXRT) observations shows that the observed hard X-ray luminosity above 4 keV decreased by a factor of approximately 3 (intrinsic luminosity by almost a factor of 6) in a period of approximately 3.6 yr. On the other hand, the soft luminosity has not varied significantly in approximately 13 yr, lending support to the extended nature of the soft emission, perhaps dominated by scattering of the nuclear X-rays. ASCA resolves the Fe K line emission into at least two components for the first time. The dominant component at 6.4 keV has an equivalent width of approximately 860 eV and full width at half maximum (FWHM) approximately equals 10(exp 4) km/sec, while the second component has an equivalent width of approximately 190 eV and appears to be narrower than the first. The total intensity of the Fe K emission decreased by factor of over 3 in response to the decrease in the continuum level, implying that a substantial part of the dominant Fe K emission must originate in a region smaller than that responsible for the soft emission. The variability provides direct evidence that the hard X-ray continuum and Fe K line in this Seyfert 2 are being observed directly through the nuclear obscuring material, not in scattered light.

Description: We present preliminary results of an ASCA observation of the classic soft-excess quasar PG 1211+143. The overall ASCA spectrum can be characterized by a blackbody with a temperature of approximately 125 eV (quasar frame) and a power law with photon index of approximately 2. Simultaneous ROSAT data are suggestive of further steepening of the spectrum just below the ASCA band. Comparison with previous observations shows that the soft flux in the 0.1-2 keV band varies by at least a factor of approximately 16, scaling roughly as the square of the hard flux in the 2-10 keV band over a timescale of approximately 13.5 yr. We also find evidence of short-term amplitude variability of up to a factor of approximately 2 on a timescale of approximately 2 x 10(exp 4) sec, in both the soft and hard flux so that the soft and hard photons are likely to originate from the same, compact, region. The data rule out variable absorption (cold or ionized) as the origin of the soft excess, favoring an intrinsic emission component. However, we argue against optically thin emission for the 'blue bump' in PG 1211+143. The large amplitude soft X-ray variability may be indicative of variations in the effective temperature, or peak, of the soft component. There is only marginal evidence for Fe K line emission between 6-7 keV in the quasar frame.

Description: The overlapping plates method has been applied to crossing-point Charge Coupled Device (CCD) observations of minor planet 243 Ida to produce absolute position measurements precise to better than 0.1 sec and differential position measurements precise to better than 0.06 sec. Although these observations numbered only 17 out of the 520 that produced the final ground-based Ida ephemeris for the Galileo spacecraft flyby, their inclusion decreased Ida's downtrack error from 78 to 60 km and its out-of-plane error from 58 to 44 km.

Description: In the present work, the high Reynolds number flow past an inclined plate with a splitter plate placed in its wake is considered numerically. A numerical conformal mapping technique is employed to transform the two-plate system into the same number of cylinders: the flow field is assumed to be two-dimensional. The vortex shedding from the inclined plate is modelled using the discrete vortex method. It is shown that the splitter plate has a profound effect on the development of the flow over a range of values of a suitably defined offset parameter and for a range of positions of the leading edge of the splitter plate. The acoustic field is also calculated and the spectrum reflects the flow results.

Description: Among the gamma-ray bursts (GRBs) observed by the Burst and Transient Source Experiment (BATSE) on board the Compton Gamma Ray Observatory we define a subclass of bursts based on similar morphology: a sharp rise followed by a longer decay time. About 7% of all the gamma-ray bursts observed by BATSE fall into this subclass. We study the spectral evolution of these bursts by fitting models to time-segmented burst spectra and find no clear distinction between the spectral evolutionary properties of this subclass and those of other bursts. Further, we study the high time resolution spectral evolution of this subclass of GRBs using their spectral hardness ratios. A majority of the bursts show hardness ratio leading the counting rate and also display a continuous hard to soft evolution. The time lag between the counting rate and the hardness ratio is found to be directly correlated with the rise time of the counting rate profile. We also find, for the first time, evidence for spectral variation in a timescale of 64 ms.

Description: In this paper we discuss a variety of issues concerning the exciting and mysterious Galactic center gamma-ray sources 1E 1740.7-2942 and GRS 1758-258. We discuss the problem associated with the highly uncertain X-ray absorption column toward 1E 1740.7-2942 and use the recent Roentgen Satellite (ROSAT) results to narrow its range to 0.5-1 x 10(exp 23)/sq cm. Then the current upper limits from deep optical and near-IR searches of stellar objects at these source locations are plotted on an H-R diagram, from which we find the mass of a potential companion star of the (supposed) black hole in GRS 1758-258 to be less than 4 solar mass and in 1E 1740.7-2942 to be less than 9 solar mass. The observed well-collimated radio jets in 1E 1740.7-2942 require the existence of a stable accretion disk (presumably from binary accretion). The apparent association of 1E 1740.7-2942 with a high-density molecular cloud, on the other hand, points to possible accretion directly from the interstellar medium (ISM). We present an analysis of the energetics and kinematics of the radio jets in 1E 1740.7-2942. We present the long-term X-ray light curves of the two sources which include both the Granat/SIGMA's 3 yr monitoring data and all the data from previous imaging balloon and satellite observations over the last decade. The possible physical mechanisms responsible for producing both the long-term X-ray variations and the radio jets are postulated. We also consider Roche lobe-overflowing, low-mass X-ray binaries and Bondi-Hoyle accretion directly from a high-density surrounding medium. We propose a plausible scenario in which both sources are binary systems with a black hole primary and a low-mass companion and they are accreting mainly from the ISM at a rate self-regulated by the interaction between the accretion flow and the emerging hard X-ray flux.

Description: In this panel discussion contributions were made by K. Strom, L. Nordh and H. Zinnecker on the contributions of surveys to the study of star formation regions, by B. Burton on a survey of galactic H I and by E. Dwek on the detection of galactic supernovae by infrared surveys. The contributions of K. Strom, L. Nordh and E. Dwek are summarized here.

Description: Data with the Two Micron All Sky Survey (2MASS) prototype camera were obtained in a 2.3 sq. deg region in Taurus containing Heiles Cloud 2, a region known from Infrared Astronomy Satellite (IRAS) observations to contain a number of very young solar type stars. Data at 1.25 (J), 1.65 (H), and 2.2 (K(sub s)) micrometers are presented. These data are representative of the type and quality of data expected from the planned near-IR surveys, 2MASS and Deep Near-Infrared Survey (DENIS). Near-IR surveys will be useful for determining the large scale variation of extinction with clouds, for determining the luminosity function in nearby clouds down to ranges of 0.1-1.0 solar luminosity, and for finding highly extincted T Tauri stars missed by IRAS because the bulk of their luminosity is emitted shortward of 12 micrometers.

Description: Red Giants on the Asymptotic Giant Branch typically are more luminous than M(sub K) = -7 mag. Therefore, a new Two Micron Sky Survey (2MASS) which will go as faint as m(sub K) = 10 mag will be able to observe most of these stars in the Milky Way. Such a complete census will enable us to develop a much better understanding of Galactic Structure. It will be important to separate the luminous red giants into their different subclasses because these subclasses trace different Galactic Populations. For example, Miras with periods less than 300 days can be used to study the 'thick disk', while Miras with periods greater than 300 days belong to the 'thin disk'.

Description: Infrared Astronomy Satellite (IRAS), ISO, Space Infrared Telescope Facility (SIRTF), WIRE, Deep Near-Infrared Survey (DENIS), and Two Micron All-Sky Survey (2MASS) observations were used to compute the maximum number of observable brown dwarfs for various infrared surveys by combining the maximum possible Oort limit (0.1 'missing' solar mass p/cu c) with all possible brown dwarf mass and age distributions. This approach shows what limits will be placed on the contribution of brown dwarfs to any possible 'missing mass' if no brown dwarfs are observed. I consider brown dwarfs with masses of 0.01-0.08 solar mass and ages of 10(exp 9)-10(exp 10) years. The full range of predicted numbers of brown dwarfs above approx. 6 times the noise of each of the below surveys is: IRAS Point Source Catalog, 0.02-6; IRAS Faint Source Catalog absolute value of b greater than 10 deg, 0.05-16; ISO (2 week 12 micrometer survey), 0.15-80; SIRTF (2 week 12 micrometer survey), 2.50-1600; WIRE (4 month 12 micrometer survey), 21.80-6000; DENIS(half sky) absolute value of b greater than 10 deg, 0.00-2000; and 2MASS(full sky) absolute value of b greater than 10 deg, 0.00-8800. A failure to find brown dwarfs in the IRAS FSC would just barely rule out about half of the mass-age range for Oort limit total masses. A failure to find brown dwarfs in 2MASS/DENIS would rule out roughly the same mass-age range, but would set a limit of 0.1-0.01 times the Oort mass in that mass-age region. No limits would be set for the other half of the mass-age range since both IRAS and 2MASS/DENIS have insufficient sensitivity for brown dwarfs with T less than 750 K. A failure to find brown dwarfs with ISO would rule out almost all of the mass-age range for Oort limit total masses, but would not set a significantly lower limit to the brown dwarf mass limit. A failure to find brown dwarfs with SIRTF or WIRE would rule out the entire mass-age range for Oort limit total masses and set an upper limit of 0.1-0.001 times the Oort mass. To date, about 18% of the IRAS FSC has been searched down to 6 sigma, and no brown dwarfs have been found. This sets a 95% upper limit of 3 in 18% of the sky, or 13 in the entire FSC for absolute value b greater than 10 deg. To begin to set useful limits from 2MASS or DENIS, approximately 400 square degrees needs to be analyzed. To date, only a few square degrees of results from the 2MASS prototype camera have been examined, with no brown dwarfs found so far.

Description: Multiple linear regression analysis was used to derive the effective solar flare contributions of each of the McIntosh classification parameters. The best fits to the combined average number of M- and X-class X-ray flares per day were found when the flare contributions were assumed to be multiplicative rather than additive. This suggests that nonlinear processes may amplify the effects of the following different active-region properties encoded in the McIntosh classifications: the length of the sunspot group, the size and shape of the largest spot, and the distribution of spots within the group. Since many of these active-region properties are correlated with magnetic field strengths and fluxes, we suggest that the derived correlations reflect a more fundamental relationship between flare production and the magnetic properties of the region. The derived flare contributions for the individual McIntosh parameters can be used to derive a flare rate for each of the three-parameter McIntosh classes. These derived flare rates can be interpreted as smoothed values that may provide better estimates of an active region's expected flare rate when rare classes are reported or when the multiple observing sites report slightly different classifications.

Description: We report Very Large Array (VLA) A-configuration studies of a sample of 49 radio galaxies at redshift less than 1. These were selected with no prior knowledge of their morphology and were chosen to match the redshift and luminosity distribution of a previously studied sample of radio-loud quasars. We compare the radio galaxies with the quasar sample and also with a sample of 29 radio galaxies selected for steep spectrum and double-lobe structure. We find that the radio galaxies have more luminous lobes and mostly weaker cores, and there is no population of one-sided sources associated with the galaxies. The radio galaxies' lobe length ratios and lobe power ratios differ from quasars. The overall sizes of the two types of sources are similar, but the radio galaxies have a 3 times larger upper envelope. The distribution of bend angles is similar but the radio galaxies have fewer very bent and straight sources. We discuss these and other comparisons in detail and suggest that while quasars appear to be viewed within a cone and radio galaxies outside it, the two types of source also have intrinsic differences, and both have individual growth and evolution scenarios. This is supported by previously observed differences in optical properties between the two source types.

Description: The Astrometric Imaging Telescope (AIT) is designed to probe the circumstellar environment by both direct imaging and indirect astrometric measurements. The Circumstellar Imager (CI) is a coronagraphic camera and is the direct imaging component of the AIT. The CI is designed to obtain high-sensitivity images of the circumstellar region. It provides crucial non-inferential information relating to the frequency, origin, and evolution of planetary systems and all forms of circumstellar matter. Such imaging is usually limited by the scattered and diffracted light halos of the star itself, which are greatly suppressed in the CI by mating a novel high-efficiency coronagraph with a phase-compensated optical system. For faint point sources in the circumstellar region, the CI will have a sensitivity in excess of 5 magnitudes fainter than the as-designed Hubble Space Telescope (HST). Laboratory data are shown for the coronagraph, which, in a diffraction-limited environment, is capable of suppressing the stellar diffraction sidelobes by several orders of magnitude without significant sacrifice of field of view. In order to realize the high rejection levels inherent in the coronagraph design, it is necessary to limit scatter in the optical systems, imposing a mid-spatial frequency figure error requirement an order of magnitude smaller than that of the HST. Experimental data directed toward meeting this requirement are also shown.

Description: The Space Infrared Telescope Facility, to be launched into a near-Earth heliocentric orbit in the year 2001, will open broad new vistas for the study, at infrared wavelengths, of the objects in the Solar System and planetary systems around other stars. This paper focuses on the study of Kuiper-belt comets and circumstellar planetary debris disks.

Description: Narrow-angle astrometry with long-baseline infrared interferometers can provide extremely high accuracies as required for indirect planet detection. Narrow-angle astrometric interferometry exploits the properties of atmospheric turbulence over fields smaller than the interferometer baseline divided by the atmospheric scale height. For such fields, accuracy is linear with star separation, and nearly inversely proportional to baseline length. To exploit these properties, the interferometer observes a relatively bright (less than 13 mag(sub K)) target in the near infrared at 2.2 micrometers, and uses phase referencing to find a reference star within the 2.2-micrometers isoplanatic patch. With this technique faint references can be found for most targets. With baselines greater than 100 m, which also minimize photon-noise errors, and with careful control of systematic errors by using laser metrology, accuracies of tens of microarcseconds/square root of (hour) should be possible.

Description: We have utilized the integral moment analysis technique of Horack & Emslie to extract information on the allowable form of the luminosity function for gamma-ray bursts observed by Burst and Transient Source Experiment (BATSE). Using the general properties of moments, we are able to derive constraints on the range of luminosity from which the gamma-ray bursts must be sampled. These constraints are independent of the form of the radial distribution of the gamma-ray bursts, and depend only on the assumptions that space is Euclidean and that the luminosity function phi(L) is distance independent. For power-law luminosity functions of the form phi(L) = A(sub 0)L(exp -alpha), we find that the range of luminosity from which 80% of the gamma-ray bursts must be sampled cannot exceed approximately 6.5, with a 3 sigma upper limit of 12-15, regardless of the value of alpha.

Description: We have analyzed the existing spectra of seven comets which show an emission feature at 7.8-13 micrometers. Most have been converted to a common calibration, taking into account the SiO feature in late-type standard stars. The spectra are compared with spectra of the Trapezium, interplanetary dust particles (IDPs), laboratory mineral samples, and small particle emission models. The emission spectra show a variety of shapes; there is no unique 'cometary silicate'. A peak at 11.20-11.25 micrometers, indicative of small crystalline olivine particles, is seen in only three comets of this sample, P/Halley, Bradfield 1987 XXIX, and Levy 1990 XX. The widths of the emission features range from 2.6 to 4.1 micrometers (FWHM). To explain the differing widths and the broad 9.8 micrometers maximum, glassy silicate particles, including both pyroxene and olivine compositions, are the most plausible candidates. Calculations of emission models confirm that small grains of glassy silicate well mixed with carbonaceous material are plausible cometary constituents. No single class of chondritic aggregate IDPs exhibits spectra closely matching the comet spectra. A mixture of IDP spectra, particularly the glass-rich aggregates, approximately matches the spectra of comets P/Halley, Levy, and Bradfield 1987 XXIX. Yet, if comets are simply a mix of IDP types, it is puzzling that the classes of IDPs are so distinct. None of the comet spectra match the spectrum of the Trapezium. Thus, the mineralogy of the cometary silicates is not the same as that of the interstellar medium. The presence of a component of crystalline silicates in comets may be evidence of mixing between high- and low-temperature regions in the solar nebula.

Description: Numerical results obtained with direct simulation Monte Carlo and Navier-Stokes methods are presented for a Mach-20 nitrogen flow about a 70-deg blunted cone. The flow conditions simuulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how the wake structure changes as a function of rarefaction, what the afterbody levels of heating are, and to what limits the continuum models are realistic as rarefaction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the after body sting are emphasizes in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

Description: The redetermination of the luni-solar precession is performed by using Very Long Baseline Interferometry (VLBI) catalogues of extragalactic radio sources containing positions at observation epochs which cover 9 years. The positions of the sources have been determined by the VLBI network of the Crustal Dynamics Project of NASA from measurements covering consecutive half year and one year intervals between 1981 and 1989. In the course of their reduction the International Astronomical Union (IAU) recommended precession terms were applied. On the assumption of an imperfection of the luni-solar precession it is expected that the positions of each source obtained at different epochs but uniformly reduced to the reference system J2000.0 exhibit an apparent proper motion. We attributed this motion tentatively to precession and solved for a correction of the luni-solar precession by making a least squares fit to the position differences between homonymous sources in the set of catalogues. The 18.6 yr terms of nutation entered the process as parameters the values of which were set by choice from the reservoir of recent determinations. Twenty-five sources contributed to the analysis with the effect that the annual and semi-annual catalogues contain 16 sources on the average. The analysis confirmed the slightly overestimated IAU value of the luni-solar precession yielding a mean correction of -3.59 +/- 1.14 mas/yr from the annual and -3.84 +/- 1.16 mas/yr from the semi-annual catalogues. We discuss these secular terms and their dependence on the adopted nutation terms.

Description: SOFIA, (Stratospheric Observatory for Infrared Astronomy) is a planned 2.5 meter telescope to be installed in a Boeing 747 aircraft and operated at altitudes from 41,000 to 46,000 feet. It will permit routine measurement of infrared radiation inaccessible from the ground-based sites, and observation of astronomical objects and transient events from anywhere in the world. The concept is based on 18 years of experience with NASA's Kuiper Airborne Observatory (KAO), which SOFIA would replace.

Description: This paper describes the status of NASA's Space Infrared Telescope Facility (SIRTF) program. SIRTF will be a cryogenically cooled observatory for infrared astronomy from space and is planned for launch early in the next decade. We discuss a newly modified baseline SIRTF mission and review the state of the detector technology which will provide the heart of SIRTF's scientific advances.

Description: Phase-resolved ROSAT observations of the soft X-ray flux from V444 Cygni confirm the orbital dependence of the flux suggested by analysis of imaging proportional counter (IPC) observations. The X-ray behavior suggests that a region of X-ray emitting gas exists between the 2 stars, probably produced by a collision between the WR and O star winds.

Description: The Broad Band X-ray Telescope (BBXRT) was designed to perform sensitive, moderate resolution spectroscopy of cosmic X-ray sources in the 0.3-10 keV band from the Space Shuttle. During its nine-day flight in December, 1990, the BBXRT observed a variety of supernova remnants and related objects. We present results from some of these observations, emphasizing the ability of the BBXRT to perform spatially-resolved spectroscopy. The improved specral resolution and efficiency over previous instruments makes possible measurements of previously undetectable lines, and the broad bandpass allows simultaneous measurements of lines from oxygen through iron.

Description: New 4 to 8 micron infrared spectroscopic observations of two oxygen-rich stars are presented and combined with IRAS LRS data to span the 4 to 24 micron wavelength range. In the 4 to 8 micron wavelength range, we observe a 7.15 micron (1400/cm) emission feature. The new feature at 7.15 microns is not uniquely correlated with any of the sharply defined 10, 11, 13.1, and 19.7 micron emission features that are known to be present in this class of circumstellar shells, but it does not appear to be correlated with the spectrally broad dust emission in the 10 to 20 micron spectral region. The feature has not been reported previously in any other astronomical environment. A reinterpretation of prior 4 to 8 micron spectroscopy of alpha Ori and R Cas reveals the presence of the 7.15 micron emission in alpha Ori and possibly in R Cas. The spectrally narrow 19.7 micron emission, that is distinctly diffeent than the relatively broad silicate 18 micron emmision feature in oxygen-rich dust shells, is also observed to be present in the LRS spectrum of SAO 197549. The implication of these observations is that a universal astronomical silicate does not exist in oxygen-rich circumstellar shells. This variety is analogous to that observed in interplanetary dust particles and may indicate an intimate relation between the classes of condensates.

Description: The visible optical power emitted from the expansion plumes from 0.4 and 2 km diameter fragments of Shoemaker-Levy are expected to be, approximately 25% and comparable to, the visible solar flux reflected from Jupiter, respectively, for several minutes, and could be easily observed by sensors on the Galileo spacecraft. Earth-based observers can detect these plumes as these expand over the SW limb of Jupiter and come into earth view some minutes after impact!

Description: The Mid-Infrared Spectrometer (MIRS) is one of the four focal plane instruments on the Infrared Telescope in Space (IRTS) mission. The instrument has been constructed, tested, and calibrated in the laboratory and is presently scheduled to be launched by a Japanese expendable launch vehicle as part of the Space Flyer Unit-1 mission in early 1995. The wavelength coverage of the MIRS ranges from 4.5 to 11.7 microns, with a spectral resolution of 0.23 to 0.36 microns. With the cryogenically cooled optics of the IRTS telescope assembly, the MIRS will be able to make an extremely sensitive survey of both point-source and extended objects over an estimated 10% of the sky.

Description: The Infrared Telescope in Space (IRTS) is a cryogenically cooled small infrared telescope that will fly aboard the small space platform Space Flyer Unit. It will survey approximately 10% of the sky with a relatively wide beam during its 20 day emission. Four focal-plane instruments will make simultaneous observations of the sky at wavelengths ranging from 1 to 1000 microns. The IRTS will provide significant information on cosmology, interstellar matter, late-type stars, and interplanetary dust. This paper describes the instrumentation and mission.

Description: We report the results of a survey of low-frequency (LF) plasma waves detected during the Ulysses Jupiter flyby. In the Jovian foreshock, two predominant wave periods are detected: 10(exp 2)-s and 5-s, as measured in the spacecraft frame. The 10(exp 2)-s waves are highly nonlinear propagate at large angles to vector-B(sub 0) (typically 50 deg), are steepened, and sometimes have attached whistler packets. For the interval analyzed the 10(exp 2)-s waves had mixed right-and left-hand polarizations. We argue that these are all consistent with being right-hand magnetosonic waves in the solar wind frame. The 10(exp 2)-s waves with attached whistler are similar to cometary waves. The trailing portions are linearly polaraized and the whistler portions circularly polarized with amplitudes decreasing linearly with time. The emissions are generated by approximately 2-keV protons flowing from the Jovian bow shock/magnetosheath into the upstream region. The instability is the ion beam instability. Higher Z ions were considered as a source of the waves but have been ruled out because of the low sunward velocities needed for their resonance. The 5-s waves have delta vector-B/B(sub 0 approximately = 0.5, are compressive and are left-hand polarized in the spacecraft frame. Local generation by three different resonant interactions were considered and have been ruled out. One possibility is that these waves are whistler mode by-products of the steepened lower-frequency magnetosonic waves. Mirror mode structures were detected throughout the outbound magnetosheath passes. For these structures, the theta(sub kB) values were consistently in the range of 80 deg to 90 deg, exceptionally high values.

Description: The discovery of two objects beyond the orbit of Pluto has extended the heliocentric range of the planetary system and provided tantalizing hints that a large swarm of comets may exist in similarly distant orbits. This discovery has important implications for understanding both the origin of the solar system and the origin of the short-period comets. Subjects covered include the following: the big break -- finding 1992 QB; Kuiper's hypothesis -- an idea takes hold; the dynamicists debate; clues in a disk; Neptune the perturber; discovery of 1993 FW; and unlocking the secrets in primordial ice and dust.

Description: We present spectroscopic observations of some 115 stars of the cluster Blanco 1, extending from the Ca II(H,K) region to the Ca II(I-R) triplet, supporting an age similar to that of the young cluster alpha Persei. The H-alpha absorption equivalent with vs (B - V) diagram forms a well-defined locus, with decreasing absorption equivalent width for decreasing effective temperature, akin to solar neighborhood dwarfs. A large spread in the Ca II surface flux, as a function of (B - V), also indicates the presence of a high degree of surface inhomogeneity, synonymous with high magnetic activity in young stars. A drop-off in the Ca II flux at (B - V) = 1.0 is also similar to the solar neighborhood stars, and shows that the primary chromospheric cooling changes from the Ca II and Mg II lines to the Balmer lines. The mean chromospheric temperature for stars at 4800 K lies between 8000 K and 10,000 K, based on theoretical models, which is somewhat higher than the older solar neighborhood dwarfs. The high mean Ca II surface flux of the sample is also consistent with that of other young clusters. We were able to measure the equivalent width of the Li(6708) line, whose strength as a function of (B - V) indicates an age similar to the young cluster alpha Persei. The lithium abundance decreases with decreasing effective temperature, consistent with the premise of lithium depletion in stars with larger convection zones. Using published photometry and a recent Zero Age Main Sequence (ZAMS) fitting method, we also re-define the distance to the cluster to be 246 pc.

Description: We present three sets of thermal-infrared observations of (2060) Chiron, obtained in 1991, 1993, and 1994. These observations allow the first estimates of the color temperature of Chiron as well as refined estimates of the radius and albedo of its nucleus. 10/20 micrometer color temperatures of 126(sub -6 sup +11) and 137(sub -9 sup +14) K are obtained from the 1993 and 1994 observations, respectively. These temperatures are consistent with the Standard Thermal Model (STM; Lebofsky & Spencer, 1989), but significantly higher than those predicted by the Isothermal Latitude Model. Our estimates of Chiron's radius based on the STM are in agreement with each other, with the observations of Lebofsky et al. (1984), and with recent occultation results (Buie et al., (1993). We obtained values for the radius of 74 +/- 11 km in 1991, 88 +/- 10 and 104 +/- 10 km in 1993, and, 94 +/- 6 and 91 +/- 13 km in 1994.

Description: The Crab and the Vela pulsars were observed to emit pulsed radiation in the radio, optical, X-ray and gamma ray windows of the electromagnetic spectrum. In all cases the emission appears in two peaks, except in the case of radio emission from the Vela pulsar which is single peaked. With accurate peak position determinations made possible by recent observations, it is interesting to see if the intra-peak separations at various wavelengths across the electromagnetic spectrum bear any relation to each other. We report here an emerging trend of a monotonic decrease of the intra-peak separation with energy. The rate of decrease is faster in the case of the Vela pulsar than in that of the Crab pulsar. Even the case of single peaked radio emission by the Vela pulsar can be viewed as being consistent with this trend. These trends provide both an opportunity and a challenge to realistic modelings of pulsed emissions by these objects.

Description: Formulas are developed for the transformation of ecliptical orbital elements from B 1950 to J 2000. The results are compared with those recommended by IAU Commission 20. Some drawbacks to the Commission 20 formulation are pointed out and we develop procedures which are consistent with standard precessional formulations.

Description: In this paper we study the influence of line-merging regions at the intermediate long-wavelength side of a continuum threshold on the computed model atmosphere structure and predicted spectrum. In order to model these regions sufficiently accurately, we have developed two concepts. First, we have extended the occupation probability formalism of Hummer and Mihalas to non-local thermodynamic equilibrium (non-LTE) plasmas. Second, in order to treat the very complicated opacity in the line merging region, we have generalized the concept of opacity distribution functions to treat non-LTE situations. All Rydberg states are consistently included within this framework, so that no arbitrary cutoff of high (LTE) levels is made. We have calculated several pure hydrogen models atmospheres for two effective temperatures, T(sub eff) = 20000 and 35000 K, and discussed the differences between models calculated with various treatments of the line merging. In particular, we have shown that the error in the predicted profiles of Balmer lines resulting from the neglect of line merging is typically of the order of 3-4%, while the errors in the far-UV portion of the Balmer continuum reaches 15-35%. The errors generally decrease with increasing effective temperature. At the same time, the internal accuracy of the models is shown to be about or below 0.5% for all predicted spectral features. We conclude that for interpreting current high-accuracy spectrophotometric observations models including the line merging are necessary, and that the formalism developed in this paper is capable of providing a sufficiently accurate and robust modeling technique.

Description: We develop a new technique for extracting three-dimensional information from multiday solar Very Large Array (VLA) observations. While standard stereoscopic methods provide a three-dimensional view of an object by combining simultaneous observations from two different aspect angles, we relax the condition of simultaneity and exploit solar rotation to vary the aspect angle. The solar radio images are decomposed into Gaussian source components, which are then cross-correlated in maps from preceding and following days. This provides measurements of the three-dimensional position of correlated source centroids. In this first paper, we describe the stereoscopic method and perform tests with simulated and real radio maps (from the VLA at 20 cm), in order to study the accuracy of altitude measurements, and the limitations introduced by (i) source confusion, (ii) source motion, and (iii) the assumed differential rotation rate. The tests demonstrate that (i) the information content of a VLA map relevant for stereoscopic correlation can be conveniently represented in terms of a small number of Gaussian components; (ii) the fitting of the three-dimensional source position is stable within a numerical accuracy of less than or approximately equal to 0.02 map pixels, (iii) the relative accuracy of the altitude determination is uniform over the solar disk, and (iv) source confusion does not affect the accuracy of stereoscopic position measurements for sources with a signal-to-noise ratio of greater than or approximately equal to 36.

Description: We have obtained 3.2-3.6 micron spectra, with a resolution lambda/delta-lambda approximately 750, of the protostar Mon R2/IRS-3 and of Elias 16, a background K giant behind the Taurus molecular cloud. A feature at 3,482 microns (2872/cm), with a full width at half-maximum of 0.09 microns (76/cm), is clearly seen in Mon R2/IRS-3. This feature is not detected in Elias 16. The 3.482 micron feature in Mon R2/IRS-3 is similar to a feature at 3.466-3.478 microns (2875-2885/cm) detected by Allamandola et al. in four protostars and attributed by these authors to a CH stretch in hydrocarbons dominated by sp3-bonded carbon. Neither Mon R2/IRS-3 nor Elias 16 shows absorption at 3.540 microns (2825/cm), which has been detected in two of the four protostars observed by Allamandola et al. and attributed by them to CH3OH ice. Our limit on CH3OH ice toward Elias 16 is compared to models of gas-grain chemistry in dark clouds. Our results confirm those of Allamandola et al. that at this resolution the 3.4 micron absorption due to dust in molecular clouds has very different spectral structure than that due to dust in the diffuse interstellar medium.

Description: ASCA observations of the gravitational lens and Butcher-Oemler cluster Abell 370 (z = 0.37) give kT = 8.8 +/- 0.8 keV and A = 0.5 +/- 0.1 cosmic. If the gas were isothermal the implied cluster mass would be M(sub vir) = (1.5 +/- 0.4) x 10(exp 15) solar masses, a value consistent with the optically-determined virial mass. We detect iron K line emission with high confidence. This measurement increases, by a large factor, the lookback time at which the presence of iron in the intracluster medium has been established. The iron abundance is marginally higher than that of low-redshift clusters of similar temperature, so our results are consistent with models in which all enrichment occurs before the epoch corresponding to z = 0.37.

Description: The transformation validity question utilizing resulting data from direct numerical simulations (DNS) of supersonic, isothermal cold wall channel flow was investigated. The DNS results stood for a wide scope of parameter and were suitable for the purpose of examining the generality of Van Driest transformation. The Van Driest law of the wall can be obtained from the inner-layer similarity arguments. It was demonstrated that the Van Driest transformation cannot be incorporated to collapse the sublayer and log-layer velocity profiles simultaneously. Velocity and temperature predictions according to the preceding composite mixing-length model were presented. Despite satisfactory congruity with the DNS data, the model must be perceived as an engineering guide and not as a rigorous analysis.

Description: A Monte Carlo solution technique has been formulated to predict the radiative heat transfer in three-dimensional, inhomogeneous participating media which exhibit spectrally dependent emission and absorption and anisotropic scattering. Details of the technique and selected numerical sensitivities are discussed. The technique was applied to a problem involving a medium composed of a gas mixture of carbon dioxide and nitrogen and suspended carbon particles. A homogeneous medium was modeled to examine the effect of total pressure and carbon-particle concentration on radiative heat transfer. Variation in total pressure, over the range studied, had minimal effect on the amount of heat radiated to the enclosure walls and on the radiative-flux distribution within the medium. Increases in the carbon particle concentration produced significantly higher heat fluxes at the boundaries and altered the radiative flux distribution. The technique was then applied to an inhomogeneous medium to examine effects of specific temperature and carbon particle concentration distributions on radiative heat transfer. For the inhomogeneous conditions examined, the largest radiative flux divergence occurs near the center of the medium and the regions near some enclosure walls act as energy sinks.

Description: We have conducted timing observations of the eclipsing millisecond binary pulsar PSR B1957+20, extending the span of data on this pulsar to more than five years. During this time the orbital period of the system has varied by roughly Delta P(sub b)/P(sub b) = 1.6 x 10(exp -7), changing quardratically with time and displaying with time and displaying an orbital period second derivative of P(sub b) = (1.43 +/- 0.08) x 10(exp -18)/sec. The previous measurement of a large negative orbital period derivative reflected only the short-term behavior of the system during the early observations; the orbital period derivative is now positive. If, as we suspect, the PSR B1957+20 system is undergoing quasi-cyclic orbital period variations similar to those found in other close binaries such as Algol and RS CVn, then the 0.025 solar mass companion to PSR B1957+20 is most likely non-degenerate, convective, and magnetically active.

Description: This Letter describes the results of an initial study of Galactic extinction and the colors of Galactic stellar populations in the near-IR using the Diffuse Infrared Background Experiment (DIRBE) aboard the Cosmic Background Explorer (COBE) spacecraft. The near-IR reddening observed by DIRBE is consistent with the extinction law tabulated by Rieke & Lebofsky (1985). The distribution of dust and stars in most of the first and fourth quadrants of the Galactic plane (0 deg less than l less than 90 deg, and 270 deg less than l less than 360 deg, respectively) can be modeled as a stellar background source seen through up to approximately 4 mag of extinction at 1.25 micrometers. The unreddened near-IR colors of the Galactic disk are similar to those of late-K and M giants. The Galactic bulge exhibits slightly bluer colors in the 2.2-3.5 micrometers range, as noted by Terndrup et al. (1991). Star-forming regions exhibit colors that indicate the presence of a approximately 900 K continuum produced by hot dust or polycyclic aromatic hydrocarbons (PAHs) contributing at wavelengths as short as 3.5 micrometers.

Description: Numerical results obtained with direct simulation Monte Carlo and Navier-Stokes methods are presented for a Mach-20 nitrogen flow about a 70-deg blunted cone. The flow conditions simulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how the wake structure changes as a function of rare faction, what the afterbody levels of heating are, and to what limits the continuum models are realistic as rarefunction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the afterbody sting are emphasized in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

Description: An Advanced Satellite for Cosmology and Astrophysics (ASCA) observation of Cygnus X-3 on 1993 June 11, in its X-ray high intensity state, has for the first time resolved the broad iron K line emission into three components: a He-like line at 6.67 +/- 0.01 keV, a H-like line at 6.96 +/- 0.02 keV, and a neutral line at 6.37 +/- 0.03 keV. The line intensities of the 6.67 keV and 6.96 keV lines are modulated with the 4.8 hr orbital period and are maximum when the continuum intensity is minimum. There is a sharp minimum of the line intensity on the rising phase of the continuum intensity. An iron absorption edge is observed at 7.19 +/- 0.02 keV. The optical depth of the absorption edge varies from 0.3 to 0.5 and is in anti-phase with the overall X-ray continuum modulation. The observed complexity of the iron K line region is greater than that had been assumed in previous spectral modeling based on observations with lower resolution detectors.

Description: We have observed the intermediate polar EX Hya with Advanced Satellite for Cosmology and Astrophysics (ASCA), and have clearly detected He-like and H-like K alpha emission lines from Mg to Fe. The intensity ratios of these pairs of lines are not compatible with an isothermal plasma, and a temperature distribution can no longer be parameterized by the conventional two emission component model. We have successfully decomposed iron line emission into thermal plasma component and flourescent component. The equivalent width of the flourescent component is approximately 80 eV.

Description: With the high sensitivity and spectral resolution of the Advanced Satellite for Cosmology and Astrophysics (ASCA) satellite, we have discovered strong emission lines from the H-like and/or He-like ions of Ne, Mg, Si, and S as well as Fe L and confirmed the complex structure of Fe K line emission in the Seyfert II galaxy NGC 1068. The continuum emission above 3 keV exhibits rather flat shape with no evidence of low energy absorption. The overall X-ray spectrum can be well explained with a model involving starburst activity plus an obscured active galactic nucleus.

Description: The symbiotic star EG Andromedae has recently been the subject of several studies investigating its wind properties. Late-type giants are usually considered to have winds driven by radiation pressure on dust. Indeed, the derived wind velocity for EG Andromedae is consistent with this model. We point out here that there is no appreciable dust opacity in the wind of EG Andromedae using constraints on extinction limits from International Ultraviolet Explorer (IUE) and far infrared fluxes from Infrared Astronomy Satellite (IRAS). An alternate mechanism must operate in this star. We suggest that the wind can be driven by radiation pressure on molecular lines.

Description: This atlas contains sample spectra from the far-ultraviolet observations of 32 symbiotic stars obtained with the International Ultraviolet Explorer (IUE) satellite. In all, 394 low-resolution spectra from the short-wavelength primary (SWP) camera covering the range 1200-2000 A have been extracted from the IUE archive, calibrated, and measured. Absolute line fluxes and wavelengths for the prominent emission lines have been tabulated. Tables of both the general properties of these symbiotics and of features specific to the spectrum of each are included. The spectra shown are representative of the different classes of symbiotic stars that are currently in the IUE archive. These include known eclipsing systems and those that have been observed in outburst (as well as quiescence).

Description: The displacement of a radio-emitting star around the barycenter of a possible planetary system can be measured by astrometric very long baseline interferometry (VLBI) observations. We have observed the radio-emitting star sigma(sup 2) CrB at 8 epochs over 5 years by VLBI and fitted its 5 astrometric parameters to the observed coordinates. The post-fit coordinate residuals have and rms scatter of 0.22 milliarcseconds and show no systematic behavior. We use this result to set a limit on the presence of planets around sigma(sup 2) CrB and conclude that our present VLBI astrometric precision corresponds to the threshold to detect a Jupiter-like planet around this star. We also discuss the astrometric monitoring program of 11 radio-emitting stars that we are conducting for the Hipparcos space mission and its possible contribution to a long-term planet search program.

Description: We investigate the question of disk formation during the protostar phase. We build on the results of Keene and Masson (1990) whose analysis of L1551 showed the millimeter continuum emission comes from both an unresolved circumstellar component, i.e., a disk and an extended cloud core. We model the dust continuum emission from the cloud core and show how it is important at 1.3 mm but negligible at 2.7 mm. Combining new 2.7 mm Owens Valley Interferometer data of IRAS-Dense cores with data from the literature we conclude that massive disks are also seen toward a number of other sources. However, 1.3 mm data from the IRAM 30 m telescope for a larger sample shows that massive disks are relatively rare, occurring around perhaps 5% of young embedded stars. This implies that either massive disks occur briefly during the embedded phase or that relatively few young stars form massive disks. At 1.3 mm the median flux of IRAS-Dense cores is nearly the same as T Tauri stars in the sample of Beckwith et al. (1990). We conclude that the typical disk mass during the embedded phase is nearly the same or less than the typical disk mass during the T Tauri phase.

Description: The accomplishments of the Space Shuttle Endeavor's mission to repair the Hubble Space Telescope (HST) are discussed. Images are presented that were taken with the HST Faint Object Camera before and after the installation of the Corrective Optics Space Telescope and Axial Replacement. Images are also presented that were taken with the old Wide Field and Planetary Camera (WFPC) and its replacement, WFPC 2. Comparisons between before-repair and after-repair images clearly indicate that the mission was a success.

Description: We have analyzed two data sets of Jupiter's limb darkening taken before (March 1990) and after (March 1991) a South Equatorial Belt (SEB) disturbance. The vertical aerosol structure has been determined for the SEB, the Equatorial Zone (EZ), and the North Equatorial Belt (NEB). We find that all regions can be well represented with a two-cloud layer model, with slightly higher NH3-cloud top altitude in the EZ than in the two belts. The cloud single-scattering albedos tilda-omega(sub cl) (determined within plus or minus 0.003) show little difference in the EZ (from 0.9956 to 0.9952) and the NEB (from 0.9925 to 0.9926) during a 1-year period, while a relatively larger decrease is found in the SEB from 0.9975 to 0.9957, indicating the recovery of the belt through the disturbance. The NH3-cloud top altitude in the SEB shows the most prominent change, lower in 1990 and higher in 1991. We suggest that the brightening and the darkening of the SEB may originate not only from the changes of tilda-omega(sub cl) but also from changes of the ammonia condensation level in the atmosphere, probably due to a decrease of the tropospheric temperature and/or an increase of the ammonia vapor.

Description: The Simplified Shuttle Payload Thermal Analyzer program (SSPTA) was developed to aid in the evaluation of thermal design concepts of instruments to be flown in the Space Shuttle cargo bay. SSPTA consists of a collection of programs that are currently used in the thermal analysis of spacecraft and have been modified for quick, preliminary analysis of payloads. SSPTA includes a reduced math model of the Shuttle cargo bay to simplify use of the program for payload analysis. One of the prime objectives in developing SSPTA was to create a program which was easy to use. With SSPTA, the user required input is simple and the user is free from many of the concerns of computer usage such as disk space handling, tape usage, and complicated program control. Although SSPTA was designed primarily to analyze Shuttle payloads, it can easily be used to perform thermal analysis in other situations. SSPTA is comprised of a system of data files called 'bins', a master program, and a set of thermal subprograms. The bin system is a collection of disk files which contain data required by or computed by the thermal subprograms. SSPTA currently has the capability of handling 50 bins. The master program serves primarily as a manager for the bin system and its interaction with the thermal subprograms. Input to the master program consists of simple user commands which direct the data manipulation procedures, prepare the data for these procedures, and call the appropriate thermal subprograms. The subprograms of SSPTA are all based on programs which have been used extensively in the analysis of orbiting spacecraft and space hardware. Subprogram CONSHAD uses the user supplied geometric radiation model to compute black body view factors, shadow factors, and a description of the surface model. The subprogram WORKSHEET uses the surface model description, optical property data, and node assignment data to prepare input for SCRIPTF. Subprogram SCRIPTF computes the inverses of the infrared (IR) and ultraviolet (UV) radiation transfer equations; it also computes the radiation coupling between nodes in the thermal model. Subprogram ORBITAL uses the shadow tables to compute incident flux intensities on each surface in the geometric model. Subprogram ABSORB uses these flux intensities combined with the IR and UV inverses to compute the IR and UV fluxes absorbed by each surface. The radiation couplings from SCRIPTF and the absorbed fluxes from ABSORB are used by subprogram TTA to compute the temperature and power balance for each node in the thermal model. Output consists of tabulated data from each of the subprograms executed during a particular analysis. Due to the modular form of SSPTA, analyses may be run in whole or in part, and new subprograms may be added by the user. SSPTA is written in FORTRAN for use on a DEC VAX-11/780. SSPTA was originally developed in 1977 for use on IBM 370 series computers. This version is an update which was ported to the VAX in 1980.

Description: Remote Interactive Particle-tracing (RIP) is a distributed-graphics program which computes particle traces for computational fluid dynamics (CFD) solution data sets. A particle trace is a line which shows the path a massless particle in a fluid will take; it is a visual image of where the fluid is going. The program is able to compute and display particle traces at a speed of about one trace per second because it runs on two machines concurrently. The data used by the program is contained in two files. The solution file contains data on density, momentum and energy quantities of a flow field at discrete points in three-dimensional space, while the grid file contains the physical coordinates of each of the discrete points. RIP requires two computers. A local graphics workstation interfaces with the user for program control and graphics manipulation, and a remote machine interfaces with the solution data set and performs time-intensive computations. The program utilizes two machines in a distributed mode for two reasons. First, the data to be used by the program is usually generated on the supercomputer. RIP avoids having to convert and transfer the data, eliminating any memory limitations of the local machine. Second, as computing the particle traces can be computationally expensive, RIP utilizes the power of the supercomputer for this task. Although the remote site code was developed on a CRAY, it is possible to port this to any supercomputer class machine with a UNIX-like operating system. Integration of a velocity field from a starting physical location produces the particle trace. The remote machine computes the particle traces using the particle-tracing subroutines from PLOT3D/AMES, a CFD post-processing graphics program available from COSMIC (ARC-12779). These routines use a second-order predictor-corrector method to integrate the velocity field. Then the remote program sends graphics tokens to the local machine via a remote-graphics library. The local machine interprets the graphics tokens and draws the particle traces. The program is menu driven. RIP is implemented on the silicon graphics IRIS 3000 (local workstation) with an IRIX operating system and on the CRAY2 (remote station) with a UNICOS 1.0 or 2.0 operating system. The IRIS 4D can be used in place of the IRIS 3000. The program is written in C (67%) and FORTRAN 77 (43%) and has an IRIS memory requirement of 4 MB. The remote and local stations must use the same user ID. PLOT3D/AMES unformatted data sets are required for the remote machine. The program was developed in 1988.

Description: This software package includes two programs, the KPD12 and the KPD12P. Both programs utilizes the vortex-blob method to simulate flow around solid bodies, in an unbounded domain using the KPD12, with periodicity in one direction using the KPD12P. The main advantage of the vortex-blob method is the ability to handle situations involving arbitrary shapes including multiple bodies. The user just supplies points on the solid boundaries; there is no grid. The KPD12 program has worked successfully on bluff bodies, stalled wings, and multiple-element airfoils. The KPD12P program has been used successfully on high-solidity separated cascades and on cases of rotating stall in cascades of thin airfoils. However, they do not capture subtle viscous effects such as incipient separation and friction drag. The KPD12 and the KPD12P programs apply the vortex-blob method to time-dependent, high-Reynolds-number flows around solid bodies. Both programs solve the two-dimensional incompressible Navier-Stokes equations, neglecting the viscous effects away from the walls. By creating new vortices along the wall at every time step, they treat the no-penetration and no-slip boundary conditions while using an influence matrix. The code automatically controls the number of vortices. Furthermore, the code has the option of treating the boundary layers by simple integral methods to determine the separation points. The KPD12 outputs forces, moments, and pressure distributions on the bodies. The KPD12P also outputs the turning angle and loss of total pressure. The source code is in Cray FORTRAN and contains a few calls to Cray vector functions which are vectorized with the Cray compiler. However, substitutes for these vector functions are provided. The code is set up to plot the bodies, vortex positions, and streamlines using the DISSPLA graphics software. The software requires a mainframe computer with at least 589k of memory available running under COS 1.16. KPD12 was developed in 1988.

Description: The Multimission Interactive Picture Planner, MIP, is a scientifically accurate and fast, 3D animation program for deep space. MIP is also versatile, reasonably comprehensive, portable, and will run on microcomputers. New techniques were developed to rapidly perform the calculations and transformations necessary to animate scientifically accurate 3D space. At the same time, portability is maintained, as the transformations and clipping have been written in FORTRAN 77 code. MIP was primarily designed to handle Voyager, Galileo, and the Space Telescope. It can, however, be adapted to handle other missions. The space simulation consists of a rotating body (usually a planet), any natural satellites, a spacecraft, the sun, stars, descriptive labelling, and field of view boxes. The central body and natural satellites are tri-axial wireframe representations with terminators, limbs, and landmarks. Hidden lines are removed for the central body and natural satellites, but not for the scene as a whole so that bodies may be seen behind one another. The program has considerable flexibility in its step time, observer position, viewed object, field of view, etc. Most parameters may be changed from the keyboard while the simulation is running. When MIP is executed it will ask the user for a control file, which should be prepared before execution. The control file identifies which mission MIP should simulate, the star catalog files, the ephemerides files to be used, the central body, planets, asteroids, and comets, and solar system landmarks and constants such as planets, asteroids, and comets. The control file also describes the fields of view. Control files are included to simulate the Voyager 1 encounter at Jupiter and the Giotto spacecraft's flyby of Halley's comet. Data is included for Voyager 1 and 2 (all 6 planetary encounters) and Giotto. MIP was written for an IBM PC or compatibles. It requires 512K of RAM, a CGA or compatible graphics adapter, and DOS 2.0 or higher. Users must supply their own graphics primitives to clear the screen, change the color, and connect 2D points with straight lines. Also, the users must tie in the graphics primitives along with their ephemeris readers. (MIP does everything else including clipping.) MIP was developed in 1988.

Description: The Steady State Thermal Analysis Program (STEADY) provides the thermal designer with a quick and convenient method for calculating heat loads and temperatures. STEADY can be used on small nodal networks for conceptual or preliminary thermal design and analysis. STEADY will accept up to 20 nodes of fixed or variable temperature, with constant or temperature-dependent thermal conductivities, and any set of consistent units. In a steady state thermal network, the heat balance on each variable temperature node must sum to zero. The general heat transfer equations are solved with a Newton-Raphson technique and refined by a fourth order quartic solution. Input data includes the number of nodes, number of boundary nodes, the fixed temperatures at all boundary nodes, initial temperature guesses for variable nodes, impressed heat loads, conduction and radiation coefficients, and control parameters such as convergence criteria, maximum iterations, and damping factors. The output is stored in a print file and tabulates final temperatures and heat flows for all nodes. STEADY is menu driven and allows the user to save files for future modification. STEADY is written in FORTRAN 77 (Ryan McFarland's RMFORTRAN) for interactive execution and has been implemented on the IBM PC computer series under DOS with a central memory requirement of approximately 92K of 8 bit bytes using a math coprocessor, and 103K bytes without the coprocessor. This program was developed in 1987.

Description: The Thermal Radiation Analyzer System, TRASYS, is a computer software system with generalized capability to solve the radiation related aspects of thermal analysis problems. TRASYS computes the total thermal radiation environment for a spacecraft in orbit. The software calculates internode radiation interchange data as well as incident and absorbed heat rate data originating from environmental radiant heat sources. TRASYS provides data of both types in a format directly usable by such thermal analyzer programs as SINDA/FLUINT (available from COSMIC, program number MSC-21528). One primary feature of TRASYS is that it allows users to write their own driver programs to organize and direct the preprocessor and processor library routines in solving specific thermal radiation problems. The preprocessor first reads and converts the user's geometry input data into the form used by the processor library routines. Then, the preprocessor accepts the user's driving logic, written in the TRASYS modified FORTRAN language. In many cases, the user has a choice of routines to solve a given problem. Users may also provide their own routines where desirable. In particular, the user may write output routines to provide for an interface between TRASYS and any thermal analyzer program using the R-C network concept. Input to the TRASYS program consists of Options and Edit data, Model data, and Logic Flow and Operations data. Options and Edit data provide for basic program control and user edit capability. The Model data describe the problem in terms of geometry and other properties. This information includes surface geometry data, documentation data, nodal data, block coordinate system data, form factor data, and flux data. Logic Flow and Operations data house the user's driver logic, including the sequence of subroutine calls and the subroutine library. Output from TRASYS consists of two basic types of data: internode radiation interchange data, and incident and absorbed heat rate data. The flexible structure of TRASYS allows considerable freedom in the definition and choice of solution method for a thermal radiation problem. The program's flexible structure has also allowed TRASYS to retain the same basic input structure as the authors update it in order to keep up with changing requirements. Among its other important features are the following: 1) up to 3200 node problem size capability with shadowing by intervening opaque or semi-transparent surfaces; 2) choice of diffuse, specular, or diffuse/specular radiant interchange solutions; 3) a restart capability that minimizes recomputing; 4) macroinstructions that automatically provide the executive logic for orbit generation that optimizes the use of previously completed computations; 5) a time variable geometry package that provides automatic pointing of the various parts of an articulated spacecraft and an automatic look-back feature that eliminates redundant form factor calculations; 6) capability to specify submodel names to identify sets of surfaces or components as an entity; and 7) subroutines to perform functions which save and recall the internodal and/or space form factors in subsequent steps for nodes with fixed geometry during a variable geometry run. There are two machine versions of TRASYS v27: a DEC VAX version and a Cray UNICOS version. Both versions require installation of the NASADIG library (MSC-21801 for DEC VAX or COS-10049 for CRAY), which is available from COSMIC either separately or bundled with TRASYS. The NASADIG (NASA Device Independent Graphics Library) plot package provides a pictorial representation of input geometry, orbital/orientation parameters, and heating rate output as a function of time. NASADIG supports Tektronix terminals. The CRAY version of TRASYS v27 is written in FORTRAN 77 for batch or interactive execution and has been implemented on CRAY X-MP and CRAY Y-MP series computers running UNICOS. The standard distribution medium for MSC-21959 (CRAY version without NASADIG) is a 1600 BPI 9-track magnetic tape in UNIX tar format. The standard distribution medium for COS-10040 (CRAY version with NASADIG) is a set of two 6250 BPI 9-track magnetic tapes in UNIX tar format. Alternate distribution media and formats are available upon request. The DEC VAX version of TRASYS v27 is written in FORTRAN 77 for batch execution (only the plotting driver program is interactive) and has been implemented on a DEC VAX 8650 computer under VMS. Since the source codes for MSC-21030 and COS-10026 are in VAX/VMS text library files and DEC Command Language files, COSMIC will only provide these programs in the following formats: MSC-21030, TRASYS (DEC VAX version without NASADIG) is available on a 1600 BPI 9-track magnetic tape in VAX BACKUP format (standard distribution medium) or in VAX BACKUP format on a TK50 tape cartridge; COS-10026, TRASYS (DEC VAX version with NASADIG), is available in VAX BACKUP format on a set of three 6250 BPI 9-track magnetic tapes (standard distribution medium) or a set of three TK50 tape cartridges in VAX BACKUP format. TRASYS was last updated in 1993.

Description: SUNTRACKER computes the azimuth and elevation angles of the sun over a user defined period of time for a given position on the moon. The program gets the selenographic (moon centered) position of the sun at a given Julian date, then converts the selenographic position of the sun into azimuth and elevation at a given position on the moon. The selenographic coordinate system is based on the equatorial plane of the moon. The origin of this system is referenced to the mean center of the apparent lunar disk. This is the point of the surface of the moon intersected by the lunar radius directed towards the center of the earth when the moon is at its mean ascending node. Selenographic longitudes are measured positive in the direction towards Mare Crisium from the lunar meridian passing through the origin. Selenographic latitudes are measured positive towards the northern hemisphere containing Mare Serenitatis, from the lunar equator. The selenographic co-longitude is obtained by subtracting the selenographic longitude from either 90 or 450 degrees. SUNTRACKER performs two main operations. The first handles the Julian and calendar date calculations. The second operation calculates the right ascension and declination of the sun and moon on the Earth's equatorial coordinate system from the adjusted Julian date. These coordinates are then transferred into the ecliptic coordinate system. From here the moon's position is configured to the heliocentric ecliptic coordinate system. It is from this system that the sun's selenographic position is determined. Algorithms are then used to compute both the physical and optical librations of the moon's motion. The limitations, restrictions and assumptions for SUNTRACKER are as follows: 1) The orbital elements used in this program do not account for nutation, aberration and precession. 2) The selenographic coordinates computed by this program are based on the 1961 Astronomical Ephemeris algorithms. In 1981, a new analytic theory on the librations of the moon was adopted by The Astronomical Almanac. This improved the method of calculating selenographic coordinates. The selenographic coordinates computed by SUNTRACKER are identical to the pre-1981 Astronomical Almanac values. A comparison between the program's values and the 1993 Astronomical Almanac's values found the maximum deviations in longitude and latitude to be .030 deg. and .034 deg. respectively. The average deviations were .013 deg. in selenographic co-longitude and .017 deg. in selenographic latitude. SUNTRACKER is written in FORTRAN 77 for IBM PC compatible computers running MS-DOS. The sample executable included on the media requires 64K of RAM and the Lahey FORTRAN 77 run time library for execution. The standard distribution medium for this program is a 5.25 inch 360K MS-DOS format diskette.

Description: TDIGG is a fast and versatile program for generating two-dimensional computational grids for use with finite-difference flow-solvers. Both algebraic and elliptic grid generation systems are included. The method for grid generation by algebraic transformation is based on an interpolation algorithm and the elliptic grid generation is established by solving the partial differential equation (PDE). Non-uniform grid distributions are carried out using a hyperbolic tangent stretching function. For algebraic grid systems, interpolations in one direction (univariate) and two directions (bivariate) are considered. These interpolations are associated with linear or cubic Lagrangian/Hermite/Bezier polynomial functions. The algebraic grids can subsequently be smoothed using an elliptic solver. For elliptic grid systems, the PDE can be in the form of Laplace (zero forcing function) or Poisson. The forcing functions in the Poisson equation come from the boundary or the entire domain of the initial algebraic grids. A graphics interface procedure using the Silicon Graphics (GL) Library is included to allow users to visualize the grid variations at each iteration. This will allow users to interactively modify the grid to match their applications. TDIGG is written in FORTRAN 77 for Silicon Graphics IRIS series computers running IRIX. This package requires either MIT's X Window System, Version 11 Revision 4 or SGI (Motif) Window System. A sample executable is provided on the distribution medium. It requires 148K of RAM for execution. The standard distribution medium is a .25 inch streaming magnetic IRIX tape cartridge in UNIX tar format. This program was developed in 1992.

Description: The Systems Improved Numerical Fluids Analysis Code, SINFAC, consists of additional routines added to the April 1983 revision of SINDA, a general thermal analyzer program. The purpose of the additional routines is to allow for the modeling of active heat transfer loops. The modeler can simulate the steady-state and pseudo-transient operations of 16 different heat transfer loop components including radiators, evaporators, condensers, mechanical pumps, reservoirs and many types of valves and fittings. In addition, the program contains a property analysis routine that can be used to compute the thermodynamic properties of 20 different refrigerants. SINFAC can simulate the response to transient boundary conditions. SINFAC was first developed as a method for computing the steady-state performance of two phase systems. It was then modified using CNFRWD, SINDA's explicit time-integration scheme, to accommodate transient thermal models. However, SINFAC cannot simulate pressure drops due to time-dependent fluid acceleration, transient boil-out, or transient fill-up, except in the accumulator. SINFAC also requires the user to be familiar with SINDA. The solution procedure used by SINFAC is similar to that which an engineer would use to solve a system manually. The solution to a system requires the determination of all of the outlet conditions of each component such as the flow rate, pressure, and enthalpy. To obtain these values, the user first estimates the inlet conditions to the first component of the system, then computes the outlet conditions from the data supplied by the manufacturer of the first component. The user then estimates the temperature at the outlet of the third component and computes the corresponding flow resistance of the second component. With the flow resistance of the second component, the user computes the conditions down stream, namely the inlet conditions of the third. The computations follow for the rest of the system, back to the first component. On the first pass, the user finds that the calculated outlet conditions of the last component do not match the estimated inlet conditions of the first. The user then modifies the estimated inlet conditions of the first component in an attempt to match the calculated values. The user estimated values are called State Variables. The differences between the user estimated values and calculated values are called the Error Variables. The procedure systematically changes the State Variables until all of the Error Variables are less than the user-specified iteration limits. The solution procedure is referred to as SCX. It consists of two phases, the Systems phase and the Controller phase. The X is to imply experimental. SCX computes each next set of State Variables in two phases. In the first phase, SCX fixes the controller positions and modifies the other State Variables by the Newton-Raphson method. This first phase is the Systems phase. Once the Newton-Raphson method has solved the problem for the fixed controller positions, SCX next calculates new controller positions based on Newton's method while treating each sensor-controller pair independently but allowing all to change in one iteration. This phase is the Controller phase. SINFAC is available by license for a period of ten (10) years to approved licensees. The licenced program product includes the source code for the additional routines to SINDA, the SINDA object code, command procedures, sample data and supporting documentation. Additional documentation may be purchased at the price below. SINFAC was created for use on a DEC VAX under VMS. Source code is written in FORTRAN 77, requires 180k of memory, and should be fully transportable. The program was developed in 1988.

Description: QUIKVIS computes the times during an Earth orbit when geometric requirements are satisfied for observing celestial objects. The observed objects may be fixed (stars, etc.) or moving (sun, moon, planets). QUIKVIS is useful for preflight analysis by those needing information on the availability of celestial objects to be observed. Two types of analyses are performed by QUIKVIS. One is used when specific objects are known, the other when targets are unknown and potentially useful regions of the sky must be identified. The results are useful in selecting candidate targets, examining the effects of observation requirements, and doing gross assessments of the effects of the orbit's right ascension of the ascending node (RAAN). The results are not appropriate when high accuracy is needed (e.g. for scheduling actual mission operations). The observation duration is calculated as a function of date, orbit node, and geometric requirements. The orbit right ascension of the ascending node can be varied to account for the effects of an uncertain launch time of day. The orbit semimajor axis and inclination are constant throughout the run. A circular orbit is assumed, but a simple program modification will allow eccentric orbits. The geometric requirements that can be processed are: 1) minimum separation angle between the line of sight to the object and the earth's horizon; 2) minimum separation angle between the line of sight to the object and the spacecraft velocity vector; 3) maximum separation angle between the line of sight to the object and the zenith direction; and 4) presence of the spacecraft in the earth's shadow. The user must supply a date or date range, the spacecraft orbit and inclination, up to 700 observation targets, and any geometric requirements to be met. The primary output is the time per orbit that conditions are satisfied, with options for sky survey maps, time since a user-specified orbit event, and bar graphs illustrating overlapping requirements. The output is printed in visually convenient lineprinter form but is also available on data files for use by postprocessors such as external XY plotters. QUIKVIS is written in FORTRAN 77 for batch or interactive execution and has been implemented on a DEC VAX 11/780 operating under VMS with a central memory requirement of approximately 500K of 8 bit bytes. QUIKVIS was developed in 1986 and revised in 1987.

Description: The Thermal Radiation Analyzer System, TRASYS, is a computer software system with generalized capability to solve the radiation related aspects of thermal analysis problems. TRASYS computes the total thermal radiation environment for a spacecraft in orbit. The software calculates internode radiation interchange data as well as incident and absorbed heat rate data originating from environmental radiant heat sources. TRASYS provides data of both types in a format directly usable by such thermal analyzer programs as SINDA/FLUINT (available from COSMIC, program number MSC-21528). One primary feature of TRASYS is that it allows users to write their own driver programs to organize and direct the preprocessor and processor library routines in solving specific thermal radiation problems. The preprocessor first reads and converts the user's geometry input data into the form used by the processor library routines. Then, the preprocessor accepts the user's driving logic, written in the TRASYS modified FORTRAN language. In many cases, the user has a choice of routines to solve a given problem. Users may also provide their own routines where desirable. In particular, the user may write output routines to provide for an interface between TRASYS and any thermal analyzer program using the R-C network concept. Input to the TRASYS program consists of Options and Edit data, Model data, and Logic Flow and Operations data. Options and Edit data provide for basic program control and user edit capability. The Model data describe the problem in terms of geometry and other properties. This information includes surface geometry data, documentation data, nodal data, block coordinate system data, form factor data, and flux data. Logic Flow and Operations data house the user's driver logic, including the sequence of subroutine calls and the subroutine library. Output from TRASYS consists of two basic types of data: internode radiation interchange data, and incident and absorbed heat rate data. The flexible structure of TRASYS allows considerable freedom in the definition and choice of solution method for a thermal radiation problem. The program's flexible structure has also allowed TRASYS to retain the same basic input structure as the authors update it in order to keep up with changing requirements. Among its other important features are the following: 1) up to 3200 node problem size capability with shadowing by intervening opaque or semi-transparent surfaces; 2) choice of diffuse, specular, or diffuse/specular radiant interchange solutions; 3) a restart capability that minimizes recomputing; 4) macroinstructions that automatically provide the executive logic for orbit generation that optimizes the use of previously completed computations; 5) a time variable geometry package that provides automatic pointing of the various parts of an articulated spacecraft and an automatic look-back feature that eliminates redundant form factor calculations; 6) capability to specify submodel names to identify sets of surfaces or components as an entity; and 7) subroutines to perform functions which save and recall the internodal and/or space form factors in subsequent steps for nodes with fixed geometry during a variable geometry run. There are two machine versions of TRASYS v27: a DEC VAX version and a Cray UNICOS version. Both versions require installation of the NASADIG library (MSC-21801 for DEC VAX or COS-10049 for CRAY), which is available from COSMIC either separately or bundled with TRASYS. The NASADIG (NASA Device Independent Graphics Library) plot package provides a pictorial representation of input geometry, orbital/orientation parameters, and heating rate output as a function of time. NASADIG supports Tektronix terminals. The CRAY version of TRASYS v27 is written in FORTRAN 77 for batch or interactive execution and has been implemented on CRAY X-MP and CRAY Y-MP series computers running UNICOS. The standard distribution medium for MSC-21959 (CRAY version without NASADIG) is a 1600 BPI 9-track magnetic tape in UNIX tar format. The standard distribution medium for COS-10040 (CRAY version with NASADIG) is a set of two 6250 BPI 9-track magnetic tapes in UNIX tar format. Alternate distribution media and formats are available upon request. The DEC VAX version of TRASYS v27 is written in FORTRAN 77 for batch execution (only the plotting driver program is interactive) and has been implemented on a DEC VAX 8650 computer under VMS. Since the source codes for MSC-21030 and COS-10026 are in VAX/VMS text library files and DEC Command Language files, COSMIC will only provide these programs in the following formats: MSC-21030, TRASYS (DEC VAX version without NASADIG) is available on a 1600 BPI 9-track magnetic tape in VAX BACKUP format (standard distribution medium) or in VAX BACKUP format on a TK50 tape cartridge; COS-10026, TRASYS (DEC VAX version with NASADIG), is available in VAX BACKUP format on a set of three 6250 BPI 9-track magnetic tapes (standard distribution medium) or a set of three TK50 tape cartridges in VAX BACKUP format. TRASYS was last updated in 1993.

Description: INS3D computes steady-state solutions to the incompressible Navier-Stokes equations. The INS3D approach utilizes pseudo-compressibility combined with an approximate factorization scheme. This computational fluid dynamics (CFD) code has been verified on problems such as flow through a channel, flow over a backwardfacing step and flow over a circular cylinder. Three dimensional cases include flow over an ogive cylinder, flow through a rectangular duct, wind tunnel inlet flow, cylinder-wall juncture flow and flow through multiple posts mounted between two plates. INS3D uses a pseudo-compressibility approach in which a time derivative of pressure is added to the continuity equation, which together with the momentum equations form a set of four equations with pressure and velocity as the dependent variables. The equations' coordinates are transformed for general three dimensional applications. The equations are advanced in time by the implicit, non-iterative, approximately-factored, finite-difference scheme of Beam and Warming. The numerical stability of the scheme depends on the use of higher-order smoothing terms to damp out higher-frequency oscillations caused by second-order central differencing. The artificial compressibility introduces pressure (sound) waves of finite speed (whereas the speed of sound would be infinite in an incompressible fluid). As the solution converges, these pressure waves die out, causing the derivation of pressure with respect to time to approach zero. Thus, continuity is satisfied for the incompressible fluid in the steady state. Computational efficiency is achieved using a diagonal algorithm. A block tri-diagonal option is also available. When a steady-state solution is reached, the modified continuity equation will satisfy the divergence-free velocity field condition. INS3D is capable of handling several different types of boundaries encountered in numerical simulations, including solid-surface, inflow and outflow, and far-field boundaries. Three machine versions of INS3D are available. INS3D for the CRAY is written in CRAY FORTRAN for execution on a CRAY X-MP under COS, INS3D for the IBM is written in FORTRAN 77 for execution on an IBM 3090 under the VM or MVS operating system, and INS3D for DEC RISC-based systems is written in RISC FORTRAN for execution on a DEC workstation running RISC ULTRIX 3.1 or later. The CRAY version has a central memory requirement of 730279 words. The central memory requirement for the IBM is 150Mb. The memory requirement for the DEC RISC ULTRIX version is 3Mb of main memory. INS3D was developed in 1987. The port to the IBM was done in 1990. The port to the DECstation 3100 was done in 1991. CRAY is a registered trademark of Cray Research Inc. IBM is a registered trademark of International Business Machines. DEC, DECstation, and ULTRIX are trademarks of the Digital Equipment Corporation.

Description: SAGE, Self Adaptive Grid codE, is a flexible tool for adapting and restructuring both 2D and 3D grids. Solution-adaptive grid methods are useful tools for efficient and accurate flow predictions. In supersonic and hypersonic flows, strong gradient regions such as shocks, contact discontinuities, shear layers, etc., require careful distribution of grid points to minimize grid error and produce accurate flow-field predictions. SAGE helps the user obtain more accurate solutions by intelligently redistributing (i.e. adapting) the original grid points based on an initial or interim flow-field solution. The user then computes a new solution using the adapted grid as input to the flow solver. The adaptive-grid methodology poses the problem in an algebraic, unidirectional manner for multi-dimensional adaptations. The procedure is analogous to applying tension and torsion spring forces proportional to the local flow gradient at every grid point and finding the equilibrium position of the resulting system of grid points. The multi-dimensional problem of grid adaption is split into a series of one-dimensional problems along the computational coordinate lines. The reduced one dimensional problem then requires a tridiagonal solver to find the location of grid points along a coordinate line. Multi-directional adaption is achieved by the sequential application of the method in each coordinate direction. The tension forces direct the redistribution of points to the strong gradient region. To maintain smoothness and a measure of orthogonality of grid lines, torsional forces are introduced that relate information between the family of lines adjacent to one another. The smoothness and orthogonality constraints are direction-dependent, since they relate only the coordinate lines that are being adapted to the neighboring lines that have already been adapted. Therefore the solutions are non-unique and depend on the order and direction of adaption. Non-uniqueness of the adapted grid is acceptable since it makes possible an overall and local error reduction through grid redistribution. SAGE includes the ability to modify the adaption techniques in boundary regions, which substantially improves the flexibility of the adaptive scheme. The vectorial approach used in the analysis also provides flexibility. The user has complete choice of adaption direction and order of sequential adaptions without concern for the computational data structure. Multiple passes are available with no restraint on stepping directions; for each adaptive pass the user can choose a completely new set of adaptive parameters. This facility, combined with the capability of edge boundary control, enables the code to individually adapt multi-dimensional multiple grids. Zonal grids can be adapted while maintaining continuity along the common boundaries. For patched grids, the multiple-pass capability enables complete adaption. SAGE is written in FORTRAN 77 and is intended to be machine independent; however, it requires a FORTRAN compiler which supports NAMELIST input. It has been successfully implemented on Sun series computers, SGI IRIS's, DEC MicroVAX computers, HP series computers, the Cray YMP, and IBM PC compatibles. Source code is provided, but no sample input and output files are provided. The code reads three datafiles: one that contains the initial grid coordinates (x,y,z), one that contains corresponding flow-field variables, and one that contains the user control parameters. It is assumed that the first two datasets are formatted as defined in the plotting software package PLOT3D. Several machine versions of PLOT3D are available from COSMIC. The amount of main memory is dependent on the size of the matrix. The standard distribution medium for SAGE is a 5.25 inch 360K MS-DOS format diskette. It is also available on a .25 inch streaming magnetic tape cartridge in UNIX tar format or on a 9-track 1600 BPI ASCII CARD IMAGE format magnetic tape. SAGE was developed in 1989, first released as a 2D version in 1991 and updated to 3D in 1993.

Description: INS3D computes steady-state solutions to the incompressible Navier-Stokes equations. The INS3D approach utilizes pseudo-compressibility combined with an approximate factorization scheme. This computational fluid dynamics (CFD) code has been verified on problems such as flow through a channel, flow over a backwardfacing step and flow over a circular cylinder. Three dimensional cases include flow over an ogive cylinder, flow through a rectangular duct, wind tunnel inlet flow, cylinder-wall juncture flow and flow through multiple posts mounted between two plates. INS3D uses a pseudo-compressibility approach in which a time derivative of pressure is added to the continuity equation, which together with the momentum equations form a set of four equations with pressure and velocity as the dependent variables. The equations' coordinates are transformed for general three dimensional applications. The equations are advanced in time by the implicit, non-iterative, approximately-factored, finite-difference scheme of Beam and Warming. The numerical stability of the scheme depends on the use of higher-order smoothing terms to damp out higher-frequency oscillations caused by second-order central differencing. The artificial compressibility introduces pressure (sound) waves of finite speed (whereas the speed of sound would be infinite in an incompressible fluid). As the solution converges, these pressure waves die out, causing the derivation of pressure with respect to time to approach zero. Thus, continuity is satisfied for the incompressible fluid in the steady state. Computational efficiency is achieved using a diagonal algorithm. A block tri-diagonal option is also available. When a steady-state solution is reached, the modified continuity equation will satisfy the divergence-free velocity field condition. INS3D is capable of handling several different types of boundaries encountered in numerical simulations, including solid-surface, inflow and outflow, and far-field boundaries. Three machine versions of INS3D are available. INS3D for the CRAY is written in CRAY FORTRAN for execution on a CRAY X-MP under COS, INS3D for the IBM is written in FORTRAN 77 for execution on an IBM 3090 under the VM or MVS operating system, and INS3D for DEC RISC-based systems is written in RISC FORTRAN for execution on a DEC workstation running RISC ULTRIX 3.1 or later. The CRAY version has a central memory requirement of 730279 words. The central memory requirement for the IBM is 150Mb. The memory requirement for the DEC RISC ULTRIX version is 3Mb of main memory. INS3D was developed in 1987. The port to the IBM was done in 1990. The port to the DECstation 3100 was done in 1991. CRAY is a registered trademark of Cray Research Inc. IBM is a registered trademark of International Business Machines. DEC, DECstation, and ULTRIX are trademarks of the Digital Equipment Corporation.

Description: The ability to treat arbitrary boundary shapes is one of the most desirable characteristics of a method for generating grids, including those about airfoils. In a grid used for computing aerodynamic flow over an airfoil, or any other body shape, the surface of the body is usually treated as an inner boundary and often cannot be easily represented as an analytic function. The GRAPE computer program was developed to incorporate a method for generating two-dimensional finite-difference grids about airfoils and other shapes by the use of the Poisson differential equation. GRAPE can be used with any boundary shape, even one specified by tabulated points and including a limited number of sharp corners. The GRAPE program has been developed to be numerically stable and computationally fast. GRAPE can provide the aerodynamic analyst with an efficient and consistent means of grid generation. The GRAPE procedure generates a grid between an inner and an outer boundary by utilizing an iterative procedure to solve the Poisson differential equation subject to geometrical restraints. In this method, the inhomogeneous terms of the equation are automatically chosen such that two important effects are imposed on the grid. The first effect is control of the spacing between mesh points along mesh lines intersecting the boundaries. The second effect is control of the angles with which mesh lines intersect the boundaries. Along with the iterative solution to Poisson's equation, a technique of coarse-fine sequencing is employed to accelerate numerical convergence. GRAPE program control cards and input data are entered via the NAMELIST feature. Each variable has a default value such that user supplied data is kept to a minimum. Basic input data consists of the boundary specification, mesh point spacings on the boundaries, and mesh line angles at the boundaries. Output consists of a dataset containing the grid data and, if requested, a plot of the generated mesh. The GRAPE program is written in FORTRAN IV for batch execution and has been implemented on a CDC 6000 series computer with a central memory requirement of approximately 135K (octal) of 60 bit words. For plotted output the commercially available DISSPLA graphics software package is required. The GRAPE program was developed in 1980.

Description: The Thermal Radiation Analyzer System, TRASYS, is a computer software system with generalized capability to solve the radiation related aspects of thermal analysis problems. TRASYS computes the total thermal radiation environment for a spacecraft in orbit. The software calculates internode radiation interchange data as well as incident and absorbed heat rate data originating from environmental radiant heat sources. TRASYS provides data of both types in a format directly usable by such thermal analyzer programs as SINDA/FLUINT (available from COSMIC, program number MSC-21528). One primary feature of TRASYS is that it allows users to write their own driver programs to organize and direct the preprocessor and processor library routines in solving specific thermal radiation problems. The preprocessor first reads and converts the user's geometry input data into the form used by the processor library routines. Then, the preprocessor accepts the user's driving logic, written in the TRASYS modified FORTRAN language. In many cases, the user has a choice of routines to solve a given problem. Users may also provide their own routines where desirable. In particular, the user may write output routines to provide for an interface between TRASYS and any thermal analyzer program using the R-C network concept. Input to the TRASYS program consists of Options and Edit data, Model data, and Logic Flow and Operations data. Options and Edit data provide for basic program control and user edit capability. The Model data describe the problem in terms of geometry and other properties. This information includes surface geometry data, documentation data, nodal data, block coordinate system data, form factor data, and flux data. Logic Flow and Operations data house the user's driver logic, including the sequence of subroutine calls and the subroutine library. Output from TRASYS consists of two basic types of data: internode radiation interchange data, and incident and absorbed heat rate data. The flexible structure of TRASYS allows considerable freedom in the definition and choice of solution method for a thermal radiation problem. The program's flexible structure has also allowed TRASYS to retain the same basic input structure as the authors update it in order to keep up with changing requirements. Among its other important features are the following: 1) up to 3200 node problem size capability with shadowing by intervening opaque or semi-transparent surfaces; 2) choice of diffuse, specular, or diffuse/specular radiant interchange solutions; 3) a restart capability that minimizes recomputing; 4) macroinstructions that automatically provide the executive logic for orbit generation that optimizes the use of previously completed computations; 5) a time variable geometry package that provides automatic pointing of the various parts of an articulated spacecraft and an automatic look-back feature that eliminates redundant form factor calculations; 6) capability to specify submodel names to identify sets of surfaces or components as an entity; and 7) subroutines to perform functions which save and recall the internodal and/or space form factors in subsequent steps for nodes with fixed geometry during a variable geometry run. There are two machine versions of TRASYS v27: a DEC VAX version and a Cray UNICOS version. Both versions require installation of the NASADIG library (MSC-21801 for DEC VAX or COS-10049 for CRAY), which is available from COSMIC either separately or bundled with TRASYS. The NASADIG (NASA Device Independent Graphics Library) plot package provides a pictorial representation of input geometry, orbital/orientation parameters, and heating rate output as a function of time. NASADIG supports Tektronix terminals. The CRAY version of TRASYS v27 is written in FORTRAN 77 for batch or interactive execution and has been implemented on CRAY X-MP and CRAY Y-MP series computers running UNICOS. The standard distribution medium for MSC-21959 (CRAY version without NASADIG) is a 1600 BPI 9-track magnetic tape in UNIX tar format. The standard distribution medium for COS-10040 (CRAY version with NASADIG) is a set of two 6250 BPI 9-track magnetic tapes in UNIX tar format. Alternate distribution media and formats are available upon request. The DEC VAX version of TRASYS v27 is written in FORTRAN 77 for batch execution (only the plotting driver program is interactive) and has been implemented on a DEC VAX 8650 computer under VMS. Since the source codes for MSC-21030 and COS-10026 are in VAX/VMS text library files and DEC Command Language files, COSMIC will only provide these programs in the following formats: MSC-21030, TRASYS (DEC VAX version without NASADIG) is available on a 1600 BPI 9-track magnetic tape in VAX BACKUP format (standard distribution medium) or in VAX BACKUP format on a TK50 tape cartridge; COS-10026, TRASYS (DEC VAX version with NASADIG), is available in VAX BACKUP format on a set of three 6250 BPI 9-track magnetic tapes (standard distribution medium) or a set of three TK50 tape cartridges in VAX BACKUP format. TRASYS was last updated in 1993.

Description: INS3D computes steady-state solutions to the incompressible Navier-Stokes equations. The INS3D approach utilizes pseudo-compressibility combined with an approximate factorization scheme. This computational fluid dynamics (CFD) code has been verified on problems such as flow through a channel, flow over a backwardfacing step and flow over a circular cylinder. Three dimensional cases include flow over an ogive cylinder, flow through a rectangular duct, wind tunnel inlet flow, cylinder-wall juncture flow and flow through multiple posts mounted between two plates. INS3D uses a pseudo-compressibility approach in which a time derivative of pressure is added to the continuity equation, which together with the momentum equations form a set of four equations with pressure and velocity as the dependent variables. The equations' coordinates are transformed for general three dimensional applications. The equations are advanced in time by the implicit, non-iterative, approximately-factored, finite-difference scheme of Beam and Warming. The numerical stability of the scheme depends on the use of higher-order smoothing terms to damp out higher-frequency oscillations caused by second-order central differencing. The artificial compressibility introduces pressure (sound) waves of finite speed (whereas the speed of sound would be infinite in an incompressible fluid). As the solution converges, these pressure waves die out, causing the derivation of pressure with respect to time to approach zero. Thus, continuity is satisfied for the incompressible fluid in the steady state. Computational efficiency is achieved using a diagonal algorithm. A block tri-diagonal option is also available. When a steady-state solution is reached, the modified continuity equation will satisfy the divergence-free velocity field condition. INS3D is capable of handling several different types of boundaries encountered in numerical simulations, including solid-surface, inflow and outflow, and far-field boundaries. Three machine versions of INS3D are available. INS3D for the CRAY is written in CRAY FORTRAN for execution on a CRAY X-MP under COS, INS3D for the IBM is written in FORTRAN 77 for execution on an IBM 3090 under the VM or MVS operating system, and INS3D for DEC RISC-based systems is written in RISC FORTRAN for execution on a DEC workstation running RISC ULTRIX 3.1 or later. The CRAY version has a central memory requirement of 730279 words. The central memory requirement for the IBM is 150Mb. The memory requirement for the DEC RISC ULTRIX version is 3Mb of main memory. INS3D was developed in 1987. The port to the IBM was done in 1990. The port to the DECstation 3100 was done in 1991. CRAY is a registered trademark of Cray Research Inc. IBM is a registered trademark of International Business Machines. DEC, DECstation, and ULTRIX are trademarks of the Digital Equipment Corporation.

Description: The ability to treat arbitrary boundary shapes is one of the most desirable characteristics of a method for generating grids. 3DGRAPE is designed to make computational grids in or about almost any shape. These grids are generated by the solution of Poisson's differential equations in three dimensions. The program automatically finds its own values for inhomogeneous terms which give near-orthogonality and controlled grid cell height at boundaries. Grids generated by 3DGRAPE have been applied to both viscous and inviscid aerodynamic problems, and to problems in other fluid-dynamic areas. 3DGRAPE uses zones to solve the problem of warping one cube into the physical domain in real-world computational fluid dynamics problems. In a zonal approach, a physical domain is divided into regions, each of which maps into its own computational cube. It is believed that even the most complicated physical region can be divided into zones, and since it is possible to warp a cube into each zone, a grid generator which is oriented to zones and allows communication across zonal boundaries (where appropriate) solves the problem of topological complexity. 3DGRAPE expects to read in already-distributed x,y,z coordinates on the bodies of interest, coordinates which will remain fixed during the entire grid-generation process. The 3DGRAPE code makes no attempt to fit given body shapes and redistribute points thereon. Body-fitting is a formidable problem in itself. The user must either be working with some simple analytical body shape, upon which a simple analytical distribution can be easily effected, or must have available some sophisticated stand-alone body-fitting software. 3DGRAPE does not require the user to supply the block-to-block boundaries nor the shapes of the distribution of points. 3DGRAPE will typically supply those block-to-block boundaries simply as surfaces in the elliptic grid. Thus at block-to-block boundaries the following conditions are obtained: (1) grids lines will match up as they approach the block-to-block boundary from either side, (2) grid lines will cross the boundary with no slope discontinuity, (3) the spacing of points along the line piercing the boundary will be continuous, (4) the shape of the boundary will be consistent with the surrounding grid, and (5) the distribution of points on the boundary will be reasonable in view of the surrounding grid. 3DGRAPE offers a powerful building-block approach to complex 3-D grid generation, but is a low-level tool. Users may build each face of each block as they wish, from a wide variety of resources. 3DGRAPE uses point-successive-over-relaxation (point-SOR) to solve the Poisson equations. This method is slow, although it does vectorize nicely. Any number of sophisticated graphics programs may be used on the stored output file of 3DGRAPE though it lacks interactive graphics. Versatility was a prominent consideration in developing the code. The block structure allows a great latitude in the problems it can treat. As the acronym implies, this program should be able to handle just about any physical region into which a computational cube or cubes can be warped. 3DGRAPE was written in FORTRAN 77 and should be machine independent. It was originally developed on a Cray under COS and tested on a MicroVAX 3200 under VMS 5.1.

Description: The recent discovery of an object bound to the sun by gravity but situated beyond the orbit of Pluto is discussed. The discovery of this object, designated 1992 QB1, may confirm the existence of the Kuiper belt, envisioned to be a disk of far-flung comets left over from the origin of the solar system. The implications of the new discovery for the future search for distant comets beyond Pluto are examined.

Description: In the present paper, direct numerical methods by which to simulate the spatially developing free shear flows in the transitional region are described and the numerical results of a spatially developing plane wake are presented. The incompressible time-dependent Navier-Stokes equations were solved using Pade finite difference approximations in the streamwise direction, a mapped pseudospectral Fourier method in the cross-stream direction, and a third-order compact Runge-Kutta scheme for time advancement. The unstable modes of the Orr-Sommerfeld equations were used to perturb the inlet of the wake. Statistical analyses were performed and some numerical results were compared with experimental measurements. When only the fundamental mode is forced, the energy spectra show amplification of the fundamental and its higher harmonics. In this case, unperturbed alternate vortices develop in the saturation region of the wake. The phase jitter around the fundamental frequency plays a critical role in generating vortices of random shape and spacing. Large- and small-scale distortions of the fundamental structure are observed. Pairing of vortices of the same sign is observed, as well as vortex coupling of vortices of the opposite sign.

Description: This book first reviews the overall aspects and background information related to thermal radiation heat transfer and incorporates new general information, advances in analytical and computational techniques, and new reference material. Coverage focuses on radiation from opaque surfaces, radiation interchange between various types of surfaces enclosing a vacuum or transparent medium, and radiation including the effects of partially transmitting media, such as combustion gases, soot, or windows. Boundary conditions and multiple layers are discussed with information on radiation in materials with nonunity refractive indices.

Description: The present approach to the prediction of instability generation that is due to the interaction of freestream disturbances with regions of subscale variations in surface boundary conditions can account for the finite Reynolds number effects, while furnishing a framework for the study of receptivity in compressible flow and in 3D boundary layers. The approach is illustrated for the case of Tollmien-Schlichting wave generation in a Blasius boundary layer, due to the interaction of a freestream acoustic wave with a localized wall inhomogeneity. Results are presented for the generation of viscous and inviscid instabilities in adverse pressure-gradient boundary layers, supersonic boundary layer instabilities, and cross-flow vortex instabilities.

Description: Flaws of the depth analysis technique of distinguishing photographic plate defects are discussed. The method is shown to be unreliable in confirming or disproving the astrophysical nature of optical transients. The arguments of Greiner and of Zytkow to the contrary are criticized.