ALBERT

Description: The scientific requirements and implications for the instruments and telescope design for the Next Generation Space Telescope (NGST) are described. A candidate concept is a deployable, 8 m diameter telescope, optimized for the near infrared region, but featuring instruments capable of observing up to 30 micrometers. The observatory is radiatively cooled to approximately 30 K.

Description: Detailed understanding of the properties of the point-spread function (PSF) of Space Telescope Imaging Spectrograph (STIS) imaging modes is necessary to separate real effects of geometrically complex astronomical sources, such as protoplanetary disks surrounding bright stars, from instrumental effects, such as scattering and diffraction. In order to investigate STIS imaging properties we have numerically simulated broadband stellar PSFs generated by STIS in CCD clear imaging mode, including the effects of the Lyot stop and the coronagraphic wedges. The input spectrum is a stellar model atmosphere of the appropriate spectral type, convolved with the pre-flight STIS CCD response function. The PSF modeler generates broadband PSFs by co-adding weighted monochromatic PSFs across the waveband.

Description: Using the Early Release Observation of 9 Comae, we demonstrate an iterative method for correcting Space Telescope Imaging Spectrograph (STIS) echelle spectra for the effects of the echelle ripple. This analytic approach allows the actual spectrum of interest to be used in the determination of its calibration. The late F star 9 Comae is not an ideal candidate for this method, due to the many absorption lines present in its E230M spectrum, yet, given this difficulty, the method still works quite well.

Description: The installation of the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope (HST) allows for the first time two-dimensional optical and ultraviolet slitless spectroscopy of faint objects from space. The STIS Parallel Survey (SPS) routinely obtains broad band images and slitless spectra of random fields in parallel with HST observations using other instruments. The SPS is designed to study a wide variety of astrophysical phenomena, including the rate of star formation in galaxies at intermediate to high redshift through the detection of emission-line galaxies. We present the first results of the SPS, which demonstrate the capability of STIS slitless spectroscopy to detect and identify high-redshift galaxies.

Description: An analysis of the fine pointing errors of the Hubble Space Telescope (HST), in the range from 1 to 15 milliarcsecs, is reported on. The emphasis is on the study of the performance characteristics of the pointing control system, fine guidance sensors and the optical telescope assembly, which produce relative attitude and astrometric measurement errors. Since the first servicing mission in December 1993, the HST relative short term position stability is of the order of 3 milliarcsecs RMS when averaged over 1 min intervals. At this level of accuracy, longer term systematic attitude errors in this range can have a noticeable impact on the telescope's observations. The various error sources are described, including: internal temperature effects; spacecraft structure temperature effects; calibration procedures, and computational inaccuracies.

Description: The prospects for global astrometric measurements with the space interferometry mission (SIM) are discussed. The SIM mission will perform four microarcsec astrometric measurements on objects as faint as 20 mag using the optical interferometry technique with a 10 m baseline. The SIM satellite will perform narrow angle astrometry and global astrometry by means of an astrometric grid. The sensitivities of the SIM global astrometric performance and the grid accuracy versus instrumental parameters and sky coverage schemes are reported on. The problems in finding suitable astrometric grid objects to support microarcsec astrometry, and related ground-based observation programs are discussed.

Description: The space interferometry mission (SIM), for launch in 2004, is expected to determine the positions of point sources to a global accuracy of 4 microarcsec, and an accuracy of 1 microarcsec over small angles. The instrument will reach 200 mag in 3 x 10(exp 4) s and has a nominal lifetime of five years. A nulling capability will be provided and synthesis imaging and near-infrared capabilities are being considered. These capabilities will extend those of the Hipparcos mission by almost three orders of magnitude in every sense. The NASA's expectation for the scientific return from the SIM mission is outlined in terms of the areas of stellar physics, galactic structure and extragalactic astrophysics, including the distance scale problem.

Description: The burst and transient source experiment (BATSE) survey of the high energy emission from supernova remnants is reported on. The ability of BATSE to continuously monitor the entire sky in the 20 keV to 2 MeV energy range enables a large group of remnants to by studied at high energies. Preliminary analysis indicates the likely detection of several supernova remnants other than the Crab nebula. Among these are MSH 15-52, Vela, Cas A and possibly HB 9. The techniques employed are discussed together with the status of the survey and its limitations.

Description: A methodology for the simulation of turbine cooling flows is being developed. The methodology seeks to combine numerical techniques that optimize both accuracy and computational efficiency. Key components of the methodology include the use of multiblock grid systems for modeling complex geometries, and multigrid convergence acceleration for enhancing computational efficiency in highly resolved fluid flow simulations. The use of the methodology has been demonstrated in several turbo machinery flow and heat transfer studies. Ongoing and future work involves implementing additional turbulence models, improving computational efficiency, adding AMR.

Description: The following are conclusions and recomendations from the study. Primary wake effect is linear reduction in (eta) with St. Secondary wake effect is skewing of suction/pressure side cooling. Steady computations match experimental Nu, but overpredict (eta). Unsteady computations elucidate wake/film interaction. Model may be used to estimate wake passing effect. Need boundary layer and full stage experiments. Need resolved film hole and full stage unsteady computations. Need validated turbulence models for film cooling.

Description: The following topics are discussed: (1) The Transient Liquid-Crystal Heat-Transfer Technique; (2) 2-D Film-Cooling Heat-Transfer on an AlliedSignal Vane; and (3) Effects of Tab Vortex Generators on Surface Heat Transfer. Downstream of a Jet in Crossflow.

Description: A tab placed on the leeward side of the nozzle was expected to increase jet penetration into the cross-flow. An experiment at UTRC showed insignificant effect. The primary objective of the present study was to confirm and explain the ineffectiveness. The overall approach of the study was to conduct experiments in a low-speed wind tunnel and to conduct hot-wire measurements for mean velocity and streamwise vorticity fields.

Description: Testing is currently in progress to acquire simple geometry surface heat transfer data for internal channels with trips and bleed holes which can be used in the development and validation of models. The transient liquid crystal technique is used on a simple multipass model with rectangular channels and normal ribs. Normal bleed holes are located on the floor of the model in the first channel. Each hole is attached to a flow meter, allowing various bleed flow rates to simulate external pressures on the blade.

Description: Reasonable heat transfer prediction can be achieved in complex geometries. Multi-block grid allows efficient placement of grid points, and efficient use of computer resources. Wilcox k-(omega) turbulence model predicts heat transfer well, and has good numerical behavior.

Description: Data from the transient gamma ray germanium detector taken during 1995 provide a limited ability to study the high resolution spectrum of the X-ray transient GRO J1655-40. During a two-day period around the peak intensity, the soft spectrum was measured and found to be consistent with a power law spectrum with photon index -3.0 +/- 0.2. No evidence was found for narrow spectral features in the 50 to 100 keV band or around 511 keV. The 3 sigma upper limit for a narrow emission line at 511 keV is 1.2 x 10(exp -2) phot/cu cm s.

Description: We describe the STIS autonomous target acquisition capabilities. We also present the results of dedicated tests executed as part of Cycle 7 calibration, following post-launch improvements to the Space Telescope Imaging Spectrograph (STIS) flight software. The residual pointing error from the acquisitions are 〈 0.5 CCD pixels, which is better than preflight estimates. Execution of peakups show clear improvement of target centering for slits of width 0.1 sec or smaller. These results may be used by Guest Observers in planning target acquisitions for their STIS programs.

Description: Various tests have been done of the Space Telescope Imaging Spectrograph (STIS) using internal wavecals to measure thermal motion of the spectral format on the detectors. In most cases, the spectral format moves less than the specification not to exceed 0.2 pixels per hour. Primary causes of the motion are (1) changes to the thermal design dictated by the warmer Aft Shroud environment and (2) on-orbit power cycling of Multi-Anode Microchannel Arrays (MAMA) electronics to minimize the effects of radiation hits on the MAMA detectors. The rear portion of the STIS optical bench is too warm to be held at a constant temperature by internal heaters. Electronics swing in temperature with an orbital and daily frequency. The thermal drift of the optical formats is not negligible, but is well behaved in most circumstances. The observer is advised to examine the trade-off between the most accurate wavelengths with best spectral/spatial resolutions versus increased overheads that directly affect the observing times. A long term concern is that the Aft Shroud thermal environment is predicted to heat up as much as one Centigrade degree per year. Progressively more of the bench would move out of thermal control. Thus the external cooler for STIS, being considered for the Third Servicing Mission is of major importance to the long term operation of STIS.

Description: The astrophysical implications of the large width of the 1809 keV gamma ray line from the decay of radioactive Al-26, recently observed with the gamma ray imaging spectrometer (GRIS), are discussed. While there may be no apparent single mechanism that can explain the observed broadening, high speed dust grains, extremely hot superbubbles and a large, low density gaseous halo were identified in the Galactic center region as the possible origins. The implications of these observations on the hot gas content in the interstellar medium are discussed.

Description: The phenomenon of broad gamma ray line splitting is discussed, together with a scenario for the suppression of the blue wings of these broad lines due to the geometry of the accelerated particle interaction region. The broad and narrow gamma ray line emissions are compared taking into account the line splitting effect. It is concluded that the observed gamma ray lines from Orion are most likely broad, implying that the low energy cosmic rays which produce this line emission consist mostly of C and heavier ions. The suppression of the proton and alpha particle abundances requires acceleration conditions such as the acceleration of the ejecta of the supernovae before mixing with the interstellar medium. Similar conditions are implied by observations of the B and Be in low metallicity stars formed during the first Gyr of galactic evolution.

Description: The SPAcecraft SIMulator (SPASIM) simulates the functions and resources of a spacecraft to quickly perform conceptual design (Phase A) trade-off and sensitivity analyses and uncover any operational bottlenecks during any part of the mission. Failure modes and operational contingencies can be evaluated allowing operational planning (what-if scenarios) and optimization for a range of mission scenarios. The payloads and subsystems are simulated, using a hierarchy of graphical models, in terms of how their functions affect resources such as propellant, power, and data. Any of the inputs and outputs of the payloads and subsystems can be plotted during the simulation or stored in a file so they can be used by other programs. Most trade-off analyses, including those that compare current versus advanced technology, can be performed by changing values in the parameter menus. However, when a component is replaced by one with a different functional architecture, its graphical model can also be modified or replaced by drawing from a component library. SPASIM has been validated using several spacecraft designs that were at least at the Critical Design Review level. The user and programmer guide, including figures, is available on line as a hypertext document. This is an easy-to-use and expandable tool which is based on MATLAB(R) and SIMULINK(R). It runs on Silicon Graphics Inc. workstations and personal computers with Windows 95(TM) or NT(TM).

Description: It is proposed to use data from the burst and transient source experiment (BATSE) onboard the Compton Gamma Ray Observatory (CGRO) as a sensitive, rapid, all sky monitor for the International Gamma Ray Astrophysics Laboratory (INTEGRAL). The all sky earth occultation monitoring and imaging capabilities of BATSE will complement and extend the Galactic plane scans of INTEGRAL as well as provide high latitude Galactic coverage which is inaccessible to INTEGRAL. Examples of transients previously observed by BATSE, its sensitivity, and the response time for these observations are described. Examples of the variability of several of the brighter sources observed by BATSE are shown.

Description: The spectrometer for INTEGRAL (SPI) is a germanium spectrometer with a wide field of view and will provide the International Gamma Ray Astrophysics Laboratory (INTEGRAL) mission with the opportunity of studying gamma ray bursts. Simulations carried out to assess the response of the instrument using data from real burst data as input are reported on. It is shown that, despite the angular resolution of 3 deg, it is possible to locate the direction of bursts with an accuracy of a few arcmin, while offering the high spectral resolution of the germanium detectors. It is remarked that the SPI field of view is similar to the size of the halo of bursts expected around M 31 on galactic models. The detectability of bursts with such a halo is discussed.

Description: The results from observations of a sample of BL Lac objects by the Compton telescope (COMPTEL) and energetic gamma ray experiment telescope (EGRET) onboard the Compton Gamma Ray Observatory (CGRO) are presented. The main targets in the sample were selected on the basis of their X-ray brightness and apparent spectral hardening at hard X-ray energies. One of the targets, PKS 0521-365 was detected by EGRET, and these data are presented together with contemporaneous ground-based radio data. Subsequent X-ray observations of this source with the Advanced Satellite for Cosomology and Astrophysics (ASCA) demonstrated that, although originally classified as a BL Lac object, it is better described as a flat spectrum radio quasar. Upper limits on the gamma ray flux from all the non-detected BL lacs are presented.

Description: Observations of the Galactic black hole candidate binary system GX 339-4 at radio, optical and X-ray wavelengths are reviewed. The radio observations reveal a compact, persistent and variable source. On the average, the radio intensity is less than 10 mJy. These data, when compared with previous data, indicate no set correlation between the radio and X-ray emission. The X-ray intensity suggests a 14.5 month cycle. The latest radio observational data are presented together with X-ray data and are discussed in the context of similar correlated observations from other X-ray binaries.

Description: The hard X-ray behavior of several X-ray binary systems containing a neutron star or a black hole candidate is analyzed in an attempt to determine the specific signature of these categories of compact objects. Limiting the consideration to two subclasses of neutron stars, Atoll sources and non-pulsating Z sources, it appears that only the Atoll sources have a spectral behavior similar to black holes. It is proposed that Atoll sources are weakly magnetized neutron stars, whereas Z sources are small radius moderate magnetized neutron stars. Large magnetic fields funnel the accreting matter, thus preventing spherical accretion and free fall if the neutron star radius is smaller than the last stable accreting orbit. Weak magnetic fields do not have this effect, and blackbody soft photons from the stellar surface are upscattered on the relativistic infalling matter, leading to excess hard X-rays. This excess is visible in two of the observed Atoll sources and in the spectrum of a black hole candidate. In the case of a Z source, a lack of photons was remarked, providing a possible signature to distinguish between these classes of objects.

Description: The submillimeter wave astronomy satellite (SWAS) mission is dedicated to the investigation of star formation and interstellar chemistry. In order to perform the mission, SWAS will survey dense molecular clouds within the Milky Way Galaxy in either the ground state or a low-lying transition of five astrophysically-significant species: H2O, H2(18)O, O2, C I and (13)CO. The observation of these lines will: test theories that predict that these species are dominant coolants of molecular clouds during early stages of their collapse to form stars and planets, and supply information concerning the abundance of species central to the chemical models of dense interstellar gas. The SWAS will use two independent Schottky barrier diode mixers and a 53 x 68 sq cm, off-axis Cassegrain antenna.

Description: Dramatic torque reversals between spin-up and spin-down have been observed in half of the persistent X-ray pulsars monitored by the Burst and Transient Space Experiment (BATSE) all-sky monitor on the Compton Gamma Ray Observatory. Theoretical models developed to explain early pulsar timing data can explain spin-down torques via a disk-magnetosphere interaction if the star nearly corotates with the inner accretion disk. To produce the observed BATSE torque reversals, however, these equilibrium models require the disk to alternate between two mass accretion rates, with M+/- producing accretion torques of similar magnitude but always of opposite sign. Moreover, in at least one pulsar (GX 1+4) undergoing secular spin-down, the neutron star spins down faster during brief (approximately 20 day) hard X-ray flares-this is opposite the correlation expected from standard theory, assuming that BATSE pulsed flux increases with mass accretion rate. The 10 day to 10 yr intervals between torque reversals in these systems are much longer than any characteristic magnetic or viscous timescale near the inner disk boundary and are more suggestive of a global disk phenomenon. We discuss possible explanations of the observed torque behavior. Despite the preferred sense of rotation defined by the binary orbit, the BATSE observations are surprisingly consistent with an earlier suggestion for GX 1+4: the disks in these systems somehow alternate between episodes of prograde and retrograde rotation. We are unaware of any mechanism that could produce a stable retrograde disk in a binary undergoing Roche lobe overflow, but such flip-flop behavior does occur in numerical simulations of wind-fed systems. One possibility is that the disks in some of these binaries are fed by an X-ray-excited wind.

Description: The HTXS mission concept combines large effective area (approximately 15,000 sq cm at 1 keV), high spectral resolution (E/Delta(E) approximately 300-3000), and broad energy bandpass (0.25-40 keV and possibly up to 100 keV) by using replicated optics together with a complement of spectroscopic instrumentation including reflection gratings readout by charge-coupled device detectors (CCDs), quantum micro-calorimeters, and cadmium zinc telluride (CZT) or comparable high energy detectors. An essential feature of this concept involves minimization of cost (approximately $350M for development and approximately $500-600M including launches) and risk by building six identical modest satellites to achieve the large area. Current mission and technology studies are targeted towards a new start in the 2002 timeframe, with first launch around 2005-2006. The HTXS mission represents a major advance, providing as much as a factor of 100 increase in sensitivity over currently planned high resolution X ray spectroscopy missions. HTXS will mark the start of a new era when high quality X ray spectra will be obtained for all classes of X ray sources, over a wide range of luminosity and distance. With its increased capabilities, HTXS will address many fundamental astrophysics questions such as the origin and distribution of the elements from carbon to zinc, the formation and evolution of clusters of galaxies, the validity of general relativity in the strong gravity limit, the evolution of supermassive black holes in active galactic nuclei, the details of supernova explosions and their aftermath, and the mechanisms involved in the heating of stellar coronae and driving of stellar winds.

Description: Detailed flow field measurements have been carried out for a turbulent circular jet perturbed by tabs and artificial excitation. Two "delta tabs" were placed at the nozzle exit at diametrically opposite locations. The excitation condition involved subharmonic resonance that manifested in a periodic vortex pairing in the near flow field. While the excitation and the tabs independently increased jet spreading, a combination of the two diminished the effect. The jet spreading was most pronounced with the tabs but was reduced when excitation was applied to the tabbed jet. The tabs generated streamwise vortex pairs that caused a lateral spreading of the jet in a direction perpendicular to the plane containing the tabs. The excitation, on the other band, organized the azimuthal vorticity into coherent ring structures whose evolution and pairing also increased entrainment by the jet. In the tabbed case, the excitation produced coherent azimuthal structures that were distorted and asymmetric in shape. The self-induction of these structures produced an effect that opposed the tendency for the lateral spreading of the streamwise vortex pairs. The passage of the distorted vortices, and their pairing, also had a cancellation effect on the time-averaged streamwise vorticity field. These led to the reduction in jet spreading.

Description: From a sample of 7 MSFC vector magnetograms,of active regions and 17 Yohkoh SXT soft X-ray images of these active regions, we have found that the total x-ray brightness of an entire active region is correlated with the total length of neutral lines on which the magnetic field is both strong (less than 250 G) and strongly sheared (shear angle greater than 75 deg) in the same active region. This correlation, if not fortuitous, is additional evidence of the importance of strong-shear strong-field neutral lines to strong heating in active regions.

Description: This rapporteur contribution attempts to summarize the discussions at this conference in the context of current funding realities.

Description: The interaction between the main passage flow through a high pressure turbine and the secondary flows due to leakage through the wheelspace rim seals is reviewed. Various experimental and analytical studies of these interactions are outlined and a numerical investigation of the interaction between the main passage flow and a wheelspace cavity seal flow is described. The numerical investigation uses a structured grid method to study the overall interaction between the turbine stage components and the wheelspace seal flow, and an unstructured grid method to resolve the detailed flow features within the geometrically complex cavity seal. The numerical results agree with various observations from experimental studies under similar flow conditions. As the flow rate through the rim cavity seal is increased, the ingestion of fluid from the main passage flow into the rim seal area decreases rapidly. A small amount of main flow gas is ingested into the rim seal area, even at high seal flow rates.

Description: Cartesian methods for computational fluid dynamics are considered to offer an accurate and robust approach for the simulation of aerodynamic flows around geometrically complex bodies. A part of the ongoing research in this domain is reviewed with the aim of providing insight into the fundamental challenges faced by the practitioners of this approach, and a guide to further research. The integration schemes used in Cartesian solvers are similar to those used in other approaches. The emphasis is on the geometric algorithms, surface modeling and boundary conditions required to design a successful Cartesian mesh scheme.

Description: In the past forty years much progress has been made in computational methods applied to the solution of problems in spacecraft hypervelocity flow and heat transfer. Although the basic thermochemical and physical modeling techniques have changed little in this time, several orders of magnitude increase in the speed of numerically solving the Navier-Stokes and associated energy equations have been achieved. The extent to which this computational power can be applied to the design of spacecraft heat shields is dependent on the proper coupling of the external flow equations to the boundary conditions and governing equations representing the thermal protection system in-depth conduction, pyrolysis and surface ablation phenomena. A discussion of the techniques used to do this in past problems as well as the current state-of-art is provided. Specific examples, including past missions such as Galileo, together with the more recent case studies of ESA/Rosetta Sample Comet Return, Mars Pathfinder and X-33 will be discussed. Modeling assumptions, design approach and computational methods and results are presented.

Description: A loosely coupled two-phase vacuum water plume model has been developed. This model consists of a source flow model to describe the expansion of water vapor, and the Lagrangian equations of motion for particle trajectories. Gas/Particle interaction is modeled through the drag force induced by the relative velocities. Particles are assumed traveling along streamlines. The equations of motion are integrated to obtain particle velocity along the streamline. This model has been used to predict the mass flux in a 5 meter radius hemispherical domain resulting from the burst of a water jet of 1.5 mm in diameter, mass flow rate of 24.2 g/s, and stagnation pressure of 21.0 psia, which is the nominal Orbiter water dump condition. The result is compared with an empirical water plume model deduced from a video image of the STS-29 water dump. To further improve the model, work has begun to numerically simulate the bubble formation and bursting present in a liquid stream injected into a vacuum. The technique of smoothed particle hydrodynamics was used to formulate this simulation. A status and results of the on-going effort are presented and compared to results from the literature.

Description: We present the results of timing and spectral studies of the dipping X-ray source X1916-053, observed by ASCA during its Performance Verification phase. The detected dipping activity is consistent with previous observations, with a period of 3008s and an intermittent secondary dip observed roughly 0.4 out of phase with the primary dip. The energy spectra of different intensity states are fitted with a power law with partial covering fraction absorption and interstellar absorption. The increase in the hardness ratio during the primary and secondary dips, and the increase in the covering fraction and column density with decreasing X-ray intensity, all imply that the dipping is caused by the photo-absorbing materials which have been suggested to be where the accreted flow hits the outer edge of the disk materials. The spectra at all intensity levels show no apparent evidence for Fe or Ne emission lines. This may be due to the low metal abundance in the accretion flow. Alternatively, the X-ray luminosity of the central source may be too weak to excite emission lines, which are assumed to be produced by X-ray photoionization of the disk materials.

Description: We report submilliarsecond-precise astrometric measurement for the late-type star AB Doradus via a combination of VLBI (very long baseline interferometry) and Hipparos satellite data. Our astrometric analysis results in the precise determination of the kinematics of this star, that reveals an orbital motion readily explained as caused by the gravitational interaction with a low-mass companion.

Description: We present 8-13 mue spectra of 23 young Stars acquired with the UKIRT CGS3 spectrometer, including T Tauri, Herbig Ae/Be, and FU Ori Stars.

Description: We presen 8-13 meu spectra of 23 young stars acquired with the UKIRT CGS3 spectromere, including T Tauri, Herbig Ae/Be, and FU Ori stars. Silicate emission and absorption features can generally be matched with the Trapezium emissivity, by employing simple models to account for optical depth effects.

Description: The illuminating source of the photodissociation region associated with the reflection nebula NGC7023 is HD200775. We probed the foreground atomic and molecular material through diagnostics observed in absorption against the background star.

Description: Recent spectral observations of asteriods have identified a plausible origin for the ordinary chondrite meteorites (OCs). Owing to the paucity of main-belt asteriods with matching spectra [1], it had previously been difficult to explain their origin.

Description: In this paper, a pointing control performance criteria is established to suppport high resolution space spectroscopy. Results indicate that these pointing control requirements are very stringent, and would typically be difficult to meet using standard 3-axis spacecraft control. To resolve this difficulty, it is shown that performance can be significantly improved using a reconfigurable control architecture that switches among a small bankof detuned Kalman filters. The effectiveness of the control reconfiguration approach is demonstrated by example on the Space Infra-Red Telescope Facility (SIRTF) pointing system, in support of the InfraRed Spectrograph (IRS) payload.

Description: For this effort we requested ROSAT images of two small gamma-ray burst (GRB) error boxes. Our goal was to search for sources that might be associated with the quiescent site of a GRB. More than 1000 GRBs have been detected in the twenty five years since their discovery, yet their origin remains a mystery. No real-time or quiescent counterparts at any wavelength have been identified as sources of GRBs despite considerable follow-up efforts. Ground based campaigns to examine GRB error boxes shortly after the bursts have revealed no transient, or highly variable objects at optical and radio wavelengths in the time period of more than 7 hours after the burst. Due to the heavy demand on X-ray satellite time, and the difficulty of re-scheduling observations, rapid follow-up observations at high energies have not been obtained as quickly as the ground based-efforts. In fact, X-ray images of GRB error boxes are normally obtained months-years after burst detection. The current fastest X-ray response time is over two weeks. Deep imaging of GRB error boxes at X-ray wavelengths is an additional observational approach that can be used to constrain models of the origin of GRBs. We combined the results of our investigation of the GRB920525 error box with that of a prior ROSAT PSPC observation of GRB910814 and submitted them for publication. We conclude that based upon our efforts and those of others, deep X-ray observations obtained long after the GRB are not likely to yield further productive results until smaller GRB error boxes become available. It may be possible to obtain such small error boxes in the near future from XTE, HETE and an augmented IPN.

Description: A distributed algorithm for a high-order-accurate finite-difference approach to the direct numerical simulation (DNS) of transition and turbulence in compressible flows is described. This work has two major objectives. The first objective is to demonstrate that parallel and distributed-memory machines can be successfully and efficiently used to solve computationally intensive and input/output intensive algorithms of the DNS class. The second objective is to show that the computational complexity involved in solving the tridiagonal systems inherent in the DNS algorithm can be reduced by algorithm innovations that obviate the need to use a parallelized tridiagonal solver.

Description: The behavior of velocity fluctuations subjected to rapid rotation is examined. The rapid rotation considered is any arbitrary combination of two basic forms of rotation, reference frame rotation and mean flow rotation. It is recognized that the two types of rotating flows differ in the manner in which the fluctuating fields are advected. The first category is comprised of flows in rotating systems of which synoptic scale geophysical flows are a good example. In this class of flows the fluctuating velocity field advects and rotates with the mean flow. In the rapid rotation limit, the Taylor-Proudman theorem describes the behavior of this class of fluctuations. Velocity fluctuations that are advected without rotation by the mean flow constitute the second category which includes vortical flows of aerodynamic interest. The Taylor-Proudman theorem is not pertinent to I his class flows and a new result appropriate to this second category of fluctuations is derived. The present development demonstrates that the fluctuating velocity fields are rendered two-dimensional and horizontally non-divergent in the limit of any large combination of reference frame rotation and mean-flow rotation. The concommitant 'geostrophic' balance of the momentum equation is, however, dependent upon the form of rapid rotation. It is also demonstrated that the evolution equations of a two-dimensional fluctuating velocity fields are frame-indifferent with any imposed mean-flow rotation. The analyses and results of this paper highlight many fundamental aspects of rotating flows and have important consequences for their turbulence closures in inertial and non-inertial frames.

Description: A free-floating droplet in microgravity is ideal for scientific observation since it is free of confounding factors such as wetting and nonsymmetrical heat transfer introduced by contact with surfaces. However, the technology to reliably deploy in microgravity has not yet been developed. In some recent fluid deployment experiments, droplets are either shaken off the dispenser or the dispenser is quickly retracted from the droplet. These solutions impart random residual motion to deployed droplet, which can be undesirable for certain investigations. In the present study, two new types of droplet injectors were built and tested. Testing of the droplet injectors consisted of neutral buoyancy tank tests, 5-sec drop tower tests at the NASA Lewis Zero Gravity Facility, and DC-9 tests. One type, the concentric injector, worked well in the neutral buoyancy tank but did not do well in low-gravity. However, it appeared that it makes a fine apparatus for constructing bubbles in low-gravity conditions. The other type, the T-injector, showed the most promise for future development. In both neutral buoyancy and DC-9 tests, water droplets were formed and deployed with some control and repeatability, although in low-gravity the residual velocities were higher than desirable. Based on our observations, further refinements are suggested for future development work.

Description: The term renormalization and renormalization group are explained by reference to various physical systems. The extension of renormalization group to turbulence is then discussed; first as a comprehensive review and second concentrating on the technical details of a few selected approaches. We conclude with a discussion of the relevance and application of renormalization group to turbulence modelling.

Description: The results of oriented scintillation spectrometer experiment (OSSE) observations of the ultraluminous infrared galaxies Arp 220 and Mrk 273 are reported. The pointings of Arp 220 and Mrk 273 concentrated on their upper limits. The gamma ray luminosities from these sources were found to be between one and two orders of magnitude smaller than the infrared luminosities. Multiwavelength luminosity spectra are produced from the radio to the gamma ray regime, and are compared with the typical multiwavelength spectra of active galactic nuclei. The lack of measured gamma ray emission provides no evidence for the existence of buried active galactic nuclei in these ultraluminous infrared galaxies, but is consistent with an origin of the infrared luminosity from starburst activity.

Description: The Rossi X-ray Timing Explorer (RXTE) observed three active galaxies during its in-orbit verification phase: NGC 4151; NGC 4945, and MCG 8-11-11. All three were detected from 2 keV to more than 100 keV by a combination of the proportional counter array (PCA) and the high energy X-ray timing experiment (HEXTE). The PCA contains five, xenon/methane, multilayer, multiwire, gas proportional counters covering the 2 to 60 keV range, while HEXTE is an array of eight NaI/CsI phoswich scintillation counters covering the 15 to 250 keV range. The three active galaxies represent the classes of Seyfert 1, Seyfert 2 and intermediate Seyfert galaxies. The results of the fitting of various models containing partial covering fractions, Compton reflection components and high energy spectral breaks are discussed.

Description: The Compton telescope (COMPTEL) onboard the Compton Gamma Ray Observatory (CGRO) detected the pulsar PSR B1951-32 at MeV energies, and found indications of a signal from PSR B0656+14. In the combined spectra from COMPTEL and the energetic gamma ray experiment telescope (EGRET) onboard CGRO, it can be seen that the maximum luminosities of these objects are reached in the COMPTEL energy range. These spectra can be compared with those from four other pulsars observed in MeV energies with COMPTEL. The spectral properties of five of the six pulsars, Vela, PSR B1509-58, PSR B1951-32 and PSR B0656+14, require breaks and bends at MeV energies. The sixth pulsar, the Crab pulsar, approximately follows a power law flux relation from keV to GeV energies. It is concluded that this spectral behavior may play a role in the discrimination between current gamma ray emission models.

Description: Observations of the 511 keV emission in the direction of 1E 1740.7-2942 (1E) using the earth burst and transient source experiment (BATSE) onboard the Compton Gamma Ray Observatory (CGRO), are presented. The CGRO phase 1 average spectrum of 1E is calculated using a method which assumes that a given source spectrum is the sum of the flux coming directly from the object and the contribution from the surrounding diffuse emission. The 1E light curve is calculated in the 40 to 150 keV range. It presents a constant flux excess of 70 mCrab in comparison with observations from the SIGMA gamma ray telescope onboard the GRANAT observatory. By removing this contribution, the 1E spectral transition from the low state to the high standard state observed by SIGMA is confirmed, and it is shown that the 511 keV flux is independent of the 1E long term evolution from low state to high standard state. It is concluded that the 511 keV emission of (4.2 +/- 1.3) x 140(exp -4) photons/sq cm s observed in the direction of 1E is mainly diffuse and spatially extended.

Description: We report on the observation of the HI filament LVC 88+36-2 with the Deep Survey/Spectrometers of EUVE. Detecting the shadow of this cloud in absorption might have proven the existence of a diffuse EUV background behind the cloud, and constrain the emission geometries of the hot plasma in the Local Interstellar Medium (LISM). By detecting diffuse emission in the EUVE spectrometers, it would also be possible to constrain the temperature and pressure of the emitting regions as well as determine whether or not the hot phase of the LISM is in equilibrium with cosmic elemental abundances. Unfortunately, we did not detect a cloud shadow with this set of observations. This is due to a combination of circumstances, but mostly due to the fact that EUV diffuse emission is weak compared to the background levels of EUVE. If a shadow exists in the EUVE, then it is less than 1% of the EUVE detector background of 1e-3 cts arcmin(exp -2) in the Deep Survey detector.

Description: This report describes the use of novel techniques for heat containment that could be used to build a high temperature transparent furnace. The primary objective of the work was to experimentally demonstrate transparent furnace operation at 1200 C. Secondary objectives were to understand furnace operation and furnace component specification to enable the design and construction of a low power prototype furnace for delivery to NASA in a follow-up project. The basic approach of the research was to couple high temperature component design with simple concept demonstration experiments that modify a commercially available transparent furnace rated at lower temperature. A detailed energy balance of the operating transparent furnace was performed, calculating heat losses through the furnace components as a result of conduction, radiation, and convection. The transparent furnace shells and furnace components were redesigned to permit furnace operation at at least 1200 C. Techniques were developed that are expected to lead to significantly improved heat containment compared with current transparent furnaces. The design of a thermal profile in a multizone high temperature transparent furnace design was also addressed. Experiments were performed to verify the energy balance analysis, to demonstrate some of the major furnace improvement techniques developed, and to demonstrate the overall feasibility of a high temperature transparent furnace. The important objective of the research was achieved: to demonstrate the feasibility of operating a transparent furnace at 1200 C.

Description: Turbulent wall-bounded complex flows are commonly encountered in engineering practice and are of considerable interest in a variety of industrial applications. The presence of a wall significantly affects turbulence characteristics. In addition to the wall effects, turbulent wall-bounded flows become more complicated by the presence of additional body forces (e.g. centrifugal force and Coriolis force) and complex geometry. Most near-wall Reynolds stress models are developed from a high-Reynolds-number model which assumes turbulence is homogenous (or quasi-homogenous). Near-wall modifications are proposed to include wall effects in near-wall regions. In this process, wall normals are introduced. Good predictions could be obtained by Reynolds stress models with wall normals. However, ambiguity arises when the models are applied in flows with multiple walls. Many models have been proposed to model turbulent flows. Among them, Reynolds stress models, in which turbulent stresses are obtained by solving the Reynolds stress transport equations, have been proved to be the most successful ones. To apply the Reynolds stress models to wall-bounded flows, near-wall corrections accounting for the wall effects are needed, and the resulting models are called near-wall Reynolds stress models. In most of the existing near-wall models, the near-wall corrections invoke wall normals. These wall-dependent near-wall models are difficult to implement for turbulent flows with complex geometry and may give inaccurate predictions due to the ambiguity of wall normals at corners connecting multiple walls. The objective of this study is to develop a more general and flexible near-wall Reynolds stress model without using any wall-dependent variable for wall-bounded turbulent flows. With the aid of near-wall asymptotic analysis and results of direct numerical simulation, a new near-wall Reynolds stress model (NNWRS) is formulated based on Speziale et al.'s high-Reynolds-stress model with wall-independent near-wall corrections. Moreover, only one damping function is used for flows with a wide range of Reynolds numbers to ensure that the near-wall modifications diminish away from the walls.

Description: We present the results of an HST and ground-based imaging study of a complete 3C sample of zeta approx. equal to 1 sources, comprising 5 quasars and 5 radio galaxies. We have observed all of the sample in essentially line-free bands at rest-frame 0.33 micrometers with WFPC2 and in rest-frame 1 micrometer images from the ground; we have also observed most of the sample in narrow-band filters centered on [O II]. We resolve continuum structure around all of our quasars in the high-resolution WFPC2 images, and in four of the five ground-based K' images. All of the quasars have some optical continuum structure that is aligned with the radio axis. In at least 3 of these cases, some of this optical structure is directly coincident with a portion of the radio structure, including optical counterparts to radio jets in 3C212 and 3C245 and an optical counterpart to a radio lobe in 3C2. These are most likely due to optical synchrotron radiation, and the radio and optical spectral indices in the northern lobe of 3C2 are consistent with this interpretation. The fact that we see a beamed optical synchotron component in the quasars but not in the radio galaxies complicates both the magnitude and the alignment comparisons. Nonetheless, the total optical and K' flux densities of the quasar hosts are consistent with those of the radio galaxies within the observed dispersion in our sample. The distributions of K' flux densities of both radio galaxies and quasar hosts exhibit similar mean and dispersion to that found for other radio galaxies at this redshift, and the average host galaxy luminosity is equivalent to, or a little fainter than, L*. The formal determination of the alignment in the optical and infrared in the two subsamples yields no significant difference between the radio galaxy and quasar subsamples, and the quasars 3C 196 and 3C 336 have aligned continuum and emission-line structure that is probably not due to beamed optical synchrotron emission. Very blue and/or edge-brightened structures are present in some objects within the probable quasar opening angle; these are possibly the result of illumination effects from the active nucleus, i.e., scattered quasar light or photoionization. In 3C 212, we see an optical object that lies 3 min. beyond the radio lobe, but which looks morphologically quite similar to the radio lobe itself. This object is bright in the infrared and has a steep spectral gradient along its length. A striking, semi-circular arc seen associated with 3C 280 may possibly be a tidal tail from a companion, enhanced in brightness by scattering or photoionization. In the near-infrared, most of the radio galaxies have elliptical morphologies with profiles that are well-fit by de Vaucouleurs r(exp 1/4)-laws and colors that are consistent with an old stellar population. All components around the quasars have optical-infrared colors that are redder than or similar to the colors of their respective nuclei; this is more consistent with a stellar origin for the emission than with a dominant scattering contribution. From the correspondence between the total magnitudes in the galaxies and quasars and the detection of aligned components in the quasars, we conclude that this study provides general support for the unification of FR II radio galaxies and quasars. Some of the objects in the sample (e.g, 3C 212) have properties that may be difficult to explain with our current understanding of the nature of FR II radio sources and the alignment effect.

Description: Despite significant achievements in computational fluid dynamics, there still remain many fluid flow phenomena not well understood. For example, the prediction of temperature distributions is inaccurate when temperature gradients are high, particularly in shock wave turbulent boundary layer interactions close to the wall. Complexities of fluid flow phenomena include transition to turbulence, relaminarization separated flows, transition between viscous and inviscid incompressible and compressible flows, among others, in all speed regimes. The purpose of this paper is to introduce a new approach, called the Flowfield-Dependent Mixed Explicit-Implicit (FDMEI) method, in an attempt to resolve these difficult issues in Computational Fluid Dynamics (CFD). In this process, a total of six implicitness parameters characteristic of the current flowfield are introduced. They are calculated from the current flowfield or changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkoehler numbers (if reacting) at each nodal point and time step. This implies that every nodal point or element is provided with different or unique numerical scheme according to their current flowfield situations, whether compressible, incompressible, viscous, inviscid, laminar, turbulent, reacting, or nonreacting. In this procedure, discontinuities or fluctuations of an variables between adjacent nodal points are determined accurately. If these implicitness parameters are fixed to certain numbers instead of being calculated from the flowfield information, then practically all currently available schemes of finite differences or finite elements arise as special cases. Some benchmark problems to be presented in this paper will show the validity, accuracy, and efficiency of the proposed methodology.

Description: Three equilibrium-air numerical solutions are presented for the Reentry-F flight-test vehicle at Mach 20, 80,000 Ft. conditions, including turbulent flow predictions. The three solutions are from a thin-layer Navier-Stokes code, coupled thin-layer and parabolized Navier-Stokes codes, and an approximate viscous shock-layer code. Boundary-layer and shock-layer profiles are presented and compared between the three solutions, revealing close agreement between the three solution methods. Notable exceptions to the close agreement, with 7-10 percent discrepancies, occur in the density profiles at the boundary-layer edge, in the boundary-layer velocity profiles, and in the shock-layer profiles in regions influenced by the nose bluntness.

Description: Since its discovery in 1992, the superluminal X-ray transient GRS 1915+105 has been extensively observed in an attempt to understand its behaviour. We present here preliminary results from a multi-wavelength campaign undertaken from July to September 1996. This study includes X-ray data from the RXTE All Sky Monitor and BATSE, two-frequency data from the Nancay radio telescope, and infrared photometry from the 1.8m Perkins telescope at Lowell Observatory. The K-band data presented herein provide the first long-term well-sampled IR light curve of GRS 1915+105. We compare the various light curves, searching for correlations in the behaviour of the source at differing wavelengths and for possible periodicities.

Description: In this paper, we examine the effectiveness of absorbing layers as non-reflecting computational boundaries for the Euler equations. The absorbing-layer equations are simply obtained by splitting the governing equations in the coordinate directions and introducing absorption coefficients in each split equation. This methodology is similar to that used by Berenger for the numerical solutions of Maxwell's equations. Specifically, we apply this methodology to three physical problems shock-vortex interactions, a plane free shear flow and an axisymmetric jet- with emphasis on acoustic wave propagation. Our numerical results indicate that the use of absorbing layers effectively minimizes numerical reflection in all three problems considered.

Description: The objective of this investigation is experimental and computational study of three dimensional viscous flow field in the nozzle passage of an axial flow turbine stage. The nozzle passage flow field has been measured using a two sensor hot-wire probe at various axial and radial stations. In addition, two component LDV measurements at one axial station (x/c(sum m) = 0.56) were performed to measure the velocity field. Static pressure measurements and flow visualization, using a fluorescent oil technique, were also performed to obtain the location of transition and the endwall limiting streamlines. A three dimensional boundary layer code, with a simple intermittency transition model, was used to predict the viscous layers along the blade and endwall surfaces. The boundary layers on the blade surface were found to be very thin and mostly laminar, except on the suction surface downstream of 70% axial chord. Strong radial pressure gradient, especially close to the suction surface, induces strong cross flow components in the trailing edge regions of the blade. On the end-walls the boundary layers were much thicker, especially near the suction corner of the casing surface, caused by secondary flow. The secondary flow region near the suction-casing surface corner indicates the presence of the passage vortex detached from the blade surface. The corner vortex is found to be very weak. The presence of a closely spaced rotor downstream (20% of the nozzle vane chord) introduces unsteadiness in the blade passage. The measured instantaneous velocity signal was filtered using FFT square window to remove the periodic unsteadiness introduced by the downstream rotor and fans. The filtering decreased the free stream turbulence level from 2.1% to 0.9% but had no influence on the computed turbulence length scale. The computation of the three dimensional boundary layers is found to be accurate on the nozzle passage blade surfaces, away from the end-walls and the secondary flow region. On the nozzle passage endwall surfaces the presence of strong pressure gradients and secondary flow limit the validity of the boundary layer code.

Description: This is a continuing program to extend and refine the ephemeris of the Geminga pulsar with annual observations for the remaining lifetime of EGRET. The data show that every revolution of Geminga is accounted for during the EGRET epoch, and that a coherent timing solution linking the phase between EGRET, COS-B, amd SAS-2, observations has now been achieved. The accuracy of the gamma-ray timing is such that the proper motion of the pulsar can now be detected, consistent with the optical determination. The measured braking index over the 24.2 yr baseline is 17 +/- 1. Further observation is required to ascertain whether this very large braking index truly represents the energy loss mechanism, perhaps related to the theory in which Geminga is near its gamma-ray death line, or whether it is a manifestation of timing noise. Statistically significant timing residuals are detected in the EGRET data; they depart from the cubic ephemeris at a level of 23 milliperiods. The residuals appear to have a sinusoidal modulation with a period of about 5.1 yr. This could simply be a manifestation of timing noise, or it could be consistent with a planet of mass 1.7/sin i solar mass orbiting Geminga at a radius of 3.3/sin i AU.

Description: The interaction between a shock wave (attached to a wedge) and small amplitude, three-dimensional disturbances of a uniform, supersonic, freestream flow are investigated. The paper extends the two-dimensional study of Duck et al, through the use of vector potentials, which render the problem tractable by the same techniques as in the two-dimensional case, in particular by expansion of the solution by means of a Fourier-Bessel series, in appropriately chosen coordinates. Results are presented for specific classes of freestream disturbances, and the study shows conclusively that the shock is stable to all classes of disturbances (i.e. time periodic perturbations to the shock do not grow downstream), provided the flow downstream of the shock is supersonic (loosely corresponding to the weak shock solution). This is shown from our numerical results and also by asymptotic analysis of the Fourier-Bessel series, valid far downstream of the shock.

Description: The motivation for this work came from a NASA Headquarters interest in investigating design concepts for a large space telescope employing active optics technology. The development of telescope optical requirements and potential optical design configurations is reported.

Description: Measurements were made in a large scale single stage turbine facility. Within the nozzle passage measurements were made using a five hole probe, a two-component Laser Doppler Velocimeter (LDV), and a single sensor hot wire probe. These measurements showed weak secondary flows at midchord, and two secondary flow loss cores at the nozzle exit. The casing vortex loss core was the larger of the two. At the exit radial inward flow was found over the entire passage, and was more pronounced in the wake. Nozzle wake decay was found to be more rapid than for an isolated vane row due to the rotor's presence. The midspan rotor flow field was measured using a two-component LDV. Measurements were made from upstream of the rotor to a chord behind the rotor. The distortion of the nozzle wake as it passed through the rotor blade row was determined. The unsteadiness in the rotor flow field was determined. The decay of the rotor wake was also characterized.

Description: Three-dimensional, incompressible turbulent jets with rectangular and elliptical cross-sections are simulated with a finite-difference numerical method. The full Navier- Stokes equations are solved at low Reynolds numbers, whereas at high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate the effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporary discretization and a fourth-order compact scheme is used for spatial discretization. Although such methods are widely used in the simulation of compressible flows, the lack of an evolution equation for pressure or density presents particular difficulty in incompressible flows. The pressure-velocity coupling must be established indirectly. It is achieved, in this study, through a Poisson equation which is solved by a compact scheme of the same order of accuracy. The numerical formulation is validated and the dispersion and dissipation errors are documented by the solution of a wide range of benchmark problems. Three-dimensional computations are performed for different inlet conditions which model the naturally developing and forced jets. The experimentally observed phenomenon of axis-switching is captured in the numerical simulation, and it is confirmed through flow visualization that this is based on self-induction of the vorticity field. Statistical quantities such as mean velocity, mean pressure, two-point velocity spatial correlations and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stress equations are presented to aid in the turbulence modeling of complex jets. Simulations of circular jets are used to quantify the effect of the non-uniform curvature of the non-circular jets.

Description: The flow characteristics in the base region of a two-dimensional supersonic compression ramp are investigated. A stream-wise oriented air jet, M = 1.75, is injected through a thin horizontal slot into a supersonic air main flow, M = 2.3, at the end of a two-dimensional compression ramp. The velocity profile and basic characteristics of the flow in the base region immediately following the ramp are determined. Visualization of the flowfield for qualitative observations is accomplished via Dark Central Ground Interferometry (DCGI). Two-dimensional velocity profiles are obtained using Laser Doppler Velocimetry (LDV). The study is the initial phase of a four-year investigation of base flow mixing. The current study is to provide more details of the flowfield.

Description: The current compute environment that most researchers are using for the calculation of 3D unsteady Computational Fluid Dynamic (CFD) results is a super-computer class machine. The Massively Parallel Processors (MPP's) such as the 160 node IBM SP2 at NAS and clusters of workstations acting as a single MPP (like NAS's SGI Power-Challenge array and the J90 cluster) provide the required computation bandwidth for CFD calculations of transient problems. If we follow the traditional computational analysis steps for CFD (and we wish to construct an interactive visualizer) we need to be aware of the following: (1) Disk space requirements. A single snap-shot must contain at least the values (primitive variables) stored at the appropriate locations within the mesh. For most simple 3D Euler solvers that means 5 floating point words. Navier-Stokes solutions with turbulence models may contain 7 state-variables. (2) Disk speed vs. Computational speeds. The time required to read the complete solution of a saved time frame from disk is now longer than the compute time for a set number of iterations from an explicit solver. Depending, on the hardware and solver an iteration of an implicit code may also take less time than reading the solution from disk. If one examines the performance improvements in the last decade or two, it is easy to see that depending on disk performance (vs. CPU improvement) may not be the best method for enhancing interactivity. (3) Cluster and Parallel Machine I/O problems. Disk access time is much worse within current parallel machines and cluster of workstations that are acting in concert to solve a single problem. In this case we are not trying to read the volume of data, but are running the solver and the solver outputs the solution. These traditional network interfaces must be used for the file system. (4) Numerics of particle traces. Most visualization tools can work upon a single snap shot of the data but some visualization tools for transient problems require dealing with time.

Description: A code has been developed to automatically simplify full chemical mechanisms. The method employed is based on the Intrinsic Low Dimensional Manifold (ILDM) method of Maas and Pope. The ILDM method is a dynamical systems approach to the simplification of large chemical kinetic mechanisms. By identifying low-dimensional attracting manifolds, the method allows complex full mechanisms to be parameterized by just a few variables; in effect, generating reduced chemical mechanisms by an automatic procedure. These resulting mechanisms however, still retain all the species used in the full mechanism. Full and skeletal mechanisms for various fuels are simplified to a two dimensional manifold, and the resulting mechanisms are found to compare well with the full mechanisms, and show significant improvement over global one step mechanisms, such as those by Westbrook and Dryer. In addition, by using an ILDM reaction mechanism in a CID code, a considerable improvement in turn-around time can be achieved.

Description: The viscous driven-cavity problem is solved using a stream-function and vorticity formulation for the incompressible Navier-Stokes equations. This report provides the user's manual and FORTRAN code for the set of governing equations presented in NASA TM-110262.

Description: Researchers are developing the technology of 'Ballistic Particle Manufacturing' (BPM) in which individual drops are precisely layered onto a substrate, and the drops are deposited so as to prevent splatting. These individual drops will ultimately be combined to form a net-shape, three-dimensional object. Our understanding of controlled drop deposition as applied to BPM is far from complete. Process parameters include the size and temperature of the liquid metal drop, its impact velocity and trajectory, and the condition and temperature of the substrate. Quantitative knowledge of the fluid mechanics and heat transfer of drop deposition and solidification are necessary to fully optimize the manufacturing process and to control the material microstructure of the final part. The object of this study is to examine the dynamics of liquid metal drops as they impinge upon a solid surface and solidify under conditions consistent with BPM (i.e. conditions which produce non-splatting drops). A program of both numerical simulations and experiments will be conducted. Questions this study will address include the following: How do the deformation and solidification of the drop depend on the properties of the fluid drop and the solid substrate? How does the presence of previously deposited drops affect the impingement and solidification process? How does the impingement of the new drop affect already deposited material? How does the cooling rate and solidification of the drops influence the material microstructure?

Description: The transient black hole candidate GRO J1655-40 was shown elsewhere to eject matter at apparently superluminal velocities. The oriented scintillation spectrometer experiment (OSSE) spectra, measured during the X-ray outbursts, are consistent with a single power law up to 600 keV, with no evidence of a high energy cutoff. The detection of the flaring behavior of this source up to 10 MeV was claimed elsewhere. The Compton telescope (COMPTEL) observations in the 0.75 to 30 MeV band of this source are analyzed. The upper limits on the time averaged emission are presented and gamma ray flares are searched for. It is concluded that the COMPTEL data can neither confirm nor refute the claimed detection of a large gamma ray flare from this source. The gamma ray flare search in five years' worth of data did not reveal evidence for the occurrence of strong flares on timescales of two days.

Description: Soft gamma ray repeaters appear to be a new class of neutron stars. While a counterpart to SGR 0525-66 was detected uniquely in the X-ray band, SGR 1806-20 and SGR 1900+14 have unusual stellar counterparts whose spectra peak in the infrared. The infrared spectra appear to contain several components: the photospheric emission from stars dominates at shorter wavelengths; a bright point source dominates at 25 micrometers, and an extended source dominates at 60 micrometers. The longer wavelength spectra are inconsistent with mono-energetic synchrotron and black body radiation models. Recent millimeter, submillimeter and infrared observations are reviewed. A preliminary analysis of the higher resolution infrared spectra of SGR 1806-20 and SGR 1900+14 is outlined. These confirm previous observations suggesting that SGR1806-20 has an outflow and that the stars comprising the counterpart to SGR 1900+14 have very similar spectra.

Description: Results from recent observations of the Orion/Monoceros region, acquired with the Compton telescope (COMPTEL), are summarized, confirming previously reported results of excessive 3 to 7 MeV emission which is probably due to nuclear de-excitation lines from C-12 and O-16 nuclei following energetic nuclear interactions. The emphasis is on findings concerning the extent of the emission region. It is found that the emission extends over the entire cloud complex, although it may result from a few localized source regions. A preliminary spectrum is presented which indicates that the emission is widely spread over the 3 to 7 MeV range, suggesting a dominantly broadline origin from energetic C and O nuclei. The spectrum shows features that may be due to the splitting of the C and O lines.

Description: The current study expands the application of computational fluid dynamics to three-dimensional multi-element high-lift systems by investigating the flow dynamics created by a slat edge. Flow is computed over a three-element high-lift configuration using an incompressible Navier-Stokes solver with structured, overset grids processed assuming full turbulence with the one-equation Baldwin-Barth turbulence model. The geometry consists of an unswept wing, which spans the wind tunnel test section, a single element half-span Fowler flap, and a three-quarter span slat. Results are presented for the wing configured for landing with a chord based Reynolds number of 3.7 million. Results for the three-quarter span slat case are compared to the full-span slat and two-dimensional investigations.

Description: This document is a manual describing how to use the Volume Grid Manipulation (VGM) software. The code is specifically designed to alter or manipulate existing surface and volume structured grids to improve grid quality through the reduction of grid line skewness, removal of negative volumes, and adaption of surface and volume grids to flow field gradients. The software uses a command language to perform all manipulations thereby offering the capability of executing multiple manipulations on a single grid during an execution of the code. The command language can be input to the VGM code by a UNIX style redirected file, or interactively while the code is executing. The manual consists of 14 sections. The first is an introduction to grid manipulation; where it is most applicable and where the strengths of such software can be utilized. The next two sections describe the memory management and the manipulation command language. The following 8 sections describe simple and complex manipulations that can be used in conjunction with one another to smooth, adapt, and reuse existing grids for various computations. These are accompanied by a tutorial section that describes how to use the commands and manipulations to solve actual grid generation problems. The last two sections are a command reference guide and trouble shooting sections to aid in the use of the code as well as describe problems associated with generated scripts for manipulation control.

Description: On 1999 August 11, a total eclipse of the Sun will be visible from within a narrow corridor which traverses the Eastern Hemisphere. The path of the Moon's umbral shadow begins in the Atlantic and crosses central Europe, the Middle East, and India, where it ends at sunset in the Bay of Bengal. A partial eclipse will be seen within the much broader path of the Moon's penumbral shadow, which includes northeastern North America, all of Europe, northern Africa, and the western half of Asia. Detailed predictions for this event are presented and include besselian elements, geographic coordinates of the path of totality, physical ephemeris of the umbra, topocentric limb profile corrections, local circumstances for approximately 1400 cities, maps of the eclipse path, weather prospects, the lunar limb profile, and the sky during totality. Tips and suggestions are also given on how to safely view and photograph the eclipse.

Description: Jet noise is a major concern in the design of commercial aircraft. Studies by various researchers suggest that aerodynamic noise is a major contributor to jet noise. Some of these studies indicate that most of the aerodynamic jet noise due to turbulent mixing occurs when there is a rapid variation in turbulent structure, i.e. rapidly growing or decaying vortices. The objective of this research was to simulate a compressible round jet to study the non-linear evolution of vortices and the resulting acoustic radiations. In particular, to understand the effect of turbulence structure on the noise. An ideal technique to study this problem is Direct Numerical Simulations(DNS), because it provides precise control on the initial and boundary conditions that lead to the turbulent structures studied. It also provides complete 3-dimensional time dependent data. Since the dynamics of a temporally evolving jet are not greatly different from those, of a spatially evolving jet, a temporal jet problem was solved, using periodicity ill the direction of the jet axis. This enables the application of Fourier spectral methods in the streamwise direction. Physically this means that turbulent structures in the jet are repeated in successive downstream cells instead of being gradually modified downstream into a jet plume. The DNS jet simulation helps us understand the various turbulent scales and mechanisms of turbulence generation in the evolution of a compressible round jet. These accurate flow solutions will be used in future research to estimate near-field acoustic radiation by computing the total outward flux across a surface and determine how it is related to the evolution of the turbulent solutions. Furthermore, these simulations allow us to investigate the sensitivity of acoustic radiations to inlet/boundary conditions, with possible application to active noise suppression. In addition, the data generated can be used to compute various turbulence quantities such as mean velocities, turbulent stresses, etc. which will aid in turbulence modeling. This report will be presented in two chapters. The first chapter describes some work on the linear stability of a supersonic round jet and the implications of this for the jet noise problem. The second chapter is an extensive discussion of numerical work using the spectral method which we use to solve the compressible Navier-Stokes equations to study turbulent jet flows. The method uses Fourier expansions in the azimuthal and streamwise direction and a 1-D B-spline basis representation in the radial direction. The B-spline basis is locally supported and this ensures block diagonal matrix equations which can be solved in O(N) steps. This is a modification of a boundary layer code developed by Robert Moser. A very accurate highly resolved Direct Numerical Simulation (DNS) of a turbulent jet flow is produced.

Description: Jet noise is a major concern in the design of commercial aircraft. Studies by various researchers suggest that aerodynamic noise is a major contributor to jet noise. Some of these studies indicate that most of the aerodynamic jet noise due to turbulent mixing occurs when there is a rapid variation in turbulent structure, i.e. rapidly growing or decaying vortices. The objective of this research was to simulate a compressible round jet to study the non-linear evolution of vortices and the resulting acoustic radiations. In particular, to understand the effect of turbulence structure on the noise. An ideal technique to study this problem is Direct Numerical Simulations (DNS), because it provides precise control on the initial and boundary conditions that lead to the turbulent structures studied. It also provides complete 3-dimensional time dependent data. Since the dynamics of a temporally evolving jet are not greatly different from those of a spatially evolving jet, a temporal jet problem was solved, using periodicity in the direction of the jet axis. This enables the application of Fourier spectral methods in the streamwise direction. Physically this means that turbulent structures in the jet are repeated in successive downstream cells instead of being gradually modified downstream into a jet plume. The DNS jet simulation helps us understand the various turbulent scales and mechanisms of turbulence generation in the evolution of a compressible round jet. These accurate flow solutions will be used in future research to estimate near-field acoustic radiation by computing the total outward flux across a surface and determine how it is related to the evolution of the turbulent solutions. Furthermore, these simulations allow us to investigate the sensitivity of acoustic radiations to inlet/boundary conditions, with possible appli(,a- tion to active noise suppression. In addition, the data generated can be used to compute, various turbulence quantities such as mean velocities, turbulent stresses, etc. which will aid in turbulence modeling. This report will be presented in two chapters. The first chapter describes some work on the linear stability of a supersonic round jet and the implications of this for the jet noise problem. The second chapter is an extensive discussion of numerical work using the spectral method which we use to solve the compressible Navier-Stokes equations to study turbulent jet flows. The method uses Fourier expansions in the azimuthal and streamwise direction and a 1-D B-spline basis representation in the radial direction. The B-spline basis is locally supported and this ensures block diagonal matrix equations which can be solved in O(N) steps. This is a modification of a boundary layer code developed by Robert Moser. A very accurate highly resolved DNS of a turbulent jet flow is produced.

Description: A numerical method for the convective heat transfer problem is developed for low speed flow at mild temperatures. A simplified energy equation is added to the incompressible Navier-Stokes formulation by using Boussinesq approximation to account for the buoyancy force. A pseudocompressibility method is used to solve the resulting set of equations for steady-state solutions in conjunction with an approximate factorization scheme. A Neumann-type pressure boundary condition is devised to account for the interaction between pressure and temperature terms, especially near a heated or cooled solid boundary. It is shown that the present method is capable of predicting the temperature field in an incompressible flow.

Description: Boundary element algorithms for the solution of steady-state and transient heat conduction are presented. The algorithms are designed for efficient coupling with computational fluid dynamic discretizations and feature piecewise linear elements with offset nodal points. The steady-state algorithm employs the fundamental solution approach; the integration kernels are computed analytically based on linear shape functions, linear elements, and variably offset nodal points. The analytic expressions for both singular and nonsingular integrands are presented. The transient algorithm employs the transient fundamental solution; the temporal integration is performed analytically and the nonsingular spatial integration is performed numerically using Gaussian quadrature. A series solution to the integration is derived for the instance of a singular integrand. The boundary-only character of the algorithm is maintained by integrating the influence coefficients from initial time. Numerical results are compared to analytical solutions to verify the current boundary element algorithms. The steady-state and transient algorithms are numerically shown to be second-order accurate in space and time, respectively.

Description: A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Description: We analyze the spatial clustering properties of a new catalog of very rich galaxy clusters selected from the APM Galaxy Survey. These clusters are of comparable richness and space density to Abell Richness Class greater than or equal to 1 clusters, but selected using an objective algorithm from a catalog demonstrably free of artificial inhomogeneities. Evaluation of the two-point correlation function xi(sub cc)(r) for the full sample and for richer subsamples reveals that the correlation amplitude is consistent with that measured for lower richness APM clusters and X-ray selected clusters. We apply a maximum likelihood estimator to find the best fitting slope and amplitude of a power law fit to x(sub cc)(r), and to estimate the correlation length r(sub 0) (the value of r at which xi(sub cc)(r) is equal to unity). For clusters with a mean space density of 1.6 x 10(exp -6) h(exp 3) MpC(exp -3) (equivalent to the space density of Abell Richness greater than or equal to 2 clusters), we find r(sub 0) = 21.3(+11.1/-9.3) h(exp -1) Mpc (95% confidence limits). This is consistent with the weak richness dependence of xi(sub cc)(r) expected in Gaussian models of structure formation. In particular, the amplitude of xi(sub cc)(r) at all richnesses matches that of xi(sub cc)(r) for clusters selected in N-Body simulations of a low density Cold Dark Matter model.

Description: A number of wavelet-based techniques for the analysis of experimental data are developed and illustrated. A multiscale analysis based on the Mexican hat wavelet is demonstrated as a tool for acquiring physical and quantitative information not obtainable by standard signal analysis methods. Experimental data for the analysis came from simultaneous hot-wire velocity traces in a bypass transition of the boundary layer on a heated flat plate. A pair of traces (two components of velocity) at one location was excerpted. A number of ensemble and conditional statistics related to dominant time scales for energy and momentum transport were calculated. The analysis revealed a lack of energy-dominant time scales inside turbulent spots but identified transport-dominant scales inside spots that account for the largest part of the Reynolds stress. Momentum transport was much more intermittent than were energetic fluctuations. This work is the first step in a continuing study of the spatial evolution of these scale-related statistics, the goal being to apply the multiscale analysis results to improve the modeling of transitional and turbulent industrial flows.

Description: This grant was in support of a guest observation using the ASCA satellite of the young, spin-powered pulsar PSR B1706-44. The pulsar is interesting for several reasons: 1) it is young and shares many similar characteristics with the Vela pulsar, 2) it is one of a few pulsars which has been detected by the EGRET detector aboard the CGRO satellite, and 3) it is one of the confirmed sources of TeV gamma-rays discovered with ground based telescopes. The goals of the observation were to search for pulsations in the X-ray domain and to study the near stellar environment to determine if the pulsar is embedded within a compact nebula as in the case of the Vela pulsar.

Description: Loop heat pipes (LHPs) can transport very large thermal power loads over long distances, through flexible, small diameter tubes against gravitational heads. In order to overcome the evaporator limit of LHPs, which is of about 0.07 MW/sq m, work was carried out to improve the efficiency by threefold to tenfold. The vapor passage geometry for the high heat flux conditions is shown. A bidisperse wick material within the circumferential vapor passages was used. Along with heat flux enhancement, several underlying issues were demonstrated, including the fabrication of bidisperse powder with controlled properties and the fabrication of a device geometry capable of replacing vapor passages with bidisperse powder.

Description: We present GHRS/HST observations of the RS CVn binary system HD 155555. Several key UV emission lines (Fe XXI, Si IV, O V, C IV) have been analyzed to provide information about the heating rate throughout the atmosphere from the chromosphere to the corona. We show that both the G and K components reveal features of a chromosphere, transition region and corona. The emission measure distribution as a function of temperature for both components is derived and compared with the RS Cvn system, HR 1099, and the Sun. The transition region and coronal lines of both stars show nonthermal broadenings of approx. 20-30 km/s. Possible physical implications for coronal heating mechanisms are discussed.

Description: With the support of NASA's Astrophysics Data Program (NRA 92-OSSA-15), we have developed the Astronomical Software Directory Service (ASDS): a distributed, searchable, WWW-based database of software packages and their related documentation. ASDS provides integrated access to 56 astronomical software packages, with more than 16,000 URLs indexed for full-text searching. Users are performing about 400 searches per month. A new aspect of our service is the inclusion of telescope and instrumentation manuals, which prompted us to change the name to the Astronomical Software and Documentation Service. ASDS was originally conceived to serve two purposes: to provide a useful Internet service in an area of expertise of the investigators (astronomical software), and as a research project to investigate various architectures for searching through a set of documents distributed across the Internet. Two of the co-investigators were then installing and maintaining astronomical software as their primary job responsibility. We felt that a service which incorporated our experience in this area would be more useful than a straightforward listing of software packages. The original concept was for a service based on the client/server model, which would function as a directory/referral service rather than as an archive. For performing the searches, we began our investigation with a decision to evaluate the Isite software from the Center for Networked Information Discovery and Retrieval (CNIDR). This software was intended as a replacement for Wide-Area Information Service (WAIS), a client/server technology for performing full-text searches through a set of documents. Isite had some additional features that we considered attractive, and we enjoyed the cooperation of the Isite developers, who were happy to have ASDS as a demonstration project. We ended up staying with the software throughout the project, making modifications to take advantage of new features as they came along, as well as influencing the software development. The Web interface to the search engine is provided by a gateway program written in C++ by a consultant to the project (A. Warnock).

Description: A great deal of progress has been made in recent years in decomposing the 2-D structure in the atmospheres of late-type stars. Doppler images of many photospheres single stars, T Tauri stars, Algols, RS CV(sub n) binaries to name a few - are regularly published (Strassmeier 1996; Richards and Albright 1996; Rice and Strassmeier 1996; Kuerster et al. 1994). Ultraviolet spectral images of chromospheres appear in the literature (e.g., Walter et al. 1987; Neff et al. 1989) but are less common owing to the difficult nature of obtaining complete phase coverage. Zeeman doppler images of magnetic fields are now feasible (e.g., Donati et al. 1992). Performing Doppler imaging of the same targets over many seasons has also been accomplished (e.g, Vogt et al. 1997). Even when a true image reconstruction is not possible due to poor spectral resolution, we can still infer a great deal about spatial structure if enough phases are observed. However, it is increasingly apparent that to make sense of recent results, many different spectral features spanning a range of formation temperature and density must be observed simultaneously for a coherent picture to emerge. Here we report on one such campaign. In 1996, we observed the southern hemisphere RS CV(sub n) binary V824 Ara (P=1(sup d).68, G5IV+K0V-IV-IV) over one complete stellar rotation with the Hubble Space Telescope and EUVE. In conjunction, radio and optical photometry and spectroscopy were obtained from the ground. Unique to this campaign is the complete phase coverage of a number of activity proxy indicators that cover source temperatures ranging from the photosphere to the corona.

Description: The present report concerns the following projects 'High energy emission from hot accretion disks from active galactic nuclei'; 'OSSE spectral timing and monitoring observations of Cygnus X- 1'; 'OSSE observations of blazars'; and 'Search for correlated time variations of bright EGRET blazars'. The most favored model for the central engine inAGNs is the supermassive black hole hypothesis (Rees 1978). The observed luminosities range from approx. 10 logical and 44-45 erg /s for objects such as Mk 421 and Seyferts to as much as approx. 10logical and 47 for powerful QSOs such as 3C 273 and 3C 279, a large fraction of the observed bolometric luminosity being, in all likelihood, beamed. As such, these objects are strong X-ray emitters (cf. Makino et al. 1987) and often undergo gamma-ray flaring detected at the EGRET range (Hartman et al. 1992) and in the case of one of two known nearby BL Lacs, Mk 421, at TeV energies as well (Punch et al. 1992, Macomb et al. 1995). Previous campaigns emphasizing radio through X-ray and even gamma-ray observations have generally found that the multiwavelength spectrum is adequately fit by a standard synchrotron self-Compton (SSC-cf. Jones et al. 1974) model of a relativistic jet (e.g. Makino et al. 1987 and Macomb et al. 1995 for the BL Lac object Mk 421) or inhomogeneous relativistic jet (Mufson et al 1990). It also was examined the gamma-gamma transperency constraints inblazars.

Description: This is the first report of a new contract. However, this project represents ongoing work, so there are completed tasks as well as newly started tasks to report. The project involves the completion of the IRAS Galaxy Atlas (IGA), a large image database produced using data from the Infrared Astronomical Satellite (IRAS). In this phase, the project switches from the production and characterization of the IGA to its use in astronomical research studies of massive star formation. The research utilizes the IGA as well as two other large data sets being produced by research partners.

Description: EUVE, ASCA, and XTE observed the eclipsing binary Algol (Beta Per) from 1-7 Feb. 96. The coordinated observation covered approximately 2 binary orbits of the system, with a net exposure of approximately 160 ksec for EUVE, 40 ksec for ASCA (in 4 pointing), and 90 ksec for XTE (in 45 pointings). We discuss results of modeling the combined EUVE, ASCA, and XTE data using continuous differential emission measure distributions, and provide constraints on the Fe abundance in the Algol system.

Description: We obtained ASCA spectra of the Seyfert 1 galaxy NGC 3516 in March 1995. Simultaneous far-UV observations were obtained with the Hopkins Ultraviolet Telescope on the Astro-2 shuttle mission. The ASCA spectrum shows a lightly absorbed power law of energy index 0.78. The low energy absorbing column is significantly less than previously seen. Prominent 0 VII and 0 VIII absorption edges are visible, but, consistent with the much lower total absorbing column, no Fe K absorption edge is detectable. A weak, narrow Fe K(alpha) emission line from cold material is present as well as a broad Fe K(alpha) line. These features are similar to those reported in other Seyfert 1 galaxies. A single warm absorber model provides only an imperfect description of the low energy absorption. In addition to a highly ionized absorber with ionization parameter U = 1.66 and a total column density of 1.4 x 10(exp 22)/sq cm, adding a lower ionization absorber with U = 0.32 and a total column of 6.9 x 10(exp 21)/sq cm significantly improves the fit. The contribution of resonant line scattering to our warm absorber models limits the Doppler parameter to less than 160 km/s at 90% confidence. Turbulence at the sound speed of the photoionized gas provides the best fit. None of the warm absorber models fit to the X-ray spectrum can match the observed equivalent widths of all the UV absorption lines. Accounting for the X-ray and UV absorption simultaneously requires an absorbing region with a broad range of ionization parameters and column densities.

Description: This grant has supported my GHRS-related activities since 1990. This included both instrumental calibration activities and independent scientific research using the Goddard High Resolution Spectrograph on the Hubble Space Telescope. The activities under this grant are essentially complete. Publications to date which have resulted in whole or in part from this grant are included.

Description: We have been obtained high-resolution data (20 in at 50 microns and 30 in at 100 microns) on the KAO using Paul Harvey's 2 x 10 element photometer in both scanning and nodding modes. The practical flux limit for scanning is about 100 Jy. For fainter sources, a nodding (beam-switching) mode, which spends more time on the source is used. This technique has been used successfully on objects as faint as 10 Jy; the 1 sigma noise for a 1 hour integration is about 1 Jy. Although not as sensitive as space-based instruments, the higher spatial resolution afforded by the KAO is essential in studying the far-infrared emission associated with young stars; in several cases we have been able to distinguish emission from multiple sources which were blended in the IRAS beam. In addition comparison of fluxes in the KAO beam to those in the much larger IRAS beam provides information on the extended low-level emission arising from the surrounding region.

Description: For many years the work at MIT aimed at the development of new concepts and technologies for space experiments in high-energy astrophysics, but not explicitly supported by flight programs, has been supported. This work has yielded new devices and techniques for X-ray astronomy, primarily low-noise, deep-depletion charge-coupled devices (CCDS) for spectrally-resolved X-ray imaging, and high-performance transmission gratings for high-resolution X-ray spectroscopy. Among the most significant recent achievements have been the development by G. Ricker and associates of the X-ray CCD camera flying on ASCA, and currently in development for AXAF and Astro-E, and the development by C. Canizares and associates of thick, 200 nm-period transmission gratings employing the phenomenon of phase shifting for high-resolution X-ray spectroscopy up to energies of 8- 1 0 keV that is essential for the operation of the AXAF High Energy Transmission Grating Spectrometer (HETGS). Through the current SR&T grant, the latter technology is now being extended successfully to the fabrication of 100 nm-period transmission gratings, which have twice the dispersion of the AXAF gratings. We note that, among other outcomes, the modest investments of past SR&T Grants at MIT resulted in the development of the key technologies for fully one-half of the scientific instrumentation on AXAF. In addition, NASA flight programs that have benefited from previous SR&T support at MIT include the SAS 3 X-ray Observatory, which carried the first rotation modulation collimator, the Focal Plane Crystal Spectrometer (FPCS) on the Einstein Observatory, the CCD cameras on ASCA and planned for Astro-E, the High Energy Transient Experiment (HETE), the Solar EUV Monitor on the Solar and Heliospheric Observatory (SOHO), the Medium Energy Neutral Atom imager (MENA) on the Image for Magnetopause-to-aurora Global Exploration (IMAGE) mission, and the recently-approved Two Wide-Angle Imaging Neutral-atom Spectrometers (TWINS) Mission of Opportunity.

Description: A multiblock approach is presented for solving two-dimensional incompressible turbulent flows on unstructured grids. The artificial compressibility form of the governing equations is solved by a vertex-centered, finite-volume implicit scheme which uses a backward Euler time discretization. Point Gauss-Seidel relaxations are used to solve the linear system of equations at each time step. This work introduces a multiblock strategy to the solution procedure, which greatly improves the efficiency of the algorithm by significantly reducing the memory requirements while not increasing the CPU time. Results presented in this work shows that the current multiblock algorithm requires 70% less memory than the single block algorithm.