ALBERT

Description: As part of both the Early Release Observations fromthe Hubble Space Telescope and the Key PRoject on the Extragalctic Distance Scale, we have obtained multi-wavelength BVR WFPC2 images for the face-on Virgo cluster spiral galaxy M11 = NGC 4321.

Description: A detailed study was made of face and annular seals under conditions where boiling, i.e., phase change of the leaking fluid, occurs within the seal. Many seals operate in this mode because of flashing due to pressure drop and/or heat input from frictional heating. High pressure, water pumps, industrial chemical pumps, and cryogenic pumps are mentioned as a few of many applications. The initial motivation was the LOX-GOX seals for the space shuttle main engine, but the study was expanded to include any face or annular seal where boiling occurs. Some of the distinctive behavior characteristics of two-phase seals were discussed, particularly their axial stability. While two-phase seals probably exhibit instability to disturbances of other degrees of freedom such as wobble, etc., under certain conditions, such analyses are too complex to be treated at present. Since an all liquid seal (with parallel faces) has a neutral axial stiffness curve, and is stabilized axially by convergent coning, other degrees of freedom stability analyses are necessary. However, the axial stability behavior of the two-phase seal is always a consideration no matter how well the seal is aligned and regardless of the speed. Hence, axial stability is thought of as the primary design consideration for two-phase seals and indeed the stability behavior under sub-cooling variations probably overshadows other concerns. The main thrust was the dynamic analysis of axial motion of two-phase face seals, principally the determination of axial stiffness, and the steady behavior of two-phase annular seals. The main conclusions are that seals with two-phase flow may be unstable if improperly balanced. Detailed theoretical analyses of low (laminar) and high (turbulent) leakage seals are presented along with computer codes, parametric studies, and in particular a simplified PC based code that allows for rapid performance prediction. A simplified combined computer code for the performance prediction over the laminar and turbulent ranges of a two-phase seal is described and documented. The analyses, results, and computer codes are summarized.

Description: The present investigation is concerned with a fully discrete accuracy and stability analysis of the one-dimensional heat equation, taking into account the evaluation of two-pass explicit schemes which simultaneously employ lumped and coupled capacity matrices. Schemes of the considered characteristics are not amenable to uncoupled semidiscrete and ordinary differential equation analyses. The obtained results illustrate that superior behavior may be achieved by schemes of the employed type when compared with the performance of the standard one-pass explicit schemes. The key idea in the considered approach is related to the utilization of a reduced-quadrature capacity matrix in the evaluation of the right-hand-side residual.

Description: We have searched for unresolved X-ray sources in the vicinity of two rich clusters of galaxies: Abell 1995 (A1995) and MS 0451.6-0305 (MS0451), using the Chandra X-ray observatory. We detected significantly more unresolved sources around A1995 than expected based on the number of X-ray sources to the same flux limit detected in deep Chandra observations of blank fields. Previous studies have also found excess X-ray sources in the vicinity of several nearby clusters of galaxies using ROSAT (Roentgen Satellite), and recently in more distant (z is approximately 0.5) clusters (RXJ0030 and 3C295) using Chandra. In contrast, we detect only 14 unresolved X-ray sources near MS0451, which is consistent with the number expected from a cluster-free background. We determine the luminosity functions of the extra sources under the assumption that they are at the distance of their respective clusters. The characteristic luminosity of the extra sources around A1995 must be an order of magnitude fainter than that of the extra sources around RXJ0030 and 3C295. The apparent lack of extra sources around MS0451 is consistent with its greater distance and the same characteristic luminosity as the A1995 sources. Hardness ratios suggest that, on average, the extra sources in A1995 may have harder spectra than those of RXJ0030 and 3C295. These results indicate that different classes of objects may dominate in different clusters, perhaps depending on the formation history and/or dynamical state of the accompanying cluster.

Description: Following the detection of strong TeV gamma ray flares from the BL Lac object 1 ES 1959+650 with the Whipple 10 m Cerenkov telescope on 2002 May 16 and 17, we performed intensive target of opportunity radio, optical, X-ray, and TeV ?ray observations from 2002 May 18 to August 14. Observations with the X-ray telescope RossiX-Ray Timing Explorer and the Whipple and HEGRA gamma-ray telescopes revealed several strong flares, enabling us to sensitively test the X-ray--gamma-ray flux correlation properties. Although the X-ray and gamma-ray fluxes seemed to be correlated in general, we found an orphan gamma-ray flare that was not accompanied by an X-ray flare. While we detected optical flux variability with the Boltwood and Abastumani observatories, the data did not give evidence for a correlation of the optical flux variability with the observed X-ray and 7-ray flares. Within statistical errors of about 0.03 Jy at 14.5 GHz and 0.05 Jy at 4.8 GHz, the radio fluxes measured with the University of Michigan Radio Astronomy Observatory stayed constant throughout the campaign; the mean values agreed well with the values measured on 2002 May 7 and June 7 at 4.9 and 15 GHz with the Very Large Array and at 4.8 GHz with archival flux measurements. After describing in detail the radio, optical, X-ray and gamma-ray light curves, and spectral energy distributions (SEDs), we present initial modeling of the SED with a simple synchrotron self- Compton model. With the addition of another TeV blazar with good broadband data, we consider the set of all TeV blazars, to begin to look for a connection of the jet properties to the properties of the central accreting black hole thought to drive the jet. Remarkably, the temporal and special X-ray and gamma-ray emission characteristics of TeV blazars are very similar, even though the mass estimates of their central black holes differ by up to 1 order of magnitude.

Description: Recent Advanced Spacecraft for Cosmology Astrophysics (ASCA) observations of G347.3-0.5, a supernova remnant (SNR) discovered in the ROSAT All-Sky Survey, reveal nonthermal emission from a region along the northwestern shell. Here we report on new pointed ASCA observations of G347.3-0.5 that confirm this result for all the bright shell regions and also reveal similar emission, although with slightly different spectral properties, from the remainder of the SNR. Curiously, no thermal X-ray emission is detected anywhere in the remnant. We derive limits on the amount of thermal emitting material present in G347.3-0.5 and present new radio continuum, CO, and infrared results that indicate that the remnant is distant and of moderate age. We show that our observations are broadly consistent with a scenario that has most of the supernova remnant shock wave still within the stellar wind bubble of its progenitor star, while part of it appears to be interacting with denser material. A point source at the center of the remnant has spectral properties similar to those expected for a neutron star and may represent the compact relic of the supernova progenitor.

Description: We determine the distances to the z approximately equal to 0.55 galaxy clusters MS 0451.6-0305 and CL 0016+16 from a maximum likelihood joint fit to interferometric Sunyaev-Zel'dovich effect (SZE) and X-ray observations. We model the intracluster medium (ICM) using a spherical isothermal beta-model. We quantify the statistical and systematic uncertainties inherent to these direct distance measurements, and we determine constraints on the Hubble parameter for three different cosmologies. For an OmegaM = 0.3, OmegaL = 0.7 cosmology, these distances imply a Hubble constant of 63(exp 12)(sub -9)(exp +21)(sub -21) km/s/Mpc, where the uncertainties correspond to statistical followed by systematic at 68% confidence. The best fit H(sub o) is 57 km/sec/Mpc for an open OmegaM = 0.3 universe and 52 km/s/Mpc for a flat Omega = 1 universe.

Description: A solution procedure for solving nonlinear time-marching problems is presented. The nonsymmetric systems of equations arising from a Newton-type linearization of these time-marching problems are solved using an iterative strategy based on the generalized minimal residual (GMRES) algorithm. Matrix-free techniques leading to reduction in storage are presented. Incorporation of a linesearch algorithm in the Newton-GMRES scheme is discussed. An automatic time-increment control strategy is developed to increase the stability of the time-marching process. High-speed flow computations demonstrate the effectiveness of these algorithms.

Description: Much of the grant's support (and associated time) was used in developmental activity, building infrastructure for the core of the work that the grant supports. Though infrastructure development was the bulk of the activity supported this year, important progress was made in research as well. The two most important "infrastructure" items were in computing hardware and personnel. Research activities were primarily focused on improving and extending. Hughes' Teukolsky-equation-based gravitational-wave generator. Several improvements have been incorporated into this generator.

Description: The current status of streamline-upwind/Petrov-Galerkin (SUPG) methods for the analysis of flow problems is surveyed in an analytical review. Problem areas addressed include classical Galerkin, upwind, artificial-diffusion, SUPG, discontinuous Galerkin, space-time FEM, and discontinuity-capturing approaches to the scalar advection-diffusion equation; incompressible flows; advective-diffusive systems; and the compressible Euler and Navier-Stokes equations. Graphs and diagrams are provided, and the good stability properties of state-of-the-art SUPG methods are pointed out.