ALBERT

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: Ground-based solar observational programs are discussed. The Birmingham (U.K.) solar oscillation network (BiSON) program, the Crimean (Ukraine) low degree program, the Stanford (California) low degree program, the Tenerife (Spain) group program, the South Pole program, the Mount Wilson (California) program, the international research on the interior of the sun (IRIS) network program, the high-l helioseismometer instrument, the ground luminosity oscillations imager (LOI) program, and the Taiwanese oscillation network (TON) program, are outlined.

Description: The UBV variability of classical T Tauri stars is investigated using a large data set compiled by Herbst et al. (1994). The analysis is based on the covariance of the UBV colors, which is independent of the obscuration and the temporal sampling, and is a powerful tool for finding trends that otherwise would be hidden in the complex time-dependence of the data. I find that there are two types of UBV variability and introduce the nomenclature 'ortho-' and 'para-UBV' variations, or in short 'type-O' and 'type-P' variations. Objects with a predominance of ortho-UBV variability show strong variability in both the Balmer and Paschen continua, whereas objects dominated by para-UBV changes vary mostly in the Paschen continuum. I show that type-O variations can be explained in terms of rotating spots, variable obscuration, or changes intrinsic to the source that obey the law Delta T(sub S)/T(sub S) varies as Delta n(sub H)/n(sub H), where T(sub S) and n(sub H) are, respectively, the gas temperature and the number density of hydrogen nuclei. Type-P variations can be explained only by changes intrinsic to the source with Delta T(sub S)/T(sub S) varies as - Delta n(sub H)/n(sub H), and I speculate on a possible origin for these changes. The covariance ellipse of the UBV colors is a versatile tool for studying the connection between variability and other spectral signatures of classical T Tauri stars; its use for this purpose is illustrated by studying the correlation between the presence of CO first-overtone bands in emission and the type of variability. I find that objects with CO band emission are also those with type-P variability and propose an observational test to validate this conclusion.

Description: We present the results of a first attempt to employ multiaperture masks to obtain spectrophotometry of H II regions in nearby galaxies. A total of 97 H II regions in six southern spiral galaxies were observed using a combination of multiaperture masks and conventional long-slit spectrophotometry. The oxygen abundances derived from the multiaperture mask observations using the empirical abundance diagnostic R(sub 23) are shown to be consistent with those from long-slit spectra and generally show better reproducibility and object definition. Although the number of objects that can be observed simultaneously with this particular system is still quite limited compared with either imaging spectrophotometry or fiber-fed spectrographs, the spectral resolution offered and high throughput in the blue help make multiaperture spectrophotometry a competitive technique for increasing the sampling of H II regions in both radial distance and luminosity. There is still no clear trend of abundance gradient with either the galaxy's luminosity or its Hubble type, although the extrapolated central abundance does appear to correlate with galaxy luminosity/mass. In order to avoid difficulty in choosing an appropriate normalizing radius, we instead plot the oxygen abundance against the underlying I-band surface brightness at the radial distance of the H II region and confirm the existence of a local metallicity-surface brightness reltaionship within the disks of spiral galaxies. Although the simple closed-boc model of galaxy evolution predicts almost the right form of this relationship, a more realistic multizone model employing expnentially decreasing gas infall provides a more satisfactory fit to the observational data, provided the expected enriched gas return from dying low-mass stars shedding their envelopes at late epochs is properly taken into account. This same model, with a star formation law based upon self-regulating star formation in a three-dimensional disk (Dopita & Ryder 1994), is equally capable of accounting for the observed relationship between recent massive star formation and stellar surface brightness (Ryder & Dopita 1994).

Description: We report new high-resolution near-infrared imaging observations of the planetary nebula J900, in broadband J, H, and K, and narrowband 2.122 micrometers (for H2 1-0 S(1)), Brackett-gamma, and continuum filters. The H2 images reveal a remarkable set of outward-extending features not detected in any previous observations of theis object. In the northwest lobe of the inner bipolar region, the H2 emission, Brackett-gamma emission, and underlying continuum are spatially resolved. The relative positions of the peaks of the line emission suggest an ionization front seen edge-on. The continuum emission, apparently located nearer to the central star, is likely due to transiently heated small dust grains. The extended H2 emission features are probably excited either by oblique shocks caused by material ejected from the nebula, or by radiation leaking from holes in the opacity surrounding the ionized shell.

Description: Recent ROSAT Position Sensitive Proportional Counter (PSPC) observations for the first time unequivocally reveal the presence of a compact source of hard X-ray emission centered on the peculiar star eta Car. These observations also show a dramatic change in the hard-band (E greater than 1.6 keV) counting rate by about a factor of 2 in a 4 month interval. Thus, strong variability, which is a characteristic of eta Car in radio through IR and visible-band wavelengths, is also observed at X-ray energies. The increase in hard X-ray emission could be the result of a tripling of the mass-loss rate in less than 4 months.

Description: Using data extracted from the ROSAT archive we have discovered a recurrent supersoft X-ray transient RX J0045.4+4154 in M31. The first outburst began on 1992 February 2 and continued for at least 4 days, until the end of the observation sequence. A second outburst that lasted more than 6 days was seen to begin on 1993 January 7. The X-ray spectrum on both occasions yields a characteristic blackbody temperature of approximately 90 eV. For a range of plausible continuum models, the hydrogen column density is (0.8-1.5) x 10(exp 21)/sq cm and is consistent with the source being located in M31. This implies an unabsorbed 0.1-2.0 keV peak luminosity of approximately 10(exp 38) ergs/sec. This is the first recurrent X-ray transient to be found in M31 and is particularly notable because it is much softer than the bright X-ray transients seen in our Galaxy. The spectrum is characteristic of the supersoft class of X-ray sources, which are thought to be accreting white dwarfs that have a hydrogen-burning surface layer. A fit to a white dwarf model atmosphere gives a temperature of 10(exp 6)K, the hottest found so far. This high temperature is consistent with a white dwarf mass of 1.3-1.4 solar mass, approaching the Chandrasekhar limit, and burning close to the nuclear stability limit.

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 life of 18th century astronomer, craftsman, and partriot David Rittenhouse is detailed. As a craftsman, he distinguished himself as one of the foremost builders of clocks. He also built magnetic compasses and surveying instruments. The finest examples of his craftsmanship are considered two orreries, mechanical solar systems. In terms of astronomical observations, his best-known contribution was his observation of the transit of Venus in 1769. Rittenhouse constructed the first diffraction grating. Working as Treasurer of Pennsylvania throughout the Revolution, he became the first director of the Mint in 1792. Astronomical observations in later life included charting the position of Uranus after its discovery.

Description: Using high-resolution spectroscopic data taken with the Goddard High Resolution Spectrographic (GHRS) aboard the Hubble Space Telescope (HST) and with the International Ultraviolet Expolorer (IUE) satellite, we compare the profiles of the Mg II h and k lines seen in stars with spectral types ranging from early K through mid-M and luminosities from giants to supergiants. For all of these stars the lines are broad emission features with a central absorption. When plotted on a velocity scale the absorption features of the h and k lines agree very well in both shape and position, as do the blue wings of the emission component. The red wings of the emission, however, show a pronounced difference, with the k line wing consistently shifted to the red of the h line wing. At present the reason for this discrepancy is unknown, but we suggest several possibilities, including radiative transfer effects and high-speed stellar winds.

Description: The X-ray observatory Advanced Satellite for Cosmology and Astrophysics (ASCA) (Astro-D) was launched in Feb 1993 and is now providing excellent spectro-imaging observations of objects from virtually every astronomical class. This paper reviews the scientific capabilities of ASCA using some of the first results from X-ray binaries, cataclysmic variables, and stars.

Description: A close examination of the 0.7- to 23-micron infrared data base acquired by Gehrz and Ney (1992), suggests that the nucleus of Comet P/Halley 1986 III emitted a burst of small dust grains during a 3-day period commencing within hours of perihelion passage on 1986 February 9.46 UT. The outburst was characterized by significant increases in the coma's grain color temperature T(sub obs), temperature excess (superheat: S = T(sub obs)/T(sub BB)), infrared luminosity, albedo, and 10-micron silicate emission feature strength. These changes are all consistent with the sudden ejection from the nucleus of a cloud of grains with radii of approximately 0.5 micron. This outburst may have produced the dust that was responsible for some of the tail streamers photographed on 1986 February 22 UT. The peak of the dust outburst occurred about 3 days before a pronounced increase in the water production rate measured by the Pioneer Venus Orbiter Ultraviolet Spectrometer. We suggest that jets that release large quantities of small particles may be largely responsible for some of the variable infrared behavior that has been reported for P/Halley and other comets during the past two decades. Such jets may also account for some of the differences IR Type I and IR Type II comets.

Description: The age relations between 36 impact craters with dark paraboloids and other geologic units and structures at these localities have been studied through photogeologic analysis of Magellan SAR images of the surface of Venus. Geologic settings in all 36 sites, about 1000 x 1000 km each, could be characterized using only 10 different terrain units and six types of structures. Mapping of such units and structures in 36 randomly distributed large regions shows evidence for a distinctive regional and global stratigraphic and geologic sequence. On the basis of this sequence we have developed a model that illustrates several major themes in the history of Venus. Most of the history of Venus (that of its first 80% or so) is not preserved in the surface geomorphological record. The major deformation associated with tessera formation in the period sometime between 0.5-1.0 b.y. ago (Ivanov and Basilevsky, 1993) is the earliest event detected. Our stratigraphic analyses suggest that following tessera formation, extensive volcanic flooding resurfaced at least 85% of the planet in the form of the presently-ridged and fractured plains. Several lines of evidence favor a high flux in the post-tessera period but we have no independent evidence for the absolute duration of ridged plains emplacement. During this time, the net state of stress in the lithosphere apparently changed from extensional to compressional, first in the form of extensive ridge belt development, followed by the formation of extensive wrinkle ridges on the flow units. Subsequently, there occurred local emplacement of smooth and lobate plains units which are presently essentially undeformed. The major events in the latest 10% of the presently preserved history of Venus are continued rifting and some associated volcanism, and the redistribution of eolian material largely derived from impact crater deposits. Detailed geologic mapping and stratigraphic synthesis are necessary to test this sequence and to address many of the outstanding problems raised by this analysis.

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: We obtained the first high-dispersion IUE spectrum, and an additional low-dispersion spectrum of the very young planetary nebula He 3-1357 and its central star SDAO 244567. In 1988, it showed a very strong C IV lambda-1549 P Cygni profile which diminished and has disappeared by 1994, but C IV absorptions are seen in a high-resolution spectrum. The electron density from the C III F(lambda 1907)/F(lambda 1909) emission ratio is log N(sub e) = 4.1/cc. The older IUE data were reprocessed, and emission-line fluxes vary greatly for the period 1988-1994. The high-dispersion data show a rich Fe V and Fe VI absorption spectrum.

Description: We have observed seven narrow-linedd Seyfert 1 (NLS1) galaxies and one high-ionization Seyfert 2 galaxy with the Very Large Array (VLA). Combining these observations with published data, we summarize the radio properties of the NLS1 galaxies for which spectropolarimetry was reported by Goodrich. Fifteen of these 17 objects now have published radio observations of high sensitivity, and only nine of those have been detected. For a Hubble parameter of 75 km/s/Mpc, the 6 cm radio powers range from 10(exp 20) to 10(exp 23) W/Hz, within the range previously found for other types of Seyfert galaxy. The median radio size of the nine VLA-detected galaxies is no larger than 300 pc, similar to the median size found by Ulvestad & Wilson for a distance-limited sample of Seyfert galaxies. Of the six NLS1 galaxies known to have significant intrinsic optical polarization, three have measurable radio axes. Two of those three galaxies have radio major axes close to 90 deg from their polarization position angles, while the third has an inner radio axis that may be nearly parallel to the polarization position angle. The former relationship is expected for a Seyfert 1 in a unified model of Seyfert galaxies, assuming no intrinsic continuum polarization.

Description: We have performed multiband photopolarimetry toward stars behind the molecular cloud L1457 (MBM 12). This cloud is the nearest known molecular cloud (65 pc) and thought to be contained within the local 'hot bubble.' The polarization shows a regular structure, indicating that the cloud is threaded by an ordered magnetic field. The wavelength dependence of the polarization seems to indicate that the grains in L1457 have higher indices of refraction than normal for interstellar clouds. However, the wavelength of maximum polarization indicates that their size distribution is close to normal.

Description: Using Viking Orbiter images, detailed photoclinometric profiles were obtained across 10 irregular depressions, 32 fretted fractures, 40 troughs and pits, 124 solitary scarps, and 370 simple grabens in the north Tharsis region of Mars. These data allow inferences to be made on the shallow crustal structure of this region. The frequency modes of measured scarp heights correspond with previous general thickness estimates of the heavily cratered and rigded plains units. The depths of the flat-floored irregular depressions (55-175 m), fretted fractures (85-890 m), and troughs and pits (60-1620 m) are also similar to scarp heights (thicknesses) of the geologic units in which these depressions occur, which suggests that the depths of these flat-floored features were controlled by erosional base levels created by lithologic contacts. Although the features have a similar age, both their depths and their observed local structural control increase in the order listed above, which suggests that the more advanced stages of associated fracturing facilitated the development of these depressions by increasing permeability. If a ground-ice zone is a factor in development of these features, as has been suggested, our observation that the depths of these features decrease with increasing latitude suggests that either the thickness of the ground-ice zone does not increase poleward or the depths of the depressions were controlled by the top of the ground-ice zone whose depth may decrease with latitude.

Description: We report the detection of a broad absorption feature near 2166/cm in the spectrum of the Taurus cloud cource Elias 18. This pre-main-sequence source is the second in Taurus, the third in our survey, and the fifth known in the sky to show the broad 2166/cm absorption feature. Of equal importance, this feature is not seen toward several other embedded sources in our survey, nor is it seen toward the source Elias 16, located behind the Taurus cloud. Laboratory experiments with interstellar ice analogs show that such a feature is associated with a complex C triple bonded to N containing compound (called X(C triple bonded to N)) that results from high-energy processing (ultraviolet irradiation or ion bombardment) of simple ice components into more complex, organic components. We find a nonlinear anticorrelation between the abundance of X(C triple bonded to N) and frozen CO components in nonpolar lattices. We find no correlation between the abundance of X(C triple bonded to N) and frozen CO in polar lattices. Because the abundances of frozen CO and H2O are strongly correlated with each other and with visual extinction toward sources embedded in and located behind the Taurus molecular cloud, these ice components usually are associated with intracloud material. Our results indicate that X(C triple bonded to N) molecules result from chemical processing of dust grains dominated by nonpolar icy mantles in the local environments of pre-main-sequence stars. Such processing of icy grains in the early solar system may be an important source of organic compounds observed in minor solar system bodies. The delivery of these organic compounds to the surface of the primitive Earth through comet impacts may have provided the raw materials for prebiotic chemistry.

Description: Theoretical electron density sensitive emission-line ratios, determined using electron impact excitation rates calculated with the R-matrix code, are presented for R = I(2s(sup 2) (1)S-2s2p 3P(sub 2))/I(2s(sup 2) (1)S-2s2p 3P(sub 1) = 1(1483 A)/I(1486 A) in N IV. These are found to be up to an order of magnitude different from those deduced by previous authors, principally due to the inclusion of excitation rates for transitions among the 2s2p (3)p fine-structure levels. The observed values of R for several planetary nebulae, symbiotic stars and the Cygnus Loop supernova remnant, measured from spectra obtained with the International Ultraviolet Explorer (IUE) satellite and the Hopkins Ultraviolet Explorer (HUT), lead to electron densities which are in excellent agreement with those deduced from line ratios in other species. This provides observational support for the accuracy of the atomic data adopted in the present calculations.

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 Advanced Satellite for Cosmology and Astrophysics (ASCA) instrument uses identical conical foil X-ray mirrors for its four telescopes. One of the major advantages of ASCA's telescopes is the first time ever use of X-ray imaging over a broad energy band and high throughput for conducting astronomical spectroscopy. Nested thin foil reflectors make possible a large effective area up to 10 keV, even under the tight weight restriction of the ASCA spacecraft. The expected characteristics of the ASCA mirrors are summarized based on ray tracings and pre-flight calibrations. The total effective area of four telescopes at 1 and 7 keV are approximately 1200 and 600 sq cm, respectively. The image size (half power diameter: HPD) is about 3 min. The field of view is 24 min and 16 min full width at half maximum for 1 and 7 keV, respectively. Observations of several X-ray sources from space during the performance verification phase have confirmed those performances as well as the alignment of the optical systems.

Description: Data from the Diffuse Infrared Background Experiment (DIRBE) instrument aboard the Cosmic Background Explorer Satellite (COBE) spacecraft have been used to examine the near and far infrared signatures of the interplanetary dust (IPD) bands. Images of the dust band pairs at ecliptic latitudes of +/- 1.4 deg and +/- 10 deg have been produced at DIRBE wavelengths from 1.25 to 100 micrometers. The observations at the shorter wavelengths provide the first evidence of scattered sunlight from particles responsible for the dust bands. It is found that the grains in the bands and those in the smooth IPD cloud have similar spectral energy distributions, suggesting similar compositions and possibly a common origin. The scattering albedos from 1.25 to 3.5 micrometers for the grains in the dust bands and those in the IPD cloud are 0.22 and 0.29, respectively. The 10 deg band pair is cooler (185 +/- 10 K) than the smooth interplanetary dust cloud (259 +/- 10 K). From both parallactic and thermal analyses, the implied location of the grains responsible for the peak brightness of the 10 deg band pair is 2.1 +/- 0.1 AU the Sun A parallactic distance of 1.4 +/- 0.2 AU is found for the peak of the 1.4 deg band pair.

Description: We report a new component of Jovian radio emission observed by the Ulysses spacecraft when Ulysses was at high Jovigraphic latitudes (greater than or approximately = 30 deg north or south of the Jovian magnetic equator). This bursty high-latitude emission is elliptically polarized in the right-hand sense when observed from northern latitudes and in the left-hand sense when observed from southern latitudes, consistent with extraordinary mode. The orientation of the polarization ellipse is observed to systematically vary with time relative to the observer. It is argued that the elliptically-polarized nature of the emission is intrinsic to the source region.

Description: The high-energy diffuse gamma-ray emission from the interstellar gas in Orion was studied using observations from the Energetic Gamma-Ray Experiment Telescope (EGRET) on the Compton Gamma-Ray Observatory (CGRO) and radio surveys of the H I and CO emission. The good correlation of the gamma-ray emission with the atomic and molecular gas permits determination of the gamma-ray emissivity per nucleon in the interstellar medium and the molecular mass calibrating ratio N(H2)/W(sub co) in Orion. The integral gamma-ray emissivity is (1.7 +/- 0.1) x 10(exp -26)/s/sr for E greater than 100 MeV, in good agreement with expectations from studies of the diffuse emission on larger scales. The N(H2)/W(sub co) ratio is (1.06 +/- 0.14) x 10(exp 20)/sq cm/(K km/s), approximately 40% less than the commonly adopted Galactic average. We find no evidence for variations of the cosmic-ray density or N(H2)/W(sub co) ratio in Orion at the sensitivity and resolution of EGRET.

Description: We compare observations of BL Lac megaparsec-scale clustering environments with those for F-R I radio galaxies V and R image data are used to show that the strengths of the clustering measured via spatial cross-correlation amplitude are similar for the two types of galaxies and are consistent with Abel richness class 0. We discuss the implications this has for the origin of the BL Lac phenomena. We also examine the evolution of clustering environments of F-R I radio galaxies and BL Lac objects with redshift.

Description: We present a grid of stellar classification spectra of moderate resolution (R approximately 1500) in the range lambda lambda 5600-9000 A, compiled from high signal-to noise spectra of 275 stars, most in the open clusters Praesepe and M67. The grid covers dwarfs from types B8 through M5, giants from G8 through M7, and subgiants from F5 through K0. We catalog atomic and molecular absorption features useful for stellar classification, and demonstrate the use of luminosity-sensitive features to distinguish between late-type dwarf and giant stars. The entire database is made available in digital format on anonymous ftp and through the World Wide Web.

Description: In the past two centuries, alert amateur and professional meteor astronomers have documented 35 outbursts of 17 individual meteor streams well enough to allow the construction of a homogeneous set of activity curves. These curves add to similar profiles of the annual streams in a previous paper (Paper 1). This paper attempts to define the type and range of phenomena that classify as meteor outbursts from which the following is concluded: Outbursts are associated with the return of the comet to perihelion (near-comet type outbursts), but occur also when the parent comet is far from perihelion and far from the Earth (far-comet type). All outbursts of a given type only, depending on encounter geometry. The activity curves, expressed in terms of Zenith Hourly Rates (ZHR), have a shape that is generally well described by: ZHR = ZHR(sub max) 10(sup(-B (the absolute value of lambda (sub dot in a circle) - lambda (sup max) (sub dot in a circle))). The steepness of the slopes varies from an exponent of B = 7 to B = 220 per degree of solar longitude, with a typical value of B = 30. In addition, most near-comet type outbursts have a broader component underlying the main peak with B approximately 1 - 7.The duration Delta t is approximately 1/B of the main peak is almost independent of location near the comet, while the background component varies considerably in duration and relative intensity from one return to another. The two components in the activity curve are due to two distinct structures in the dust distribution near the parent comet, where the main component can be due to a sheet of dust that emanates from the IRAS dust trail. This brings the total number of distinct structures in meteor streams to four, including the two structures from the annual stream activity in Paper 1.

Description: On five occasions between 1992 June 29 and 1994 May 3, we have used the Australia Telescope Compact Array (ATCA) to image Eta Carinae at a wavelength of 3 cm and a resolution of 1 arcsec. These observations have revealed remarkable activity. Since 1992 June, the total flux density has increase from 0.8 to 2.2 Jy, and the original single compact source has grown to a complex of sources spread over an area of about 16 sq arcsec. Strong hydrogen recombination-line spectral emission has appeared at the site of the strongest of these new sources. This recombination emission has the largest spectral width ever observed from a star, +/- 250 km/s, and reveals gas with turbulent velocities as great as 250 km/s approaching us at an average velocity of about 200 km/s. We believe that this radio outburst has been caused by a more than threefold increase of ultraviolet luminosity, and consequent ionization of previously neutral gas clouds.

Description: The bright quasar Kaz 102, which lies in the vicinity of the North Ecliptic Pole, was monitored during the ROSAT All Sky Survey for 121.5 days from 1990 July 30 to 1991 January 25. In the course of the survey, optical photometry with various filters was peformed at several epochs, together with UV (IUE) and optical spectrophotometry. The spectral energy distribution in the 3 x 10(exp 14) -3 x 10(exp 17) Hz range is obtained simultaneously among the various frequencies to less than or = 1 day. No clear case of variability can be made in the X-rays, while in the optical and UV variability of 10%-20% is apparent. An analysis of IUE and Einstein archives indicates a doubling timescale of years for the UV and soft X-ray flux. The X-ray photon index, which in 1979 was rather flat (Gamma = 0.8(+0.6 -0.4), in 1990/1991 was found to be Gamma = 2.22 +/- 0.13, a typical value for radio-quiet quasars in this energy range. The overall energy distribution and the variability are discussed.

Description: This paper describes the calibration of the two solid-state Si(Li) X-ray detectors and the X-ray telescopes that flew as part of Goddard Space Flight Center's Broad Band X-Ray Telescope (BBXRT) experiment on board the space shuttle Columbia in 1990 December. During the 9 day shuttle mission, BBXRT performed approximately 150 observations of 82 celestial sources. The content of the archive is summarized here. Although BBXRT had a relatively short life, it stands as a milestone in X-ray astronomy as being the first instrument to offer moderate spectral resolution over a wide bandpass (0.3-12.0 keV). Among other things, this paper discusses the effective area calibration of the instrument, the flux calibration and flux corrections for off-axis observations, the detector background, and optimal background subtraction techniques. The on-axis effective area calibration for the central detector elements was performed using data from the Crab Nebula, while other carefully selected targets were used for calibration of the outer detector elements. The remaining systematic uncertainties in the effective area calibration for point sources observed both on and off axis are generally less than 5%-10%. The energy scale is known to better than 0.5% at 6.0 keV for both detectors. The results presented here have an impact on the calibration of other medium resolution X-ray experiments such as the CCDs on board ASCA (Astro-D) as well as the Advanced X-Ray Astrophysics Facility detectors.

Description: The BL Lacertae object OJ 287 was intensively monitored with the IUE satellite from 1993 March 15 to 20 in the 2000-3000 A wavelength region. The very low emission state of the source hampered the detection during part of the 23 performed exposures. The light curve at 2650 A constructed with the 11 best images exhibits a variable trend with a factor is less than or approximately 2 enhancement of the flux in about 3 hours. Simultaneous ground based observations show an optical flux variability of smaller amplitude, but correlated with the UV light curve without any apparent lag. The radio light curve is nearly constant. The UV emission state is the weakest observed in 15 years for this flaring blazard, being a factor of less than orapproximately 20 lower than the maximum recorded one.

Description: The twin Voyager spacecraft, launched nearly two decades ago, continue to operate and are now searching for the edge of our solar system, the heliopause. Voyager's giant-planet flybys of Jupiter, Saturn, Uranus, and Neptune have provided data that are likely to remain the definitive data set for the foreseeable future and have led to many ongoing discoveries. As the spacecraft move toward the heliopause, they are also providing data on the structure of the heliosphere. This article discusses the discoveries resulting from the flyby and heliosphere data that have been made within the past five years.

Description: Echelle mode spectra obtained with the Goddard High Resolution Spectrograph (GHRS) onboard the Hubble Space Telescope (HST) have been used to determine the abundances of the heavy elements Pt, Au, and Hg in the chemically peculiar HgMn-type stars kappa Cancri and chi Lupi. The abundances were determined by fitting observed line profiles with synthetically generated spectra and are found to be enhanced relative to solar system values by between three and five orders of magnitude in both stars. The Hg isotope mixture in kappa Cancri is found to resemble the terrestrial mixture while that of chi Lupi is dominated by the heaviest isotope. As determined from multiple ionization states, the abundances place constraints upon theories attempting to explain the large superficial abundances of heavy elements.

Description: The Goddard High Resoulution Spectrograph (GHRS) on the Hubble Space Telescope (HST) has been used to obtain medium (R = 20,000) and high (R = 85,000) resoultion UV spectra of chromosphere emission features for the M3.4 III star gamma Cru. Small Science Aperture (SSA) G270M and Echelle-B spectra of selected regions in the 2300-2850 A range were obtained to determine the kinematics of the chromosphere using lines of C2), Fe2, Co2, Si1/2), Ni2, Mn2, and Mg2. Profiles of C2) (UV 0.01) lines and fluorescently excited lines of low optical depth indicate average turbulent velocities (Doppler FWHM) of 30.2 +/- 1.3 and 28.8 +/- 1.3 km/s, respectively. The fluorescent emission lines (mean RV = 21.3 +/- 0.9 km/s) and the wings of the emission components of Fe2 lines (mean RV = 22.8 +/- 0.4 km/s) are approximately at rest relative to the radial velocity of the star (21 km/s), while the C2) lines show a modest inflow (mean RV = 23.1 +/- 0.9 km/s). The more opaque lines of Fe2 and Mg2 exhibit complex profiles resulting from line formation in an optically thick, extended expanding atmosphere. The emission wings of these lines are broadened by multiple scattering, and they are centered near the photospheric radial velocity. Closer to line center, these strong lines show a strong blueshifted self-absorption feature (already seen in IUE data), indicative of formation in an expanding chromosphere, and a previously unseen dip in the profiles on the red side of line center. The absorption components, when extracted using simple Gaussian fits, show strong correlations with the relative optical depths of the lines. The derived absorption flow velocities converge to the photospheric velocity as one examines spectra features formed deeper in the atmosphere. The blueward abosrption velocity increases in magnitude from about 7 to 14 km/s with increasing line optical depth - the strong absorptions directly map the acceleration of the outflowing stellar wind, while the interpretation of the weaker redshifted abosrptions is more ambiguous, indicating either an inflow of material or formation in an extended, spherically expanding outflow. The Mg2 and Fe2 profiles, taken together, imply that the wind speed decreases between the atmospheric layers where the Mg2 and Fe2 self-absorption components are formed. Interstellar absorptions are seen in the resonance lines of Mg2 (UV 1) and Fe2 (UV1) with zero-volt lower levels, at about -3 km/s, consistent with models of the interstellar medium in the direction of gamma Cru. Finally, we have detected the Mg2 'satellite lines' seen in solar spectra obtained above the limb. In gamma Cru these lines are probably fluorescently excited by H Ly beta.

Description: Scientists are still trying to piece together what exactly occurred when comet fragments from Shoemaker-Levy 9 impacted with Jupiter last year. Several theories have been from the data received from Galileo and the Hubble Space Telescope (HST) together with the modeling of the impacts using supercomputers, scientists hope to discover the answers to their questions and perhaps to resolve the differences that have arisen between the different model simulations.

Description: A brief discussion on the characteristics of pulsars is given followed by a review of millisecond pulsar discoveries including the very first, PRS B1937+21, discovered in 1982. Methods of timing millisecond pulsars and the accuracy of millisecond pulsars as clocks are discussed. Possible reasons for the pulse residuals, or differences between the observed and predicted pulse arrival times for millisecond pulsars, are given.

Description: We present data on the galactic X-ray source 1E 1740.7-2942 from the Oriented Scintillation Spectrometer Experiment (OSSE) on board NASA's Compton Gamma Ray Observatory (CGRO). Episodes of increased low-energy gamma radiation have been reported from this source, including 1-day events in 1990 October and 1992 September. These events, of intensity 7 x 10(exp -3) and 4 x 10(exp -3) photons/sq cm/s, respectively, have been interpreted as broadened and redshifted positron annihilation radiation. OSSE conducted observations of the Galactic Center region during a 21-day interval from 1992 September 17 thru 1992 October 8. This includes the time of increased 200-450 keV emission from 1E 1740.7-2942 reported by SIGMA. The OSSE observations do not confirm this event. For the specific outburst recorded by SIGMA, 1992 Sep. 19.42-20.58 (UT), OSSE data provide an upper limit (3 sigma) of 2.4 x 10(exp -3) photons/sq cm/s.

Description: X-ray observations of Abell 548 reveal that the hot gas in this cluster is at least as complex as the galaxy distribution. Our ROSAT Position Sensitive Proportional Counter (PSPC) image is used in conjunction with optical data from the Minnesota plate-scanning project and redshift data in the literature to investigate the degree of substructure in the intracluster medium and the galaxy distribution. A548 has several X-ray components: hot gas associated with clumps of galaxies, individual sources, and a newly discovered diffuse component. This diffuse component may be promordial gas falling into the cluster for the first time. The optical data suggest that this cluster consists of four major components, not three (as indicated when velocity data are ignored). Simple models of the orbital dynamics suggest that the cluster may not yet have reached its turnaround radius.

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: For the past 20 years, the 91 cm telescope in NASA's Kuiper Airborne Observatory (KAO) has enabled scientists to observe infrared sources which are obscured by the earth's atmosphere at ground-based sites, and to observe transient astronomical events from anywhere in the world. To augment this capability, the United States and German Space Agencies (NASA and DARA) are collaborating in plans to replace the KAO with a 2.5 meter telescope installed in a Boeing 747 aircraft: SOFIA - The Stratospheric Observatory for Infrared Astronomy. SOFIA's large aperture, wide wavelength coverage, mobility, accessibility, and sophisticated instruments will permit a broad range of scientific studies, some of which are described here. Its unique features complement the capabilities of other future space missions. In addition, SOFIA has important potential as a stimulus for development of new technology and as a national resource for education of K-12 teachers. If started in 1996, SOFIA will be flying in the year 2000.

Description: We have obtained images of the Orion Bar region through narrow-band filters at 6.2 and 7.8 mu m to determine whether the emission observed at these wavelengths arises from C-C modes in polycyclic aromatic hydrocarbons (PAH's). Morphologically, the distribution of the emissions are similar, but appear different in detail. Quantitative comparisons of these band intensities with previous data at 8.4 and 11.3 mu m indicates that they are indeed consistent with emission from fluorescently excited PAH molecules.

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: The nearby barred spiral galaxy NGC 1672 has been observed with the Position Sensitive Proportional Counter (PSPC) and High Resolution Imager (HRI) instruments on board the ROSAT X-ray satellite. NGC 1672 is thought to have an obscured Seyfert nucleus, and it has strong starburst activity as well. Three bright X-ray sources with luminosities 1-2 x 10(exp 40) erg/s are clearly identified with NGC 1672. The strongest lies at the nucleus, and the other two lie at the ends of NGC 1672's prominent bar, locations that are also bright in H alpha and near-infrared images. The nuclear source is resolved by the HRI on about the scale of the recently identified nuclear ring, and one of the sources at the ends of the bar is also probably resolved. The X-ray spectrum of the nuclear source is quite soft, having a Raymond-Smith plasma temperature of approximately equals 0.7 keV and little evidence for intrinsic absorption. The ROSAT band X-ray flux of the nuclear source appears to be dominated not by X-ray binary emission but rather by diffuse gas emission. The absorption and emission properties of the sources, as well as their spatial extents, lead us to models of superbubbles driven by supernovae. However, the large density and emission measure of the nuclear X-ray source stretch the limits that can be comfortably accommodated by these models. We do not detect direct emission from the putative Seyfert nucleus, although an alternative model for the nuclear source is thermal emission from gas that is photoionized by a hidden Seyfert nucleus. The spectra of the other two X-ray sources are harder than that of the nuclear source, and have similar difficulties with regard to superbubble models.

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

Description: A two dimensional finite volume method is used to predict the film coefficients in the transitional flow region (laminar or turbulent) for the radiator panel tubes. The code used to perform this analysis is CAST (Computer Aided Simulation of Turbulent Flows). The information gathered from this code is then used to augment a Sinda85 model that predicts overall performance of the radiator. A final comparison is drawn between the results generated with a Sinda85 model using the Sinda85 provided transition region heat transfer correlations and the Sinda85 model using the CAST generated data.

Description: This paper presents a probabilistic one-dimensional finite element model for heat transfer processes in porous heat exchangers. The Galerkin approach is used to develop the finite element matrices. Some of the submatrices are asymmetric due to the presence of the flow term. The Neumann expansion is used to write the temperature distribution as a series of random variables, and the expectation operator is applied to obtain the mean and deviation statistics. To demonstrate the feasibility of the formulation, a one-dimensional model of heat transfer phenomenon in superfluid flow through a porous media is considered. Results of this formulation agree well with the Monte-Carlo simulations and the analytical solutions. Although the numerical experiments are confined to parametric random variables, a formulation is presented to account for the random spatial variations.

Description: This paper presents a numerical scheme, based on the finite element method, to solve strongly coupled fluid flow and heat transfer problems. The surface radiation effect for gray, diffuse and isothermal surfaces is considered. A procedure for obtaining the view factors between the radiating surfaces is discussed. The overall solution strategy is verified by comparing the available results with those obtained using this approach. An analysis of a thermosyphon is undertaken and the effect of considering the surface radiation is clearly explained.

Description: An algorithm has been developed for the forced convective diffusion-reaction problem for convection inside and outside a droplet by a recirculating flow field hydrodynamically coupled at the droplet interface with an external flow field that at infinity becomes a uniform 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 results are shown here for the case of no reaction 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 critical state of vortex cores downstream of vortex breakdown has been studied. Base vortical flows were computed using the Reynolds-averaged, axisymmetric Navier-Stokes equations. Standard K - epsilon, RNG and second-order Reynolds stress models were employed. Results indicate that the return to supercriticality is highly dependent on the turbulence model. The K - epsilon model predicted a rapid return of the vortex to supercritical conditions, the location of which showed little sensitivity to changes in the swirl ratio. The Reynolds stress model predicted that the vortex remains subcritical to the end of the domain for each of the swirl ratios employed, and provided results in qualitative agreement with experimental work. The RNG model produced intermediate results, with a downstream movement in the critical location with increasing swirl. Calculations for which area reductions were introduced at the exit in a subcritical flow were also performed using the Reynolds stress model. The structure of the resulting recirculation zone was altered significantly. However, when area reductions were employed within supercritical flows as predicted using the two-equation models, no significant influence on the recirculation zone was noted.

Description: The view factors which are used in diffuse-gray radiation enclosure calculations are often computed by approximate numerical integrations. These approximately calculated view factors will usually not satisfy the important physical constraints of reciprocity and closure. In this paper several view-factor rectification algorithms are reviewed and a rectification algorithm based on a least-squares numerical filtering scheme is proposed with both weighted and unweighted classes. A Monte-Carlo investigation is undertaken to study the propagation of view-factor and surface-area uncertainties into the heat transfer results of the diffuse-gray enclosure calculations. It is found that the weighted least-squares algorithm is vastly superior to the other rectification schemes for the reduction of the heat-flux sensitivities to view-factor uncertainties. In a sample problem, which has proven to be very sensitive to uncertainties in view factor, the heat transfer calculations with weighted least-squares rectified view factors are very good with an original view-factor matrix computed to only one-digit accuracy. All of the algorithms had roughly equivalent effects on the reduction in sensitivity to area uncertainty in this case study.

Description: A transient, one-dimensional numerical code is developed to model the liquid motion in an axial groove with square cross section. Axial variation in liquid level, shear stress and heat transfer between the groove wall and the liquid, evaporation and transient body forces are accounted for in the model. Dryout and rewet of the groove are allowed; the front location is determined numerically using conservation of mass and linear extrapolation. Several numerical test results are presented and discussed.

Description: A higher-order finite-difference technique is developed to calculate the developing-flow field of steady incompressible laminar flows in the entrance regions of circular pipes. Navier-Stokes equations governing the motion of such a flow field are solved by using this new finite-difference scheme. This new technique can increase the accuracy of the finite-difference approximation, while also providing the option of using unevenly spaced clustered nodes for computation such that relatively fine grids can be adopted for regions with large velocity gradients. The velocity profile at the entrance of the pipe is assumed to be uniform for the computation. The velocity distribution and the surface pressure drop of the developing flow then are calculated and compared to existing experimental measurements reported in the literature. Computational results obtained are found to be in good agreement with existing experimental correlations and therefore, the reliability of the new technique has been successfully tested.

Description: The Systems Improved Numerical Differencing Analyzer and Fluid Integrator (SINDA/FLUINT) code has often been used to determine the transient and steady-state response of various thermal and fluid flow networks. While this code is an often used design and analysis tool, the validation of this program has been limited to a few simple studies. For the current study, the SINDA/FLUINT code was compared to four different analytical solutions. The thermal analyzer portion of the code (conduction and radiative heat transfer, SINDA portion) was first compared to two separate solutions. The first comparison examined a semi-infinite slab with a periodic surface temperature boundary condition. Next, a small, uniform temperature object (lumped capacitance) was allowed to radiate to a fixed temperature sink. The fluid portion of the code (FLUINT) was also compared to two different analytical solutions. The first study examined a tank filling process by an ideal gas in which there is both control volume work and heat transfer. The final comparison considered the flow in a pipe joining two infinite reservoirs of pressure. The results of all these studies showed that for the situations examined here, the SINDA/FLUINT code was able to match the results of the analytical solutions.

Description: An IBM Personal Computer (PC) version of the Groove Analysis program (GAP) was developed to predict the steady state heat transport capability of an axially grooved heat pipe for a specified groove geometry and working fluid. In the model, the capillary limit is determined by the numerical solution of the differential equation for momentum conservation with the appropriate boundary conditions. This governing equation accounts for the hydrodynamic losses due to friction in liquid and vapor flows and due to liquid/vapor shear interaction. Back-pumping in both 0-g and 1-g is accounted for in the boundary condition at the condenser end. Slug formation in 0-g and puddle flow in 1-g are also considered in the model. At the user's discretion, the code will perform the analysis for various fluid inventories (undercharge, nominal charge, overcharge, or a fixed fluid charge) and heat pipe elevations. GAP will also calculate the minimum required heat pipe wall thickness for pressure containment at design temperatures that are greater than or lower than the critical temperature of the working fluid. This paper discusses the theory behind the development of the GAP model. It also presents the many useful and powerful capabilities of the model. Furthermore, a correlation of flight test performance data and the predictions using GAP are presented and discussed.

Description: The paper presents the numerical solution of heat and mass transfer during cross-flow (orthogonal) mixed convection. In this class of flow, a buoyancy-driven transport in the vertical direction and a forced convective flow in the horizontal direction results in a three-dimensional boundary layer structure adjacent to the plate. The rates of heat and mass transfer are determined by a combined influence of the two transport processes. The equations for the conservation of mass, momentum, energy, and species concentration were solved along with appropriate boundary conditions to determine the distributions of velocity components, temperature, and concentration across the thickness of the boundary layer at different locations on the plate. Results were expressed in dimensionless form using Reynolds number, Richardson number for heat transfer, Richardson number for mass transfer, Prandtl number, and Schmidt number as parameters. It was found that the transport is dominated by buoyancy at smaller vertical locations and at larger distances away from the forced convection leading edge. Effects of forced convection appeared to be very strong at smaller horizontal distances from the leading edge. The cross stream forced convection enhanced the rate of heat and mass transfer by a very significant amount.