ALBERT

Description: NASA/GRO grant NAG 5-2081, at the University of Chicago, has provided support for a broad program of theoretical research in nuclear astrophysics and related areas, with regard to gamma-ray and hard X-ray emission from classical nova explosions. This research emphasized the possible detection of 22Na gamma-ray line emission from nearby novae involving ONeMg white dwarfs, the detailed examination of 26Al production in novae, and the possible detection of the predicted early gamma ray emission from novae that arises from the decay of the short lived, positron emitting isotopes of CNO elements. Studies of nova related problems have consumed an increasing fraction of the Principal Investigator's research efforts over the past decade. Current research addresses problems associated with the standard model for the outbursts of the classical novae: the occurrence of thermonuclear runaways (TNR) in the accreted hydrogen rich envelopes on white dwarfs in close binary systems (see, e.g., the reviews by Truran 1982; and Shara 1989). Research in progress and planned for the next three years has three main objectives: (1) to gain an improved understanding of the early evolution of the light curves of, particularly, the fastest novae; (2) to gain an improved understanding of the relative importance of the various possible mechanisms of envelope hydrogen depletion (e.g. winds, common envelope driven mass loss, and nuclear burning) to the long term evolution of novae in outburst; and (3) to seek to provide a somewhat more definitive statement of the role of classical novae in nucleosynthesis. Our proposed 2-D studies of convection during the early phases of the TNR and our systematic attempt to incorporate an improved treatment of radiation hydrodynamics into the hydrodynamic code utilized in our calculations, are particularly relevant to the first of these objectives. Further 2-D studies of the effects of common envelope evolution are intended to provide more realistic constraints on the mass depletion mechanisms. Finally, detailed calculations of the thermonuclear history of the matter ejected in novae will be carried out for representative nova configurations involving both carbon-oxygen (CO) and oxygen-neon-magnesium (ONeMg) white dwarfs.

Description: The funds from this grant were used to support observations and analysis with the International Ultraviolet Explorer (IUE) satellite telescope. The main area of scientific research concerned the variability analyses of ultraviolet spectra of Active Galactic Nuclei, primarily quasars, Seyfert galaxies, and BL Lacertae objects. The Colorado group included, at various times, the P.I. (J.M. Shull), Research Associate Dr. Rick Edelson, and graduate students Jon Saken, Elise Sachs, and Steve Penton. A portion of the work was also performed by CU undergraduate student Cheong-ming Fu. A major product of the effort was a database of all IUE spectra of active galactic nuclei. This database is being analyzed to obtain spectral indices, line fluxes, and continuum fluxes for over 500 AGN. As a by-product of this project, we implemented a new, improved technique of spectral extraction of IUE spectra, which has been used in several AGN-WATCH campaigns (on the Seyfert galaxy NGC 4151 and on the BL Lac object PKS 2155-304).

Description: The NPARC Alliance is a partnership between the NASA Lewis Research Center (LeRC) and the USAF Arnold Engineering Development Center (AEDC) dedicated to the establishment of a national CFD capability, centered on the NPARC Navier-Stokes computer program. The three main tasks of the Alliance are user support, code development, and validation. The present paper is a status report on the validation effort. It describes the validation approach being taken by the Alliance. Representative results are presented for laminar and turbulent flat plate boundary layers, a supersonic axisymmetric jet, and a glancing shock/turbulent boundary layer interaction. Cases scheduled to be run in the future are also listed. The archive of validation cases is described, including information on how to access it via the Internet.

Description: The progress made during 1995 on the Monte-Carlo gamma-ray spectrum simulation program BSIMUL is discussed. Several features have been added, including the ability to model shields that are tapered cylinders. Several simulations were made on the Near Earth Asteroid Rendezvous detector.

Description: The second guide star catalog (GSC-2) is a project that aims to create a complete catalog of stars and galaxies to about the 18th magnitude, and would contain colors, magnitudes, positions and proper motions. The catalog would provide an object list for the construction of an input catalog for possible future astrometric satellites, such as the global astrometric interferometer for astrophysics (GAIA) satellite. With a schedule that is compatible with the projected timeframe of GAIA, the GSC-2 could be available in time to support the preparatory astrophysical observations for an input catalog. The availability of the digitalized images for the preparation of finding charts would improve the efficiency of this task.

Description: We present radio observations of the gravitational lens PKS 1830-211 at 8.4 and 15 GHz acquired using the Very Large Array. The observations were made over a 13 month period. Significant flux density changes over this period provide strong constraints on the time delay between the two lensed images and suffest a value of 44 +/- 9 days. This offers new direct evidence that this source is indeed a gravitational lens. The lens distance is dependent upon the model chosen, but reasonable limits on the mass of the lensing galaxy suggest that it is unlikely to be at a redshift less than a few tenths, and may well be significantly more distant.

Description: The power laws are approximately f(exp -1.9), f(exp -1.9), and f(exp -2.1) respectively for the Grigg-Shjellerup (GS), Giacobini-Zinner (GZ), and Halley (H) comets. Other than similarities in the power spectra, the magnetic field turbulence is considerably different at the three comets. Phase steepening is demonstrated to occur at the trailing edges of the GS waves. This is probably due to nonlinear steepening plus dispersion of the left-hand mode components, i.e., the turbulence is whistler-mode. This too can be explained by nonlinear steepening plus dispersion of the magnetosonic waves. At the level of GS and GZ turbulence development when the spacecraft measurements were made, classical three-wave processes, such as the decay or modulation instabilities do not appear to play important roles. It is most likely that the nonlinear steepening and dispersive time scales are more rapid than three-wave processes, and the latter had not had time to develop for the relatively new turbulence. The wave turbulence at Halley is linearly polarized. The exact nature of this turbulence is still not well understood. Several possibilities are suggested, based on a preliminary analyses.

Description: A study of X-ray emission from five short-period Algol-type binaries based on observations with Advanced Satellite for Cosmology and Astrophysics (ASCA) and ROSAT is presented. We have observed RZ Cas with both satellites, and beta Per, U Cep, delta Lib, and TW Dra with ROSAT. Significant intensity variations are seen in the X-ray emission from RZ Cas, U Cep, TW Dra, and delta Lib. These variations seem unrelated to the eclipsing behavior of these systems and are probably due to either rotational modulation of compact active regions on the surfaces of the chromospherically active secondary components or to flaring activity in the systems. The spectra of all but one of the systems require the presence of at least two discrete plasma components with different temperatures (0.6 - 0.7 keV, and approximately 2 keV) and the abundances of the medium-Z elements 20% - 50% of the solar photospheric values. The high resolving power and signal-to-noise ratio of the ASCA spectra allow us to individually constrain the coronal abundances of O, Ne, Mg, Si, S, and Fe in RZ Cas. We demonstrate that, if we use the elemental abundances and temperatures obtained from the analysis of their ASCA spectra as (fixed) inputs, to fit the ROSAT PSPC spectra well requires the presence of a third component (kT approximately 0.2 - 0.3 keV) in RZ Cas and beta Per. A continuous emission measure model of the power-law type (EM(T) variesas (T/T(sub max)(sup alpha)) generally gives a poor fit to the ASCA and ROSAT data on most sources. Circumstellar or circumbinary absorbing matter seems to be present in some of these systems, as indicated by the variable total column density needed to fit their X-ray spectra.

Description: We study the onset of a pure Marangoni convection in a liquid layer with two deformable interfaces in the no-gravity environment. Both oscillatory and stationary instabilities are considered for a wide range of parameters. It is shown that only stationary instability is possible when surface tension at the colder interface is lower than that at the hotter one. Oscillatory instability tends to disappear and to be replaced by the stationary instability with increase of the Prandtl number and decrease of surface tension at the colder interface.

Description: The 183-GHz water vapor line was tentatively detected on Mars in January 1991, with the IRAM 30-m millimeter antenna, under extremely dry atmospheric conditions. The measurement refers to the whole disk. The spectral line, although marginally detected, can be fit with a constant H2O mixing ratio of 1.0 x 10(exp -5), which corresponds to a water abundance of 1 pr-microns; in any case, an upper limit of 3 pr-microns is inferred. This value is comparable to the very small abundances measured by Clancy (1992) 5 weeks before our observation and seems to imply both seasonal and long-term variations in the martian water cycle.

Description: Titan's large free eccentricity results in significant tidal dissipation. This can be used to constrain the existence and depth of hydrocarbon oceans. A hydrodynamical-numerical 2 deg ocean tide model has been constructed to investigate this connection. The model allows some simple land configurations. The results indicate that existence of such an ocean over the age of the solar system is hard to explain, as is the existence of the large eccentricity itself. If such an ocean exists, it is likely to be more than 500 m deep, ignoring the influence of land masses.

Description: A consistent solution of the radiative transfer equation characterizing photon transport in a semi-infinite medium of refractive index greater than or equal to one is obtained following the method of Sobolev. Fresnel specular reflection, Snell's law and isotropic scattering are assumed. An algorithm is developed and its accuracy is demonstrated. A numerical Laplace transform inversion leads to an efficient evaluation for the interior flux and source function distributions.

Description: 24 years of Lunar Laser Ranging (LLR) observations and 16 years of Very Long Baseline Interferometry (VLBI) observations are combined in a global analysis to yield improved estimates of the Earth's precession and nutation. The correction to the International Astronomical Union (IAU) (1976) precession constant inferred from this joint VLBI/LLR analysis is -3.00 +/- 0.20 milliarcsec/yr (mas/yr). A significant obliquity rate correction of -0.20 +/- 0.08 mas/yr is also found. In all, 32 forced nutation coefficients are estimated. These coefficients confirm that the IAU (1980) nutation theory is in error by several mas. The estimated nutation coeficients are found to vary by as much as several tenths of mas, depending on the a priori nutation model used to analyze the VLBI and LLR data. Forced circular nutations derived from this analysis agree with the ZMOA-1990-2 nutation theory at the 0.2 mas level for the 18.6 yr terms, and at the 0.05 mas level for the other terms (periods less than or = 1 yr). A retrograde free core nutation with an amplitude of 0.20 mas is also detected. Its phase is found to be very sensitive to the precise value of the free core nutation period used in the solution. Separate analyses of four independent subsets of the LVBI data indicate no significant variations of the free core nutation since 1988. The pre-1988 estimates of the free core nutation are consistent with the post-1988 estimates but are not accurate enough to rule out possible variations of the free core nutation at these earlier epochs.

Description: The Ulysses spacecraft has gathered data from within flows from the Sun's southern polar coronal hole, the first in situ measurement of this region. We present a brief analysis of the heliospheric magnetic field data from this region, using a fractal method. As is the case near the ecliptic, estimated spectral exponents are near 5/3 on spacecraft scales of seconds to minutes. On longer time scales, however, there appears to be a significantly different population in polar flows, which is similar to that found by the Helios spacecraft in fast solar wind flows at 0.3 AU.

Description: A calculation, employing a detailed model of neutral oxygen, is carried out to give fluorescent line intensities expected in a long-proposed photoexcitation by accidental resonance (PAR) process in which hydrogen Lyman-beta photoexcites the oxygen spectrum. The results pertain to the optically thin case but provide an upper limit to the fluorescent intensities which can be attained. They are applied to analyze line ratios involving the strong 8446 A line observed in classical novae during the diffusion-enhanced and Orion phases. Operation of the PAR process in the novae is verified. It is found that photoexcitation rates in the ejecta reach values greater than 0.1/sec, corresponding to hydrogen Lyman-beta radiation field intensities greater than 1250 ergs/cm/sec/sr.

Description: We have investigated the properties of the OMC-2 and OMC-3 cores in the Orion giant molecular cloud using high spatial spectral resolution observations of several transitions of the (13)CO, C(18)O, C(S-32) and C(S-34) molecules taken with the SEST telescope. The OMC-2 core consists of one clump (22 solar mass) with a radius of 0.11 pc surrounded by a cluster of 11 discrete infrared sources. The H2 column density and volume density in the center of this clump are 2 x 10(exp 22)/sq cm and 9 x 10(exp 5)/cu cm respectively. From a comparison between physical parameters derived from C(18)O and C(S-32) observations we conclude that the molecular envelope around the core has been completely removed by these sources and that only the very dense gas is left. OMC-3 shows a more complex elongated structure in C(18)O and CS than OMC-2. The C(S-32) and C(S-34) maps show that the denser region can be separated into at least sub-cores of roughly equal sizes (radius approximately equals 0.13 pc), with n(H2) = 6 x 10(exp 5)/cu cm, and a mass of 10 solar mass (from C(S-32)). The very different masses obtained for the central core from C(18)O and C(S-32) (55 and 12 solar mass respectively) indicate that a massive envelope is still present around the very dense sub-cores. We report the first detection of several molecular outflows in OMC-3. The presence of an IRAS source and the first detection of these outflows confirm that star formation is going on in OMC-3. Based on the different physical properties of these regions compared with OMC-1, OMC-2 appears to be in an intermediate evolutionary stage between OMC-1 and OMC-3.

Description: The depletion of condensable elements onto grains in gaseous nebulae can provide evidence that dust is well mixed with the ionized gas. Al and CA are two of the most depleted elements in the general interstellar medium, and it is therefore important to measure their abundances within the ionized region of nebulae. We compute a large grid of photoionization models and identify sets of line ratios which are relatively insensitive to stellar and nebular parameters, and are thus excellent diagnostics for determining relative abundances. Based on the absence of the ((Ca II) lambda lambda 291, 7324 doublet and the detection of Al II) lambda lambda 2660, 2669 in the ultraviolet, we determine the extent of aluminum and calcium depletion onto grains in NGC 7027 and the Orion Nebula. Our results show a approximately 0.3 dex depletion for Al, but a depletion of more than two and a half orders of magnitude for Ca. A similar calculation based on Mg II lambda 2798 yields roughly a 0.8 dex depletion for Mg. This reaffirms the discrepancy between depletion determined from high and low ionization Mg lines. We also find evidence for a 'depletion gradient' in Ca in NGC 7027, since the calcium depletion we infer for the outer, more neutral regions using (Ca II) is somewhat higher than that inferred for the inner high-ionization region, using (Ca v). This gradient can test current models of the survival of grains within hot ionized gas.

Description: Results are reported of the Surface Tension Driven Convection Experiment (STDCE) aboard USML-1 Spacelab. Steady and transient thermocapillary flows were investigated in a 10 cm dia. circular container filled with 10 Cs silicone oil. The velocity and temperature fields were studied in detail under various conditions. It is shown in this paper how the Marangoni number affects the velocity field. A numerical analysis of the flows was also conducted and its results were compared to the experimental data. Good agreement is shown.

Description: The modern scenario of evolution of massive binary systems predicts the existence of a subclass of binary radio pulsars (PSRs) with black holes (BHs). Their Galactic number was evaluated as approximately 1 per 1000 single pulsars (Lipunov et al. 1994b). Distinctive properties of such binaries would be (1) mass of the unseen companion M(sub c) greater than 3-4 solar mass and (2) absence of eclipses of the pulsar radiation with no distinctive variance of the dispersion measure along the pulsar orbit. The pulsars themselves must be similar to standard isolated ones. The recently discovered binary 1 s pulsar PSR B0042-73 = PSR J0045-7319 in the Small Magellanic Cloud (SMC) with a massive companion in a highly elongated (eccentricity e =0.8) 51 day orbit (Kaspi et al. 1994) may be the first such pulsar with a BH. The paradoxical fact that the first pulsar discovered in the SMC proved to be in a binary system can be naturally understood if its companion actually is a 10-30 solar mass black hole. We illustrate this fact by the numerical calculation of evolution of radio pulsars after a star formation burst.

Description: We calculate the structure of an effectively optically thin and geometrically thin accretion disk in the Kerr geometry, including electron-positron pairs. We show that the properties of the disk solutions are strongly dependent on the angular momentum of the central black hole. We find that close to a rapidly rotating hole there can be an appreciable pair density even for modest accretion rates. Pair critical accretion rates recently discovered in Newtonian disk models are also shown to be present in the general relativistic models, and we show that the geometrically thin disk approximations easily break down for rapidly rotating holes.

Description: We examine the hydrodynamic origin of relativistic outflows in active galactic nuclei (AGN). Specifically, we propose that the presence of a population of relativistic hadrons in the AGN 'central engine' and the associated neutron production suffices to produce outflows which under rather general conditions could be relativistic. The main such condition is that the size of the neutron production region be larger than the neutron flight path tau(sub n) approximately 3 x 10(exp 13) cm. This condition guarantees that the mean energy per particle in the proton fluid, resulting from the decay of the neutrons outside their production region, be greater than the proton rest mass. The expansion of this fluid can then lead naturally to a relativistic outflow by conversion of its internal energy to directed motion. We follow the development of such flows by solving the mass, energy as well as the kinetic equation for the proton gas in steady state, taking into account the source terms due to compute accurately the adiabatic index of the expanding gas, and in conjunction with Bernoulli's equation the detailed evolution of the bulk Lorentz factor. We further examine the role of large-scale magnetic fields in confining these outflows to produce the jets observed at larger scales.

Description: The topology of the magnetic field in the heliosheath is illustrated using plots of the field lines. It is shown that the Archimedean spiral inside the terminal shock is rotated back in the heliosheath into nested spirals that are advected in the direction of the interstellar wind. The 22-year solar magnetic cycle is imprinted onto these field lines in the form of unipolar magnetic envelopes surrounded by volumes of strongly mixed polarity. Each envelope is defined by the changing tilt of the heliospheric current sheet, which is in turn defined by the boundary of unipolar high-latitude regions on the Sun that shrink to the pole at solar maximum and expand to the equator at solar minimum. The detailed shape of the envelopes is regulated by the solar wind velocity structure in the heliosheath.

Description: SN 1987A hard X-ray continuum spectra obtained on 1987 October 29, 1988 April 9-10, and 1988 November 11 from balloon-flight measurments are presented. The spectra, spanning the energy range from 25 keV to 300 keV, have been analyzed using a detector response matrix inversion technique that converts the spectra form counts/s/sq cm keV to photons/s/sq cm keV allowing direct comparison with theoretical models. The results indicate that the bulk of the (56)Co is mixed moderately through the inner regions of the supernova envelope but they do not preclude the mixing of a small amount of the (56)Co farther out into the envelope necessary to account for the observed (56)Co line fluxes. The effect of the ratio (57)Co to (56)Co on the 1988 November 11 continuum spectrum is discussed.

Description: The UV line profile structure of high-ionization resonance lines found with the International Ultraviolet Explorer (IUE) in the brightest of four multiply imaged sources (images-A) in the candidate gravitational lens UM 425 = QSO 1120+019 indicates broad absorption line (BAL) structure. The deep-broad trough associated with the O IV line extends to velocities approiximately -12,000 km/s, and contains disrete features that suggest multicomponent velocity structure. This structure may include contributions from C IV absorption from the early-type galaxy that is believed to lens UM 425. A strong absorption feature in the blue wing of the Lyman-alpha lambda 1216 emission line may be a Lyman alpha absorption system at a Z(sub Ly alpha) = 1.437 +/- 0.003, or it may be formed by the superposition of the broad N V lambda lambda 1238, 1242 absorption trough on the extended blue emission wing of the QSO Lyman-alpha line. We obtained a redshift of Z(sub QSO) = 1.471 +/- 0.003 from Lyman-alpha lambda 1215, consistent with the redshift found by Meylan and Djorgovski in the optical. The Lyman-alpha line appears unusally weak due to the presence of N V lambda 1240 BAL absorption. A Lyman-limit absorption system at lambda 912 was not observed in the QSO rest frame. The detection of BAL structure in the other weaker ground-state resonance lines of N II (l) and S IV (l) was not found, suggesting these lines are formed in a region that is distinct from the BAL component. Detection of BAL structure in the other fainter images in this system with Hubble Space Telescope (HST) instrumentation, similar to structure observed here in image A, could provide evidence that UM 425 is a gravitational lens.

Description: We investigate the physical basis for the timescale of impulsive-phase, redshifted Lyman-alpha emission in stellar flares on the assumption that it is determined by energy losses in a nonthermal proton beam that is penetrating the chromosphere from above. The temporal evolution of ionization and heating in representative model chromospheres subjected to such beams is calculated. The treatment of 'stopping' of beam protons takes into account their interactions with (1) electrons bound in neutral hydrogen, (2) nuclei of neutral hydrogen, (3) free electrons, and (4) ambient thermal protons. We find that, for constant incident beam flux, the system attains an equilibrium with the beam energy input to the chromosphere balanced by radiative losses. In equilibrium, the beam penetration depth is constant, and erosion of the chromosphere ceases. If the redshifted, impulsive-phase stellar flare Lyman-alpha emission is produced by downstreaming hydrogen formed through charge exchange between beam protons and ambient hydrogen, then the emission should end when the beam no longer reaches neutral hydrogen. The durations of representative emission events calculated on this assumption range from 0.1 to 14 s. The stronger the beam, the shorter the timescale over which the redshifted Lyman-alpha emission can be observed.

Description: Through analysis of available optical spectrophotometric data and radio flux density measurements in the literature, it is demonstrated that a good correlation exists between the radio power and bolometric luminosity of the optically-selected OSOs in the Bright Quasar Sample (BOS) of Schmidt and Green (1983). We have recently used VLBI measurements of a sample of ultraluminous infrared galaxies to infer the likely existence of radio-quiet Active Galactic Nuclei (AGNs) deeply enshrouded in dust within their nuclei (Lonsdale, Smith, and Lonsdale 1993). We employ the radio-bolometric luminosity correlation for the BQS quasars to test whether these hypothetical buried AGNs can be energetically responsible for the observed far-infrared luminosities of the ultraluminous infrared galaxies. The ultraluminous infrared galaxies are shown to follow the same relation between radio core power and bolometric luminosity as the radio-quiet QSOs, suggesting that buried AGNs can account for essentially all the observed infrared luminosity, and raising the possibility that any starburst which may be in progress may not be energetically dominant. The broader implications of the radio-optical correlation in quasars for AGNs and luminous infrared galaxy models and the use of radio astronomy as a probe of the central powerhouse in radio quiet AGNs and luminous infrared galaxies are briefly discussed.

Description: The proposed design and construction of the Fizeau astrometric mapping explorer (FAME), a small astrometric instrument for use on an artificial satellite, is reported on. The instrument and spacecraft are designed to slowly spin and will repeatedly scan great circles on the sky so that, over a period of time, it will cover the complete sky and repeat in a manner similar to that of the Hipparcos satellite. The instrument will use the two fixed dilute aperture telescopes to measure a fixed angle between stars and detect the positions, magnitude and color of all stars crossing its field of view to a visual magnitude of approximately 15 mag. The aim of the instrument is to obtain a catalog of positions, proper motions and parallaxes of all stars down to about 15 mag, with a magnitude dependent accuracy of positions of 20 to 800 micro arcsec, proper motions of 20 to 800 micro arcsec per year and parallaxes of 20 to 800 micro arcsec.

Description: The concepts related to the operation and design of the global astrometric interferometer for astrophysics (GAIA) bring together solutions chosen for the astrometry satellite and interferometric techniques. Like the Hipparcos satellite, GAIA is a continuously scanning instrument for which the integration time on any observed object is limited by the field of view of the detector. If a final astrometric accuracy of 10 microarcsec is aimed at, a field of 1 deg in diameter is needed. A design is presented for the proposed 2.6 m baseline Fizeau interferometer with two 40 cm apertures and overall dimensions compatible with the size of the Ariane 5 payload shroud. It has a 0.9 deg diffraction limited field of view. The response of the optical system to small perturbations on each optical element is given in terms of the fringe visibility, which is shown to be dependent on the sub-aperture spot separation. The robustness of the design to thermal, mechanical and manufacturing errors is discussed. The unavoidable distortion present in wide field optical systems is analyzed in terms of displacement of the interference fringes.

Description: A direct link to an extragalactic reference system is considered as being a principle aim of the global astrometric interferometer for astrophysics (GAIA) mission. The data available from an extragalactic data base and a quasi stellar object (QSO) catalog were used to obtain an estimation of the number of QSO link candidates. The quality of presently available data and the expected accuracy of the extragalactic link are discussed. It is concluded that at least 150 QSO's must be observed by GAIA in order to guarantee an accuracy of better than 1 microarcsec/year for the link. New observations will be needed before the GAIA launch in order to reduce uncertainties in the positions, magnitudes and redshifts for some known quasars. The variability of QSO's with magnitudes near the GAIA observation limit can raise a potential problem. The motions of nearby QSO's are expected to be much smaller than 2 microarcsec/year, and therefore, will not affect the accuracy of the link in the proposed GAIA mission.

Description: We report on large-scale ab initio multiconfiguration Hartree-Fock calculations for the UV0.01 multiplet, 2s(sup 2)2p(sup 2)P(sub J) - 2s2p(sup 2 4)P(sub J prime), in N III. The resulting transition probabilities agree very well with recent semiempirical calculations, and the lifetimes for two of the three upper levels agree with experiments. The deviation for the third level is discussed. Comparisons made with the highest quality IUE echelle spectra available -- those of RR Tel and V1016 Cyg (both photoionized sources with electron densities below 10(exp 8)/cu cm) -- show that computed branching ratios of lines sharing a common upper level are in agreement with observations to within uncertainties of +/- 10%. High-quality solar limb data or stellar data from the Hubble Space Telescope (HST) could, in principle, be used to determine whether the theoretical or measured lifetimes for the discrepant level are in error. Unfortunately, stellar data for high-density plasmas (N(sub e) greater than 10(exp 11)/cu cm are needed) do not yet exist, and existing solar data lack the photometric precision to address this problem.

Description: Thermal drag, a variant of the Yarkovsky effect, may act on small asteroids with sizes from a few meters to a few tens of meters. Yarkovsky thermal drag comes from an asteroid's absorbing sunlight in the visible and reradiating it in the infrared. Since the infrared photons have momentum, by action-reaction, they kick the asteroid when they leave its surface. The reradiation, which is asymmetric in latitude over the asteroid, gives a net force along the asteroid's pole. Due to the asteroid's thermal inertia, averaging this force over one orbital period produces a net drag if the spin axis has a component in the orbital plane. Thermal drag tends to circularize orbits. It can increase or decrease orbital inclinations. An object whose spin axis points in random directions over its lifetime displays little change in orbital inclination. Thermal drag appears to have little to do with the delivery of chondrites from the asteroid belt; the thermal drag timescale (10(exp 8) years for meter-sizzed objects) is long compared with their cosmic ray exposure ages, and aphelia in the asteroid belt are not expected for mature thermal drag orbits. However, Yarkovsky thermal drag may act on the recently discovered near-Earth asteroids, which have radii of 10-30 m. Asteroid 1992 DA, for instance, might have its orbit shrunk by 0.1 AU in 3 x 10(exp 7) years, removing it from an Earth-crossing orbit. The near-Earth asteroids also tend to have small to moderate orbital eccentricities, as expected for highly evolved thermal drag objects. However, the time needed to bring them in from the asteroid belt (about 10(exp 9) years) is long compared with the collisonal and dynamical lifetimes (both about 10(exp 8) years) for Earth-crossing objects, arguing against their emplacement by thermal drag.

Description: High-contrast peaks in the cosmic microwave background (CMB) anisotropy can appear as unresolved sources to observers. We fit simluated CMB maps generated with a cold dark matter model to a set of unresolved features at instrumental resolution 0.5 deg-1.5 deg to derive the integral number density per steradian n (greater than absolute value of T) of features brighter than threshold temperature absolute value of T and compare the results to recent experiments. A typical medium-scale experiment observing 0.001 sr at 0.5 deg resolution would expect to observe one feature brighter than 85 micro-K after convolution with the beam profile, with less than 5% probability to observe a source brighter than 150 micro-K. Increasing the power-law index of primordial density perturbations n from 1 to 1.5 raises these temperature limits absolute value of T by a factor of 2. The MSAM features are in agreement with standard cold dark matter models and are not necessarily evidence for processes beyond the standard model.

Description: We have simulated full-sky maps of the cosmic microwave background (CMB) anisotropy expected from cold dark matter (CDM) models at 0.5 deg and 1.0 deg angular resolution. Statistical properties of the maps are presented as a function of sky coverage, angular resolution, and instrument noise, and the implications of these results for observability of the Doppler peak are discussed. The rms fluctuations in a map are not a particularly robust probe of the existence of a Doppler peak; however, a full correlation analysis can provide reasonable sensitivity. We find that sensitivity to the Doppler peak depends primarily on the fraction of sky covered, and only secondarily on the angular resolution and noise level. Color plates of the simulated maps are presented to illustrate the anisotropies.

Description: Soft gamma repeaters are characterized by recurrent activity consisting of short-duration outbursts of high-energy emission that is typically of temperature less than 40 keV. One recent model of repeaters is that they originate in the environs of neutron stars with superstrong magnetic fields, perhaps greater than 10(exp 14) G. In such fields, the exotic process of magnetic photon splitting gamma yields gamma gamma acts very effectively to reprocess gamma-ray radiation down to hard X-ray energies. In this Letter, the action of photon splitting is considered in some detail, via the solution of photon kinetic equations, determining how it limits the hardness of emission in strongly magnetized repeaters, and thereby obtaining observational constraints to the field in SGR 1806-20.

Description: The COMPTEL sky map of the 1.8 MeV line emission from Al-26 shows an extended diffuse distribution along the Galactic plane with a peculiar large scale asymmetry about the Galactic center (GC) and a clumpy structure with several noticeable hot spots. The most prominent hot spot at the GC appears shifted to positive longitude by about 2 deg. Nearby supernova remnants or Wolf-Rayet stars are plausible explanations for individual hot spots such as in the Vela region. We show that the global asymmetry and most hot spots can be understood by a more general model: the Al-26 sources are confined in the spiral arms of our Galaxy, and coagent star formation is responsible for additional clumpiness. The off-centered GC hot spot is probably due to the observed asymmetrically distributed circumnuclear molecular gas interacting with a central stellar bar, whose existence has been inferred from IR observations.

Description: Two dimensional spectral types for each of the stars observed in the global astrometric interferometer for astrophysics (GAIA) mission would provide additional information for the galactic structure and stellar evolution studies, as well as helping in the identification of unusual objects and populations. The classification of the large quantity generated spectra requires that automated techniques are implemented. Approaches for the automatic classification are reviewed, and a metric-distance method is discussed. In tests, the metric-distance method produced spectral types with mean errors comparable to those of human classifiers working at similar resolution. Data and equipment requirements for an automated classification survey, are discussed. A program of auxiliary observations is proposed to yield spectral types and radial velocities for the GAIA-observed stars.

Description: We consider the evolution and heating of dust embedded in the hot interstellar medium of isolated elliptical galaxies. We first construct a new set of galaxy models spanning a decade in luminosity which we use to study the evolution of dust as it is ejected by stellar sources, merges and flows with the interstellar gas, and is sputtered away due to the presence of the hot gas. We find that although grains can flow a considerable distance from the parent star in some cases before being sputtering away, the grain size distribution at a given location is accurately determined by assuming in situ sputtering of dust grains as they are ejected by the parent star. We find that dust heating is dominated by absorption of ambient starlight. Heating due to collisions with energetic electrons in the hot gas and absorption of thermal X-rays is smaller by orders of magnitude. We also find that for the largest galaxy considered (L(sub B) approximately 10 (exp 11) solar luminosity) the energy lost from the hot gas due to electronic collisions with the dust is about an order of magnitude less than that lost due to thermal emission in X-rays. The gas in smaller galaxies loses even less of its energy via heating of grains. In deriving this result, we find that the inclusion of grain sputtering is a crucial ingredient of the calculation. The presence of dust grains is therefore energetically unimportant for elliptical galaxy cooling flow models.

Description: Near-infrared images of the Galactic bulge at 1.25, 2.2, 3.5, and 4.9 microns obtained by the Diffuse Infrared Background Experiment (DIRBE) onboard the Cosmic Background Explorer (COBE) satellite are used to characterize its morphology and to determine its infrared luminosity and mass. Earlier analysis of the DIRBE observations (Weiland et al. 1994) provided supporting evidence for the claim made by Blitz & Spergel (1991) that the bulge is bar-shaped with its near end in the first Galactic quadrant. Adopting various triaxial analytical functions to represent the volume emissivity of the source, we confirm the barlike nature of the bulge and show that triaxial Gaussian-type functions provide a better fit to the data than other classes of functions, including an axisymmetric spheroid. The introduction of a `boxy' geometry, such as the one used by Kent, Dame, & Fazio (1991) improves the fit to the data. Our results show that the bar is rotated in the plane with its near side in the first Galactic quadrant creating an angle of 20 deg +/- 10 deg between its major axis and the line of sight to the Galactic center. Typical axis ratios of the bar are (1:0.33 +/- 0.11:0.23 +/- 0.08), resembling the geometry of prolate spheroids. There is no statistically significant evidence for an out-of-plane tilt of the bar at 2.2 microns, and marginal evidence for a tilt of approximately equal 2 deg at 4.9 microns. The introduction of a roll around the intrinsic major axis of the bulge improves the `boxy' appearance of some functions. A simple integration of the observed projected intensity of the bulge gives a bulge luminosity of 1.2 x 10(exp 9), 4.1 x 10(exp 8), 2.3 x 10(exp 8), and 4.3 x 10(exp 7) solar luminosity, respectively, at 1.25, 2.2, 3.5, and 4.9 microns wavelength for a Galactocentric distance of 8.5 kpc. The 2.2 microns luminosity function of the bulge population in the direction of Baade's window yields a bolometric luminosity of L(sub bol) = 5.3 x 10(exp 9) solar luminosity. Stellar evolutionary models relate this luminosity to the number of main-sequence progenitor stars that currently populate the red giant branch. Combined with the recent determination of the main-sequence turnoff mass for the bulge by the Hubble Space Telescope (Holtzman et al. 1993) we derive a photometrically determined bulge mass of approximately equal to 1.3 x 10(exp 10) solar mass for a Salpeter initial mass function extended down to 0.1 solar mass.

Description: We report the results from observations of National Radio Astronomy Observatory (NRAO) 140 carried out in the X-ray band using ROSAT and Advanced Satellite for Cosmology and Astrophysics (ASCA) in 1992 and 1994. We find the source to be absorbed by an effective hydrogen column density of approximately 3 x 10(exp 21) atoms/sq cm during both epochs, consistent with the combined atomic and molecular column inferred from radio measurements along this line of sight through the Perseus cloud complex. We compare these results with previous X-ray absorption measurements and briefly comment upon the origin of the excess absorption which has been seen toward this source. We find the ASCA spectrum of NRAO 140 is well described by a power law of energy index alpha = 0.73 +/- 0.03 and also yields the tightest constraint to date on Fe K-shell emission, with 90% confidence upper limits of 38 and 31 eV for a narrow line at a rest frame energy of 6.4 and 6.7 keV respectively. This, along with a lack of hardening towards higher energies, suggests that either NRAO 140 is devoid of cold reprocessing material, the reprocessing material has a geometry in which the imprinted features are weak, and/or the X-ray emission is relativistically beamed towards us. We also report the detection of a serendipitous source in both the Advanced Satellite for Cosmology and Astrophysics Gas Imaging Spectrometer (ASCA GIS) and ROSAT fields of view. We identify this source as the binary star system IX Per and find its spectrum to be well fitted by a two-temperature Raymond-Smith plasma.

Description: We have produced two-dimensional maps of the intensity ratio, Q(sub 60), of 60 micron infrared to 20 cm radio continuum emission, for a set of 25 nearby galaxies, mostly spirals. The ratio maps were obtained from infrared images made using IRAS data with the maximum correlation method, and radio images made using VLA data. Before taking the ratio, the radio images were processed so as to have the same resolution properties as the infrared images; the final spatial resolution in all cases is approximately 1 min, corresponding to 1 - 2 kpc for most galaxies. This resolution represents a significant improvement over previous studies. Our new high-resolution maps confirm the slow decrease of Q(sub 60) with increasing radial distance from the nucleus, but show additional structure which is probably associated with separate sites of active star formation in the spiral arms. The maps show Q(sub 60) to be more closely related to infrared surface brightness than to the radial distance r in the galaxy disk. We note also that the Q(sub 60) gradients are absent (or at least reduced) for the edge-on galaxies, a property which can be attributed to the dilution of contrast due to the averaging of the additional structure along the line of sight. The results are all in qualitative agreement with the suggestion that the radio image represents a smeared version of the infrared image, as would be expected on the basis of current models in which the infrared-radio correlation is driven by the formation of massive stars, and the intensity distribution of radio emission is smeared as a result of the propagation of energetic electrons accelerated during the supernova phase.

Description: The U-B color has been suggested as a predictor of the presence of water of hydration on asteroids. Photometry from the Eight-Color Asteroid Survey (ECAS) was used to test this concept. An overlap in U-B color prevents this magnitude difference from distinguishing between surface material that was thermally processed at higher temperatures and surface material that was aqueously altered. Two tests of the presence of water of hydration using visible spectral region photometry failed to flag those few higher albedo M- and E-class asteroids having photometry that shows a 3.0-micrometers water of hydration absorption. These asteroids probably contain little or no oxidized iron in their surface material.

Description: The Monte Carlo technique of simulating diffusive particle acceleration at shocks has made spectral predictions that compare extremely well with particle distributions observed at the quasi-parallel region of the earth's bow shock. The current extension of this work to compare simulation predictions with particle spectra at oblique interplanetary shocks has required the inclusion of significant cross-field diffusion (strong scattering) in the simulation technique, since oblique shocks are intrinsically inefficient in the limit of weak scattering. In this paper, we present results from the method we have developed for the inclusion of cross-field diffusion in our simulations, namely model predictions of particle spectra downstream of oblique subluminal shocks. While the high-energy spectral index is independent of the shock obliquity and the strength of the scattering, the latter is observed to profoundly influence the efficiency of injection of cosmic rays into the acceleration process.

Description: A standard problem in gamma-ray astronomy data analysis is the decomposition of a set of observed counts, described by Poisson statistics, according to a given multicomponent linear model, with underlying physical count rates or fluxes which are to be estimated from the data. Despite its conceptual simplicity, the linear least-squares (LLSQ) method for solving this problem has generally been limited to situations in which the number n(sub i) of counts in each bin i is not too small, conventionally more than 5-30. It seems to be widely believed that the failure of the LLSQ method for small counts is due to the failure of the Poisson distribution to be even approximately normal for small numbers. The cause is more accurately the strong anticorrelation between the data and the wieghts w(sub i) in the weighted LLSQ method when square root of n(sub i) instead of square root of bar-n(sub i) is used to approximate the uncertainties, sigma(sub i), in the data, where bar-n(sub i) = E(n(sub i)), the expected value of N(sub i). We show in an appendix that, avoiding this approximation, the correct equations for the Poisson LLSQ (PLLSQ) problems are actually identical to those for the maximum likelihood estimate using the exact Poisson distribution. We apply the method to solve a problem in high-resolution gamma-ray spectroscopy for the JPL High-Resolution Gamma-Ray Spectrometer flown on HEAO 3. Systematic error in subtracting the strong, highly variable background encountered in the low-energy gamma-ray region can be significantly reduced by closely pairing source and background data in short segments. Significant results can be built up by weighted averaging of the net fluxes obtained from the subtraction of many individual source/background pairs. Extension of the approach to complex situations, with multiple cosmic sources and realistic background parameterizations, requires a means of efficiently fitting to data from single scans in the narrow (approximately = 1.2 keV, HEAO 3) energy channels of a Ge spectrometer, where the expected number of counts obtained per scan may be very low. Such an analysis system is discussed and compared to the method previously used.

Description: Four spectral emission features of the N(sub KK) = 4(sub 04) -3(sub 13) rotational transition of methylene (CH2) have been detected at signal levels 5-7 sigma above noise toward the hot core of the Orion-KL nebula and the molecular cloud in proximity to the continuum source W51 M. Specifically, in both sources we have resolved the F = 6-5, 5-4, and 4-3 hyperfine transitions of the J = 5-4 fine-structure levels and detected the blended hyperfine structure of the J = 4-3 fine structure levels. At the J = 3-2 fine-structure levels, we have observed new transitions of NS, a known interstellar molecule, which severely contaminates the search for CH2 hyperfine transitions. These new sensitive observations finally confirm the existence of interstellar CH2 which was tentatively reported by us some years ago.

Description: The star, Altair (A7 IV-V), is clearly shown to have Lyman-alpha emission of chromospheric origin, while no evidence is found for the Mg II emission reported in previous investigations. We present non-Local Thermodymanic Equilibrium (non-LTE) semiempirical models incorporating partial redistribution of the chromosphere of Altair that reproduce the observed Lyman-alpha emission and the Mg II resonance absorption at 2800 A. We unambiguously establihed that chromospheres exist at spectral types as early as A7 on the main sequence, and we also demonstrate that it very unlikely that the observed emission originates in a corotating expanding wind. This result represents a new challenge for chromospheric heating theories. It may indicate that both differential rotation and convection layers, at least near the equator, exist in this fast rotating (v sin i = 220 km/s) star.

Description: We use the two-point correlation function of the extrema points (peaks and valleys) in the Cosmic Background Explorer (COBE) Differential Microwave Radiometers (DMR) 2 year sky maps as a test for non-Gaussian temperature distribution in the cosmic microwave background anisotropy. A maximum-likelihood analysis compares the DMR data to n = 1 toy models whose random-phase spherical harmonic components a(sub lm) are drawn from either Gaussian, chi-square, or log-normal parent populations. The likelihood of the 53 GHz (A+B)/2 data is greatest for the exact Gaussian model. There is less than 10% chance that the non-Gaussian models tested describe the DMR data, limited primarily by type II errors in the statistical inference. The extrema correlation function is a stronger test for this class of non-Gaussian models than topological statistics such as the genus.

Description: During its full-sky survey, the Energetic Gamma Ray Experiment Telescope (EGRET) high-energy instrument aboard the Compton Gamma-Ray Observatory detected gamma-ray emission in the energy range above 30 MeV from a source identified as the blazar PKS 0420-014. This object was observed during two separate viewing periods in 1992 February/March and May/June. An intensity decrease above 100 MeV of a factor of at least 1.5 from a maximum of (5.0 +/- 1.4) x 10(exp -7) photons 1 sq cm/s was observed during that time interval indicating extensive variability. The photon spectrum in the range between 30 and 10,000 MeV at the time of the maximum intensity is well represented by a power law with an exponent of -1.9 +/- 0.3. Some similarities with other EGRET detected blazars are briefly discussed.

Description: Eleven tests were carried out in DLR's high enthalpy tunnel in Goettingen (HEG), Germany, for reservoir conditions ranging from 10 to 23 MJ/kg and a free stream Mach number of approximately 10. A blunted cone model with a cylindrical afterbody (sting) was investigated. To obtain information on the influence of defined parameters on the body back flow, especially of the reacting gas, the heat transfer rate along the body contour was measured with fast response surface thermocouples on the forebody and sensitive thin film heat transfer gauges on the base and sting of the model. Flow visualization with a holographic interferometry system was provided. Tests are described, and a preliminary interpretation of the observed flow effects are given.

Description: Get Away Special (GAS) payload G-093, also called VORTEX (VOrtex Ring Transit EXperiment), is an investigation of the propagation of a vortex ring through a liquid-gas interface in microgravity. This process results in the formation of one or more liquid droplets similar to earth based liquid atomization systems. In the absence of gravity, surface tension effects dominate the drop formation process. The Shuttle's microgravity environment allows the study of the same fluid atomization processes as using a larger drop size than is possible on Earth. This enables detailed experimental studies of the complex flow processes encountered in liquid atomization systems. With VORTEX, deformations in both the vortex ring and the fluid surface will be measured closely for the first time in a parameters range that accurately resembles liquid atomization. The experimental apparatus will record images of the interactions for analysis after the payload has been returned to earth. The current design of the VORTEX payload consists of a fluid test cell with a vortex ring generator, digital imaging system, laser illumination system, computer based controller, batteries for payload power, and an array of housekeeping and payload monitoring sensors. It is a self-contained experiment and will be flown on board the Space Shuttle in a 5 cubic feet GAS canister. The VORTEX Project is entirely run by students at the University of Michigan but is overseen by a faculty advisor acting as the payload customer and the contact person with NASA. This paper summarizes both the technical and programmatic aspects of the VORTEX Project.

Description: This paper focuses on the flight results of the Cryogenic Two-Phase Flight Experiment (CRYOTP), which was a Hitchhiker based experiment that flew on the space shuttle Columbia in March of 1994 (STS-62). CRYOTP tested two new technologies for advanced cryogenic thermal control; the Space Heat Pipe (SHP), which was a constant conductance cryogenic heat pipe, and the Brilliant Eyes Thermal Storage Unit (BETSU), which was a cryogenic phase-change thermal storage device. These two devices were tested independently during the mission. Analysis of the flight data indicated that the SHP was unable to start in either of two attempts, for reasons related to the fluid charge, parasitic heat leaks, and cryocooler capacity. The BETSU test article was successfully operated with more than 250 hours of on-orbit testing including several cooldown cycles and 56 freeze/thaw cycles. Some degradation was observed with the five tactical cryocoolers used as thermal sinks, and one of the cryocoolers failed completely after 331 hours of operation. Post-flight analysis indicated that this problem was most likely due to failure of an electrical controller internal to the unit.

Description: The Capillary Pumped Loop Flight Experiment (CAPL) employs a passive two-phase thermal control system that uses the latent heat of vaporization of ammonia to transfer heat over long distances. CAPL was designed as a prototype of the Earth Observing System (EOS) instrument thermal control systems. The purpose of the mission was to provide validation of the system performance in micro-gravity, prior to implementation on EOS. CAPL was flown on STS-60 in February, 1994, with some unexpected results related to gravitational effects on two-phase systems. Flight test results and post flight investigations will be addressed, along with a brief description of the experiment design.

Description: This paper presents an overview of the Cryogenic Test Bed (CTB) experiments including experiment results, integration techniques used, and lessons learned during integration, test and flight phases of the Cryogenic Heat Pipe Flight Experiment (STS-53) and the Cryogenic Two Phase Flight Experiment (OAST-2, STS-62). We will also discuss the Cryogenic Flexible Diode Heat Pipe (CRYOFD) experiment which will fly in the 1996/97 time frame and the fourth flight of the CTB which will fly in the 1997/98 time frame. The two missions tested two oxygen axially grooved heat pipes, a nitrogen fibrous wick heat pipe and a 2-methylpentane phase change material thermal storage unit. Techniques were found for solving problems with vibration from the cryo-collers transmitted through the compressors and the cold heads, and mounting the heat pipe without introducing parasitic heat leaks. A thermally conductive interface material was selected that would meet the requirements and perform over the temperature range of 55 to 300 K. Problems are discussed with the bi-metallic thermostats used for heater circuit protection and the S-Glass suspension straps originally used to secure the BETSU PCM in the CRYOTP mission. Flight results will be compared to 1-g test results and differences will be discussed.

Description: Comets are recognized as among the most scientifically important objects in the solar system. They are presumed relics of the early primitive material in the solar nebula and are believed to have provided a general enrichment of volatiles to the inner solar system. The Cometary Coma Chemical Composition (C4) Mission, a proposed Discovery-Class Mission, will analyze materials released into the coma, providing information leading to the understanding of the chemical composition and make-up of the cometary nucleus. As one of two scientific instruments in the C4 spacecraft, an advanced and streamlined version of the Cometary Ice and Dust Experiment (CIDEX), a mini-CIDEX, will employ an X-Ray Fluorescence (XRF) spectrometer to determine bulk elemental composition of cometary dust grains and a Gas Chromatograph/Ion Mobility Spectrometer (GC/IMS) for determination of the molecular composition of dust and ices following stepwise pyrolysis and combustion. A description of the mini-CIDEX IMS will be provided as well as data from analyses conducted using the mini-CIDEX breadboard instrument.

Description: With the recent focus on the needs of design and applications CFD, research groups have begun to address the traditional bottlenecks of grid generation and surface modeling. Now, a host of emerging technologies promise to shortcut or dramatically simplify the simulation process. This paper discusses the current status of these emerging technologies. It will argue that some tools are already available which can have positive impact on portions of the design cycle. However, in most cases, these tools need to be integrated into specific engineering systems and process cycles to be used effectively. The rapidly maturing status of unstructured and Cartesian approaches for inviscid simulations makes suggests the possibility of highly automated Euler-boundary layer simulations with application to loads estimation and even preliminary design. Similarly, technology is available to link block structured mesh generation algorithms with topology libraries to avoid tedious re-meshing of topologically similar configurations. Work in algorithmic based auto-blocking suggests that domain decomposition and point placement operations in multi-block mesh generation may be properly posed as problems in Computational Geometry, and following this approach may lead to robust algorithmic processes for automatic mesh generation.

Description: The last decade has witnessed a vigorous and sustained research effort on unstructured methods for computational fluid dynamics. Unstructured mesh generators and flow solvers have evolved to the point where they are now in use for design purposes throughout the aerospace industry. In this paper we survey the various mesh types, structured as well as unstructured, and examine their relative strengths and weaknesses. We argue that unstructured methodology does offer the best prospect for the next generation of computational fluid dynamics algorithms.

Description: A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree data structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: A gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment and other accepted computational results for a series of low and moderate Reynolds number flows.

Description: Grid related issues of the Chimera overset grid method are discussed in the context of a method of solution and analysis of unsteady three-dimensional viscous flows. The state of maturity of the various pieces of support software required to use the approach is considered. Current limitations of the approach are identified.

Description: This paper presents 'A view from the trenches' on CFD grid generation from a Pratt & Whitney perspective. We anticipate that other organizations have similar views. We focus on customer expectations and the consequent requirements. We enunciate a vision for grid generation, discuss issues that developers must recognize.

Description: This paper presents a perspective on the requirements that Computational Fluid Dynamics (CFD) technology must meet for its effective use in aerospace design. General observations are made on current aerospace design practices and deficiencies are noted that must be rectified for the U.S. aerospace industry to maintain its leadership position in the global marketplace. In order to rectify deficiencies, industry is transitioning to an integrated product and process development (IPPD) environment and design processes are undergoing radical changes. The role of CFD in producing data that design teams need to support flight vehicle development is briefly discussed. An overview of the current state of the art in CFD is given to provide an assessment of strengths and weaknesses of the variety of methods currently available, or under development, to produce aerodynamic data. Effectiveness requirements are examined from a customer/supplier view point with design team as customer and CFD practitioner as supplier. Partnership between the design team and CFD team is identified as an essential requirement for effective use of CFD. Rapid turnaround, reliable accuracy, and affordability are offered as three key requirements that CFD community must address if CFD is to play its rightful role in supporting the IPPD design environment needed to produce high quality yet affordable designs.

Description: The efforts in geometry modeling and grid generation at the NASA Lewis Research Center, as applied to the computational fluid dynamic (CFD) analysis of aeropropulsion systems, are presented. The efforts are mainly characterized by a focus on the analysis of components of an aeropropulsion system, which involve turbulent viscous flow with heat transfer and chemistry. Thus, this discussion will follow that characterization and will sequence through the components of typical propulsion systems consisting of inlets, compressors, combustors, turbines, and nozzles. For each component, some applications of CFD analysis will be presented to show how CFD is used to compute the desired performance information, how geometry modeling and grid generation are performed, and what issues have developed related to geometry modeling and grid generation. The discussion will illustrate the following needs related to geometry modeling and grid generation as observed in aeropropulsion analysis: (1) accurate and efficient resolution of turbulent viscous and chemically-reacting flowfields; (2) easy-to-use interfaces with CAD data for automated grid generation about complex geometries; and (3) automated batch grid generation software for use with design and optimization software.

Description: This viewgraph presentation discusses key features of current combustion system CFD modeling capabilities at GE Aircraft Engines provided by the CONCERT code; CONCERT development history; modeling applied for designing engine combustion systems; modeling applied to improve fundamental understanding; CONCERT3D results for current production combustors; CONCERT3D model of NASA/GE E3 combustor; HYBRID CONCERT CFD/Monte-Carlo modeling approach; and future modeling directions.

Description: A procedure was developed to improve the turn-around time for computational fluid dynamics (CFD) simulations of an inlet-bleed problem involving oblique shock-wave/boundary-layer interactions on a flat plate with bleed into a plenum through one or more circular holes. This procedure is embodied in a preprocessor called AUTOMAT. With AUTOMAT, once data for the geometry and flow conditions have been specified (either interactively or via a namelist), it will automatically generate all input files needed to perform a three-dimensional Navier-Stokes simulation of the prescribed inlet-bleed problem by using the PEGASUS and OVERFLOW codes. The input files automatically generated by AUTOMAT include those for the grid system and those for the initial and boundary conditions. The grid systems automatically generated by AUTOMAT are multi-block structured grids of the overlapping type. Results obtained by using AUTOMAT are presented to illustrate its capability.

Description: The role of Computational Fluid Dynamics (CFD) at Ames Research Center has expanded to address a broad range of aeronautical problems, including wind tunnel support, flight test support, design, and analysis. Balancing the requirements of each new problem against the available resources - software, hardware, time, and expertise - is critical to the effective use of CFD. Several case studies of recent applications highlight the depth of CFD capability at Ames, the tradeoffs involved in various approaches, and lessons learned in the use of CFD as an engineering tool.

Description: When computing the flow around complex three dimensional configurations, the generation of the mesh is the most time consuming part of any calculation. With some meshing technologies this can take of the order of a man month or more. The requirement for a number of design iterations coupled with ever decreasing time allocated for design leads to the need for a significant acceleration of this process. Of the two competing approaches, block-structured and unstructured, only the unstructured approach will allow fully automatic mesh generation directly from a CAD model. Using this approach coupled with the techniques described in this paper, it is possible to reduce the mesh generation time from man months to a few hours on a workstation. The desire to closely couple a CFD code with a design or optimization algorithm requires that the changes to the geometry be performed quickly and in a smooth manner. This need for smoothness necessitates the use of Bezier polynomials in place of the more usual NURBS or cubic splines. A two dimensional Bezier polynomial based design system is described.

Description: The results of a first order perfect gas correction for the effects of the boundary layer formation within expansion tubes with nozzles are presented. The analytical model developed to describe the boundary layer formation within the expansion tube and an expansion nozzle located at the end of the acceleration tube is based on the Karman integral equations. The results of this analytical model are compared with experimental data from an expansion diffuser. The model provides a useful tool for the preliminary design of nozzles for such facilities.

Description: In this paper the capabilities of an automated CFD system which is currently available at NLR are demonstrated. Transonic flow around the AS28G wing/body configuration and hypersonic flow through a generic three-dimensional mixed-compression airbreathing inlet are simulated. An assessment of the level of automation of the current CFD-system is made. The problem-turnaround time lies within the order of a week for both applications.

Description: Detailed simulations of viscous flows in complicated geometries pose a significant challenge to current capabilities of Computational Fluid Dynamics (CFD). To enable routine application of CFD to this class of problems, advanced methodologies are required that employ (a) automated grid generation, (b) adaptivity, (c) accurate discretizations and efficient solvers, and (d) advanced software techniques. Each of these ingredients contributes to increased accuracy, efficiency (in terms of human effort and computer time), and/or reliability of CFD software. In the long run, methodologies employing structured grid systems will remain a viable choice for routine simulation of flows in complex geometries only if genuinely automatic grid generation techniques for structured grids can be developed and if adaptivity is employed more routinely. More research in both these areas is urgently needed.

Description: A discussion of the strengths and weaknesses of overset composite grid and solution technology is given, along with a sampling of current work in the area. Major trends are identified, and the observation is made that generalized and hybridized overset methods provide a natural framework for combining disparate mesh types and physics models. Because of this, the author concludes that overset methods will be the foundation for the general purpose computational fluid dynamics programs of the future.

Description: The status of CFD methods based on the use of block-structured grids for analyzing viscous flows over complex configurations is examined. The objective of the present study is to make a realistic assessment of the usability of such grids for routine computations typically encountered in the aerospace industry. It is recognized at the very outset that the total turnaround time, from the moment the configuration is identified until the computational results have been obtained and postprocessed, is more important than just the computational time. Pertinent examples will be cited to demonstrate the feasibility of solving flow over practical configurations of current interest on block-structured grids.

Description: This viewgraph presentation discusses turbulence modeling requirements, development philosophy, and approach; two major areas of concentration (high speed and low speed turbulence modeling); high speed turbulence modeling; compressibility effects; turbulence models adapted to USA code; M = 9.2 flat plate flow; Mach 7.05 flow over axisymmetric flare; Mach 8.6 flow over cold wall edge; low speed turbulence modeling; turbulence models being assessed; turbulence model deck structure and integration with Navier-Stokes solver; nonlinear algebraic-stress model; rotation modified k-epsilon model; and Reynolds stress model.

Description: This viewgraph presentation covers gas turbine combustor flow physics, turbulence model investigations, turbulent combustion modeling, and present status and future needs.

Description: This viewgraph presentation discusses geometry and flow configuration, effect of y+ on heat transfer computations, standard and extended k-epsilon turbulence model results with wall function, low-Re model results (the Lam-Bremhorst model without wall function), a criterion for flow reversal in a radially rotating square duct, and a summary.

Description: This viewgraph presentation demonstrates that computationally efficient k-l and k-kl turbulence models have been developed and implemented at Lockheed Fort Worth Company. Many years of experience have been gained applying two equation turbulence models to complex three-dimensional flows for design and analysis.

Description: The objective of this viewgraph presentation is to evaluate turbulence models for integrated aircraft components such as the forebody, wing, inlet, diffuser, nozzle, and afterbody. The one-equation models have replaced the algebraic models as the baseline turbulence models. The Spalart-Allmaras one-equation model consistently performs better than the Baldwin-Barth model, particularly in the log-layer and free shear layers. Also, the Sparlart-Allmaras model is not grid dependent like the Baldwin-Barth model. No general turbulence model exists for all engineering applications. The Spalart-Allmaras one-equation model and the Chien k-epsilon models are the preferred turbulence models. Although the two-equation models often better predict the flow field, they may take from two to five times the CPU time. Future directions are in further benchmarking the Menter blended k-w/k-epsilon and algorithmic improvements to reduce CPU time of the two-equation model.

Description: This viewgraph presentation summarizes some CFD experience at GE Aircraft Engines for flows in the primary gaspath of a gas turbine engine and in turbine blade cooling passages. It is concluded that application of the standard k-epsilon turbulence model with wall functions is not adequate for accurate CFD simulation of aerodynamic performance and heat transfer in the primary gas path of a gas turbine engine. New models are required in the near-wall region which include more physics than wall functions. The two-layer modeling approach appears attractive because of its computational complexity. In addition, improved CFD simulation of film cooling and turbine blade internal cooling passages will require anisotropic turbulence models. New turbulence models must be practical in order to have a significant impact on the engine design process. A coordinated turbulence modeling effort between NASA centers would be beneficial to the gas turbine industry.

Description: Program goals at the Center for Modeling of Turbulence and Transition (CMOTT), NASA Lewis Research Center, are (1) to develop reliable turbulence (including bypass transition) and combustion models for complex flows in propulsion systems and (2) to integrate developed models into deliverable CFD tools for propulsion systems in collaboration with industry. This viewgraph presentation covers the following topics: development of turbulence and combustion models; collaboration with industry and technology transfer; isotropic eddy viscosity models; algebraic Reynolds stress models; scalar turbulence models; second order closure models; multiple scale k-epsilon models; and PDF modeling of turbulent reacting flows.

Description: This study is a further investigation of space-exposed samples recovered from the LDEF satellite and the Franco-Russian 'Aragatz' dust collection experiment on the Mir Space Station. Impact craters with diameters ranging from 1 to 900 micron were found on the retrieved samples. Elemental analysis of residues found in the impact craters was carried out using Energy Dispersive X-ray spectrometry (EDX). The analyses show evidence of micrometeoroid and orbital debris origins for the impacts. The proportions of these two components vary according to particle size and experimental position with respect to the leading edge of the spacecraft. On the LDEF leading edge 17 percent of the impacts were apparently caused by micrometeoroids and 11 percent by debris; on the LDEF trailing edge 23 percent of the impacts are apparently caused by micrometeoroids and 4 percent consist of debris particles - mostly larger than 3 micron in diameter - in elliptical orbits around the Earth. For Mir, the analyses indicate that micrometeoroids form 23 percent of impacts and debris 9 percent. However, we note that 60-70 percent of the craters are unidentifiable, so the definitive proportions of natural v. man-made particles are yet to be determined. Experiments carried out using a light gas gun to accelerate glass spheres and fragments demonstrate the influence of particle shape on crater morphology. The experiments also show that it is more difficult to analyze the residues produced by an irregular fragment than those produced by a spherical projectile. If the particle is travelling above a certain velocity, it vaporizes upon impact and no residues are left. Simulation experiments carried out with an electrostatic accelerator indicate that this limit is about 14 km/s for Fe particles impacting Al targets. This chemical analysis cut-off may bias interpretations of the relative populations of meteoroid and orbital debris. Oblique impacts and multiple foil detectors provide a higher likelihood of detection of residues as the velocities involved are lower.

Description: Approximately 1000 impact craters on the Chemistry of Meteoroid Experiment (CME) have been analyzed by means of Scanning Electron Microscopy (SEM) and Energy Dispersive X-ray Analysis (EDXA) to determine the compositional make-up of projectile residues. This report completes our systematic survey of gold and aluminum surfaces exposed at the trailing-edge (A03) and forward-facing (A11) LDEF sites, respectively. The major categories for the projectile residues were (1) natural, with diverse subgroups such as chondritic, monomineralic silicates, and sulfides, and (2) man made, that were classified into aluminum (metallic or oxide) and miscellaneous materials (such as stainless steel, paint flakes, etc). On CME gold collectors on LDEF's trailing edge approximately 11 percent of all craters greater than 100 micron in diameter were due to man-made debris, the majority (8.6 percent) caused by pure aluminum, approximately 31.4 percent were due to cosmic dust, while the remaining 58 percent were indeterminate via the analytical techniques utilized in this study. The aluminum surfaces located at the A11 forward-facing site did not permit analysis of aluminum impactors, but approximately 9.4 percent of all craters were demonstratably caused by miscellaneous debris materials and approximately 39.2 percent were the result of natural particles, leaving approximately 50 percent which were indeterminate. Model considerations and calculations are presented that focus on the crater-production rates for features greater than 100 micron in diameter, and on assigning the intermediate crater population to man-made or natural particles. An enhancement factor of 6 in the crater-production rate of natural impactors for the 'forward-facing' versus the 'trailing-edge' CME collectors was found to best explain all observations (i.e., total crater number(s), as well as their computational characteristics). Enhancement factors of 10 and 4 are either too high or too low. It is also suggested that approximately 45 percent of all craters greater than 100 micron in diameter are caused by man-made impactors on the A11 surfaces. This makes the production rate for craters greater than 100 micron in diameter, resulting from orbital debris, a factor of 40 higher on the forward-facing sides as opposed to the trailing-edge direction.

Description: Examination of LDEF's various surfaces shows numerous craters and holes due to hypervelocity impacts of meteoroids and man-made orbital debris. In this paper, the crater numbers as reported by Humes have been analyzed in an effort to understand the orbital debris and natural meteoroid environment in LEO. To determine the fraction of man-made to natural impacts, the side to top ratio of impacts and results of the Chemistry of Micrometeoroids Experiment are used. For craters in the 100 micron to 500 micron size range, about 25 percent to 30 percent of the impacts on the forward-facing surfaces and about 10 percent of the impacts on the trailing surfaces were estimated due to man-made orbital debris. A technique has been developed to convert crater numbers to particle fluxes, taking the fact into account that the distributions of impact velocity and incidence angle vary over the different surfaces of LDEF, as well as the ratio of the surface area flux to the cross-sectional area flux. Applying this technique, Humes' data concerning craters with limiting lip diameters of 100 micron, 200 micron and 500 micron have been converted into orbital debris and meteoroid fluxes ranging from about 20 micron to 200 micron particle diameter. The results exhibit good agreement with orbital debris model and meteoroid model. The converted meteoroid flux is slightly larger than Grun's model (by 40 to 70 percent). The converted orbital debris flux is slightly lower than Kessler's model for particle diameter smaller than about 30 micron and slightly larger than the model for particle diameter larger than about 40 micron. Taking also into account the IDE data point at about 0.8 micron particle diameter, it suggests to change the slope log (flux) versus log (diameter) of orbital debris flux in the 1 micron to 100 micron particle diameter range from 2.5 to 1.9.

Description: We discuss the Shannon-Wehrl entropy within the squeezing vocabulary for the cosmological and black hole particle production.

Description: An assessment of two unstructured methods is presented in this paper. A tetrahedral unstructured method USM3D, developed at NASA Langley Research Center is compared to a Cartesian unstructured method, SPLITFLOW, developed at Lockheed Fort Worth Company. USM3D is an upwind finite volume solver that accepts grids generated primarily from the Vgrid grid generator. SPLITFLOW combines an unstructured grid generator with an implicit flow solver in one package. Both methods are exercised on three test cases, a wing, and a wing body, and a fully expanded nozzle. The results for the first two runs are included here and compared to the structured grid method TEAM and to available test data. On each test case, the set up procedure are described, including any difficulties that were encountered. Detailed descriptions of the solvers are not included in this paper.

Description: The past ten years have seen steady progress in surface modeling procedures, and wholesale changes in grid generation technology. Today, it seems fair to state that a satisfactory grid can be developed to model nearly any configuration of interest. The issues at present focus on operational concerns such as cost and quality. Continuing evolution of the engineering process is placing new demands on the technologies of surface modeling and grid generation. In the evolution toward a multidisciplinary analysis-bascd design environment, methods developed for Computational Fluid Dynamics are finding acceptance in many additional applications. These two trends, the normal evolution of the process and a watershed shift toward concurrent and multidisciplinary analysis, will be considered in assessing current capabilities and needed technological improvements.

Description: This viewgraph presentation discusses what models are used in this package and what their advantages and disadvantages are, how the probability density function (PDF) model is implemented and the features of the program, and what can be expected in the future from the NASA Lewis PDF code.

Description: This viewgraph presentation discusses an extension of the probability density function (PDF) method to the modeling of spray flames to evaluate the limitations and capabilities of this method in the modeling of gas-turbine combustor flows. The comparisons show that the general features of the flowfield are correctly predicted by the present solution procedure. The present solution appears to provide a better representation of the temperature field, particularly, in the reverse-velocity zone. The overpredictions in the centerline velocity could be attributed to the following reasons: (1) the use of k-epsilon turbulence model is known to be less precise in highly swirling flows and (2) the swirl number used here is reported to be estimated rather than measured.

Description: Viewgraphs are presented on computation of turbulent combustion, governing equations, closure problem, PDF modeling of turbulent reactive flows, validation cases, current projects, and collaboration with industry and technology transfer.

Description: Viewgraphs are presented on modeling turbulent reacting flows, regimes of turbulent combustion, regimes of premixed and regimes of non-premixed turbulent combustion, chemical closure models, flamelet model, conditional moment closure (CMC), NO(x) emissions from turbulent H2 jet flames, probability density function (PDF), departures from chemical equilibrium, mixing models for PDF methods, comparison of predicted and measured H2O mass fractions in turbulent nonpremixed jet flames, experimental evidence of preferential diffusion in turbulent jet flames, and computation of turbulent reacting flows.

Description: This viewgraph presentation discusses project description, turbulence models, and computational engine and results for second-order closures for compressible turbulence.

Description: This viewgraph presentation gives a profile of Advanced Scientific Computing (ASC) Ltd., applications, clients and clients' needs, ASC's directions, and how the Center for Modeling of Turbulence and Transition (CMOTT) can help.

Description: This viewgraph presentation discusses the following: STAR-CD computational features; STAR-CD turbulence models; common features of industrial complex flows; industry-specific CFD development requirements; applications and experiences of industrial complex flows, including flow in rotating disc cavities, diffusion hole film cooling, internal blade cooling, and external car aerodynamics; and conclusions on turbulence modeling needs.

Description: This viewgraph presentation provides a brief review of two-equation eddy-viscosity models (TEM's) from the perspective of applied CFD. It provides objective assessment of both well-known and newer models, compares model predictions from various TEM's with experiments, identifies sources of modeling error and gives historical perspective of their effects on model performance and assessment, and recommends directions for future research on TEM's.

Description: This viewgraph presentation discusses the following: introduction to CFD Research Corporation; experiences with two-equation models - models used, numerical difficulties, validation and applications, and strengths and weaknesses; and answers to three questions posed by the workshop organizing committee - what are your customers telling you, what are you doing in-house, and how can NASA-CMOTT (Center for Modeling of Turbulence and Transition) help.

Description: This viewgraph presentation concludes that a Monte Carlo probability density function (PDF) solution successfully couples with an existing finite volume code; PDF solution method applied to turbulent reacting flows shows good agreement with data; and PDF methods must be run on parallel machines for practical use.

Description: This viewgraph presentation discusses (1) turbulence modeling: challenges in turbulence modeling, desirable attributes of turbulence models, turbulence models in FLUENT, and examples using FLUENT; and (2) combustion modeling: turbulence-chemistry interaction and FLUENT equilibrium model. As of now, three turbulence models are provided: the conventional k-epsilon model, the renormalization group model, and the Reynolds-stress model. The renormalization group k-epsilon model has broadened the range of applicability of two-equation turbulence models. The Reynolds-stress model has proved useful for strongly anisotropic flows such as those encountered in cyclones, swirlers, and combustors. Issues remain, such as near-wall closure, with all classes of models.

Description: This viewgraph presentation discusses joint velocity-scalar PDF method; turbulent combustion modeling issues for gas turbine combustors; PDF calculations for a recirculating flow; stochastic dissipation model; joint PDF calculations for swirling flows; spray calculations; reduced kinetics/manifold methods; parallel processing; and joint PDF focus areas.

Description: The Orbiting Meteoroid and Debris Counting Experiment (OMDC) flew for approximately 90 days in a highly elliptical earth orbit onboard the Clementine Interstage Adapter (ISA) Spacecraft. This experiment obtained data on the impact flux of natural micrometeoroids and it provided limited information on the population of small mass man-made debris as a function of altitude in near earth space. The flight of the OMDC experiment on the ISA spacecraft also demonstrated that the ultra-lightweight, low-power, particle impact detector system that was used is a viable system for flights on future spacecraft to monitor the population of small mass man-made debris particles and to map the cosmic dust environment encountered on interplanetary missions. An overview of the ISA spacecraft mission, the approach to the OMDC experiment, and the data obtained by the experiment are presented.

Description: An algorithm has been developed for time-dependent forced convective diffusion-reaction having convection by a recirculating flow field within the drop that is hydrodynamically coupled at the interface with a convective external flow field that at infinity becomes a uniform free-streaming flow. The concentration field inside the droplet is likewise coupled with that outside by boundary conditions at the interface. A chemical reaction can take place either inside or outside the droplet, or reactions can take place in both phases. The algorithm has been implemented, and for comparison results are shown here for the case of no reaction in either phase and for the case of an external first order reaction, both for unsteady behavior. For pure interphase mass transfer, concentration isocontours, local and average Sherwood numbers, and average droplet concentrations have been obtained as a function of the physical properties and external flow field. For mass transfer enhanced by an external reaction, in addition to the above forms of results, we present the enhancement factor, with the results now also depending upon the (dimensionless) rate of reaction.

Description: The application of spectral methods, using a Chebyshev collocation scheme, to solve hydrodynamic stability problems is demonstrated on the Benard problem. Implementation of the Chebyshev collocation formulation is described. The performance of the spectral scheme is compared with that of a 2nd order finite difference scheme. An exact solution to the Marangoni-Benard problem is used to evaluate the performance of both schemes. The error of the spectral scheme is at least seven orders of magnitude smaller than finite difference error for a grid resolution of N = 15 (number of points used). The performance of the spectral formulation far exceeded the performance of the finite difference formulation for this problem. The spectral scheme required only slightly more effort to set up than the 2nd order finite difference scheme. This suggests that the spectral scheme may actually be faster to implement than higher order finite difference schemes.

Description: The effects of different turbulence boundary conditions were examined for two classical flows: a turbulent plane free shear layer and a flat plate turbulent boundary layer with zero pressure gradient. The flow solver used was DTNS, an incompressible Reynolds averaged Navier-Stokes solver with k-epsilon turbulence modeling, developed at the U.S. Navy David Taylor Research Center. Six different combinations of turbulence boundary conditions at the inflow boundary were investigated: In case 1, 'exact' k and epsilon profiles were used; in case 2, the 'exact' k profile was used, and epsilon was extrapolated upstream; in case 3, both k and epsilon were extrapolated; in case 4, the turbulence intensity (I) was 1 percent, and the turbulent viscosity (mu(sub t)) was equal to the laminar viscosity; in case 5, the 'exact' k profile was used and mu(sub t) was equal to the laminar viscosity; in case 6, the I was 1 percent, and epsilon was extrapolated. Comparisons were made with experimental data, direct numerical simulation results, or theoretical predictions as applicable. Results obtained with DTNS showed that turbulence boundary conditions can have significant impacts on the solutions, especially for the free shear layer.

Description: The exchange of models is one of the most serious problems currently encountered in the practice of spacecraft thermal analysis. Essentially, the problem originates in the diversity of computing environments that are used across different sites, and the consequent proliferation of native tool formats. Furthermore, increasing pressure to reduce the development's life cycle time has originated a growing interest in the so-called spacecraft concurrent engineering. In this context, the realization of the interdependencies between different disciplines and the proper communication between them become critical issues. The use of a neutral format represents a step forward in addressing these problems. Such a means of communication is adopted by consensus. A neutral format is not directly tied to any specific tool and it is kept under stringent change control. Currently, most of the groups promoting exchange formats are contributing with their experience to STEP, the Standard for Exchange of Product Model Data, which is being developed under the auspices of the International Standards Organization (ISO 10303). This paper presents the different efforts made in Europe to provide the spacecraft thermal analysis community with a Thermal Neutral Format (TNF) based on STEP. Following an introduction with some background information, the paper presents the characteristics of the STEP standard. Later, the first efforts to produce a STEP Spacecraft Thermal Application Protocol are described. Finally, the paper presents the currently harmonized European activities that follow up and extend earlier work on the area.

Description: Using global interpolation functions (GIF's) boundary element solutions are obtained for two-dimensional laminar flows. Two schemes are proposed for handling the convective terms. The first treats convection as a forcing function, and converts the flow equations to pseudo-Poisson equations. In the second scheme, some convective effect is incorporated into the fundamental solution used in constructing the pertinent integral equations. The lid-driven cavity flow is selected as the benchmark problem.